29 #include "surfaceInterpolate.H" 46 template<
class Type,
class GType>
93 template<
class Type,
class GType>
109 "gammaSnGradCorr("+vf.
name()+
')',
123 tgammaSnGradCorr.
ref().replace
130 return tgammaSnGradCorr;
136 template<
class Type,
class GType>
150 tLaplacian.
ref().rename(
"laplacian(" + vf.
name() +
')');
156 template<
class Type,
class GType>
180 this->tsnGradScheme_().deltaCoeffs(vf),
185 tmp<SType> tfaceFluxCorrection = gammaSnGradCorr(SfGammaCorr, vf);
187 if (this->tsnGradScheme_().corrected())
189 tfaceFluxCorrection.
ref() +=
190 SfGammaSn*this->tsnGradScheme_().correction(vf);
193 const word corrName(tfaceFluxCorrection().
name());
197 const word oldName(corrName +
"_0");
198 const scalar
relax(vf.
mesh().equationRelaxationFactor(oldName));
201 if (obr.foundObject<SType>(oldName))
203 SType& oldCorrection = obr.lookupObjectRef<SType>(oldName);
205 trelaxedCorrection.ref() *=
relax;
206 trelaxedCorrection.ref() += (1.0-
relax)*oldCorrection;
208 oldCorrection = tfaceFluxCorrection;
212 SType*
s =
new SType(oldName, tfaceFluxCorrection);
221 )().primitiveField();
232 template<
class Type,
class GType>
254 SfGammaSn*this->tsnGradScheme_().
snGrad(vf)
255 + gammaSnGradCorr(SfGammaCorr, vf)
259 tLaplacian.
ref().rename
261 "laplacian(" +
gamma.name() +
',' + vf.
name() +
')' virtual bool coupled() const
True if the patch field is coupled.
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
const Internal::FieldType & primitiveField() const noexcept
Return a const-reference to the internal field values.
static scalar Sn(const scalar a, const scalar x)
const word & name() const noexcept
Return the object name.
T & ref() const
Return non-const reference to the contents of a non-null managed pointer.
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
static tmp< GeometricField< typename innerProduct< vector, Type >::type, fvsPatchField, surfaceMesh > > dotInterpolate(const surfaceVectorField &Sf, const GeometricField< Type, fvPatchField, volMesh > &tvf)
Interpolate field onto faces.
tmp< GeometricField< cmptType, PatchField, GeoMesh > > component(const direction) const
Return a component of the field.
Abstract base class with a fat-interface to all derived classes covering all possible ways in which t...
::Foam::direction nComponents(const expressions::valueTypeCode) noexcept
The number of components associated with given valueTypeCode.
Ignore writing from objectRegistry::writeObject()
Basic second-order laplacian using face-gradients and Gauss' theorem.
#define forAll(list, i)
Loop across all elements in list.
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for expressions::valueTypeCode::INVALID.
Calculate the gradient of the given field.
A class for handling words, derived from Foam::string.
virtual tmp< Field< Type > > gradientInternalCoeffs() const
Return the matrix diagonal coefficients corresponding to the.
const objectRegistry & db() const noexcept
Return the local objectRegistry.
static tmp< T > New(Args &&... args)
Construct tmp with forwarding arguments.
A special matrix type and solver, designed for finite volume solutions of scalar equations. Face addressing is used to make all matrix assembly and solution loops vectorise.
tmp< GeometricField< Type, fvPatchField, volMesh > > fvcLaplacian(const GeometricField< Type, fvPatchField, volMesh > &)
Calculate the divergence of the given field.
const Mesh & mesh() const noexcept
Return mesh.
const fileName & instance() const noexcept
Read access to instance path component.
faceFluxFieldPtrType & faceFluxCorrectionPtr()
Return pointer to face-flux non-orthogonal correction field.
Field< Type > & source() noexcept
const scalarField & upper() const
virtual tmp< Field< Type > > gradientBoundaryCoeffs() const
Return the matrix source coefficients corresponding to the.
Mesh data needed to do the Finite Volume discretisation.
A special matrix type and solver, designed for finite volume solutions of scalar equations.
const FieldField< Field, Type > & internalCoeffs() const noexcept
fvBoundary scalar field containing pseudo-matrix coeffs for internal cells
const FieldField< Field, Type > & boundaryCoeffs() const noexcept
fvBoundary scalar field containing pseudo-matrix coeffs for boundary cells
A class for managing temporary objects.
static tmp< fvMatrix< Type > > fvmLaplacianUncorrected(const surfaceScalarField &gammaMagSf, const surfaceScalarField &deltaCoeffs, const GeometricField< Type, fvPatchField, volMesh > &)
Registry of regIOobjects.
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))
Defines the attributes of an object for which implicit objectRegistry management is supported...
tmp< fvMatrix< Type > > fvmLaplacian(const GeometricField< GType, fvsPatchField, surfaceMesh > &, const GeometricField< Type, fvPatchField, volMesh > &)
An abstract base class with a fat-interface to all derived classes covering all possible ways in whic...
const Boundary & boundaryField() const noexcept
Return const-reference to the boundary field.
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > snGrad(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
const dimensionSet & dimensions() const noexcept
Return dimensions.