Conjugate heat transfer with real gas EOS

aalzoubi06 · October 9, 2019, 04:23

Hello All,
I am wondering if there is a possibility to run conjugate heat transfer such as chtMultiRegionFoam using a real gas Equation of State such as Peng–Robinson.

Many thanks in advance.

dasith0001 · August 18, 2022, 03:58

Very keen to know if this is possible at all, adapting Peng-Robinson into chtmultiRegionFoam

jherb · August 25, 2022, 12:08

Yes, I used Peng-Robinson in a simulation with chtMultiRegionFoam and it worked.

dasith0001 · September 20, 2022, 04:13

Quote:

Originally Posted by jherb

Yes, I used Peng-Robinson in a simulation with chtMultiRegionFoam and it worked.

Which version are you working with? Can you provide more details of the thermophysical file ?

Thank you
Dasith

Yann · September 20, 2022, 04:26

Hi,

As long as Peng-Robinson is implemented in the thermophysical models, it should be possible to use it with chtMultiRegionFoam : https://www.openfoam.com/documentati...onOfState.html

Searching the forum might also help to find answers: Usage of PengRobinsonGas

I hope this helps,
Yann

dasith0001 · September 20, 2022, 21:16

Quote:

Originally Posted by Yann

Hi,

As long as Peng-Robinson is implemented in the thermophysical models, it should be possible to use it with chtMultiRegionFoam : https://www.openfoam.com/documentati...onOfState.html

Searching the forum might also help to find answers: Usage of PengRobinsonGas

I hope this helps,
Yann

I used the following thernophysical file

PHP Code:


			
thermoType

{

    type            heRhoThermo;

    mixture         pureMixture;

    transport       const;

    thermo          hConst;

    equationOfState PengRobinsonGas;

    specie          specie;

    energy          sensibleEnthalpy;

}





mixture

{

    specie

    {

        molWeight       4.002602; //g/mol

    }



    equationOfState

    {

        Tc 126.19;

        Vc 0.089500;

        Pc 3397800;

        omega 0.04;

    //    rhoCoeffs<8>    ( 3.8936 -0.016463 3.2101e-05 -2.9174e-08 9.9889e-12 0 0 0 );

    }



    thermodynamics

    {

        Cp              5190; //2700; //J/(kg*K)

        Hf              0;

    }

    transport

    {

        mu              1.9e-05; //1.8e-05

        Pr              0.7;

    }

}

and following error came up

Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState PengRobinsonGas;
specie specie;
energy sensibleEnthalpy;
}

--> FOAM FATAL IO ERROR: (openfoam-2206)
Unknown rhoReactionThermo type

Valid rhoReactionThermo types :

69
(
heRhoThermo<homogeneousMixture<const<hConst<incomp ressiblePerfectGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<homogeneousMixture<const<hConst<perfec tGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<homogeneousMixture<sutherland<janaf<in compressiblePerfectGas<specie>>,sensibleEnthalpy>> >
heRhoThermo<homogeneousMixture<sutherland<janaf<pe rfectGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<inhomogeneousMixture<const<hConst<inco mpressiblePerfectGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<inhomogeneousMixture<const<hConst<perf ectGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<inhomogeneousMixture<sutherland<janaf< incompressiblePerfectGas<specie>>,sensibleEnthalpy >>>
heRhoThermo<inhomogeneousMixture<sutherland<janaf< perfectGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<multiComponentMixture<const<eConst<adi abaticPerfectFluid<specie>>,sensibleInternalEnergy >>>
heRhoThermo<multiComponentMixture<const<eConst<inc ompressiblePerfectGas<specie>>,sensibleInternalEne rgy>>>
heRhoThermo<multiComponentMixture<const<eConst<per fectFluid<specie>>,sensibleInternalEnergy>>>
heRhoThermo<multiComponentMixture<const<eConst<per fectGas<specie>>,sensibleInternalEnergy>>>
heRhoThermo<multiComponentMixture<const<eConst<rho Const<specie>>,sensibleInternalEnergy>>>
heRhoThermo<multiComponentMixture<const<hConst<adi abaticPerfectFluid<specie>>,sensibleEnthalpy>>>
heRhoThermo<multiComponentMixture<const<hConst<inc ompressiblePerfectGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<multiComponentMixture<const<hConst<per fectFluid<specie>>,sensibleEnthalpy>>>
heRhoThermo<multiComponentMixture<const<hConst<per fectGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<multiComponentMixture<const<hConst<rho Const<specie>>,sensibleEnthalpy>>>
heRhoThermo<multiComponentMixture<polynomial<hPoly nomial<icoPolynomial<specie>>,sensibleEnthalpy>>>
heRhoThermo<multiComponentMixture<polynomial<hPoly nomial<icoPolynomial<specie>>,sensibleInternalEner gy>>>
heRhoThermo<multiComponentMixture<polynomial<janaf <perfectGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<multiComponentMixture<polynomial<janaf <perfectGas<specie>>,sensibleInternalEnergy>>>
heRhoThermo<multiComponentMixture<sutherland<janaf <incompressiblePerfectGas<specie>>,sensibleEnthalp y>>>
heRhoThermo<multiComponentMixture<sutherland<janaf <incompressiblePerfectGas<specie>>,sensibleInterna lEnergy>>>
heRhoThermo<multiComponentMixture<sutherland<janaf <perfectGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<multiComponentMixture<sutherland<janaf <perfectGas<specie>>,sensibleInternalEnergy>>>
heRhoThermo<pureMixture<WLF<eConst<rhoConst<specie >>,sensibleInternalEnergy>>>
heRhoThermo<pureMixture<const<eConst<adiabaticPerf ectFluid<specie>>,sensibleInternalEnergy>>>
heRhoThermo<pureMixture<const<eConst<incompressibl ePerfectGas<specie>>,sensibleInternalEnergy>>>
heRhoThermo<pureMixture<const<eConst<perfectFluid< specie>>,sensibleInternalEnergy>>>
heRhoThermo<pureMixture<const<eConst<perfectGas<sp ecie>>,sensibleInternalEnergy>>>
heRhoThermo<pureMixture<const<eConst<rhoConst<spec ie>>,sensibleInternalEnergy>>>
heRhoThermo<pureMixture<const<hConst<Boussinesq<sp ecie>>,sensibleEnthalpy>>>
heRhoThermo<pureMixture<const<hConst<Boussinesq<sp ecie>>,sensibleInternalEnergy>>>
heRhoThermo<pureMixture<const<hConst<adiabaticPerf ectFluid<specie>>,sensibleEnthalpy>>>
heRhoThermo<pureMixture<const<hConst<incompressibl ePerfectGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<pureMixture<const<hConst<perfectFluid< specie>>,sensibleEnthalpy>>>
heRhoThermo<pureMixture<const<hConst<perfectGas<sp ecie>>,sensibleEnthalpy>>>
heRhoThermo<pureMixture<const<hConst<rhoConst<spec ie>>,sensibleEnthalpy>>>
heRhoThermo<pureMixture<polynomial<hPolynomial<ico Polynomial<specie>>,sensibleEnthalpy>>>
heRhoThermo<pureMixture<polynomial<hPolynomial<ico Polynomial<specie>>,sensibleInternalEnergy>>>
heRhoThermo<pureMixture<sutherland<janaf<Boussines q<specie>>,sensibleEnthalpy>>>
heRhoThermo<pureMixture<sutherland<janaf<Boussines q<specie>>,sensibleInternalEnergy>>>
heRhoThermo<pureMixture<sutherland<janaf<incompres siblePerfectGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<pureMixture<sutherland<janaf<incompres siblePerfectGas<specie>>,sensibleInternalEnergy>>>
heRhoThermo<pureMixture<sutherland<janaf<perfectGa s<specie>>,sensibleEnthalpy>>>
heRhoThermo<pureMixture<sutherland<janaf<perfectGa s<specie>>,sensibleInternalEnergy>>>
heRhoThermo<reactingMixture<const<eConst<adiabatic PerfectFluid<specie>>,sensibleInternalEnergy>>>
heRhoThermo<reactingMixture<const<eConst<incompres siblePerfectGas<specie>>,sensibleInternalEnergy>>>
heRhoThermo<reactingMixture<const<eConst<perfectFl uid<specie>>,sensibleInternalEnergy>>>
heRhoThermo<reactingMixture<const<eConst<perfectGa s<specie>>,sensibleInternalEnergy>>>
heRhoThermo<reactingMixture<const<eConst<rhoConst< specie>>,sensibleInternalEnergy>>>
heRhoThermo<reactingMixture<const<hConst<adiabatic PerfectFluid<specie>>,sensibleEnthalpy>>>
heRhoThermo<reactingMixture<const<hConst<incompres siblePerfectGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<reactingMixture<const<hConst<perfectFl uid<specie>>,sensibleEnthalpy>>>
heRhoThermo<reactingMixture<const<hConst<perfectGa s<specie>>,sensibleEnthalpy>>>
heRhoThermo<reactingMixture<const<hConst<rhoConst< specie>>,sensibleEnthalpy>>>
heRhoThermo<reactingMixture<polynomial<hPolynomial <icoPolynomial<specie>>,sensibleEnthalpy>>>
heRhoThermo<reactingMixture<polynomial<hPolynomial <icoPolynomial<specie>>,sensibleInternalEnergy>> >
heRhoThermo<reactingMixture<sutherland<janaf<incom pressiblePerfectGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<reactingMixture<sutherland<janaf<incom pressiblePerfectGas<specie>>,sensibleInternalEnerg y>>>
heRhoThermo<reactingMixture<sutherland<janaf<perfe ctGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<reactingMixture<sutherland<janaf<perfe ctGas<specie>>,sensibleInternalEnergy>>>
heRhoThermo<singleStepReactingMixture<sutherland<j anaf<perfectGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<singleStepReactingMixture<sutherland<j anaf<perfectGas<specie>>,sensibleInternalEnergy>>>
heRhoThermo<veryInhomogeneousMixture<const<hConst< incompressiblePerfectGas<specie>>,sensibleEnthalpy >>>
heRhoThermo<veryInhomogeneousMixture<const<hConst< perfectGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<veryInhomogeneousMixture<sutherland<ja naf<incompressiblePerfectGas<specie>>,sensibleEnth alpy>>>
heRhoThermo<veryInhomogeneousMixture<sutherland<ja naf<perfectGas<specie>>,sensibleEnthalpy>>>
)

type mixture transport thermo equationOfState specie energy

heRhoThermo homogeneousMixture const hConst incompressiblePerfectGas specie sensibleEnthalpy
heRhoThermo homogeneousMixture const hConst perfectGas specie sensibleEnthalpy
heRhoThermo homogeneousMixture sutherland janaf incompressiblePerfectGas specie sensibleEnthalpy
heRhoThermo homogeneousMixture sutherland janaf perfectGas specie sensibleEnthalpy
heRhoThermo inhomogeneousMixture const hConst incompressiblePerfectGas specie sensibleEnthalpy
heRhoThermo inhomogeneousMixture const hConst perfectGas specie sensibleEnthalpy
heRhoThermo inhomogeneousMixture sutherland janaf incompressiblePerfectGas specie sensibleEnthalpy
heRhoThermo inhomogeneousMixture sutherland janaf perfectGas specie sensibleEnthalpy
heRhoThermo multiComponentMixture const eConst adiabaticPerfectFluid specie sensibleInternalEnergy
heRhoThermo multiComponentMixture const eConst incompressiblePerfectGas specie sensibleInternalEnergy
heRhoThermo multiComponentMixture const eConst perfectFluid specie sensibleInternalEnergy
heRhoThermo multiComponentMixture const eConst perfectGas specie sensibleInternalEnergy
heRhoThermo multiComponentMixture const eConst rhoConst specie sensibleInternalEnergy
heRhoThermo multiComponentMixture const hConst adiabaticPerfectFluid specie sensibleEnthalpy
heRhoThermo multiComponentMixture const hConst incompressiblePerfectGas specie sensibleEnthalpy
heRhoThermo multiComponentMixture const hConst perfectFluid specie sensibleEnthalpy
heRhoThermo multiComponentMixture const hConst perfectGas specie sensibleEnthalpy
heRhoThermo multiComponentMixture const hConst rhoConst specie sensibleEnthalpy
heRhoThermo multiComponentMixture polynomial hPolynomial icoPolynomial specie sensibleEnthalpy
heRhoThermo multiComponentMixture polynomial hPolynomial icoPolynomial specie sensibleInternalEnergy
heRhoThermo multiComponentMixture polynomial janaf perfectGas specie sensibleEnthalpy
heRhoThermo multiComponentMixture polynomial janaf perfectGas specie sensibleInternalEnergy
heRhoThermo multiComponentMixture sutherland janaf incompressiblePerfectGas specie sensibleEnthalpy
heRhoThermo multiComponentMixture sutherland janaf incompressiblePerfectGas specie sensibleInternalEnergy
heRhoThermo multiComponentMixture sutherland janaf perfectGas specie sensibleEnthalpy
heRhoThermo multiComponentMixture sutherland janaf perfectGas specie sensibleInternalEnergy
heRhoThermo pureMixture WLF eConst rhoConst specie sensibleInternalEnergy
heRhoThermo pureMixture const eConst adiabaticPerfectFluid specie sensibleInternalEnergy
heRhoThermo pureMixture const eConst incompressiblePerfectGas specie sensibleInternalEnergy
heRhoThermo pureMixture const eConst perfectFluid specie sensibleInternalEnergy
heRhoThermo pureMixture const eConst perfectGas specie sensibleInternalEnergy
heRhoThermo pureMixture const eConst rhoConst specie sensibleInternalEnergy
heRhoThermo pureMixture const hConst Boussinesq specie sensibleEnthalpy
heRhoThermo pureMixture const hConst Boussinesq specie sensibleInternalEnergy
heRhoThermo pureMixture const hConst adiabaticPerfectFluid specie sensibleEnthalpy
heRhoThermo pureMixture const hConst incompressiblePerfectGas specie sensibleEnthalpy
heRhoThermo pureMixture const hConst perfectFluid specie sensibleEnthalpy
heRhoThermo pureMixture const hConst perfectGas specie sensibleEnthalpy
heRhoThermo pureMixture const hConst rhoConst specie sensibleEnthalpy
heRhoThermo pureMixture polynomial hPolynomial icoPolynomial specie sensibleEnthalpy
heRhoThermo pureMixture polynomial hPolynomial icoPolynomial specie sensibleInternalEnergy
heRhoThermo pureMixture sutherland janaf Boussinesq specie sensibleEnthalpy
heRhoThermo pureMixture sutherland janaf Boussinesq specie sensibleInternalEnergy
heRhoThermo pureMixture sutherland janaf incompressiblePerfectGas specie sensibleEnthalpy
heRhoThermo pureMixture sutherland janaf incompressiblePerfectGas specie sensibleInternalEnergy
heRhoThermo pureMixture sutherland janaf perfectGas specie sensibleEnthalpy
heRhoThermo pureMixture sutherland janaf perfectGas specie sensibleInternalEnergy
heRhoThermo reactingMixture const eConst adiabaticPerfectFluid specie sensibleInternalEnergy
heRhoThermo reactingMixture const eConst incompressiblePerfectGas specie sensibleInternalEnergy
heRhoThermo reactingMixture const eConst perfectFluid specie sensibleInternalEnergy
heRhoThermo reactingMixture const eConst perfectGas specie sensibleInternalEnergy
heRhoThermo reactingMixture const eConst rhoConst specie sensibleInternalEnergy
heRhoThermo reactingMixture const hConst adiabaticPerfectFluid specie sensibleEnthalpy
heRhoThermo reactingMixture const hConst incompressiblePerfectGas specie sensibleEnthalpy
heRhoThermo reactingMixture const hConst perfectFluid specie sensibleEnthalpy
heRhoThermo reactingMixture const hConst perfectGas specie sensibleEnthalpy
heRhoThermo reactingMixture const hConst rhoConst specie sensibleEnthalpy
heRhoThermo reactingMixture polynomial hPolynomial icoPolynomial specie sensibleEnthalpy
heRhoThermo reactingMixture polynomial hPolynomial icoPolynomial specie sensibleInternalEnergy
heRhoThermo reactingMixture sutherland janaf incompressiblePerfectGas specie sensibleEnthalpy
heRhoThermo reactingMixture sutherland janaf incompressiblePerfectGas specie sensibleInternalEnergy
heRhoThermo reactingMixture sutherland janaf perfectGas specie sensibleEnthalpy
heRhoThermo reactingMixture sutherland janaf perfectGas specie sensibleInternalEnergy
heRhoThermo singleStepReactingMixture sutherland janaf perfectGas specie sensibleEnthalpy
heRhoThermo singleStepReactingMixture sutherland janaf perfectGas specie sensibleInternalEnergy
heRhoThermo veryInhomogeneousMixture const hConst incompressiblePerfectGas specie sensibleEnthalpy
heRhoThermo veryInhomogeneousMixture const hConst perfectGas specie sensibleEnthalpy
heRhoThermo veryInhomogeneousMixture sutherland janaf incompressiblePerfectGas specie sensibleEnthalpy
heRhoThermo veryInhomogeneousMixture sutherland janaf perfectGas specie sensibleEnthalpy

file: constant/midAir/thermophysicalProperties.thermoType at line 19 to 25.

From static typename ThermoConstructTable::mapped_type Foam::basicThermo::getThermoOrDie(const Foam::dictionary&, ThermoConstructTable&, const Foam::word&, const wordList&) [with Thermo = Foam::rhoReactionThermo; ThermoConstructTable = Foam::HashTable<Foam::autoPtr<Foam::rhoReactionThe rmo> (*)(const Foam::fvMesh&, const Foam::word&), Foam::word, Foam::Hash<Foam::word> >; typename ThermoConstructTable::mapped_type = Foam::autoPtr<Foam::rhoReactionThermo> (*)(const Foam::fvMesh&, const Foam::word&); Foam::wordList = Foam::List<Foam::word>]
in file ./src/thermophysicalModels/basic/lnInclude/basicThermoTemplates.C at line 50.

FOAM exiting

I would appreciate your help, thank you
Dasith

dasith0001 · September 21, 2022, 02:57

Update;

I tried to update src/thermophysicalModel/basic/rhoThermo/rhoThermos.C

as

PHP Code:


			
makeThermos

(

    rhoThermo,

    heRhoThermo,

    pureMixture,

    constTransport,

    sensibleEnthalpy,

    hConstThermo,

    PengRobinsonGas,

    specie

);

But still the error is

--> FOAM FATAL IO ERROR: (openfoam-2206)
Unknown rhoReactionThermo type

Valid rhoReactionThermo types :

70
(
heRhoThermo<homogeneousMixture<const<hConst<PengRo binsonGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<homogeneousMixture<const<hConst<incomp ressiblePerfectGas<specie>>,sensibleEnthalpy>>>
heRhoThermo<homo..............

Thank you for any help.
Dasith

Yann · September 21, 2022, 04:56

Hi Dasith,

OK I understand, this error usually happens when trying to use an invalid thermo package: it's not possible to use any combination of models (thermo type, transport, thermo, equationOfState, etc)

After a quick search I found this thread: PengRobinsonGas crashing with high pressures on fine mesh
where Peng-Robinson is used with this thermo package:

Code:

thermoType
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       polynomial;
    thermo          hPolynomial;
    equationOfState PengRobinsonGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

So I guess you have 2 options: using a suitable thermo package to be able to use Peng-Robinson, which will require to adapt your case setup (no idea if this will be compatible with chtMultiRegionFoam) or compile your own thermo package to be able to use it as you wish. But it might require to adjust a bunch of stuff in the code and you need to be sure your changes make sense, physically speaking.

Yann

dasith0001 · September 28, 2022, 08:12

Hi All,

Thank you for all the ideas, finally figure it out.

OpenFoam Version 10 facilitates PengRobinsonGas model with chtMultiRegionFoam solver, my model works fine now.

Cheers,
Dasith

Yann · September 28, 2022, 09:26

Hi Dasith,

Thanks for your feedback.
Would you mind sharing your thermophysicalProperties setup? It could help other people who want to use PengRobinsonGas model.

Thanks,
Yann

dasith0001 · September 30, 2022, 05:11

Hi,

Not a problem, this is what I used. no guarantee for the values though.

PHP Code:


			
/*--------------------------------*- C++ -*----------------------------------*\

  =========                 |

  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox

   \\    /   O peration     | Website:  https://openfoam.org

    \\  /    A nd           | Version:  10

     \\/     M anipulation  |

\*---------------------------------------------------------------------------*/

FoamFile

{

    format      ascii;

    class       dictionary;

    location    "constant";

    object      physicalProperties;

}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //



thermoType

{

    type            heRhoThermo;     //hePsiThermo; //heRhoThermo

    mixture         pureMixture;     //homogeneousMixture pureMixture

    transport       const; 

    thermo          hConst;

    equationOfState PengRobinsonGas;

    specie          specie;

    energy          sensibleEnthalpy;

}



mixture

{

    specie

    {

        molWeight       39.948; //g/mol

    }



    equationOfState

    {

        Tc 151.15;       //- Critical Temperature [K]

        Vc 0.07430328;   //- Critical volume [m^3/kmol]

        Pc 4.875e+6;     //- Critical Pressure [Pa]

        #scalar Zc_;     //- Critical compression factor [-]

        omega 0.001;     //- Acentric factor [-]

    

    }



    thermodynamics

    {

        Cp              580;  //J/(kg*K) 

        Hf              0;

    }

    transport

    {

        mu              5.05E-05;

        Pr              0.65951;

    }

}



// ************************************************************************* //

Thanks,
Dasith

October 9, 2019, 04:23	Conjugate heat transfer with real gas EOS	#1
aalzoubi06 New Member Ahmad Al-Zoubi Join Date: Nov 2013 Posts: 14 Rep Power: 11	Hello All, I am wondering if there is a possibility to run conjugate heat transfer such as chtMultiRegionFoam using a real gas Equation of State such as Peng–Robinson. Many thanks in advance.

September 21, 2022, 02:57		#7
dasith0001 Member Dasith Desh Join Date: Mar 2021 Location: Sydney Posts: 34 Rep Power: 3	Update; I tried to update src/thermophysicalModel/basic/rhoThermo/rhoThermos.C as PHP Code: `makeThermos ( rhoThermo, heRhoThermo, pureMixture, constTransport, sensibleEnthalpy, hConstThermo, PengRobinsonGas, specie );` But still the error is --> FOAM FATAL IO ERROR: (openfoam-2206) Unknown rhoReactionThermo type Valid rhoReactionThermo types : 70 ( heRhoThermo<homogeneousMixture<const<hConst<PengRo binsonGas<specie>>,sensibleEnthalpy>>> heRhoThermo<homogeneousMixture<const<hConst<incomp ressiblePerfectGas<specie>>,sensibleEnthalpy>>> heRhoThermo<homo.............. Thank you for any help. Dasith

September 21, 2022, 04:56		#8
Yann Senior Member Yann Join Date: Apr 2012 Location: France Posts: 572 Rep Power: 19	Hi Dasith, OK I understand, this error usually happens when trying to use an invalid thermo package: it's not possible to use any combination of models (thermo type, transport, thermo, equationOfState, etc) After a quick search I found this thread: PengRobinsonGas crashing with high pressures on fine mesh where Peng-Robinson is used with this thermo package: Code: thermoType { type hePsiThermo; mixture pureMixture; transport polynomial; thermo hPolynomial; equationOfState PengRobinsonGas; specie specie; energy sensibleEnthalpy; } So I guess you have 2 options: using a suitable thermo package to be able to use Peng-Robinson, which will require to adapt your case setup (no idea if this will be compatible with chtMultiRegionFoam) or compile your own thermo package to be able to use it as you wish. But it might require to adjust a bunch of stuff in the code and you need to be sure your changes make sense, physically speaking. Yann

September 30, 2022, 05:11		#11
dasith0001 Member Dasith Desh Join Date: Mar 2021 Location: Sydney Posts: 34 Rep Power: 3	Hi, Not a problem, this is what I used. no guarantee for the values though. PHP Code: /--------------------------------- C++ -----------------------------------\ ========= \| \\ / F ield \| OpenFOAM: The Open Source CFD Toolbox \\ / O peration \| Website: https://openfoam.org \\ / A nd \| Version: 10 \\/ M anipulation \| \---------------------------------------------------------------------------/ FoamFile { format ascii; class dictionary; location "constant"; object physicalProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // thermoType { type heRhoThermo; //hePsiThermo; //heRhoThermo mixture pureMixture; //homogeneousMixture pureMixture transport const; thermo hConst; equationOfState PengRobinsonGas; specie specie; energy sensibleEnthalpy; } mixture { specie { molWeight 39.948; //g/mol } equationOfState { Tc 151.15; //- Critical Temperature [K] Vc 0.07430328; //- Critical volume [m^3/kmol] Pc 4.875e+6; //- Critical Pressure [Pa] #scalar Zc_; //- Critical compression factor [-] omega 0.001; //- Acentric factor [-] } thermodynamics { Cp 580; //J/(kgK) Hf 0; } transport { mu 5.05E-05; Pr 0.65951; } } // ************************************************************************ // Thanks, Dasith Yann likes this.

Thread Tools	Search this Thread
Show Printable Version Email this Page	Search this Thread: Advanced Search
Display Modes
Linear Mode Switch to Hybrid Mode Switch to Threaded Mode

August 18, 2022, 03:58		#2
dasith0001 Member Dasith Desh Join Date: Mar 2021 Location: Sydney Posts: 34 Rep Power: 3	Very keen to know if this is possible at all, adapting Peng-Robinson into chtmultiRegionFoam

August 25, 2022, 12:08		#3
jherb Senior Member Joachim Herb Join Date: Sep 2010 Posts: 648 Rep Power: 20	Yes, I used Peng-Robinson in a simulation with chtMultiRegionFoam and it worked.

September 20, 2022, 04:26		#5
Yann Senior Member Yann Join Date: Apr 2012 Location: France Posts: 572 Rep Power: 19	Hi, As long as Peng-Robinson is implemented in the thermophysical models, it should be possible to use it with chtMultiRegionFoam : https://www.openfoam.com/documentati...onOfState.html Searching the forum might also help to find answers: Usage of PengRobinsonGas I hope this helps, Yann

September 28, 2022, 08:12		#9
dasith0001 Member Dasith Desh Join Date: Mar 2021 Location: Sydney Posts: 34 Rep Power: 3	Hi All, Thank you for all the ideas, finally figure it out. OpenFoam Version 10 facilitates PengRobinsonGas model with chtMultiRegionFoam solver, my model works fine now. Cheers, Dasith

September 28, 2022, 09:26		#10
Yann Senior Member Yann Join Date: Apr 2012 Location: France Posts: 572 Rep Power: 19	Hi Dasith, Thanks for your feedback. Would you mind sharing your thermophysicalProperties setup? It could help other people who want to use PengRobinsonGas model. Thanks, Yann

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
UDF for thermal conductivity relation of nano fluid	ngonbadipoor@yahoo.com	Fluent UDF and Scheme Programming	29	April 1, 2017 04:02
Erron: FLUENT received fatal signal (ACCESS_VIOLATION) in DPM Heat and Mass Transfer	Melvin	FLUENT	2	November 17, 2016 02:06
Error - Solar absorber - Solar Thermal Radiation	MichaelK	CFX	12	September 1, 2016 06:15
Heat Flux at wall in a conjugate heat transfer problem	Chander	CFX	2	July 9, 2011 23:22
Modelling the heat transfer during compression and cooling of natural gas	pano	Main CFD Forum	0	December 10, 2010 16:53