[xiaoyoyo]

Hi Everyone,

Have anyone tried the chtMultiRegionFoam solver application for high velocity& High Pressure?
I was using the cht heater tutorial.
when I increased the bottom air velocity (~100m/s) and pressure (~10 atm), I got the error massage after several iterations.

--> FOAM FATAL ERROR:
Maximum number of iterations exceeded

From function thermo<Thermo, Type>::T(scalar f, scalar T0, scalar (thermo<Thermo, Type>::*F)(const scalar) const, scalar (thermo<Thermo, Type>::*dFdT)(const scalar) const, scalar (thermo<Thermo, Type>::*limit)(const scalar) const) const
in file /home/opencfd/OpenFOAM/OpenFOAM-2.2.0/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::species::thermo<Foam::hConstThermo<Foam:er fectGas<Foam::specie> >, Foam::sensibleEnthalpy>::T(double, double, double, double (Foam::species::thermo<Foam::hConstThermo<Foam:e rfectGas<Foam::specie> >, Foam::sensibleEnthalpy>::*)(double, double) const, double (Foam::species::thermo<Foam::hConstThermo<Foam:e rfectGas<Foam::specie> >, Foam::sensibleEnthalpy>::*)(double, double) const, double (Foam::species::thermo<Foam::hConstThermo<Foam:e rfectGas<Foam::specie> >, Foam::sensibleEnthalpy>::*)(double) const) const at ??:?
#3 Foam::heRhoThermo<Foam::rhoThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleEnthalpy> > > >::calculate() at ??:?
#4 Foam::heRhoThermo<Foam::rhoThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleEnthalpy> > > >::correct() at ??:?
#5
at ??:?
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
at ??:?
Aborted (core dumped)

I will appreciate it if you can give me any suggestion or show me the example of your solver and case.

Anyone can help?
thanks a lot!

[laplacian]

Originally Posted by xiaoyoyo

Hi Everyone,
Have anyone tried the chtMultiRegionFoam solver application for high velocity& High Pressure?
I was using the cht heater tutorial.
when I increased the bottom air velocity (~100m/s) and pressure (~10 atm), I got the error massage after several iterations.

I will appreciate it if you can give me any suggestion or show me the example of your solver and case.

Anyone can help?
thanks a lot!

Hi YuanLiu,
your topic is really interesting also for me.
I think one of the problems could be connected with the Mach 0.3 and the formulation of energy equation:
openfoam-solving/116578-big-problem-predicting-adiabatic-wall-temperature-wing-mach-0-33-rhosimplecfoam.html
Which version of OF are you using? (the last release has quite big differences in energy balance)
You could try several options:
1) replace hRhoThermo with hPsiThermo
2) modify the SIMPLE algorithm to SIMPLEC
3) check which h-Eqn is implemented: you could include the term of viscous dissipation

BTW which is your Re number:laminar or turbulent?

[xiaoyoyo]

Originally Posted by laplacian

Hi YuanLiu,
your topic is really interesting also for me.
I think one of the problems could be connected with the Mach 0.3 and the formulation of energy equation:
openfoam-solving/116578-big-problem-predicting-adiabatic-wall-temperature-wing-mach-0-33-rhosimplecfoam.html
Which version of OF are you using? (the last release has quite big differences in energy balance)
You could try several options:
1) replace hRhoThermo with hPsiThermo
2) modify the SIMPLE algorithm to SIMPLEC
3) check which h-Eqn is implemented: you could include the term of viscous dissipation

BTW which is your Re number:laminar or turbulent?

Thanks a lot! laplacian!
I am using the newest version of OF. And I tried both turbulent and laminar case.
I will try as you suggested. Thank you again.