All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
limitedCubicV.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-2022 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::limitedCubicVLimiter
26 
27 Description
28  Class with limiter function which returns the limiter for the
29  limitedCubicV interpolation scheme based on r obtained from the LimiterFunc
30  class.
31 
32  Used in conjunction with the template class LimitedScheme.
33 
34 SourceFiles
35  limitedCubicV.C
36 
37 \*---------------------------------------------------------------------------*/
38 
39 #ifndef limitedCubicV_H
40 #define limitedCubicV_H
41 
42 #include "vector.H"
43 
44 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
45 
46 namespace Foam
47 {
48 
49 /*---------------------------------------------------------------------------*\
50  Class limitedCubicVLimiter Declaration
51 \*---------------------------------------------------------------------------*/
52 
53 template<class LimiterFunc>
55 :
56  public LimiterFunc
57 {
58  scalar k_;
59  scalar twoByk_;
60 
61 public:
62 
64  :
65  k_(readScalar(is))
66  {
67  if (k_ < 0 || k_ > 1)
68  {
70  << "coefficient = " << k_
71  << " should be >= 0 and <= 1"
72  << exit(FatalIOError);
73  }
74 
75  // Avoid the /0 when k_ = 0
76  twoByk_ = 2.0/max(k_, small);
77  }
78 
79  scalar limiter
80  (
81  const scalar cdWeight,
82  const scalar faceFlux,
83  const typename LimiterFunc::phiType& phiP,
84  const typename LimiterFunc::phiType& phiN,
85  const typename LimiterFunc::gradPhiType& gradcP,
86  const typename LimiterFunc::gradPhiType& gradcN,
87  const vector& d
88  ) const
89  {
90  scalar twor = twoByk_*LimiterFunc::r
91  (
92  faceFlux, phiP, phiN, gradcP, gradcN, d
93  );
94 
95  vector fV = cdWeight*phiP + (1.0 - cdWeight)*phiN;
96 
97  scalar fVphiP = fV & phiP;
98  scalar fVphiN = fV & phiN;
99 
100  scalar fVphiU;
101 
102  if (faceFlux > 0)
103  {
104  fVphiU = fVphiP;
105  }
106  else
107  {
108  fVphiU = fVphiN;
109  }
110 
111  // Calculate the face value using cubic interpolation
112  scalar fVphif =
113  cdWeight*(fVphiP - 0.25*(fV & (d & gradcN)))
114  + (1 - cdWeight)*(fVphiN + 0.25*(fV & (d & gradcP)));
115 
116  scalar fVphiCD = cdWeight*fVphiP + (1 - cdWeight)*fVphiN;
117 
118  // Calculate the effective limiter for the cubic interpolation
119  scalar cubicLimiter =
120  (fVphif - fVphiU)/stabilise(fVphiCD - fVphiU, small);
121 
122  // Limit the limiter to obey the TVD constraint
123  return max(min(min(twor, cubicLimiter), 2), 0);
124  }
125 };
126 
127 
128 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
129 
130 } // End namespace Foam
131 
132 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
133 
134 #endif
135 
136 // ************************************************************************* //
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)
errorManipArg< error, int > exit(error &err, const int errNo=1)
Definition: errorManip.H:124
An Istream is an abstract base class for all input systems (streams, files, token lists etc)...
Definition: Istream.H:57
limitedCubicVLimiter(Istream &is)
Definition: limitedCubicV.H:62
scalar limiter(const scalar cdWeight, const scalar faceFlux, const typename LimiterFunc::phiType &phiP, const typename LimiterFunc::phiType &phiN, const typename LimiterFunc::gradPhiType &gradcP, const typename LimiterFunc::gradPhiType &gradcN, const vector &d) const
Definition: limitedCubicV.H:79
Class with limiter function which returns the limiter for the limitedCubicV interpolation scheme base...
Definition: limitedCubicV.H:53
layerAndWeight min(const layerAndWeight &a, const layerAndWeight &b)
bool readScalar(const char *buf, doubleScalar &s)
Read whole of buf as a scalar. Return true if successful.
Definition: doubleScalar.H:75
tmp< DimensionedField< scalar, GeoMesh > > stabilise(const DimensionedField< scalar, GeoMesh > &dsf, const dimensioned< scalar > &ds)
#define FatalIOErrorInFunction(ios)
Report an error message using Foam::FatalIOError.
Definition: error.H:318
Namespace for OpenFOAM.
IOerror FatalIOError