hEqn.H

Go to the documentation of this file.

{
    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)
    );
 
    if (meshFluxCorr)
    {
        surfaceScalarField phihMesh
        (
            fvc::interpolate(betav*rho*h)*mesh.phi()
        );
 
       hEqn -= fvc::div(phihMesh);
    }
 
    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;
 
    radiation->correct();
}