[brbbhatti]

Hello Everyone,

I am using Liquid Evaporation Model for Evaporation of Urea and Water. I have worked on OF 2.1.0 version and have problem with the Exhaust temperature drop and lot of colleagues here on forum ask me to use OF 2.3.0.

So now i have also same or even more problems and errors in OpenFoam 2.3.0.

1. When i use Liquid Evaporation Model for H20 and SingleKineticRateDevolitilization Model for Urea together the output results are good but they are only running in serial when i try to run it in parallel it crash. I am using simpleReactingParcelFoam solver.
Below you find the error i recieved

[4] #0 Foam::error:rintStack(Foam::Ostream&) in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[4] #1 Foam::sigFpe::sigHandler(int) in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[4] #2 in "/lib64/libc.so.6"
[4] #3 Foam::tmp<Foam::Field<double> > Foam:perator/<double>(Foam::UList<double> const&, double const&) in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[4] #4 Foam::liquidMixtureProperties::X(Foam::Field<doubl e> const&) const in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libliquidMixtureProperties.so"
[4] #5 Foam::LiquidEvaporation<Foam::ReactingCloud<Foam:: ThermoCloud<Foam::KinematicCloud<Foam::Cloud<Foam: :ReactingMultiphaseParcel<Foam::ReactingParcel<Foa m::ThermoParcel<Foam::KinematicParcel<Foam:artic le> > > > > > > > >::Tvap(Foam::Field<double> const&) const in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/liblagrangianIntermediate.so"
[4] #6
[4] in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/simpleReactingParcelFoam"
[4] #7
[4] in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/simpleReactingParcelFoam"
[4] #8
[4] in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/simpleReactingParcelFoam"
[4] #9
[4] in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/simpleReactingParcelFoam"
[4] #10
[4] at simpleReactingParcelFoam.C:0
[4] #11
[4] in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/simpleReactingParcelFoam"
[4] #12 __libc_start_main in "/lib64/libc.so.6"
[4] #13
[4] at /home/abuild/rpmbuild/BUILD/glibc-2.14.1/csu/../sysdeps/x86_64/elf/start.S:116
[fe-z0as9:01363] *** Process received signal ***
[fe-z0as9:01363] Signal: Floating point exception (8)
[fe-z0as9:01363] Signal code: (-6)
[fe-z0as9:01363] Failing at address: 0x12a6800000553
[fe-z0as9:01363] [ 0] /lib64/libc.so.6(+0x34e10) [0x2afbeb957e10]
[fe-z0as9:01363] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x2afbeb957d95]
[fe-z0as9:01363] [ 2] /lib64/libc.so.6(+0x34e10) [0x2afbeb957e10]
[fe-z0as9:01363] [ 3] /usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4FoamdvIdEENS_3tmpINS_5Field IT_EEEERKNS_5UListIS3_EERKd+0x65) [0x2afbe4350085]
[fe-z0as9:01363] [ 4] /usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libliquidMixtureProperties.so(_ZNK4Foam23liquidMix tureProperties1XERKNS_5FieldIdEE+0x18c) [0x2afbe797b2bc]
[fe-z0as9:01363] [ 5] /usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/liblagrangianIntermediate.so(_ZNK4Foam17LiquidEvap orationINS_13ReactingCloudINS_11ThermoCloudINS_14K inematicCloudINS_5CloudINS_24ReactingMultiphasePar celINS_14ReactingParcelINS_12ThermoParcelINS_15Kin ematicParcelINS_8particleEEEEEEEEEEEEEEEEEE4TvapER KNS_5FieldIdEE+0x1c) [0x2afbe677926c]
[fe-z0as9:01363] [ 6] simpleReactingParcelFoam() [0x479d7d]
[fe-z0as9:01363] [ 7] simpleReactingParcelFoam() [0x487207]
[fe-z0as9:01363] [ 8] simpleReactingParcelFoam() [0x4a02b8]
[fe-z0as9:01363] [ 9] simpleReactingParcelFoam() [0x4acfa8]
[fe-z0as9:01363] [10] simpleReactingParcelFoam() [0x442824]
[fe-z0as9:01363] [11] simpleReactingParcelFoam() [0x42b4b6]
[fe-z0as9:01363] [12] /lib64/libc.so.6(__libc_start_main+0xed) [0x2afbeb94423d]
[fe-z0as9:01363] [13] simpleReactingParcelFoam() [0x431c01]
[fe-z0as9:01363] *** End of error message ***

.

The error shows that there is a problem with the liquidMixtureProperties but i did not find the solution of that so can you explain me what is wrong here why its not running??

2. Second problem when i use only Liquid Evaporation Model for evaporation of both water and urea it also crashes and gives error just after 1 simulation. it is giving the error showing that there is a problem with the calculation of the Temperature as the temperature is not in the certain range. The error is given below

FOAM Warning :
From function janafThermo<EquationOfState>::limit(const scalar T) const
in file /usr/local/openfoam/OpenFOAM-2.3.0/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 108
attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = -2.62401293e+114
--> FOAM Warning :
From function janafThermo<EquationOfState>::limit(const scalar T) const
in file /usr/local/openfoam/OpenFOAM-2.3.0/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 108
attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = -2.85584986e+114

But when i run the same simulation in OF 2.1.0 it is running OK! but in OF 2.1.0 i am using LTSReactingParcelFoam.

I didn't understand why this is not working properly in OF 2.3.0 ??

3. Third problem When i run the same simulation of liquid evaporation model only taking the water then simulation run for 10 steps and then crashes. The error in log file is different this time now it is related to pressure. The error shows that my pressure is not in the desired range it is exceeding the triple point limit and i did not understand why is it so.

on other hand the same simulation is running in OF 2.1.0 without any error but not in OF 2.3.0. The error details are given below

FOAM Warning :
From function Foam::scalar Foam::liquidMixtureProperties:vInvert( const scalar, const scalarField&) const
in file liquidMixtureProperties/liquidMixtureProperties.C at line 142
Pressure below triple point pressure: p = 115313.434 < Pt = 14554275

and on terminal it show error

#0 Foam::error:rintStack(Foam::Ostream&) in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib64/libc.so.6"
#3 pow in "/lib64/libm.so.6"
#4 Foam::CH4N2O:(double, double) const in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libliquidProperties.so"
#5 Foam::LiquidEvaporation<Foam::ReactingCloud<Foam:: ThermoCloud<Foam::KinematicCloud<Foam::Cloud<Foam: :ReactingMultiphaseParcel<Foam::ReactingParcel<Foa m::ThermoParcel<Foam::KinematicParcel<Foam:artic le> > > > > > > > >::calculate(double, int, double, double, double, double, double, double, double, double, Foam::Field<double> const&, Foam::Field<double>&) const in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/liblagrangianIntermediate.so"
#6
in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/simplePorousReactingParcelFoam"
#7
in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/simplePorousReactingParcelFoam"
#8
in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/simplePorousReactingParcelFoam"
#9
in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/simplePorousReactingParcelFoam"
#10
at simplePorousReactingParcelFoam.C:0
#11
at simplePorousReactingParcelFoam.C:0
#12
in "/usr/local/openfoam/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/bin/simplePorousReactingParcelFoam"
#13 __libc_start_main in "/lib64/libc.so.6"
#14
at /home/abuild/rpmbuild/BUILD/glibc-2.14.1/csu/../sysdeps/x86_64/elf/start.S:116
Floating point exception

So what should i do because OF 2.3.0 is giving more errors and problem.
And OF 2.1.0 is also giving the some problem with Temperature drop calculation.
Temperature drop: means that i am analytically making the analysis of the change in the temperature at outlet. I calculate the amount of energy by taking in account the mass of species evaporate during simulation and using some theoretical Formula. Then calculate the temperature at outlet and compare it with the results i got from simulation but the results are not same.

If you need some more details of model i am using please check the link below where i posted early all details
http://www.cfd-online.com/Forums/ope...sel-spray.html

Thanks a lot

Regards,
Bilal

[dkxls]

Which evaporation model are you using?
I hope 'liquidEvaporationBoil', the other one is not suited for spray evaporation. If you search the forum, you will find more explanation on this by Tommaso Lucchini.

[brbbhatti]

Originally Posted by dkxls

Which evaporation model are you using?
I hope 'liquidEvaporationBoil', the other one is not suited for spray evaporation. If you search the forum, you will find more explanation on this by Tommaso Lucchini.

Hi Armin,

Thanks. i am using Liquid Evaporation Model.
What you mean is that LiquidEvaporationBoil Model i have to use instead of the other one right?
I check the details you mention and Tommaso Lucchini just mention that there is a problem with the boiling conditions nothing more.

And one more thing if i change the Model what you suggest do i run it in OF 2.1.0 or OF 2.3.0?

Best Regards,

Bilal

[dkxls]

Originally Posted by brbbhatti

Thanks. i am using Liquid Evaporation Model.

Which one?

Originally Posted by brbbhatti

What you mean is that LiquidEvaporationBoil Model i have to use instead of the other one right?

You should use this one for spray simulation:

    phaseChangeModel liquidEvaporationBoil;

Originally Posted by brbbhatti

And one more thing if i change the Model what you suggest do i run it in OF 2.1.0 or OF 2.3.0?

2.3.x use a git version, there are bug fixes in already and more to come...

Also, for spray simulation I would recommend you the 'sprayFoam' solver.

[brbbhatti]

Hi Armin,

Thanks once again.

I am using

PhaseChangeModel liquidEvaporation;

But as you suggest i will run the simulation on liquidEvaporationBoil Model and then i will let you know what i got.

But one thing you suggest me to use 2.3.x i will run it again with new changes to see temperature behavior but as i mentioned there are lot of other failures coming on especially i can not run it parallel .....anything on that

[brbbhatti]

Originally Posted by dkxls

Which one?


You should use this one for spray simulation:

    phaseChangeModel liquidEvaporationBoil;

2.3.x use a git version, there are bug fixes in already and more to come...

Also, for spray simulation I would recommend you the 'sprayFoam' solver.

Hi Armin,

I run the simulation with changes as you suggest on OF 2.3.0

phaseChangeModel liquidEvaporationBoil;

But i got the same error message when i use liquidEvaporationBoil for mixture

FOAM Warning :
From function Foam::scalar Foam::liquidMixtureProperties:vInvert( const scalar, const scalarField&) const
in file liquidMixtureProperties/liquidMixtureProperties.C at line 142
Pressure below triple point pressure: p = 115313.434 < Pt = 14554275

And

FOAM Warning :
From function janafThermo<EquationOfState>::limit(const scalar T) const
in file /usr/local/openfoam/OpenFOAM-2.3.0/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 108
attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = -2.62401293e+114
--> FOAM Warning :
From function janafThermo<EquationOfState>::limit(const scalar T) const
in file /usr/local/openfoam/OpenFOAM-2.3.0/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 108
attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = -2.85584986e+114

So any suggestion now what is the problem??

Best Regards,
Bilal

[dkxls]

Sorry, cannot help you with this one.

There is most likely something wrong in your case setup, which leads to these wrong temperatures.

If you are sure that all your boundary and initial conditions are right, I would probably check on the time step and on the maxCo number for the droplets.

-Armin

[brbbhatti]

Hi Armin,

I did some changes with the time step and track time in OF 2.1.0. As i increase the Time step (DeltaT) it give the desired results but only when i use both Liquid Evaporation and SingleKineticRateDevolitilization model. As i use only LiquidEvaporation Model the temperature loss at exhaust is more and i did not understand what's the reason.
When i run the evaporation model with water only in OF 2.1.0 the temperature value is ok but when i use mixture of water and urea then temperature is not right. Maybe there is some problem with the liquid properties of Urea or Enthalpy calculation by using NSRDS functions.

Thanks

Best Regards,
Bilal

[Polli]

Hi!

I have exactly the same problem as mentioned in the first part.

"When i use Liquid Evaporation Model for H20 and SingleKineticRateDevolitilization Model for Urea together the output results are good but they are only running in serial when i try to run it in parallel it crash. I am using simpleReactingParcelFoam solver.
Below you find the error i recieved"

do you already fixed the problem with the breakdown in the parallel simulation?

To your problem:

FOAM Warning :
From function janafThermo<EquationOfState>::limit(const scalar T) const
in file /usr/local/openfoam/OpenFOAM-2.3.0/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 108
attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = -2.62401293e+114
--> FOAM Warning :
From function janafThermo<EquationOfState>::limit(const scalar T) const
in file /usr/local/openfoam/OpenFOAM-2.3.0/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 108
attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = -2.85584986e+114

You have to reduce the timestep, dont take the adjusted timestep with the corant number. I had the same problem. If this dont work than reduce the ralaxation factors from 1 to 0.85 than it works. Mabe this have also to do with radiation. I got this problem several times but now it runs perfect.

[brbbhatti]

Hi Werner,

No, still i can not run that in OF 2.3.0 still having same problem.
So its better you use OF 2.1.0 it is working ok in series and parallel.

And i tried to reduce the time step but results are also not good. its better if you increase a track time you will get better results.

Regards,

Bilal

[dkxls]

Are you guys sure you are running the very latest 2.3.x git-version?

There have been several bug-fixes in the 2.3.x branch which also affected parallel runs (bugs that have also been present in 2.1, but might not show due to all the other things that go wrong there).

See e.g. here:
http://www.openfoam.org/mantisbt/view.php?id=1304
http://www.openfoam.org/mantisbt/view.php?id=1298

[brbbhatti]

Hi Armin,

Thanks for the updates, i will try it and let you know.