57 template<
class Type,
class NVDweight>
100 if (edgeFlux_[
edge] > 0)
102 d.removeCollinear(
n[owner[
edge]]);
106 d.removeCollinear(
n[neighbour[
edge]]);
109 d *=
mesh.edgeInterpolation::lPN().internalField()[
edge];
114 d = vector::uniform(SMALL);
125 gradc[neighbour[
edge]],
135 if (bWeights[patchI].
coupled())
139 const scalarField& pEdgeFlux = edgeFlux_.boundaryField()[patchI];
157 mesh.areaCentres().boundaryField()[patchI].patchInternalField()
162 mesh.areaCentres().boundaryField()[patchI]
163 .patchNeighbourField()
168 mesh.faceAreaNormals().boundaryField()[patchI]
169 .patchInternalField()
174 mesh.faceAreaNormals().boundaryField()[patchI]
175 .patchNeighbourField()
180 mesh.edgeInterpolation::lPN().boundaryField()[patchI]
185 vector d(CN[edgeI] - CP[edgeI]);
187 if (pEdgeFlux[edgeI] > 0)
189 d.removeCollinear(nP[edgeI]);
193 d.removeCollinear(nN[edgeI]);
201 d = vector::uniform(SMALL);
219 return tWeightingFactors;
Finite area mesh (used for 2-D non-Euclidian finite area method) defined using a patch of faces on a ...
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
T & ref() const
Return non-const reference to the contents of a non-null managed pointer.
bool coupled(solutionDict.getOrDefault("coupledEnergyField", false))
tmp< areaScalarField > limiter(const areaScalarField &phi)
#define forAll(list, i)
Loop across all elements in list.
Calculate the gradient of the given field.
An edge is a list of two point labels. The functionality it provides supports the discretisation on a...
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
GeometricField< scalar, faePatchField, edgeMesh > edgeScalarField
virtual tmp< edgeScalarField > weights(const GeometricField< Type, faPatchField, areaMesh > &) const
Return the interpolation weighting factors.
A Vector of values with scalar precision, where scalar is float/double depending on the compilation f...
Internal::FieldType & primitiveFieldRef(const bool updateAccessTime=true)
Return a reference to the internal field values.
const dimensionedScalar c
Speed of light in a vacuum.
Boundary & boundaryFieldRef(const bool updateAccessTime=true)
Return a reference to the boundary field.
Field< vector > vectorField
Specialisation of Field<T> for vector.
A class for managing temporary objects.
const Boundary & boundaryField() const noexcept
Return const-reference to the boundary field.
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
tmp< GeometricField< typename outerProduct< vector, Type >::type, faPatchField, areaMesh >> grad(const GeometricField< Type, faePatchField, edgeMesh > &ssf)