thermophysicalProperties

thermoType dictionary

There are mixture, transport, thermo, equationOfState, specie and energy in thermoType.

For fluid

thermoType
{
    type            heRhoThermo;
    mixture         pureMixture; 
    transport       const;
    thermo          hConst; 
    equationOfState rhoConst; 
    specie          specie;
    energy          sensibleEnthalpy;
}

for solid

thermoType
{
    type            heSolidThermo;
    mixture         pureMixture; 
    transport       constIso; 
    thermo          hConst; 
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleEnthalpy;
}

• type : heRhoThermo for fluid, heSolidThermo for solid
• mixture : If species is not included, pureMixture, otherwise multiComponentMixture. For solid, pureMixture.
• transport : If viscosity and thermal conductivity is constant, const, if polynimial, polynomial. For solid, if thermal conductivity is constant constIso.
• thermo : hConst for specific heat is constant, hPolynomial for Polynomial.
• equationOfState : rhoConst for density is constant, perfectGas for Perfect gas, incompressiblePerfectGas for Incompressible perfect gas. icoPolynomial for Polynomial
• specie : always specie
• energy : always sensibleEnthalpy

mixture dictionary

If species is not included, material properties are set in mixture dictionary. There are thermodynamics, transport, specie, equationOfState in mixture. Example is as follows

...
mixture
{
    thermodynamics
    {
        Cp      1006.0;
        Hf      0;
    }
    transport
    {
        mu      1.79e-05;
        Pr      0.7349959183673469;
    }
    specie
    {
        nMoles      1;
        molWeight   28.966;
    }
    equationOfState
    {
        rho     1.225;
    }
}
...

thermodynamics

Set Cp and Hf. Cp is specific heat capacity and Hf is heat of formation. BaramFlow has no capability of chemical reaction yet, so Hf is not used and set as 0.

If Cp is constant, dictionary is as follows

thermodynamics
{
    Cp  <value>;
    Hf  0;
}

If Cp is Polynomial, dictionary is like this

thermodynamics
{
    Hf          0;
    Sf          0;
    CpCoeffs<8> (<a0> <a1> <a2> ... <a7>)
}

Sf is standard entropy and it also not use and set as 0.

As polynomial equation is expressed as 7th equation as follows, 8 coefficients are required. The coefficients not given are set as 0.

$Cp = a_0 + a_1 T + a_2 T^2 + a_3 T^3 + a_4 T^4 + a_5 T^5 + a_6 T^6 + a_7 T^7$

transport

Set viscosity, mu and Prandtl Number, Pr.

If viscosity and thermal conductivity is constannt,

transport
{
    mu  <value>;
    Pr  <value>;
}

If viscosity and thermal conductivity is polynomial,

transport
{
    muCoeffs<8>     (<a0> <a1> <a2> ... <a7>)
    kappaCoeffs<8>  (<a0> <a1> <a2> ... <a7>)
}

If viscosity is Sutherland,

transport
{
    As  <value>;
    Ts  <value>;
}

As is Sutherland coefficient, Ts is Sutherland temperature

specie

Set molecular weight at specie

specie
{
    nMoles      1;
    molWeight   <value>;
}

equationOfState

If density is constant,

equationOfState
{
    rho <value>;
}

If density is polynomial

equationOfState
{
    rhoCoeffs<8>     (<a0> <a1> <a2> ... <a7>)
}

If density is perfect gas, equationOfState dictionary is not needed.

if density is incomprssible perect gas,

equationOfState
{
    pRef    <value>;
}

Pressure at [Reference] is used at pRef.

transportProperties

Define material properties for multi-phase flow. There must be phases dictionary and additional dictionaries for each phase.

Surface tension and Cavitation is set here.

phases

For 2 pahse

phases  (<phase1> <phase2>);

For more than 2 pahse

phases
(
    <phase1> 
    {
        transportModel  Newtonian;
        nu              <value>;
        rho             <value>;
    } 
    <phase2>
    {
        transportModel  Newtonian;
        nu              <value>;
        rho             <value>;
    }
    ... 
);

Properties for each phase

For 2 phase

<phase1>
{
    transportModel  Newtonian;
    nu              <value>;
    rho             <value>;
}
<phase2>
{
    transportModel  Newtonian;
    nu              <value>;
    rho             <value>;
}

For more than 2 phase, noting is needed because properties are set in phases.

Surface tension

For 2 phase

sigma   <value>;

For more than 2 phase, values for each other are needed as follows

sigmas
(
    (<phase1> <phase2>) <value>
    (<phase1> <phase3>) <value>
    (<phase2> <phase3>) <value>
    ...
);

Cavitation

Set vapor pressure at pSat and cavitation model at phaseChangeTwoPhaseMixture

pSat                        <value>;
phaseChangeTwoPhaseMixture  SchnerrSauer; //Kunz, Merkle, Zwart

Coefficients for each model must be set as dictionary as follows

Schnerr-Sauer model

SchnerrSauerCoeffs
{
    n       <value>;
    dNuc    <value>;
    Cc      <value>;
    Cv      <value>;
}

Kunz, Merkle model

KunzCoeffs // MerkleCoeffs
{
    UInf    <value>;
    tInf    <value>;
    Cc      <value>;
    Cv      <value>;
}

Zwart-Gerber-Belamri model

ZwartCoeffs
{
    aNuc    <value>;
    dNuc    <value>;
    Cc      <value>;
    Cv      <value>;
}