zhyang-dev/CFDEMcoupling-PFM
1
0
Fork 0
mirror of https://github.com/ParticulateFlow/CFDEMcoupling-PFM.git synced 2025-12-08 06:37:44 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

convert tabs to spaces

Browse Source
This commit is contained in:
danielque
2023-04-07 14:25:59 +02:00
parent 57f672e0e0
commit e4de205ef5
1 changed files with 20 additions and 20 deletions
Download Patch File Download Diff File Expand all files Collapse all files

40
applications/solvers/cfdemSolverPimple/pEqn.H
View File

@ -10,22 +10,22 @@ phi = voidfractionf*phiByVoidfraction;
surfaceScalarField phiHbyA surfaceScalarField phiHbyA
( (
"phiHbyA", "phiHbyA",
( (
(fvc::interpolate(HbyA) & mesh.Sf() ) (fvc::interpolate(HbyA) & mesh.Sf() )
+ phicForces //explicit contribution + phicForces //explicit contribution
+ rAUfvoidfraction*fvc::ddtCorr(U, phiByVoidfraction) //correction + rAUfvoidfraction*fvc::ddtCorr(U, phiByVoidfraction) //correction
) )
); );
if (modelType=="A") if (modelType=="A")
rAUvoidfraction = volScalarField("(voidfraction2|A(U))",rAU*voidfraction*voidfraction); rAUvoidfraction = volScalarField("(voidfraction2|A(U))",rAU*voidfraction*voidfraction);
// Update the fixedFluxPressure BCs to ensure flux consistency // Update the fixedFluxPressure BCs to ensure flux consistency
if (modelType=="A") if (modelType=="A")
{ {
volScalarField rUsed = rAU*voidfraction; volScalarField rUsed = rAU*voidfraction;
constrainPressure(p, U, phiHbyA, rUsed,MRF); constrainPressure(p, U, phiHbyA, rUsed,MRF);
} }
else constrainPressure(p, U, phiHbyA, rAU,MRF); else constrainPressure(p, U, phiHbyA, rAU,MRF);
@ -37,23 +37,23 @@ while (pimple.correctNonOrthogonal())
fvm::laplacian(rAUvoidfraction, p) == fvc::div(voidfractionf*phiHbyA) + particleCloud.ddtVoidfraction() fvm::laplacian(rAUvoidfraction, p) == fvc::div(voidfractionf*phiHbyA) + particleCloud.ddtVoidfraction()
); );
pEqn.setReference(pRefCell, pRefValue); pEqn.setReference(pRefCell, pRefValue);
pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter()))); pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter())));
if (pimple.finalNonOrthogonalIter()) if (pimple.finalNonOrthogonalIter())
{ {
phiByVoidfraction = phiHbyA - pEqn.flux()/voidfractionf; phiByVoidfraction = phiHbyA - pEqn.flux()/voidfractionf;
phi = voidfractionf*phiByVoidfraction; phi = voidfractionf*phiByVoidfraction;
#include "continuityErrorPhiPU.H" #include "continuityErrorPhiPU.H"
// Explicitly relax pressure for momentum corrector // Explicitly relax pressure for momentum corrector
p.relax(); p.relax();
U = fvc::reconstruct(phiHbyA) U = fvc::reconstruct(phiHbyA)
- rAU*fvc::reconstruct(pEqn.flux()/voidfractionf/rAUf); - rAU*fvc::reconstruct(pEqn.flux()/voidfractionf/rAUf);
U.correctBoundaryConditions(); U.correctBoundaryConditions();
fvOptions.correct(U); fvOptions.correct(U);
} }
} }