ENH: patchCloud: return pTraits<Type>::max for unfound points
[OpenFOAM-1.7.x.git] / applications / solvers / multiphase / multiphaseInterFoam / MRFMultiphaseInterFoam / UEqn.H
blobb47e7111922fdca962813649538ff88c5e633789
1     surfaceScalarField muEff
2     (
3         "muEff",
4         mixture.muf()
5       + fvc::interpolate(rho*turbulence->nut())
6     );
8     fvVectorMatrix UEqn
9     (
10         fvm::ddt(rho, U)
11       + fvm::div(mixture.rhoPhi(), U)
12       - fvm::laplacian(muEff, U)
13       - (fvc::grad(U) & fvc::grad(muEff))
14     //- fvc::div(muEff*(fvc::interpolate(dev(fvc::grad(U))) & mesh.Sf()))
15     );
16     mrfZones.addCoriolis(rho, UEqn);
18     UEqn.relax();
20     if (momentumPredictor)
21     {
22         solve
23         (
24             UEqn
25          ==
26             fvc::reconstruct
27             (
28                 (
29                     mixture.surfaceTensionForce()
30                   - ghf*fvc::snGrad(rho)
31                   - fvc::snGrad(p_rgh)
32                 ) * mesh.magSf()
33             )
34         );
35     }