reaction_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: reaction_source_terms.h
//    - Description: This class is used to assemble cell matrix and cell residual for reaction
//      source terms in the catalyst layer for various equation classes.
//    - Developers: Madhur Bhaiya and M. Secanell
//    - $Id: reaction_source_terms.h 1459 2013-08-30 22:27:10Z secanell $
//
//---------------------------------------------------------------------------

#ifndef _FCST_FUELCELLSHOP_EQUATION_REACTION_SOURCE_TERMS_H_
#define _FCST_FUELCELLSHOP_EQUATION_REACTION_SOURCE_TERMS_H_

#include "utils/fcst_constants.h"

#include "equations/equation_base.h"

#include "layers/catalyst_layer.h"

#include "reactions/tafel_kinetics.h"
#include "reactions/dual_path_kinetics.h"


namespace FuelCellShop
{
    namespace Equation
    {
        


        DeclException2(VariableNotFoundForKinetics,
                       std::string,
                       std::string,
                       << "For " << arg1 << " kinetics source terms, \"" << arg2 << "\" is not found as one of the solution variables.");
        
        
        
        template<int dim>
        class ReactionSourceTerms : public EquationBase<dim>
        {
        public:
            

            
            ReactionSourceTerms(FuelCell::SystemManagement& system_management);
            
            virtual ~ReactionSourceTerms();
            
            virtual void declare_parameters(ParameterHandler& param) const;
            
            virtual void initialize(ParameterHandler& param);
            
            inline void set_cathode_kinetics(FuelCellShop::Kinetics::BaseKinetics* kinetics)
            {
                cathode_kinetics = kinetics;
            }
            
            inline void set_anode_kinetics(FuelCellShop::Kinetics::BaseKinetics* kinetics)
            {
                anode_kinetics = kinetics;
            }
            
            

            
            virtual void assemble_cell_matrix(AppFrame::MatrixVector& cell_matrices,
                                              const typename AppFrame::DoFApplication<dim>::CellInfo& cell_info,
                                              FuelCellShop::Layer::BaseLayer<dim>* const layer);
            
            virtual void assemble_cell_residual(AppFrame::FEVector& cell_residual,
                                                const typename AppFrame::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;

            
        protected:
            


            FuelCellShop::Kinetics::BaseKinetics* cathode_kinetics;
            
            FuelCellShop::Kinetics::BaseKinetics* anode_kinetics;
            
            

            
            virtual void make_assemblers_generic_constant_data();
            
            virtual void make_assemblers_cell_constant_data(const typename AppFrame::DoFApplication<dim>::CellInfo& cell_info);
            
            virtual void make_assemblers_cell_variable_data(const typename AppFrame::DoFApplication<dim>::CellInfo& cell_info,
                                                            FuelCellShop::Layer::BaseLayer<dim>* const layer);
            
            virtual void assemble_matrix_for_equation(AppFrame::MatrixVector& cell_matrices,
                                                      const typename AppFrame::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_oxygen;
            
            VariableInfo x_water;
            
            VariableInfo phi_s;
            
            VariableInfo phi_m;
            
            VariableInfo t_rev;
            
            double F;
            
            

            
            std::vector< std::vector<double> > phi_T_cell;
            
            std::vector< std::vector<double> > phi_xOxygen_cell;
            
            std::vector< std::vector<double> > phi_phiS_cell;
            
            std::vector< std::vector<double> > phi_phiM_cell;
            
            std::vector< std::vector<double> > phi_xWater_cell;
            
            std::vector<double> current_cell;
            
            std::vector<double> dcurrent_dphiS_cell;
            
            std::vector<double> dcurrent_dphiM_cell;
            
            std::vector<double> dcurrent_dT_cell;
            
            std::vector<double> dcurrent_dxOxygen_cell;
            
            std::vector<double> dcurrent_dxWater_cell;
            
            double active_area_cell;
            
            double factor_protontranseq_cell;
            
            double factor_electrontranseq_cell;
            
            double factor_oxygentranseq_cell;
            
            double factor_watertranseq_cell;
            
            

            
            bool cell_matrix_counter;
            
            bool cell_residual_counter;
            
            
             unsigned int last_iter_cell;            
            
        };
             
    } // Equation
             
} // FuelCellShop

#endif