a884a44e47
Patch contributed by Tobias Holzmann Resolves patch request https://bugs.openfoam.org/view.php?id=2524
39 lines
810 B
C
39 lines
810 B
C
if (finalIter)
|
|
{
|
|
mesh.data::add("finalIteration", true);
|
|
}
|
|
|
|
{
|
|
for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
|
|
{
|
|
fvScalarMatrix hEqn
|
|
(
|
|
fvm::ddt(betav*rho, h)
|
|
- (
|
|
thermo.isotropic()
|
|
? fvm::laplacian(betav*thermo.alpha(), h, "laplacian(alpha,h)")
|
|
: fvm::laplacian(betav*taniAlpha(), h, "laplacian(alpha,h)")
|
|
)
|
|
==
|
|
fvOptions(rho, h)
|
|
);
|
|
|
|
hEqn.relax();
|
|
|
|
fvOptions.constrain(hEqn);
|
|
|
|
solvPerfh = hEqn.solve(mesh.solver(h.select(finalIter)));
|
|
|
|
fvOptions.correct(h);
|
|
}
|
|
}
|
|
|
|
thermo.correct();
|
|
|
|
Info<< "Min/max T:" << min(thermo.T()).value() << ' '
|
|
<< max(thermo.T()).value() << endl;
|
|
|
|
if (finalIter)
|
|
{
|
|
mesh.data::remove("finalIteration");
|
|
}
|