src/lagrangian/intermediate/submodels/Thermodynamic/InjectionModel/ThermoLookupTableInjection/ThermoLookupTableInjection.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright held by original author
   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 "ThermoLookupTableInjection.H"
  28 #include "scalarIOList.H"
  29
  30 // * * * * * * * * * * * * Protected Member Functions  * * * * * * * * * * * //
  31
  32 template<class CloudType>
  33 Foam::label Foam::ThermoLookupTableInjection<CloudType>::parcelsToInject
  34 (
  35     const scalar time0,
  36     const scalar time1
  37 ) const
  38 {
  39     if ((time0 >= 0.0) && (time0 < duration_))
  40     {
  41         return round(injectorCells_.size()*(time1 - time0)*nParcelsPerSecond_);
  42     }
  43     else
  44     {
  45         return 0;
  46     }
  47 }
  48
  49
  50 template<class CloudType>
  51 Foam::scalar Foam::ThermoLookupTableInjection<CloudType>::volumeToInject
  52 (
  53     const scalar time0,
  54     const scalar time1
  55 ) const
  56 {
  57     scalar volume = 0.0;
  58     if ((time0 >= 0.0) && (time0 < duration_))
  59     {
  60         forAll(injectors_, i)
  61         {
  62             volume += injectors_[i].mDot()/injectors_[i].rho()*(time1 - time0);
  63         }
  64     }
  65
  66     return volume;
  67 }
  68
  69
  70 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  71
  72 template<class CloudType>
  73 Foam::ThermoLookupTableInjection<CloudType>::ThermoLookupTableInjection
  74 (
  75     const dictionary& dict,
  76     CloudType& owner
  77 )
  78 :
  79     InjectionModel<CloudType>(dict, owner, typeName),
  80     inputFileName_(this->coeffDict().lookup("inputFile")),
  81     duration_(readScalar(this->coeffDict().lookup("duration"))),
  82     nParcelsPerSecond_
  83     (
  84         readScalar(this->coeffDict().lookup("parcelsPerSecond"))
  85     ),
  86     injectors_
  87     (
  88         IOobject
  89         (
  90             inputFileName_,
  91             owner.db().time().constant(),
  92             owner.db(),
  93             IOobject::MUST_READ,
  94             IOobject::NO_WRITE
  95         )
  96     ),
  97     injectorCells_(0)
  98 {
  99     // Set/cache the injector cells
 100     injectorCells_.setSize(injectors_.size());
 101     forAll(injectors_, i)
 102     {
 103         this->findCellAtPosition(injectorCells_[injectorI], injectors_[i].x());
 104     }
 105
 106     // Determine volume of particles to inject
 107     this->volumeTotal_ = 0.0;
 108     forAll(injectors_, i)
 109     {
 110         this->volumeTotal_ += injectors_[i].mDot()/injectors_[i].rho();
 111     }
 112     this->volumeTotal_ *= duration_;
 113 }
 114
 115
 116 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 117
 118 template<class CloudType>
 119 Foam::ThermoLookupTableInjection<CloudType>::~ThermoLookupTableInjection()
 120 {}
 121
 122
 123 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 124
 125 template<class CloudType>
 126 bool Foam::ThermoLookupTableInjection<CloudType>::active() const
 127 {
 128     return true;
 129 }
 130
 131
 132 template<class CloudType>
 133 Foam::scalar Foam::ThermoLookupTableInjection<CloudType>::timeEnd() const
 134 {
 135     return this->SOI_ + duration_;
 136 }
 137
 138
 139 template<class CloudType>
 140 void Foam::ThermoLookupTableInjection<CloudType>::setPositionAndCell
 141 (
 142     const label parcelI,
 143     const label nParcels,
 144     const scalar time,
 145     vector& position,
 146     label& cellOwner
 147 )
 148 {
 149     label injectorI = parcelI*injectorCells_.size()/nParcels;
 150
 151     position = injectors_[injectorI].x();
 152     cellOwner = injectorCells_[injectorI];
 153 }
 154
 155
 156 template<class CloudType>
 157 void Foam::ThermoLookupTableInjection<CloudType>::setProperties
 158 (
 159     const label parcelI,
 160     const label nParcels,
 161     const scalar,
 162     typename CloudType::parcelType* pPtr
 163 )
 164 {
 165     label injectorI = parcelI*injectorCells_.size()/nParcels;
 166
 167     // set particle velocity
 168     parcel.U() = injectors_[injectorI].U();
 169
 170     // set particle diameter
 171     parcel.d() = injectors_[injectorI].d();
 172
 173     // set particle density
 174     parcel.rho() = injectors_[injectorI].rho();
 175
 176     // set particle temperature
 177     parcel.T() = injectors_[injectorI].T();
 178
 179     // set particle specific heat capacity
 180     parcel.cp() = injectors_[injectorI].cp();
 181 }
 182
 183
 184 template<class CloudType>
 185 bool Foam::ThermoLookupTableInjection<CloudType>::fullyDescribed() const
 186 {
 187     return true;
 188 }
 189
 190
 191 template<class CloudType>
 192 bool Foam::ThermoLookupTableInjection<CloudType>::validInjection(const label)
 193 {
 194     return true;
 195 }
 196
 197
 198 // ************************************************************************* //