FreeFOAM The Cross-Platform CFD Toolkit
fanFvPatchFields.C
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 \*---------------------------------------------------------------------------*/
25 
26 #include "fanFvPatchFields.H"
28 #include <finiteVolume/volFields.H>
30 
31 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
32 
33 namespace Foam
34 {
35  makePatchTypeField(fvPatchScalarField, fanFvPatchScalarField);
36 }
37 
38 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
39 
40 //- Specialisation of the jump-condition for the pressure
41 template<>
43 {
44  if (updated())
45  {
46  return;
47  }
48 
49  jump_ = f_[0];
50 
51  if (f_.size() > 1)
52  {
53  const surfaceScalarField& phi =
54  db().lookupObject<surfaceScalarField>("phi");
55 
56  const fvsPatchField<scalar>& phip =
57  patch().patchField<surfaceScalarField, scalar>(phi);
58 
59  scalarField Un = max
60  (
61  scalarField::subField(phip, size()/2)
62  /scalarField::subField(patch().magSf(), size()/2),
63  scalar(0)
64  );
65 
67  {
68  Un /=
70  (
71  patch().lookupPatchField<volScalarField, scalar>("rho"),
72  size()/2
73  );
74  }
75 
76  for(label i=1; i<f_.size(); i++)
77  {
78  jump_ += f_[i]*pow(Un, i);
79  }
80 
81  jump_ = max(jump_, scalar(0));
82  }
83 
85 }
86 
87 
88 // ************************ vim: set sw=4 sts=4 et: ************************ //