turbulent_channel_flow_skin_friction_check.h

#ifndef __TURBULENT_CHANNEL_FLOW_SKIN_FRICTION_CHECK__
#define __TURBULENT_CHANNEL_FLOW_SKIN_FRICTION_CHECK__

#include "tests.h"
#include "parameters/parameters_flow_solver.h"
#include "physics/navier_stokes.h"

namespace PHiLiP {
namespace Tests {

template <int dim, int nspecies, int nstate>
class TurbulentChannelFlowSkinFrictionCheck: public TestsBase
{
private:
    using XVelocityInitialConditionEnum = Parameters::FlowSolverParam::XVelocityInitialConditionType;
public:
    TurbulentChannelFlowSkinFrictionCheck(
        const Parameters::AllParameters *const parameters_input,
        const dealii::ParameterHandler &parameter_handler_input);

    ~TurbulentChannelFlowSkinFrictionCheck() {};

    const dealii::ParameterHandler &parameter_handler;
    
    const double half_channel_height;

    const XVelocityInitialConditionEnum xvelocity_initial_condition_type;

    const double y_top_wall; 
    const double y_bottom_wall; 
    const double normal_vector_top_wall; 
    const double normal_vector_bottom_wall; 

    std::shared_ptr< Physics::NavierStokes_ChannelFlowConstantSourceTerm_WallModel<dim,nspecies,dim+2,double> > navier_stokes_channel_flow_constant_source_term_wall_model_physics;
    bool check_wall_model; 

    int run_test () const override;
private:
    double get_x_velocity(const double y) const;
    double get_x_velocity_gradient(const double y) const;
    double get_wall_shear_stress() const;
    double get_bulk_velocity() const;
    double get_skin_friction_coefficient() const;
    double get_integral_of_x_velocity(const double y_plus) const;
    double get_wall_shear_stress_from_friction_reynolds_number() const;
    double get_wall_shear_stress_from_wall_model() const;
};

} // End of Tests namespace
} // End of PHiLiP namespace

#endif