Skip to content
Snippets Groups Projects
Commit 3b2cd0b1 authored by Henry Weller's avatar Henry Weller Committed by Andrew Heather
Browse files

DPMDyMFoam, DPMDyMFoam: Corrected support for closed-domain simulations

Also added support for extrapolated pressure boundary conditions.
parent b39b90a9
Branches
Tags
No related merge requests found
{
volVectorField HbyA("HbyA", Uc);
HbyA = rAUc*UcEqn.H();
volVectorField HbyA(constrainHbyA(rAUc*UcEqn.H(), Uc, p));
surfaceScalarField phiHbyA
(
......@@ -8,7 +7,6 @@
(
fvc::flux(HbyA)
+ alphacf*rAUcf*fvc::ddtCorr(Uc, Ucf)
+ phicForces
)
);
......@@ -19,6 +17,8 @@
fvc::makeAbsolute(phiHbyA, Uc);
}
phiHbyA += phicForces;
// Update the pressure BCs to ensure flux consistency
constrainPressure(p, Uc, phiHbyA, rAUcf);
......
{
volVectorField HbyA("HbyA", Uc);
HbyA = rAUc*UcEqn.H();
volVectorField HbyA(constrainHbyA(rAUc*UcEqn.H(), Uc, p));
surfaceScalarField phiHbyA
(
......@@ -8,10 +7,16 @@
(
fvc::flux(HbyA)
+ alphacf*rAUcf*fvc::ddtCorr(Uc, phic)
+ phicForces
)
);
if (p.needReference())
{
adjustPhi(phiHbyA, Uc, p);
}
phiHbyA += phicForces;
// Update the pressure BCs to ensure flux consistency
constrainPressure(p, Uc, phiHbyA, rAUcf);
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment