42namespace areaSurfaceFilmModels
57 const word& modelType,
66 this->
solution().subDict(
"PIMPLE").get<
bool>(
"momentumPredictor")
70 this->
solution().subDict(
"PIMPLE").get<label>(
"nOuterCorr")
72 nCorr_(this->
solution().subDict(
"PIMPLE").get<label>(
"nCorr")),
75 this->
solution().subDict(
"PIMPLE").get<label>(
"nFilmCorr")
86 pRef_(
dict.get<scalar>(
"pRef")),
130 "ppf_" + regionName_,
143 "phif_" + regionName_,
156 "phi2s_" + regionName_,
180 g_(meshObjects::gravity::
New(primaryMesh().time())),
236 faOptions_(
Foam::fa::options::
New(
p))
244 Info <<
"No finite area options present" <<
endl;
275 Info<<
"Max film Courant Number: " <<
CoNum
276 <<
" Film velocity magnitude: " <<
velMag <<
endl;
306 && isA<movingWallVelocityFvPatchVectorField>(pp)
310 refCast<const movingWallVelocityFvPatchVectorField>(pp);
320 Uw.primitiveFieldRef() = UsWall() - nHat*(UsWall() & nHat);
388 Up.primitiveFieldRef() =
vsmPtr_->mapToSurface(Ust);
391 Up.primitiveFieldRef() -= nHat*(
Up.primitiveField() & nHat);
460 const scalar massSource,
461 const vector& momentumSource,
462 const scalar pressureSource,
463 const scalar energySource
Internal::FieldType & primitiveFieldRef(const bool updateAccessTime=true)
Return a reference to the internal field.
Boundary & boundaryFieldRef(const bool updateAccessTime=true)
Return a reference to the boundary field.
const Boundary & boundaryField() const
Return const-reference to the boundary field.
const Internal & internalField() const
Return a const-reference to the dimensioned internal field.
Defines the attributes of an object for which implicit objectRegistry management is supported,...
dimensionedScalar deltaT() const
Return time step.
This boundary condition is not designed to be evaluated; it is assmued that the value is assigned via...
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.
const edgeScalarField & deltaCoeffs() const
Return reference to difference factors array.
const areaVectorField & faceAreaNormals() const
Return face area normals.
A finiteVolume patch using a polyPatch and a fvBoundaryMesh.
virtual label size() const
Return size.
label index() const
Return the index of this patch in the fvBoundaryMesh.
This boundary condition provides a velocity condition for cases with moving walls.
tmp< vectorField > Uwall() const
Return wall velocity field.
const Type & lookupObject(const word &name, const bool recursive=false) const
const polyBoundaryMesh & boundaryMesh() const
Return boundary mesh.
A patch is a list of labels that address the faces in the global face list.
virtual bool write(const bool valid=true) const
Write using setting from DB.
const edgeScalarField & phif() const
Access to momentum flux.
areaVectorField Uf_
Film velocity.
word pName_
Name of the pressure field.
areaScalarField h_
Film hight.
virtual void postEvolveRegion()
Post-evolve film.
word UName_
Name of the velocity field.
virtual void addSources(const label patchi, const label facei, const scalar massSource, const vector &momentumSource, const scalar pressureSource, const scalar energySource=0)
Add sources.
areaScalarField gn_
Normal gravity field.
tmp< areaScalarField > alpha() const
Wet indicator using h0.
const regionFaModel & region() const
Access to this region.
uniformDimensionedVectorField g_
Gravity.
edgeScalarField phif_
Film momentum flux.
Foam::fa::options & faOptions()
Return faOptions.
Foam::fa::options & faOptions_
faOptions
tmp< areaScalarField > pg() const
Map primary static pressure.
volScalarField pnSource_
Normal pressure by particles.
dimensionedScalar h0_
Smallest numerical thickness.
tmp< areaVectorField > Uw() const
Wall velocity.
const areaVectorField & Uf() const
Access const reference Uf.
virtual void preEvolveRegion()
Pre-evolve film.
const areaScalarField & h() const
Access const reference h.
tmp< areaVectorField > Up() const
Primary region velocity at film hight. Assume the film to be.
volScalarField massSource_
Mass.
const edgeScalarField & phi2s() const
Access continuity flux.
virtual ~liquidFilmBase()
Destructor.
scalar pRef()
Access to pRef.
word UName() const
Name of the U field.
const dimensionedScalar & h0() const
Return h0.
volVectorField momentumSource_
Momentum.
edgeScalarField phi2s_
Film height flux.
dimensionedScalar deltaWet_
Delta wet for sub-models.
virtual scalar CourantNumber() const
Courant number evaluation.
tmp< areaVectorField > Us() const
Film surface film velocity.
scalar pRef_
Reference absolute pressure.
const uniformDimensionedVectorField & g() const
Gravity.
const areaScalarField & gn() const
Access const reference gn.
Base class for area region models.
virtual void postEvolveRegion()
Post-evolve region.
const Time & time() const
Return the reference to the time database.
autoPtr< volSurfaceMapping > vsmPtr_
Volume-to surface mapping.
const fvMesh & primaryMesh() const
Return the reference to the primary mesh database.
virtual void preEvolveRegion()
Pre-evolve region.
const fvPatch & patch_
Reference to fvPatch.
const faMesh & regionMesh() const
Return the region mesh database.
Selector class for relaxation factors, solver type and solution.
A class for managing temporary objects.
A class for handling words, derived from Foam::string.
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
const Foam::word liquidFilmName("liquidFilm")
const dimensionSet dimPressure
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
const dimensionSet dimless
Dimensionless.
dimensionedScalar pos0(const dimensionedScalar &ds)
const dimensionSet dimEnergy
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
const dimensionSet dimVelocity
messageStream Info
Information stream (stdout output on master, null elsewhere)
Ostream & endl(Ostream &os)
Add newline and flush stream.
dimensionedScalar sqrt(const dimensionedScalar &ds)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
void reduce(const List< UPstream::commsStruct > &comms, T &value, const BinaryOp &bop, const int tag, const label comm)
dimensioned< vector > dimensionedVector
Dimensioned vector obtained from generic dimensioned type.
static constexpr const zero Zero
Global zero (0)
const dimensionSet dimAcceleration
bool interpolate(const vector &p1, const vector &p2, const vector &o, vector &n, scalar l)
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh > > &tdf1, const word &name, const dimensionSet &dimensions)
Global function forwards to reuseTmpDimensionedField::New.
const dimensionSet dimMass(1, 0, 0, 0, 0, 0, 0)
Calculate the second temporal derivative.
#define defineRunTimeSelectionTable(baseType, argNames)
Define run-time selection table.