FreeFOAM The Cross-Platform CFD Toolkit
DispersionRASModel.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 | 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 \*---------------------------------------------------------------------------*/
25 
27 
28 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
29 
30 template<class CloudType>
32 (
33  const dictionary& dict,
34  CloudType& owner
35 )
36 :
37  DispersionModel<CloudType>(dict, owner),
38  turbulence_
39  (
40  owner.mesh().objectRegistry::lookupObject<compressible::RASModel>
41  (
42  "RASProperties"
43  )
44  ),
45  kPtr_(NULL),
46  ownK_(false),
47  epsilonPtr_(NULL),
48  ownEpsilon_(false)
49 {}
50 
51 
52 // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
53 
54 template<class CloudType>
56 {
57  cacheFields(false);
58 }
59 
60 
61 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
62 
63 template<class CloudType>
65 {
66  if (store)
67  {
68  tmp<volScalarField> tk = this->turbulence().k();
69  if (tk.isTmp())
70  {
71  kPtr_ = tk.ptr();
72  ownK_ = true;
73  }
74  else
75  {
76  kPtr_ = tk.operator->();
77  ownK_ = false;
78  }
79 
80  tmp<volScalarField> tepsilon = this->turbulence().epsilon();
81  if (tepsilon.isTmp())
82  {
83  epsilonPtr_ = tepsilon.ptr();
84  ownEpsilon_ = true;
85  }
86  else
87  {
88  epsilonPtr_ = tepsilon.operator->();
89  ownEpsilon_ = false;
90  }
91  }
92  else
93  {
94  if (ownK_ && kPtr_)
95  {
96  delete kPtr_;
97  kPtr_ = NULL;
98  ownK_ = false;
99  }
100  if (ownEpsilon_ && epsilonPtr_)
101  {
102  delete epsilonPtr_;
103  epsilonPtr_ = NULL;
104  ownEpsilon_ = false;
105  }
106  }
107 }
108 
109 
110 // ************************ vim: set sw=4 sts=4 et: ************************ //