34 namespace regionModels
36 namespace surfaceFilmModels
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;
#define forAll(list, i)
Loop across all elements in list.
Base class for surface film models.
Kinematic form of single-cell layer surface film model.
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 dimensionedScalar sigma
Stefan-Boltzmann constant: default SI units: [W/m^2/K^4].
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.
addToRunTimeSelectionTable(surfaceFilmRegionModel, kinematicSingleLayer, mesh)
Macros for easy insertion into run-time selection tables.
const surfaceFilmRegionModel & film() const
Return const access to the film surface film model.
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.
const dimensionedVector & g() const
Return the acceleration due to gravity.
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 > &)
virtual ~BrunDrippingInjection()
Destructor.
virtual const volScalarField & magSf() const
Return the face area magnitudes / [m^2].
const volScalarField & sigma() const
Return const access to the surface tension [kg/s^2].
fileType type(const fileName &, const bool checkVariants=true, const bool followLink=true)
Return the file type: directory or file.
BrunDrippingInjection(surfaceFilmRegionModel &film, const dictionary &dict)
Construct from surface film model.
dimensioned< scalar > mag(const dimensioned< Type > &)
virtual const volScalarField & rho() const
Return the film density [kg/m^3].
A class for managing temporary objects.
defineTypeNameAndDebug(kinematicSingleLayer, 0)