src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcelI.H

   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 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  28
  29 template<class ParcelType>
  30 inline Foam::ReactingMultiphaseParcel<ParcelType>::constantProperties::
  31 constantProperties
  32 (
  33     const dictionary& parentDict
  34 )
  35 :
  36     ReactingParcel<ParcelType>::constantProperties(parentDict),
  37     LDevol_(dimensionedScalar(this->dict().lookup("LDevol")).value()),
  38     hRetentionCoeff_
  39     (
  40         dimensionedScalar(this->dict().lookup("hRetentionCoeff")).value()
  41     )
  42 {
  43     if ((hRetentionCoeff_ < 0) || (hRetentionCoeff_ > 1))
  44     {
  45         FatalErrorIn
  46         (
  47             "ReactingMultiphaseParcel<ParcelType>::constantProperties::"
  48             "constantProperties"
  49         )   << "hRetentionCoeff must be in the range 0 to 1" << nl
  50             << exit(FatalError) << endl;
  51     }
  52
  53     hRetentionCoeff_ = max(1e-06, hRetentionCoeff_);
  54 }
  55
  56
  57 template<class ParcelType>
  58 inline Foam::ReactingMultiphaseParcel<ParcelType>::trackData::trackData
  59 (
  60     ReactingMultiphaseCloud<ParcelType>& cloud,
  61     const constantProperties& constProps,
  62     const interpolation<scalar>& rhoInterp,
  63     const interpolation<vector>& UInterp,
  64     const interpolation<scalar>& muInterp,
  65     const interpolation<scalar>& TInterp,
  66     const interpolation<scalar>& CpInterp,
  67     const interpolation<scalar>& pInterp,
  68     const vector& g
  69 )
  70 :
  71     ReactingParcel<ParcelType>::trackData
  72     (
  73         cloud,
  74         constProps,
  75         rhoInterp,
  76         UInterp,
  77         muInterp,
  78         TInterp,
  79         CpInterp,
  80         pInterp,
  81         g
  82     ),
  83     cloud_(cloud),
  84     constProps_(constProps)
  85 {}
  86
  87
  88 template <class ParcelType>
  89 inline Foam::ReactingMultiphaseParcel<ParcelType>::ReactingMultiphaseParcel
  90 (
  91     ReactingMultiphaseCloud<ParcelType>& owner,
  92     const vector& position,
  93     const label cellI
  94 )
  95 :
  96     ReactingParcel<ParcelType>(owner, position, cellI),
  97     YGas_(0),
  98     YLiquid_(0),
  99     YSolid_(0),
 100     canCombust_(false)
 101 {}
 102
 103
 104 template<class ParcelType>
 105 inline Foam::ReactingMultiphaseParcel<ParcelType>::ReactingMultiphaseParcel
 106 (
 107     ReactingMultiphaseCloud<ParcelType>& owner,
 108     const vector& position,
 109     const label cellI,
 110     const label typeId,
 111     const scalar nParticle0,
 112     const scalar d0,
 113     const vector& U0,
 114     const scalarField& Y0,
 115     const scalarField& YGas0,
 116     const scalarField& YLiquid0,
 117     const scalarField& YSolid0,
 118     const constantProperties& constProps
 119 )
 120 :
 121     ReactingParcel<ParcelType>
 122     (
 123         owner,
 124         position,
 125         cellI,
 126         typeId,
 127         nParticle0,
 128         d0,
 129         U0,
 130         Y0,
 131         constProps
 132     ),
 133     YGas_(YGas0),
 134     YLiquid_(YLiquid0),
 135     YSolid_(YSolid0),
 136     canCombust_(false)
 137 {}
 138
 139
 140 // * * * * * * * * * constantProperties Member Functions * * * * * * * * * * //
 141
 142 template<class ParcelType>
 143 inline Foam::scalar
 144 Foam::ReactingMultiphaseParcel<ParcelType>::constantProperties::LDevol() const
 145 {
 146     return LDevol_;
 147 }
 148
 149
 150 template<class ParcelType>
 151 inline Foam::scalar
 152 Foam::ReactingMultiphaseParcel<ParcelType>::constantProperties::
 153 hRetentionCoeff() const
 154 {
 155     return hRetentionCoeff_;
 156 }
 157
 158
 159 // * * * * * * * * * * * trackData Member Functions  * * * * * * * * * * * * //
 160
 161 template<class ParcelType>
 162 inline Foam::ReactingMultiphaseCloud<ParcelType>&
 163 Foam::ReactingMultiphaseParcel<ParcelType>::trackData::cloud()
 164 {
 165     return cloud_;
 166 }
 167
 168
 169 template<class ParcelType>
 170 inline const typename Foam::ReactingMultiphaseParcel<ParcelType>::
 171 constantProperties&
 172 Foam::ReactingMultiphaseParcel<ParcelType>::trackData::constProps() const
 173 {
 174     return constProps_;
 175 }
 176
 177
 178 // * * * * * * * * * * ThermoParcel Member Functions * * * * * * * * * * * * //
 179
 180 template<class ParcelType>
 181 inline const Foam::scalarField& Foam::ReactingMultiphaseParcel<ParcelType>::
 182 YGas() const
 183 {
 184     return YGas_;
 185 }
 186
 187
 188 template<class ParcelType>
 189 inline const Foam::scalarField& Foam::ReactingMultiphaseParcel<ParcelType>::
 190 YLiquid() const
 191 {
 192     return YLiquid_;
 193 }
 194
 195
 196 template<class ParcelType>
 197 inline const Foam::scalarField& Foam::ReactingMultiphaseParcel<ParcelType>::
 198 YSolid() const
 199 {
 200     return YSolid_;
 201 }
 202
 203
 204 template<class ParcelType>
 205 inline bool Foam::ReactingMultiphaseParcel<ParcelType>::canCombust() const
 206 {
 207     return canCombust_;
 208 }
 209
 210
 211 template<class ParcelType>
 212 inline Foam::scalarField& Foam::ReactingMultiphaseParcel<ParcelType>::YGas()
 213 {
 214     return YGas_;
 215 }
 216
 217
 218 template<class ParcelType>
 219 inline Foam::scalarField& Foam::ReactingMultiphaseParcel<ParcelType>::YLiquid()
 220 {
 221     return YLiquid_;
 222 }
 223
 224
 225 template<class ParcelType>
 226 inline Foam::scalarField& Foam::ReactingMultiphaseParcel<ParcelType>::YSolid()
 227 {
 228     return YSolid_;
 229 }
 230
 231
 232 template<class ParcelType>
 233 inline bool& Foam::ReactingMultiphaseParcel<ParcelType>::canCombust()
 234 {
 235     return canCombust_;
 236 }
 237
 238 // ************************************************************************* //