From 17afa7d79b09f7d2dbf71f55925c041e0f2bd32f Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Mon, 17 Feb 2020 14:04:45 +0000 Subject: [PATCH] reactingEulerFoam::AnisothermalPhaseModel: Added a continuity error compensation term to the internal energy pressure work Reduced the accumulation of error for incompressible and low compressibility cases. Partly resolves report https://bugs.openfoam.org/view.php?id=3442 --- .../AnisothermalPhaseModel/AnisothermalPhaseModel.C | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/AnisothermalPhaseModel/AnisothermalPhaseModel.C b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/AnisothermalPhaseModel/AnisothermalPhaseModel.C index 8cef201487..5959c0fafc 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/AnisothermalPhaseModel/AnisothermalPhaseModel.C +++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/AnisothermalPhaseModel/AnisothermalPhaseModel.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org - \\ / A nd | Copyright (C) 2015-2018 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2015-2020 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -134,7 +134,7 @@ Foam::AnisothermalPhaseModel::heEqn() tEEqn.ref() += filterPressureWork ( fvc::div(fvc::absolute(alphaPhi, alpha, U), this->thermo().p()) - + this->thermo().p()*fvc::ddt(alpha) + + (fvc::ddt(alpha) - contErr/this->rho())*this->thermo().p() ); } else if (this->thermo_->dpdt())