setRDeltaT.H
Go to the documentation of this file.
1 {
2  volScalarField& rDeltaT = trDeltaT.ref();
3 
4  const dictionary& pimpleDict = pimple.dict();
5 
6  scalar maxCo
7  (
8  pimpleDict.lookupOrDefault<scalar>("maxCo", 0.8)
9  );
10 
12  (
13  pimpleDict.lookupOrDefault<scalar>("rDeltaTSmoothingCoeff", 0.02)
14  );
15 
16  scalar rDeltaTDampingCoeff
17  (
18  pimpleDict.lookupOrDefault<scalar>("rDeltaTDampingCoeff", 1.0)
19  );
20 
21  scalar maxDeltaT
22  (
23  pimpleDict.lookupOrDefault<scalar>("maxDeltaT", GREAT)
24  );
25 
26  volScalarField rDeltaT0("rDeltaT0", rDeltaT);
27 
28  // Set the reciprocal time-step from the local Courant number
29  rDeltaT.ref() = max
30  (
31  1/dimensionedScalar("maxDeltaT", dimTime, maxDeltaT),
32  fvc::surfaceSum(mag(phi))()()
33  /((2*maxCo)*mesh.V()*rho())
34  );
35 
36  if (pimple.transonic())
37  {
39  (
40  "phid",
42  );
43 
44  rDeltaT.ref() = max
45  (
46  rDeltaT(),
47  fvc::surfaceSum(mag(phid))()()
48  /((2*maxCo)*mesh.V()*psi())
49  );
50  }
51 
52  // Update tho boundary values of the reciprocal time-step
53  rDeltaT.correctBoundaryConditions();
54 
55  Info<< "Flow time scale min/max = "
56  << gMin(1/rDeltaT.primitiveField())
57  << ", " << gMax(1/rDeltaT.primitiveField()) << endl;
58 
59  if (rDeltaTSmoothingCoeff < 1.0)
60  {
62  }
63 
64  Info<< "Smoothed flow time scale min/max = "
65  << gMin(1/rDeltaT.primitiveField())
66  << ", " << gMax(1/rDeltaT.primitiveField()) << endl;
67 
68  // Limit rate of change of time scale
69  // - reduce as much as required
70  // - only increase at a fraction of old time scale
71  if
72  (
74  && runTime.timeIndex() > runTime.startTimeIndex() + 1
75  )
76  {
77  rDeltaT =
78  rDeltaT0
79  *max(rDeltaT/rDeltaT0, scalar(1) - rDeltaTDampingCoeff);
80 
81  Info<< "Damped flow time scale min/max = "
82  << gMin(1/rDeltaT.primitiveField())
83  << ", " << gMax(1/rDeltaT.primitiveField()) << endl;
84  }
85 }
scalar maxCo(readScalar(pimpleDict.lookup("maxCo")))
surfaceScalarField & phi
scalar maxDeltaT(pimpleDict.lookupOrDefault< scalar >("maxDeltaT", GREAT))
U
Definition: pEqn.H:83
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
Type gMin(const FieldField< Field, Type > &f)
const dictionary & pimple
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:253
void smooth(volScalarField &field, const scalar coeff)
Definition: fvcSmooth.C:35
volScalarField rDeltaT0("rDeltaT0", rDeltaT)
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
dynamicFvMesh & mesh
scalar rDeltaTDampingCoeff(pimpleDict.lookupOrDefault< scalar >("rDeltaTDampingCoeff", 1.0))
Type gMax(const FieldField< Field, Type > &f)
static tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > interpolate(const GeometricField< Type, fvPatchField, volMesh > &tvf, const surfaceScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
surfaceScalarField phid("phid", fvc::interpolate(psi) *(fvc::flux(HbyA)+rhorAUf *fvc::ddtCorr(rho, U, phi)/fvc::interpolate(rho)))
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
Definition: dimensionSets.H:51
messageStream Info
dimensioned< scalar > mag(const dimensioned< Type > &)
tmp< GeometricField< Type, fvPatchField, volMesh > > surfaceSum(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
const volScalarField & psi
tmp< surfaceScalarField > flux(const volVectorField &vvf)
Return the face-flux field obtained from the given volVectorField.
Definition: fvcFlux.C:32
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
tmp< volScalarField > trDeltaT
Definition: createRDeltaT.H:3
scalar rDeltaTSmoothingCoeff(pimpleDict.lookupOrDefault< scalar >("rDeltaTSmoothingCoeff", 0.1))