|March 31, 2013, 12:43||
density calculation in cold flow and reacting flow solver
Join Date: Nov 2012
Posts: 168Rep Power: 4
I noticed there is a small difference about the density calculation in rhoPimpleFoam and fireFoam solvers:
For rhoPimpleFoam, in pEqn.H, density is updated before and after solving the pressure equations: rho = thermo.rho(); And then in the file rhoPimpleFoam.C, after the turbulence is updated turbulence->correct(), the density will not be caluclated in this current time loop.
For fireFoam, in pEqn.H, rho = thermo.rho() is enfored before the pressure Poison equation is solved. However, after the pressure in pEq.H, no density update. And then calculate turbulence in fireFoam.C, after that the density is calculated in that source file.
I not know the reason why this changed was made. Can anybody say something about this problem?
|Thread||Thread Starter||Forum||Replies||Last Post|
|Laminar Stiff reacting flow won't converge||aarongoodwin||FLUENT||1||March 13, 2012 20:13|
|fortran code for reacting flow||mohamad||Main CFD Forum||0||December 21, 2011 02:05|
|MutiSpecies Reacting Flow||sanjibdsharma||OpenFOAM Running, Solving & CFD||1||August 2, 2009 08:14|
|PISO for reacting flow||newcomer_ to_ CFD||Main CFD Forum||0||September 25, 2005 18:49|
|regarding reacting flow with coupled solver||kiran||FLUENT||2||April 15, 2005 04:57|