applications/solvers/coupled/conjugateHeatFoam/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::ddt(T)
  14               + fvm::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         //Done automatically in 1.6.1-ext
  24         //TFluidEqn->boundaryManipulate(T.boundaryField());
  25
  26         fvScalarMatrix* TSolidEqn = new fvScalarMatrix
  27         (
  28             fvm::ddt(rhoCpsolid, Tsolid)
  29           - fvm::laplacian(ksolidf, Tsolid, "laplacian(k,T)")
  30           + fvm::SuSp(-solidThermo.S()/Tsolid, Tsolid)
  31         );
  32
  33         //Done automatically in 1.6.1-ext
  34         //TSolidEqn->boundaryManipulate(Tsolid.boundaryField());
  35
  36
  37         // Add fluid equation
  38         TEqns.set(0, TFluidEqn);
  39
  40         // Add solid equation
  41         TEqns.set(1, TSolidEqn);
  42
  43         TEqns.solve();
  44     }
  45 }