lanl/singularity-eos

Add extrapolation to rho = 0 for spiner tables

Aperta

#541 aperta il 15 ago 2025

 (0 commenti) (1 reazione) (0 assegnatari)C++ (21 fork)auto 404
enhancementgood first issuehelp wanted

Metriche repository

Star
 (36 stelle)
Metriche merge PR
 (Metriche PR in attesa)

Descrizione

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.

Guida contributor