twoPhaseMixtureThermo: Corrected mixture molecular weight

Resolves bug-report https://bugs.openfoam.org/view.php?id=3645
This commit is contained in:
Henry Weller
2021-03-21 10:09:46 +00:00
parent 7dc7956327
commit 8ec5eac35e

View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org \\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2013-2020 OpenFOAM Foundation \\ / A nd | Copyright (C) 2013-2021 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -398,7 +398,7 @@ Foam::tmp<Foam::scalarField> Foam::twoPhaseMixtureThermo::Cpv
Foam::tmp<Foam::volScalarField> Foam::twoPhaseMixtureThermo::W() const Foam::tmp<Foam::volScalarField> Foam::twoPhaseMixtureThermo::W() const
{ {
return 1/(Alpha1_/thermo1_->W() + Alpha2_/thermo1_->W()); return 1/(Alpha1_/thermo1_->W() + Alpha2_/thermo2_->W());
} }
@ -410,7 +410,7 @@ Foam::tmp<Foam::scalarField> Foam::twoPhaseMixtureThermo::W
return return
1/( 1/(
Alpha1_.boundaryField()[patchi]/thermo1_->W(patchi) Alpha1_.boundaryField()[patchi]/thermo1_->W(patchi)
+ Alpha2_.boundaryField()[patchi]/thermo1_->W(patchi) + Alpha2_.boundaryField()[patchi]/thermo2_->W(patchi)
); );
} }