sorption_source_terms.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: sorption_source_terms.h
//    - Description: This class is used to assemble cell matrix and cell residual
//      corresponding to sorption/desorption of water inside the catalyst layer.
//    - Developers: Marc Secanell, Madhur Bhaiya
//
//---------------------------------------------------------------------------

#ifndef _FCST_FUELCELLSHOP_EQUATION_SORPTION_SOURCE_TERMS_H_
#define _FCST_FUELCELLSHOP_EQUATION_SORPTION_SOURCE_TERMS_H_

#include <equations/equation_base.h>
#include <layers/catalyst_layer.h>


namespace FuelCellShop
{
    namespace Equation
    {  


        DeclException2(VariableNotFoundForSorption,
                       std::string,
                       std::string,
                       << arg1 << " should be one of the solution variables, in order to account for " << arg2 << " .");
        
        
        template<int dim>
        class SorptionSourceTerms : public EquationBase<dim>
        {
        public:
            

            
            SorptionSourceTerms(FuelCell::SystemManagement& system_management,boost::shared_ptr< FuelCell::ApplicationCore::ApplicationData > data = 
            boost::shared_ptr< FuelCell::ApplicationCore::ApplicationData >());
            
            virtual ~SorptionSourceTerms();
            
            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 adjust_internal_cell_couplings(std::vector< couplings_map >& ) const;
            
            virtual void print_equation_info() const;
            
            double get_time_constant() const
            {
                return time_constant;
            }
            
            inline bool get_flag_sorp_heat_cl() const
            {
                return flag_sorp_heat_cl;
            }
            
            
        protected:
            

            
            bool flag_sorp_heat_cl;
            
            

            
            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_cell_variable_data(const typename FuelCell::ApplicationCore::DoFApplication<dim>::CellInfo& cell_info,
                                                            FuelCellShop::Layer::BaseLayer<dim>* const layer);
            
            virtual void assemble_matrix_for_equation(FuelCell::ApplicationCore::MatrixVector& cell_matrices,
                                                      const typename FuelCell::ApplicationCore::DoFApplication<dim>::CellInfo& cell_info,
                                                      const std::string& eq_name,
                                                      const FEValuesBase<dim>& test_fe,
                                                      const std::vector< std::vector<double> >& test_shape_functions,
                                                      const double& sourceterm_factor);
            

            
            VariableInfo x_water;
            
            VariableInfo lambda;
            
            VariableInfo t_rev;
            
            double time_constant;
            
            

                        
            double rho_dry_cell;
            
            double EW_cell;
            
            std::vector<double> lambda_eq_cell;
            
            std::vector<double> dlambdaEq_dxWater_cell;
            
            std::vector<double> dlambdaEq_dT_cell;
            
            std::vector<double> h_sorp_cell;
            
            std::vector<double> dhsorp_dT_cell;
            
            std::vector< std::vector<double> > phi_T_cell;
            
            std::vector< std::vector<double> > phi_lambda_cell;
            
            std::vector< std::vector<double> > phi_xWater_cell;
            
            
            bool cell_residual_counter;
            
             unsigned int last_iter_cell;            
            
        };
             
    } // Equation
             
} // FuelCellShop

#endif