PlessisMasliyahDragForce.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) 2013-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 
27 #include "volFields.H"
28 
29 // * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
30 
31 template<class CloudType>
33 (
34  const scalar Re
35 ) const
36 {
37  if (Re > 1000.0)
38  {
39  return 0.44*Re;
40  }
41  else
42  {
43  return 24.0*(1.0 + 0.15*pow(Re, 0.687));
44  }
45 }
46 
47 
48 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
49 
50 template<class CloudType>
52 (
53  CloudType& owner,
54  const fvMesh& mesh,
55  const dictionary& dict
56 )
57 :
58  ParticleForce<CloudType>(owner, mesh, dict, typeName, true),
59  alphac_
60  (
61  this->mesh().template lookupObject<volScalarField>
62  (
63  this->coeffs().lookup("alphac")
64  )
65  )
66 {}
67 
68 
69 template<class CloudType>
71 (
73 )
74 :
76  alphac_
77  (
78  this->mesh().template lookupObject<volScalarField>
79  (
80  this->coeffs().lookup("alphac")
81  )
82  )
83 {}
84 
85 
86 // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
87 
88 template<class CloudType>
90 {}
91 
92 
93 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
94 
95 template<class CloudType>
97 (
98  const typename CloudType::parcelType& p,
99  const typename CloudType::parcelType::trackingData& td,
100  const scalar dt,
101  const scalar mass,
102  const scalar Re,
103  const scalar muc
104 ) const
105 {
106  scalar alphac(alphac_[p.cell()]);
107 
108  scalar cbrtAlphap(pow(1.0 - alphac, 1.0/3.0));
109 
110  scalar A =
111  26.8*pow3(alphac)
112  /(
113  sqr(cbrtAlphap)
114  *(1.0 - cbrtAlphap)
115  *sqr(1.0 - sqr(cbrtAlphap))
116  + small
117  );
118 
119  scalar B =
120  sqr(alphac)
121  /sqr(1.0 - sqr(cbrtAlphap));
122 
123  return forceSuSp
124  (
125  Zero,
126  (mass/p.rho())
127  *(A*(1.0 - alphac)/alphac + B*Re)*muc/(alphac*sqr(p.d()))
128  );
129 }
130 
131 
132 // ************************************************************************* //
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)
Abstract base class for particle forces.
Definition: ParticleForce.H:53
Helper container for force Su and Sp terms.
Definition: forceSuSp.H:61
dynamicFvMesh & mesh
static const zero Zero
Definition: zero.H:97
PlessisMasliyahDragForce(CloudType &owner, const fvMesh &mesh, const dictionary &dict)
Construct from mesh.
virtual ~PlessisMasliyahDragForce()
Destructor.
ParcelType parcelType
Type of parcel the cloud was instantiated for.
Definition: DSMCCloud.H:218
PlessisMasliyahDragForce drag model for solid spheres.
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 coupled force.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar pow3(const dimensionedScalar &ds)
Mesh data needed to do the Finite Volume discretisation.
Definition: fvMesh.H:78
Templated base class for dsmc cloud.
Definition: DSMCCloud.H:69
scalarField Re(const UList< complex > &cf)
Definition: complexFields.C:97