38#ifndef Foam_FA_regionFaIBMList_H
39#define Foam_FA_regionFaIBMList_H
125 #include "regionFaIBMList.txx"
Info<< nl;Info<< "Write faMesh in vtk format:"<< nl;{ vtk::uindirectPatchWriter writer(aMesh.patch(), fileName(aMesh.time().globalPath()/vtkBaseFileName));writer.writeGeometry();globalIndex procAddr(aMesh.nFaces());labelList cellIDs;if(UPstream::master()) { cellIDs.resize(procAddr.totalSize());for(const labelRange &range :procAddr.ranges()) { auto slice=cellIDs.slice(range);slice=identity(range);} } writer.beginCellData(4);writer.writeProcIDs();writer.write("cellID", cellIDs);writer.write("area", aMesh.S().field());writer.write("normal", aMesh.faceAreaNormals());writer.beginPointData(1);writer.write("normal", aMesh.pointAreaNormals());Info<< " "<< writer.output().name()<< nl;}{ vtk::lineWriter writer(aMesh.points(), aMesh.edges(), fileName(aMesh.time().globalPath()/(vtkBaseFileName+"-edges")));writer.writeGeometry();writer.beginCellData(4);writer.writeProcIDs();{ Field< scalar > fld(faMeshTools::flattenEdgeField(aMesh.magLe(), true))
IOdictionary(const IOobject &io, const dictionary *fallback=nullptr)
Construct given an IOobject and optional fallback dictionary content.
constexpr PtrList() noexcept
Class to control time during OpenFOAM simulations that is also the top-level objectRegistry.
UPtrList(Detail::PtrListDetail< regionFaIBM > &&ptrs) noexcept
Finite area mesh (used for 2-D non-Euclidian finite area method) defined using a patch of faces on a ...
virtual ~regionFaIBMList()=default
Destructor.
areaVectorField & Ustar() noexcept
Ustar velocity.
const areaVectorField & Ustar() const noexcept
Ustar velocity.
areaVectorField Ustar_
Interpolated velocity.
regionFaIBMList(const Time &runTime, const faMesh &aMesh)
Constructor.
void updateMesh()
Update mesh.
void zeroFilter(Type &fld) const
Apply the zero filter to an area field.
bool correct(const areaVectorField &Us)
Correct the regionFaIBM model.
GeometricField< vector, faPatchField, areaMesh > areaVectorField