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

Requesting help in LES set up for SandiaD flame

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

Like Tree1Likes
  • 1 Post By arvindpj

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 27, 2018, 01:49
Default Requesting help in LES set up for SandiaD flame
  #1
New Member
 
Divyanshu Solanki
Join Date: May 2018
Posts: 1
Rep Power: 0
Divyanshu01 is on a distinguished road
Hello all !
I have been trying to simulate Sandia D flame using LES, but an error message always shows up, which is,


--> FOAM Warning :
From function Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::sc alar) const [with EquationOfState = Foam::perfectGas<Foam::specie>; Foam::scalar = double]
in file /home/ubuntu/OpenFOAM/OpenFOAM-5.x/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 117
attempt to use janafThermo<EquationOfState> out of temperature range 250 -> 3000; T = -184.276


and a series of such messages with different values of T. I can see that T is diverging but I can not figure out the cause.


The program then terminates with the message,


#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::combustionModels::EDC<Foam::combustionModels ::psiChemistryCombustion>::correct() at ??:?
#4 ? at ??:?
#5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6 ? at ??:?
Floating point exception (core dumped)




Please help me understand the cause of divergence and meaning of this last message.


Thanks.
Divyanshu01 is offline   Reply With Quote

Old   May 31, 2018, 10:27
Default
  #2
Member
 
Arvind Jay
Join Date: Sep 2012
Posts: 96
Rep Power: 12
arvindpj is on a distinguished road
Quote:
Originally Posted by Divyanshu01 View Post
Hello all !
I have been trying to simulate Sandia D flame using LES, but an error message always shows up, which is,


--> FOAM Warning :
From function Foam::scalar Foam::janafThermo<EquationOfState>::limit(Foam::sc alar) const [with EquationOfState = Foam:erfectGas<Foam::specie>; Foam::scalar = double]
in file /home/ubuntu/OpenFOAM/OpenFOAM-5.x/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 117
attempt to use janafThermo<EquationOfState> out of temperature range 250 -> 3000; T = -184.276


and a series of such messages with different values of T. I can see that T is diverging but I can not figure out the cause.


The program then terminates with the message,


#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::combustionModels::EDC<Foam::combustionModels :siChemistryCombustion>::correct() at ??:?
#4 ? at ??:?
#5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6 ? at ??:?
Floating point exception (core dumped)




Please help me understand the cause of divergence and meaning of this last message.


Thanks.

Use a temperature limiter in fvObjects to constrain the temperature between Tmin and Tmax. But this is usually a fix, the main cause may be due to large time step or some of the reaction ODE settings.

Cheers,
Jay
Zhiheng Wang likes this.
arvindpj is offline   Reply With Quote

Old   July 19, 2018, 14:17
Default Check combustion model and and boundary Condition
  #3
Member
 
Zhiheng Wang
Join Date: Mar 2016
Posts: 72
Rep Power: 8
Zhiheng Wang is on a distinguished road
Quote:
Originally Posted by arvindpj View Post
Use a temperature limiter in fvObjects to constrain the temperature between Tmin and Tmax. But this is usually a fix, the main cause may be due to large time step or some of the reaction ODE settings.

Cheers,
Jay
Hello,
Please give a brief about combustion model, constant/reaction file , chemistryProperties file( in case of Arrehinius rate) and boundary condition.
it seems your temperature is going down compared to the minimum temperature in thermo.gasproperties may be you are using wrong combustion parameters.
Zhiheng Wang is offline   Reply With Quote

Old   February 1, 2021, 01:51
Default Issues in LES of Sandia Flame D with FGMFoam
  #4
New Member
 
Sandeep
Join Date: Apr 2017
Location: IIT Delhi, New Delhi, India
Posts: 10
Rep Power: 7
Sandeep lamba is on a distinguished road
Hi all,

I am trying to simulate Sandia Flame D with FGMFoam, I have already done RANS of the same and got desired results. Now as I am trying to simulate it with LES I have done some fixes such as Turbulent inlet conditions for Velocity and changing turbulent model to LES (kEqn and Smagorinsky ). Also I have worked on the mesh refinement. But, still I am not getting acceptable results infact getting very poor results. (some contours are attached ). The simulation starts very well but only up to a small time (0.045 Sec). But later on the dissipation stops as can be seen in the contours attached.


Someone please comment on the problem.

Thanks in advance.


Kind regards
T_nut_Ch4_0.4.jpg

T_nut_Ch4_0.03.jpg

T_nut_Ch4_0.05.jpg

T_nut_Ch4_0.06png.jpg
Sandeep lamba is offline   Reply With Quote

Old   February 28, 2021, 00:37
Default
  #5
New Member
 
Zahid Hussain
Join Date: Jan 2016
Posts: 1
Rep Power: 0
STHX is on a distinguished road
Hi Sandeep,

I tried similar LES this weekend. It seems that LES results are not correct for following reasons.
1) Grid is not refined for LES (in case you are using similar grid size as RANS tutorial)
2) LES is strictly for 3D simulations.

Because of second reason, 2D LES will be wrong irrespective of how refined your grid is.

Hope it helps. I have not tried 3D LES yet, for it takes more computation resources.
STHX is offline   Reply With Quote

Old   March 1, 2021, 00:16
Default
  #6
New Member
 
Sandeep
Join Date: Apr 2017
Location: IIT Delhi, New Delhi, India
Posts: 10
Rep Power: 7
Sandeep lamba is on a distinguished road
Thanks for the reply and suggestion. Now I am moving to 3D and first trying RANS. I have made 45 Degree segment of Sandia flame (mesh pic link https://drive.google.com/file/d/1Z_o...ew?usp=sharing). Now I am confused which type of boundary condition should be given for the symmetric cut sections. I have used CyclicAMI but when simulation run new error of ''No MRF model persent'' appears. Any hint regarding this please.
Sandeep lamba is offline   Reply With Quote

Old   July 9, 2021, 07:16
Default
  #7
New Member
 
Sandeep
Join Date: Apr 2017
Location: IIT Delhi, New Delhi, India
Posts: 10
Rep Power: 7
Sandeep lamba is on a distinguished road
I did some fixes with my mesh and now using 3D geometry. I am getting some improved results. But, there seems some problems with Pressure fields (pressure wave are visible in domain). the boundary conditions for U is fixed value at inlet and inletoutlet at outlet. For pressure it is zerogradient at inlet and totalPressure at outlet.
The combustion model is FGM and LES model is kEqn.
Anyone please comment.
the link for contours is here:

https://drive.google.com/drive/folde...0y?usp=sharing

kind regards


Sandeep Lamba
Sandeep lamba is offline   Reply With Quote

Old   July 9, 2021, 09:21
Default
  #8
Senior Member
 
piu58's Avatar
 
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 727
Rep Power: 13
piu58 is on a distinguished road
The inletoutlet b.c. may be problematic. I recommend at least make a few experiments with zero gradient first.

LES is very demanding in all of the model: Mesh, boundary conditions, schemes. So I recommend further starting with conservative schemes firs.

It is always a good idea of using relaxation factors.

You need first stabilize the simulation. At the next step you should get the results in the expected range of values. In your case, this is near the inlet.

If you managed this you may start to get reliable results for the more problematic regions of your simulation.
__________________
Uwe Pilz
--
Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950)
piu58 is offline   Reply With Quote

Old   July 10, 2021, 03:04
Default
  #9
New Member
 
Sandeep
Join Date: Apr 2017
Location: IIT Delhi, New Delhi, India
Posts: 10
Rep Power: 7
Sandeep lamba is on a distinguished road
Thank you Uwe Pilz for quick reply. I used zeroGradient at outlet and also relaxation factor as you told and it worked. There are no more pressure wave, other fields are also improved.
Again thanks a lot.


kind regards

Sandeep Lamba
Sandeep lamba is offline   Reply With Quote

Old   January 25, 2022, 07:59
Default
  #10
Member
 
Amirreza Niazmehr
Join Date: Nov 2018
Posts: 40
Rep Power: 5
Amirreza_pro is on a distinguished road
Hi guys, I'm new to openFoam and I have a question about simulation of Sandia Flame D.
I have already simulated Sandia_D with RAS in openFoam.
Now I'm going to simulate Sandia_D with LES.
But I don't know how to change the set up to get the results for LES.
Should I change the setting for the RAS model or there is alreay a test case in toturial folders.
I would be appreciated if someone can help me by this case.
Amirreza_pro is offline   Reply With Quote

Reply

Tags
flame model, les, sandia

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
[ICEM] ICEM Scripting Issues tylerplowright ANSYS Meshing & Geometry 33 September 27, 2021 16:35
Remeshing for Rotation in CFX via ICEMCFD gera CFX 3 September 23, 2017 12:29
Ansys Licence Serve on Ubuntu 16.04 LTS david.pasquale ANSYS 2 January 20, 2017 11:52
[ICEM] Question about the use of ICEM tcl scripting lnk ANSYS Meshing & Geometry 9 January 16, 2013 09:24
Help with GNUPlot Renato. Main CFD Forum 6 June 6, 2007 19:51


All times are GMT -4. The time now is 00:55.