FreeFOAM The Cross-Platform CFD Toolkit
laminar.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) 1991-2010 OpenCFD Ltd.
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 Class
25  Foam::compressible::laminar
26 
27 Description
28  Turbulence model for laminar compressible flow.
29 
30 SourceFiles
31  laminar.C
32 
33 \*---------------------------------------------------------------------------*/
34 
35 #ifndef compressibleLaminar_H
36 #define compressibleLaminar_H
37 
39 
40 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
41 
42 namespace Foam
43 {
44 namespace compressible
45 {
46 
47 /*---------------------------------------------------------------------------*\
48  Class laminar Declaration
49 \*---------------------------------------------------------------------------*/
50 
51 class laminar
52 :
53  public turbulenceModel
54 {
55 
56 public:
57 
58  //- Runtime type information
59  TypeName("laminar");
60 
61  // Constructors
62 
63  //- Construct from components
64  laminar
65  (
66  const volScalarField& rho,
67  const volVectorField& U,
68  const surfaceScalarField& phi,
69  const basicThermo& thermophysicalModel
70  );
71 
72 
73  // Selectors
74 
75  //- Return a reference to the selected turbulence model
76  static autoPtr<laminar> New
77  (
78  const volScalarField& rho,
79  const volVectorField& U,
80  const surfaceScalarField& phi,
81  const basicThermo& thermophysicalModel
82  );
83 
84 
85  //- Destructor
86  virtual ~laminar()
87  {}
88 
89 
90  // Member Functions
91 
92  //- Return the turbulence viscosity, i.e. 0 for laminar flow
93  virtual tmp<volScalarField> mut() const;
94 
95  //- Return the effective viscosity, i.e. the laminar viscosity
96  virtual tmp<volScalarField> muEff() const
97  {
98  return tmp<volScalarField>(new volScalarField("muEff", mu()));
99  }
100 
101  //- Return the effective turbulent thermal diffusivity,
102  // i.e. the laminar thermal diffusivity
104  {
105  return tmp<volScalarField>(new volScalarField("alphaEff", alpha()));
106  }
107 
108  //- Return the turbulence kinetic energy, i.e. 0 for laminar flow
109  virtual tmp<volScalarField> k() const;
110 
111  //- Return the turbulence kinetic energy dissipation rate,
112  // i.e. 0 for laminar flow
113  virtual tmp<volScalarField> epsilon() const;
114 
115  //- Return the Reynolds stress tensor, i.e. 0 for laminar flow
116  virtual tmp<volSymmTensorField> R() const;
117 
118  //- Return the effective stress tensor, i.e. the laminar stress
119  virtual tmp<volSymmTensorField> devRhoReff() const;
120 
121  //- Return the source term for the momentum equation
122  virtual tmp<fvVectorMatrix> divDevRhoReff(volVectorField& U) const;
123 
124  //- Correct the laminar viscosity
125  virtual void correct();
126 
127  //- Read turbulenceProperties dictionary
128  virtual bool read();
129 };
130 
131 
132 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
133 
134 } // End namespace incompressible
135 } // End namespace Foam
136 
137 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
138 
139 #endif
140 
141 // ************************ vim: set sw=4 sts=4 et: ************************ //