[alimea]

Hi foamers,

I want to add a non-linear equation of state (JWL) to foam to use in sonicFoam.
In this equation, I need to calculate rho, Cp, h, s, psi, Z, and cpMcv and some functions such as:

correct (which is used in EEqn.H: thermo.correct() )
pThermo(which is used in createFields.H: thermo.pThermo() )
p (which is used in createFields.H: thermo.p() )
rho (which is used in createFields.H: thermo.rho() )
he (which is used in createFieldRefs.H: thermo.he() )
psi (which is used in createFieldRefs.H: thermo.psi() )

But JWL is non-linear and calculation of parameters are not as easy as perfectGas! However, I have calculated rho using Newton method in rho() function. But I don't know how to calculate the others!

Thanks