lanl/singularity-eos

Add extrapolation to rho = 0 for spiner tables

Offen

#541 geöffnet am 15.08.2025

 (0 Kommentare) (1 Reaktion) (0 zugewiesene Personen)C++ (21 Forks)auto 404
enhancementgood first issuehelp wanted

Repository-Metriken

Stars
 (36 Sterne)
PR-Merge-Metriken
 (PR-Metriken ausstehend)

Beschreibung

Right now spiner's logarithmic interpolation means that it can never extrapolate to exactly rho = 0. An easy fix would be to check if rho < the minimum density on the table, and if it is compute pressure as $P(rho, T) = \rho \Gamma(\rho_m, T)$ where here $\rho_m$ is the minimum density and $\Gamma$ is the Gruneisen coefficient. And one could do a similar thing for temperature.

Thanks @jhp-lanl for reminding me of this issue.

Contributor Guide