zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

denseParticleFoam: Use Dcf*phid rather than Fdf for consistency

Browse Source
This commit is contained in:
Henry Weller
2023-04-08 16:04:42 +01:00
parent 0c08c2888c
commit 59218933a3
2 changed files with 3 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
applications/solvers/lagrangian/denseParticleFoam/denseParticleFoam.C
View File

@ -164,13 +164,10 @@ int main(int argc, char *argv[])
// Face continuous-dispersed phase drag coefficient
const surfaceScalarField Dcf(fvc::interpolate(Dc));
// Face dispersed phase drag force
const surfaceScalarField Fdf(fvc::flux(Fd));
// Effective flux of the dispersed phase
const surfaceScalarField phid
(
Fdf/(Dcf + dimensionedScalar(Dc.dimensions(), small))
fvc::flux(Fd)/(Dcf + dimensionedScalar(Dc.dimensions(), small))
);
// Face buoyancy force

4
applications/solvers/lagrangian/denseParticleFoam/pEqn.H
View File

@ -3,7 +3,7 @@
const volScalarField r1ADUc(1/(1 + rAUc*Dc));
const surfaceScalarField rAUcf(fvc::interpolate(rAUc));
const surfaceScalarField r1ADUcf(1/(1 + rAUcf*fvc::interpolate(Dc)));
const surfaceScalarField r1ADUcf(1/(1 + rAUcf*Dcf));
const surfaceScalarField rADUcf("Dp", r1ADUcf*rAUcf);
volVectorField HbyA(constrainHbyA(rAUc*UcEqn.H(), Uc, p));
@ -27,7 +27,7 @@
fvc::makeAbsolute(phiHbyAD, Uc);
}
phiHbyAD += rADUcf*(Fgf + Fdf);
phiHbyAD += rADUcf*(Fgf + Dcf*phid);
// Update the pressure BCs to ensure flux consistency
constrainPressure(p, Uc, phiHbyAD, rADUcf);