|July 9, 2015, 10:10||
phi == pEqn.flux() rhoSimpleFoam
Join Date: Jul 2015
Posts: 3Rep Power: 3
Does anyone know why for the transonic case in rhoSimpleFoam the flux is calculated as:
phi == pEqn.flux(); ?
If transonic is turned off, then the flux is calculated in the 'normal' way, i.e. first calculting the flux without the contribution of the pressure gradient:
phi = fvc::interpolate(rho)*(fvc::interpolate(U) & mesh.Sf() )
And after the calulation of the pressure, the old flux value gets corrected using the new pressure:
phi -= pEqn.flux();
So I do not understand why the flux in the transonic case is not calculated as:
phi == fvc::interpolate(rho)*(fvc::interpolate(U) & mesh.Sf() ) - pEqn.flux()
I hope somebody can help me with this.
Thank you in advance!
|flux, peqn, phi, rhosimplefoam, transonic|
|Thread||Thread Starter||Forum||Replies||Last Post|
|phi -= pEqn.flux() vs. linearInterpolate(U) & mesh.Sf()||santiagomarquezd||OpenFOAM Programming & Development||32||June 12, 2014 01:50|
|mixerVesselAMI2D's mass is not balancing||sharonyue||OpenFOAM Running, Solving & CFD||6||June 10, 2013 09:34|
|pEqn.flux()||cheng1988sjtu||OpenFOAM Running, Solving & CFD||4||November 21, 2012 16:08|
|Problem with rhoSimpleFoam : exploding enthalpy and density at the walls||david39||OpenFOAM Running, Solving & CFD||6||January 18, 2011 12:49|
|Turbulence Model phi vs phi_||doug||OpenFOAM Running, Solving & CFD||4||November 10, 2009 05:33|