RedlichKwongGasLaw

more accurate gas law for modeling real gas behavior

For the use of the Redlich-Kwong gas law define the PhysicalProperties density, cv, lambda, eta and the initial absolute pressure for the material with index $Material$ by using

density($Material$) = (%MED_REDLICH_KWONG%, MolarMass, PressureCritical, TemperatureCritical) cv($Material$) = (%MED_REDLICH_KWONG%, MolarMass, PressureCritical, TemperatureCritical) lambda($Material$) = (%MED_REDLICH_KWONG%, MolarMass, PressureCritical, TemperatureCritical) eta($Material$) = (%MED_REDLICH_KWONG%, MolarMass, PressureCritical, TemperatureCritical) absolute_pressure($Material$) = InitAbsolutePressure

The parameters are:

MolarMass [g/mol] of the material, e.g. Hydrogen: 2.01588
PressureCritical [Pa]: pressure from the critical point data of the material, e.g. Hydrogen: 1.3152*10^6
TemperatureCritical [K]: temperature from the critical point data of the material, e.g. Hydrogen: 33.19

Example:

begin_alias{} "TCRIT" = "33.19" # [K] "PCRIT" = "1.3152e6" # [Pa] "Mw" = "2.01588" # [g/mol] "p0" = "93.6" # [bar] end_alias ... density($GAS$) = (%MED_REDLICH_KWONG%, [&Mw&], [&PCRIT&], [&TCRIT&]) # density in [kg/(m^3)] cv($GAS$) = (%MED_REDLICH_KWONG%, [&Mw&], [&PCRIT&], [&TCRIT&]) # heat capacity in [Nm/(kg*K)] lambda($GAS$) = (%MED_REDLICH_KWONG%, [&Mw&], [&PCRIT&], [&TCRIT&]) # heat conductivity in [W/(m*K)] eta($GAS$) = (%MED_REDLICH_KWONG%, [&Mw&], [&PCRIT&], [&TCRIT&]) # viscosity in [Pa*s] absolute_pressure($GAS$) = [&p0&*100000.0] # initial pressure in [Pa] ... begin_alias{} "wall" = " BC$...$ ACTIVE$...$ IDENT%...% MAT$GAS$ TOUCH%...% MOVE$...$ LAYER0 CHAMBER1 " end_alias

Do not forget the absolute_pressure ( see also COEFF_p_divV )!!!