incompressible_single_component_NS_equations.h

Go to the documentation of this file.

// ----------------------------------------------------------------------------
//
// FCST: Fuel Cell Simulation Toolbox
//
// Copyright (C) 2006-2015 by Energy Systems Design Laboratory, University of Alberta
//
// This software is distributed under the MIT license
// For more information, see the README file in /doc/LICENSE
//
// - Class: incompressible_single_component_NS_equations.h
// - Description: This class describes steady-state incompressible and isothermal Navier-Stokes
//   fluid transport equations for a single-phase single-component case
// - Developers: Valentin N. Zingan, University of Alberta
//
// ----------------------------------------------------------------------------

#ifndef _FCST_FUELCELLSHOP_EQUATION_INCOMPRESSIBLE_SINGLE_COMPONENT_NS_EQUATIONS_H_
#define _FCST_FUELCELLSHOP_EQUATION_INCOMPRESSIBLE_SINGLE_COMPONENT_NS_EQUATIONS_H_

#include <equations/equation_base.h>
#include <layers/channel.h>
#include <layers/experimental_porous_layer.h>
#include <utils/fcst_units.h>

namespace FuelCellShop
{
    namespace Equation
    {

        template<int dim>
        class IncompressibleSingleComponentNSEquations : public EquationBase<dim>
        {
            public:



                IncompressibleSingleComponentNSEquations(FuelCell::SystemManagement& system_management,boost::shared_ptr< FuelCell::ApplicationCore::ApplicationData > data = 
                            boost::shared_ptr< FuelCell::ApplicationCore::ApplicationData >());

                virtual ~IncompressibleSingleComponentNSEquations();

                virtual void declare_parameters(ParameterHandler& param) const;

                virtual void initialize(ParameterHandler& param);




                virtual void assemble_cell_matrix(FuelCell::ApplicationCore::MatrixVector&                                 cell_matrices,
                                                    const typename FuelCell::ApplicationCore::DoFApplication<dim>::CellInfo& cell_info,
                                                    FuelCellShop::Layer::BaseLayer<dim>* const                               layer);

                virtual void assemble_cell_residual(FuelCell::ApplicationCore::FEVector&                                     cell_residual,
                                                    const typename FuelCell::ApplicationCore::DoFApplication<dim>::CellInfo& cell_info,
                                                    FuelCellShop::Layer::BaseLayer<dim>* const                               layer);

                virtual void assemble_bdry_matrix(FuelCell::ApplicationCore::MatrixVector&                                 bdry_matrices,
                                                    const typename FuelCell::ApplicationCore::DoFApplication<dim>::FaceInfo& bdry_info,
                                                    FuelCellShop::Layer::BaseLayer<dim>* const                               layer);

                virtual void assemble_bdry_residual(FuelCell::ApplicationCore::FEVector&                                     bdry_residual,
                                                    const typename FuelCell::ApplicationCore::DoFApplication<dim>::FaceInfo& bdry_info,
                                                    FuelCellShop::Layer::BaseLayer<dim>* const                               layer);




                const double& get_inlet_outlet_velocity_max() const
                {
                    return inlet_outlet_velocity_max;
                }
                
                const std::string get_inlet_outlet_equation() const
                {
                    return inlet_outlet_equation;
                }
                
                const unsigned int get_inlet_outlet_boundary_ID() const
                {
                    return inlet_outlet_boundary_ID;
                }
                
                const std::string get_press_vel_comp_apply_to() const
                {
                    return press_vel_comp_apply_to;
                }

                virtual void print_equation_info() const;


            protected:



                virtual void make_internal_cell_couplings();

                virtual void make_matrix_block_indices();

                virtual void make_residual_indices();


                // DATA //



                bool inertia_in_channels;

                bool shear_stress_in_channels;

                bool gravity_in_channels;

                bool inertia_in_porous_media;

                bool shear_stress_in_porous_media;

                bool gravity_in_porous_media;

                std::string drag_in_porous_media;

                bool normal_velocity_is_suppressed_weakly;




                Tensor<1,dim> gravity_acceleration;

                SymmetricTensor<2,dim> unit;

                double theta;

                double eta;

                double       inlet_outlet_velocity_max;
                std::string  inlet_outlet_equation;
                unsigned int inlet_outlet_boundary_ID;
                std::string  press_vel_comp_apply_to;


        };

    } // Equation

} // FuelCellShop

#endif