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

updates for phase change modelling

Browse Source
This commit is contained in:
andy
2009-05-15 13:14:59 +01:00
parent f27c3cff74
commit 3f11699fce
3 changed files with 31 additions and 23 deletions
Download Patch File Download Diff File Expand all files Collapse all files

20
src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C
View File

@ -108,18 +108,17 @@ Foam::scalar Foam::ReactingMultiphaseParcel<ParcelType>::updateMassFractions
{
scalarField& YMix = this->Y_;
scalar dMassGasTot = sum(dMassGas);
scalar dMassLiquidTot = sum(dMassLiquid);
scalar dMassSolidTot = sum(dMassSolid);
scalar massGas =
this->updateMassFraction(mass0*YMix[GAS], dMassGas, YGas_);
scalar massLiquid =
this->updateMassFraction(mass0*YMix[LIQ], dMassLiquid, YLiquid_);
scalar massSolid =
this->updateMassFraction(mass0*YMix[SLD], dMassSolid, YSolid_);
this->updateMassFraction(mass0*YMix[GAS], dMassGas, YGas_);
this->updateMassFraction(mass0*YMix[LIQ], dMassLiquid, YLiquid_);
this->updateMassFraction(mass0*YMix[SLD], dMassSolid, YSolid_);
scalar massNew = max(massGas + massLiquid + massSolid, ROOTVSMALL);
scalar massNew = mass0 - (dMassGasTot + dMassLiquidTot + dMassSolidTot);
YMix[GAS] = (mass0*YMix[GAS] - dMassGasTot)/massNew;
YMix[LIQ] = (mass0*YMix[LIQ] - dMassLiquidTot)/massNew;
YMix[GAS] = massGas/massNew;
YMix[LIQ] = massLiquid/massNew;
YMix[SLD] = 1.0 - YMix[GAS] - YMix[LIQ];
return massNew;
@ -189,6 +188,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::calc
d0,
T0,
U0,
mass0,
idL,
YMix[LIQ],
YLiquid_,

29
src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
View File

@ -67,19 +67,25 @@ void Foam::ReactingParcel<ParcelType>::updateCellQuantities
template<class ParcelType>
void Foam::ReactingParcel<ParcelType>::updateMassFraction
Foam::scalar Foam::ReactingParcel<ParcelType>::updateMassFraction
(
const scalar mass0,
const scalarField& dMass,
scalarField& Y
)
) const
{
scalar mass1 = mass0 + sum(dMass);
scalar mass1 = mass0 - sum(dMass);
forAll(Y, i)
// only update the mass fractions if the new particle mass is finite
if (mass1 > ROOTVSMALL)
{
Y[i] = (Y[i]*mass0 - dMass[i])/mass1;
forAll(Y, i)
{
Y[i] = (Y[i]*mass0 - dMass[i])/mass1;
}
}
return mass1;
}
@ -116,13 +122,10 @@ void Foam::ReactingParcel<ParcelType>::calc
// Return enthalpy source and calc mass transfer due to phase change
scalar ShPC =
calcPhaseChange(td, dt, cellI, d0, T0, U0, 0, 1.0, Y_, dMassPC);
calcPhaseChange(td, dt, cellI, d0, T0, U0, mass0, 0, 1.0, Y_, dMassPC);
// Update particle component mass fractions
updateMassFraction(mass0, dMassPC, Y_);
// Update mass
scalar mass1 = mass0 - sum(dMassPC);
// Update particle component mass and mass fractions
scalar mass1 = updateMassFraction(mass0, dMassPC, Y_);
// Heat transfer
@ -218,6 +221,7 @@ Foam::scalar Foam::ReactingParcel<ParcelType>::calcPhaseChange
const scalar d,
const scalar T,
const vector& U,
const scalar mass,
const label idPhase,
const scalar YPhase,
const scalarField& YComponents,
@ -248,6 +252,9 @@ Foam::scalar Foam::ReactingParcel<ParcelType>::calcPhaseChange
dMassPC
);
// Limit phase change mass by availability of each specie
dMassPC = min(mass*YPhase*YComponents, dMassPC);
scalar dMassTot = sum(dMassPC);
// Add to cumulative phase change mass

5
src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.H
View File

@ -196,6 +196,7 @@ protected:
const scalar d, // diameter
const scalar T, // temperature
const vector& U, // velocity
const scalar mass, // mass
const label idPhase, // id of phase involved in phase change
const scalar YPhase, // total mass fraction
const scalarField& YComponents, // component mass fractions
@ -203,12 +204,12 @@ protected:
);
//- Update mass fraction
void updateMassFraction
scalar updateMassFraction
(
const scalar mass0,
const scalarField& dMass,
scalarField& Y
);
) const;
public: