FreeFOAM The Cross-Platform CFD Toolkit
NVDTVD.H
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 1991-2010 OpenCFD Ltd.
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;
51 
52 
53  // Null Constructor
54 
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  // Stabilise for division
85  gradcf = stabilise(gradcf, VSMALL);
86 
87  return 1 - 0.5*gradf/gradcf;
88  }
89 
90 
91  scalar r
92  (
93  const scalar faceFlux,
94  const scalar phiP,
95  const scalar phiN,
96  const vector& gradcP,
97  const vector& gradcN,
98  const vector& d
99  ) const
100  {
101  scalar gradf = phiN - phiP;
102 
103  scalar gradcf;
104 
105  if (faceFlux > 0)
106  {
107  gradcf = d & gradcP;
108  }
109  else
110  {
111  gradcf = d & gradcN;
112  }
113 
114  // Stabilise for division
115  gradf = stabilise(gradf, VSMALL);
116 
117  return 2*(gradcf/gradf) - 1;
118  }
119 };
120 
121 
122 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
123 
124 } // End namespace Foam
125 
126 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
127 
128 #endif
129 
130 // ************************ vim: set sw=4 sts=4 et: ************************ //