eos_cell

eos_cell simply calls the equation of state on an input density, temperature, and composition, and then outputs the full thermodynamic state. This is mainly used to understand the thermodynamics one might encounter in a simulation when using a particular EOS.

Getting Started

The eos_cell code is located in Microphysics/unit_test/eos_cell. An inputs file which sets the default parameters for your thermodynamic state is needed to run the test.

Setting the thermodynamics

The parameters that affect the thermodynamics are:

  • unit_test.density : the initial density

  • unit_test.temperature : the initial temperature

The composition can be set either by specifying individual mass fractions or setting unit_test.uniform_xn as described in Defining composition.

Building and Running the Code

The code can be built simply as:

make

Note

Even though there are no reactions, a network is still required, and can be set via the NETWORK_DIR build variable. By default, the aprox13 network is used.

The network choice serves only to set the composition, and a general_null network may also be used.

The build process will automatically create links in the build directory to any required EOS table.

To run the code, in the eos_cell directory run:

./main3d.gnu.ex inputs_eos

where inputs_eos is the provided inputs file. You may edit the thermodynamic state in that file prior to running.

Output

All output is directed to stdout and simply lists the entries in the full eos_t datatype. Example output is shown below:

AMReX (22.10-75-g0a3ee1486cc3) initialized
calling the EOS on a single zone state...
rho = 1000000
T =   1000000000
xn = 1 0 0 0 0 0 0 0 0 0 0 0 0 
p = 7.33195768e+22
e = 1.284393683e+17
h = 2.017589451e+17
s = 431635027.3
dpdT = 6.581921221e+13
dpdr = 7.850159761e+16
dedT = 134065803.7
dedr = 7500364598
dhdT = 199885015.9
dhdr = 1.26823854e+10
dsdT = 0.1340658037
dsdr = -65.81921221
dpde = 490947.0602
dpdr_e = 7.481931566e+16
cv = 134065803.7
cp = 189251543
xne = 3.011083784e+29
xnp = 0
eta = 0.8073999904
pele = 5.008681351e+22
ppos = 0
mu = 1.333333333
mu_e = 2
y_e = 0.5
gam1 = 1.511401882
cs = 332889390.6
abar = 4
zbar = 2
AMReX (22.10-75-g0a3ee1486cc3) finalized