|
[Sponsors] |
Bug in energy equation of reactingFoam (openfoam 6) |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
September 6, 2018, 06:28 |
Bug in energy equation of reactingFoam (openfoam 6)
|
#1 |
Member
benoit favier
Join Date: Jun 2017
Posts: 64
Rep Power: 9 |
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 |
|
September 6, 2018, 06:41 |
|
#2 |
Member
benoit favier
Join Date: Jun 2017
Posts: 64
Rep Power: 9 |
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 ); " |
|
September 6, 2018, 11:23 |
|
#3 |
Member
benoit favier
Join Date: Jun 2017
Posts: 64
Rep Power: 9 |
This modification in the code of the solver doesnt resolve the problem
|
|
January 14, 2019, 08:38 |
I encountered the same problem use OpenFOAM-v1612+
|
#4 | |
New Member
haixiang zeng
Join Date: Sep 2018
Posts: 1
Rep Power: 0 |
Quote:
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? |
||
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Question regarding Energy equation in reactingFoam | Dan1788 | OpenFOAM Programming & Development | 6 | August 5, 2014 01:22 |
Energy equation in OpenFoam | nwpukaka | OpenFOAM Programming & Development | 2 | June 24, 2014 03:40 |
problem on energy and flow equation solving separately | preetam69 | FLUENT | 2 | July 30, 2013 22:08 |
error message | cuteapathy | CFX | 14 | March 20, 2012 07:45 |
Why FVM for high-Re flows? | Zhong Lei | Main CFD Forum | 23 | May 14, 1999 14:22 |