zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

TwoResistanceHeatTransferPhaseSystem: Restored asymmetric stabilisation term

Browse Source 
See commit 7e799e9f for details.
This commit is contained in:
Will Bainbridge
2019-09-16 11:33:48 +01:00
parent df683520c8
commit fc2e2b54da
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
applications/solvers/multiphase/reactingEulerFoam/phaseSystems/PhaseSystems/TwoResistanceHeatTransferPhaseSystem/TwoResistanceHeatTransferPhaseSystem.C
View File

@ -479,10 +479,10 @@ heatTransfer() const
const volScalarField HEff(H1*H2/(H1 + H2));
*eqns[phase1.name()] +=
HEff*(phase2.thermo().T() - phase1.thermo().T())
+ HEff/Cpv1*he1 - fvm::Sp(HEff/Cpv1, he1);
+ H1/Cpv1*he1 - fvm::Sp(H1/Cpv1, he1);
*eqns[phase2.name()] +=
HEff*(phase1.thermo().T() - phase2.thermo().T())
+ HEff/Cpv2*he2 - fvm::Sp(HEff/Cpv2, he2);
+ H2/Cpv2*he2 - fvm::Sp(H2/Cpv2, he2);
}
return eqnsPtr;