agglomerate_ionomer_sun.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
// C++ Interface: agglomerate_ionomer_sun.h
//
// Description: Used to solve a system of equations representing a
//              spherical ionomer-filled agglomerate.
//
// Author: Marc Secanell, Phil Wardlaw
//         University of Alberta
//
// Copyright: See COPYING file that comes with this distribution
//
//---------------------------------------------------------------------------

#ifndef FUEL_CELL__IONOMER_AGGLOMERATE_SUN__H
#define FUEL_CELL__IONOMER_AGGLOMERATE_SUN__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 IonomerAgglomerateSun : public AgglomerateBase
        {
            public:


                static const std::string concrete_name;


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

                static IonomerAgglomerateSun const* PROTOTYPE;

                virtual boost::shared_ptr<FuelCellShop::MicroScale::MicroScaleBase> create_replica ()
                {
                    return boost::shared_ptr<FuelCellShop::MicroScale::MicroScaleBase> (new FuelCellShop::MicroScale::IonomerAgglomerateSun ());
                }

                IonomerAgglomerateSun ();
                IonomerAgglomerateSun(std::string concrete_name);

                /*
                 * Protected virtual member function for declaring parameters, pure
                 * in MicroScaleBase, implemented here in IonomerAgglomerateSun.
                 * 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 IonomerAgglomerateSun.
                 * 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:

                /* 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;


                double compute_Er (const double k_c, const double D);

                double compute_dEr (const double k_c,const double dk_c, const double D);
                
                double epsilon_V;


        };// class IonomerAgglomerateSun

    } // namespace Layer

} // namespace FuelCellShop

#endif