37 sizes[0] =
mesh.nCells();
38 sizes[1] =
mesh.nPoints();
54 volumeCell_(
mesh.V()),
68 volumeDual_[triIs[0]] += v;
69 volumeDual_[triIs[1]] += v;
70 volumeDual_[triIs[2]] += v;
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.
const fvMesh & mesh_
The mesh on which the averaging is to be done.
virtual void average()
Calculate the average.
Dual-mesh lagrangian averaging procedure.
tmp< Field< Type > > primitiveField() const
Return an internal field of the average.
TypeGrad interpolateGrad(const barycentric &coordinates, const tetIndices &tetIs) const
Interpolate gradient.
void average()
Calculate the average.
A field of fields is a PtrList of fields with reference counting.
Generic templated field type.
Defines the attributes of an object for which implicit objectRegistry management is supported,...
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Sums a given list of (at least two or more) fields and outputs the result into a new field,...
Mesh data needed to do the Finite Volume discretisation.
const volVectorField & C() const
Return cell centres as volVectorField.
void syncPointData(List< Type > &pointData, const CombineOp &cop, const TransformOp &top) const
Helper to synchronise coupled patch point data.
static List< tetIndices > cellTetIndices(const polyMesh &mesh, label cI)
Return the tet decomposition of the given cell, see.
const globalMeshData & globalData() const
Return parallel info.
bool interpolate() const noexcept
Same as isPointData()
Storage and named access for the indices of a tet which is part of the decomposition of a cell.
triFace faceTriIs(const polyMesh &mesh, const bool warn=true) const
Return the indices corresponding to the tri on the face for.
tetPointRef tet(const polyMesh &mesh) const
Return the geometry corresponding to this tet.
label cell() const noexcept
Return the cell index.
scalar mag() const
Return volume.
A class for managing temporary objects.
A triangular face using a FixedList of labels corresponding to mesh vertices.
type
Volume classification types.
PtrList< coordinateSystem > coordinates(solidRegions.size())
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))
IOobject io("surfaceFilmProperties", mesh.time().constant(), mesh, IOobject::READ_IF_PRESENT, IOobject::NO_WRITE, false)
List< label > labelList
A List of labels.
dimensionedSphericalTensor inv(const dimensionedSphericalTensor &dt)
#define forAll(list, i)
Loop across all elements in list.