basicXiSubG.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 | www.openfoam.com
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8  Copyright (C) 2011-2013 OpenFOAM Foundation
9 -------------------------------------------------------------------------------
10 License
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 Class
27  Foam::XiGModels::basicSubGrid
28 
29 Description
30 
31  Basic sub-grid obstacle flame-wrinkling generation rate coefficient model.
32  Details supplied by J Puttock 2/7/06.
33 
34  \f$ G_{sub} \f$ denotes the generation coefficient and it is given by
35 
36  \f[
37  G_{sub} = k_{1} /frac{\vert \dwea{\vec{U}} \vert}{L_{obs}}
38  \frac{/Xi_{{sub}_{eq}}-1}{/Xi_{sub}}
39  \f]
40 
41  and the removal:
42 
43  \f[
44  - k_{1} /frac{\vert \dwea{\vec{U}} \vert}{L_{sub}}
45  \frac{\Xi_{sub}-1}{\Xi_{sub}}
46  \f]
47 
48  Finally, \f$ G_{sub} \f$ is added to generation rate \f$ G_{in} \f$
49  due to the turbulence.
50 
51 SourceFiles
52  basicSubGrid.C
53 
54 \*---------------------------------------------------------------------------*/
55 
56 #ifndef basicSubGrid_H
57 #define basicSubGrid_H
58 
59 #include "XiGModel.H"
60 
61 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
62 
63 namespace Foam
64 {
65 namespace XiGModels
66 {
67 
68 /*---------------------------------------------------------------------------*\
69  Class basicSubGrid Declaration
70 \*---------------------------------------------------------------------------*/
71 
72 class basicSubGrid
73 :
74  public XiGModel
75 {
76  // Private data
77 
78  //- Sub-grid generation rate coefficient
79  scalar k1;
80 
81  //- Xi generation rate model due to turbulence
82  autoPtr<XiGModel> XiGModel_;
83 
84 
85  // Private Member Functions
86 
87  //- No copy construct
88  basicSubGrid(const basicSubGrid&) = delete;
89 
90  //- No copy assignment
91  void operator=(const basicSubGrid&) = delete;
92 
93 
94 public:
95 
96  //- Runtime type information
97  TypeName("basicSubGridG");
98 
99 
100  // Constructors
101 
102  //- Construct from components
104  (
105  const dictionary& XiGProperties,
106  const psiuReactionThermo& thermo,
108  const volScalarField& Su
109  );
110 
111 
112  //- Destructor
113  virtual ~basicSubGrid();
114 
115 
116  // Member Functions
117 
118  //- Return the flame-wrinkling generation rate
119  virtual tmp<volScalarField> G() const;
120 
121  //- Return the flame diffusivity
122  virtual tmp<volScalarField> Db() const;
123 
124  //- Update properties from given dictionary
125  virtual bool read(const dictionary& XiGProperties);
126 };
127 
128 
129 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
130 
131 } // End namespace XiGModels
132 } // End namespace Foam
133 
134 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
135 
136 #endif
137 
138 // ************************************************************************* //
turbulence
Info<< "Reading field U\n"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho());volVectorField rhoU(IOobject("rhoU", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *U);volScalarField rhoE(IOobject("rhoE", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *(e+0.5 *magSqr(U)));surfaceScalarField pos(IOobject("pos", runTime.timeName(), mesh), mesh, dimensionedScalar("pos", dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar("neg", dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\n"<< endl;autoPtr< compressible::turbulenceModel > turbulence(compressible::turbulenceModel::New(rho, U, phi, thermo))
Definition: createFields.H:94
Foam::tmp
A class for managing temporary objects.
Definition: PtrList.H:61
Foam::psiuReactionThermo
Foam::psiuReactionThermo.
Definition: psiuReactionThermo.H:55
thermo
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
Foam::XiGModels::basicSubGrid::read
virtual bool read(const dictionary &XiGProperties)
Update properties from given dictionary.
Foam::XiGModel
Base-class for all Xi generation models used by the b-Xi combustion model. See Technical Report SH/RE...
Definition: XiGModel.H:56
Su
zeroField Su
Definition: alphaSuSp.H:1
Foam::XiGModels::basicSubGrid::Db
virtual tmp< volScalarField > Db() const
Return the flame diffusivity.
Foam::dictionary
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Definition: dictionary.H:123
Foam
Namespace for OpenFOAM.
Definition: atmBoundaryLayer.C:33
Foam::XiGModels::basicSubGrid::G
virtual tmp< volScalarField > G() const
Return the flame-wrinkling generation rate.
Foam::XiGModels::basicSubGrid::~basicSubGrid
virtual ~basicSubGrid()
Destructor.
Foam::autoPtr
Pointer management similar to std::unique_ptr, with some additional methods and type checking.
Definition: HashPtrTable.H:53
Foam::compressible::RASModel
RASModel< EddyDiffusivity< turbulenceModel > > RASModel
Definition: turbulentFluidThermoModel.H:66
Foam::XiGModels::basicSubGrid
Basic sub-grid obstacle flame-wrinkling generation rate coefficient model. Details supplied by J Putt...
Definition: basicXiSubG.H:71
Foam::GeometricField< scalar, fvPatchField, volMesh >
Foam::XiGModels::basicSubGrid::TypeName
TypeName("basicSubGridG")
Runtime type information.
XiGModel.H