applications/solvers/coupled/conjugateHeatSimpleFoam/solveEnergy.H

   1 {
   2     // Solid side
   3 #   include "readSolidControls.H"
   4
   5     for (int nonOrth = 0; nonOrth <= nNonOrthCorr; nonOrth++)
   6     {
   7         coupledFvScalarMatrix TEqns(2);
   8
   9         fvScalarMatrix* TFluidEqn = new fvScalarMatrix
  10         (
  11             rho*Cp*
  12             (
  13                 fvm::div(phi, T)
  14               + fvm::SuSp(-fvc::div(phi), T)
  15             )
  16           - fvm::laplacian(kappaEff, T)
  17             ==
  18             radiation->Ru()
  19           - fvm::Sp(4.0*radiation->Rp()*pow3(T), T)
  20           + 3.0*radiation->Rp()*pow4(T)
  21         );
  22
  23         TFluidEqn->relax();
  24
  25         fvScalarMatrix* TSolidEqn = new fvScalarMatrix
  26         (
  27           - fvm::laplacian(kSolidf, Tsolid, "laplacian(k,T)")
  28           + fvm::SuSp(-solidThermo.S()/Tsolid, Tsolid)
  29         );
  30
  31         TSolidEqn->relax();
  32
  33         // Add fluid equation
  34         TEqns.set(0, TFluidEqn);
  35
  36         // Add solid equation
  37         TEqns.set(1, TSolidEqn);
  38
  39         TEqns.solve();
  40     }
  41 }