research/OpenFOAM-5.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

rhoPimpleDyMFoam: Updated transonic formulation for consistency with sonicFoam

Browse source
This commit is contained in:
Henry Weller 2017-05-12 18:34:00 +01:00
parent 3aec2ad4ac
commit e3dbff06e6
1 changed files with 9 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

13
applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/pEqn.H
View file

@ -1,3 +1,5 @@
rho = thermo.rho();
volScalarField rAU(1.0/UEqn.A()); volScalarField rAU(1.0/UEqn.A());
surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho*rAU)); surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho*rAU));
volVectorField HbyA(constrainHbyA(rAU*UEqn.H(), U, p)); volVectorField HbyA(constrainHbyA(rAU*UEqn.H(), U, p));
@ -27,7 +29,8 @@ if (pimple.transonic())
"phid", "phid",
(fvc::interpolate(psi)/fvc::interpolate(rho))*phiHbyA (fvc::interpolate(psi)/fvc::interpolate(rho))*phiHbyA
); );
phiHbyA -= fvc::interpolate(p)*phid;
phiHbyA -= fvc::interpolate(psi*p)*phiHbyA/fvc::interpolate(rho);
while (pimple.correctNonOrthogonal()) while (pimple.correctNonOrthogonal())
{ {
@ -86,9 +89,11 @@ U.correctBoundaryConditions();
fvOptions.correct(U); fvOptions.correct(U);
K = 0.5*magSqr(U); K = 0.5*magSqr(U);
pressureControl.limit(p); if (pressureControl.limit(p))
p.correctBoundaryConditions(); {
rho = thermo.rho(); p.correctBoundaryConditions();
rho = thermo.rho();
}
if (!pimple.transonic()) if (!pimple.transonic())
{ {