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

Xoodles PitzDaily out of temperature rangepersisten error

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 10, 2007, 04:56
Default Hi in a domain of the scale fo
  #1
Member
 
O R
Join Date: Mar 2009
Posts: 50
Rep Power: 8
arkangel is on a distinguished road
Hi in a domain of the scale found in pitzDaily, What are the parameters to be chosen in ignites sites in order to ensure a complete combustion.

All is running stable and then an error
FOAM FATAL ERROR : attempt to use janafThermo<equationofstate> out of temperature range 200 -> 6000;

I saw a lot of people having this error , so far noone has a satisfactory answer

the Courant number is normally less than 1.5
and b -min(b)- is bound to be close to 0 (O(e-12)) with a combustion progress of 50 or 60 %

What else can be wrong
arkangel is offline   Reply With Quote

Old   April 11, 2007, 04:05
Default No one?
  #2
Member
 
O R
Join Date: Mar 2009
Posts: 50
Rep Power: 8
arkangel is on a distinguished road
No one?
arkangel is offline   Reply With Quote

Old   April 11, 2007, 09:10
Default To avoid an out of range error
  #3
Member
 
E. David Huckaby
Join Date: Mar 2009
Posts: 57
Rep Power: 8
david_h is on a distinguished road
To avoid an out of range error you might need to reduce your courant number below 0.5. Weller et. al. (Proc. of Comb. Inst. 1998) mention that a CFL < 0.5 is required for stability, CFL < 0.2 for accuracy. This has also been mentioned on the message board.

Are you using a fixed time-step or having the solver select a time-step based on a CFL condition ?

To alter the ignition you might try adding additional ignition sites and/or changing the parameters of the existing ignition sites.

Dave
david_h is offline   Reply With Quote

Old   April 26, 2007, 03:18
Default Hi Thanks one last que
  #4
Member
 
O R
Join Date: Mar 2009
Posts: 50
Rep Power: 8
arkangel is on a distinguished road
Hi

Thanks

one last question how can i do this
"having the solver select a time-step based on a CFL condition "

till a fixed time step
arkangel is offline   Reply With Quote

Old   April 26, 2007, 08:54
Default To satisfy CFL < 0.2 you can m
  #5
Member
 
E. David Huckaby
Join Date: Mar 2009
Posts: 57
Rep Power: 8
david_h is on a distinguished road
To satisfy CFL < 0.2 you can make the following changes to "system/controlDict":

adjustTimeStep yes;
maxCo 0.2;
david_h is offline   Reply With Quote

Old   April 27, 2007, 10:59
Default hi Thanks it is not working
  #6
Member
 
O R
Join Date: Mar 2009
Posts: 50
Rep Power: 8
arkangel is on a distinguished road
hi Thanks

it is not working it is like i havent done anything , maybe i mistype something
here is my controlDict

application Xoodles;

startFrom startTime;

startTime 0;

stopAt endTime;

endTime 0.3;

deltaT 1e-05;

writeControl runTime;

writeInterval 0.005;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression uncompressed;

timeFormat general;

timePrecision 6;

graphFormat raw;

runTimeModifiable yes;

adjustTimeStep yes;

maxCo 0.2;
arkangel is offline   Reply With Quote

Old   April 27, 2007, 11:04
Default Hi, what happens by using
  #7
Member
 
Tommaso Lucchini
Join Date: Mar 2009
Posts: 82
Rep Power: 8
lucchini is on a distinguished road
Hi,

what happens by using maxCo 0.1? Is it crashing at the same point?

Try also to tighten the tolerances of the fvSolution by an order or two of magnitude and see what happens.

Bye
lucchini is offline   Reply With Quote

Old   April 27, 2007, 11:22
Default Hi Tommaso , David I mean I
  #8
Member
 
O R
Join Date: Mar 2009
Posts: 50
Rep Power: 8
arkangel is on a distinguished road
Hi Tommaso , David

I mean I am not having control over the Time step, it is the same fiexed tiem step

but i am trying all your reomendations during this weekend
arkangel is offline   Reply With Quote

Old   April 27, 2007, 11:33
Default Hi again, the timeStep rem
  #9
Member
 
Tommaso Lucchini
Join Date: Mar 2009
Posts: 82
Rep Power: 8
lucchini is on a distinguished road
Hi again,

the timeStep remains fixed because the Xoodles solver does not include the following lines in the time loop.


# include "readTimeControls.H"
# include "readPISOControls.H"

and after

# include "compressibleCourantNo.H"

also this one is needed

# include "setDeltaT.H"

Furthermore,

# include "setInitialDeltaT.H"

Should be pasted before the start of the time loop.

Try to have a look to the XiFoam application and then to the Xoodles application and try to put the correct files where needed to set the deltaT.

I hope I have been useful.

Bye

Tommaso
lucchini is offline   Reply With Quote

Old   April 30, 2007, 07:07
Default Hi Tommaso, thanks for your
  #10
Member
 
O R
Join Date: Mar 2009
Posts: 50
Rep Power: 8
arkangel is on a distinguished road
Hi Tommaso,
thanks for your help i am testing now seems to be working

for someone else who might be interested
# include "readTimeControls.H"
should be placed before
# include "setInitialDeltaT.H" too
arkangel is offline   Reply With Quote

Old   May 10, 2007, 12:41
Default Hi me again ! I want to kno
  #11
Member
 
O R
Join Date: Mar 2009
Posts: 50
Rep Power: 8
arkangel is on a distinguished road
Hi me again !

I want to know what is the relationship between these parameters in ignitions sites field:

duration 0.01;
strength 4;


I pretty sure that duration is the time in sec. during which the sparkle takes place. but if I am wrong (please correct me) . If I increase the duration time i also have to increase the strength in order to get the same result , the point is that i have no idea what is the relationship , what is the meaning of strength anyways

can someone enlighten me ?
arkangel is offline   Reply With Quote

Old   May 10, 2007, 16:04
Default Hi O R (?!?!) the ignition
  #12
Member
 
Tommaso Lucchini
Join Date: Mar 2009
Posts: 82
Rep Power: 8
lucchini is on a distinguished road
Hi O R (?!?!)

the ignition treatment works as follows:

for its duration, a source term which is proportional to the strength is added to the equation for b. Have a look at the ignite.H file to know more about that.

Bye

Tommaso
lucchini is offline   Reply With Quote

Old   May 14, 2007, 07:28
Default Hi I am still having the s
  #13
Member
 
O R
Join Date: Mar 2009
Posts: 50
Rep Power: 8
arkangel is on a distinguished road
Hi

I am still having the same error

FOAM FATAL ERROR : attempt to use janafThermo<equationofstate> out of temperature range 200 -> 5000; T = -736.126

the maxCo is 0.1 (in fact the last max Courant was 0.099034) the last residual for p, Ux, Uy,and k were 2.55e-08,2.44e-08,6.55e-09, 7.5e-10.

Something I can not still understand. In the previous time step the min temp either in T or Tu field were bigger than 200 , and in the next time step,(with a deltaT = 8.14815e-08 ) in one point the temperature simply dropped to -736.126 (how could this change happened in such small deltaT and Courant number)

and 3 time steps before the error i get Solution singularity for Xi and b , not sure if this is directly related , because I got the same singularity solution warning for Xi in the first Time Step
arkangel is offline   Reply With Quote

Old   June 7, 2007, 09:16
Default BUMP
  #14
Member
 
O R
Join Date: Mar 2009
Posts: 50
Rep Power: 8
arkangel is on a distinguished road
BUMP
arkangel is offline   Reply With Quote

Old   August 13, 2007, 06:07
Default BUMP !! no one
  #15
Member
 
O R
Join Date: Mar 2009
Posts: 50
Rep Power: 8
arkangel is on a distinguished road
BUMP !! no one
arkangel is offline   Reply With Quote

Old   August 14, 2007, 02:28
Default Hello, I had the same error
  #16
Senior Member
 
Markus Rehm
Join Date: Mar 2009
Location: Erlangen (Germany)
Posts: 176
Rep Power: 8
markusrehm is on a distinguished road
Hello,

I had the same error using reactingFoam. I limited the enthalpy as suggested by Markus Hartinger (http://www.cfd-online.com/OpenFOAM_D...tml?1177099819). When I was at a certain stage I could continue my calculation with the (unlimited) standard solver.

Don't know if that might help.

Regards Markus
markusrehm is offline   Reply With Quote

Old   February 11, 2011, 05:32
Default
  #17
fox
New Member
 
sebastian fox
Join Date: Feb 2011
Posts: 1
Rep Power: 0
fox is on a distinguished road
Hi Foamers,
im a new foamer and wondering if anybody has found a satisfying answer to the problem with the error from janafThermo:

Im currently working on a combustion induced vortex breakdown with XiFoam respectively myXiFoam (within this solver i have changed the calculation of the desity to be only depending on the compressibility) . Im using a 2D mesh (Axi-symmetric geometry using the wedge patch type) and running with low Courrant numbers 0.2<Co<0.5.

Anyway i get for both solvers the error message

-> FOAM FATAL ERROR: attempt to use janafThermo<equationOfState> out of temperature range 100 -> 5000; T = 99.3763

i have tryed to change Co number, fvSchemes (Crank-Nicholson, Upwind), time step and mesh density. So far, all attempts failed and I dont find any satisfying answer within the posts been made.

So please I would warmly and deeply appreciate any kind of help to solve this problem

greets
fox is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
PitzDaily example and simpleFoam csomerlot OpenFOAM Running, Solving & CFD 2 November 6, 2013 12:20
Run pitzdaily with xifoam tavasoly OpenFOAM Running, Solving & CFD 0 July 8, 2008 04:25
PitzDaily tutorial with Xoodles christianvhoersten OpenFOAM Running, Solving & CFD 0 January 11, 2008 11:17
How to run the cases Xoodles%5cpitzDaily3D and Xoodles%5cpitzDaily cfdfans OpenFOAM Running, Solving & CFD 4 October 18, 2007 03:47
Xoodles temperature change without combustion hannes OpenFOAM Running, Solving & CFD 5 October 17, 2005 16:31


All times are GMT -4. The time now is 17:54.