36 namespace regionModels
38 namespace areaSurfaceFilmModels
47 BrunDrippingInjection,
54 BrunDrippingInjection::BrunDrippingInjection
63 coeffDict_.getCheckOrDefault<scalar>
70 dCoeff_(coeffDict_.getOrDefault<scalar>(
"dCoeff", 3.3)),
71 deltaStable_(coeffDict_.getOrDefault<scalar>(
"deltaStable", 0)),
72 diameter_(film.regionMesh().nFaces(), -1.0)
92 const scalar magg =
mag(
film().
g().value());
96 bool dripping =
false;
100 const scalar rhoc =
rho[facei];
101 const scalar lc =
sqrt(
sigma[facei]/(rhoc*magg));
102 const scalar deltaStable =
max 104 3*lc*
sqrt(1 -
sqr(sinAlpha[facei]))
109 if (
delta[facei] > deltaStable)
111 const scalar massDrip =
112 availableMass[facei]*(
delta[facei] - deltaStable);
116 const scalar diam =
dCoeff_*lc;
119 massToInject[facei] += massDrip;
120 availableMass[facei] -= massDrip;
122 diameterToInject[facei] = diam;
132 diameterToInject[facei] = 0;
133 massToInject[facei] = 0;
void addToInjectedMass(const scalar dMass)
Add to injected mass.
scalar deltaStable_
Stable film thickness - drips only formed if thickness exceeds this threshold value.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
const areaScalarField & h() const
Access const reference h.
A list of keyword definitions, which are a keyword followed by a number of values (eg...
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
virtual const areaScalarField & sigma() const =0
Access const reference sigma.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar sqrt(const dimensionedScalar &ds)
scalar dCoeff_
Coefficient relating the diameter of the drops formed to the capillary length.
defineTypeNameAndDebug(kinematicThinFilm, 0)
dimensionedScalar sigma("sigma", dimMass/sqr(dimTime), transportProperties)
Base class for film injection models, handling mass transfer from the film.
Macros for easy insertion into run-time selection tables.
#define forAll(list, i)
Loop across all elements in list.
fileName::Type type(const fileName &name, const bool followLink=true)
Return the file type: DIRECTORY or FILE, normally following symbolic links.
scalarList diameter_
Diameters of particles to inject into the dripping.
addToRunTimeSelectionTable(liquidFilmBase, kinematicThinFilm, dictionary)
scalar ubarStar_
Critical non-dimensional interface velocity.
const uniformDimensionedVectorField & g
const areaScalarField & gn() const
Access const reference gn.
virtual const areaScalarField & rho() const =0
Access const reference rho.
A class for managing temporary objects.
const liquidFilmBase & film() const
Return const access to the film surface film model.