src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/KinematicLookupTableInjection/KinematicLookupTableInjection.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2008-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 "KinematicLookupTableInjection.H"
  28 #include "scalarIOList.H"
  29
  30 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  31
  32 template<class CloudType>
  33 Foam::label Foam::KinematicLookupTableInjection<CloudType>::INPUT_FILE_COLS = 9;
  34
  35 // * * * * * * * * * * * * Protected Member Functions  * * * * * * * * * * * //
  36
  37 template<class CloudType>
  38 Foam::label Foam::KinematicLookupTableInjection<CloudType>::parcelsToInject
  39 (
  40     const scalar time0,
  41     const scalar time1
  42 ) const
  43 {
  44     if ((time0 >= 0.0) && (time0 < duration_))
  45     {
  46         return round(injectorCells_.size()*(time1 - time0)*nParcelsPerSecond_);
  47     }
  48     else
  49     {
  50         return 0;
  51     }
  52 }
  53
  54
  55 template<class CloudType>
  56 Foam::scalar Foam::KinematicLookupTableInjection<CloudType>::volumeToInject
  57 (
  58     const scalar time0,
  59     const scalar time1
  60 ) const
  61 {
  62     scalar volume = 0.0;
  63     if ((time0 >= 0.0) && (time0 < duration_))
  64     {
  65         forAll(mDot_, injectorI)
  66         {
  67             volume += mDot_[injectorI]/rho_[injectorI]*(time1 - time0);
  68         }
  69     }
  70
  71     return volume;
  72 }
  73
  74
  75 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  76
  77 template<class CloudType>
  78 Foam::KinematicLookupTableInjection<CloudType>::KinematicLookupTableInjection
  79 (
  80     const dictionary& dict,
  81     CloudType& owner
  82 )
  83 :
  84     InjectionModel<CloudType>(dict, owner, typeName),
  85     inputFileName_(this->coeffDict().lookup("inputFile")),
  86     duration_(readScalar(this->coeffDict().lookup("duration"))),
  87     nParcelsPerSecond_
  88     (
  89         readScalar(this->coeffDict().lookup("parcelsPerSecond"))
  90     ),
  91     x_(0),
  92     U_(0),
  93     d_(0),
  94     rho_(0),
  95     mDot_(0),
  96     injectorCells_(0)
  97 {
  98     scalarListIOList injectorData
  99     (
 100         IOobject
 101         (
 102             inputFileName_,
 103             owner.db().time().constant(),
 104             owner.db(),
 105             IOobject::MUST_READ,
 106             IOobject::NO_WRITE
 107         )
 108     );
 109
 110     x_.setSize(injectorData.size());
 111     U_.setSize(injectorData.size());
 112     d_.setSize(injectorData.size());
 113     rho_.setSize(injectorData.size());
 114     mDot_.setSize(injectorData.size());
 115
 116     // Populate lists
 117     forAll(injectorData, injectorI)
 118     {
 119         if (injectorData[injectorI].size() != INPUT_FILE_COLS)
 120         {
 121             FatalErrorIn
 122             (
 123                 "KinematicLookupTableInjection"
 124                 "("
 125                     "const dictionary&,"
 126                     "CloudType& owner"
 127                 ")"
 128             )   << "Incorrect number of entries in injector specification "
 129                 << "- found " << injectorData[injectorI].size()
 130                 << ", expected " << INPUT_FILE_COLS << ":" << nl
 131                 << "    x0 x1 x2 u0 u1 u2 d rho mDot " << nl
 132                 << exit(FatalError);
 133         }
 134         x_[injectorI].component(0) = injectorData[injectorI][0];
 135         x_[injectorI].component(1) = injectorData[injectorI][1];
 136         x_[injectorI].component(2) = injectorData[injectorI][2];
 137         U_[injectorI].component(0) = injectorData[injectorI][3];
 138         U_[injectorI].component(1) = injectorData[injectorI][4];
 139         U_[injectorI].component(2) = injectorData[injectorI][5];
 140         d_[injectorI] = injectorData[injectorI][6];
 141         rho_[injectorI] = injectorData[injectorI][7];
 142         mDot_[injectorI] = injectorData[injectorI][8];
 143     }
 144
 145     // Set/cache the injector cells
 146     injectorCells_.setSize(injectorData.size());
 147     forAll(x_, injectorI)
 148     {
 149         this->findCellAtPosition(injectorCells_[injectorI], x_[injectorI]);
 150     }
 151
 152     // Determine volume of particles to inject
 153     this->volumeTotal_ = 0.0;
 154     forAll(mDot_, injectorI)
 155     {
 156         this->volumeTotal_ += mDot_[injectorI]/rho_[injectorI];
 157     }
 158     this->volumeTotal_ *= duration_;
 159 }
 160
 161
 162 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 163
 164 template<class CloudType>
 165 Foam::KinematicLookupTableInjection<CloudType>::~KinematicLookupTableInjection()
 166 {}
 167
 168
 169 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 170
 171 template<class CloudType>
 172 bool Foam::KinematicLookupTableInjection<CloudType>::active() const
 173 {
 174     return true;
 175 }
 176
 177
 178 template<class CloudType>
 179 Foam::scalar Foam::KinematicLookupTableInjection<CloudType>::timeEnd() const
 180 {
 181     return this->SOI_ + duration_;
 182 }
 183
 184
 185 template<class CloudType>
 186 void Foam::KinematicLookupTableInjection<CloudType>::setPositionAndCell
 187 (
 188     const label parcelI,
 189     const label nParcels,
 190     const scalar time,
 191     vector& position,
 192     label& cellOwner
 193 )
 194 {
 195     label injectorI = parcelI*injectorCells_.size()/nParcels;
 196
 197     position = x_[injectorI];
 198     cellOwner = injectorCells_[injectorI];
 199 }
 200
 201
 202 template<class CloudType>
 203 void Foam::KinematicLookupTableInjection<CloudType>::setProperties
 204 (
 205     const label parcelI,
 206     const label nParcels,
 207     const scalar,
 208     typename CloudType::parcelType& parcel
 209 )
 210 {
 211     label injectorI = parcelI*injectorCells_.size()/nParcels;
 212
 213     // set particle velocity
 214     parcel.U() = U_[injectorI];
 215
 216     // set particle diameter
 217     parcel.d() = d_[injectorI];
 218
 219     // set particle density
 220     parcel.rho() = rho_[injectorI];
 221 }
 222
 223
 224 template<class CloudType>
 225 bool Foam::KinematicLookupTableInjection<CloudType>::fullyDescribed() const
 226 {
 227     return true;
 228 }
 229
 230
 231 template<class CloudType>
 232 bool Foam::KinematicLookupTableInjection<CloudType>::validInjection
 233 (
 234     const label
 235 )
 236 {
 237     return true;
 238 }
 239
 240
 241 // ************************************************************************* //