applications/solvers/coupled/MRFPorousFoam/couplingTerms.H

   1 {
   2     // Calculate grad p coupling matrix. Needs to be here if one uses
   3     // gradient schemes with limiters.  VV, 9/June/2014
   4     BlockLduSystem<vector, vector> pInU(fvm::grad(p));
   5
   6     // Calculate div U coupling.  Could be calculated only once since
   7     // it is only geometry dependent.  VV, 9/June/2014
   8     BlockLduSystem<vector, scalar> UInp(fvm::UDiv(U));
   9
  10     // Last argument in insertBlockCoupling says if the column direction
  11     // should be incremented. This is needed for arbitrary positioning
  12     // of U and p in the system. This could be better. VV, 30/April/2014
  13     UpEqn.insertBlockCoupling(0, 3, pInU, true);
  14     UpEqn.insertBlockCoupling(3, 0, UInp, false);
  15 }