54 const GeometricField<Type, fvsPatchField, surfaceMesh>& ssf
59 const fvMesh&
mesh = ssf.mesh();
67 tmp<GeometricField<GradType, fvPatchField, volMesh>> treconField
69 new GeometricField<GradType, fvPatchField, volMesh>
73 "reconstruct("+ssf.name()+
')',
80 dimensioned<GradType>(ssf.dimensions()/
dimArea,
Zero),
85 Field<GradType>& rf = treconField();
89 label own = owner[facei];
90 label nei = neighbour[facei];
92 rf[own] += (Cf[facei] -
C[own])*ssf[facei];
93 rf[nei] -= (Cf[facei] -
C[nei])*ssf[facei];
96 const typename GeometricField<Type, fvsPatchField, surfaceMesh>::
97 Boundary& bsf = ssf.boundaryField();
101 const fvsPatchField<Type>& psf = bsf[patchi];
104 const vectorField& pCf = Cf.boundaryField()[patchi];
108 label own = pOwner[pFacei];
109 rf[own] += (pCf[pFacei] -
C[own])*psf[pFacei];
115 treconField().correctBoundaryConditions();
131 const tmp<GeometricField<Type, fvsPatchField, surfaceMesh>>& tssf
135 tmp<GeometricField<GradType, fvPatchField, volMesh>> tvf
typeOfRank< typenamepTraits< arg1 >::cmptType, direction(pTraits< arg1 >::rank)+direction(pTraits< arg2 >::rank)>::type type
Reconstruct volField from a face flux field.
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh > > reconstruct(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
GeometricField< vector, fvsPatchField, surfaceMesh > surfaceVectorField
GeometricField< vector, fvPatchField, volMesh > volVectorField
const dimensionSet dimArea(sqr(dimLength))
Field< vector > vectorField
Specialisation of Field<T> for vector.
static constexpr const zero Zero
Global zero (0)
UList< label > labelUList
A UList of labels.
#define forAll(list, i)
Loop across all elements in list.
static const char *const typeName
The type name used in ensight case files.