tafel_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: tafel_kinetics.h
//    - Description: Header file for Tafel Kinetics model class.
//    - Developers: M. Moore, M. Bhaiya, M. Secanell and V. Zingan
//    - $Id: tafel_kinetics.h 2605 2014-08-15 03:36:44Z secanell $
//
//---------------------------------------------------------------------------

#ifndef _FUELCELLSHOP__TAFEL_KINETICS_H
#define _FUELCELLSHOP__TAFEL_KINETICS_H

// Include openFCST routines:
#include "base_kinetics.h"
#include "fcst_constants.h"

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

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

using namespace dealii;

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



            TafelKinetics();

            TafelKinetics(const bool);

            ~TafelKinetics();

            virtual void declare_parameters(ParameterHandler& param) const;

            virtual void initialize(ParameterHandler& param);




            static const std::string concrete_name;



            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::TafelKinetics ( ));
            }
            static TafelKinetics const* PROTOTYPE;



            virtual void init_kin_param()
            {
                Assert( !kin_param_initialized, ExcInternalError() );
                Assert( catalyst != NULL, ExcMessage("Catalyst object not initialized in the TafelKinetics object.") );
                Assert( phi_m.is_initialized() && phi_s.is_initialized() && T.is_initialized(), ExcMessage("Either phi_m/phi_s/T is not set in the TafelKinetics object.") );
                Assert( reactants_map.size()>0, ExcMessage("Atleast one reactant should be set in the TafelKinetics object, using set_reactant_concentrations method.") );

                std::vector<VariableNames> names;
                for ( std::map< VariableNames, SolutionVariable >::const_iterator iter=reactants_map.begin(); iter!=reactants_map.end(); ++iter )
                    names.push_back(iter->first);

                catalyst->reference_concentration(names, ref_conc);
                catalyst->reaction_order(names, gamma);
                catalyst->alpha_cathodic(alpha_c);

                kin_param_initialized = true;
            }

            double alpha_c;

            std::map< VariableNames, double > ref_conc;

            std::map< VariableNames, double > gamma;

        };

    } //Kinetics

} //FuelCellShop

#endif //_FUELCELLSHOP__TAFEL_KINETICS_H