NVDTVD.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) 2011-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::NVDTVD
26 
27 Description
28  Foam::NVDTVD
29 
30 \*---------------------------------------------------------------------------*/
31 
32 #ifndef NVDTVD_H
33 #define NVDTVD_H
34 
35 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
36 
37 namespace Foam
38 {
39 
40 /*---------------------------------------------------------------------------*\
41  Class LimitedScheme Declaration
42 \*---------------------------------------------------------------------------*/
43 
44 class NVDTVD
45 {
46 
47 public:
48 
49  typedef scalar phiType;
50  typedef vector gradPhiType;
51 
52 
53  // Null Constructor
54 
55  NVDTVD()
56  {}
57 
58 
59  // Member Functions
60 
61  scalar phict
62  (
63  const scalar faceFlux,
64  const scalar phiP,
65  const scalar phiN,
66  const vector& gradcP,
67  const vector& gradcN,
68  const vector& d
69  ) const
70  {
71  scalar gradf = phiN - phiP;
72 
73  scalar gradcf;
74 
75  if (faceFlux > 0)
76  {
77  gradcf = d & gradcP;
78  }
79  else
80  {
81  gradcf = d & gradcN;
82  }
83 
84  if (mag(gradf) >= 1000*mag(gradcf))
85  {
86  return 1 - 0.5*1000*sign(gradcf)*sign(gradf);
87  }
88  else
89  {
90  return 1 - 0.5*gradf/gradcf;
91  }
92  }
93 
94 
95  scalar r
96  (
97  const scalar faceFlux,
98  const scalar phiP,
99  const scalar phiN,
100  const vector& gradcP,
101  const vector& gradcN,
102  const vector& d
103  ) const
104  {
105  scalar gradf = phiN - phiP;
106 
107  scalar gradcf;
108 
109  if (faceFlux > 0)
110  {
111  gradcf = d & gradcP;
112  }
113  else
114  {
115  gradcf = d & gradcN;
116  }
117 
118  if (mag(gradcf) >= 1000*mag(gradf))
119  {
120  return 2*1000*sign(gradcf)*sign(gradf) - 1;
121  }
122  else
123  {
124  return 2*(gradcf/gradf) - 1;
125  }
126  }
127 };
128 
129 
130 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
131 
132 } // End namespace Foam
133 
134 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
135 
136 #endif
137 
138 // ************************************************************************* //
Foam::NVDTVD.
Definition: NVDTVD.H:44
vector gradPhiType
Definition: NVDTVD.H:49
scalar phict(const scalar faceFlux, const scalar phiP, const scalar phiN, const vector &gradcP, const vector &gradcN, const vector &d) const
Definition: NVDTVD.H:61
scalar phiType
Definition: NVDTVD.H:48
scalar r(const scalar faceFlux, const scalar phiP, const scalar phiN, const vector &gradcP, const vector &gradcN, const vector &d) const
Definition: NVDTVD.H:95
Namespace for OpenFOAM.
dimensionedScalar sign(const dimensionedScalar &ds)
dimensioned< scalar > mag(const dimensioned< Type > &)