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

Bug in energy equation of reactingFoam (openfoam 6)

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 6, 2018, 05:28
Default Bug in energy equation of reactingFoam (openfoam 6)
  #1
Member
 
benoit favier
Join Date: Jun 2017
Posts: 64
Rep Power: 8
Benben is on a distinguished road
Hello,

I found a weird comportement of reactingFoam in openfoam 6 (maybe is it similar in older versions).

I setup a case with H2O and N2 as species, reaction off, and a closed box for the mesh.
On every walls of the box, boundary conditions are "no slip" for U, "zero gradient" for P and T.
The only thing is that i initialized the interior of the box to full N2 (Y[H2O] = 0.), and set the walls to a fixed fraction of H2O (Y[H2O] = 0.2).
This way the H2O should diffuse inside of the box, and there should be no convective flux.

The H2O is diffusing correctly but the temperature goes crazy for no reason.

My guess is that there should be a molecular diffusivity term in the energy equation, and there is none. So the species that diffuse dont bring with them there enthalpy.
I attached the (very light) test case for people to be able to test it. You just have to run "reactingFoam" inside, the mesh is already generated. It will create only 5 time steps.
Box.zip
Benben is offline   Reply With Quote

Old   September 6, 2018, 05:41
Default
  #2
Member
 
benoit favier
Join Date: Jun 2017
Posts: 64
Rep Power: 8
Benben is on a distinguished road
I found this thread :
having trouble using reactingFoam with reactions turned off
That talks about a similar problem.

"reactingFoam algorithm is an efficient, approximate method for turbulent reacting flow. However, in the laminar limit the approximation breaks down, leading to temperature inconsistencies.

In the laminar case, if you have a composition gradient, the only way to ensure that the diffusion does not give rise to temperature errors is by enforcing a Lewis no = 1 constraint, e.g. update the sensible enthalpy eqn to:

fvScalarMatrix hsEqn
(
fvm::ddt(rho, hs)
+ mvConvection->fvmDiv(phi, hs)
// - fvm::laplacian(turbulence->alphaEff(), hs)
- fvm::laplacian(turbulence->muEff(), hs) // unit lewis no.
==
DpDt
+ chemistrySh
);
"
Benben is offline   Reply With Quote

Old   September 6, 2018, 10:23
Default
  #3
Member
 
benoit favier
Join Date: Jun 2017
Posts: 64
Rep Power: 8
Benben is on a distinguished road
This modification in the code of the solver doesnt resolve the problem
Benben is offline   Reply With Quote

Old   January 14, 2019, 07:38
Default I encountered the same problem use OpenFOAM-v1612+
  #4
New Member
 
haixiang zeng
Join Date: Sep 2018
Posts: 1
Rep Power: 0
zeng-ujs is on a distinguished road
Quote:
Originally Posted by Benben View Post
This modification in the code of the solver doesnt resolve the problem
Hi,Benben.

Recently i met the same problem when i ignition with pilot flame . the temperature dropped crazy and total reached less than 200K. Firstly my thought was there maybe something wrong with boundary, so i set a 2D case to exam it , and amazing everything is fine . i have no idea what should i do, do you have solved this problem?
zeng-ujs is offline   Reply With Quote

Reply

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
Question regarding Energy equation in reactingFoam Dan1788 OpenFOAM Programming & Development 6 August 5, 2014 00:22
Energy equation in OpenFoam nwpukaka OpenFOAM Programming & Development 2 June 24, 2014 02:40
problem on energy and flow equation solving separately preetam69 FLUENT 2 July 30, 2013 21:08
error message cuteapathy CFX 14 March 20, 2012 06:45
Why FVM for high-Re flows? Zhong Lei Main CFD Forum 23 May 14, 1999 13:22


All times are GMT -4. The time now is 10:05.