reaction_source_terms.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: reaction_source_terms.h
// - Description: This class is used to assemble both cell matrix and cell residual
//                for reaction source terms in the catalyst layers for various equation classes
// - Developers: Madhur Bhaiya, Marc Secanell, Valentin N. Zingan
//
// ----------------------------------------------------------------------------

#ifndef _FCST_FUELCELLSHOP_EQUATION_REACTION_SOURCE_TERMS_H_
#define _FCST_FUELCELLSHOP_EQUATION_REACTION_SOURCE_TERMS_H_

#include <equations/equation_base.h>
#include <equations/reaction_heat.h>
#include <equations/reaction_source_terms_base.h>

namespace FuelCellShop
{
    namespace Equation
    {

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



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

            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(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 >& dst);

            virtual void print_equation_info() const;

            inline bool get_irrev_heat_ccl() const
            {
                return irrev_heat_ccl;
            }

            inline bool get_irrev_heat_acl() const
            {
                return irrev_heat_acl;
            }

            inline bool get_rev_heat() const
            {
                return rev_heat;
            }

            inline double get_factor_rev_heat_ccl() const
            {
                return factor_rev_heat_ccl;
            }

            inline bool get_water_vap_heat_ccl() const
            {
                return water_vap_heat_ccl;
            }

            inline FuelCellShop::Kinetics::BaseKinetics* get_cathode_kinetics() const
            {
                return cathode_kinetics;
            }

            inline FuelCellShop::Kinetics::BaseKinetics* get_anode_kinetics() const
            {
                return anode_kinetics;
            }


        protected:



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




            bool irrev_heat_ccl;

            bool irrev_heat_acl;

            bool rev_heat;

            double factor_rev_heat_ccl;

            bool water_vapour_phase;

            bool water_vap_heat_ccl;
            
            std::string equation_name_liquid_water;




            FuelCellShop::Kinetics::BaseKinetics* cathode_kinetics;

            FuelCellShop::Kinetics::BaseKinetics* anode_kinetics;




            FuelCellShop::Equation::ReactionHeat* cathode_reactionheat;

            FuelCellShop::Equation::ReactionHeat* anode_reactionheat;




            VariableInfo x_oxygen;

            VariableInfo x_hydrogen;
            
            VariableInfo x_water;

            VariableInfo phi_s;

            VariableInfo lambda;

            VariableInfo phi_m;

            VariableInfo t_rev;

            VariableInfo s_liquid_water;
            
            VariableInfo p_liquid_water;

            double F;

            double M_water;




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

            std::vector< std::vector<double> > phi_xOxygen_cell;

            std::vector< std::vector<double> > phi_xHydrogen_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< std::vector<double> > phi_s_cell;
            
            std::vector< std::vector<double> > phi_p_cell;

            std::vector<double> current_cell;

            std::vector<double> heat_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_dxHydrogen_cell;
            
            std::vector<double> dcurrent_dxWater_cell;

            std::vector<double> dheat_dphiS_cell;

            std::vector<double> dheat_dphiM_cell;

            std::vector<double> dheat_dT_cell;

            std::vector<double> dheat_dxOxygen_cell;

            std::vector<double> dheat_dxHydrogen_cell;
            
            std::vector<double> dheat_dxWater_cell;

            double active_area_cell;

            double factor_protontranseq_cell;

            double factor_electrontranseq_cell;

            double factor_oxygentranseq_cell;
            
            double factor_hydrogentranseq_cell;
            
            double factor_watertranseq_cell;

            double factor_saturationtranseq_cell;
            
            double factor_capillarytranseq_cell;




            bool cell_matrix_counter;

            bool cell_residual_counter;

            unsigned int last_iter_cell;
            
            
             };
             
        } // Equation
             
} // FuelCellShop
        
#endif