physics_model.h

#ifndef __PHYSICS_MODEL__
#define __PHYSICS_MODEL__

#include "physics.h"
#include "navier_stokes.h"
#include "model.h"

namespace PHiLiP {
namespace Physics {

template <int dim, int nstate, typename real, int nstate_baseline_physics>
class PhysicsModel : public PhysicsBase <dim, nstate, real>
{
public:
    PhysicsModel(
        const Parameters::AllParameters                              *const parameters_input,
        Parameters::AllParameters::PartialDifferentialEquation       baseline_physics_type,
        std::shared_ptr< ModelBase<dim,nstate,real> >                model_input,
        std::shared_ptr< ManufacturedSolutionFunction<dim,real> >    manufactured_solution_function,
        const bool                                                   has_nonzero_diffusion,
        const bool                                                   has_nonzero_physical_source);

    const int n_model_equations;

    std::shared_ptr< PhysicsBase<dim,nstate_baseline_physics,real> > physics_baseline;

    std::shared_ptr< ModelBase<dim,nstate,real> > model;

    std::array<dealii::Tensor<1,dim,real>,nstate> convective_flux (
        const std::array<real,nstate> &conservative_soln) const;

    std::array<dealii::Tensor<1,dim,real>,nstate> dissipative_flux (
        const std::array<real,nstate> &conservative_soln,
        const std::array<dealii::Tensor<1,dim,real>,nstate> &solution_gradient,
        const dealii::types::global_dof_index cell_index) const;

    std::array<real,nstate> physical_source_term (
        const dealii::Point<dim,real> &pos,
        const std::array<real,nstate> &conservative_soln,
        const std::array<dealii::Tensor<1,dim,real>,nstate> &solution_gradient,
        const dealii::types::global_dof_index cell_index) const;

    std::array<real,nstate> source_term (
        const dealii::Point<dim,real> &pos,
        const std::array<real,nstate> &conservative_soln,
        const real current_time,
        const dealii::types::global_dof_index cell_index) const;

    //===========================================================================================
    // All other functions required by PhysicsBase:
    //===========================================================================================
    std::array<dealii::Tensor<1,dim,real>,nstate> convective_numerical_split_flux (
        const std::array<real,nstate> &conservative_soln1,
        const std::array<real,nstate> &conservative_soln2) const;

    std::array<real,nstate> compute_entropy_variables (
                const std::array<real,nstate> &conservative_soln) const;

    std::array<real,nstate> compute_conservative_variables_from_entropy_variables (
                const std::array<real,nstate> &entropy_var) const;

    std::array<real,nstate> convective_eigenvalues (
        const std::array<real,nstate> &/*solution*/,
        const dealii::Tensor<1,dim,real> &/*normal*/) const;

    real max_convective_eigenvalue (const std::array<real,nstate> &soln) const;

    real max_convective_normal_eigenvalue (
        const std::array<real,nstate> &soln,
        const dealii::Tensor<1,dim,real> &normal) const override;

    real max_viscous_eigenvalue (const std::array<real,nstate> &soln) const;

    void boundary_face_values (
        const int /*boundary_type*/,
        const dealii::Point<dim, real> &/*pos*/,
        const dealii::Tensor<1,dim,real> &/*normal*/,
        const std::array<real,nstate> &/*soln_int*/,
        const std::array<dealii::Tensor<1,dim,real>,nstate> &/*soln_grad_int*/,
        std::array<real,nstate> &/*soln_bc*/,
        std::array<dealii::Tensor<1,dim,real>,nstate> &/*soln_grad_bc*/) const;


    dealii::Vector<double> post_compute_derived_quantities_vector (
        const dealii::Vector<double>              &uh,
        const std::vector<dealii::Tensor<1,dim> > &duh,
        const std::vector<dealii::Tensor<2,dim> > &dduh,
        const dealii::Tensor<1,dim>               &normals,
        const dealii::Point<dim>                  &evaluation_points) const;
    

    std::vector<std::string> post_get_names () const;
    

    std::vector<dealii::DataComponentInterpretation::DataComponentInterpretation> post_get_data_component_interpretation () const;
    

    dealii::UpdateFlags post_get_needed_update_flags () const; 

protected:
    const MPI_Comm mpi_communicator; 
    const int mpi_rank; 
    const int n_mpi; 

    dealii::ConditionalOStream pcout;
};

} // Physics namespace
} // PHiLiP namespace

#endif