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

multiphaseEulerFoam: Corrected handling of K for energy equation

Browse Source 
Resolves bug-report https://bugs.openfoam.org/view.php?id=3634
This commit is contained in:
Henry Weller
2021-02-24 15:37:28 +00:00
parent 9c82a94e31
commit a060c63235
2 changed files with 16 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

21
applications/solvers/multiphase/multiphaseEulerFoam/phaseSystems/phaseModel/AnisothermalPhaseModel/AnisothermalPhaseModel.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2015-2020 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2015-2021 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -102,12 +102,20 @@ Foam::AnisothermalPhaseModel<BasePhaseModel>::heEqn()
{
const volScalarField& alpha = *this;
const volVectorField U(this->U());
const surfaceScalarField alphaPhi(this->alphaPhi());
const surfaceScalarField alphaRhoPhi(this->alphaRhoPhi());
const tmp<volVectorField> tU(this->U());
const volVectorField& U(tU());
const volScalarField contErr(this->continuityError());
const volScalarField K(this->K());
const tmp<surfaceScalarField> talphaPhi(this->alphaPhi());
const surfaceScalarField& alphaPhi(talphaPhi());
const tmp<surfaceScalarField> talphaRhoPhi(this->alphaRhoPhi());
const surfaceScalarField& alphaRhoPhi(talphaRhoPhi());
const tmp<volScalarField> tcontErr(this->continuityError());
const volScalarField& contErr(tcontErr());
tmp<volScalarField> tK(this->K());
const volScalarField& K(tK());
volScalarField& he = this->thermo_->he();
@ -119,6 +127,7 @@ Foam::AnisothermalPhaseModel<BasePhaseModel>::heEqn()
+ fvc::ddt(alpha, this->rho(), K) + fvc::div(alphaRhoPhi, K)
- contErr*K
+ this->divq(he)
==
alpha*this->Qdot()

3
applications/solvers/multiphase/multiphaseEulerFoam/phaseSystems/phaseModel/MovingPhaseModel/MovingPhaseModel.C
View File

@ -264,8 +264,7 @@ void Foam::MovingPhaseModel<BasePhaseModel>::correctKinematics()
if (K_.valid())
{
K_.clear();
K();
K_.ref() = 0.5*magSqr(this->U());
}
}