44 scalar boundedBackwardFaDdtScheme::deltaT_()
const 50 scalar boundedBackwardFaDdtScheme::deltaT0_()
const 76 scalar deltaT = deltaT_();
77 scalar deltaT0 = deltaT0_();
79 scalar coefft = 1 + deltaT/(deltaT + deltaT0);
80 scalar coefft00 = deltaT*deltaT/(deltaT0*(deltaT + deltaT0));
81 scalar coefft0 = coefft + coefft00;
124 "ddt("+dt.
name()+
')',
131 scalar deltaT = deltaT_();
132 scalar deltaT0 = deltaT0_();
134 scalar coefft = 1 + deltaT/(deltaT + deltaT0);
135 scalar coefft00 = deltaT*deltaT/(deltaT0*(deltaT + deltaT0));
136 scalar coefft0 = coefft + coefft00;
144 -rDeltaT*(coefft0 - coefft00)*dt
150 tdtdt0.
ref().primitiveFieldRef() = (-rDeltaT.
value()*dt.value())*
169 "ddt("+vf.
name()+
')',
176 scalar deltaT = deltaT_();
177 scalar deltaT0 = deltaT0_(vf);
240 return tmp<areaScalarField>
266 "ddt0("+vf.
name()+
')',
273 scalar deltaT = deltaT_();
274 scalar deltaT0 = deltaT0_(vf);
335 return tmp<areaScalarField>
361 "ddt("+
rho.name()+
','+vf.
name()+
')',
368 scalar deltaT = deltaT_();
369 scalar deltaT0 = deltaT0_(vf);
432 return tmp<areaScalarField>
458 "ddt0("+
rho.name()+
','+vf.
name()+
')',
465 scalar deltaT = deltaT_();
466 scalar deltaT0 = deltaT0_(vf);
527 return tmp<areaScalarField>
553 "ddt("+
rho.name()+
','+vf.
name()+
')',
560 scalar deltaT = deltaT_();
561 scalar deltaT0 = deltaT0_(vf);
603 rho.oldTime().primitiveField()*
606 rho.oldTime().oldTime().primitiveField()
615 rho.oldTime().boundaryField()*
618 rho.oldTime().oldTime().boundaryField()*
627 return tmp<areaScalarField>
654 "ddt0("+
rho.name()+
','+vf.
name()+
')',
661 scalar deltaT = deltaT_();
662 scalar deltaT0 = deltaT0_(vf);
703 rho.oldTime().primitiveField()*
706 rho.oldTime().oldTime().primitiveField()*
714 rho.oldTime().boundaryField()*
717 rho.oldTime().oldTime().boundaryField()*
726 return tmp<areaScalarField>
758 scalar rDeltaT = 1.0/deltaT_();
760 scalar deltaT = deltaT_();
761 scalar deltaT0 = deltaT0_(vf);
789 fam.diag() = (coefft*rDeltaT)*
mesh().
S();
793 fam.source() = rDeltaT*
829 scalar rDeltaT = 1.0/deltaT_();
831 scalar deltaT = deltaT_();
832 scalar deltaT0 = deltaT0_(vf);
860 fam.diag() = (coefft*rDeltaT*
rho.value())*
mesh().
S();
864 fam.source() = rDeltaT*
rho.value()*
900 scalar rDeltaT = 1.0/deltaT_();
902 scalar deltaT = deltaT_();
903 scalar deltaT0 = deltaT0_(vf);
912 rho.oldTime().oldTime().internalField()*
914 -
rho.oldTime().oldTime().oldTime().internalField()*
920 rho.oldTime().internalField()*
922 -
rho.oldTime().oldTime().internalField()*
935 fam.diag() = (coefft*rDeltaT)*
rho.primitiveField()*
mesh().
S();
939 fam.source() = rDeltaT*
941 coefft0*
rho.oldTime().primitiveField()
943 - coefft00*
rho.oldTime().oldTime().primitiveField()
951 coefft0*
rho.oldTime().primitiveField()
953 - coefft00*
rho.oldTime().oldTime().primitiveField()
const Type & value() const noexcept
Return const reference to value.
const GeometricField< Type, PatchField, GeoMesh > & oldTime() const
Return old time field.
const Internal::FieldType & primitiveField() const noexcept
Return a const-reference to the internal field values.
dimensionedScalar deltaT0() const
Return old time step.
tmp< areaScalarField > facDdt0(const dimensionedScalar)
const DimensionedField< scalar, areaMesh > & S00() const
Return old-old-time face areas.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
tmp< areaScalarField > facDdt(const dimensionedScalar)
const Time & time() const
Return reference to time.
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.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
const Internal & internalField() const noexcept
Return a const-reference to the dimensioned internal field.
faDdtScheme< scalar >::addIstreamConstructorToTable< boundedBackwardFaDdtScheme > addboundedBackwardFaDdtSchemeIstreamConstructorToTable_
const DimensionedField< scalar, areaMesh > & S() const
Return face areas.
Ignore writing from objectRegistry::writeObject()
virtual const objectRegistry & thisDb() const
Return reference to the mesh database.
dimensionedScalar pos(const dimensionedScalar &ds)
tmp< areaScalarField > limiter(const areaScalarField &phi)
tmp< faScalarMatrix > famDdt(const areaScalarField &)
defineTypeNameAndDebug(limitHeight, 0)
static tmp< T > New(Args &&... args)
Construct tmp with forwarding arguments.
const dimensionSet & dimensions() const noexcept
Return const reference to dimensions.
const faMesh & mesh() const
Return mesh reference.
static word timeName(const scalar t, const int precision=precision_)
Return time name of given scalar time formatted with the given precision.
Calculate the divergence of the given field.
static const word & calculatedType() noexcept
The type name for calculated patch fields.
const word & name() const noexcept
Return const reference to name.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
bool moving() const
Is mesh moving.
Template specialisation for scalar faMatrix.
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
A class for managing temporary objects.
Defines the attributes of an object for which implicit objectRegistry management is supported...
Calculate the matrix for the second temporal derivative.
const Boundary & boundaryField() const noexcept
Return const-reference to the boundary field.
dimensionedScalar deltaT() const
Return time step.
const dimensionSet dimArea(sqr(dimLength))
const DimensionedField< scalar, areaMesh > & S0() const
Return old-time face areas.
const dimensionSet & dimensions() const noexcept
Return dimensions.
GeometricField< scalar, faPatchField, areaMesh > areaScalarField
static constexpr const zero Zero
Global zero (0)