7859083246
The deprecated non-const tmp functionality is now on the compiler switch NON_CONST_TMP which can be enabled by adding -DNON_CONST_TMP to EXE_INC in the Make/options file. However, it is recommended to upgrade all code to the new safer tmp by using the '.ref()' member function rather than the non-const '()' dereference operator when non-const access to the temporary object is required. Please report any problems on Mantis. Henry G. Weller CFD Direct.
39 lines
798 B
C
39 lines
798 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);
|
|
|
|
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");
|
|
}
|