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

Is rhoCentralFoam correctly implemented?

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

Like Tree2Likes
  • 2 Post By randolph

LinkBack Thread Tools Search this Thread Display Modes
Old   March 10, 2017, 07:04
Exclamation Is rhoCentralFoam correctly implemented?
Di Cheng
Join Date: May 2010
Location: Beijing, China
Posts: 47
Rep Power: 13
chengdi is on a distinguished road
I compared the variables in rhoCentralFoam.C and "Implementation of semi-discrete, non-staggered central schemes in a colocated, polyhedral, finite volume framework, for high-speed viscous flows"(Which is a theory guide of rhoCentralFoam.C)

However, I cannot follow the code and paper at the same time. The variables are different. Especially the aphiv_pos and aphiv_neg in the code. I think these variables are corresponded to alpha*phi_f_plus in equation 7 of the paper. However, the formula is not consisitent from code to paper. Can anyone figure out what is wrong?
chengdi is offline   Reply With Quote

Old   January 2, 2019, 10:26
Senior Member
Reviewer #2
Join Date: Jul 2015
Location: Gainesville,FL
Posts: 133
Rep Power: 8
randolph is on a distinguished road
I know this is an old post. But just in case someone else comes across in the future.

The aphiv_pos in the code is "alpha*Phi_f_postive-Omega_f" in equation 7 of the paper.

For the continuity equation (when the variable is rho), the mass flux is calculated as equation 7, which interprets to the following in the code.
phi = aphiv_pos*rho_pos + aphiv_neg*rho_neg;

For the momentum equation, the flux that calculated in the code also includes the contribution of the pressure gradient as the following

surfaceVectorField phiUp
(aphiv_pos*rhoU_pos + aphiv_neg*rhoU_neg)
+ (a_pos*p_pos + a_neg*p_neg)*mesh.Sf()

I hope this help
emjay and febriyan91 like this.
randolph is offline   Reply With Quote


openfoam, rhocentralfoam

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
How to apply a turbulence model to the rhoCentralFoam solver? arussell92 OpenFOAM Pre-Processing 16 December 19, 2020 09:12
Modify rhoCentralFoam: other equations of state fivos OpenFOAM Programming & Development 5 July 29, 2020 14:17
InterDyMFoam+simpleFunctionObject Elham OpenFOAM Running, Solving & CFD 5 July 10, 2017 12:59
rhoCentralFoam Runge Kutta Henning86 OpenFOAM Running, Solving & CFD 24 December 13, 2016 11:12
rhoCentralFoam solver with Slip BCs fails in Parallel Only JLight OpenFOAM Running, Solving & CFD 2 October 11, 2012 22:08

All times are GMT -4. The time now is 02:22.