Go to the documentation of this file.
35 #include "surfaceInterpolate.H"
45 namespace regionModels
47 namespace surfaceFilmModels
87 const scalar rMin = 1
e-6;
90 forAll(definedPatchRadii_, i)
92 label patchi = definedPatchRadii_[i].first();
93 scalar definedInvR1 = 1.0/
max(rMin, definedPatchRadii_[i].second());
98 const scalar rMax = 1e6;
101 if (
mag(invR1[i]) < 1/rMax)
130 label cellO = own[facei];
131 label cellN = nbr[facei];
133 if (
phi[facei] > phiMax[cellO])
135 phiMax[cellO] =
phi[facei];
136 cosAngle[cellO] = -gHat_ & nf[facei];
138 if (-
phi[facei] > phiMax[cellN])
140 phiMax[cellN] = -
phi[facei];
141 cosAngle[cellN] = -gHat_ & -nf[facei];
154 if (phip[i] > phiMax[celli])
156 phiMax[celli] = phip[i];
157 cosAngle[celli] = -gHat_ & nf[i];
205 mesh.time().timeName(),
211 zeroGradientFvPatchScalarField::typeName
214 volCosAngle.correctBoundaryConditions();
218 return max(
min(cosAngle, scalar(1)), scalar(-1));
224 curvatureSeparation::curvatureSeparation
232 deltaByR1Min_(coeffDict_.getOrDefault<scalar>(
"deltaByR1Min", 0)),
233 definedPatchRadii_(),
237 if (magG_ < ROOTVSMALL)
240 <<
"Acceleration due to gravity must be non-zero"
244 gHat_ = film.
g().
value()/magG_;
258 const label patchi = patchIDs[j];
260 if (!uniquePatchIDs.found(patchi))
262 const scalar radius = prIn[i].second();
265 uniquePatchIDs.insert(patchi);
270 definedPatchRadii_.
transfer(prData);
290 refCast<const kinematicSingleLayer>(this->film());
304 const scalar Fthreshold = 1
e-10;
309 if ((invR1[i] > 0) && (
delta[i]*invR1[i] > deltaByR1Min_))
311 scalar R1 = 1.0/(invR1[i] + ROOTVSMALL);
312 scalar R2 = R1 +
delta[i];
315 scalar Fi = -
delta[i]*
rho[i]*magSqrU[i]*72.0/60.0*invR1[i];
319 - 0.5*
rho[i]*magG_*invR1[i]*(
sqr(R1) -
sqr(R2))*cosAngle[i];
322 scalar Fs =
sigma[i]/R2;
324 Fnet[i] = Fi +
Fb + Fs;
326 if (Fnet[i] + Fthreshold < 0)
334 massToInject = separated*availableMass;
335 diameterToInject = separated*
delta;
336 availableMass -= separated*availableMass;
338 addToInjectedMass(
sum(separated*availableMass));
347 mesh.time().timeName(),
353 zeroGradientFvPatchScalarField::typeName
356 volFnet.correctBoundaryConditions();
int debug
Static debugging option.
Defines the attributes of an object for which implicit objectRegistry management is supported,...
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
virtual ~curvatureSeparation()
Destructor.
A polyBoundaryMesh is a polyPatch list with additional search methods and registered IO.
A class for managing temporary objects.
static constexpr const zero Zero
Global zero (0)
const dimensionSet dimVelocity
A 1D vector of objects of type <T> that resizes itself as necessary to accept the new objects.
virtual const volVectorField & U() const
Return the film velocity [m/s].
tmp< volScalarField > calcInvR1(const volVectorField &U) const
Calculate local (inverse) radius of curvature.
An abstract base class with a fat-interface to all derived classes covering all possible ways in whic...
Calculate the divergence of the given field.
Kinematic form of single-cell layer surface film model.
void append(const T &val)
Append an element at the end of the list.
const polyBoundaryMesh & boundaryMesh() const
Return boundary mesh.
const Type & value() const
Return const reference to value.
const dimensionSet dimForce
addToRunTimeSelectionTable(surfaceFilmRegionModel, kinematicSingleLayer, mesh)
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
wordList names() const
Return a list of patch names.
const dimensionedVector & g() const
Return the acceleration due to gravity.
#define forAll(list, i)
Loop across all elements in list.
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
An indexed form of CGAL::Triangulation_face_base_2<K> used to keep track of the vertices in the trian...
tmp< vectorField > nf() const
Return face normals.
Base class for surface film models.
const fvMesh & regionMesh() const
Return the region mesh database.
A finiteVolume patch using a polyPatch and a fvBoundaryMesh.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
virtual const surfaceScalarField & phi() const
Return the film flux [kg.m/s].
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Base class for film injection models, handling mass transfer from the film.
Macros for easy insertion into run-time selection tables.
Mesh data needed to do the Finite Volume discretisation.
errorManipArg< error, int > exit(error &err, const int errNo=1)
Internal::FieldType & primitiveFieldRef(const bool updateAccessTime=true)
Return a reference to the internal field.
const volScalarField & sigma() const
Return const access to the surface tension [kg/s2].
defineTypeNameAndDebug(kinematicSingleLayer, 0)
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
const fvPatch & patch() const
Return patch.
Calculate the gradient of the given field.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
fileName::Type type(const fileName &name, const bool followLink=true)
Return the file type: DIRECTORY or FILE, normally following symbolic links.
virtual const volVectorField & nHat() const
Return the patch normal vectors.
const volScalarField & delta() const
Return const access to the film thickness [m].
const dimensionedScalar e
Elementary charge.
dimensioned< Type > sum(const DimensionedField< Type, GeoMesh > &df)
void transfer(List< T > &list)
Transfer contents of the argument List into this.
virtual const labelUList & faceCells() const
Return faceCells.
#define forAllReverse(list, i)
Reverse loop across all elements in list.
A List with indirect addressing.
Operations on lists of strings.
labelList findIndices(const ListType &input, typename ListType::const_reference val, label start=0)
Linear search to find all occurrences of given element.
dimensionedScalar sigma("sigma", dimMass/sqr(dimTime), transportProperties)
A 2-tuple for storing two objects of dissimilar types. The container is similar in purpose to std::pa...
virtual const volScalarField & rho() const
Return the film density [kg/m3].
Smooth ATC in cells next to a set of patches supplied by type.
const dimensionSet dimless
Dimensionless.
tmp< scalarField > calcCosAngle(const surfaceScalarField &phi) const
Calculate the cosine of the angle between gravity vector and.