halbarath / eoslib Goto Github PK

EOSsolveBC is not working for ANEOSmaterial, as the lower bound for the bisection rootfinder is outside the density range of the ANEOSmaterial table. This can be easily fixed by changing the lower bound to a value inside the table, or we can switch cases and set it dynamically depending on the material. This also goes for the upper bound.

During the implementation of EOSlib in our IC code ballic I found that the following functions are missing:

• Calculation of u_cold(rho)

• A function that returns an EOS material string for a given material id

• A function that outputs the EOS material parameters or table size

We have to add the SCvH EOS (Saumon et al. 1995, Vazan et al. 2013) to EOSlib.

The following functions have to be modified:

• EOSinitMaterial()
• EOSinitIsentropicLookup()
• EOSfinalizeMaterial()
• EOSPrintMat()
• EOSPofRhoU()
• EOSPofRhoT()
• EOSCofRhoU()
• EOSPCofRhoU()
• EOSPCTofRhoU()
• EOSIsentropic()
• EOSTofRhoU()
• EOSUofRhoT()
• EOSRhoofPT()
• EOSRhoofUT() (***)
• EOSIsInTable()
• EOSdPdRho()
• EOSdPdU()
• EOSdUdRho()
• EOSdPdRhoofRhoT() (**)
• EOSdPdT() (**)
• EOSUCold()

**: These functions are currently not used but could be useful in future versions of ballic
***: These functions are not used

EOSlib.h:

• Add header file in EOSlib.h
• Define material id for the eos in EOSlib.h
• Define range of iMat for the eos in EOSlib.h
• Add material to struct EOSmaterial in EOSlib.h

Makefile:

• Add object files to objects
• Add eos and path to the library to INCLUDES