Removed unnecessary return statement
[foam-extend-3.2.git] / applications / solvers / multiphase / cavitatingFoam / UEqn.H
blob01911faaa3dbd1b23e49004c2ba1573226050510
1     surfaceScalarField muEff
2     (
3         "muEff",
4         twoPhaseProperties.muf()
5       + fvc::interpolate(rho*turbulence->nut())
6     );
8     fvVectorMatrix UEqn
9     (
10         fvm::ddt(rho, U)
11       + fvm::div(phi, U)
12       - fvm::laplacian(muEff, U)
13     //- (fvc::grad(U) & fvc::grad(muf))
14       - fvc::div(muEff*(fvc::interpolate(dev(fvc::grad(U))) & mesh.Sf()))
15     );
17     UEqn.relax();
19     if (momentumPredictor)
20     {
21         solve(UEqn == -fvc::grad(p));
22     }
24     Info<< "max(U) " << max(mag(U)).value() << endl;