geometryScheme.hpp

Go to the documentation of this file.

// SPDX-License-Identifier: MIT
// SPDX-FileCopyrightText: 2023 NeoFOAM authors
#pragma once
 
 
#include "NeoFOAM/core/primitives/vector.hpp"
#include "NeoFOAM/core/primitives/scalar.hpp"
#include "NeoFOAM/core/executor/executor.hpp"
#include "NeoFOAM/finiteVolume/cellCentred/fields/surfaceField.hpp"
#include "NeoFOAM/mesh/unstructured.hpp"
 
namespace NeoFOAM::finiteVolume::cellCentred
{
 
class GeometrySchemeFactory
{
 
public:
 
    GeometrySchemeFactory(const UnstructuredMesh& mesh);
 
    virtual ~GeometrySchemeFactory() = default;
 
    virtual void updateWeights(const Executor& exec, SurfaceField<scalar>& weights) = 0;
 
    virtual void updateDeltaCoeffs(const Executor& exec, SurfaceField<scalar>& deltaCoeffs) = 0;
 
    virtual void
    updateNonOrthDeltaCoeffs(const Executor& exec, SurfaceField<scalar>& nonOrthDeltaCoeffs) = 0;
 
    virtual void
    updateNonOrthDeltaCoeffs(const Executor& exec, SurfaceField<Vector>& nonOrthDeltaCoeffs) = 0;
};
 
/* @class GeometryScheme
 * @brief Implements a method to compute deltaCoeffs
 */
class GeometryScheme
{
public:
 
    GeometryScheme(
        const Executor& exec,
        std::unique_ptr<GeometrySchemeFactory> kernel,
        const SurfaceField<scalar>& weights,
        const SurfaceField<scalar>& deltaCoeffs,
        const SurfaceField<scalar>& nonOrthDeltaCoeffs,
        const SurfaceField<Vector>& nonOrthCorrectionVectors
    );
 
    GeometryScheme(
        const Executor& exec,
        const UnstructuredMesh& mesh,
        std::unique_ptr<GeometrySchemeFactory> kernel
    );
 
    GeometryScheme(const UnstructuredMesh& mesh // will lookup the kernel
    );
 
    virtual ~GeometryScheme() = default;
 
    const SurfaceField<scalar>& weights() const;
 
    const SurfaceField<scalar>& deltaCoeffs() const;
 
    const SurfaceField<scalar>& nonOrthDeltaCoeffs() const;
 
    const SurfaceField<Vector>& nonOrthCorrectionVectors() const;
 
    void update();
 
    std::string name() const;
 
    // add selection mechanism via dictionary later
    static const std::shared_ptr<GeometryScheme> readOrCreate(const UnstructuredMesh& mesh);
 
private:
 
    const Executor exec_;
    const UnstructuredMesh& mesh_;
    std::unique_ptr<GeometrySchemeFactory> kernel_;
 
    SurfaceField<scalar> weights_;
    SurfaceField<scalar> deltaCoeffs_;
    SurfaceField<scalar> nonOrthDeltaCoeffs_;
    SurfaceField<Vector> nonOrthCorrectionVectors_;
};
 
} // namespace NeoFOAM