Giesekus.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) 2019-2020 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::laminarModels::Giesekus
26 
27 Description
28  Giesekus model for viscoelasticity using the upper-convected time
29  derivative of the stress tensor with support for multiple modes.
30 
31  Reference:
32  \verbatim
33  Giesekus, H., 1982.
34  A simple constitutive equation for polymer fluids based on the
35  concept of deformation-dependent tensional mobility.
36  J. Non-Newton. Fluid. 11, 69–109.
37  \endverbatim
38 
39 See also
40  Foam::laminarModels::Maxwell
41 
42 SourceFiles
43  Giesekus.C
44 
45 \*---------------------------------------------------------------------------*/
46 
47 #ifndef Giesekus_H
48 #define Giesekus_H
49 
50 #include "Maxwell.H"
51 
52 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
53 
54 namespace Foam
55 {
56 namespace laminarModels
57 {
58 
59 /*---------------------------------------------------------------------------*\
60  Class Giesekus Declaration
61 \*---------------------------------------------------------------------------*/
62 
63 template<class BasicMomentumTransportModel>
64 class Giesekus
65 :
66  public Maxwell<BasicMomentumTransportModel>
67 {
68 protected:
69 
70  // Protected data
71 
72  // Mode coefficients
73 
75 
76 
77  // Protected Member Functions
78 
80  (
81  const label modei,
83  ) const;
84 
85 
86 public:
87 
88  typedef typename BasicMomentumTransportModel::alphaField alphaField;
89  typedef typename BasicMomentumTransportModel::rhoField rhoField;
90  typedef typename BasicMomentumTransportModel::transportModel transportModel;
91 
92 
93  //- Runtime type information
94  TypeName("Giesekus");
95 
96 
97  // Constructors
98 
99  //- Construct from components
100  Giesekus
101  (
102  const alphaField& alpha,
103  const rhoField& rho,
104  const volVectorField& U,
105  const surfaceScalarField& alphaRhoPhi,
106  const surfaceScalarField& phi,
107  const transportModel& transport,
108  const word& type = typeName
109  );
110 
111  //- Disallow default bitwise copy construction
112  Giesekus(const Giesekus&) = delete;
113 
114 
115 
116  //- Destructor
117  virtual ~Giesekus()
118  {}
119 
120 
121  // Member Functions
122 
123  //- Re-read model coefficients if they have changed
124  virtual bool read();
125 
126 
127  // Member Operators
128 
129  //- Disallow default bitwise assignment
130  void operator=(const Giesekus&) = delete;
131 };
132 
133 
134 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
135 
136 } // End namespace laminarModels
137 } // End namespace Foam
138 
139 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
140 
141 #ifdef NoRepository
142  #include "Giesekus.C"
143 #endif
144 
145 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
146 
147 #endif
148 
149 // ************************************************************************* //
PtrList< dimensionedScalar > alphaGs_
Definition: Giesekus.H:73
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Definition: label.H:59
Generalised Maxwell model for viscoelasticity using the upper-convected time derivative of the stress...
Definition: Maxwell.H:76
virtual bool read()
Re-read model coefficients if they have changed.
Definition: Giesekus.C:92
volScalarField alpha(IOobject("alpha", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), lambda *max(Ua &U, zeroSensitivity))
BasicMomentumTransportModel::transportModel transportModel
Definition: Giesekus.H:89
virtual tmp< fvSymmTensorMatrix > sigmaSource(const label modei, volSymmTensorField &sigma) const
Definition: Giesekus.C:41
Giesekus(const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport, const word &type=typeName)
Construct from components.
Definition: Giesekus.C:59
BasicMomentumTransportModel::alphaField alphaField
Definition: Giesekus.H:87
virtual ~Giesekus()
Destructor.
Definition: Giesekus.H:116
A class for handling words, derived from string.
Definition: word.H:59
Giesekus model for viscoelasticity using the upper-convected time derivative of the stress tensor wit...
Definition: Giesekus.H:63
phi
Definition: correctPhi.H:3
virtual tmp< volSymmTensorField > sigma() const
Return the stress tensor [m^2/s^2].
Definition: Maxwell.C:269
void operator=(const Giesekus &)=delete
Disallow default bitwise assignment.
U
Definition: pEqn.H:72
A templated 1D list of pointers to objects of type <T>, where the size of the array is known and used...
Definition: List.H:70
TypeName("Giesekus")
Runtime type information.
fileType type(const fileName &, const bool checkVariants=true, const bool followLink=true)
Return the file type: directory or file.
Definition: POSIX.C:488
A class for managing temporary objects.
Definition: PtrList.H:53
BasicMomentumTransportModel::rhoField rhoField
Definition: Giesekus.H:88
Namespace for OpenFOAM.