FreeFOAM The Cross-Platform CFD Toolkit
cellCuts Class Reference

Description of cuts across cells. More...

#include <dynamicMesh/cellCuts.H>


Detailed Description

Description of cuts across cells.

Description of cut is given as list of vertices and list of edges to be cut (and position on edge).

Does some checking of correctness/non-overlapping of cuts. 2x2x2 refinement has to be done in three passes since cuts can not overlap (would make addressing too complicated)

Introduces concept of 'cut' which is either an existing vertex or a edge.

Input can either be

  1. list of cut vertices and list of cut edges. Constructs cell circumference walks ('cellLoops').
  1. list of cell circumference walks. Will filter them so they don't overlap.
  1. cellWalker and list of cells to refine (refineCell). Constructs cellLoops and does B. cellWalker is class which can cut a single cell using a plane through the cell centre and in a certain normal direction

CellCuts constructed from cellLoops (B, C) can have multiple cut-edges and/or cut-point as long as there is per face only one (or none) cut across a face, i.e. a face can only be split into two faces.

The information available after construction:

  • pointIsCut, edgeIsCut.
  • faceSplitCut : the cross-cut of a face.
  • cellLoops : the loop which cuts across a cell
  • cellAnchorPoints : per cell the vertices on one side which are considered the anchor points.

AnchorPoints: connected loops have to be oriented in the same way to be able to grow new internal faces out of the same bottom faces. (limitation of the mapping procedure). The loop is cellLoop is oriented such that the normal of it points towards the anchorPoints. This is the only routine which uses geometric information.

Limitations:

  • cut description is very precise. Hard to get right.
  • do anchor points need to be unique per cell? Very inefficient.
  • is orientation guaranteed to be correct? Uses geometric info so can go wrong on highly distorted meshes.
  • is memory inefficient. Probably need to use Maps instead of labelLists.
Source files

Definition at line 109 of file cellCuts.H.

+ Inheritance diagram for cellCuts:
+ Collaboration diagram for cellCuts:

List of all members.

Public Member Functions

 ClassName ("cellCuts")
 Runtime type information.
 cellCuts (const polyMesh &mesh, const labelList &cutCells, const labelList &meshVerts, const labelList &meshEdges, const scalarField &meshEdgeWeights)
 Construct from cells to cut and pattern of cuts.
 cellCuts (const polyMesh &mesh, const labelList &meshVerts, const labelList &meshEdges, const scalarField &meshEdgeWeights)
 Construct from pattern of cuts. Detect cells to cut.
 cellCuts (const polyMesh &mesh, const labelList &cellLabels, const labelListList &cellLoops, const List< scalarField > &cellEdgeWeights)
 Construct from complete cellLoops through specified cells.
 cellCuts (const polyMesh &mesh, const cellLooper &cellCutter, const List< refineCell > &refCells)
 Construct from list of cells to cut and direction to cut in.
 cellCuts (const polyMesh &mesh, const cellLooper &cellCutter, const labelList &cellLabels, const List< plane > &cutPlanes)
 Construct from list of cells to cut and plane to cut with and.
 cellCuts (const polyMesh &mesh, const boolList &pointIsCut, const boolList &edgeIsCut, const scalarField &edgeWeight, const Map< edge > &faceSplitCut, const labelListList &cellLoops, const label nLoops, const labelListList &cellAnchorPoints)
 Construct from components.
 ~cellCuts ()
void clearOut ()
 Clear out demand driven storage.
const boolListpointIsCut () const
 Is mesh point cut.
const boolListedgeIsCut () const
 Is edge cut.
const scalarFieldedgeWeight () const
 If edge is cut gives weight (ratio between start() and end())
const labelListListfaceCuts () const
 Cuts per existing face (includes those along edge of face)
const Map< edge > & faceSplitCut () const
 Gives for split face the two cuts that split the face into two.
const labelListListcellLoops () const
 For each cut cell the cut along the circumference.
label nLoops () const
 Number of valid cell loops.
const labelListListcellAnchorPoints () const
 For each cut cell the points on the 'anchor' side of the cell.
pointField loopPoints (const label cellI) const
 Returns coordinates of points on loop for given cell.
labelList nonAnchorPoints (const labelList &cellPoints, const labelList &anchorPoints, const labelList &loop) const
 Invert anchor point selection.
void flip (const label cellI)
 Flip loop for cellI. Updates anchor points as well.
void flipLoopOnly (const label cellI)
 Flip loop for cellI. Does not update anchors. Use with care.
void writeOBJ (Ostream &os, const pointField &loopPoints, label &vertI) const
 debugging:Write list of cuts to stream in OBJ format
void writeOBJ (Ostream &os) const
 debugging:Write all of cuts to stream in OBJ format
void writeCellOBJ (const fileName &dir, const label cellI) const
 debugging:Write edges of cell and loop
- Public Member Functions inherited from edgeVertex
 edgeVertex (const polyMesh &mesh)
 Construct from mesh.
const polyMeshmesh () const
bool isEdge (const label eVert) const
label getEdge (const label eVert) const
label getVertex (const label eVert) const
label vertToEVert (const label vertI) const
label edgeToEVert (const label edgeI) const
point coord (const label cut, const scalar weight) const
label cutPairToEdge (const label cut0, const label cut1) const
OstreamwriteCut (Ostream &os, const label cut, const scalar) const
 Write cut description to Ostream.
OstreamwriteCuts (Ostream &os, const labelList &, const scalarField &) const
 Write cut descriptions to Ostream.

Additional Inherited Members

- Static Public Member Functions inherited from edgeVertex
static void updateLabels (const labelList &map, List< refineCell > &)
 Update refine list from map. Used to update cell/face labels.
static void updateLabels (const labelList &map, Map< label > &)
 Update map from map. Used to update cell/face labels.
static void updateLabels (const labelList &map, labelHashSet &)
 Update map from map. Used to update cell/face labels.
static bool isEdge (const primitiveMesh &mesh, const label eVert)
 is eVert an edge?
static label getEdge (const primitiveMesh &mesh, const label eVert)
 convert eVert to edge label
static label getVertex (const primitiveMesh &mesh, const label eVert)
 convert eVert to vertex label
static label vertToEVert (const primitiveMesh &mesh, const label vertI)
 Convert pointI to eVert.
static label edgeToEVert (const primitiveMesh &mesh, const label edgeI)
 Convert edgeI to eVert.
static point coord (const primitiveMesh &, const label cut, const scalar weight)
 Return coordinate of cut (uses weight if edgeCut)
static label cutPairToEdge (const primitiveMesh &, const label cut0, const label cut1)
 Find mesh edge (or -1) between two cuts.

Constructor & Destructor Documentation

cellCuts ( const polyMesh mesh,
const labelList cutCells,
const labelList meshVerts,
const labelList meshEdges,
const scalarField meshEdgeWeights 
)

Construct from cells to cut and pattern of cuts.

Definition at line 2698 of file cellCuts.C.

References Foam::endl(), and Foam::Pout.

cellCuts ( const polyMesh mesh,
const labelList meshVerts,
const labelList meshEdges,
const scalarField meshEdgeWeights 
)

Construct from pattern of cuts. Detect cells to cut.

Definition at line 2744 of file cellCuts.C.

References Foam::endl(), Foam::identity(), primitiveMesh::nCells(), and Foam::Pout.

cellCuts ( const polyMesh mesh,
const labelList cellLabels,
const labelListList cellLoops,
const List< scalarField > &  cellEdgeWeights 
)

Construct from complete cellLoops through specified cells.

Checks for consistency. Constructs cut-cut addressing and cellAnchorPoints.

Definition at line 2788 of file cellCuts.C.

References Foam::endl(), and Foam::Pout.

cellCuts ( const polyMesh mesh,
const cellLooper cellCutter,
const List< refineCell > &  refCells 
)

Construct from list of cells to cut and direction to cut in.

(always through cell centre) and celllooper.

Definition at line 2833 of file cellCuts.C.

References Foam::endl(), and Foam::Pout.

cellCuts ( const polyMesh mesh,
const cellLooper cellCutter,
const labelList cellLabels,
const List< plane > &  cutPlanes 
)

Construct from list of cells to cut and plane to cut with and.

celllooper. (constructor above always cuts through cell centre)

Definition at line 2877 of file cellCuts.C.

References Foam::endl(), and Foam::Pout.

cellCuts ( const polyMesh mesh,
const boolList pointIsCut,
const boolList edgeIsCut,
const scalarField edgeWeight,
const Map< edge > &  faceSplitCut,
const labelListList cellLoops,
const label  nLoops,
const labelListList cellAnchorPoints 
)

Construct from components.

Definition at line 2924 of file cellCuts.C.

References Foam::endl(), and Foam::Pout.

~cellCuts ( )

Definition at line 2955 of file cellCuts.C.


Member Function Documentation

ClassName ( "cellCuts"  )

Runtime type information.

void clearOut ( )

Clear out demand driven storage.

Definition at line 2961 of file cellCuts.C.

References Foam::deleteDemandDrivenData().

const boolList& pointIsCut ( ) const
inline

Is mesh point cut.

Definition at line 536 of file cellCuts.H.

Referenced by meshCutAndRemove::setRefinement().

const boolList& edgeIsCut ( ) const
inline

Is edge cut.

Definition at line 542 of file cellCuts.H.

Referenced by meshCutAndRemove::setRefinement(), and meshCutter::setRefinement().

const scalarField& edgeWeight ( ) const
inline

If edge is cut gives weight (ratio between start() and end())

Definition at line 548 of file cellCuts.H.

Referenced by meshCutAndRemove::setRefinement(), and meshCutter::setRefinement().

const labelListList& faceCuts ( ) const
inline

Cuts per existing face (includes those along edge of face)

Cuts in no particular order

Definition at line 555 of file cellCuts.H.

const Map<edge>& faceSplitCut ( ) const
inline

Gives for split face the two cuts that split the face into two.

Definition at line 565 of file cellCuts.H.

Referenced by meshCutAndRemove::setRefinement(), and meshCutter::setRefinement().

const labelListList& cellLoops ( ) const
inline

For each cut cell the cut along the circumference.

Definition at line 571 of file cellCuts.H.

Referenced by meshCutAndRemove::setRefinement(), and meshCutter::setRefinement().

label nLoops ( ) const
inline

Number of valid cell loops.

Definition at line 577 of file cellCuts.H.

Referenced by refinementIterator::setRefinement(), meshCutAndRemove::setRefinement(), and meshCutter::setRefinement().

const labelListList& cellAnchorPoints ( ) const
inline

For each cut cell the points on the 'anchor' side of the cell.

Definition at line 583 of file cellCuts.H.

Referenced by meshCutAndRemove::setRefinement(), and meshCutter::setRefinement().

Foam::pointField loopPoints ( const label  cellI) const

Returns coordinates of points on loop for given cell.

Uses cellLoops_ and edgeWeight_

Definition at line 2969 of file cellCuts.C.

References forAll, and List< T >::size().

Foam::labelList nonAnchorPoints ( const labelList cellPoints,
const labelList anchorPoints,
const labelList loop 
) const

Invert anchor point selection.

Definition at line 1188 of file cellCuts.C.

References Foam::findIndex(), forAll, List< T >::setSize(), and List< T >::size().

void flip ( const label  cellI)

Flip loop for cellI. Updates anchor points as well.

Definition at line 2995 of file cellCuts.C.

References mesh, and Foam::reverse().

void flipLoopOnly ( const label  cellI)

Flip loop for cellI. Does not update anchors. Use with care.

(only if you're sure loop orientation is wrong)

Definition at line 3013 of file cellCuts.C.

References Foam::reverse().

void writeOBJ ( Ostream os,
const pointField loopPoints,
label &  vertI 
) const

debugging:Write list of cuts to stream in OBJ format

Definition at line 3022 of file cellCuts.C.

References Foam::endl(), forAll, Vector< Cmpt >::x(), Vector< Cmpt >::y(), and Vector< Cmpt >::z().

void writeOBJ ( Ostream os) const

debugging:Write all of cuts to stream in OBJ format

Definition at line 3048 of file cellCuts.C.

References forAll, and Foam::meshTools::writeOBJ().

void writeCellOBJ ( const fileName dir,
const label  cellI 
) const

debugging:Write edges of cell and loop

Definition at line 3059 of file cellCuts.C.

References Foam::meshTools::writeOBJ().


The documentation for this class was generated from the following files: