34 namespace regionModels
36 namespace surfaceFilmModels
64 return 0.037*
pow(Re, 0.8)*
cbrt(Sc);
71 standardPhaseChange::standardPhaseChange
78 deltaMin_(
readScalar(coeffDict_.lookup(
"deltaMin"))),
80 TbFactor_(coeffDict_.lookupOrDefault<scalar>(
"TbFactor", 1.1))
105 const label vapId = thermo.carrierId(filmThermo.
name());
130 const scalar pc = pInf[celli];
133 const scalar Tb = filmThermo.
Tb(pc);
139 const scalar pSat = filmThermo.
pv(pc, Tloc);
142 const scalar hVap = filmThermo.
hl(pc, Tloc);
148 const scalar Cp = filmThermo.
Cp(pc, Tloc);
149 const scalar Tcorr =
max(0.0, T[celli] - Tb);
150 const scalar qCorr = limMass[celli]*Cp*(Tcorr);
156 const scalar rhoInfc = rhoInf[celli];
159 const scalar muInfc = muInf[celli];
162 const scalar Re = rhoInfc*
mag(dU[celli])*
L_/muInfc;
165 const scalar Wvap = thermo.carrier().W(vapId);
168 const scalar Wliq = filmThermo.
W();
171 const scalar Ys = Wliq*pSat/(Wliq*pSat + Wvap*(pc - pSat));
174 const scalar Dab = filmThermo.
D(pc, Tloc);
177 const scalar Sc = muInfc/(rhoInfc*(Dab + ROOTVSMALL));
180 const scalar
Sh = this->
Sh(Re, Sc);
183 const scalar hm = Sh*Dab/(
L_ + ROOTVSMALL);
186 dm = dt*magSf[celli]*rhoInfc*hm*(Ys - YInf[celli])/(1.0 - Ys);
189 dMass[celli] +=
min(limMass[celli],
max(dm, 0));
193 dEnergy[celli] += dm*hs[celli];
#define forAll(list, i)
Loop across all elements in list.
virtual const word & name() const =0
Return the specie name.
virtual void correctModel(const scalar dt, scalarField &availableMass, scalarField &dMass, scalarField &dEnergy)
Correct.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
A list of keyword definitions, which are a keyword followed by any number of values (e...
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
virtual scalar Cp(const scalar p, const scalar T) const =0
Return specific heat capacity [J/kg/K].
const scalar L_
Length scale / [m].
virtual scalar hl(const scalar p, const scalar T) const =0
Return latent heat [J/kg].
dimensionedScalar sqrt(const dimensionedScalar &ds)
virtual const volVectorField & Us() const
Return the film surface velocity [m/s].
const volScalarField & rhoPrimary() const
Density [kg/m3].
Macros for easy insertion into run-time selection tables.
const volScalarField & muPrimary() const
Viscosity [Pa.s].
const filmThermoModel & filmThermo() const
Film thermo.
addToRunTimeSelectionTable(surfaceFilmModel, kinematicSingleLayer, mesh)
const volVectorField & UPrimary() const
Velocity [m/s].
scalar Sh(const scalar Re, const scalar Sc) const
Return Sherwood number as a function of Reynolds and Schmidt numbers.
virtual ~standardPhaseChange()
Destructor.
const scalar TbFactor_
Boiling temperature factor / [].
dimensionedScalar cbrt(const dimensionedScalar &ds)
const PtrList< volScalarField > & YPrimary() const
Specie mass fractions [0-1].
Thermo package for (S)olids (L)iquids and (G)ases Takes reference to thermo package, and provides:
bool readScalar(const char *buf, doubleScalar &s)
Read whole of buf as a scalar. Return true if succesful.
virtual scalar D(const scalar p, const scalar T) const =0
Return diffusivity [m2/s].
virtual scalar Tb(const scalar p) const =0
Return boiling temperature [K].
const volScalarField & pPrimary() const
Pressure [Pa].
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
const volScalarField & delta() const
Return const access to the film thickness [m].
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
virtual scalar W() const =0
Return molecular weight [kg/kmol].
virtual const volScalarField & magSf() const
Return the face area magnitudes / [m2].
Base class for surface film models.
virtual const volScalarField & hs() const
Return the film sensible enthalpy [J/kg].
dimensioned< scalar > mag(const dimensioned< Type > &)
virtual const volScalarField & rho() const
Return the film density [kg/m3].
const scalar deltaMin_
Minimum film height for model to be active.
defineTypeNameAndDebug(kinematicSingleLayer, 0)
const SLGThermo & thermo() const
Return const reference to the SLGThermo object.
virtual scalar pv(const scalar p, const scalar T) const =0
Return vapour pressure [Pa].
virtual const volScalarField & T() const
Return the film mean temperature [K].