agglomerate_ionomer_analytical.h

Go to the documentation of this file.

// ----------------------------------------------------------------------------
//
// FCST: Fuel Cell Simulation Toolbox
//
// Copyright (C) 2006-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: agglomerate_ionomer_analytical.h
// - Description: Used to solve a system of equations representing a spherical ionomer-filled agglomerate.
// - Developers: Peter Dobson <pdobson@ualberta.ca>
//               Marc Secanell Gallart, University of Alberta
// - $Id: agglomerate_ionomer_analytical.h 2605 2014-08-15 03:36:44Z secanell $
//
// ----------------------------------------------------------------------------

#ifndef FUEL_CELL__IONOMER_AGGLOMERATE_ANALYTICAL__H
#define FUEL_CELL__IONOMER_AGGLOMERATE_ANALYTICAL__H


//------------------------------
// FUEL CELL DECLARATIONS
//-----------------------------

#include <microscale/agglomerate_base.h>
#include <reactions/tafel_kinetics.h>
#include <application_core/fcst_variables.h>
#include <reactions/dual_path_kinetics.h>
#include <materials/catalyst_base.h>


namespace FuelCellShop
{
    namespace MicroScale
    {
        class IonomerAgglomerateAnalytical : public AgglomerateBase
        {
        public:


            static const std::string concrete_name;
            
            /*
             * Set the composition and structure of the agglomerate
             */

            
            virtual SolutionMap compute_current ( );
            
            virtual std::vector<double> compute_derivative_current ();
            
            virtual std::string get_name(){
                    return concrete_name;
                }

            virtual double aux_volume_fraction(){
                return 0;
            }



        protected:
            /*
             * Virtual function for returning film thickness in nano meters.
             *
             */
            virtual double get_film_thickness(){
                return delta_agg*1e7;
            }

            /*
             * Virtual function for returning agglomerate radius in nano meters.
             *
             */
            virtual double get_radius(){
                return r_agg*1e7;
            }

            /*
             * Set the composition and structure of the agglomerate
             */
            virtual void set_structure ();

            //name Instance Delivery (Prototype)
            static IonomerAgglomerateAnalytical const* PROTOTYPE;
            
            virtual boost::shared_ptr<FuelCellShop::MicroScale::MicroScaleBase> create_replica ()
            {
                return boost::shared_ptr<FuelCellShop::MicroScale::MicroScaleBase> (new FuelCellShop::MicroScale::IonomerAgglomerateAnalytical ());
            }
            
            IonomerAgglomerateAnalytical ();
            IonomerAgglomerateAnalytical(std::string concrete_name);

            /*
             * Protected virtual member function for declaring parameters, pure
             * in MicroScaleBase, implemented here in IonomerAgglomerateAnalytical.
             * Calls parent AgglomerateBase
             */
            virtual void declare_parameters (ParameterHandler &param) const
            {
                AgglomerateBase::declare_parameters(param);
                param.enter_subsection(concrete_name);{
                    //No parameters at the moment
                }
                param.leave_subsection();
            }

            /*
             * Protected virtual member function for initializing parameters, pure
             * in MicroScaleBase, implemented here in IonomerAgglomerateAnalytical.
             * Calls parent AgglomerateBase
             */
            virtual void initialize (ParameterHandler &param)
            {
                AgglomerateBase::initialize(param);
                param.enter_subsection(concrete_name);{
                    //No parameters at the moment
                }
                param.leave_subsection();
            }

        private:
            
            double compute_Er (const double k_c, const double D);
            
            double compute_dEr (const double k_c,const double dk_c, const double D);
            
            /* Private member function to check kinetics are appropriate for analytical formulation.
             * Throws exception in event of inappropriate kinetic conditions.
             */
            void check_kinetics();

            /*bool to monitor if we have checked kinetic conditions */
            bool checked_kinetics;


        };// class IonomerAgglomerate2
        
    } // namespace Layer
    
} // namespace FuelCellShop

#endif