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

ENH: Changing normalization with YMix for the ReactingMultiPhaseParcel

Browse Source 
The parcel mass fractions transfer of ReactingMultiPhaseParcel was modified
in order to be consistent between processors.
This commit is contained in:
sergio
2020-12-15 14:49:02 -08:00
parent 373d88a4c4
commit dc0372858d
1 changed files with 16 additions and 22 deletions
Download Patch File Download Diff File Expand all files Collapse all files

38
src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcelIO.C
View File

@ -6,7 +6,7 @@
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
Copyright (C) 2011-2017 OpenFOAM Foundation Copyright (C) 2011-2017 OpenFOAM Foundation
Copyright (C) 2016-2019 OpenCFD Ltd. Copyright (C) 2016-2020 OpenCFD Ltd.
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
This file is part of OpenFOAM. This file is part of OpenFOAM.
@ -71,12 +71,6 @@ Foam::ReactingMultiphaseParcel<ParcelType>::ReactingMultiphaseParcel
YGas_.transfer(Yg); YGas_.transfer(Yg);
YLiquid_.transfer(Yl); YLiquid_.transfer(Yl);
YSolid_.transfer(Ys); YSolid_.transfer(Ys);
// scale the mass fractions
const scalarField& YMix = this->Y_;
YGas_ /= YMix[GAS] + ROOTVSMALL;
YLiquid_ /= YMix[LIQ] + ROOTVSMALL;
YSolid_ /= YMix[SLD] + ROOTVSMALL;
} }
is.check(FUNCTION_NAME); is.check(FUNCTION_NAME);
@ -136,7 +130,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::readFields
label i = 0; label i = 0;
for (ReactingMultiphaseParcel<ParcelType>& p : c) for (ReactingMultiphaseParcel<ParcelType>& p : c)
{ {
p.YGas_[j] = YGas[i]/(p.Y()[GAS] + ROOTVSMALL); p.YGas_[j] = YGas[i]/(max(p.Y()[GAS], SMALL));
++i; ++i;
} }
} }
@ -156,7 +150,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::readFields
label i = 0; label i = 0;
for (ReactingMultiphaseParcel<ParcelType>& p : c) for (ReactingMultiphaseParcel<ParcelType>& p : c)
{ {
p.YLiquid_[j] = YLiquid[i]/(p.Y()[LIQ] + ROOTVSMALL); p.YLiquid_[j] = YLiquid[i]/(max(p.Y()[LIQ], SMALL));
++i; ++i;
} }
} }
@ -176,7 +170,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::readFields
label i = 0; label i = 0;
for (ReactingMultiphaseParcel<ParcelType>& p : c) for (ReactingMultiphaseParcel<ParcelType>& p : c)
{ {
p.YSolid_[j] = YSolid[i]/(p.Y()[SLD] + ROOTVSMALL); p.YSolid_[j] = YSolid[i]/(max(p.Y()[SLD], SMALL));
++i; ++i;
} }
} }
@ -224,7 +218,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::writeFields
label i = 0; label i = 0;
for (const ReactingMultiphaseParcel<ParcelType>& p0 : c) for (const ReactingMultiphaseParcel<ParcelType>& p0 : c)
{ {
YGas[i] = p0.YGas()[j]*p0.Y()[GAS]; YGas[i] = p0.YGas()[j]*max(p0.Y()[GAS], SMALL);
++i; ++i;
} }
@ -248,7 +242,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::writeFields
label i = 0; label i = 0;
for (const ReactingMultiphaseParcel<ParcelType>& p0 : c) for (const ReactingMultiphaseParcel<ParcelType>& p0 : c)
{ {
YLiquid[i] = p0.YLiquid()[j]*p0.Y()[LIQ]; YLiquid[i] = p0.YLiquid()[j]*max(p0.Y()[LIQ], SMALL);
++i; ++i;
} }
@ -272,7 +266,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::writeFields
label i = 0; label i = 0;
for (const ReactingMultiphaseParcel<ParcelType>& p0 : c) for (const ReactingMultiphaseParcel<ParcelType>& p0 : c)
{ {
YSolid[i] = p0.YSolid()[j]*p0.Y()[SLD]; YSolid[i] = p0.YSolid()[j]*max(p0.Y()[SLD], SMALL);
++i; ++i;
} }
@ -358,7 +352,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::readObjects
label i = 0; label i = 0;
for (ReactingMultiphaseParcel<ParcelType>& p0 : c) for (ReactingMultiphaseParcel<ParcelType>& p0 : c)
{ {
p0.YGas()[j]*p0.Y()[GAS] = YGas[i]; p0.YGas()[j]*max(p0.Y()[GAS], SMALL) = YGas[i];
++i; ++i;
} }
} }
@ -373,7 +367,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::readObjects
label i = 0; label i = 0;
for (ReactingMultiphaseParcel<ParcelType>& p0 : c) for (ReactingMultiphaseParcel<ParcelType>& p0 : c)
{ {
p0.YLiquid()[j]*p0.Y()[LIQ] = YLiquid[i]; p0.YLiquid()[j]*max(p0.Y()[LIQ], SMALL) = YLiquid[i];
++i; ++i;
} }
} }
@ -388,7 +382,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::readObjects
label i = 0; label i = 0;
for (ReactingMultiphaseParcel<ParcelType>& p0 : c) for (ReactingMultiphaseParcel<ParcelType>& p0 : c)
{ {
p0.YSolid()[j]*p0.Y()[SLD] = YSolid[i]; p0.YSolid()[j]*max(p0.Y()[SLD], SMALL) = YSolid[i];
++i; ++i;
} }
} }
@ -424,7 +418,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::writeObjects
label i = 0; label i = 0;
for (const ReactingMultiphaseParcel<ParcelType>& p0 : c) for (const ReactingMultiphaseParcel<ParcelType>& p0 : c)
{ {
YGas[i] = p0.YGas()[j]*p0.Y()[GAS]; YGas[i] = p0.YGas()[j]*max(p0.Y()[GAS], SMALL);
++i; ++i;
} }
} }
@ -439,7 +433,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::writeObjects
label i = 0; label i = 0;
for (const ReactingMultiphaseParcel<ParcelType>& p0 : c) for (const ReactingMultiphaseParcel<ParcelType>& p0 : c)
{ {
YLiquid[i] = p0.YLiquid()[j]*p0.Y()[LIQ]; YLiquid[i] = p0.YLiquid()[j]*max(p0.Y()[LIQ], SMALL);
++i; ++i;
} }
} }
@ -454,7 +448,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::writeObjects
label i = 0; label i = 0;
for (const ReactingMultiphaseParcel<ParcelType>& p0 : c) for (const ReactingMultiphaseParcel<ParcelType>& p0 : c)
{ {
YSolid[i] = p0.YSolid()[j]*p0.Y()[SLD]; YSolid[i] = p0.YSolid()[j]*max(p0.Y()[SLD], SMALL);
++i; ++i;
} }
} }
@ -471,9 +465,9 @@ Foam::Ostream& Foam::operator<<
const ReactingMultiphaseParcel<ParcelType>& p const ReactingMultiphaseParcel<ParcelType>& p
) )
{ {
scalarField YGasLoc(p.YGas()*p.Y()[0]); scalarField YGasLoc(p.YGas());
scalarField YLiquidLoc(p.YLiquid()*p.Y()[1]); scalarField YLiquidLoc(p.YLiquid());
scalarField YSolidLoc(p.YSolid()*p.Y()[2]); scalarField YSolidLoc(p.YSolid());
if (os.format() == IOstream::ASCII) if (os.format() == IOstream::ASCII)
{ {