CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Peng Robinson into chtmultiRegionFoam

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 18, 2022, 04:14
Default Peng Robinson into chtmultiRegionFoam
  #1
Senior Member
 
Desh
Join Date: Mar 2021
Location: Sydney
Posts: 118
Rep Power: 5
dasith0001 is on a distinguished road
Hi Foamers,

I would greatly appreciate if one could direct me how to incorporate Pen Robinson gas model into chtmultiRegionFoam solver.

I am currently running on OpenFOAM-v2012 and chtMultiRegionFoam solver does not seem to recognise Pen Robinson gas model, any ideas?

Thank you.
Desh
dasith0001 is offline   Reply With Quote

Old   August 18, 2022, 04:23
Default
  #2
Senior Member
 
Kumaresh
Join Date: Oct 2016
Posts: 355
Rep Power: 12
Kummi is on a distinguished road
Send a message via Yahoo to Kummi
Peng Robinson equation of state ? Have you tried fixing through rho model
Kummi is offline   Reply With Quote

Old   September 20, 2022, 04:15
Default
  #3
Senior Member
 
Desh
Join Date: Mar 2021
Location: Sydney
Posts: 118
Rep Power: 5
dasith0001 is on a distinguished road
Quote:
Originally Posted by Kummi View Post
Peng Robinson equation of state ? Have you tried fixing through rho model
Sorry I don't think I follow you, can you please elaborate how to fix it through rho model?

Thank you.
Dasith
dasith0001 is offline   Reply With Quote

Old   September 20, 2022, 05:43
Default
  #4
Senior Member
 
Kumaresh
Join Date: Oct 2016
Posts: 355
Rep Power: 12
Kummi is on a distinguished road
Send a message via Yahoo to Kummi
I believe you are trying simulate real gas scenario by Peng Robinson model. Have you tried fixing something like this in thermoPhysicalProperties of your model,
Quote:
thermoType
{
* * type * * * * * *hePsiThermo;
* * mixture * * * * pureMixture;
* * transport * * * sutherland;
* * thermo * * * * *hConst;
* * energy * * * * *sensibleEnthalpy;
* * equationOfState PengRobinsonGas;
* * specie * * * * *specie;
}
There is no such information in your previous question either.
There are some threads related to Peng Robinson in CFD forum. Have you double-checked it ?
Usage of PengRobinsonGas

Post your errors please.

Thank you
Kummi is offline   Reply With Quote

Old   September 20, 2022, 21:09
Default
  #5
Senior Member
 
Desh
Join Date: Mar 2021
Location: Sydney
Posts: 118
Rep Power: 5
dasith0001 is on a distinguished road
Hi,

Sorry for inadequate info. my fluid region is under immense pressure and temperature and that's why I am moving away from perfectgas model. And I already followed the link you have mentioned. my thermophysical file for the fluid is as follows.

PHP Code:
{
    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 I am getting this following error when using chtmultiRegionFoam solver


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

Thank you.

Best Regards,
Dasith
dasith0001 is offline   Reply With Quote

Old   September 21, 2022, 02:59
Default
  #6
Senior Member
 
Desh
Join Date: Mar 2021
Location: Sydney
Posts: 118
Rep Power: 5
dasith0001 is on a distinguished road
please find an update in this link

HTML Code:
https://www.cfd-online.com/Forums/openfoam-pre-processing/221218-conjugate-heat-transfer-real-gas-eos.html#post836197
Thank you
dasith0001 is offline   Reply With Quote

Old   September 21, 2022, 03:05
Default
  #7
Senior Member
 
Desh
Join Date: Mar 2021
Location: Sydney
Posts: 118
Rep Power: 5
dasith0001 is on a distinguished road
Quote:
Originally Posted by Kummi View Post
I believe you are trying simulate real gas scenario by Peng Robinson model. Have you tried fixing something like this in thermoPhysicalProperties of your model,

There is no such information in your previous question either.
There are some threads related to Peng Robinson in CFD forum. Have you double-checked it ?
Usage of PengRobinsonGas

Post your errors please.

Thank you
I run the same as you have mentioned and the error is same .

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

Valid rhoReactionThermo types :

70 
thank you
dasith0001 is offline   Reply With Quote

Old   September 22, 2022, 01:28
Default
  #8
Senior Member
 
Kumaresh
Join Date: Oct 2016
Posts: 355
Rep Power: 12
Kummi is on a distinguished road
Send a message via Yahoo to Kummi
Hello Dasith,
The solver you are using doesn't support customized rho model. So, you didn't construct PengRobinsonGas properly.
Let me look into your work by arranging a meeting.
Thank you
Kummi is offline   Reply With Quote

Old   September 28, 2022, 08:09
Default
  #9
Senior Member
 
Desh
Join Date: Mar 2021
Location: Sydney
Posts: 118
Rep Power: 5
dasith0001 is on a distinguished road
Hi Kummi,

Thank you very much for offering the help, but I figure it out.

OpenFoam Version 10 facilitates the PengRobinson model with chtmultiRegionFoam. My model works fine now.

Thanks
Dasith
dasith0001 is offline   Reply With Quote

Reply

Tags
chtmultiregionfoam, peng robinson

« Previous Thread | Next Thread »

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts
BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Setting the height of the stream in the free channel kevinmccartin CFX 12 October 13, 2022 22:43
Domain Reference Pressure and mass flow inlet boundary AdidaKK CFX 75 August 20, 2018 06:37
mass flow in is not equal to mass flow out saii CFX 12 March 19, 2018 06:21
error message cuteapathy CFX 14 March 20, 2012 07:45
Constant velocity of the material Sas CFX 15 July 13, 2010 09:56


All times are GMT -4. The time now is 13:06.

Contact Us - CFD Online - Privacy Statement - Top