64 "cellDisplacement" + cmptName_,
72 cellMotionBoundaryTypes<scalar>(pointDisplacement_.boundaryField())
74 pointLocation_(nullptr),
81 coeffDict().
found(
"frozenPointsZone")
82 ? fvMesh_.pointZones().findZoneID(coeffDict().lookup(
"frozenPointsZone"))
97 Info<<
"displacementComponentLaplacianFvMotionSolver:" <<
nl
98 <<
" diffusivity : " << diffusivityPtr_().type() <<
nl
99 <<
" frozenPoints zone : " << frozenPointsZone_ <<
endl;
122 Info<<
"displacementComponentLaplacianFvMotionSolver :"
123 <<
" Read pointVectorField "
124 << pointLocation_().name()
125 <<
" to be used for boundary conditions on points."
127 <<
"Boundary conditions:"
128 << pointLocation_().boundaryField().types() <<
endl;
152 if (pointLocation_.valid())
156 Info<<
"displacementComponentLaplacianFvMotionSolver : applying "
157 <<
" boundary conditions on " << pointLocation_().name()
158 <<
" to new point location."
164 pointLocation_().primitiveFieldRef() = fvMesh_.points();
166 pointLocation_().primitiveFieldRef().replace
169 points0_ + pointDisplacement_.primitiveField()
172 pointLocation_().correctBoundaryConditions();
175 if (frozenPointsZone_ != -1)
177 const pointZone& pz = fvMesh_.pointZones()[frozenPointsZone_];
181 label pointi = pz[i];
183 pointLocation_()[pointi][cmpt_] = points0_[pointi];
199 points0_ + pointDisplacement_.primitiveField()
203 if (frozenPointsZone_ != -1)
205 const pointZone& pz = fvMesh_.pointZones()[frozenPointsZone_];
209 label pointi = pz[i];
211 curPoints[pointi][cmpt_] = points0_[pointi];
215 twoDCorrectPoints(curPoints);
226 movePoints(fvMesh_.points());
228 diffusivityPtr_->correct();
229 pointDisplacement_.boundaryFieldRef().updateCoeffs();
235 diffusivityPtr_->operator()(),
237 "laplacian(diffusivity,cellDisplacement)"
252 diffusivityPtr_.reset(
nullptr);
256 coeffDict().lookup(
"diffusivity")
270 diffusivityPtr_.reset(
nullptr);
274 coeffDict().lookup(
"diffusivity")
#define forAll(list, i)
Loop across all elements in list.
Macros for easy insertion into run-time selection tables.
static pointMesh & New(const polyMesh &mesh)
void replace(const direction, const UList< cmptType > &)
Replace a component field of the field.
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
Virtual base class for displacement motion solver.
virtual void topoChange(const polyTopoChangeMap &)
Update local data for topology changes.
virtual void mapMesh(const polyMeshMap &)
Update from another mesh using the given map.
A list of keyword definitions, which are a keyword followed by any number of values (e....
const word & name() const
Return const reference to name.
Mesh motion solver for an fvMesh. Based on solving the cell-centre Laplacian for the given component ...
virtual tmp< pointField > curPoints() const
Return point location obtained from the current motion field.
virtual void topoChange(const polyTopoChangeMap &)
Update corresponding to the given map.
displacementComponentLaplacianFvMotionSolver(const word &name, const polyMesh &, const dictionary &)
Construct from polyMesh and dictionary.
~displacementComponentLaplacianFvMotionSolver()
Destructor.
virtual void mapMesh(const polyMeshMap &)
Update from another mesh using the given map.
virtual void solve()
Solve for motion.
const Time & time() const
Return the top-level database.
Base class for fvMesh based motionSolvers.
const fvMesh & fvMesh_
The fvMesh to be moved.
Abstract base class for cell-centre mesh motion diffusivity.
static autoPtr< motionDiffusivity > New(const fvMesh &mesh, Istream &mdData)
Select null constructed.
Virtual base class for mesh motion solver.
A subset of mesh points. The labels of points in the zone can be obtained from the addressing() list.
Class containing mesh-to-mesh mapping information.
Mesh consisting of general polyhedral cells.
Class containing mesh-to-mesh mapping information after a change in polyMesh topology.
A class for managing temporary objects.
T & ref() const
Return non-const reference or generate a fatal error.
Templated form of IOobject providing type information for file reading and header type checking.
bool headerOk()
Read header (uses typeGlobalFile to find file) and check.
tmp< PointField< Type > > interpolate(const VolField< Type > &) const
Interpolate volField using inverse distance weighting.
A class for handling words, derived from string.
Calculate the matrix for the laplacian of the field.
autoPtr< CompressibleMomentumTransportModel > New(const volScalarField &rho, const volVectorField &U, const surfaceScalarField &phi, const viscosity &viscosity)
tmp< fvMatrix< Type > > laplacian(const VolField< Type > &vf, const word &name)
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Ostream & endl(Ostream &os)
Add newline and flush stream.
vectorField pointField
pointField is a vectorField.
defineTypeNameAndDebug(combustionModel, 0)
word name(const complex &)
Return a string representation of a complex.
fileType type(const fileName &, const bool checkVariants=true, const bool followLink=true)
Return the file type: directory or file.
SolverPerformance< Type > solve(fvMatrix< Type > &, const word &)
Solve returning the solution statistics given convergence tolerance.
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
conserve primitiveFieldRef()+