strong_dg_les.hpp

#ifndef PHILIP_DG_STRONG_LES_HPP
#define PHILIP_DG_STRONG_LES_HPP

#include "strong_dg.hpp"
#include "physics/large_eddy_simulation.h"
#include "physics/navier_stokes_model.h"

namespace PHiLiP {

/*  Contains the functions that need to be templated on the number of state variables.
 */
#if PHILIP_DIM==1 // dealii::parallel::distributed::Triangulation<dim> does not work for 1D
template <int dim, int nspecies, int nstate, typename real, typename MeshType = dealii::Triangulation<dim>>
#else
template <int dim, int nspecies, int nstate, typename real, typename MeshType = dealii::parallel::distributed::Triangulation<dim>>
#endif
class DGStrongLES: public DGStrong<dim, nspecies, nstate, real, MeshType>
{
protected:
    using Triangulation = typename DGStrong<dim,nspecies,nstate,real,MeshType>::Triangulation;

public:
    DGStrongLES(
        const Parameters::AllParameters *const parameters_input,
        const unsigned int degree,
        const unsigned int max_degree_input,
        const unsigned int grid_degree_input,
        const std::shared_ptr<Triangulation> triangulation_input);

    ~DGStrongLES();

    std::shared_ptr < Physics::LargeEddySimulationBase<dim, nspecies, nstate, real > > pde_model_les_double;

    virtual void allocate_model_variables() override;

    void update_model_variables() override;

protected:
    void update_cellwise_volume_and_poly_degree();

    virtual void update_cellwise_mean_quantities();

    using DGBase<dim,nspecies,real,MeshType>::pcout; 
    
}; // end of DGStrongLES class

/*  Contains the functions that need to be templated on the number of state variables.
 */
#if PHILIP_DIM==1 // dealii::parallel::distributed::Triangulation<dim> does not work for 1D
template <int dim, int nspecies, int nstate, typename real, typename MeshType = dealii::Triangulation<dim>>
#else
template <int dim, int nspecies, int nstate, typename real, typename MeshType = dealii::parallel::distributed::Triangulation<dim>>
#endif
class DGStrongLES_ShearImproved: public DGStrongLES<dim, nspecies, nstate, real, MeshType>
{
protected:
    using Triangulation = typename DGStrongLES<dim,nspecies,nstate,real,MeshType>::Triangulation;

public:
    DGStrongLES_ShearImproved(
        const Parameters::AllParameters *const parameters_input,
        const unsigned int degree,
        const unsigned int max_degree_input,
        const unsigned int grid_degree_input,
        const std::shared_ptr<Triangulation> triangulation_input);

    ~DGStrongLES_ShearImproved();

    void allocate_model_variables() override;

protected:
    void update_cellwise_mean_quantities() override;

    using DGBase<dim,nspecies,real,MeshType>::pcout; 
    
}; // end of DGStrongLES class

/*  Contains the functions that need to be templated on the number of state variables.
 */
#if PHILIP_DIM==1 // dealii::parallel::distributed::Triangulation<dim> does not work for 1D
template <int dim, int nspecies, int nstate, typename real, typename MeshType = dealii::Triangulation<dim>>
#else
template <int dim, int nspecies, int nstate, typename real, typename MeshType = dealii::parallel::distributed::Triangulation<dim>>
#endif
class DGStrongLES_DynamicSmagorinsky: public DGStrongLES<dim, nspecies, nstate, real, MeshType>
{
protected:
    using Triangulation = typename DGStrongLES<dim,nspecies,nstate,real,MeshType>::Triangulation;

public:
    DGStrongLES_DynamicSmagorinsky(
        const Parameters::AllParameters *const parameters_input,
        const unsigned int degree,
        const unsigned int max_degree_input,
        const unsigned int grid_degree_input,
        const std::shared_ptr<Triangulation> triangulation_input);

    ~DGStrongLES_DynamicSmagorinsky();

    const double dynamic_smagorinsky_model_constant_clipping_limit;

    void allocate_model_variables() override;

protected:
    void update_cellwise_mean_quantities() override;

    using DGBase<dim,nspecies,real,MeshType>::pcout; 
    
}; // end of DGStrongLES_ShearImproved class

/*  Contains the functions that need to be templated on the number of state variables.
 */
#if PHILIP_DIM==1 // dealii::parallel::distributed::Triangulation<dim> does not work for 1D
template <int dim, int nspecies, int nstate, typename real, typename MeshType = dealii::Triangulation<dim>>
#else
template <int dim, int nspecies, int nstate, typename real, typename MeshType = dealii::parallel::distributed::Triangulation<dim>>
#endif
class DGStrong_ChannelFlow: public DGStrong<dim, nspecies, nstate, real, MeshType>
{
protected:
    using Triangulation = typename DGStrong<dim,nspecies,nstate,real,MeshType>::Triangulation;

public:
    DGStrong_ChannelFlow(
        const Parameters::AllParameters *const parameters_input,
        const unsigned int degree,
        const unsigned int max_degree_input,
        const unsigned int grid_degree_input,
        const std::shared_ptr<Triangulation> triangulation_input);

    ~DGStrong_ChannelFlow();

protected:
    const double channel_height; 
    const double half_channel_height; 
    const double channel_friction_velocity_reynolds_number; 
    const int number_of_cells_x_direction; 
    const int number_of_cells_y_direction; 
    const int number_of_cells_z_direction; 
    const double pi_val; 
    const double domain_length_x; 
    const double domain_length_y; 
    const double domain_length_z; 
    const double domain_volume; 

    const double channel_bulk_velocity_reynolds_number;

    const double channel_centerline_velocity_reynolds_number;

    const double total_wall_area; 
public:
    std::shared_ptr < Physics::NavierStokesWithModelSourceTerms<dim, nspecies, nstate, real > > pde_model_navier_stokes_double;

    void allocate_model_variables() override;

    void update_model_variables() override;

protected:
    using DGBase<dim,nspecies,real,MeshType>::pcout; 

    void set_bulk_flow_quantities(); 
    double get_average_wall_shear_stress() const; 
    
}; // end of DGStrong_ChannelFlow class

} // PHiLiP namespace

#endif