Is rhoCentralFoam correctly implemented?
 

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?
Attachment 54540

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
Rdf