34 namespace regionModels
36 namespace surfaceFilmModels
46 BrunDrippingInjection::BrunDrippingInjection
53 ubarStar_(coeffDict_.lookupOrDefault(
"ubarStar", 1.62208)),
54 dCoeff_(coeffDict_.lookupOrDefault(
"dCoeff", 3.3)),
55 deltaStable_(coeffDict_.lookupOrDefault(
"deltaStable", scalar(0))),
76 refCast<const kinematicSingleLayer>(this->
film());
86 const scalar magg =
mag(film.
g().
value());
90 bool dripping =
false;
92 if (sinAlpha[celli] > SMALL && delta[celli] >
deltaStable_)
94 const scalar rhoc = rho[celli];
95 const scalar lc =
sqrt(sigma[celli]/(rhoc*magg));
96 const scalar deltaStable =
max 98 3*lc*
sqrt(1 -
sqr(sinAlpha[celli]))
103 if (delta[celli] > deltaStable)
105 const scalar ddelta =
max(delta[celli] - deltaStable, 0);
107 const scalar massDrip =
108 min(availableMass[celli],
max(ddelta*rhoc*magSf[celli], 0));
112 const scalar diam =
dCoeff_*lc;
115 massToInject[celli] += massDrip;
116 availableMass[celli] -= massDrip;
118 diameterToInject[celli] = diam;
128 diameterToInject[celli] = 0;
129 massToInject[celli] = 0;
const surfaceFilmModel & film() const
Return const access to the film surface film model.
#define forAll(list, i)
Loop across all elements in list.
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 > &)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar sqrt(const dimensionedScalar &ds)
const dimensionedVector & g() const
Return the accleration due to gravity.
const dimensionedScalar sigma
Stefan-Boltzmann constant: default SI units: [W/m2/K4].
tmp< volScalarField > gNorm() const
Return the gravity normal-to-patch component contribution.
scalar deltaStable_
Stable film thickness - drips only formed if thickness.
scalar ubarStar_
Critical non-dimensional interface velocity.
scalar dCoeff_
Coefficient relating the diameter of the drops formed to.
Macros for easy insertion into run-time selection tables.
addToRunTimeSelectionTable(surfaceFilmModel, kinematicSingleLayer, mesh)
const Type & value() const
Return const reference to value.
const fvMesh & regionMesh() const
Return the region mesh database.
scalarList diameter_
Diameters of particles to inject into the dripping.
void addToInjectedMass(const scalar dMass)
Add to injected mass.
Base class for film injection models, handling mass transfer from the film.
const volScalarField & delta() const
Return const access to the film thickness [m].
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
fileName::Type type(const fileName &, const bool followLink=true)
Return the file type: DIRECTORY or FILE.
virtual ~BrunDrippingInjection()
Destructor.
virtual const volScalarField & magSf() const
Return the face area magnitudes / [m2].
const volScalarField & sigma() const
Return const access to the surface tension [kg/s2].
Base class for surface film models.
dimensioned< scalar > mag(const dimensioned< Type > &)
virtual const volScalarField & rho() const
Return the film density [kg/m3].
A class for managing temporary objects.
defineTypeNameAndDebug(kinematicSingleLayer, 0)