37 #ifndef averageNeighbourFvGeometryScheme_H 38 #define averageNeighbourFvGeometryScheme_H 86 const scalar minRatio,
tmp< fvMatrix< Type > > correction(const fvMatrix< Type > &)
Return the correction form of the given matrix by subtracting the matrix multiplied by the current fi...
A list of keyword definitions, which are a keyword followed by a number of values (eg...
virtual ~averageNeighbourFvGeometryScheme()=default
Destructor.
tmp< pointField > averageCentres(const pointField &cellCentres, const pointField &faceCentres, const vectorField &faceNormals) const
TypeName("averageNeighbour")
Runtime type information.
Geometry calculation scheme to minimise non-orthogonality/.
const label nIters_
Number of averaging iterations.
tmp< pointField > averageNeighbourCentres(const pointField &cellCentres, const vectorField &faceNormals, const scalarField &faceWeights) const
Average neighbouring cell centres to minimise non-ortho.
const scalar relax_
Blending between old-iteration cell centres and current average.
const fvMesh & mesh() const
Return mesh reference.
label clipFaceTet(const scalar minRatio, const vectorField &faceCentres, const vectorField &faceNormals, vectorField &faceCorrection) const
Clip face-centre correction vector if new triangle area would get below min. Return number of clipped...
void makeNonOrthoWeights(const pointField &cellCentres, const vectorField &faceNormals, scalarField &cosAngles, scalarField &faceWeights) const
Make weights based on non-orthogonality.
Geometry calculation scheme with automatic stabilisation for high-aspect ratio cells.
label clipPyramids(const pointField &cellCentres, const vectorField &faceCentres, const vectorField &faceNormals, const scalarField &minOwnHeight, const scalarField &minNeiHeight, vectorField &correction) const
Clip correction vector if new pyramid height would get below min. Return number of clipped cells...
Mesh data needed to do the Finite Volume discretisation.
A class for managing temporary objects.
virtual void movePoints()
Do what is necessary if the mesh has moved.
const scalar minRatio_
Clipping for pyramid heights - allowable shrinkage as fraction.
void makePyrHeights(const pointField &cellCentres, const vectorField &faceCentres, const vectorField &faceNormals, scalarField &ownHeight, scalarField &neiHeight) const
Calculate pyramid heights.