diff --git a/applications/solvers/multiphase/interFoam/interMixingFoam/alphaEqns.H b/applications/solvers/multiphase/interFoam/interMixingFoam/alphaEqns.H index d720aa97f654aa5d4f835efc4fe9f128b8f59a86..3f75765dd730917fc023a73eb9243fa27bf8db3b 100644 --- a/applications/solvers/multiphase/interFoam/interMixingFoam/alphaEqns.H +++ b/applications/solvers/multiphase/interFoam/interMixingFoam/alphaEqns.H @@ -108,6 +108,11 @@ ); } + alphaPhi1 = alphaPhi1BD + lambda*alphaPhi1; + + // Reset allLambda to 1.0 + allLambda = 1.0; + // Create the complete flux for alpha2 surfaceScalarField alphaPhi2 ( @@ -172,7 +177,6 @@ } // Construct the limited fluxes - alphaPhi1 = alphaPhi1BD + lambda*alphaPhi1; alphaPhi2 = alphaPhi2BD + lambda*alphaPhi2; // Solve for alpha1