HarrisCrighton.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) 2013-2021 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::ParticleStressModels::HarrisCrighton
26 
27 Description
28  Inter-particle stress model of Harris and Crighton
29 
30  The stress value takes the following form:
31  \f[
32  \frac{P_s \alpha^\beta}{ \mathrm{max} \left( \alpha_{pack} - \alpha ,
33  \epsilon ( 1 - \alpha ) \right) }
34  \f]
35  Here, \f$\alpha\f$ is the volume fraction of the dispersed phase, and the
36  other values are modelling constants. A small value \f$\epsilon\f$ is used
37  to limit the denominator to ensure numerical stability.
38 
39  Reference:
40  \verbatim
41  "Solitons, solitary waves, and voidage disturbances in gas-fluidised
42  beds"
43  S Harris and D Crighton,
44  Journal of Fluid Mechanics
45  Volume 266, Pages 243-276, 1994
46  \endverbatim
47 
48 SourceFiles
49  HarrisCrighton.C
50 
51 \*---------------------------------------------------------------------------*/
52 
53 #ifndef HarrisCrighton_H
54 #define HarrisCrighton_H
55 
56 #include "ParticleStressModel.H"
57 
58 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
59 
60 namespace Foam
61 {
62 namespace ParticleStressModels
63 {
64 
65 /*---------------------------------------------------------------------------*\
66  Class HarrisCrighton Declaration
67 \*---------------------------------------------------------------------------*/
68 
69 class HarrisCrighton
70 :
71  public ParticleStressModel
72 {
73  // Private Data
74 
75  //- Solid pressure coefficient
76  scalar pSolid_;
77 
78  //- Exponent of the volume fraction
79  scalar beta_;
80 
81  //- Smallest allowable difference from the packed volume fraction
82  scalar eps_;
83 
84 
85  // Private Member Functions
86 
87  //- Return the limited denominator of the radial distribution function
88  tmp<Field<scalar>> denominator(const Field<scalar>& alpha) const;
89 
90 
91 public:
92 
93  //- Runtime type information
94  TypeName("HarrisCrighton");
95 
96 
97  //- Constructors
98 
99  //- Construct from components
101 
102  //- Construct copy
103  HarrisCrighton(const HarrisCrighton& hc);
104 
105  //- Clone
106  virtual autoPtr<ParticleStressModel> clone() const
107  {
109  (
110  new HarrisCrighton(*this)
111  );
112  }
113 
114 
115  //- Destructor
116  virtual ~HarrisCrighton();
117 
118 
119  //- Member Functions
120 
121  //- Collision stress
123  (
124  const Field<scalar>& alpha,
125  const Field<scalar>& rho,
126  const Field<scalar>& uRms
127  ) const;
128 
129  //- Collision stress derivative w.r.t. the volume fraction
131  (
132  const Field<scalar>& alpha,
133  const Field<scalar>& rho,
134  const Field<scalar>& uRms
135  ) const;
136 };
137 
138 
139 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
140 
141 } // End namespace ParticleStressModels
142 } // End namespace Foam
143 
144 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
145 
146 #endif
147 
148 // ************************************************************************* //
Base class for inter-particle stress models.
Inter-particle stress model of Harris and Crighton.
tmp< Field< scalar > > tau(const Field< scalar > &alpha, const Field< scalar > &rho, const Field< scalar > &uRms) const
Member Functions.
HarrisCrighton(const dictionary &dict)
Constructors.
virtual autoPtr< ParticleStressModel > clone() const
Clone.
tmp< Field< scalar > > dTaudTheta(const Field< scalar > &alpha, const Field< scalar > &rho, const Field< scalar > &uRms) const
Collision stress derivative w.r.t. the volume fraction.
TypeName("HarrisCrighton")
Runtime type information.
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
Definition: autoPtr.H:51
A list of keyword definitions, which are a keyword followed by any number of values (e....
Definition: dictionary.H:162
A class for managing temporary objects.
Definition: tmp.H:55
volScalarField alpha(IOobject("alpha", runTime.name(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), lambda *max(Ua &U, zeroSensitivity))
Namespace for OpenFOAM.
dictionary dict