Hi,

I am working with the ESCRS(7_gas combustion model), in my opinion is the temperature calculation wrong. The coding is provided in Subroutine GXSRCS. I found this subroutine, here the temperature is calculated by the use of its ancillary subroutines, named: CKHBMS and CKCPBS. Where can I find them??? I will be really grateful if anybody can help me.

Regards Marc

Marc,

I don't know how the energy equation is solved in Pheonics but I do know the two subroutines you're trying to find.

They are CHEMKIN subroutines and are originally developed by Kee et al in Sandia NL some 20 years ago. Now CHEMKIN is managed by a private company called Reaction Design. (www.reactiondesign.com)

CKHBMS returns mass-based mixture (bulk) enthalpy given mass fractions and temperature.

CKCPBS returns mass-based mixture specific heat capacity (Cp) given mixture mass fractions and temperature.

The value obtained from CKHBMS (the last entry of the call list) should be used in the convective (advective) terms. The value from CKCPBS is normally used for transient term.

CP

Hi Marc,

You may find these routines source in the CHEMKIN library, residing in c:\phoenics\D_EARTH\D_OPT\D_CHEM\GXCHKI.FOR

Rami