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.
tmp< areaScalarField > facDdt0(const dimensionedScalar)
scalar deltaTValue() const noexcept
Return time step value.
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()
tmp< areaScalarField > limiter(const areaScalarField &phi)
virtual const objectRegistry & thisDb() const
Return reference to the mesh database.
dimensionedScalar pos(const dimensionedScalar &ds)
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 a time name for the given scalar time value formatted with the given precision.
Calculate the divergence of the given field.
scalar deltaT0Value() const noexcept
Return old time step value.
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)