69 const cell&
c = mesh_.cells()[cellI_];
71 bool fullyBelow =
true;
72 bool fullyAbove =
true;
74 label nFaceBelowInterface = 0;
77 for (
const label facei :
c)
79 const label faceStatus = cutFace_.calcSubFace(facei, cutValue_);
83 cutFaceCentres_.append(cutFace_.subFaceCentre());
84 cutFaceAreas_.append(cutFace_.subFaceArea());
85 isoFaceEdges_.append(cutFace_.surfacePoints());
89 else if (faceStatus == -1)
91 cutFaceCentres_.append(cutFace_.subFaceCentre());
92 cutFaceAreas_.append(cutFace_.subFaceArea());
94 nFaceBelowInterface++;
102 if (!fullyBelow && !fullyAbove)
118 if (
mag(faceArea_) < ROOTVSMALL)
120 if (nFaceBelowInterface == 0)
124 subCellCentre_ =
Zero;
133 subCellCentre_ = mesh_.C()[cellI_];
134 subCellVolume_ = mesh_.V()[cellI_];
140 cutFaceCentres_.append(faceCentre_);
141 cutFaceAreas_.append(faceArea_);
152 VOF_ = subCellVolume_ / mesh_.V()[cellI_];
157 subCellCentre_ =
Zero;
164 subCellCentre_ = mesh_.C()[cellI_];
165 subCellVolume_ = mesh_.V()[cellI_];
175 if (facePoints_.empty())
178 calcIsoFacePointsFromEdges
195 cutFaceCentres_.clear();
196 cutFaceAreas_.clear();
197 isoFaceEdges_.clear();
201 subCellCentre_ =
Zero;
202 subCellVolume_ = -10;
dimensioned< Type > average(const DimensionedField< Type, GeoMesh > &f1)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
cutCellIso(const fvMesh &mesh, scalarField &f)
Construct from fvMesh and a scalarField.
Class for cutting a face, faceI, of an fvMesh, mesh_, at its intersection with an isosurface defined ...
Service routines for cutting a cell, celli, of an fvMesh, mesh_, at its intersection with a surface...
A 1D vector of objects of type <T> that resizes itself as necessary to accept the new objects...
void clearStorage()
Resets internal values.
A cell is defined as a list of faces with extra functionality.
Mesh data needed to do the Finite Volume discretisation.
const dimensionedScalar c
Speed of light in a vacuum.
label calcSubCell(const label cellI, const scalar cutValue)
Sets internal values and returns face status.
const DynamicList< point > & facePoints()
Returns the points of the cutting PLICface.
static constexpr const zero Zero
Global zero (0)