lanl/singularity-eos

Add extrapolation to rho = 0 for spiner tables

Aberta

#541 aberto em 15 de ago. de 2025

 (0 comentário) (1 reação) (0 responsável)C++ (21 forks)auto 404
enhancementgood first issuehelp wanted

Métricas do repositório

Stars
 (36 estrelas)
Métricas de merge de PR
 (Métricas PR pendentes)

Description

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.

Guia do colaborador