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.getOrDefault<scalar>("maxCo", 0.8)
9 );
10
12 (
13 pimpleDict.getOrDefault<scalar>("rDeltaTSmoothingCoeff", 0.02)
14 );
15
17 (
18 pimpleDict.getOrDefault<scalar>("rDeltaTDampingCoeff", 1.0)
19 );
20
21 scalar maxDeltaT
22 (
23 pimpleDict.getOrDefault<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 {
38 surfaceScalarField phid
39 (
40 "phid",
41 fvc::interpolate(psi)*fvc::flux(U)
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 the 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 {
61 fvc::smooth(rDeltaT, rDeltaTSmoothingCoeff);
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 =
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}
Y[inertIndex] max(0.0)
surfaceScalarField & phi
pimpleControl & pimple
U
Definition: pEqn.H:72
const volScalarField & psi
scalar rDeltaTDampingCoeff(pimpleDict.getOrDefault< scalar >("rDeltaTDampingCoeff", 1.0))
const dictionary & pimpleDict
Definition: setRDeltaT.H:32
scalar rDeltaTSmoothingCoeff(pimpleDict.getOrDefault< scalar >("rDeltaTSmoothingCoeff", 0.1))
volScalarField rDeltaT0("rDeltaT0", rDeltaT)
surfaceScalarField phid("phid", fvc::interpolate(psi) *(fvc::flux(HbyA)+MRF.zeroFilter(rhorAUf *fvc::ddtCorr(rho, U, phi)/fvc::interpolate(rho))))
dynamicFvMesh & mesh
engineTime & runTime
tmp< volScalarField > trDeltaT
Definition: createRDeltaT.H:3
scalar maxCo
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:82
messageStream Info
Information stream (stdout output on master, null elsewhere)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:372
Type gMin(const FieldField< Field, Type > &f)
Type gMax(const FieldField< Field, Type > &f)