NumericalStabilization.h

Go to the documentation of this file.

 
#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()
    {
#ifdef USE_PETSC
        OGS_FATAL(
            "FluxCorrectedTransport scheme is not implemented to work with MPI "
            "parallelization.");
#endif
    }
 
    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