39inline void resizeMatrix(SquareMatrix<T>& mat,
const label
n)
49inline void resizeField
67 pdrBlock_(std::cref<PDRblock>(PDRblock::null()))
83 pdrBlock_ = std::cref<PDRblock>(pdrBlock);
92 const labelVector faceDims(cellDims + labelVector::one);
97 const label maxDim =
cmptMax(pdrBlock.sizes()) + 2;
100 resizeField(
surf, cellDims);
116 resizeField(
drag_s, cellDims);
117 resizeField(
drag_r, cellDims);
134 resizeMatrix(
c_drag, maxDim);
136 resizeField(
face_patch, faceDims, labelVector::uniform(-1));
Info<< nl<< "Wrote faMesh in vtk format: "<< writer.output().name()<< nl;}{ vtk::lineWriter writer(aMesh.points(), aMesh.edges(), fileName(aMesh.mesh().time().globalPath()/"finiteArea-edges"));writer.writeGeometry();writer.beginCellData(4);writer.writeProcIDs();{ Field< scalar > fld(faMeshTools::flattenEdgeField(aMesh.magLe(), true))
IjkField< vector > drag_r
Directional drag from round obstacles.
IjkField< vector > betai_inv1
IjkField< vector > along_block
IjkField< scalar > v_block
Volume blockage.
SquareMatrix< scalar > ac_lblock
IjkField< vector > area_block_r
Summed area blockage (directional) from round obstacles.
Vector< List< scalar > > overlap_1d
SquareMatrix< scalar > b_lblock
IjkField< vector > sub_count
Number of obstacles parallel to specified direction.
IjkField< Vector< bool > > hole_in_face
Face field for (directional) hole in face.
SquareMatrix< symmTensor2D > c_drag
Cell-centred drag.
IjkField< vector > area_block_s
Summed area blockage (directional) from sharp obstacles.
SquareMatrix< scalar > abperim
SquareMatrix< scalar > bc_lblock
IjkField< scalar > obs_size
Obstacle size in cell.
SquareMatrix< scalar > aboverlap
IjkField< Vector< bool > > dirn_block
A total directional blockage in the cell.
IjkField< scalar > obs_count
Number of obstacles in cell.
IjkField< vector > grating_count
IjkField< symmTensor > drag_s
Tensorial drag from sharp obstacles.
PDRarrays()
Construct null.
IjkField< vector > face_block
IjkField< labelVector > face_patch
Face field for (directional) for patch Id.
SquareMatrix< scalar > a_lblock
IjkField< scalar > surf
Surface area in cell.
SquareMatrix< scalar > c_count
void reset()
Reset to defaults.
Vector< label > labelVector
Vector of labels.
void cmptMax(FieldField< Field, typename FieldField< Field, Type >::cmptType > &cf, const FieldField< Field, Type > &f)
static constexpr const zero Zero
Global zero (0)