SESAME Table subroutines for AstroBEAR

S2GETI, S2EOSI - These routines are to be incorporated into a hydro code which uses the inverted form of an equation of state. Density and internal energy are the independent variables, and pressure and temperature are the dependent variables.

  1. S2GETI is used to get data from the library 
CALL S2GETI (IR, IDS2, TBLS, LCNT, LU, IFL)

IR material region number

IDS2 Sesame material number

TBLS name of array designated for storage of tables

LCNT current word in array TBLS

LU unit number for library

IFL error flag
  1. Subroutine S2EOSI is used to compute an EOS point. That is, it computes the pressure, temperature, and their derivatives for a given density and internal energy. 
CALL S2EOSI (IR, TBLS, R, E, P, T)

IR material region number

TBLS name of array which contains the EOS tables

R density in Mg/m3

E internal energy in MJ/kg

P, T pressure, temperature vectors

P(1), T(1) pressure in GPa, temperature in Kelvins,

P(2), T(2) density derivatives, (P/r)E, (T/r)E

P(3), T(3) energy derivatives, (P/E)r, (T/E)r.

For certain materials, the library also has tables of the pressure, temperature, density, and internal energy along the vapor-liquid coexistence curve. This information is needed in reactor safety problems. Routines S2GET and S2GETI can be modified to access the coexistence data, and routine LA401A can be used to compute the thermodynamic quantities.

Comments

No comments.