applications/solvers/incompressible/simpleFoam/windSimpleFoam/UEqn.H

   1     // Solve the Momentum equation
   2
   3     tmp<fvVectorMatrix> UEqn
   4     (
   5         fvm::div(phi, U)
   6       + turbulence->divDevReff(U)
   7     );
   8
   9     // Add resistance on the actuation disks
  10     actuationDisks.addSu(UEqn());
  11
  12     UEqn().relax();
  13
  14     solve(UEqn() == -fvc::grad(p));