Go to the documentation of this file.
37 sizes[0] =
mesh.nCells();
38 sizes[1] =
mesh.nPoints();
54 volumeCell_(
mesh.V()),
59 forAll(this->mesh_.C(), celli)
62 polyMeshTetDecomposition::cellTetIndices(this->mesh_, celli);
67 const scalar v = tetIs.
tet(this->mesh_).
mag();
68 volumeDual_[triIs[0]] += v;
69 volumeDual_[triIs[1]] += v;
70 volumeDual_[triIs[2]] += v;
74 mesh.globalData().syncPointData
90 volumeCell_(am.volumeCell_),
91 volumeDual_(am.volumeDual_),
102 this->mesh_.globalData().syncPointData
123 dataCell_[tetIs.
cell()] +=
125 / (0.25*volumeCell_[tetIs.
cell()]);
127 for (label i = 0; i < 3; ++i)
129 dataDual_[triIs[i]] +=
131 / (0.25*volumeDual_[triIs[i]]);
163 const label celli(tetIs.
cell());
171 this->mesh_.points()[triIs[0]] - this->mesh_.C()[celli],
172 this->mesh_.points()[triIs[1]] - this->mesh_.C()[celli],
173 this->mesh_.points()[triIs[2]] - this->mesh_.C()[celli]
178 const vector t( -
T.T().x() -
T.T().y() -
T.T().z());
187 const Type
s(dataCell_[celli]);
189 return (
T & S) + (t*
s);
Base class for lagrangian averaging methods.
List< label > labelList
A List of labels.
Defines the attributes of an object for which implicit objectRegistry management is supported,...
tetPointRef tet(const polyMesh &mesh) const
Return the geometry corresponding to this tet.
A field of fields is a PtrList of fields with reference counting.
gmvFile<< "tracers "<< particles.size()<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().x()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().y()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().z()<< " ";}gmvFile<< nl;forAll(lagrangianScalarNames, i){ word name=lagrangianScalarNames[i];IOField< scalar > s(IOobject(name, runTime.timeName(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
label cell() const
Return the cell.
A class for managing temporary objects.
static constexpr const zero Zero
Global zero (0)
AveragingMethod< Type >::TypeGrad TypeGrad
Public typedefs.
void add(const barycentric &coordinates, const tetIndices &tetIs, const Type &value)
Member Functions.
#define forAll(list, i)
Loop across all elements in list.
tmp< GeometricField< Type, faPatchField, areaMesh > > average(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
Area-weighted average a edgeField creating a areaField.
dimensionedSphericalTensor inv(const dimensionedSphericalTensor &dt)
tmp< Field< Type > > primitiveField() const
Return an internal field of the average.
PtrList< coordinateSystem > coordinates(solidRegions.size())
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Dual(const IOobject &io, const dictionary &dict, const fvMesh &mesh)
Constructors.
triFace faceTriIs(const polyMesh &mesh, const bool warn=true) const
Return the indices corresponding to the tri on the face for.
Mesh data needed to do the Finite Volume discretisation.
Type interpolate(const barycentric &coordinates, const tetIndices &tetIs) const
Interpolate.
TypeGrad interpolateGrad(const barycentric &coordinates, const tetIndices &tetIs) const
Interpolate gradient.
Storage and named access for the indices of a tet which is part of the decomposition of a cell.
A triangular face using a FixedList of labels corresponding to mesh vertices.
Dual-mesh lagrangian averaging procedure.
void average()
Calculate the average.
scalar mag() const
Return volume.
outerProduct< vector, Type >::type TypeGrad
Protected typedefs.