diff --git a/applications/solvers/dfLowMachFoam/EEqn.H b/applications/solvers/dfLowMachFoam/EEqn.H index d72cc335a..73b972cc3 100644 --- a/applications/solvers/dfLowMachFoam/EEqn.H +++ b/applications/solvers/dfLowMachFoam/EEqn.H @@ -3,29 +3,27 @@ fvScalarMatrix EEqn ( - turbName == "laminar" - ? - ( - fvm::ddt(rho, he) + mvConvection->fvmDiv(phi, he) - + fvc::ddt(rho, K) + fvc::div(phi, K) - - dpdt - - fvm::laplacian(turbulence->alphaEff(), he) - + diffAlphaD - == - fvc::div(hDiffCorrFlux) - ) - : - ( + fvm::ddt(rho, he) + mvConvection->fvmDiv(phi, he) + fvc::ddt(rho, K) + fvc::div(phi, K) - dpdt - - fvm::laplacian(turbulence->alphaEff(), he) - ) + == + ( + turbName == "laminar" + ? + ( + fvm::laplacian(turbulence->alpha(), he) + - diffAlphaD + + fvc::div(hDiffCorrFlux) + ) + : + ( + fvm::laplacian(turbulence->alphaEff(), he) + ) + ) ); EEqn.relax(); EEqn.solve("ha"); - - }