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());
89 const scalar magg =
mag(film.
g().
value());
93 bool dripping =
false;
95 if (sinAlpha[celli] > small && delta[celli] >
deltaStable_)
97 const scalar rhoc = rho[celli];
98 const scalar lc =
sqrt(sigma[celli]/(rhoc*magg));
99 const scalar deltaStable =
max 101 3*lc*
sqrt(1 -
sqr(sinAlpha[celli]))
106 if (delta[celli] > deltaStable)
108 const scalar ddelta =
max(delta[celli] - deltaStable, 0);
110 const scalar massDrip =
111 min(availableMass[celli],
max(ddelta*rhoc*magSf[celli], 0));
115 const scalar diam =
dCoeff_*lc;
118 massToEject[celli] += massDrip;
119 availableMass[celli] -= massDrip;
121 diameterToEject[celli] = diam;
131 diameterToEject[celli] = 0;
132 massToEject[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 > &)
BrunDrippingEjection(surfaceFilmRegionModel &film, const dictionary &dict)
Construct from surface film model.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Base class for film ejection models, handling mass transfer from the film.
scalarList diameter_
Diameters of particles to eject into the dripping.
dimensionedScalar sqrt(const dimensionedScalar &ds)
scalar ubarStar_
Critical non-dimensional interface velocity.
Macros for easy insertion into run-time selection tables.
const surfaceFilmRegionModel & film() const
Return const access to the film surface film model.
scalar dCoeff_
Coefficient relating the diameter of the drops formed to.
const dimensionedScalar sigma
Stefan-Boltzmann constant: default SI units: [W/m^2/K^4].
const Type & value() const
Return const reference to value.
const volVectorField & nHat() const
Return the patch normal vectors.
const fvMesh & regionMesh() const
Return the region mesh database.
const dimensionedVector & g() const
Return the acceleration due to gravity.
void addToEjectedMass(const scalar dMass)
Add to ejected mass.
virtual ~BrunDrippingEjection()
Destructor.
const volScalarField & delta() const
Return const access to the film thickness [m].
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
tmp< volScalarField > sigma() const
Return the surface tension coefficient [kg/s^2].
const volScalarField::Internal & magSf() const
Return the face area magnitudes [m^2].
scalar deltaStable_
Stable film thickness - drips only formed if thickness.
fileType type(const fileName &, const bool checkVariants=true, const bool followLink=true)
Return the file type: directory or file.
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
dimensioned< scalar > mag(const dimensioned< Type > &)
const volScalarField & rho() const
Return the film density [kg/m^3].
A class for managing temporary objects.
defineTypeNameAndDebug(kinematicSingleLayer, 0)
addToRunTimeSelectionTable(surfaceFilmModel, noFilm, mesh)