applications/solvers/multiphase/cavitatingFoam/rhoEqn.H

   1 {
   2     fvScalarMatrix rhoEqn
   3     (
   4         fvm::ddt(rho)
   5       + fvm::div(phiv, rho)
   6     );
   7
   8     rhoEqn.solve();
   9
  10     phi = rhoEqn.flux();
  11
  12     Info<< "max-min rho: " << max(rho).value()
  13         << " " << min(rho).value() << endl;
  14
  15     rho == max(rho, rhoMin);
  16 }