35 #include "surfaceInterpolate.H" 45 namespace regionModels
47 namespace surfaceFilmModels
84 scalarField& invR1 = tinvR1.ref().primitiveFieldRef();
87 const scalar rMin = 1
e-6;
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];
148 const fvPatch& pp = phip.patch();
149 const labelList& faceCells = pp.faceCells();
153 label celli = faceCells[i];
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
231 gradNHat_(fvc::
grad(film.nHat())),
232 deltaByR1Min_(coeffDict_.getOrDefault<scalar>(
"deltaByR1Min", 0)),
233 definedPatchRadii_(),
234 magG_(
mag(film.
g().value())),
237 if (
magG_ < ROOTVSMALL)
240 <<
"Acceleration due to gravity must be non-zero" 258 const label patchi = patchIDs[j];
260 if (!uniquePatchIDs.
found(patchi))
262 const scalar radius = prIn[i].second();
265 uniquePatchIDs.insert(patchi);
290 refCast<const kinematicSingleLayer>(this->
film());
304 const scalar Fthreshold = 1
e-10;
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];
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;
347 mesh.time().timeName(),
353 zeroGradientFvPatchScalarField::typeName
356 volFnet.correctBoundaryConditions();
const polyBoundaryMesh & boundaryMesh() const
Return boundary mesh.
List< Tuple2< label, scalar > > definedPatchRadii_
List of radii for patches - if patch not defined, radius.
fvsPatchField< scalar > fvsPatchScalarField
ITstream & lookup(const word &keyword, enum keyType::option matchOpt=keyType::REGEX) const
Find and return an entry data stream. FatalIOError if not found, or not a stream. ...
const Type & value() const noexcept
Return const reference to value.
void size(const label n)
Older name for setAddressableSize.
vector gHat_
Direction of gravity vector.
tmp< scalarField > calcCosAngle(const surfaceScalarField &phi) const
Calculate the cosine of the angle between gravity vector and.
Base class for surface film models.
errorManipArg< error, int > exit(error &err, const int errNo=1)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
Kinematic form of single-cell layer surface film model.
error FatalError
Error stream (stdout output on all processes), with additional 'FOAM FATAL ERROR' header text and sta...
A list of keyword definitions, which are a keyword followed by a number of values (eg...
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
A 2-tuple for storing two objects of dissimilar types. The container is similar in purpose to std::pa...
void append(const T &val)
Append an element at the end of the list.
A 1D array of objects of type <T>, where the size of the vector is known and used for subscript bound...
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Curvature film separation model.
bool found(const T &val, label pos=0) const
True if the value if found in the list.
labelList findIndices(const ListType &input, typename ListType::const_reference val, label start=0)
Linear search to find all occurrences of given element.
Operations on lists of strings.
const dimensionSet dimless
Dimensionless.
scalar magG_
Magnitude of gravity vector.
tmp< volScalarField > calcInvR1(const volVectorField &U) const
Calculate local (inverse) radius of curvature.
const dictionary coeffDict_
Coefficients dictionary.
dimensionedScalar sigma("sigma", dimMass/sqr(dimTime), transportProperties)
virtual const volScalarField & rho() const =0
Return the film density [kg/m3].
addToRunTimeSelectionTable(surfaceFilmRegionModel, kinematicSingleLayer, mesh)
Macros for easy insertion into run-time selection tables.
dimensioned< Type > sum(const DimensionedField< Type, GeoMesh > &df)
const surfaceFilmRegionModel & film() const
Return const access to the film surface film model.
#define forAll(list, i)
Loop across all elements in list.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
fileName::Type type(const fileName &name, const bool followLink=true)
Return the file type: DIRECTORY or FILE, normally following symbolic links.
UPtrList< const labelUList > faceCells() const
Return a list of faceCells for each patch.
const dimensionedScalar e
Elementary charge.
A 1D vector of objects of type <T> that resizes itself as necessary to accept the new objects...
Calculate the gradient of the given field.
wordList names() const
Return a list of patch names.
virtual const volScalarField & delta() const =0
Return the film thickness [m].
const fvMesh & regionMesh() const
Return the region mesh database.
virtual const volVectorField & U() const =0
Return the film velocity [m/s].
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
A polyBoundaryMesh is a polyPatch list with additional search methods and registered IO...
Calculate the divergence of the given field.
void addToInjectedMass(const scalar dMass)
Add to injected mass.
const dimensionedVector & g() const
Return the acceleration due to gravity.
const uniformDimensionedVectorField & g
const dimensionSet dimForce
int debug
Static debugging option.
Base class for film injection models, handling mass transfer from the film.
virtual ~curvatureSeparation()
Destructor.
Internal::FieldType & primitiveFieldRef(const bool updateAccessTime=true)
Return a reference to the internal field values.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
Mesh data needed to do the Finite Volume discretisation.
A List with indirect addressing. Like IndirectList but does not store addressing. ...
Field< vector > vectorField
Specialisation of Field<T> for vector.
#define forAllReverse(list, i)
Reverse loop across all elements in list.
List< label > labelList
A List of labels.
A class for managing temporary objects.
An indexed form of CGAL::Triangulation_face_base_2<K> used to keep track of the vertices in the trian...
volTensorField gradNHat_
Gradient of surface normals.
virtual const volScalarField & sigma() const =0
Return the film surface tension [N/m].
defineTypeNameAndDebug(kinematicSingleLayer, 0)
scalar deltaByR1Min_
Minimum gravity driven film thickness (non-dimensionalised delta/R1)
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
tmp< GeometricField< typename outerProduct< vector, Type >::type, faPatchField, areaMesh >> grad(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
static constexpr const zero Zero
Global zero (0)
const dimensionSet dimVelocity