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

Problem With reactingFOAM

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 24, 2010, 03:26
Default Problem With reactingFOAM
  #1
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
Hi,
I am trying to run a case using reactingFOAM. But I am getting the following error:

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

From function janafThermo<equationofstate>::checkT(const scalar T) const
in file /home/dm2/henry/OpenFOAM/OpenFOAM-1.2/src/thermophysicalModels/specie/lnInclude/ janafThermoI.H at line 73

The temperatures that I have supplied as BC are 1000K and 600K for O2 and H2 respectively. (Its H2-O2 combustion.)

My Courant No. is 0.2 and my max. cell skewness = 0.66838.
My "checkMesh" results are all OK.

Can anybody please tell me where am I going wrong?
nakul is offline   Reply With Quote

Old   November 24, 2010, 04:02
Default
  #2
Senior Member
 
Karl-Johan Nogenmyr
Join Date: Mar 2009
Location: Linköping
Posts: 249
Rep Power: 12
kalle is on a distinguished road
Hi,

this is a common problem. For me it was always caused by erroneous BC's.

Though this user does not get janaf errors, you can see how BC's can be set:

reactingFoam boundary conditions (inletOutlet BC at outlets)

Other problems may be misconfigured inlets/initial condition. Please check that initally every cell and face have a proper mixture.

K
kalle is offline   Reply With Quote

Old   November 24, 2010, 07:34
Default
  #3
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
Thanx Karl for such an early reply. I changed my BC accordingly.

But I have few doubts :
1) What should be the BC for k and epsilon as my problem also employs k-epsilon turbulence model?

2) At the inlet for air, the air is moving with certain velocity so wether the BC for pressure should be "total pressure" or the "fixedValue" of known static presure? The airspeed at inlet is known to me.

3) Moreover my case also includes compressibiltiy. So in your opinion should there be any changes in BC.

Actually I didn't specified exactly the BC as you directed in that link, rather I have tweaked them a little bit, because the solver was giving the same error earlier than the time, when it used to blow with previously specified BC. So could you please tell me the logic behind changing these BC or you may refer me to some study material if possible!!
nakul is offline   Reply With Quote

Old   November 25, 2010, 04:01
Default
  #4
Senior Member
 
Karl-Johan Nogenmyr
Join Date: Mar 2009
Location: Linköping
Posts: 249
Rep Power: 12
kalle is on a distinguished road
Hi,

1) For k and eps, I would use the same types as for other scalars (except p), that is: prescribed value at inlet and inletOutlet at outlets.

2, 3) For incompressible (in this case low-Mach number cases) or weakly compressible, I would use zeroGradient at inlet and totalPressure at outlet. If you are running hi-Mach numbers, you would need other BC's. Can't help you much there, unfortunately. If you by compressibility mean mainly density dependence on temperature, and pressure fluctuations remain low, I would go for a low-Mach number formulation, where rho(T, Yi, h) and you introduce a constant global pressure in the state equation.

Logic behind BC's: To have a well posed system which gives you what you want could be a pragmatic logic :-)

K
kalle is offline   Reply With Quote

Old   November 25, 2010, 05:44
Default
  #5
Senior Member
 
Nakul
Join Date: Apr 2010
Location: India
Posts: 147
Rep Power: 7
nakul is on a distinguished road
Thanx Karl !

Actually I do have high Mach no.flow of order of Mach No. =2.5.

Do you think "waveTransmissive" BC for pressure in super-sonic flow good?

Could you tell me any references, if possible?

If anybody else can help, it would be greatly appreciated!!!
nakul is offline   Reply With Quote

Old   October 4, 2013, 05:01
Default mass fraction of products
  #6
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 3
yash.aesi is on a distinguished road
Greetings oll,
I am using OF-2.2 and tried to solve my case using reactingFoam solver with PasR model . i am getting the temperature profile almost as required but i am not getting the mass fraction of the products as desired . it coming almost double of what i desire.

i checked my BC's many times and i don't think there is any mistake . i am attaching the files plz have a look .
so can anybody tell me then where i am doing wrong ??

Thanks in advance


Regards ,
Sonu
Attached Files
File Type: zip 0.zip (6.3 KB, 10 views)
File Type: zip constant1.zip (4.7 KB, 16 views)
yash.aesi 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
UDF compiling problem Wouter Fluent UDF and Scheme Programming 6 June 6, 2012 04:43
Incoherent problem table in hollow-fiber spinning Gianni FLUENT 0 April 5, 2008 10:33
natural convection problem for a CHT problem Se-Hee CFX 2 June 10, 2007 06:29
Adiabatic and Rotating wall (Convection problem) ParodDav CFX 5 April 29, 2007 19:13
Is this problem well posed? Thomas P. Abraham Main CFD Forum 5 September 8, 1999 14:52


All times are GMT -4. The time now is 15:20.