OpenFOAM logo
Open Source CFD Toolkit

swirlInjector.H

Go to the documentation of this file.
00001 /*---------------------------------------------------------------------------*\
00002   =========                 |
00003   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
00004    \\    /   O peration     |
00005     \\  /    A nd           | Copyright (C) 1991-2005 OpenCFD Ltd.
00006      \\/     M anipulation  |
00007 -------------------------------------------------------------------------------
00008 License
00009     This file is part of OpenFOAM.
00010 
00011     OpenFOAM is free software; you can redistribute it and/or modify it
00012     under the terms of the GNU General Public License as published by the
00013     Free Software Foundation; either version 2 of the License, or (at your
00014     option) any later version.
00015 
00016     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
00017     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
00018     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
00019     for more details.
00020 
00021     You should have received a copy of the GNU General Public License
00022     along with OpenFOAM; if not, write to the Free Software Foundation,
00023     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
00024 
00025 Class
00026     swirlInjector
00027 
00028 Description
00029 
00030 SourceFiles
00031     swirlInjectorI.H
00032     swirlInjector.C
00033     swirlInjectorIO.C
00034 
00035 \*---------------------------------------------------------------------------*/
00036 
00037 #ifndef swirlInjector_H
00038 #define swirlInjector_H
00039 
00040 #include "injectorType.H"
00041 #include "vector.H"
00042 
00043 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00044 
00045 namespace Foam
00046 {
00047 
00048 /*---------------------------------------------------------------------------*\
00049                            Class swirlInjector Declaration
00050 \*---------------------------------------------------------------------------*/
00051 
00052 class swirlInjector
00053 :
00054     public injectorType
00055 {
00056 
00057 private:
00058 
00059     // Private data
00060 
00061         typedef VectorSpace<Vector<scalar>, scalar, 2> pair;
00062 
00063         dictionary propsDict_;
00064 
00065         vector position_;
00066         vector direction_;
00067         scalar d_;
00068         scalar mass_;
00069         scalar injectionPressure_;
00070         scalar T_;
00071         label nParcels_;
00072         scalarField X_;
00073         List<pair> massFlowRateProfile_;
00074         List<pair> injectionPressureProfile_;
00075         List<pair> velocityProfile_;
00076         List<pair> CdProfile_;
00077         List<pair> TProfile_;
00078         scalar averageParcelMass_;
00079 
00080         bool pressureIndependentVelocity_;
00081 
00082         //- two orthogonal vectors that are also orthogonal
00083         //  to the injection direction
00084         vector tangentialInjectionVector1_, tangentialInjectionVector2_;
00085 
00086 
00087     // Private Member Functions
00088 
00089         //- Disallow default bitwise copy construct
00090         swirlInjector(const swirlInjector&);
00091 
00092         //- Disallow default bitwise assignment
00093         void operator=(const swirlInjector&);
00094 
00095         //- Integrate the mass-flow profile and set
00096         //  the integral value
00097         void setIntegralProfileValue();
00098 
00099         //- Create two vectors orthonoal to each other
00100         //  and the injection vector
00101         void setTangentialVectors();
00102 
00103         //- Return the fraction of the total injected liquid
00104         scalar fractionOfInjection(const scalar time) const;
00105 
00106 
00107 /*
00108 
00109         //- Return the average injection velocity
00110         scalar averageInjectionVelocityTime(const scalar liquidDensity) const;
00111 
00112         //- Return the instantaneous injection velocity
00113         scalar injectionVelocity(const scalar, const scalar liquidDensity) const;
00114 
00115 */
00116 
00117 public:
00118 
00119     //- Runtime type information
00120     TypeName("swirlInjector");
00121 
00122 
00123     // Constructors
00124 
00125         //- Construct from components
00126         swirlInjector
00127         (
00128             const Time& t,
00129             const dictionary& dict
00130         );
00131 
00132 
00133     // Destructor
00134 
00135         ~swirlInjector();
00136 
00137 
00138     // Member Functions
00139 
00140         //- Return number of particles to inject
00141         label nParcelsToInject
00142         (
00143             const scalar t0,
00144             const scalar t1
00145         ) const;
00146 
00147         //- Return the injection position
00148         const vector position() const;
00149 
00150         //- Return the injection position
00151         vector position
00152         (
00153             const scalar time,
00154             const bool twoD,
00155             const scalar angleOfWedge,
00156             const vector& axisOfSymmetry,
00157             const vector& axisOfWedge,
00158             const vector& axisOfWedgeNormal,
00159             Random& rndGen
00160         ) const;
00161     
00162         //- Return the injector diameter
00163         scalar d() const;
00164 
00165         //- Return the injection direction
00166         const vector& direction() const;
00167 
00168         //- Return the mass of the injected particle
00169         scalar mass
00170         (
00171             const scalar t0,
00172             const scalar t1,
00173             const bool twoD,
00174             const scalar angleOfWedge
00175         ) const;
00176 
00177         //- Return the mass injected by the injector
00178         scalar mass() const;
00179 
00180         //- Return the mass flow rate profile
00181         List<pair> massFlowRateProfile() const;
00182 
00183         //- Return the instantaneous mass flow rate profile
00184         scalar massFlowRate(const scalar time) const;
00185 
00186         //- Return the pressure injection profile
00187         List<pair> injectionPressureProfile() const;
00188 
00189         //- Return the instantaneous injection pressure
00190         scalar injectionPressure(const scalar time) const;
00191 
00192         //- Return the velocity injection profile
00193         List<pair> velocityProfile() const;
00194 
00195         //- Return the instantaneous velocity
00196         scalar velocity(const scalar time) const;
00197 
00198         //- Return the discharge coefficient    
00199         List<pair> CdProfile() const;
00200 
00201         //- Return the instantaneous discharge coefficient
00202         scalar Cd(const scalar time) const;
00203 
00204         //- Return the fuel mass fractions of the injected particle
00205         const scalarField& X() const;
00206 
00207         //- Return the temperatue profile of the injected parcel
00208         List<pair> T() const;
00209 
00210         //- Return the temperatue of the injected parcel
00211         scalar T(const scalar time) const;
00212 
00213         //- Return the start-of-injection time
00214         scalar tsoi() const;
00215 
00216         //- Return the end-of-injection time
00217         scalar teoi() const;
00218 
00219         //- Return the injected liquid mass
00220         scalar injectedMass(const scalar t) const;
00221 
00222         void correctProfiles
00223         (
00224             const liquidMixture& fuel,
00225             const scalar referencePressure
00226         );
00227         
00228         bool pressureIndependentVelocity() const
00229         {
00230             return pressureIndependentVelocity_;
00231         }
00232 };
00233 
00234 
00235 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00236 
00237 } // End namespace Foam
00238 
00239 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00240 
00241 #endif
00242 
00243 // ************************************************************************* //
For further information go to www.openfoam.org