dual_path_kinetics.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//
//    FCST: Fuel Cell Simulation Toolbox
//
//    Copyright (C) 2006-13 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: dual_path_kinetics.h
//    - Description:  Dual Path Kinetics model for hydrogen oxidation reaction
//    - Developers: M. Secanell, M. Moore and Madhur Bhaiya
//    - $Id: dual_path_kinetics.h 2605 2014-08-15 03:36:44Z secanell $
//
//---------------------------------------------------------------------------

#ifndef _FUELCELLSHOP__DUAL_PATH_KINETICS_H
#define _FUELCELLSHOP__DUAL_PATH_KINETICS_H

// Include openFCST routines:
#include <reactions/base_kinetics.h>

// Include deal.II classes
#include <deal.II/base/parameter_handler.h>
#include <deal.II/base/point.h>
#include <deal.II/base/function.h>
#include <deal.II/lac/vector.h>
#include <deal.II/fe/fe_values.h>

//Include STL
#include<cmath>
#include<iostream>

using namespace dealii;

namespace FuelCellShop
{
    namespace Kinetics
    {
        class DualPathKinetics : 
        public  BaseKinetics
        {
        public:


            static const std::string concrete_name;         
            


            DualPathKinetics();
            
            DualPathKinetics(const bool);
            
            ~DualPathKinetics();
            
            virtual void declare_parameters(ParameterHandler &param) const;
            
            virtual void initialize(ParameterHandler &param);
            


            virtual void set_reaction_kinetics(const ReactionNames & name)
            {
                if (name == HOR)
                    name_reaction_kinetics = name;

                else
                {
                    const std::type_info& info = typeid(*this);
                    FcstUtilities::log << "Only HOR reaction is to be implemented in " << __FUNCTION__ << " called in Class " << info.name()  << std::endl;
                    exit(1);
                }
            }
            


            virtual void current_density (std::vector<double>&);
            
            virtual void derivative_current (std::map< VariableNames, std::vector<double> >&);
            
                                
        protected: 


            virtual boost::shared_ptr<FuelCellShop::Kinetics::BaseKinetics > create_replica ()
            {
                return boost::shared_ptr<FuelCellShop::Kinetics::BaseKinetics > (new FuelCellShop::Kinetics::DualPathKinetics ( ));
            }
            static DualPathKinetics const* PROTOTYPE;


            
            virtual void init_kin_param()
            {
                Assert( !kin_param_initialized, ExcInternalError() );
                Assert( catalyst != NULL, ExcMessage("Catalyst object not initialized in the DualPathKinetics object.") );
                Assert( catalyst->get_reaction_name() == HOR, ExcMessage("Catalyst object in the DualPathKinetics not set to HOR reaction name.") );
                Assert( phi_m.is_initialized() && phi_s.is_initialized() && T.is_initialized(), ExcMessage("Either phi_m/phi_s/T is not set in the DualPathKinetics object.") );
                Assert( reactants_map.find(hydrogen_concentration) != reactants_map.end(), ExcMessage("Hydrogen concentration is not set in the DualPathKinetics object.") );

                std::vector<VariableNames> names(1, hydrogen_concentration);
                std::map< VariableNames, double > cref_map;
                catalyst->reference_concentration(names, cref_map);
                ref_conc_H2 = cref_map[hydrogen_concentration];
                
                kin_param_initialized = true;
            }
            
            double ref_conc_H2;
            
            double j_0T;
            
            double j_0H;
            
            double potential_constant;
            
            double ref_potential;
            
            
            
        };
        
    } //Kinetics 
    
} //FuelCellShop

#endif //_FUELCELLSHOP__DUAL_PATH_KINETICS_H