thermal_transport_equation.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//
//    FCST: Fuel Cell Simulation Toolbox
//
//    Copyright (C) 2013 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: thermal_transport_equation.h 18-04-2013
//    - Description: Equation class for Thermal transport
//    - Developers: Madhur Bhaiya
//
//---------------------------------------------------------------------------

#ifndef _FCST_FUELCELLSHOP_EQUATION_THERMAL_TRANSPORT_EQUATION_H_
#define _FCST_FUELCELLSHOP_EQUATION_THERMAL_TRANSPORT_EQUATION_H_

#include <utils/fcst_units.h>
#include <equations/equation_base.h>
#include <layers/gas_diffusion_layer.h>
#include <layers/micro_porous_layer.h>
#include <layers/catalyst_layer.h>
#include <layers/membrane_layer.h>

//STD
#include <string>
#include <sstream>


namespace FuelCellShop
{
    namespace Equation
    {


        DeclException2(VariableNotFoundForOhmicHeat,
                       std::string,
                       std::string,
                       << "For "  << arg1 << " ohmic heating source term set as True, \"" << arg2 << "\" is not found as one of the solution variables.");

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



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

            virtual ~ThermalTransportEquation();

            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);




            virtual void print_equation_info() const;

            inline bool get_electron_ohmic_heat_gdl() const
            {   return electron_ohmic_heat_gdl;}

            inline bool get_electron_ohmic_heat_mpl() const
            {   return electron_ohmic_heat_mpl;}

            inline bool get_electron_ohmic_heat_cl() const
            {   return electron_ohmic_heat_cl;}

            inline bool get_proton_ohmic_heat_cl() const
            {   return proton_ohmic_heat_cl;}

            inline bool get_proton_ohmic_heat_ml() const
            {   return proton_ohmic_heat_ml;}


        protected:



            bool electron_ohmic_heat_gdl;

            bool electron_ohmic_heat_mpl;

            bool electron_ohmic_heat_cl;

            bool proton_ohmic_heat_cl;

            bool proton_ohmic_heat_ml;

            bool enthalpy_fickian_transport;

            bool enthalpy_lambda_transport;

            bool flag_thermoosmosis;




            std::map<unsigned int, double> const_heat_flux_map;

            std::map<unsigned int, std::vector<double>> conv_heat_flux_map;




            virtual void make_assemblers_generic_constant_data();

            virtual void make_assemblers_cell_constant_data(const typename FuelCell::ApplicationCore::DoFApplication<dim>::CellInfo& cell_info);

            virtual void make_assemblers_bdry_constant_data(const typename FuelCell::ApplicationCore::DoFApplication<dim>::FaceInfo& bdry_info);

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

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


            template <typename porelayer>
            void gas_enthalpy_transport_assemblers_cell_data(const typename FuelCell::ApplicationCore::DoFApplication<dim>::CellInfo& cell_info,
                                                             porelayer* const layer);

            template <typename polymer>
            void lambda_enthalpy_transport_assemblers_cell_data(const typename FuelCell::ApplicationCore::DoFApplication<dim>::CellInfo& cell_info,
                                                                const polymer* const layer);




            virtual void make_internal_cell_couplings();




            VariableInfo t_rev;

            VariableInfo phi_s;

            VariableInfo phi_m;

            VariableInfo lambda;

            VariableInfo s_liquid_water;
            
            VariableInfo p_liquid_water;

            std::map< std::string, VariableInfo > xi_map;




            std::vector< Tensor<2,dim> > keff_cell;

            std::vector< Tensor<2,dim> > dkeff_dT_cell;

            Tensor<2,dim> sigmaSeff_cell;

            std::vector<double> sigmaMeff_cell;

            std::vector<double> dsigmaMeff_dT_cell;

            std::vector<double> dsigmaMeff_dlambda_cell;


            std::vector< std::vector<double> > phi_T_cell;

            std::vector< std::vector< Tensor<1,dim> > > grad_phi_T_cell;

            std::vector< std::vector< Tensor<1,dim> > > grad_phi_phiS_cell;

            std::vector< std::vector< Tensor<1,dim> > > grad_phi_phiM_cell;

            std::vector< std::vector<double> > phi_lambda_cell;

            std::vector< std::vector< Tensor<1,dim> > > grad_phi_lambda_cell;

            std::vector< std::vector<double> > phi_s_cell;
            
            std::vector< std::vector<double> > phi_p_cell;


            std::vector< Tensor<1,dim> > grad_phiM_cell_old;

            std::vector< Tensor<1,dim> > grad_phiS_cell_old;

            std::vector< Tensor<1,dim> > grad_lambda_cell_old;




            std::map< std::string, std::vector< Tensor<2,dim> > > conc_Deff_dHdT_map;

            std::map< std::string, std::vector< Tensor<2,dim> > > dT_concDeffdHdT_map;

            std::map< std::string, std::vector< Tensor<2,dim> > > ds_concDeffdHdT_map;

            std::map< std::string, std::vector< Tensor<2,dim> > > dp_concDeffdHdT_map;
            std::map< std::string, std::vector< std::vector< Tensor<1,dim> > > > grad_phi_xi_map;

            std::vector<double> electroosmotic_dhdT;

            std::vector<double> delectroosmotic_dhdT_dlambda;

            std::vector<double> delectroosmotic_dhdT_dT;

            std::vector<double> backdiff_dhdT;

            std::vector<double> dbackdiff_dhdT_dlambda;

            std::vector<double> dbackdiff_dhdT_dT;

            std::vector<double> thermoosmotic_dhdT;

            std::vector<double> dthermoosmotic_dhdT_dT;




            std::vector< Tensor<2,dim> > dkeff_dT_bdry;

            std::vector< std::vector<double> > phi_T_bdry;


            bool cell_residual_counter;

            bool bdry_residual_counter;

            unsigned int last_iter_cell;

            unsigned int last_iter_bdry;

        };

    } // Equation

} // FuelCellShop

#endif