diff --git a/applications/solvers/multiphase/compressibleLesInterFoam/alphaEqns.H b/applications/solvers/multiphase/compressibleLesInterFoam/alphaEqns.H index 12eea31127e74f8a38b9ce7f0eab662ff3b1d9f4..819cd0f538b6a87d0448618482caf75c52e3c793 100644 --- a/applications/solvers/multiphase/compressibleLesInterFoam/alphaEqns.H +++ b/applications/solvers/multiphase/compressibleLesInterFoam/alphaEqns.H @@ -28,7 +28,7 @@ ), // Divergence term is handled explicitly to be // consistent with the explicit transport solution - divU*min(alpha1, 1.0) + divU*min(alpha1, scalar(1)) ); forAll(dgdt, celli) diff --git a/applications/solvers/multiphase/compressibleLesInterFoam/createFields.H b/applications/solvers/multiphase/compressibleLesInterFoam/createFields.H index c6dd6c0f48a62b428e9748e8cfee5447652150a1..6fa5049de47edfa4053792ba6015990e627e1a36 100644 --- a/applications/solvers/multiphase/compressibleLesInterFoam/createFields.H +++ b/applications/solvers/multiphase/compressibleLesInterFoam/createFields.H @@ -140,7 +140,7 @@ ); volScalarField dgdt = - pos(alpha2)*fvc::div(phi)/max(alpha2, 0.0001); + pos(alpha2)*fvc::div(phi)/max(alpha2, scalar(0.0001)); // Construct interface from alpha1 distribution interfaceProperties interface(alpha1, U, twoPhaseProperties); diff --git a/applications/solvers/multiphase/compressibleLesInterFoam/pEqn.H b/applications/solvers/multiphase/compressibleLesInterFoam/pEqn.H index b794d4df2d57dfeb3e78e2eeef2016055a080d7a..ebf24498ade0bd1d1572d6055da8b5bccd369ffa 100644 --- a/applications/solvers/multiphase/compressibleLesInterFoam/pEqn.H +++ b/applications/solvers/multiphase/compressibleLesInterFoam/pEqn.H @@ -40,7 +40,10 @@ solve ( - (max(alpha1, 0.0)*(psi1/rho1) + max(alpha2, 0.0)*(psi2/rho2)) + ( + max(alpha1, scalar(0))*(psi1/rho1) + + max(alpha2, scalar(0))*(psi2/rho2) + ) *pdEqnComp() + pdEqnIncomp );