Second-oder Crank-Nicolson implicit ddt using the current and previous time-step fields as well as the previous time-step ddt. More...
Public Types | |
typedef ddtScheme< Type >::fluxFieldType | fluxFieldType |
![]() | |
typedef GeometricField< typename flux< Type >::type, fvsPatchField, surfaceMesh > | fluxFieldType |
Additional Inherited Members | |
![]() | |
static tmp< ddtScheme< Type > > | New (const fvMesh &mesh, Istream &schemeData) |
Return a pointer to a new ddtScheme created on freestore. More... | |
![]() | |
ddtScheme (const ddtScheme &) | |
Disallow copy construct. More... | |
void | operator= (const ddtScheme &) |
Disallow default bitwise assignment. More... | |
![]() | |
refCount () | |
Construct null initializing count to 0. More... | |
![]() | |
const fvMesh & | mesh_ |
Second-oder Crank-Nicolson implicit ddt using the current and previous time-step fields as well as the previous time-step ddt.
The Crank-Nicolson scheme is often unstable for complex flows in complex geometries and it is necessary to "off-centre" the scheme to stabilize it while retaining greater temporal accuracy than the first-order Euler-implicit scheme. Off-centering is specified via the mandatory coefficient ocCoeff
in the range [0,1] following the scheme name e.g.
ddtSchemes { default CrankNicolson 0.9; }
or with an optional "ramp" function to transition from the Euler scheme to Crank-Nicolson over a initial period to avoid start-up problems, e.g.
ddtSchemes { default CrankNicolson ocCoeff { type scale; scale linearRamp; duration 0.01; value 0.9; }; }
With a coefficient of 1 the scheme is fully centred and second-order, with a coefficient of 0 the scheme is equivalent to Euler-implicit. A coefficient of 0.9 has been found to be suitable for a range of cases for which higher-order accuracy in time is needed and provides similar accuracy and stability to the backward scheme. However, the backward scheme has been found to be more robust and provides formal second-order accuracy in time.
The advantage of the Crank-Nicolson scheme over backward is that only the new and old-time values are needed, the additional terms relating to the fluxes and sources are evaluated at the mid-point of the time-step which provides the opportunity to limit the fluxes in such a way as to ensure boundedness while maintaining greater accuracy in time compared to the Euler-implicit scheme. This approach is now used with MULES in the interFoam family of solvers. Boundedness cannot be guaranteed with the backward scheme.
cnCoeff = 1.0/(1.0 + ocCoeff);
Definition at line 110 of file CrankNicolsonDdtScheme.H.
typedef ddtScheme<Type>::fluxFieldType fluxFieldType |
Definition at line 286 of file CrankNicolsonDdtScheme.H.
CrankNicolsonDdtScheme | ( | const fvMesh & | mesh | ) |
Construct from mesh.
Definition at line 281 of file CrankNicolsonDdtScheme.C.
References polyMesh::moving(), and fvMesh::V00().
CrankNicolsonDdtScheme | ( | const fvMesh & | mesh, |
Istream & | is | ||
) |
Construct from mesh and Istream.
Definition at line 297 of file CrankNicolsonDdtScheme.C.
References dict, Foam::exit(), Foam::FatalIOError, FatalIOErrorInFunction, CrankNicolsonDdtScheme< Type >::fvcDdt(), token::isNumber(), Function1< Type >::New(), token::number(), CrankNicolsonDdtScheme< Type >::ocCoeff(), and Istream::putBack().
TypeName | ( | "CrankNicolson" | ) |
Runtime type information.
|
inline |
Return mesh reference.
Definition at line 222 of file CrankNicolsonDdtScheme.H.
References ddtScheme< Type >::mesh().
Referenced by CrankNicolsonDdtScheme< Type >::fvcDdt(), CrankNicolsonDdtScheme< Type >::fvcDdtPhiCorr(), CrankNicolsonDdtScheme< Type >::fvcDdtUfCorr(), CrankNicolsonDdtScheme< Type >::fvmDdt(), CrankNicolsonDdtScheme< Type >::meshPhi(), and CrankNicolsonDdtScheme< Type >::ocCoeff().
|
inline |
Return the current off-centreing coefficient.
Definition at line 228 of file CrankNicolsonDdtScheme.H.
References Foam::constant::atomic::alpha, CrankNicolsonDdtScheme< Type >::fvcDdt(), CrankNicolsonDdtScheme< Type >::fvmDdt(), CrankNicolsonDdtScheme< Type >::mesh(), psi, and rho.
Referenced by CrankNicolsonDdtScheme< Type >::CrankNicolsonDdtScheme().
|
virtual |
Implements ddtScheme< Type >.
Definition at line 346 of file CrankNicolsonDdtScheme.C.
References dimensioned< Type >::dimensions(), Foam::dimTime, CrankNicolsonDdtScheme< Type >::mesh(), dimensioned< Type >::name(), fvMesh::time(), Time::timeName(), fvMesh::V(), fvMesh::V0(), fvMesh::V00(), and Foam::Zero.
Referenced by CrankNicolsonDdtScheme< Type >::CrankNicolsonDdtScheme(), CrankNicolsonDdtScheme< Type >::fvcDdt(), and CrankNicolsonDdtScheme< Type >::ocCoeff().
|
virtual |
Implements ddtScheme< Type >.
Definition at line 408 of file CrankNicolsonDdtScheme.C.
References GeometricField< Type, PatchField, GeoMesh >::boundaryField(), DimensionedField< Type, GeoMesh >::dimensions(), Foam::fv::ff(), CrankNicolsonDdtScheme< Type >::fvcDdt(), CrankNicolsonDdtScheme< Type >::mesh(), IOobject::name(), GeometricField< Type, PatchField, GeoMesh >::oldTime(), fvMesh::time(), Time::timeName(), fvMesh::V(), fvMesh::V0(), fvMesh::V00(), and dimensioned< Type >::value().
|
virtual |
Implements ddtScheme< Type >.
Definition at line 495 of file CrankNicolsonDdtScheme.C.
References GeometricField< Type, PatchField, GeoMesh >::boundaryField(), dimensioned< Type >::dimensions(), DimensionedField< Type, GeoMesh >::dimensions(), Foam::fv::ff(), CrankNicolsonDdtScheme< Type >::fvcDdt(), CrankNicolsonDdtScheme< Type >::mesh(), dimensioned< Type >::name(), IOobject::name(), GeometricField< Type, PatchField, GeoMesh >::oldTime(), GeometricField< Type, PatchField, GeoMesh >::primitiveField(), fvMesh::time(), Time::timeName(), fvMesh::V(), fvMesh::V0(), fvMesh::V00(), and dimensioned< Type >::value().
|
virtual |
Implements ddtScheme< Type >.
Definition at line 585 of file CrankNicolsonDdtScheme.C.
References GeometricField< Type, PatchField, GeoMesh >::boundaryField(), DimensionedField< Type, GeoMesh >::dimensions(), Foam::fv::ff(), CrankNicolsonDdtScheme< Type >::fvcDdt(), CrankNicolsonDdtScheme< Type >::mesh(), IOobject::name(), GeometricField< Type, PatchField, GeoMesh >::oldTime(), GeometricField< Type, PatchField, GeoMesh >::primitiveField(), fvMesh::time(), Time::timeName(), fvMesh::V(), fvMesh::V0(), fvMesh::V00(), and dimensioned< Type >::value().
|
virtual |
Implements ddtScheme< Type >.
Definition at line 685 of file CrankNicolsonDdtScheme.C.
References GeometricField< Type, PatchField, GeoMesh >::boundaryField(), DimensionedField< Type, GeoMesh >::dimensions(), Foam::fv::ff(), CrankNicolsonDdtScheme< Type >::fvmDdt(), CrankNicolsonDdtScheme< Type >::mesh(), IOobject::name(), GeometricField< Type, PatchField, GeoMesh >::oldTime(), GeometricField< Type, PatchField, GeoMesh >::primitiveField(), fvMesh::time(), Time::timeName(), fvMesh::V(), fvMesh::V0(), fvMesh::V00(), and dimensioned< Type >::value().
|
virtual |
Implements ddtScheme< Type >.
Definition at line 815 of file CrankNicolsonDdtScheme.C.
References lduMatrix::diag(), DimensionedField< Type, GeoMesh >::dimensions(), Foam::dimTime, Foam::dimVol, Foam::fv::ff(), CrankNicolsonDdtScheme< Type >::mesh(), IOobject::name(), GeometricField< Type, PatchField, GeoMesh >::oldTime(), fvMatrix< Type >::source(), fvMesh::V(), fvMesh::V0(), fvMesh::V00(), and dimensioned< Type >::value().
Referenced by CrankNicolsonDdtScheme< Type >::fvcDdt(), CrankNicolsonDdtScheme< Type >::fvmDdt(), and CrankNicolsonDdtScheme< Type >::ocCoeff().
|
virtual |
Implements ddtScheme< Type >.
Definition at line 897 of file CrankNicolsonDdtScheme.C.
References lduMatrix::diag(), dimensioned< Type >::dimensions(), DimensionedField< Type, GeoMesh >::dimensions(), Foam::dimTime, Foam::dimVol, Foam::fv::ff(), CrankNicolsonDdtScheme< Type >::fvmDdt(), CrankNicolsonDdtScheme< Type >::mesh(), dimensioned< Type >::name(), IOobject::name(), GeometricField< Type, PatchField, GeoMesh >::oldTime(), fvMatrix< Type >::source(), fvMesh::V(), fvMesh::V0(), fvMesh::V00(), and dimensioned< Type >::value().
|
virtual |
Implements ddtScheme< Type >.
Definition at line 979 of file CrankNicolsonDdtScheme.C.
References GeometricField< Type, PatchField, GeoMesh >::boundaryField(), lduMatrix::diag(), DimensionedField< Type, GeoMesh >::dimensions(), Foam::dimTime, Foam::dimVol, Foam::fv::ff(), CrankNicolsonDdtScheme< Type >::fvmDdt(), CrankNicolsonDdtScheme< Type >::mesh(), IOobject::name(), GeometricField< Type, PatchField, GeoMesh >::oldTime(), GeometricField< Type, PatchField, GeoMesh >::primitiveField(), fvMatrix< Type >::source(), fvMesh::V(), fvMesh::V0(), fvMesh::V00(), and dimensioned< Type >::value().
|
virtual |
Implements ddtScheme< Type >.
Definition at line 1069 of file CrankNicolsonDdtScheme.C.
References GeometricField< Type, PatchField, GeoMesh >::boundaryField(), lduMatrix::diag(), DimensionedField< Type, GeoMesh >::dimensions(), Foam::dimTime, Foam::dimVol, Foam::fv::ff(), CrankNicolsonDdtScheme< Type >::fvcDdtUfCorr(), CrankNicolsonDdtScheme< Type >::mesh(), IOobject::name(), GeometricField< Type, PatchField, GeoMesh >::oldTime(), GeometricField< Type, PatchField, GeoMesh >::primitiveField(), fvMatrix< Type >::source(), fvMesh::V(), fvMesh::V0(), fvMesh::V00(), and dimensioned< Type >::value().
|
virtual |
Implements ddtScheme< Type >.
Definition at line 1180 of file CrankNicolsonDdtScheme.C.
References DimensionedField< Type, GeoMesh >::dimensions(), ddtScheme< Type >::fvcDdtPhiCoeff(), CrankNicolsonDdtScheme< Type >::fvcDdtPhiCorr(), Foam::fvc::interpolate(), CrankNicolsonDdtScheme< Type >::mesh(), IOobject::name(), GeometricField< Type, PatchField, GeoMesh >::oldTime(), fvMesh::Sf(), fvMesh::time(), and Time::timeName().
Referenced by CrankNicolsonDdtScheme< Type >::fvcDdtPhiCorr(), CrankNicolsonDdtScheme< Type >::fvcDdtUfCorr(), and CrankNicolsonDdtScheme< Type >::fvmDdt().
|
virtual |
Implements ddtScheme< Type >.
Definition at line 1241 of file CrankNicolsonDdtScheme.C.
References DimensionedField< Type, GeoMesh >::dimensions(), Foam::fvc::dotInterpolate(), ddtScheme< Type >::fvcDdtPhiCoeff(), CrankNicolsonDdtScheme< Type >::fvcDdtUfCorr(), CrankNicolsonDdtScheme< Type >::mesh(), IOobject::name(), GeometricField< Type, PatchField, GeoMesh >::oldTime(), fvMesh::Sf(), fvMesh::time(), and Time::timeName().
Referenced by CrankNicolsonDdtScheme< Type >::fvcDdtPhiCorr(), and CrankNicolsonDdtScheme< Type >::fvcDdtUfCorr().
|
virtual |
Implements ddtScheme< Type >.
Definition at line 1303 of file CrankNicolsonDdtScheme.C.
References Foam::abort(), Foam::fvc::ddtCorr(), DimensionedField< Type, GeoMesh >::dimensions(), Foam::dimVelocity, Foam::FatalError, FatalErrorInFunction, ddtScheme< Type >::fvcDdtPhiCoeff(), CrankNicolsonDdtScheme< Type >::fvcDdtPhiCorr(), CrankNicolsonDdtScheme< Type >::fvcDdtUfCorr(), Foam::fvc::interpolate(), CrankNicolsonDdtScheme< Type >::mesh(), IOobject::name(), GeometricField< Type, PatchField, GeoMesh >::null(), GeometricField< Type, PatchField, GeoMesh >::oldTime(), fvMesh::Sf(), fvMesh::time(), and Time::timeName().
|
virtual |
Implements ddtScheme< Type >.
Definition at line 1398 of file CrankNicolsonDdtScheme.C.
References Foam::abort(), Foam::fvc::ddtCorr(), Foam::dimArea, DimensionedField< Type, GeoMesh >::dimensions(), Foam::dimVelocity, Foam::fvc::dotInterpolate(), Foam::FatalError, FatalErrorInFunction, ddtScheme< Type >::fvcDdtPhiCoeff(), CrankNicolsonDdtScheme< Type >::fvcDdtPhiCorr(), CrankNicolsonDdtScheme< Type >::mesh(), CrankNicolsonDdtScheme< Type >::meshPhi(), IOobject::name(), GeometricField< Type, PatchField, GeoMesh >::null(), GeometricField< Type, PatchField, GeoMesh >::oldTime(), fvMesh::Sf(), fvMesh::time(), and Time::timeName().
|
virtual |
Implements ddtScheme< Type >.
Definition at line 1493 of file CrankNicolsonDdtScheme.C.
References Foam::dimVolume, CrankNicolsonDdtScheme< Type >::mesh(), Foam::name(), IOobject::NO_READ, IOobject::NO_WRITE, GeometricField< Type, PatchField, GeoMesh >::oldTime(), phi, fvMesh::phi(), and timeName.
Referenced by CrankNicolsonDdtScheme< Type >::fvcDdtPhiCorr().
tmp< surfaceScalarField > fvcDdtUfCorr | ( | const GeometricField< scalar, fvPatchField, volMesh > & | U, |
const GeometricField< scalar, fvsPatchField, surfaceMesh > & | Uf | ||
) |
tmp< surfaceScalarField > fvcDdtPhiCorr | ( | const volScalarField & | U, |
const surfaceScalarField & | phi | ||
) |
tmp< surfaceScalarField > fvcDdtUfCorr | ( | const volScalarField & | rho, |
const volScalarField & | U, | ||
const surfaceScalarField & | Uf | ||
) |
tmp< surfaceScalarField > fvcDdtPhiCorr | ( | const volScalarField & | rho, |
const volScalarField & | U, | ||
const surfaceScalarField & | phi | ||
) |