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: Added API for droplet boiling

Browse Source
This commit is contained in:
andy
2011-12-13 18:35:14 +00:00
parent 188ad97333
commit fe66afddbc
7 changed files with 24 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -240,6 +240,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::calc
dt,
cellI,
Res,
Prs,
Ts,
mus/rhos,
d0,

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

@ -329,6 +329,7 @@ void Foam::ReactingParcel<ParcelType>::calc
dt,
cellI,
Res,
Prs,
Ts,
mus/rhos,
d0,
@ -464,6 +465,7 @@ void Foam::ReactingParcel<ParcelType>::calcPhaseChange
const scalar dt,
const label cellI,
const scalar Re,
const scalar Pr,
const scalar Ts,
const scalar nus,
const scalar d,
@ -500,11 +502,14 @@ void Foam::ReactingParcel<ParcelType>::calcPhaseChange
dt,
cellI,
Re,
Pr,
d,
nus,
T,
Ts,
pc_,
this->Tc_,
YComponents,
dMassPC
);

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

@ -210,6 +210,7 @@ protected:
const scalar dt, // timestep
const label cellI, // owner cell
const scalar Re, // Reynolds number
const scalar Pr, // Prandtl number
const scalar Ts, // Surface temperature
const scalar nus, // Surface kinematic viscosity
const scalar d, // diameter

10
src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C
View File

@ -140,20 +140,20 @@ void Foam::LiquidEvaporation<CloudType>::calculate
const scalar dt,
const label cellI,
const scalar Re,
const scalar Pr,
const scalar d,
const scalar nu,
const scalar T,
const scalar Ts,
const scalar pc,
const scalar Tc,
const scalarField& Yl,
scalarField& dMassPC
) const
{
// construct carrier phase species volume fractions for cell, cellI
const scalarField Xc(calcXc(cellI));
// droplet surface area
const scalar A = pi*sqr(d);
// calculate mass transfer of each specie in liquid
forAll(activeLiquids_, i)
{
@ -180,7 +180,7 @@ void Foam::LiquidEvaporation<CloudType>::calculate
// mass transfer coefficient [m/s]
const scalar kc = Sh*Dab/(d + ROOTVSMALL);
// vapour concentration at droplet surface [kmol/m3] at film temperature
// vapour concentration at surface [kmol/m3] at film temperature
const scalar Cs = pSat/(specie::RR*Ts);
// vapour concentration in bulk gas [kmol/m3] at film temperature
@ -190,7 +190,7 @@ void Foam::LiquidEvaporation<CloudType>::calculate
const scalar Ni = max(kc*(Cs - Cinf), 0.0);
// mass transfer [kg]
dMassPC[lid] += Ni*A*liquids_.properties()[lid].W()*dt;
dMassPC[lid] += Ni*pi*sqr(d)*liquids_.properties()[lid].W()*dt;
}
}

3
src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.H
View File

@ -111,11 +111,14 @@ public:
const scalar dt,
const label cellI,
const scalar Re,
const scalar Pr,
const scalar d,
const scalar nu,
const scalar T,
const scalar Ts,
const scalar pc,
const scalar Tc,
const scalarField& Yl,
scalarField& dMassPC
) const;

6
src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/PhaseChangeModel/PhaseChangeModel.C
View File

@ -135,6 +135,9 @@ void Foam::PhaseChangeModel<CloudType>::calculate
const scalar,
const scalar,
const scalar,
const scalar,
const scalar,
const scalarField&,
scalarField&
) const
{
@ -150,6 +153,9 @@ void Foam::PhaseChangeModel<CloudType>::calculate
"const scalar, "
"const scalar, "
"const scalar, "
"const scalar, "
"const scalar, "
"const scalarField&,"
"scalarField&"
") const"
);

3
src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/PhaseChangeModel/PhaseChangeModel.H
View File

@ -158,11 +158,14 @@ public:
const scalar dt,
const label cellI,
const scalar Re,
const scalar Pr,
const scalar d,
const scalar nu,
const scalar T,
const scalar Ts,
const scalar pc,
const scalar Tc,
const scalarField& Yl,
scalarField& dMassPC
) const;