39namespace surfaceFilmModels
56 ubarStar_(coeffDict_.getOrDefault<scalar>(
"ubarStar", 1.62208)),
57 dCoeff_(coeffDict_.getOrDefault<scalar>(
"dCoeff", 3.3)),
58 deltaStable_(coeffDict_.getOrDefault<scalar>(
"deltaStable", 0)),
59 diameter_(film.regionMesh().nCells(), -1.0)
79 refCast<const kinematicSingleLayer>(this->
film());
93 bool dripping =
false;
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 massToInject[celli] += massDrip;
119 availableMass[celli] -= massDrip;
121 diameterToInject[celli] = diam;
131 diameterToInject[celli] = 0;
132 massToInject[celli] = 0;
Macros for easy insertion into run-time selection tables.
#define addToRunTimeSelectionTable(baseType, thisType, argNames)
Add to construction table with typeName as the key.
virtual void correct()
Solve the turbulence equations and correct the turbulence viscosity.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
const Type & value() const
Return const reference to value.
virtual const volScalarField & magSf() const
Return the face area magnitudes / [m2].
Film Dripping mass transfer model.
scalar dCoeff_
Coefficient relating the diameter of the drops formed to.
scalar deltaStable_
Stable film thickness - drips only formed if thickness.
virtual ~BrunDrippingInjection()
Destructor.
scalar ubarStar_
Critical non-dimensional interface velocity.
scalarList diameter_
Diameters of particles to inject into the dripping.
const surfaceFilmRegionModel & film() const
Return const access to the film surface film model.
Base class for film injection models, handling mass transfer from the film.
void addToInjectedMass(const scalar dMass)
Add to injected mass.
Kinematic form of single-cell layer surface film model.
Base class for surface film models.
virtual const volScalarField & sigma() const =0
Return the film surface tension [N/m].
virtual const volScalarField & rho() const =0
Return the film density [kg/m3].
const dimensionedVector & g() const
Return the acceleration due to gravity.
virtual const volScalarField & delta() const =0
Return the film thickness [m].
A class for managing temporary objects.
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
fileName::Type type(const fileName &name, const bool followLink=true)
Return the file type: DIRECTORY or FILE, normally following symbolic links.
dimensionedScalar sqrt(const dimensionedScalar &ds)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
#define forAll(list, i)
Loop across all elements in list.