mirror of
https://develop.openfoam.com/Development/openfoam.git
synced 2025-11-28 03:28:01 +00:00
77b03e2e0c
e.g. (fvc::interpolate(HbyA) & mesh.Sf()) -> fvc::flux(HbyA) This removes the need to create an intermediate face-vector field when computing fluxes which is more efficient, reduces the peak storage and improved cache coherency in addition to providing a simpler and cleaner API.
72 lines
1.7 KiB
C
72 lines
1.7 KiB
C
{
|
|
surfaceScalarField alphaPhi
|
|
(
|
|
IOobject
|
|
(
|
|
"alphaPhi",
|
|
runTime.timeName(),
|
|
mesh
|
|
),
|
|
mesh,
|
|
dimensionedScalar("0", phi.dimensions(), 0)
|
|
);
|
|
|
|
surfaceScalarField phir(fvc::flux(UdmModel.Udm()));
|
|
|
|
if (nAlphaSubCycles > 1)
|
|
{
|
|
dimensionedScalar totalDeltaT = runTime.deltaT();
|
|
surfaceScalarField alphaPhiSum
|
|
(
|
|
IOobject
|
|
(
|
|
"alphaPhiSum",
|
|
runTime.timeName(),
|
|
mesh
|
|
),
|
|
mesh,
|
|
dimensionedScalar("0", phi.dimensions(), 0)
|
|
);
|
|
|
|
for
|
|
(
|
|
subCycle<volScalarField> alphaSubCycle(alpha1, nAlphaSubCycles);
|
|
!(++alphaSubCycle).end();
|
|
)
|
|
{
|
|
#include "alphaEqn.H"
|
|
alphaPhiSum += (runTime.deltaT()/totalDeltaT)*alphaPhi;
|
|
}
|
|
|
|
alphaPhi = alphaPhiSum;
|
|
}
|
|
else
|
|
{
|
|
#include "alphaEqn.H"
|
|
}
|
|
|
|
// Apply the diffusion term separately to allow implicit solution
|
|
// and boundedness of the explicit advection
|
|
{
|
|
fvScalarMatrix alpha1Eqn
|
|
(
|
|
fvm::ddt(alpha1) - fvc::ddt(alpha1)
|
|
- fvm::laplacian(turbulence->nut(), alpha1)
|
|
);
|
|
|
|
alpha1Eqn.solve(mesh.solver("alpha1Diffusion"));
|
|
|
|
alphaPhi += alpha1Eqn.flux();
|
|
alpha2 = 1.0 - alpha1;
|
|
|
|
Info<< "Phase-1 volume fraction = "
|
|
<< alpha1.weightedAverage(mesh.Vsc()).value()
|
|
<< " Min(" << alpha1.name() << ") = " << min(alpha1).value()
|
|
<< " Max(" << alpha1.name() << ") = " << max(alpha1).value()
|
|
<< endl;
|
|
}
|
|
|
|
rhoPhi = alphaPhi*(rho1 - rho2) + phi*rho2;
|
|
rho = mixture.rho();
|
|
}
|