applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/TEqn.H

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