00001 surfaceScalarField muEff 00002 ( 00003 "muEff", 00004 twoPhaseProperties.muf() 00005 + fvc::interpolate(rho*turbulence->nut()) 00006 ); 00007 00008 fvVectorMatrix UEqn 00009 ( 00010 fvm::ddt(rho, U) 00011 + fvm::div(rhoPhi, U) 00012 - fvm::laplacian(muEff, U) 00013 - (fvc::grad(U) & fvc::grad(muEff)) 00014 //- fvc::div(muEff*(fvc::interpolate(dev(fvc::grad(U))) & mesh.Sf())) 00015 ); 00016 00017 UEqn.relax(); 00018 00019 if (momentumPredictor) 00020 { 00021 solve 00022 ( 00023 UEqn 00024 == 00025 fvc::reconstruct 00026 ( 00027 ( 00028 fvc::interpolate(interface.sigmaK())*fvc::snGrad(alpha1) 00029 - ghf*fvc::snGrad(rho) 00030 - fvc::snGrad(p_rgh) 00031 ) * mesh.magSf() 00032 ) 00033 ); 00034 } 00035 00036 // ************************ vim: set sw=4 sts=4 et: ************************ //
