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

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<BasePhaseModel>::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())