NonSphereDragForce.H
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-2019 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 Class
25  Foam::NonSphereDragForce
26 
27 Description
28  Drag model for non-spherical particles.
29 
30  The drag coefficient is given by:
31  \f[
32  Cd = \frac{24}{Re} (1 + a Re^{b}) + Re \frac{c}{Re + d}
33  \f]
34 
35  Where \f$a\f$, \f$b\f$, \f$c\f$, and \f$d\f$ are coefficients, calculated
36  as functions of \f$phi\f$. \f$phi\f$ is the ratio of the surface area of a
37  sphere with the same volume as the particle to the actual surface area of
38  the particle, and must be between 0 and 1.
39 
40  Reference:
41  \verbatim
42  Haider, A., & Levenspiel, O. (1989).
43  Drag coefficient and terminal velocity of spherical and nonspherical
44  particles.
45  Powder technology, 58(1), 63-70.
46  \endverbatim
47 
48 \*---------------------------------------------------------------------------*/
49 
50 #ifndef NonSphereDragForce_H
51 #define NonSphereDragForce_H
52 
53 #include "ParticleForce.H"
54 
55 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
56 
57 namespace Foam
58 {
59 
60 /*---------------------------------------------------------------------------*\
61  Class NonSphereDragForce Declaration
62 \*---------------------------------------------------------------------------*/
63 
64 template<class CloudType>
66 :
67  public ParticleForce<CloudType>
68 {
69 protected:
70 
71  // Protected Data
72 
73  //- The ratio of the surface area of a sphere with the same volume as
74  // the particle to the actual surface area of the particle. Between 0
75  // and 1.
76  const scalar phi_;
77 
78  //- Coefficient
79  const scalar a_;
80 
81  //- Coefficient
82  const scalar b_;
83 
84  //- Coefficient
85  const scalar c_;
86 
87  //- Coefficient
88  const scalar d_;
89 
90 
91 public:
92 
93  //- Runtime type information
94  TypeName("nonSphereDrag");
95 
96 
97  // Constructors
98 
99  //- Construct from mesh
101  (
102  CloudType& owner,
103  const fvMesh& mesh,
104  const dictionary& dict
105  );
106 
107  //- Construct copy
109 
110  //- Construct and return a clone
111  virtual autoPtr<ParticleForce<CloudType>> clone() const
112  {
114  (
116  );
117  }
118 
119 
120  //- Destructor
121  virtual ~NonSphereDragForce();
122 
123 
124  // Member Functions
125 
126  // Evaluation
127 
128  //- Calculate the non-coupled force
129  virtual forceSuSp calcCoupled
130  (
131  const typename CloudType::parcelType& p,
132  const typename CloudType::parcelType::trackingData& td,
133  const scalar dt,
134  const scalar mass,
135  const scalar Re,
136  const scalar muc
137  ) const;
138 };
139 
140 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
141 
142 } // End namespace Foam
143 
144 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
145 
146 #ifdef NoRepository
147  #include "NonSphereDragForce.C"
148 #endif
149 
150 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
151 
152 #endif
153 
154 // ************************************************************************* //
dictionary dict
const scalar a_
Coefficient.
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:158
const fvMesh & mesh() const
Return the mesh database.
Abstract base class for particle forces.
Definition: ParticleForce.H:53
const CloudType & owner() const
Return const access to the cloud owner.
const scalar phi_
The ratio of the surface area of a sphere with the same volume as.
virtual autoPtr< ParticleForce< CloudType > > clone() const
Construct and return a clone.
Helper container for force Su and Sp terms.
Definition: forceSuSp.H:61
TypeName("nonSphereDrag")
Runtime type information.
ParcelType parcelType
Type of parcel the cloud was instantiated for.
Definition: DSMCCloud.H:215
const scalar c_
Coefficient.
virtual ~NonSphereDragForce()
Destructor.
Mesh data needed to do the Finite Volume discretisation.
Definition: fvMesh.H:78
const scalar b_
Coefficient.
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
Definition: PtrList.H:52
virtual forceSuSp calcCoupled(const typename CloudType::parcelType &p, const typename CloudType::parcelType::trackingData &td, const scalar dt, const scalar mass, const scalar Re, const scalar muc) const
Calculate the non-coupled force.
const scalar d_
Coefficient.
volScalarField & p
Drag model for non-spherical particles.
Templated base class for dsmc cloud.
Definition: DSMCCloud.H:69
scalarField Re(const UList< complex > &cf)
Definition: complexFields.C:97
NonSphereDragForce(CloudType &owner, const fvMesh &mesh, const dictionary &dict)
Construct from mesh.
Namespace for OpenFOAM.