Go to the documentation of this file.
40 namespace regionModels
42 namespace 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_,
242 if (!faOptions_.optionList::size())
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);
421 volScalarField::Boundary& pw =
460 const scalar massSource,
461 const vector& momentumSource,
462 const scalar pressureSource,
463 const scalar energySource
466 massSource_.boundaryFieldRef()[patchi][facei] += massSource;
467 pnSource_.boundaryFieldRef()[patchi][facei] += pressureSource;
468 momentumSource_.boundaryFieldRef()[patchi][facei] += momentumSource;
word UName_
Name of the velocity field.
int debug
Static debugging option.
const dimensionSet dimPressure
Defines the attributes of an object for which implicit objectRegistry management is supported,...
edgeScalarField phif_
Film momentum flux.
Selector class for relaxation factors, solver type and solution.
static options & New(const fvPatch &p)
A class for handling words, derived from Foam::string.
dimensionedScalar deltaT() const
Return time step.
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
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 dimensionSet dimVelocity
const dimensionSet dimEnergy
virtual void postEvolveRegion()
Post-evolve film.
tmp< areaScalarField > alpha() const
Wet indicator using h0.
volScalarField pnSource_
Normal pressure by particles.
const edgeScalarField & phi2s() const
Access continuity flux.
tmp< areaVectorField > Up() const
Primary region velocity at film hight. Assume the film to be.
word UName() const
Name of the U field.
const polyBoundaryMesh & boundaryMesh() const
Return boundary mesh.
areaScalarField gn_
Normal gravity field.
Ostream & endl(Ostream &os)
Add newline and flush stream.
const Type & value() const
Return const reference to value.
dimensionedScalar pos0(const dimensionedScalar &ds)
T get(const word &keyword, enum keyType::option matchOpt=keyType::REGEX) const
tmp< areaVectorField > Us() const
Film surface film velocity.
T lookupOrDefault(const word &keyword, const T &deflt, enum keyType::option matchOpt=keyType::REGEX) const
const Internal & internalField() const
Return a const-reference to the dimensioned internal field.
virtual void preEvolveRegion()
Pre-evolve film.
dimensionedScalar h0_
Smallest numerical thickness.
liquidFilmBase(const word &modelType, const fvPatch &patch, const dictionary &dict)
Construct from type name and mesh and dict.
dimensioned< vector > dimensionedVector
Dimensioned vector obtained from generic dimensioned type.
const areaVectorField & Uf() const
Access const reference Uf.
virtual label size() const
Return size.
tmp< areaScalarField > pg() const
Map primary static pressure.
const uniformDimensionedVectorField & g() const
Gravity.
void reduce(const List< UPstream::commsStruct > &comms, T &Value, const BinaryOp &bop, const int tag, const label comm)
uniformDimensionedVectorField g_
Gravity.
const edgeScalarField & deltaCoeffs() const
Return reference to difference factors array.
messageStream Info
Information stream (stdout output on master, null elsewhere)
A patch is a list of labels that address the faces in the global face list.
const areaScalarField & gn() const
Access const reference gn.
const areaVectorField & faceAreaNormals() const
Return face area normals.
volScalarField massSource_
Mass.
autoPtr< volSurfaceMapping > vsmPtr_
Volume-to surface mapping.
A finiteVolume patch using a polyPatch and a fvBoundaryMesh.
volVectorField momentumSource_
Momentum.
defineRunTimeSelectionTable(liquidFilmBase, dictionary)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
virtual void postEvolveRegion()
Post-evolve region.
const Type & lookupObject(const word &name, const bool recursive=false) const
areaVectorField Uf_
Film velocity.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
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,...
label index() const
Return the index of this patch in the fvBoundaryMesh.
This boundary condition is not designed to be evaluated; it is assmued that the value is assigned via...
const regionFaModel & region() const
Access to this region.
const dimensionSet dimMass(1, 0, 0, 0, 0, 0, 0)
const Time & time() const
Return the reference to the time database.
virtual void preEvolveRegion()
Pre-evolve region.
word pName_
Name of the pressure field.
Internal::FieldType & primitiveFieldRef(const bool updateAccessTime=true)
Return a reference to the internal field.
const dimensionSet dimAcceleration
virtual void addSources(const label patchi, const label facei, const scalar massSource, const vector &momentumSource, const scalar pressureSource, const scalar energySource=0)
Add sources.
dimensionedScalar deltaWet_
Delta wet for sub-models.
tmp< vectorField > Uwall() const
Return wall velocity field.
defineTypeNameAndDebug(kinematicThinFilm, 0)
This boundary condition provides a velocity condition for cases with moving walls.
const dimensionedScalar & h0() const
Return h0.
Boundary & boundaryFieldRef(const bool updateAccessTime=true)
Return a reference to the boundary field.
dimensionedScalar sqrt(const dimensionedScalar &ds)
Base class for area region models.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
Foam::fa::options & faOptions_
faOptions
const dimensionedScalar e
Elementary charge.
const areaScalarField & h() const
Access const reference h.
static const word null
An empty word.
scalar pRef()
Access to pRef.
Foam::fa::options & faOptions()
Return faOptions.
const fvPatch & patch_
Reference to fvPatch.
scalar pRef_
Reference absolute pressure.
areaScalarField h_
Film hight.
virtual ~liquidFilmBase()
Destructor.
edgeScalarField phi2s_
Film height flux.
const Foam::word liquidFilmName("liquidFilm")
static tmp< GeometricField< Type, faePatchField, edgeMesh > > interpolate(const GeometricField< Type, faPatchField, areaMesh > &tvf, const edgeScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
static const gravity & New(const Time &runTime)
Construct on Time.
const Boundary & boundaryField() const
Return const-reference to the boundary field.
const dimensionSet dimless
Dimensionless.
tmp< areaVectorField > Uw() const
Wall velocity.
Base class for thermal 2D shells.
virtual scalar CourantNumber() const
Courant number evaluation.
const edgeScalarField & phif() const
Access to momentum flux.