GidaspowViscosity.C
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) 2011-2018 OpenFOAM Foundation
9 -------------------------------------------------------------------------------
10 License
11  This file is part of OpenFOAM.
12 
13  OpenFOAM is free software: you can redistribute it and/or modify it
14  under the terms of the GNU General Public License as published by
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 
26 \*---------------------------------------------------------------------------*/
27 
28 #include "GidaspowViscosity.H"
29 #include "mathematicalConstants.H"
31 
32 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
33 
34 namespace Foam
35 {
36 namespace kineticTheoryModels
37 {
38 namespace viscosityModels
39 {
42 }
43 }
44 }
45 
46 
47 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
48 
50 (
51  const dictionary& dict
52 )
53 :
55 {}
56 
57 
58 // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
59 
61 {}
62 
63 
64 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
65 
68 (
69  const volScalarField& alpha1,
70  const volScalarField& Theta,
71  const volScalarField& g0,
72  const volScalarField& rho1,
73  const volScalarField& da,
74  const dimensionedScalar& e
75 ) const
76 {
77  const scalar sqrtPi = sqrt(constant::mathematical::pi);
78 
79  return da*sqrt(Theta)*
80  (
81  (4.0/5.0)*sqr(alpha1)*g0*(1.0 + e)/sqrtPi
82  + (1.0/15.0)*sqrtPi*g0*(1.0 + e)*sqr(alpha1)
83  + (1.0/6.0)*sqrtPi*alpha1
84  + (10.0/96.0)*sqrtPi/((1.0 + e)*g0)
85  );
86 }
87 
88 
89 // ************************************************************************* //
mathematicalConstants.H
Foam::tmp
A class for managing temporary objects.
Definition: PtrList.H:61
Foam::kineticTheoryModels::viscosityModel
Definition: viscosityModel.H:55
alpha1
const volScalarField & alpha1
Definition: setRegionFluidFields.H:8
Foam::kineticTheoryModels::viscosityModels::Gidaspow::~Gidaspow
virtual ~Gidaspow()
Destructor.
Definition: GidaspowViscosity.C:60
Foam::kineticTheoryModels::viscosityModels::addToRunTimeSelectionTable
addToRunTimeSelectionTable(viscosityModel, Gidaspow, dictionary)
GidaspowViscosity.H
rho1
volScalarField & rho1
Definition: setRegionFluidFields.H:27
dict
dictionary dict
Definition: searchingEngine.H:14
Foam::dictionary
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Definition: dictionary.H:123
addToRunTimeSelectionTable.H
Macros for easy insertion into run-time selection tables.
Foam::dimensioned< scalar >
Foam
Namespace for OpenFOAM.
Definition: atmBoundaryLayer.C:33
Foam::kineticTheoryModels::viscosityModels::Gidaspow::Gidaspow
Gidaspow(const dictionary &dict)
Construct from components.
Definition: GidaspowViscosity.C:50
Foam::sqr
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Definition: dimensionedSymmTensor.C:51
Foam::constant::mathematical::pi
constexpr scalar pi(M_PI)
Foam::sqrt
dimensionedScalar sqrt(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:144
Foam::constant::electromagnetic::e
const dimensionedScalar e
Elementary charge.
Definition: createFields.H:11
Foam::kineticTheoryModels::viscosityModels::defineTypeNameAndDebug
defineTypeNameAndDebug(Gidaspow, 0)
Foam::kineticTheoryModels::viscosityModels::Gidaspow
Definition: GidaspowViscosity.H:53
Foam::kineticTheoryModels::viscosityModels::Gidaspow::nu
tmp< volScalarField > nu(const volScalarField &alpha1, const volScalarField &Theta, const volScalarField &g0, const volScalarField &rho1, const volScalarField &da, const dimensionedScalar &e) const
Definition: GidaspowViscosity.C:68
Foam::GeometricField< scalar, fvPatchField, volMesh >