|
[Sponsors] |
July 11, 2019, 18:53 |
JanafThermo Warning
|
#1 |
New Member
Miriana De Rose
Join Date: Jun 2019
Posts: 8
Rep Power: 7 |
Hello,
while running my simulation i get the following Warning: --> 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/pawan/OpenFOAM/OpenFOAM-v1812/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 117 attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = -388.868 --> 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/pawan/OpenFOAM/OpenFOAM-v1812/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 117 attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = -111.631 min/max(T) = 200, 1742.55 any idea on the reason why i get that? Thanks. |
|
July 12, 2019, 05:13 |
|
#2 |
Senior Member
Yann
Join Date: Apr 2012
Location: France
Posts: 1,224
Rep Power: 28 |
Hello Miriana,
The warning you get means the temperature somewhere in your domain is out of the valid temperature range of the Janaf tables. As you can see in the warning, the Janaf table is valid for temperature ranging from 200 to 5000K. You seem to get negative temperature which is quite unphysical. It's hard to help you more without knowing what you're trying to simulate, but most of the time this is due to inappropriate boundary conditions leading to divergence. On a side note, please do not post the same question on several topics. I advise you to read the Guide: How to ask a question on the forums. Yann |
|
July 12, 2019, 05:27 |
|
#3 |
New Member
Miriana De Rose
Join Date: Jun 2019
Posts: 8
Rep Power: 7 |
Hello Yann,
thank you for your answer and advice. i know there could be some problems on the BC, but I can't figure it out where exactly. I'm simulating a gas dispersion, I set: - velocity Field value, - velocity inlet value, - Temperature Field value, - Temperature Inlet value. the simulation does not give me warning during the first Time steps, but after a certain time it comes out with the janaf warning. |
|
July 12, 2019, 05:53 |
|
#4 |
Senior Member
Yann
Join Date: Apr 2012
Location: France
Posts: 1,224
Rep Power: 28 |
What solver are you using?
Could you post the detail of the BC you are using? (inlet and outlet) What BC do you use for pressure? |
|
July 12, 2019, 06:14 |
|
#5 |
New Member
Miriana De Rose
Join Date: Jun 2019
Posts: 8
Rep Power: 7 |
I'm using SmallPoolFire3D.
The BC i set are: Velocity - U internalField uniform -5; (because i'm calculating on the opposite orientation of x-axes) sides type pressureInletOutletVelocity; value uniform (0 0 0); walls type noSlip; inlet type fixedValue; value uniform (-0.05 0 0); outlet type pressureInletOutletVelocity; value uniform (0 0 0); Temperature - T internalField uniform 298.15; sides { type inletOutlet; inletValue $internalField; value $internalField; } walls { type zeroGradient; } inlet { type fixedValue; value uniform 313.15; } outlet { type inletOutlet; inletValue $internalField; value $internalField; } Pressure - p internalField uniform 101325; sides { type calculated; value $internalField; } inlet { type calculated; value $internalField; } outlet { type calculated; value $internalField; } } in the folder Constant i changed the g value to dimensions [0 1 -2 0 0 0 0]; value (0 9.8 0); Because after translating the geometry from fluent, it is upside down, so the dispersion is directed in the opposite side of the y-axes. meaning that g is positive according to the y-axes orientation. |
|
July 12, 2019, 07:56 |
|
#6 |
Senior Member
Yann
Join Date: Apr 2012
Location: France
Posts: 1,224
Rep Power: 28 |
Just to be clear, smallPoolFire3D is a tutorial and it uses the fireFoam solver.
What have you modified on the smallPoolFire3D tutorial? As I understand, you changed the geometry and you use one coming from fluent, am I right? How does your domain look like? about your pressure BC, fireFoam use the p_rgh pressure, defined as: This is why all the BC in your p file are set to "calculated". You have to define your BC in the p_rgh file. Could you post more information about your domain and set up? Or maybe upload a light version of your case if possible, so people can give a look at it and maybe help your to solve your problem? You can have a look at this topic to see what kind of information can be useful for other users to help you. In addition, please use the CODE tags when you post your BC or other OpenFOAM files snippets in order to improve the readability of your posts. I am not very familiar with gas dispersion but is there a specific reason for choosing fireFoam over other solvers? Yann |
|
July 12, 2019, 09:47 |
|
#7 |
New Member
Miriana De Rose
Join Date: Jun 2019
Posts: 8
Rep Power: 7 |
Yes, i use smallPoolFire3D tutorial and fireFoam solver!!!
here's the geometry: https://www.dropbox.com/s/5px442yfvm...ageom.PNG?dl=0 The inlet is the top of the highest tank, the outlet is the top of the domain. I modified: Velocity; Temperature; g and of course the mesh. I haven't changed the pressure, cause i need atmospheric pressure inside my domain. I'm using fireFoam because I need further calculation on fire spreading. Thank you. |
|
July 12, 2019, 10:49 |
|
#8 |
Senior Member
Yann
Join Date: Apr 2012
Location: France
Posts: 1,224
Rep Power: 28 |
Ok thank you I understand your choice now.
About your velocity BC on the inlet, you say you set: Code:
type fixedValue; value uniform (-0.05 0 0); Code:
value uniform (0 0.05 0); |
|
July 12, 2019, 11:03 |
|
#9 |
New Member
Miriana De Rose
Join Date: Jun 2019
Posts: 8
Rep Power: 7 |
I put the value
Code:
type fixedValue; value uniform (-0.05 0 0); Do you mean it cannot work if I choose the path to be on the x-axes, while the outlet is on top of the domain, parallel to the x-axes? constant/href it refers to the top of the domain? i didn't know that! thanks a lot. |
|
July 12, 2019, 11:20 |
|
#10 |
Senior Member
Yann
Join Date: Apr 2012
Location: France
Posts: 1,224
Rep Power: 28 |
If your inlet patch is horizontal (normal to y) and your BC imposes the velocity component Uy to be 0 on the patch it means there will be not flux through this patch.
I'm not very familiar with the prghTotalHydrostaticPressure BC but this is what I guess from the description of the BC and from tutorials using this BC. (oppositeBurningPanels for instance) Maybe some other users who are more familiar with it can confirm that. |
|
April 17, 2023, 14:45 |
|
#11 |
New Member
Prakash Bhattarai
Join Date: Mar 2023
Posts: 2
Rep Power: 0 |
Hello MiriDR,
I am facing the same problem from quiet a while could you please provide some information on how did you manage to solve it if you have of course. Thank you in advance |
|
Tags |
janaf, janafthermo, janafthermo warning |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Caffa 3D code | Waliur Rahman | Main CFD Forum | 0 | May 29, 2018 01:53 |
[swak4Foam] installation problem with version 0.2.3 | Claudio87 | OpenFOAM Community Contributions | 9 | May 8, 2013 11:20 |
is internalField(U) equivalent to zeroGradient? | immortality | OpenFOAM Running, Solving & CFD | 7 | March 29, 2013 02:27 |
latest OpenFOAM-1.6.x from git failed to compile | phsieh2005 | OpenFOAM Bugs | 25 | February 9, 2010 05:37 |
[blockMesh] Axisymmetrical mesh | Rasmus Gjesing (Gjesing) | OpenFOAM Meshing & Mesh Conversion | 10 | April 2, 2007 15:00 |