applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/TEqn.H

   1 {
   2     volScalarField kappaEff
   3     (
   4         "kappaEff",
   5         turbulence->nu()/Pr + turbulence->nut()/Prt
   6     );
   7
   8     fvScalarMatrix TEqn
   9     (
  10         fvm::div(phi, T)
  11       - fvm::Sp(fvc::div(phi), T)
  12       - fvm::laplacian(kappaEff, T)
  13     );
  14
  15     TEqn.relax();
  16
  17     eqnResidual = TEqn.solve().initialResidual();
  18     maxResidual = max(eqnResidual, maxResidual);
  19
  20     rhok = 1.0 - beta*(T - TRef);
  21 }