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