BlobsSheetAtomization.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2011-2018 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 #include "BlobsSheetAtomization.H"
27 
28 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
29 
30 template<class CloudType>
32 (
33  const dictionary& dict,
34  CloudType& owner
35 )
36 :
37  AtomizationModel<CloudType>(dict, owner, typeName),
38  B_(readScalar(this->coeffDict().lookup("B"))),
39  angle_(readScalar(this->coeffDict().lookup("angle")))
40 {}
41 
42 
43 template<class CloudType>
45 (
47 )
48 :
50  B_(am.B_),
51  angle_(am.angle_)
52 {}
53 
54 
55 // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
56 
57 template<class CloudType>
59 {}
60 
61 
62 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
63 
64 template<class CloudType>
66 {
67  return 1.0;
68 }
69 
70 
71 template<class CloudType>
73 {
74  return false;
75 }
76 
77 
78 template<class CloudType>
80 (
81  const scalar dt,
82  scalar& d,
83  scalar& liquidCore,
84  scalar& tc,
85  const scalar rho,
86  const scalar mu,
87  const scalar sigma,
88  const scalar volFlowRate,
89  const scalar rhoAv,
90  const scalar Urel,
91  const vector& pos,
92  const vector& injectionPos,
93  const scalar pAmbient,
94  const scalar chi,
95  Random& rndGen
96 ) const
97 {
98  scalar lBU =
99  B_
100  * sqrt
101  (
102  rho*sigma*d*cos(angle_*constant::mathematical::pi/360.0)
103  / sqr(rhoAv*Urel)
104  );
105 
106  scalar pWalk = mag(pos - injectionPos);
107 
108  if (pWalk > lBU)
109  {
110  liquidCore = 0.0;
111  }
112 }
113 
114 
115 // ************************************************************************* //
#define readScalar
Definition: doubleScalar.C:38
dictionary dict
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:137
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar sqrt(const dimensionedScalar &ds)
Templated atomization model class.
Definition: SprayCloud.H:47
BlobsSheetAtomization(const dictionary &dict, CloudType &cloud)
Construct from dictionary.
virtual ~BlobsSheetAtomization()
Destructor.
stressControl lookup("compactNormalStress") >> compactNormalStress
dimensionedScalar cos(const dimensionedScalar &ds)
virtual scalar initLiquidCore() const
Initial value of liquidCore.
Random number generator.
Definition: Random.H:57
Primary Breakup Model for pressure swirl atomizers.
virtual void update(const scalar dt, scalar &d, scalar &liquidCore, scalar &tc, const scalar rho, const scalar mu, const scalar sigma, const scalar volFlowRate, const scalar rhoAv, const scalar Urel, const vector &pos, const vector &injectionPos, const scalar pAmbient, const scalar chi, Random &rndGen) const
dimensioned< scalar > mag(const dimensioned< Type > &)
Templated base class for dsmc cloud.
Definition: DSMCCloud.H:69
virtual bool calcChi() const
Flag to indicate if chi needs to be calculated.