real-fluid based thermophysical models OpenFOAM

Installation URL: https://github.com/danhnam11/realFluidThermophysicalModels-6

A real-fluid based thermophysicalModels library OF-6

General Information

An updated thermophysicalModels library of OpenFOAM 6.0. with real-fluid models for reacting flow simulations at high pressure. Readers are referred to documents provided in the documentation directory for the detailed implementation and extension guide. They are written for the development in OpenFOAM-6 but they can be referred for the development in other versions.

List of implemented real-fluid models in the new library

Modified Soave-Redlich-Kwong (SRK) model for the equation of state [1, 2].
Peng-Robinson (PR) model for the equation of state [3]
JANAF-based model for real-fluid thermodynamic properties.
Chung’s model (1988) for dynamic viscosity and thermal conductivity [4].
Mixture averaged model for mass diffusivity of individual species in a mixture in which the binary diffusion coefficients are based on Fuller’s model and Takahashi correction at high pressure [5].
Mixture averaged model for mass diffusivity of individual species in a mixture in which the binary diffusion coefficients are based on Standard Kinetic Theory [6].

Combinations of implemented models available in the library

The runtime names of thermophysical models need to be specified correctly in the thermoType dictionary and chemistryType dictionary based on the combinations of implemented models available in the library. For example:

In the constant/thermophysicalProperties file:

thermoType
{
	type            hePsiThermo;                   //(1)
	mixture         SRKchungTakaReactingMixture;   //(2)
	transport       chungTaka;                     //(3)
	thermo          rfJanaf;                       //(4)
	energy          sensibleEnthalpy;              //(7)
	equationOfState soaveRedlichKwong;             //(5)
	specie          rfSpecie;                      //(6)
}

In the constant/chemistryProperties file:

chemistryType
{
	solver   ode;                                  // either ode/EulerImplicit/none;
	method   SRKchungTakaStandard;                 //(8)
}

There are 20 options (combinations) available for reacting flow simulations in the library as follows:

No	type(1)	mixture(2)	transport(3)	thermo(4)	EoS(5)	specie(6)	Energy(7)	method(8)
1	Psi/Rho	SRK-C-Ta	C-Ta	rfJ	SRK	rfSp	sens/Int	SRK-C-Ta-Stand
2	Psi/Rho	SRK-C-Ki	C-Ki	rfJ	SRK	rfSp	sens/Int	SRK-C-Ki-Stand
3	Psi/Rho	PR-C-Ta	C-Ta	rfJ	PR	rfSp	sens/Int	PR-C-Ta-Stand
4	Psi/Rho	PR-C-Ki	C-Ki	rfJ	PR	rfSp	sens/Int	PR-C-Ki-Stand
5	Psi/Rho	id-Ki	su-Ki	J	per	sp	sens/Int	id-Ki-Stand

where the abbreviations of model names are:

(1) Psi: hePsiThermo; Rho: heRhoThermo. (2) SRK-C-Ta: SRKchungTakaReactingMixture; SRK-C-Ki: SRKchungKineticReactingMixture; PR-C-Ta: PRchungTakaReactingMixture; PR-C-Ki: PRchungKineticReactingMixture; id-Ki: idKineticReactingMixture. (3) C-Ta: chungTaka; C-Ki: chungKinetic; su-Ki: sutherlandKinetic. (4) rfJ: rfJanaf; J: janaf. (5) SRK: soaveRedlichKwong; PR: Peng; per: perfectGas. (6) rfSp: rfSpecie; sp: specie. (7) sens: sensibleEnthalpy; Int: sensibleInternalEnergy. (8) SRK-C-Ta-Stand: SRKchungTakaStandard; SRK-C-Ki-Stand: SRKchungKineticStandard; PR-C-Ta-Stand: PRchungTakaStandard; PR-C-Ki-Stand: PRchungKineticStandard; id-Ki-Stand: idKineticStandard.