CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

phi == pEqn.flux() rhoSimpleFoam

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

LinkBack Thread Tools Search this Thread Display Modes
Old   July 9, 2015, 10:10
Question phi == pEqn.flux() rhoSimpleFoam
New Member
Join Date: Jul 2015
Posts: 3
Rep Power: 8
David1 is on a distinguished road

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()

as well.

I hope somebody can help me with this.
Thank you in advance!
David1 is offline   Reply With Quote


flux, peqn, phi, rhosimplefoam, transonic

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
phi -= pEqn.flux() vs. linearInterpolate(U) & mesh.Sf() santiagomarquezd OpenFOAM Programming & Development 33 October 15, 2019 10:24
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 15:08
Problem with rhoSimpleFoam : exploding enthalpy and density at the walls david39 OpenFOAM Running, Solving & CFD 6 January 18, 2011 11:49
Turbulence Model phi vs phi_ doug OpenFOAM Running, Solving & CFD 4 November 10, 2009 04:33

All times are GMT -4. The time now is 00:31.