diff --git a/applications/solvers/multiphase/interFoam/alphaEqn.H b/applications/solvers/multiphase/interFoam/alphaEqn.H index 4f44b732b585bb8b0048f99917e5cdfbf7b5aaac..1f07a7b4b727d72d2073441f2bfa4934a89e10e9 100644 --- a/applications/solvers/multiphase/interFoam/alphaEqn.H +++ b/applications/solvers/multiphase/interFoam/alphaEqn.H @@ -65,16 +65,16 @@ } } - if (MULESCorr) - { - tmp<surfaceScalarField> phiCN(phi); + tmp<surfaceScalarField> phiCN(phi); - // Calculate the Crank-Nicolson off-centred volumetric flux - if (ocCoeff > 0) - { - phiCN = (1.0 - cnCoeff)*phi + cnCoeff*phi.oldTime(); - } + // Calculate the Crank-Nicolson off-centred volumetric flux + if (ocCoeff > 0) + { + phiCN = (1.0 - cnCoeff)*phi + cnCoeff*phi.oldTime(); + } + if (MULESCorr) + { fvScalarMatrix alpha1Eqn ( #ifdef LTSSOLVE