30 #include "surfaceInterpolate.H" 58 mesh().time().timeName(),
74 tdtdt.
ref().primitiveFieldRef() =
101 "ddt("+vf.
name()+
')',
102 mesh().time().timeName(),
127 tdtdt.
ref().boundaryFieldRef().
128 template evaluateCoupled<coupledFvPatch>();
158 "ddt("+
rho.name()+
','+vf.
name()+
')',
159 mesh().time().timeName(),
184 tdtdt.
ref().boundaryFieldRef().
185 template evaluateCoupled<coupledFvPatch>();
191 return tmp<GeometricField<Type, fvPatchField, volMesh>>
193 new GeometricField<Type, fvPatchField, volMesh>
215 "ddt("+
rho.name()+
','+vf.
name()+
')',
216 mesh().time().timeName(),
236 -
rho.oldTime().boundaryField()
244 tdtdt.
ref().boundaryFieldRef().
245 template evaluateCoupled<coupledFvPatch>();
251 return tmp<GeometricField<Type, fvPatchField, volMesh>>
253 new GeometricField<Type, fvPatchField, volMesh>
277 mesh().time().timeName(),
300 alpha.boundaryField()
304 -
alpha.oldTime().boundaryField()
305 *
rho.oldTime().boundaryField()
313 tdtdt.
ref().boundaryFieldRef().
314 template evaluateCoupled<coupledFvPatch>();
320 return tmp<GeometricField<Type, fvPatchField, volMesh>>
322 new GeometricField<Type, fvPatchField, volMesh>
347 "ddt("+sf.
name()+
')',
348 mesh().time().timeName(),
381 scalar rDeltaT = 1.0/
mesh().time().deltaTValue();
416 scalar rDeltaT = 1.0/
mesh().time().deltaTValue();
453 scalar rDeltaT = 1.0/
mesh().time().deltaTValue();
455 fvm.
diag() = rDeltaT*
rho.primitiveField()*
mesh().Vsc();
460 *
rho.oldTime().primitiveField()
466 *
rho.oldTime().primitiveField()
493 scalar rDeltaT = 1.0/
mesh().time().deltaTValue();
496 rDeltaT*
alpha.primitiveField()*
rho.primitiveField()*
mesh().Vsc();
501 *
alpha.oldTime().primitiveField()
502 *
rho.oldTime().primitiveField()
508 *
alpha.oldTime().primitiveField()
509 *
rho.oldTime().primitiveField()
527 fluxFieldType phiUf0(
mesh().Sf() &
Uf.oldTime());
528 fluxFieldType phiCorr
539 "ddtCorr(" +
U.name() +
',' +
Uf.name() +
')',
540 mesh().time().timeName(),
543 this->fvcDdtPhiCoeff(
U.oldTime(), phiUf0, phiCorr)
555 const fluxFieldType&
phi 560 fluxFieldType phiCorr
571 "ddtCorr(" +
U.name() +
',' +
phi.name() +
')',
572 mesh().time().timeName(),
575 this->fvcDdtPhiCoeff(
U.oldTime(),
phi.oldTime(), phiCorr)
601 rho.oldTime()*
U.oldTime()
604 fluxFieldType phiUf0(
mesh().Sf() &
Uf.oldTime());
614 +
rho.name() +
',' +
U.name() +
',' +
Uf.name() +
')',
615 mesh().time().timeName(),
618 this->fvcDdtPhiCoeff(rhoU0, phiUf0, phiCorr,
rho.oldTime())
629 fluxFieldType phiUf0(
mesh().Sf() &
Uf.oldTime());
630 fluxFieldType phiCorr
642 +
rho.name() +
',' +
U.name() +
',' +
Uf.name() +
')',
643 mesh().time().timeName(),
659 <<
"dimensions of Uf are not correct" 662 return fluxFieldType::null();
673 const fluxFieldType&
phi 686 rho.oldTime()*
U.oldTime()
689 fluxFieldType phiCorr
701 +
rho.name() +
',' +
U.name() +
',' +
phi.name() +
')',
702 mesh().time().timeName(),
721 fluxFieldType phiCorr
726 return tmp<fluxFieldType>
733 +
rho.name() +
',' +
U.name() +
',' +
phi.name() +
')',
734 mesh().time().timeName(),
750 <<
"dimensions of phi are not correct" 753 return fluxFieldType::null();
const scalarField & diag() const
const Type & value() const noexcept
Return const reference to value.
const GeometricField< Type, PatchField, GeoMesh > & oldTime() const
Return old time field.
error FatalError
Error stream (stdout output on all processes), with additional 'FOAM FATAL ERROR' header text and sta...
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
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.
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< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tf1, const word &name, const dimensionSet &dimensions, const bool initCopy=false)
Global function forwards to reuseTmpDimensionedField::New.
const dimensionSet dimVol(dimVolume)
Older spelling for dimVolume.
tmp< GeometricField< Type, fvPatchField, volMesh > > fvcDdt(const dimensioned< Type > &)
Generic GeometricField class.
Generic dimensioned Type class.
tmp< fluxFieldType > fvcDdtPhiCorr(const GeometricField< Type, fvPatchField, volMesh > &U, const fluxFieldType &phi)
tmp< surfaceScalarField > meshPhi(const GeometricField< Type, fvPatchField, volMesh > &)
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.
autoPtr< surfaceVectorField > Uf
const dimensionSet & dimensions() const noexcept
Return const reference to dimensions.
errorManip< error > abort(error &err)
tmp< fluxFieldType > fvcDdtUfCorr(const GeometricField< Type, fvPatchField, volMesh > &U, const GeometricField< Type, fvsPatchField, surfaceMesh > &Uf)
Calculate the divergence of the given field.
const word & name() const noexcept
Return const reference to name.
Field< Type > & source() noexcept
static const word & calculatedType() noexcept
The type name for calculated patch fields.
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
A special matrix type and solver, designed for finite volume solutions of scalar equations.
A class for managing temporary objects.
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
Defines the attributes of an object for which implicit objectRegistry management is supported...
const Boundary & boundaryField() const noexcept
Return const-reference to the boundary field.
tmp< fvMatrix< Type > > fvmDdt(const GeometricField< Type, fvPatchField, volMesh > &)
const dimensionSet dimArea(sqr(dimLength))
const dimensionSet & dimensions() const noexcept
Return dimensions.
static constexpr const zero Zero
Global zero (0)
const dimensionSet dimVelocity