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

Liquid Evaporation Model Error in OF 2.3.0

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 14, 2014, 12:31
Question Liquid Evaporation Model Error in OF 2.3.0
  #1
Member
 
bilal
Join Date: Mar 2014
Location: Germany
Posts: 30
Rep Power: 10
brbbhatti is on a distinguished road
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

Quote:
[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

Quote:
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


Quote:
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

Quote:
#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
brbbhatti is offline   Reply With Quote

Old   May 15, 2014, 03:43
Default
  #2
Senior Member
 
dkxls's Avatar
 
Armin
Join Date: Feb 2011
Location: Helsinki, Finland
Posts: 156
Rep Power: 17
dkxls will become famous soon enough
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.
dkxls is offline   Reply With Quote

Old   May 15, 2014, 05:38
Default
  #3
Member
 
bilal
Join Date: Mar 2014
Location: Germany
Posts: 30
Rep Power: 10
brbbhatti is on a distinguished road
Quote:
Originally Posted by dkxls View Post
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
brbbhatti is offline   Reply With Quote

Old   May 15, 2014, 08:06
Default
  #4
Senior Member
 
dkxls's Avatar
 
Armin
Join Date: Feb 2011
Location: Helsinki, Finland
Posts: 156
Rep Power: 17
dkxls will become famous soon enough
Quote:
Originally Posted by brbbhatti View Post
Thanks. i am using Liquid Evaporation Model.
Which one?

Quote:
Originally Posted by brbbhatti View Post
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:
Code:
    phaseChangeModel liquidEvaporationBoil;
Quote:
Originally Posted by brbbhatti View Post
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.
dkxls is offline   Reply With Quote

Old   May 15, 2014, 09:35
Default
  #5
Member
 
bilal
Join Date: Mar 2014
Location: Germany
Posts: 30
Rep Power: 10
brbbhatti is on a distinguished road
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 is offline   Reply With Quote

Old   May 15, 2014, 10:09
Default
  #6
Member
 
bilal
Join Date: Mar 2014
Location: Germany
Posts: 30
Rep Power: 10
brbbhatti is on a distinguished road
Quote:
Originally Posted by dkxls View Post
Which one?


You should use this one for spray simulation:
Code:
    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

Quote:
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

Quote:
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
brbbhatti is offline   Reply With Quote

Old   May 15, 2014, 10:47
Default
  #7
Senior Member
 
dkxls's Avatar
 
Armin
Join Date: Feb 2011
Location: Helsinki, Finland
Posts: 156
Rep Power: 17
dkxls will become famous soon enough
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
dkxls is offline   Reply With Quote

Old   May 15, 2014, 11:45
Default
  #8
Member
 
bilal
Join Date: Mar 2014
Location: Germany
Posts: 30
Rep Power: 10
brbbhatti is on a distinguished road
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
brbbhatti is offline   Reply With Quote

Old   June 16, 2014, 08:36
Default
  #9
New Member
 
Werner
Join Date: Apr 2014
Posts: 19
Rep Power: 10
Polli is on a distinguished road
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.
Polli is offline   Reply With Quote

Old   June 16, 2014, 09:30
Default
  #10
Member
 
bilal
Join Date: Mar 2014
Location: Germany
Posts: 30
Rep Power: 10
brbbhatti is on a distinguished road
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
brbbhatti is offline   Reply With Quote

Old   June 16, 2014, 09:46
Default
  #11
Senior Member
 
dkxls's Avatar
 
Armin
Join Date: Feb 2011
Location: Helsinki, Finland
Posts: 156
Rep Power: 17
dkxls will become famous soon enough
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
dkxls is offline   Reply With Quote

Old   June 16, 2014, 10:40
Default
  #12
Member
 
bilal
Join Date: Mar 2014
Location: Germany
Posts: 30
Rep Power: 10
brbbhatti is on a distinguished road
Hi Armin,

Thanks for the updates, i will try it and let you know.
brbbhatti is offline   Reply With Quote

Reply

Tags
evaporation model, problem with of 2.3.0

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
model evaporation tfe FLOW-3D 9 July 6, 2015 07:05
Droplet evaporation in fluidised bed. Help with the model? nb92 FLUENT 0 February 27, 2014 05:17
Locating near wall cells for evaporation model mehrdad_kbg OpenFOAM Programming & Development 2 February 4, 2014 03:56
Liquid rocket engine non-premixed combustion model Erik FLUENT 2 November 28, 2013 15:09
multiphase model for evaporation with volatility james FLUENT 0 October 6, 2004 04:55


All times are GMT -4. The time now is 14:50.