|June 20, 2011, 13:55||
Tmax in reactingFoam depends on CFL-number (?)
Join Date: Mar 2009
Posts: 29Rep Power: 8
I am using OpenFoam-1.7.1 (git: 03e7e056c215). There is a tight dependence between "T gas max" and the CFL number in the reactingFoam tutorial case. Lower Courant-numbers lead to lower maximum temperatures and vice versa.
In "applications/solvers/combustion/reactingFoam/chemistry.H" the line
// Chalmers PaSR model
kappa = (runTime.deltaT() + tc)/(runTime.deltaT() + tc + tk);
with an _absolute_ runTime.deltaT() clearly is an error, isn't it?
Also to test mesh dependence i tried this:
I refined the original mesh from the tutorial from 4000 -> 400000 cells (x10 in both directions) and ran the tutorial again. It crashes both in single and parallel at around 0.0007s. I also tried removing runTime.deltaT() from the kappa-equation in the source code, recompiled it, and ran it to the same effect. At some point epsilon becomes negative and the solver crashes.
Last edited by andras; July 5, 2011 at 08:15.
|cfl-number, reactingfoam, temperature|
|Thread||Thread Starter||Forum||Replies||Last Post|
|unable to install openfoam 1.7 using apt-get||maazhmd||OpenFOAM Installation||7||July 24, 2011 11:18|
|snappyHexMesh won't work - zeros everywhere!||sc298||OpenFOAM Native Meshers: snappyHexMesh and Others||2||March 27, 2011 21:11|
|what is CFL number?||cfdiscool||FLUENT||3||February 11, 2010 03:56|
|external flow with snappyHexMesh||chelvistero||OpenFOAM||11||January 15, 2010 20:43|
|CFL number and time accuracy for LES||Li Yang||Main CFD Forum||2||August 1, 2002 06:11|