42 displacementComponentLaplacianFvMotionSolver,
50 Foam::displacementComponentLaplacianFvMotionSolver::
51 displacementComponentLaplacianFvMotionSolver
63 "cellDisplacement" + cmptName_,
73 pointDisplacement_.dimensions(),
76 cellMotionBoundaryTypes<scalar>(pointDisplacement_.boundaryField())
78 pointLocation_(
nullptr),
85 coeffDict().
found(
"frozenPointsZone")
86 ? fvMesh_.pointZones().findZoneID(coeffDict().
lookup(
"frozenPointsZone"))
92 coeffDict().lookupOrDefault
99 if (applyPointLocation)
108 fvMesh_.time().timeName(),
119 Info<<
"displacementComponentLaplacianFvMotionSolver :" 120 <<
" Read pointVectorField " 121 << pointLocation_().name()
122 <<
" to be used for boundary conditions on points." 124 <<
"Boundary conditions:" 125 << pointLocation_().boundaryField().types() <<
endl;
149 if (pointLocation_.valid())
153 Info<<
"displacementComponentLaplacianFvMotionSolver : applying " 154 <<
" boundary conditions on " << pointLocation_().name()
155 <<
" to new point location." 161 pointLocation_().primitiveFieldRef() = fvMesh_.points();
163 pointLocation_().primitiveFieldRef().replace
166 points0_ + pointDisplacement_.primitiveField()
169 pointLocation_().correctBoundaryConditions();
172 if (frozenPointsZone_ != -1)
174 const pointZone& pz = fvMesh_.pointZones()[frozenPointsZone_];
178 label pointi = pz[i];
180 pointLocation_()[pointi][cmpt_] = points0_[pointi];
196 points0_ + pointDisplacement_.primitiveField()
200 if (frozenPointsZone_ != -1)
202 const pointZone& pz = fvMesh_.pointZones()[frozenPointsZone_];
206 label pointi = pz[i];
208 curPoints[pointi][cmpt_] = points0_[pointi];
212 twoDCorrectPoints(curPoints);
223 movePoints(fvMesh_.points());
225 diffusivityPtr_->correct();
226 pointDisplacement_.boundaryFieldRef().updateCoeffs();
232 diffusivityPtr_->operator()(),
234 "laplacian(diffusivity,cellDisplacement)" 249 diffusivityPtr_.reset(
nullptr);
253 coeffDict().
lookup(
"diffusivity")
tmp< GeometricField< Type, pointPatchField, pointMesh > > interpolate(const GeometricField< Type, fvPatchField, volMesh > &) const
Interpolate volField using inverse distance weighting.
#define forAll(list, i)
Loop across all elements in list.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
virtual void updateMesh(const mapPolyMesh &)
Update local data for topology changes.
T & ref() const
Return non-const reference or generate a fatal error.
~displacementComponentLaplacianFvMotionSolver()
Destructor.
Ostream & endl(Ostream &os)
Add newline and flush stream.
A simple wrapper around bool so that it can be read as a word: true/false, on/off, yes/no, y/n, t/f, or none.
Calculate the matrix for the laplacian of the field.
static word timeName(const scalar, const int precision=precision_)
Return time name of given scalar time.
void replace(const direction, const UList< cmptType > &)
Replace a component field of the field.
virtual void updateMesh(const mapPolyMesh &)
Update topology.
Class containing mesh-to-mesh mapping information after a change in polyMesh topology.
conserve primitiveFieldRef()+
Macros for easy insertion into run-time selection tables.
static const pointMesh & New(const polyMesh &mesh)
IOdictionary is derived from dictionary and IOobject to give the dictionary automatic IO functionalit...
Base class for fvMesh based motionSolvers.
vectorField pointField
pointField is a vectorField.
stressControl lookup("compactNormalStress") >> compactNormalStress
SolverPerformance< Type > solve(fvMatrix< Type > &, const dictionary &)
Solve returning the solution statistics given convergence tolerance.
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
const Time & time() const
Return time.
defineTypeNameAndDebug(combustionModel, 0)
virtual void solve()
Solve for motion.
fileName::Type type(const fileName &, const bool followLink=true)
Return the file type: DIRECTORY or FILE.
A subset of mesh points. The labels of points in the zone can be obtained from the addressing() list...
tmp< fvMatrix< Type > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Mesh consisting of general polyhedral cells.
A class for managing temporary objects.
static autoPtr< motionDiffusivity > New(const fvMesh &mesh, Istream &mdData)
Select null constructed.
Virtual base class for displacement motion solver.
virtual tmp< pointField > curPoints() const
Return point location obtained from the current motion field.
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...