CFD Online Discussion Forums - Development of an All mach number fractional step method

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Programming & Development (https://www.cfd-online.com/Forums/openfoam-programming-development/)

- - Development of an All mach number fractional step method (https://www.cfd-online.com/Forums/openfoam-programming-development/181549-development-all-mach-number-fractional-step-method.html)

Development of an All mach number fractional step method

Dear Foamers,

I have recently come across this paper "A Primitive Variable Central Flux Scheme for All
Mach Number Flows" that triggered me to realise this solver in OpenFOAM.
Here is reference for this paper:

Shi, Ke; Morris, Scott; Jemcov, Aleksandar
A Primitive Variable Central Flux Scheme for All Mach Number Flows (Conference)
53rd AIAA Aerospace Sciences Meeting, 2015, At Kissimmee, FL, 2015.

However, I am stuck at the central flux part:

Code:

surfaceScalarField rho_left

(

     "rho_left",

      fvc::interpolate(rho, pos, "reconstruct(rho)")

);

surfaceScalarField rho_right

(

     "rho_right",

     fvc::interpolate(rho, neg, "reconstruct(rho)")

);



surfaceVectorField U_left

(

     "U_left",

     fvc::interpolate(U, pos, "reconstruct(U)")

);



surfaceVectorField U_right

(

     "U_right",

     fvc::interpolate(U, neg, "reconstruct(U)")

);



surfaceScalarField phil((rho_left*U_left) & mesh.Sf());

surfaceScalarField phir((rho_right*U_right) & mesh.Sf());

surfaceScalarField pRML(runTime.deltaT()*fvc::snGrad(p)*mesh.magSf());

surfaceScalarField phif(0.5*(phil + phir) + pRML);

which corresponds to Eqs 11-17 of the paper. With this implementation, the continuity error is at an order of 1e-02 which is not at all acceptable. I have attached the solver and a forwardStep case for you to help me with this problem.

Many thanks,
Carlen

Can you put the paper?

Sent from my GT-I8190L using CFD Online Forum mobile app

Quote:

Originally Posted by Santiago (Post 630106)

Can you put the paper?

Sent from my GT-I8190L using CFD Online Forum mobile app

Here is the link for this paper: http://sci-hub.cc/10.2514/6.2015-1265

Thanks! I'll check it out this weekend and see how I can give you a hand

Sent from my GT-I8190L using CFD Online Forum mobile app

Anyone would like to give me a hand:)