FreeFOAM The Cross-Platform CFD Toolkit
NVDVTVDV.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::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;
51 
52  // Null Constructor
53 
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  // 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 vector& phiP,
95  const vector& phiN,
96  const tensor& gradcP,
97  const tensor& gradcN,
98  const vector& d
99  ) const
100  {
101  vector gradfV = phiN - phiP;
102  scalar gradf = gradfV & gradfV;
103 
104  scalar gradcf;
105 
106  if (faceFlux > 0)
107  {
108  gradcf = gradfV & (d & gradcP);
109  }
110  else
111  {
112  gradcf = gradfV & (d & gradcN);
113  }
114 
115  // Stabilise for division
116  gradf = stabilise(gradf, VSMALL);
117 
118  return 2*(gradcf/gradf) - 1;
119  }
120 };
121 
122 
123 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
124 
125 } // End namespace Foam
126 
127 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
128 
129 #endif
130 
131 // ************************ vim: set sw=4 sts=4 et: ************************ //