oversetPatchPhiErr.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 | www.openfoam.com
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8  Copyright (C) 2021 OpenCFD Ltd.
9 -------------------------------------------------------------------------------
10 License
11  This file is part of OpenFOAM.
12 
13  OpenFOAM is free software: you can redistribute it and/or modify it
14  under the terms of the GNU General Public License as published by
15  the Free Software Foundation, either version 3 of the License, or
16  (at your option) any later version.
17 
18  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
19  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
20  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
21  for more details.
22 
23  You should have received a copy of the GNU General Public License
24  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
25 
26 \*---------------------------------------------------------------------------*/
27 
28 #include "oversetPatchPhiErr.H"
30 #include "volFields.H"
31 #include "volFieldsFwd.H"
32 
33 // * * * * * * * * * * * * * * * Global Functions * * * * * * * * * * * * * //
34 
36 (
37  const fvScalarMatrix& m,
38  const surfaceScalarField& phi
39 )
40 {
41  const volScalarField::Boundary& bm = m.psi().boundaryField();
42 
43  forAll(bm, patchi)
44  {
45  const auto& fvp = bm[patchi];
46 
48  {
49  const oversetFvPatchField<scalar>& op =
50  refCast<const oversetFvPatchField<scalar>>(fvp);
51 
52  op.fringeFlux(m, phi);
53  }
54  }
55 }
56 
57 
58 // ************************************************************************* //
Forwards and collection of common volume field types.
void fringeFlux(const fvMatrix< Type > &m, const surfaceScalarField &phi) const
Calculate patch flux (helper function). Requires.
#define forAll(list, i)
Loop across all elements in list.
Definition: stdFoam.H:421
const GeometricField< Type, fvPatchField, volMesh > & psi(const label i=0) const
Return psi.
Definition: fvMatrix.H:485
oversetPatchPhiErr
A special matrix type and solver, designed for finite volume solutions of scalar equations. Face addressing is used to make all matrix assembly and solution loops vectorise.
Definition: fvPatchField.H:64
void oversetPatchPhiErr(const fvScalarMatrix &m, const surfaceScalarField &phi)
Prints out sum of m.flux and phi over the fringe faces.
const Type * isA(const U &obj)
Check if dynamic_cast to Type is possible.
Definition: typeInfo.H:88
Boundary condition for use on overset patches. To be run in combination with special dynamicFvMesh ty...