From c12573c79967a50b7db2439f898f3a1e6b7eaf41 Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Wed, 12 Oct 2016 18:43:18 +0100 Subject: [PATCH] twoPhaseEulerFoam, reactingTwoPhaseEulerFoam: Corrected support for implicitPhasePressure with nAlphaCorr > 1 Resolves bug-report http://bugs.openfoam.org/view.php?id=2290 --- .../reactingTwoPhaseEulerFoam/twoPhaseSystem/twoPhaseSystem.C | 2 +- .../twoPhaseEulerFoam/twoPhaseSystem/twoPhaseSystem.C | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/twoPhaseSystem/twoPhaseSystem.C b/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/twoPhaseSystem/twoPhaseSystem.C index 6c039598a..8115f05cd 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/twoPhaseSystem/twoPhaseSystem.C +++ b/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/twoPhaseSystem/twoPhaseSystem.C @@ -403,7 +403,7 @@ void Foam::twoPhaseSystem::solve() fvScalarMatrix alpha1Eqn ( fvm::ddt(alpha1) - fvc::ddt(alpha1) - - fvm::laplacian(alphaDbyA, alpha1, "bounded") + - fvm::laplacian(alphaDbyA(), alpha1, "bounded") ); alpha1Eqn.relax(); diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseSystem/twoPhaseSystem.C b/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseSystem/twoPhaseSystem.C index 011cfbcd5..7c2323a03 100644 --- a/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseSystem/twoPhaseSystem.C +++ b/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseSystem/twoPhaseSystem.C @@ -523,7 +523,7 @@ void Foam::twoPhaseSystem::solve() fvScalarMatrix alpha1Eqn ( fvm::ddt(alpha1) - fvc::ddt(alpha1) - - fvm::laplacian(alpha1alpha2f*pPrimeByA_(), alpha1, "bounded") + - fvm::laplacian(alpha1alpha2f()*pPrimeByA_(), alpha1, "bounded") ); alpha1Eqn.relax();