src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcelI.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
   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
  26 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  27
  28 template<class ParcelType>
  29 inline
  30 Foam::ReactingParcel<ParcelType>::constantProperties::constantProperties()
  31 :
  32     ParcelType::constantProperties(),
  33     pMin_(0.0),
  34     constantVolume_(false),
  35     Tvap_(0.0),
  36     Tbp_(0.0)
  37 {}
  38
  39
  40 template<class ParcelType>
  41 inline Foam::ReactingParcel<ParcelType>::constantProperties::constantProperties
  42 (
  43     const constantProperties& cp
  44 )
  45 :
  46     ParcelType::constantProperties(cp),
  47     pMin_(cp.pMin_),
  48     constantVolume_(cp.constantVolume_),
  49     Tvap_(cp.Tvap_),
  50     Tbp_(cp.Tbp_)
  51 {}
  52
  53
  54 template<class ParcelType>
  55 inline Foam::ReactingParcel<ParcelType>::constantProperties::constantProperties
  56 (
  57     const dictionary& parentDict,
  58     const bool readFields
  59 )
  60 :
  61     ParcelType::constantProperties(parentDict, readFields),
  62     pMin_(0.0),
  63     constantVolume_(false),
  64     Tvap_(0.0),
  65     Tbp_(0.0)
  66 {
  67     if (readFields)
  68     {
  69         this->dict().lookup("pMin") >> pMin_;
  70         this->dict().lookup("constantVolume") >> constantVolume_;
  71         this->dict().lookup("Tvap") >> Tvap_;
  72         this->dict().lookup("Tbp") >> Tbp_;
  73     }
  74 }
  75
  76
  77 template<class ParcelType>
  78 inline Foam::ReactingParcel<ParcelType>::constantProperties::constantProperties
  79 (
  80     const label parcelTypeId,
  81     const scalar rhoMin,
  82     const scalar rho0,
  83     const scalar minParticleMass,
  84     const scalar youngsModulus,
  85     const scalar poissonsRatio,
  86     const scalar T0,
  87     const scalar TMin,
  88     const scalar Cp0,
  89     const scalar epsilon0,
  90     const scalar f0,
  91     const scalar Pr,
  92     const scalar pMin,
  93     const Switch& constantVolume,
  94     const scalar Tvap,
  95     const scalar Tbp
  96 )
  97 :
  98     ParcelType::constantProperties
  99     (
 100         parcelTypeId,
 101         rhoMin,
 102         rho0,
 103         minParticleMass,
 104         youngsModulus,
 105         poissonsRatio,
 106         T0,
 107         TMin,
 108         Cp0,
 109         epsilon0,
 110         f0,
 111         Pr
 112     ),
 113     pMin_(pMin),
 114     constantVolume_(constantVolume),
 115     Tvap_(Tvap),
 116     Tbp_(Tbp)
 117 {}
 118
 119
 120 template<class ParcelType>
 121 inline Foam::ReactingParcel<ParcelType>::ReactingParcel
 122 (
 123     const polyMesh& mesh,
 124     const vector& position,
 125     const label cellI,
 126     const label tetFaceI,
 127     const label tetPtI
 128 )
 129 :
 130     ParcelType(mesh, position, cellI, tetFaceI, tetPtI),
 131     mass0_(0.0),
 132     Y_(0),
 133     pc_(0.0)
 134 {}
 135
 136
 137 template<class ParcelType>
 138 inline Foam::ReactingParcel<ParcelType>::ReactingParcel
 139 (
 140     const polyMesh& mesh,
 141     const vector& position,
 142     const label cellI,
 143     const label tetFaceI,
 144     const label tetPtI,
 145     const label typeId,
 146     const scalar nParticle0,
 147     const scalar d0,
 148     const scalar dTarget0,
 149     const vector& U0,
 150     const vector& f0,
 151     const vector& angularMomentum0,
 152     const vector& torque0,
 153     const scalarField& Y0,
 154     const constantProperties& constProps
 155 )
 156 :
 157     ParcelType
 158     (
 159         mesh,
 160         position,
 161         cellI,
 162         tetFaceI,
 163         tetPtI,
 164         typeId,
 165         nParticle0,
 166         d0,
 167         dTarget0,
 168         U0,
 169         f0,
 170         angularMomentum0,
 171         torque0,
 172         constProps
 173     ),
 174     mass0_(0.0),
 175     Y_(Y0),
 176     pc_(0.0)
 177 {
 178     // Set initial parcel mass
 179     mass0_ = this->mass();
 180 }
 181
 182
 183 // * * * * * * * * * constantProperties Member Functions * * * * * * * * * * //
 184
 185 template<class ParcelType>
 186 inline Foam::scalar
 187 Foam::ReactingParcel<ParcelType>::constantProperties::pMin() const
 188 {
 189     return pMin_;
 190 }
 191
 192
 193 template<class ParcelType>
 194 inline Foam::Switch
 195 Foam::ReactingParcel<ParcelType>::constantProperties::constantVolume() const
 196 {
 197     return constantVolume_;
 198 }
 199
 200
 201 template<class ParcelType>
 202 inline Foam::scalar
 203 Foam::ReactingParcel<ParcelType>::constantProperties::Tvap() const
 204 {
 205     return Tvap_;
 206 }
 207
 208
 209 template<class ParcelType>
 210 inline Foam::scalar
 211 Foam::ReactingParcel<ParcelType>::constantProperties::Tbp() const
 212 {
 213     return Tbp_;
 214 }
 215
 216
 217 // * * * * * * * * * * ThermoParcel Member Functions * * * * * * * * * * * * //
 218
 219 template<class ParcelType>
 220 inline Foam::scalar Foam::ReactingParcel<ParcelType>::mass0() const
 221 {
 222     return mass0_;
 223 }
 224
 225
 226 template<class ParcelType>
 227 inline const Foam::scalarField& Foam::ReactingParcel<ParcelType>::Y() const
 228 {
 229     return Y_;
 230 }
 231
 232
 233 template<class ParcelType>
 234 inline Foam::scalar Foam::ReactingParcel<ParcelType>::pc() const
 235 {
 236     return pc_;
 237 }
 238
 239
 240 template<class ParcelType>
 241 inline Foam::scalar& Foam::ReactingParcel<ParcelType>::pc()
 242 {
 243     return pc_;
 244 }
 245
 246
 247 template<class ParcelType>
 248 inline Foam::scalar& Foam::ReactingParcel<ParcelType>::mass0()
 249 {
 250     return mass0_;
 251 }
 252
 253
 254 template<class ParcelType>
 255 inline Foam::scalarField& Foam::ReactingParcel<ParcelType>::Y()
 256 {
 257     return Y_;
 258 }
 259
 260
 261 // ************************************************************************* //