Go to the documentation of this file.
38 namespace regionModels
40 namespace areaSurfaceFilmModels
49 void contactAngleForce::initialise()
56 contactAngleForce::contactAngleForce
64 Ccf_(coeffDict_.get<scalar>(
"Ccf")),
69 typeName +
":fContactForceMask",
99 typeName +
":contactForce",
101 film().primaryMesh(),
129 const label faceO = own[edgei];
130 const label faceN = nbr[edgei];
133 if ((talpha()[faceO] > 0.5) && (talpha()[faceN] < 0.5))
137 else if ((talpha()[faceO] < 0.5) && (talpha()[faceN] > 0.5))
142 if (facei != -1 && mask_[facei] > 0.5)
147 gradAlpha[facei]/(
mag(gradAlpha[facei]) + ROOTVSMALL)
152 Ccf_*
n*
sigma[facei]*(1 - cosTheta)/invDx/
rhof[facei];
164 label face0 = faces[edgei];
183 tfvm.
ref() += tForce;
Defines the attributes of an object for which implicit objectRegistry management is supported,...
faPatchField<Type> abstract base class. This class gives a fat-interface to all derived classes cover...
A class for handling words, derived from Foam::string.
faMatrix< vector > faVectorMatrix
const fvMesh & primaryMesh() const
Return the reference to the primary mesh database.
A class for managing temporary objects.
static constexpr const zero Zero
Global zero (0)
const Field< Type > & field() const
Return field.
const dimensionSet dimDensity
static word timeName(const scalar t, const int precision=precision_)
tmp< areaScalarField > alpha() const
Wet indicator using h0.
tmp< GeometricField< typename outerProduct< vector, Type >::type, faPatchField, areaMesh >> grad(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
Unit conversion functions.
Base class for film (stress-based) force models.
const dimensionSet dimForce
virtual bool writeTime() const
Flag to indicate when to write a property.
#define forAll(list, i)
Loop across all elements in list.
dimensioned< vector > dimensionedVector
Dimensioned vector obtained from generic dimensioned type.
const labelUList & edgeFaces() const
Return edge-face addressing.
const dimensionSet dimArea(sqr(dimLength))
const DimensionedField< scalar, areaMesh > & S() const
Return face areas.
const edgeScalarField & deltaCoeffs() const
Return reference to difference factors array.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
const faMesh & regionMesh() const
Return the region mesh database.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Macros for easy insertion into run-time selection tables.
const labelUList & neighbour() const
Internal face neighbour.
constexpr scalar degToRad(const scalar deg) noexcept
Conversion from degrees to radians.
const labelUList & owner() const
Internal face owner.
Internal::FieldType & primitiveFieldRef(const bool updateAccessTime=true)
Return a reference to the internal field.
virtual const areaScalarField & rho() const =0
Access const reference rho.
defineTypeNameAndDebug(kinematicThinFilm, 0)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
const Time & time() const
Return the top-level database.
static const Vector< scalar > zero
virtual const areaScalarField & sigma() const =0
Access const reference sigma.
dimensionedScalar sigma("sigma", dimMass/sqr(dimTime), transportProperties)
surfaceScalarField rhof(fvc::interpolate(rho, "div(phi,rho)"))
const liquidFilmBase & film() const
Return const access to the film surface film model.
const faPatch & patch() const
Return patch.
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
const dimensionSet dimless
Dimensionless.
dimensionedScalar cos(const dimensionedScalar &ds)