The open source CFD toolbox
Lun.H
Go to the documentation of this file.
1/*---------------------------------------------------------------------------*\
2 ========= |
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4 \\ / O peration |
5 \\ / A nd | www.openfoam.com
6 \\/ M anipulation |
7-------------------------------------------------------------------------------
8 Copyright (C) 2013-2016 OpenFOAM Foundation
9-------------------------------------------------------------------------------
11 This file is part of OpenFOAM.
12
13 OpenFOAM is free software: you can redistribute it and/or modify it
15 the Free Software Foundation, either version 3 of the License, or
16 (at your option) any later version.
17
18 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
19 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
20 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
21 for more details.
22
23 You should have received a copy of the GNU General Public License
24 along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
25
26Class
27 Foam::ParticleStressModels::Lun
28
29Group
30 grpLagrangianIntermediateMPPICParticleStressSubModels
31
32Description
33 Inter-particle stress model of Lun et al
34
35 The stress value takes the following form:
36 \f[
37 \left( \alpha \rho + \alpha^2 \rho (1 + e) \frac{3}{5}
38 \left( 1 - \left( \frac{\alpha}{\alpha_{pack}} \right)^\frac{1}{3}
39 \right) \right) \frac{1}{3} \sigma^2
40 \f]
41 Here, \f$\alpha\f$ is the volume fraction of the dispersed phase,
42 \f$\rho\f$ is the density of the dispersed phase, \f$e\f$ is a coefficient
43 of restitution, and \f$\sigma\f$ is the RMS velocity fluctuation.
44
45 Reference:
46 \verbatim
47 "Kinetic theories for granular flow: inelastic particles in Couette
48 flow and slightly inelastic particles in a general flowfield"
49 C Lun, S Savage, G Jeffrey, N Chepurniy
50 Journal of Fluid Mechanics
51 Volume 140, Pages 223-256, 1984
52 \endverbatim
53
54SourceFiles
55 Lun.C
56
57\*---------------------------------------------------------------------------*/
58
59#ifndef Lun_H
60#define Lun_H
61
62#include "ParticleStressModel.H"
63
64// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
65
66namespace Foam
67{
68namespace ParticleStressModels
69{
70
71/*---------------------------------------------------------------------------*\
72 Class Lun Declaration
73\*---------------------------------------------------------------------------*/
75class Lun
76:
78{
79 // Private data
80
81 //- Coefficient of restitution
82 scalar e_;
83
84 //- Smallest allowable difference from the packed volume fraction
85 scalar eps_;
86
87
88public:
89
90 //- Runtime type information
91 TypeName("Lun");
92
93
94 //- Constructors
95
96 //- Construct from components
97 Lun(const dictionary& dict);
98
99 //- Construct copy
100 Lun(const Lun& hc);
101
102 //- Clone
103 virtual autoPtr<ParticleStressModel> clone() const
104 {
106 (
107 new Lun(*this)
108 );
109 }
110
111
112 //- Destructor
113 virtual ~Lun();
114
115
116 //- Member Functions
117
118 //- Collision stress
120 (
121 const Field<scalar>& alpha,
122 const Field<scalar>& rho,
123 const Field<scalar>& uRms
124 ) const;
125
126 //- Collision stress derivative w.r.t. the volume fraction
128 (
129 const Field<scalar>& alpha,
130 const Field<scalar>& rho,
131 const Field<scalar>& uRms
132 ) const;
133};
134
135
136// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
137
138} // End namespace ParticleStressModels
139} // End namespace Foam
140
141// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
142
143#endif
144
145// ************************************************************************* //
Generic templated field type.
Definition: Field.H:82
Base class for inter-particle stress models.
Inter-particle stress model of Lun et al.
Definition: Lun.H:77
tmp< Field< scalar > > tau(const Field< scalar > &alpha, const Field< scalar > &rho, const Field< scalar > &uRms) const
Member Functions.
Definition: Lun.C:83
TypeName("Lun")
Runtime type information.
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.
Definition: Lun.C:107
virtual ~Lun()
Destructor.
Definition: Lun.C:75
virtual autoPtr< ParticleStressModel > clone() const
Clone.
Definition: Lun.H:102
Pointer management similar to std::unique_ptr, with some additional methods and type checking.
Definition: autoPtr.H:66
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Definition: dictionary.H:126
A class for managing temporary objects.
Definition: tmp.H:65
Namespace for OpenFOAM.
volScalarField & alpha
dictionary dict
#define TypeName(TypeNameString)
Declare a ClassName() with extra virtual type info.
Definition: typeInfo.H:73