zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

BUG: Fixing BC's fixed value for species in multicomponent system

Browse Source
This commit is contained in:
sergio
2021-09-01 12:24:36 -07:00
parent 851cf11c0a
commit a88d0e7181
1 changed files with 22 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

25
src/phaseSystemModels/multiphaseInter/phasesSystem/phaseModel/MultiComponentPhaseModel/MultiComponentPhaseModel.C
View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com \\ / A nd | www.openfoam.com
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
Copyright (C) 2017-2020 OpenCFD Ltd. Copyright (C) 2017-2021 OpenCFD Ltd.
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
This file is part of OpenFOAM. This file is part of OpenFOAM.
@ -90,7 +90,9 @@ MultiComponentPhaseModel
( (
IOobject::groupName("X" + species_[i], phaseName), IOobject::groupName("X" + species_[i], phaseName),
fluid.mesh().time().timeName(), fluid.mesh().time().timeName(),
fluid.mesh() fluid.mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE
), ),
Y()[i] Y()[i]
) )
@ -116,7 +118,7 @@ void Foam::MultiComponentPhaseModel<BasePhaseModel, phaseThermo>
{ {
const dimensionedScalar Wi const dimensionedScalar Wi
( (
"W", "Wi",
dimMass/dimMoles, dimMass/dimMoles,
thermo().composition().W(i) thermo().composition().W(i)
); );
@ -126,6 +128,23 @@ void Foam::MultiComponentPhaseModel<BasePhaseModel, phaseThermo>
X_[i] = W*Y()[i]/Wi; X_[i] = W*Y()[i]/Wi;
Xtotal += X_[i]; Xtotal += X_[i];
X_[i].correctBoundaryConditions(); X_[i].correctBoundaryConditions();
const volScalarField::Boundary& YBf = Y()[i].boundaryField();
forAll(YBf, patchi)
{
const fvPatchScalarField& YPf = YBf[patchi];
if (YPf.fixesValue())
{
scalarField& xbf = X_[i].boundaryFieldRef()[patchi];
const scalarField& ybf = Y()[i].boundaryField()[patchi];
const scalarField& Wbf = W.boundaryField()[patchi];
forAll(xbf, facei)
{
xbf[facei] = Wbf[facei]*ybf[facei]/Wi.value();
}
}
}
} }
} }
X_[inertIndex_] = 1.0 - Xtotal; X_[inertIndex_] = 1.0 - Xtotal;