compressible_navier_stokes_explicit_solving_strategy_runge_kutta.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Ruben Zorrilla, Eduard Gómez
//
//
 
#if !defined(KRATOS_COMPRESSIBLE_NAVIER_STOKES_EXPLICIT_SOLVING_STRATEGY_RUNGE_KUTTA)
#define KRATOS_COMPRESSIBLE_NAVIER_STOKES_EXPLICIT_SOLVING_STRATEGY_RUNGE_KUTTA
 
// System includes
 
// External includes
 
// Project includes
#include "includes/define.h"
#include "includes/model_part.h"
#include "factories/factory.h"
#include "solving_strategies/strategies/explicit_solving_strategy_runge_kutta.h"
 
// Application includes
#include "fluid_dynamics_application_variables.h"
#include "compressible_navier_stokes_explicit_solving_strategy.h"
 
namespace Kratos
{
 
 
 
 
 
 
template <class TSparseSpace, class TDenseSpace, class TButcherTableau>
class CompressibleNavierStokesExplicitSolvingStrategyRungeKutta
: public CompressibleNavierStokesExplicitSolvingStrategy<ExplicitSolvingStrategyRungeKutta<TSparseSpace, TDenseSpace, TButcherTableau>>
{
public:
 
    typedef CompressibleNavierStokesExplicitSolvingStrategy<ExplicitSolvingStrategyRungeKutta<TSparseSpace, TDenseSpace, TButcherTableau>> BaseType;
 
    typedef typename BaseType::ExplicitBuilderType ExplicitBuilderType;
 
    typedef typename TDenseSpace::VectorType LocalSystemVectorType;
 
    KRATOS_CLASS_POINTER_DEFINITION(CompressibleNavierStokesExplicitSolvingStrategyRungeKutta);
 
    KRATOS_DEFINE_LOCAL_FLAG(SHOCK_CAPTURING);
 
 
    explicit CompressibleNavierStokesExplicitSolvingStrategyRungeKutta(
        ModelPart &rModelPart,
        Parameters ThisParameters)
        : BaseType(rModelPart)
    {
        KRATOS_TRY
 
        // Validate and assign defaults
        ThisParameters = this->ValidateAndAssignParameters(ThisParameters, this->GetDefaultParameters());
        this->AssignSettings(ThisParameters);
 
        KRATOS_CATCH("")
    }
 
    explicit CompressibleNavierStokesExplicitSolvingStrategyRungeKutta(
        ModelPart &rModelPart,
        typename ExplicitBuilderType::Pointer pExplicitBuilder,
        bool MoveMeshFlag = false,
        int RebuildLevel = 0)
        : BaseType(rModelPart, pExplicitBuilder, MoveMeshFlag, RebuildLevel)
    {
    }
 
    explicit CompressibleNavierStokesExplicitSolvingStrategyRungeKutta(
        ModelPart &rModelPart,
        bool MoveMeshFlag = false,
        int RebuildLevel = 0)
        : BaseType(rModelPart, MoveMeshFlag, RebuildLevel)
    {
    }
 
    CompressibleNavierStokesExplicitSolvingStrategyRungeKutta(const CompressibleNavierStokesExplicitSolvingStrategyRungeKutta &Other) = delete;
 
    virtual ~CompressibleNavierStokesExplicitSolvingStrategyRungeKutta() = default;
 
 
 
 
    Parameters GetDefaultParameters() const override
    {
        KRATOS_TRY
 
        Parameters default_parameters {};
        default_parameters.AddString("explicit_solving_strategy", Name());
 
        // Getting base class default parameters
        const Parameters base_default_parameters = BaseType::GetDefaultParameters();
        default_parameters.RecursivelyAddMissingParameters(base_default_parameters);
        return default_parameters;
 
        KRATOS_CATCH("")
    }
 
    static std::string Name()
    {
        return "compressible_navier_stokes_explicit_solving_strategy_runge_kutta_" + TButcherTableau::Name();
    }
 
    std::string Info() const override
    {
        return "CompressibleNavierStokesExplicitSolvingStrategyRungeKutta";
    }
 
 
protected:
 
 
 
 
 
 
 
    void InitializeRungeKuttaIntermediateSubStep() override
    {
        BaseType::InitializeRungeKuttaIntermediateSubStep();
        this->InitializeEverySubstep();
    }
 
    void InitializeRungeKuttaLastSubStep() override
    {
        BaseType::InitializeRungeKuttaLastSubStep();
        this->InitializeEverySubstep();
    }
 
    void FinalizeRungeKuttaIntermediateSubStep() override
    {
        BaseType::FinalizeRungeKuttaIntermediateSubStep();
        this->FinalizeEverySubstep();
    }
 
    void FinalizeRungeKuttaLastSubStep() override
    {
        BaseType::FinalizeRungeKuttaLastSubStep();
        this->FinalizeEverySubstep();
    }
 
 
 
 
 
 
 
private:
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
}; /* Class CompressibleNavierStokesExplicitSolvingStrategyRungeKutta */
 
 
 
 
template<class TSparseSpace, class TDenseSpace>
using CompressibleNavierStokesExplicitSolvingStrategyForwardEuler = CompressibleNavierStokesExplicitSolvingStrategyRungeKutta<TSparseSpace, TDenseSpace, ButcherTableauForwardEuler>;
 
template<class TSparseSpace, class TDenseSpace>
using CompressibleNavierStokesExplicitSolvingStrategyRungeKutta4 = CompressibleNavierStokesExplicitSolvingStrategyRungeKutta<TSparseSpace, TDenseSpace, ButcherTableauRK4>;
 
template<class TSparseSpace, class TDenseSpace>
using CompressibleNavierStokesExplicitSolvingStrategyRungeKutta3TVD = CompressibleNavierStokesExplicitSolvingStrategyRungeKutta<TSparseSpace, TDenseSpace, ButcherTableauRK3TVD>;
 
} /* namespace Kratos.*/
 
#endif /* KRATOS_COMPRESSIBLE_NAVIER_STOKES_EXPLICIT_SOLVING_STRATEGY_RUNGE_KUTTA  defined */