instabilityG.C
Go to the documentation of this file.
1/*---------------------------------------------------------------------------*\
2 ========= |
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4 \\ / O peration |
5 \\ / A nd | www.openfoam.com
6 \\/ M anipulation |
7-------------------------------------------------------------------------------
8 Copyright (C) 2011-2015 OpenFOAM Foundation
9-------------------------------------------------------------------------------
10License
11 This file is part of OpenFOAM.
12
13 OpenFOAM is free software: you can redistribute it and/or modify it
14 under the terms of the GNU General Public License as published by
15 the Free Software Foundation, either version 3 of the License, or
16 (at your option) any later version.
17
18 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
19 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
20 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
21 for more details.
22
23 You should have received a copy of the GNU General Public License
24 along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
25
26\*---------------------------------------------------------------------------*/
27
28#include "instabilityG.H"
30
31// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
32
33namespace Foam
34{
35namespace XiGModels
36{
37 defineTypeNameAndDebug(instabilityG, 0);
38 addToRunTimeSelectionTable(XiGModel, instabilityG, dictionary);
39}
40}
41
42
43// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
44
46(
47 const dictionary& XiGProperties,
48 const psiuReactionThermo& thermo,
49 const compressible::RASModel& turbulence,
50 const volScalarField& Su
51)
52:
53 XiGModel(XiGProperties, thermo, turbulence, Su),
54 GIn_("GIn", dimless/dimTime, XiGModelCoeffs_),
55 lambdaIn_("lambdaIn", dimLength, XiGModelCoeffs_),
56 XiGModel_(XiGModel::New(XiGModelCoeffs_, thermo, turbulence, Su))
57{}
58
59
60// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
61
63{}
64
65
66// * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
67
69{
70 volScalarField turbXiG(XiGModel_->G());
71 return (GIn_*GIn_/(GIn_ + turbXiG) + turbXiG);
72}
73
74
76{
77 const objectRegistry& db = Su_.db();
78 const volScalarField& Xi = db.lookupObject<volScalarField>("Xi");
79 const volScalarField& rho = db.lookupObject<volScalarField>("rho");
80 const volScalarField& mgb = db.lookupObject<volScalarField>("mgb");
81
82 return XiGModel_->Db()
83 + rho*Su_*(Xi - 1.0)*mgb*(0.5*lambdaIn_)/(mgb + 1.0/lambdaIn_);
84}
85
86
87bool Foam::XiGModels::instabilityG::read(const dictionary& XiGProperties)
88{
89 XiGModel::read(XiGProperties);
90
91 XiGModelCoeffs_.readEntry("GIn", GIn_);
92 XiGModelCoeffs_.readEntry("lambdaIn", lambdaIn_);
93
94 return true;
95}
96
97
98// ************************************************************************* //
Macros for easy insertion into run-time selection tables.
#define addToRunTimeSelectionTable(baseType, thisType, argNames)
Add to construction table with typeName as the key.
virtual bool read()
Re-read model coefficients if they have changed.
Flame-surface instabilityG flame-wrinkling generation rate coefficient model used in XiModel....
Definition: instabilityG.H:58
virtual tmp< volScalarField > G() const
Return the flame-wrinkling generation rate.
virtual tmp< volScalarField > Db() const
Return the flame diffusivity.
virtual ~instabilityG()
Destructor.
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
A class for managing temporary objects.
Definition: tmp.H:65
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
Definition: className.H:121
compressible::turbulenceModel & turbulence
zeroField Su
Definition: alphaSuSp.H:1
Namespace for OpenFOAM.
const dimensionSet dimless
Dimensionless.
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
Definition: dimensionSets.H:52
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
Definition: dimensionSets.H:53
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:82
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh > > &tdf1, const word &name, const dimensionSet &dimensions)
Global function forwards to reuseTmpDimensionedField::New.