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

[fireFoam] Pyrolysis of a wood crib case, janaf temperature divergence

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 14, 2021, 06:06
Default [fireFoam] Pyrolysis of a wood crib case, janaf temperature divergence
  #1
New Member
 
Mohamed Baraket
Join Date: Mar 2021
Posts: 6
Rep Power: 5
BarakeX is on a distinguished road
Hello
I'm trying to model an experimental model-scale case that consists in burning a wood crib with heptane pools. Case description is shown in this pdf link:

https://www.iafss.org/publications/f...ss_10-1305.pdf

I have made the model of the wood crib and exported it as an stl file, then I took the tutorial case 'oppositeBurningPanels' from the LES fireFoam solver in OpenFoam-v2012 as a base, tweaked the parameters and boundary conditions to match my case, used snappyHexMesh for the wood crib and the compartment meshing.

this is the first model of the wood crib:



I was trying to extrude the wood crib stl as a pyrolysis region with the extrudeToRegionMesh utility. I don't know how to use its other options but I tried to extrude the region around the stl face zone of the wood crib that's inside the internal mesh, and I get extruded cells inside the stl (which is outside the internal mesh).
The extruded cell layers eventually intersect, but I ignored it at that time and tried to run the simulation. it was set for 20 seconds at an adjustable time step.
Eventually the simulation stopped halfway with an error from the janaf table:

Quote:
--> FOAM Warning :
From 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/OpenFOAM-v2012/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 120
attempt to use janafThermo<EquationOfState> out of temperature range 200 -> 5000; T = -1041.53982114
I got outrageous values of the temperature.
I refined the mesh, tweaked the solution schemes to gain stability over precision, and still got the same thing.

I tried to simplify the wood crib model and I made this:



I tried to extrude the panel region the same way using the stl as a face zone.
The image below shows how the intersected cells of the panel region:



The two images below show the compartment and the wood crib (second version) together:






I can't find the root of the divergence (the exact area) but I'm assuming it's the panel region that's causing the problem.

I really appreciate any advice, but I also want to know is there any way to properly extrude a panel/pyrolysis region out of such geometry inside the internal mesh without causing any problems or cell intersections.

There's a kind of a similar case published, but I can't find enough resources on that, so I can at least make a proper panel region.
here's the link:
https://www.sciencedirect.com/scienc...79711217301893

Thanks in advance.

Last edited by BarakeX; April 14, 2021 at 09:09. Reason: png links didn't work
BarakeX is offline   Reply With Quote

Old   March 26, 2022, 01:48
Default
  #2
Member
 
Francisco T
Join Date: Nov 2011
Location: Melbourne, Australia
Posts: 64
Blog Entries: 1
Rep Power: 14
frantov is on a distinguished road
Hi there
is a bit old, but I think the mesh in the depth direction is not fine enough.
What you can do is make a simpler model, where you model one or two rows and columns of the crib, this will allow you to refine the mesh in the depth direction and see if runs well.
frantov is offline   Reply With Quote

Reply

Tags
extrudetoregionmesh, firefoam, janafthermo warning, pyrolysis

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
equationOfState icoTabulated error carl_r OpenFOAM Running, Solving & CFD 3 March 14, 2022 07:05
Fail to converge when solving with a fabricated solution zizhou FLUENT 0 March 22, 2021 07:33
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 07:27
Pitz Daily Test Case temperature measurements panara OpenFOAM Running, Solving & CFD 2 March 22, 2007 04:00
Free surface boudary conditions with SOLA-VOF Fan Main CFD Forum 10 September 9, 2006 13:24


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