setRDeltaT.H
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8 License
9  This file is part of OpenFOAM.
10 
11  OpenFOAM is free software: you can redistribute it and/or modify it
12  under the terms of the GNU General Public License as published by
13  the Free Software Foundation, either version 3 of the License, or
14  (at your option) any later version.
15 
16  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
19  for more details.
20 
21  You should have received a copy of the GNU General Public License
22  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
23 
24 \*---------------------------------------------------------------------------*/
25 
26 {
27  volScalarField& rDeltaT = trDeltaT();
28 
29  const dictionary& pimpleDict = pimple.dict();
30 
31  // Maximum flow Courant number
32  scalar maxCo(readScalar(pimpleDict.lookup("maxCo")));
33 
34  // Maximum time scale
35  scalar maxDeltaT(pimpleDict.lookupOrDefault<scalar>("maxDeltaT", GREAT));
36 
37  // Smoothing parameter (0-1) when smoothing iterations > 0
39  (
40  pimpleDict.lookupOrDefault<scalar>("rDeltaTSmoothingCoeff", 0.1)
41  );
42 
43  // Damping coefficient (1-0)
44  scalar rDeltaTDampingCoeff
45  (
46  pimpleDict.lookupOrDefault<scalar>("rDeltaTDampingCoeff", 0.2)
47  );
48 
49  // Maximum change in cell temperature per iteration
50  // (relative to previous value)
51  scalar alphaTemp(pimpleDict.lookupOrDefault("alphaTemp", 0.05));
52 
53 
54  Info<< "Time scales min/max:" << endl;
55 
56  // Cache old reciprocal time scale field
57  volScalarField rDeltaT0("rDeltaT0", rDeltaT);
58 
59  // Flow time scale
60  {
61  rDeltaT.dimensionedInternalField() =
62  (
64  /((2*maxCo)*mesh.V()*rho.dimensionedInternalField())
65  );
66 
67  // Limit the largest time scale
68  rDeltaT.max(1/maxDeltaT);
69 
70  Info<< " Flow = "
71  << gMin(1/rDeltaT.internalField()) << ", "
72  << gMax(1/rDeltaT.internalField()) << endl;
73  }
74 
75  // Reaction source time scale
76  {
77  volScalarField::DimensionedInternalField rDeltaTT
78  (
79  mag
80  (
81  parcels.hsTrans()/(mesh.V()*runTime.deltaT())
82  + combustion->Sh()()
83  )
84  /(
85  alphaTemp
86  *rho.dimensionedInternalField()
88  *T.dimensionedInternalField()
89  )
90  );
91 
92  Info<< " Temperature = "
93  << gMin(1/(rDeltaTT.field() + VSMALL)) << ", "
94  << gMax(1/(rDeltaTT.field() + VSMALL)) << endl;
95 
96  rDeltaT.dimensionedInternalField() = max
97  (
98  rDeltaT.dimensionedInternalField(),
99  rDeltaTT
100  );
101  }
102 
103  // Update tho boundary values of the reciprocal time-step
104  rDeltaT.correctBoundaryConditions();
105 
106  // Spatially smooth the time scale field
107  if (rDeltaTSmoothingCoeff < 1.0)
108  {
110  }
111 
112  // Limit rate of change of time scale
113  // - reduce as much as required
114  // - only increase at a fraction of old time scale
115  if
116  (
117  rDeltaTDampingCoeff < 1.0
118  && runTime.timeIndex() > runTime.startTimeIndex() + 1
119  )
120  {
121  rDeltaT = max
122  (
123  rDeltaT,
124  (scalar(1.0) - rDeltaTDampingCoeff)*rDeltaT0
125  );
126  }
127 
128  Info<< " Overall = "
129  << gMin(1/rDeltaT.internalField())
130  << ", " << gMax(1/rDeltaT.internalField()) << endl;
131 }
132 
133 
134 // ************************************************************************* //
Info<< "Creating combustion model\n"<< endl;autoPtr< combustionModels::psiCombustionModel > combustion(combustionModels::psiCombustionModel::New( mesh ))
#define readScalar
Definition: doubleScalar.C:38
dimensioned< scalar > mag(const dimensioned< Type > &)
scalar maxCo(readScalar(pimpleDict.lookup("maxCo")))
tmp< GeometricField< Type, fvPatchField, volMesh > > surfaceSum(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
messageStream Info
dynamicFvMesh & mesh
const dictionary & pimple
const volScalarField & T
Definition: createFields.H:25
scalar rDeltaTSmoothingCoeff(pimpleDict.lookupOrDefault< scalar >("rDeltaTSmoothingCoeff", 0.1))
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
scalar alphaTemp(pimpleDict.lookupOrDefault("alphaTemp", 0.05))
volScalarField rDeltaT0("rDeltaT0", rDeltaT)
Type gMin(const FieldField< Field, Type > &f)
surfaceScalarField & phi
rDeltaT dimensionedInternalField()
tmp< volScalarField > trDeltaT
Definition: createRDeltaT.H:3
scalar rDeltaTDampingCoeff(pimpleDict.lookupOrDefault< scalar >("rDeltaTDampingCoeff", 1))
psiReactionThermo & thermo
Definition: createFields.H:32
scalar maxDeltaT(pimpleDict.lookupOrDefault< scalar >("maxDeltaT", GREAT))
Type gMax(const FieldField< Field, Type > &f)
void smooth(volScalarField &field, const scalar coeff)
Definition: fvcSmooth.C:35