PDRFoam: Added betav to ddtCorr

Resolves bug-report https://bugs.openfoam.org/view.php?id=3014
This commit is contained in:
Henry Weller
2018-08-13 14:47:30 +01:00
parent 535d1a46d2
commit a2c8fc63d9
2 changed files with 3 additions and 3 deletions

View File

@ -146,7 +146,7 @@ void PDRkEpsilon::correct()
volScalarField GR(drag.Gk());
volScalarField LI
(C4_*(Lobs + dimensionedScalar("minLength", dimLength, vSmall)));
(C4_*(Lobs + dimensionedScalar("minLength", dimLength, rootVSmall)));
// Dissipation equation
tmp<fvScalarMatrix> epsEqn

View File

@ -11,7 +11,7 @@ if (pimple.transonic())
fvc::interpolate(psi)
*(
fvc::flux(HbyA)
+ fvc::interpolate(rho*rAU)*fvc::ddtCorr(rho, U, phi)
+ fvc::interpolate(betav*rho*rAU)*fvc::ddtCorr(rho, U, phi)
/fvc::interpolate(rho)
)
);
@ -42,7 +42,7 @@ else
"phiHbyA",
(
fvc::flux(rho*HbyA)
+ fvc::interpolate(rho*rAU)*fvc::ddtCorr(rho, U, phi)
+ fvc::interpolate(betav*rho*rAU)*fvc::ddtCorr(rho, U, phi)
)
);