src/lagrangian/intermediate/submodels/Kinematic/DispersionModel/DispersionRASModel/DispersionRASModel.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 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 the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     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, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #include "DispersionRASModel.H"
  28
  29 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  30
  31 template<class CloudType>
  32 Foam::DispersionRASModel<CloudType>::DispersionRASModel
  33 (
  34     const dictionary& dict,
  35     CloudType& owner
  36 )
  37 :
  38     DispersionModel<CloudType>(dict, owner),
  39     turbulence_
  40     (
  41         owner.mesh().objectRegistry::lookupObject<compressible::RASModel>
  42         (
  43             "RASProperties"
  44         )
  45     ),
  46     kPtr_(NULL),
  47     ownK_(false),
  48     epsilonPtr_(NULL),
  49     ownEpsilon_(false)
  50 {}
  51
  52
  53 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
  54
  55 template<class CloudType>
  56 Foam::DispersionRASModel<CloudType>::~DispersionRASModel()
  57 {
  58     cacheFields(false);
  59 }
  60
  61
  62 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  63
  64 template<class CloudType>
  65 void Foam::DispersionRASModel<CloudType>::cacheFields(const bool store)
  66 {
  67     if (store)
  68     {
  69         tmp<volScalarField> tk = this->turbulence().k();
  70         if (tk.isTmp())
  71         {
  72             kPtr_ = tk.ptr();
  73             ownK_ = true;
  74         }
  75         else
  76         {
  77             kPtr_ = tk.operator->();
  78             ownK_ = false;
  79         }
  80
  81         tmp<volScalarField> tepsilon = this->turbulence().epsilon();
  82         if (tepsilon.isTmp())
  83         {
  84             epsilonPtr_ = tepsilon.ptr();
  85             ownEpsilon_ = true;
  86         }
  87         else
  88         {
  89             epsilonPtr_ = tepsilon.operator->();
  90             ownEpsilon_ = false;
  91         }
  92     }
  93     else
  94     {
  95         if (ownK_ && kPtr_)
  96         {
  97             delete kPtr_;
  98             ownK_ = false;
  99         }
 100         if (ownEpsilon_ && epsilonPtr_)
 101         {
 102             delete epsilonPtr_;
 103             ownEpsilon_ = false;
 104         }
 105     }
 106 }
 107
 108
 109 // ************************************************************************* //