Removed unnecessary return statement
[foam-extend-3.2.git] / applications / solvers / multiphase / compressibleInterFoam / UEqn.H
blob0b1a9ac029d776dbd57da20f91c14da8ee88f5a7
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(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     );
17     UEqn.relax();
19     if (momentumPredictor)
20     {
21         solve
22         (
23             UEqn
24          ==
25             fvc::reconstruct
26             (
27                 fvc::interpolate(rho)*(g & mesh.Sf())
28               + (
29                     fvc::interpolate(interface.sigmaK())*fvc::snGrad(alpha1)
30                   - fvc::snGrad(p)
31                 ) * mesh.magSf()
32             )
33         );
34     }