src/lagrangian/dieselSpray/injector/definedInjector/definedInjector.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2008 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 Class
  26     Foam::definedInjector
  27
  28 Description
  29     User specified MFR vs time and velocity vs time
  30
  31 SourceFiles
  32     definedInjectorI.H
  33     definedInjector.C
  34
  35 \*---------------------------------------------------------------------------*/
  36
  37 #ifndef definedInjector_H
  38 #define definedInjector_H
  39
  40 #include "injectorType.H"
  41 #include "vector.H"
  42
  43 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  44
  45 namespace Foam
  46 {
  47
  48 /*---------------------------------------------------------------------------*\
  49                            Class definedInjector Declaration
  50 \*---------------------------------------------------------------------------*/
  51
  52 class definedInjector
  53 :
  54     public injectorType
  55 {
  56
  57 private:
  58
  59     // Private data
  60
  61         typedef VectorSpace<Vector<scalar>, scalar, 2> pair;
  62
  63         dictionary propsDict_;
  64
  65         vector position_;
  66         vector direction_;
  67         scalar d_;
  68         scalar mass_;
  69         scalar T_;
  70         label nParcels_;
  71         scalarField X_;
  72         List<pair> massFlowRateProfile_;                         // aj
  73         List<pair> velocityProfile_;                    // aj
  74         List<pair> injectionPressureProfile_;
  75         List<pair> CdProfile_;
  76         List<pair> TProfile_;
  77         scalar averageParcelMass_;
  78
  79         bool pressureIndependentVelocity_;
  80
  81         //- two orthogonal vectors that are also orthogonal
  82         //  to the injection direction
  83         vector tangentialInjectionVector1_, tangentialInjectionVector2_;
  84
  85
  86     // Private Member Functions
  87
  88         //- Disallow default bitwise copy construct
  89         definedInjector(const definedInjector&);
  90
  91         //- Disallow default bitwise assignment
  92         void operator=(const definedInjector&);
  93
  94         //- Create two vectors orthonoal to each other
  95         //  and the injection vector
  96         void setTangentialVectors();
  97
  98         //- Return the fraction of the total injected liquid
  99         scalar fractionOfInjection(const scalar time) const;
 100
 101         //- Return the instantaneous injection velocity
 102         scalar injectionVelocity(const scalar) const;
 103
 104 public:
 105
 106     //- Runtime type information
 107     TypeName("definedInjector");
 108
 109
 110     // Constructors
 111
 112         //- Construct from components
 113         definedInjector
 114         (
 115             const Time& t,
 116             const dictionary& dict
 117         );
 118
 119
 120     // Destructor
 121
 122         ~definedInjector();
 123
 124
 125     // Member Functions
 126
 127         //- Return number of particles to inject
 128         label nParcelsToInject
 129         (
 130             const scalar t0,
 131             const scalar t1
 132         ) const;
 133
 134         //- Return the injection position
 135         const vector position() const;
 136
 137         //- Return the injection position
 138         vector position
 139         (
 140             const scalar time,
 141             const bool twoD,
 142             const scalar angleOfWedge,
 143             const vector& axisOfSymmetry,
 144             const vector& axisOfWedge,
 145             const vector& axisOfWedgeNormal,
 146             Random& rndGen
 147         ) const;
 148
 149         //- Return the injector diameter
 150         scalar d() const;
 151
 152         //- Return the injection direction
 153         const vector& direction() const;
 154
 155         //- Return the mass of the injected particle
 156         scalar mass
 157         (
 158             const scalar t0,
 159             const scalar t1,
 160             const bool twoD,
 161             const scalar angleOfWedge
 162         ) const;
 163
 164         //- Return the mass injected by the injector
 165         scalar mass() const;
 166
 167         //- Return the fuel mass fractions of the injected particle
 168         const scalarField& X() const;
 169
 170         //- Return the temperature profile of the injected particle
 171         List<pair> T() const;
 172
 173         //- Return the temperature of the injected particle
 174         scalar T(const scalar time) const;
 175
 176         //- Return the start-of-injection time
 177         scalar tsoi() const;
 178
 179         //- Return the end-of-injection time
 180         scalar teoi() const;
 181
 182         //- Return the injected liquid mass
 183         scalar injectedMass(const scalar t) const;
 184
 185         List<pair> massFlowRateProfile() const
 186         {
 187             return massFlowRateProfile_;
 188         }
 189
 190         scalar massFlowRate(const scalar time) const;
 191
 192         List<pair> injectionPressureProfile() const
 193         {
 194             return injectionPressureProfile_;
 195         }
 196
 197         scalar injectionPressure(const scalar time) const;
 198
 199         List<pair> velocityProfile() const
 200         {
 201             return velocityProfile_;
 202         }
 203
 204         scalar velocity(const scalar time) const;
 205
 206         List<pair> CdProfile() const
 207         {
 208             return CdProfile_;
 209         }
 210
 211         scalar Cd(const scalar time) const;
 212
 213         void correctProfiles
 214         (
 215             const liquidMixture& fuel,
 216             const scalar referencePressure
 217         );
 218
 219         bool pressureIndependentVelocity() const
 220         {
 221             return pressureIndependentVelocity_;
 222         }
 223
 224 };
 225
 226
 227 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 228
 229 } // End namespace Foam
 230
 231 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 232
 233 #endif
 234
 235 // ************************************************************************* //