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-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::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 
91 
92  //- Runtime type information
93  TypeName("Giesekus");
94 
95 
96  // Constructors
97 
98  //- Construct from components
99  Giesekus
100  (
101  const alphaField& alpha,
102  const rhoField& rho,
103  const volVectorField& U,
104  const surfaceScalarField& alphaRhoPhi,
105  const surfaceScalarField& phi,
106  const viscosity& viscosity,
107  const word& type = typeName
108  );
109 
110  //- Disallow default bitwise copy construction
111  Giesekus(const Giesekus&) = delete;
112 
113 
114 
115  //- Destructor
116  virtual ~Giesekus()
117  {}
118 
119 
120  // Member Functions
121 
122  //- Re-read model coefficients if they have changed
123  virtual bool read();
124 
125 
126  // Member Operators
127 
128  //- Disallow default bitwise assignment
129  void operator=(const Giesekus&) = delete;
130 };
131 
132 
133 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
134 
135 } // End namespace laminarModels
136 } // End namespace Foam
137 
138 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
139 
140 #ifdef NoRepository
141  #include "Giesekus.C"
142 #endif
143 
144 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
145 
146 #endif
147 
148 // ************************************************************************* //
PtrList< dimensionedScalar > alphaGs_
Definition: Giesekus.H:73
U
Definition: pEqn.H:72
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))
virtual tmp< fvSymmTensorMatrix > sigmaSource(const label modei, volSymmTensorField &sigma) const
Definition: Giesekus.C:41
BasicMomentumTransportModel::alphaField alphaField
Definition: Giesekus.H:87
Giesekus(const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const viscosity &viscosity, const word &type=typeName)
Construct from components.
Definition: Giesekus.C:59
virtual ~Giesekus()
Destructor.
Definition: Giesekus.H:115
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
Abstract base class for all fluid physical properties.
Definition: viscosity.H:49
void operator=(const Giesekus &)=delete
Disallow default bitwise assignment.
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.