Go to the documentation of this file.
36 namespace reconstruction
46 Foam::reconstruction::isoAlpha::isoAlpha
67 isoFaceTol_(modelDict().getOrDefault<scalar>(
"isoFaceTol", 1
e-8)),
68 surfCellTol_(modelDict().getOrDefault<scalar>(
"surfCellTol", 1
e-8)),
69 sIterIso_(mesh_, ap_, surfCellTol_)
81 if (uptodate && !forceUpdate)
90 if (ap_.size() != mesh_.
nPoints())
void clear()
Clear the addressed list, i.e. set the size to zero.
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> reconstruct(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
defineTypeNameAndDebug(isoAlpha, 0)
bool topoChanging() const
Is mesh topology changing.
static constexpr const zero Zero
Global zero (0)
A 1D vector of objects of type <T> that resizes itself as necessary to accept the new objects.
static const volPointInterpolation & New(const fvMesh &mesh, Args &&... args)
Get existing or create a new MeshObject.
volVectorField normal_
Interface area normals.
const volScalarField & alpha1
tmp< GeometricField< Type, pointPatchField, pointMesh > > interpolate(const GeometricField< Type, fvPatchField, volMesh > &) const
Interpolate volField using inverse distance weighting.
boolList interfaceCell_
Is interface cell?
bool isASurfaceCell(const scalar alpha1) const
Determine if a cell is a surface cell.
#define forAll(list, i)
Loop across all elements in list.
volVectorField centre_
Interface centres.
label nCells() const
Number of mesh cells.
Reconstructs an interface (centre and normal vectors) consisting of isosurfaces to match the internal...
const vector & surfaceArea() const
The area vector of cutting isosurface.
volScalarField & alpha1_
Reference to the VoF Field.
void resize(const label newSize)
Adjust allocated size of list.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Macros for easy insertion into run-time selection tables.
addToRunTimeSelectionTable(reconstructionSchemes, isoAlpha, components)
bool alreadyReconstructed(bool forceUpdate=true) const
Is the interface already up-to-date?
label cellStatus()
The cellStatus.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
const dimensionedScalar e
Elementary charge.
label nPoints() const
Number of mesh points.
const point & surfaceCentre() const
The centre of cutting isosurface.
Original code supplied by Henning Scheufler, DLR (2019)
virtual void reconstruct(bool forceUpdate=true)
Reconstructs the interface.
DynamicField< T, SizeMin > & append(const T &val)
Append an element at the end of the list.
label vofCutCell(const label celli, const scalar alpha1, const scalar tol, const label maxIter)
DynamicField< label > interfaceLabels_
List of cell labels that have an interface.