parameters_flow_solver.h

#ifndef __PARAMETERS_FLOW_SOLVER_H__
#define __PARAMETERS_FLOW_SOLVER_H__

#include <deal.II/base/parameter_handler.h>
#include "parameters/parameters.h"

namespace PHiLiP {

namespace Parameters {

class FlowSolverParam
{
public:

    enum FlowCaseType{
        taylor_green_vortex,
        decaying_homogeneous_isotropic_turbulence,
        burgers_viscous_snapshot,
        naca0012,
        burgers_rewienski_snapshot,
        burgers_inviscid,
        convection_diffusion,
        advection,
        periodic_1D_unsteady,
        gaussian_bump,
        isentropic_vortex,
        kelvin_helmholtz_instability,
        non_periodic_cube_flow,
        sod_shock_tube,
        low_density_2d,
        leblanc_shock_tube,
        shu_osher_problem,
        advection_limiter,
        burgers_limiter
        };
    FlowCaseType flow_case_type; 

    unsigned int poly_degree; 
    unsigned int max_poly_degree_for_adaptation; 
    

    double final_time; 
    double constant_time_step; 
    double courant_friedrichs_lewy_number; 

    std::string unsteady_data_table_filename;

    bool steady_state; 
    bool steady_state_polynomial_ramping; 

    bool error_adaptive_time_step; 
    bool adaptive_time_step; 

    std::string sensitivity_table_filename;

    std::string input_mesh_filename;
    bool use_gmsh_mesh; 
    bool mesh_reader_verbose_output;

    bool use_periodic_BC_in_x; 
    bool use_periodic_BC_in_y; 
    bool use_periodic_BC_in_z; 

    int x_periodic_id_face_1; 
    int x_periodic_id_face_2; 

    int y_periodic_id_face_1; 
    int y_periodic_id_face_2; 

    int z_periodic_id_face_1; 
    int z_periodic_id_face_2; 

    bool restart_computation_from_file; 
    bool output_restart_files; 
    std::string restart_files_directory_name; 
    unsigned int restart_file_index; 
    int output_restart_files_every_x_steps; 
    double output_restart_files_every_dt_time_intervals; 

    unsigned int grid_degree; 
    double grid_left_bound; 
    double grid_right_bound; 
    unsigned int number_of_grid_elements_per_dimension; 
    int number_of_mesh_refinements; 
    double channel_height; 
    double channel_length; 
    double bump_height; 
    int number_of_subdivisions_in_x_direction; 
    int number_of_subdivisions_in_y_direction; 
    int number_of_subdivisions_in_z_direction; 

    double expected_kinetic_energy_at_final_time;

    double expected_theoretical_dissipation_rate_at_final_time;

    enum DensityInitialConditionType{
        uniform,
        isothermal,
        };
    DensityInitialConditionType density_initial_condition_type;
    bool do_calculate_numerical_entropy;

    double atwood_number;

    enum ApplyInitialConditionMethod{
        interpolate_initial_condition_function,
        project_initial_condition_function,
        read_values_from_file_and_project
        };
    ApplyInitialConditionMethod apply_initial_condition_method; 

    std::string input_flow_setup_filename_prefix;

    bool output_velocity_field_at_fixed_times; 
    std::string output_velocity_field_times_string; 
    unsigned int number_of_times_to_output_velocity_field; 
    bool output_vorticity_magnitude_field_in_addition_to_velocity; 
    std::string output_flow_field_files_directory_name; 

    bool end_exactly_at_final_time; 

    static void declare_parameters (dealii::ParameterHandler &prm);

    void parse_parameters (dealii::ParameterHandler &prm);
};

} // Parameters namespace

} // PHiLiP namespace

#endif