CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

How to run the cases Xoodles%5cpitzDaily3D and Xoodles%5cpitzDaily

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 15, 2006, 12:45
Default Dear all I run the Xoodles\
  #1
cfdfans
Guest
 
Posts: n/a
Dear all

I run the Xoodles\pitzDaily3D, use the default 'combustionProperties' values.
and run Xoodles\pitzDaily, set 'ignite' is 'yes',but the two cases do not ignite and combustion. when I increase the value of 'strength' of 'ignitionSites' , the cases ignite,but reach error soon,
"attempt to use janafThermo<equationofstate> out of temperature range 200 -> 5000; T = 192.877 "

I wonder how to run the cases Xoodles\pitzDaily3D and Xoodles\pitzDaily with combustion? what are the meaning of the following Properties? when 'strength' is '2', means how high temperature? and what is the 'ignitionSphereFraction' , 'ignitionThickness' and 'ignitionKernelArea' ?

thanks !

qf

--------------------------------------------------


ignitionSites
(

{
location (0 0 0);
diameter 0.003;
start 0;
duration 0.001;
strength 2;
}
);

ignitionSphereFraction 1;

ignitionThickness ignitionThickness [0 1 0 0 0 0 0] 0;

ignitionCircleFraction 0;

ignitionKernelArea ignitionKernelArea [0 2 0 0 0 0 0] 0;

----------------------------------------------------
  Reply With Quote

Old   May 18, 2006, 08:57
Default Anyone can help me to get the
  #2
cfdfans
Guest
 
Posts: n/a
Anyone can help me to get the correct results of Xoodles\pitzDaily and XoodlespitzDaily3D ?

thanks
  Reply With Quote

Old   February 11, 2007, 17:46
Default BUMP
  #3
Member
 
O R
Join Date: Mar 2009
Posts: 50
Rep Power: 7
arkangel is on a distinguished road
BUMP
arkangel is offline   Reply With Quote

Old   February 16, 2007, 05:11
Default I run the case again this ti
  #4
Member
 
O R
Join Date: Mar 2009
Posts: 50
Rep Power: 7
arkangel is on a distinguished road
I run the case again this time i set the values for boundary conditions to adiabaticWall and run to the end
arkangel is offline   Reply With Quote

Old   October 18, 2007, 04:47
Default Hello qinfei, at least in m
  #5
Member
 
Hannes Kröger
Join Date: Mar 2009
Location: Rostock, Germany
Posts: 92
Rep Power: 8
hannes is on a distinguished road
Hello qinfei,

at least in my case, which had well initialised temperature and combustion progress fields from another run with a different combustion model, I found this problem (janafThermo out of temperature range) to arise in the newton solver in "specieThermo".

I suppose, that due to an unlucky shape of h(T) or because of the break in h(T) and cp(T) at Tcommon, the calculated tangent has a too low positive ascent at certain temperatures and the newton iteration gives very low or even negative temperatures.
When the newton solver tries to evaluate the janaf polynomials for the next iteration, the error occurs.

For a very quick and very dirty fix, I inserted a
if (Tnew<TMIN) return TMIN;
at line 74 in $FOAM_SRC/thermophysicalModels/specie/thermo/specieThermo/specieThermoI.H with TMIN being of course problem specific.

Best regards, Hannes
hannes is online now   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
Xoodles PitzDaily out of temperature rangepersisten error arkangel OpenFOAM Running, Solving & CFD 16 February 11, 2011 05:32
Division by zero in Xoodles hannes OpenFOAM Bugs 3 August 4, 2008 12:04
PitzDaily tutorial with Xoodles christianvhoersten OpenFOAM Running, Solving & CFD 0 January 11, 2008 11:17
Soot model Xoodles deepblue17 OpenFOAM Running, Solving & CFD 0 May 2, 2006 12:57
Xoodles temperature change without combustion hannes OpenFOAM Running, Solving & CFD 5 October 17, 2005 17:31


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