a7155a7e0a
Patch fields on cyclic patches which have overridden the cyclic constraint using a "patchType cyclic;" setting cannot be decomposed. OpenFOAM does not have processor variants of jumpCyclic, porousBafflePressure, etc... Using these conditions in a decomposed case requires the cyclic to be constrained to a single processor. This change catches this problem in decomposePar and reconstructPar and raises a fatal error, rather than continuing and silently converting these overridden boundary conditions to a standard processorCyclic patch field. Resolves bug report https://bugs.openfoam.org/view.php?id=3916
190 lines
5.2 KiB
C++
190 lines
5.2 KiB
C++
/*---------------------------------------------------------------------------*\
|
|
========= |
|
|
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
|
\\ / O peration | Website: https://openfoam.org
|
|
\\ / A nd | Copyright (C) 2011-2022 OpenFOAM Foundation
|
|
\\/ M anipulation |
|
|
-------------------------------------------------------------------------------
|
|
License
|
|
This file is part of OpenFOAM.
|
|
|
|
OpenFOAM is free software: you can redistribute it and/or modify it
|
|
under the terms of the GNU General Public License as published by
|
|
the Free Software Foundation, either version 3 of the License, or
|
|
(at your option) any later version.
|
|
|
|
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
|
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
|
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
|
for more details.
|
|
|
|
You should have received a copy of the GNU General Public License
|
|
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
|
|
|
Class
|
|
Foam::fvFieldDecomposer
|
|
|
|
Description
|
|
Finite Volume volume and surface field decomposer.
|
|
|
|
SourceFiles
|
|
fvFieldDecomposer.C
|
|
fvFieldDecomposerDecomposeFields.C
|
|
|
|
\*---------------------------------------------------------------------------*/
|
|
|
|
#ifndef fvFieldDecomposer_H
|
|
#define fvFieldDecomposer_H
|
|
|
|
#include "fvMesh.H"
|
|
#include "directFvPatchFieldMapper.H"
|
|
#include "surfaceFields.H"
|
|
|
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
|
|
|
namespace Foam
|
|
{
|
|
|
|
class IOobjectList;
|
|
|
|
/*---------------------------------------------------------------------------*\
|
|
Class fvFieldDecomposer Declaration
|
|
\*---------------------------------------------------------------------------*/
|
|
|
|
class fvFieldDecomposer
|
|
{
|
|
public:
|
|
|
|
// Public Classes
|
|
|
|
//- Patch field decomposer class
|
|
class patchFieldDecomposer
|
|
:
|
|
public labelList,
|
|
public directFvPatchFieldMapper
|
|
{
|
|
public:
|
|
|
|
// Constructors
|
|
|
|
//- Construct given addressing
|
|
patchFieldDecomposer(const labelUList& addressing);
|
|
};
|
|
|
|
|
|
private:
|
|
|
|
// Private Data
|
|
|
|
//- Reference to complete mesh
|
|
const fvMesh& completeMesh_;
|
|
|
|
//- Reference to processor mesh
|
|
const fvMesh& procMesh_;
|
|
|
|
//- Reference to face addressing
|
|
const labelList& faceAddressing_;
|
|
|
|
//- Reference to cell addressing
|
|
const labelList& cellAddressing_;
|
|
|
|
//- Reference to face addressing boundary field
|
|
const surfaceLabelField::Boundary& faceAddressingBf_;
|
|
|
|
//- List of patch field decomposers
|
|
PtrList<patchFieldDecomposer> patchFieldDecomposers_;
|
|
|
|
|
|
// Private Member Functions
|
|
|
|
//- Convert a processor patch to the corresponding complete patch index
|
|
label completePatchID(const label procPatchi) const;
|
|
|
|
//- Map cell values to faces
|
|
template<class Type>
|
|
static tmp<Field<Type>> mapCellToFace
|
|
(
|
|
const labelUList& owner,
|
|
const labelUList& neighbour,
|
|
const Field<Type>& field,
|
|
const labelUList& addressing
|
|
);
|
|
|
|
//- Map face values to faces
|
|
template<class Type>
|
|
static tmp<Field<Type>> mapFaceToFace
|
|
(
|
|
const Field<Type>& field,
|
|
const labelUList& addressing,
|
|
const bool isFlux
|
|
);
|
|
|
|
|
|
public:
|
|
|
|
// Constructors
|
|
|
|
//- Construct from components
|
|
fvFieldDecomposer
|
|
(
|
|
const fvMesh& completeMesh,
|
|
const fvMesh& procMesh,
|
|
const labelList& faceAddressing,
|
|
const labelList& cellAddressing,
|
|
const surfaceLabelField::Boundary& faceAddressingBf
|
|
);
|
|
|
|
//- Disallow default bitwise copy construction
|
|
fvFieldDecomposer(const fvFieldDecomposer&) = delete;
|
|
|
|
|
|
//- Destructor
|
|
~fvFieldDecomposer();
|
|
|
|
|
|
// Member Functions
|
|
|
|
//- Decompose volume field
|
|
template<class Type>
|
|
tmp<GeometricField<Type, fvPatchField, volMesh>>
|
|
decomposeField
|
|
(
|
|
const GeometricField<Type, fvPatchField, volMesh>& field
|
|
) const;
|
|
|
|
//- Decompose surface field
|
|
template<class Type>
|
|
tmp<GeometricField<Type, fvsPatchField, surfaceMesh>>
|
|
decomposeField
|
|
(
|
|
const GeometricField<Type, fvsPatchField, surfaceMesh>& field
|
|
) const;
|
|
|
|
//- Decompose a list of fields
|
|
template<class GeoField>
|
|
void decomposeFields(const PtrList<GeoField>& fields) const;
|
|
|
|
|
|
// Member Operators
|
|
|
|
//- Disallow default bitwise assignment
|
|
void operator=(const fvFieldDecomposer&) = delete;
|
|
};
|
|
|
|
|
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
|
|
|
} // End namespace Foam
|
|
|
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
|
|
|
#ifdef NoRepository
|
|
#include "fvFieldDecomposerDecomposeFields.C"
|
|
#endif
|
|
|
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
|
|
|
#endif
|
|
|
|
// ************************************************************************* //
|