applications/solvers/heatTransfer/buoyantBoussinesqPisoFoam/TEqn.H

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