fvsPatchField.H

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
    \\  /    A nd           | www.openfoam.com
     \\/     M anipulation  |
-------------------------------------------------------------------------------
    Copyright (C) 2011-2016 OpenFOAM Foundation
    Copyright (C) 2019-2024 OpenCFD Ltd.
-------------------------------------------------------------------------------
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::fvsPatchField
 
Description
    An abstract base class with a fat-interface to all derived classes
    covering all possible ways in which they might be used.
 
    The first level of derivation is to basic patchFields which cover
    zero-gradient, fixed-gradient, fixed-value and mixed conditions.
 
    The next level of derivation covers all the specialised typed with
    specific evaluation procedures, particularly with respect to specific
    fields.
 
SourceFiles
    fvsPatchField.C
    fvsPatchFieldNew.C
 
\*---------------------------------------------------------------------------*/
 
#ifndef Foam_fvsPatchField_H
#define Foam_fvsPatchField_H
 
#include "fvPatch.H"
#include "DimensionedField.H"
#include "fieldTypes.H"
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
namespace Foam
{
 
// Forward Declarations
class dictionary;
class objectRegistry;
class fvPatchFieldMapper;
class surfaceMesh;
 
template<class Type> class fvsPatchField;
template<class Type> class calculatedFvsPatchField;
 
template<class Type>
Ostream& operator<<(Ostream&, const fvsPatchField<Type>&);
 
 
/*---------------------------------------------------------------------------*\
                      Class fvsPatchFieldBase Declaration
\*---------------------------------------------------------------------------*/
 
//- Template invariant parts for fvsPatchField
class fvsPatchFieldBase
{
    // Private Data
 
        //- Reference to patch
        const fvPatch& patch_;
 
        //- Optional patch type
        //  Used to allow specified boundary conditions to be applied
        //  to constraint patches by providing the constraint
        //  patch type as 'patchType'
        word patchType_;
 
 
protected:
 
    // Protected Member Functions
 
        //- Read dictionary entries.
        //  Useful when initially constructed without a dictionary
        virtual void readDict(const dictionary& dict);
 
 
public:
 
    //- Debug switch to disallow the use of generic fvsPatchField
    static int disallowGenericPatchField;
 
    //- Runtime type information
    TypeName("fvsPatchField");
 
 
    // Constructors
 
        //- Construct from patch
        explicit fvsPatchFieldBase(const fvPatch& p);
 
        //- Construct from patch and dictionary (unused)
        fvsPatchFieldBase(const fvPatch& p, const dictionary& dict);
 
        //- Copy construct with new patch
        fvsPatchFieldBase(const fvsPatchFieldBase& rhs, const fvPatch& p);
 
        //- Copy construct
        fvsPatchFieldBase(const fvsPatchFieldBase& rhs);
 
 
    //- Destructor
    virtual ~fvsPatchFieldBase() = default;
 
 
    // Static Member Functions
 
        //- The type name for \c empty patch fields
        static const word& emptyType() noexcept
        {
            return Foam::fieldTypes::emptyType;
        }
 
        //- The type name for \c calculated patch fields
        static const word& calculatedType() noexcept
        {
            return Foam::fieldTypes::calculatedType;
        }
 
        //- The type name for \c zeroValue patch fields
        static const word& zeroValueType() noexcept
        {
            return Foam::fieldTypes::zeroValueType;
        }
 
 
    // Member Functions
 
    // Attributes
 
        //- True if the value of the patch field is altered by assignment
        virtual bool assignable() const
        {
            return true;
        }
 
        //- True if the patch field fixes a value.
        //  Needed to check if a level has to be specified while solving
        //  Poissons equations.
        virtual bool fixesValue() const
        {
            return false;
        }
 
        //- True if the patch field is coupled
        virtual bool coupled() const
        {
            return false;
        }
 
 
    // Access
 
        //- The associated objectRegistry
        const objectRegistry& db() const;
 
        //- Return the patch
        const fvPatch& patch() const noexcept
        {
            return patch_;
        }
 
        //- The optional patch type
        const word& patchType() const noexcept
        {
            return patchType_;
        }
 
        //- The optional patch type
        word& patchType() noexcept
        {
            return patchType_;
        }
 
        //- True if the type does not correspond to the constraint type
        virtual bool constraintOverride() const
        {
            return !patchType_.empty() && patchType_ != type();
        }
 
 
    // Solution
 
        //- True if the boundary condition has already been updated.
        //- This is always true for fvsPatchField
        bool updated() const noexcept
        {
            return true;
        }
 
        //- Set updated state. This is a no-op for fvsPatchField
        void setUpdated(bool state) noexcept
        {}
 
        //- True if the matrix has already been manipulated
        //- Always false for fvsPatchField
        bool manipulatedMatrix() const noexcept
        {
            return false;
        }
 
        //- Set matrix manipulated state.
        //- This is a no-op for fvsPatchField
        void setManipulated(bool state) noexcept
        {}
 
 
    // Check
 
        //- Check that patches are identical
        void checkPatch(const fvsPatchFieldBase& rhs) const;
};
 
 
/*---------------------------------------------------------------------------*\
                        Class fvsPatchField Declaration
\*---------------------------------------------------------------------------*/
 
template<class Type>
class fvsPatchField
:
    public fvsPatchFieldBase,
    public Field<Type>
{
public:
 
    // Public Data Types
 
        //- The patch type for the patch field
        typedef fvPatch Patch;
 
        //- The value_type for the patch field
        typedef Type value_type;
 
        //- The component type for patch field
        typedef typename pTraits<Type>::cmptType cmptType;
 
        //- The internal field type associated with the patch field
        typedef DimensionedField<Type, surfaceMesh> Internal;
 
        //- Type for a \em calculated patch
        typedef calculatedFvsPatchField<Type> Calculated;
 
 
private:
 
    // Private Data
 
        //- Reference to internal field
        const DimensionedField<Type, surfaceMesh>& internalField_;
 
 
protected:
 
    // Protected Member Functions
 
        //- Read the "value" entry into \c *this.
        //  The reading can be optional (default), mandatory etc.
        //  \returns True on success
        bool readValueEntry
        (
            const dictionary& dict,
            IOobjectOption::readOption readOpt = IOobjectOption::LAZY_READ
        );
 
        //- Write \c *this field as a "value" entry
        void writeValueEntry(Ostream& os) const
        {
            Field<Type>::writeEntry("value", os);
        }
 
 
public:
 
    // Declare run-time constructor selection tables
 
        declareRunTimeSelectionTable
        (
            tmp,
            fvsPatchField,
            patch,
            (
                const fvPatch& p,
                const DimensionedField<Type, surfaceMesh>& iF
            ),
            (p, iF)
        );
 
        declareRunTimeSelectionTable
        (
            tmp,
            fvsPatchField,
            patchMapper,
            (
                const fvsPatchField<Type>& ptf,
                const fvPatch& p,
                const DimensionedField<Type, surfaceMesh>& iF,
                const fvPatchFieldMapper& m
            ),
            (dynamic_cast<const fvsPatchFieldType&>(ptf), p, iF, m)
        );
 
        declareRunTimeSelectionTable
        (
            tmp,
            fvsPatchField,
            dictionary,
            (
                const fvPatch& p,
                const DimensionedField<Type, surfaceMesh>& iF,
                const dictionary& dict
            ),
            (p, iF, dict)
        );
 
 
    // Constructors
 
        //- Construct from patch and internal field
        fvsPatchField
        (
            const fvPatch&,
            const DimensionedField<Type, surfaceMesh>&
        );
 
        //- Construct from patch, internal field and value
        fvsPatchField
        (
            const fvPatch&,
            const DimensionedField<Type, surfaceMesh>&,
            const Type& value
        );
 
        //- Construct from patch, internal field and patch field
        fvsPatchField
        (
            const fvPatch&,
            const DimensionedField<Type, surfaceMesh>&,
            const Field<Type>& pfld
        );
 
        //- Construct from patch, internal field and patch field
        fvsPatchField
        (
            const fvPatch&,
            const DimensionedField<Type, surfaceMesh>&,
            Field<Type>&& pfld
        );
 
        //- Construct from patch, internal field and dictionary
        fvsPatchField
        (
            const fvPatch&,
            const DimensionedField<Type, surfaceMesh>&,
            const dictionary& dict,
            IOobjectOption::readOption requireValue = IOobjectOption::MUST_READ
        );
 
        //- Construct, forwarding to readOption variant
        fvsPatchField
        (
            const fvPatch& p,
            const DimensionedField<Type, surfaceMesh>& iF,
            const dictionary& dict,
            const bool needValue
        )
        :
            fvsPatchField
            (
                p, iF, dict,
                (needValue? IOobjectOption::MUST_READ : IOobjectOption::NO_READ)
            )
        {}
 
        //- Construct by mapping the given fvsPatchField onto a new patch
        fvsPatchField
        (
            const fvsPatchField<Type>&,
            const fvPatch&,
            const DimensionedField<Type, surfaceMesh>&,
            const fvPatchFieldMapper&
        );
 
        //- Construct onto new patch with internal field reference
        //- and specified value
        fvsPatchField
        (
            const fvsPatchField<Type>& pfld,
            const fvPatch& p,
            const DimensionedField<Type, surfaceMesh>& iF,
            const Type& value
        );
 
        //- Copy construct with internal field reference
        fvsPatchField
        (
            const fvsPatchField<Type>&,
            const DimensionedField<Type, surfaceMesh>&
        );
 
        //- Copy construct
        fvsPatchField(const fvsPatchField<Type>& pfld)
        :
            fvsPatchField<Type>(pfld, pfld.internalField())
        {}
 
 
        //- Clone patch field with its own internal field reference
        virtual tmp<fvsPatchField<Type>> clone() const
        {
            return tmp<fvsPatchField<Type>>
            (
                new fvsPatchField<Type>(*this, this->internalField_)
            );
        }
 
        //- Clone with an internal field reference
        virtual tmp<fvsPatchField<Type>> clone
        (
            const DimensionedField<Type, surfaceMesh>& iF
        ) const
        {
            return tmp<fvsPatchField<Type>>
            (
                new fvsPatchField<Type>(*this, iF)
            );
        }
 
 
    // Factory Methods
 
        //- Clone a patch field, optionally with internal field reference etc.
        template<class DerivedPatchField, class... Args>
        static tmp<fvsPatchField<Type>> Clone
        (
            const DerivedPatchField& pf,
            Args&&... args
        )
        {
            return tmp<fvsPatchField<Type>>
            (
                new DerivedPatchField(pf, std::forward<Args>(args)...)
            );
        }
 
        //- Return a pointer to a new patchField created on freestore given
        //  patch and internal field
        //  (does not set the patch field values)
        static tmp<fvsPatchField<Type>> New
        (
            const word& patchFieldType,
            const fvPatch&,
            const DimensionedField<Type, surfaceMesh>&
        );
 
        //- Return a pointer to a new patchField created on freestore given
        //  patch and internal field
        //  (does not set the patch field values)
        //  Allows override of constraint type
        static tmp<fvsPatchField<Type>> New
        (
            const word& patchFieldType,
            const word& actualPatchType,
            const fvPatch&,
            const DimensionedField<Type, surfaceMesh>&
        );
 
        //- Return a pointer to a new patchField created on freestore from
        //  a given fvsPatchField mapped onto a new patch
        static tmp<fvsPatchField<Type>> New
        (
            const fvsPatchField<Type>&,
            const fvPatch&,
            const DimensionedField<Type, surfaceMesh>&,
            const fvPatchFieldMapper&
        );
 
        //- Return a pointer to a new patchField created on freestore
        //  from dictionary
        static tmp<fvsPatchField<Type>> New
        (
            const fvPatch&,
            const DimensionedField<Type, surfaceMesh>&,
            const dictionary&
        );
 
        //- Return a pointer to a new calculatedFvsPatchField created on
        //  freestore without setting patchField values
        static tmp<fvsPatchField<Type>> NewCalculatedType
        (
            const fvPatch& p
        );
 
        //- Return a pointer to a new calculatedFvsPatchField created on
        //  freestore without setting patchField values
        template<class AnyType>
        static tmp<fvsPatchField<Type>> NewCalculatedType
        (
            const fvsPatchField<AnyType>& pf
        );
 
 
    //- Destructor
    virtual ~fvsPatchField() = default;
 
 
    // Member Functions
 
    // Access
 
        //- Return const-reference to the dimensioned internal field
        const DimensionedField<Type, surfaceMesh>& internalField()
        const noexcept
        {
            return internalField_;
        }
 
        //- Return const-reference to the internal field values
        const Field<Type>& primitiveField() const noexcept
        {
            return internalField_;
        }
 
 
    // Mapping Functions
 
        //- Map (and resize as needed) from self given a mapping object
        virtual void autoMap
        (
            const fvPatchFieldMapper&
        );
 
        //- Reverse map the given fvsPatchField onto this fvsPatchField
        virtual void rmap
        (
            const fvsPatchField<Type>&,
            const labelList&
        );
 
 
    // Evaluation Functions
 
        //- Initialise the evaluation of the patch field
        virtual void initEvaluate
        (
            const Pstream::commsTypes commsType = Pstream::commsTypes::buffered
        )
        {}
 
        //- Evaluate the patch field, sets Updated to false
        virtual void evaluate
        (
            const Pstream::commsTypes commsType = Pstream::commsTypes::buffered
        )
        {}
 
        //- Initialise the evaluation of the patch field after a local
        //  operation
        virtual void initEvaluateLocal
        (
            const Pstream::commsTypes commsType = Pstream::commsTypes::buffered
        )
        {}
 
        //- Evaluate the patch field after a local operation (e.g. *=)
        virtual void evaluateLocal
        (
            const Pstream::commsTypes commsType = Pstream::commsTypes::buffered
        )
        {}
 
 
    // Other
 
        //- Write the patch "type"
        virtual void write(Ostream& os) const;
 
        //- Check against given patch field
        void check(const fvsPatchField<Type>&) const;
 
 
    // Member Operators
 
        virtual void operator=(const UList<Type>&);
 
        virtual void operator=(const fvsPatchField<Type>&);
        virtual void operator+=(const fvsPatchField<Type>&);
        virtual void operator-=(const fvsPatchField<Type>&);
        virtual void operator*=(const fvsPatchField<scalar>&);
        virtual void operator/=(const fvsPatchField<scalar>&);
 
        virtual void operator+=(const Field<Type>&);
        virtual void operator-=(const Field<Type>&);
 
        virtual void operator*=(const Field<scalar>&);
        virtual void operator/=(const Field<scalar>&);
 
        virtual void operator=(const Type&);
        virtual void operator+=(const Type&);
        virtual void operator-=(const Type&);
        virtual void operator*=(const scalar);
        virtual void operator/=(const scalar);
 
 
        // Force an assignment irrespective of form of patch
 
        virtual void operator==(const fvsPatchField<Type>&);
        virtual void operator==(const Field<Type>&);
        virtual void operator==(const Type&);
 
        // Prevent automatic comparison rewriting (c++20)
        bool operator!=(const fvsPatchField<Type>&) const = delete;
        bool operator!=(const Field<Type>&) const = delete;
        bool operator!=(const Type&) const = delete;
 
 
    // Ostream Operator
 
        friend Ostream& operator<< <Type>(Ostream&, const fvsPatchField<Type>&);
};
 
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
} // End namespace Foam
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#ifdef NoRepository
    #include "fvsPatchField.C"
    #include "fvsPatchFieldNew.C"
    #include "calculatedFvsPatchField.H"
#endif
 
// Runtime selection macros
#include "fvsPatchFieldMacros.H"
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#endif
 
// ************************************************************************* //