NumericalStabilization.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include <Eigen/Core>
#include <variant>
#include <vector>
 
#include "BaseLib/Error.h"
#include "NumLib/DOF/GlobalMatrixProviders.h"
 
namespace NumLib
{
struct NoStabilization final
{
};

class IsotropicDiffusionStabilization final
{
public:
    IsotropicDiffusionStabilization(double const cutoff_velocity,
                                    double const tuning_parameter,
                                    std::vector<double>&& element_sizes);
 
    double computeArtificialDiffusion(std::size_t const element_id,
                                      double const velocity_norm) const;
 
private:
    double const cutoff_velocity_;

    double const tuning_parameter_ = 0.5;

    std::vector<double> const element_sizes_;
};

class FullUpwind final
{
public:
    explicit FullUpwind(double const cutoff_velocity)
        : cutoff_velocity_(cutoff_velocity)
    {
    }
 
    double getCutoffVelocity() const { return cutoff_velocity_; }
 
private:
    double const cutoff_velocity_;
};
 
class FluxCorrectedTransport final
{
public:
    explicit FluxCorrectedTransport() {}
 
    void calculateFluxCorrectedTransport(
        const double t, const double dt, std::vector<GlobalVector*> const& x,
        std::vector<GlobalVector*> const& x_prev, int const process_id,
        const MathLib::MatrixSpecifications& matrix_specification,
        GlobalMatrix& M, GlobalMatrix& K, GlobalVector& b);
};
 
using NumericalStabilization =
    std::variant<NoStabilization, IsotropicDiffusionStabilization, FullUpwind,
                 FluxCorrectedTransport>;
}  // namespace NumLib