ICCP.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
// C++ Interface: agglomerate_ionomer_1D.h
//
// Description: Class that solves solid carbon particle, with Pt surface loading,
//              surrounded by ionomer thin film
//
// Authors: Philip Wardlaw,
//          University of Alberta.
//
// Copyright: See COPYING file that comes with this distribution
//
//---------------------------------------------------------------------------

#include<microscale/micro_scale_base.h>

//FCST materials and kinetics classes
#include <materials/catalyst_base.h>
#include <materials/polymer_electrolyte_material_base.h>
#include <reactions/base_kinetics.h>
#include <reactions/double_trap_kinetics.h>
#include <reactions/tafel_kinetics.h>
#include <reactions/dual_path_kinetics.h>
#include <microscale/agglomerate_base.h> //For SphericalAgglomerateGeometry


#ifndef ICCO_H_
#define ICCO_H_


namespace FuelCellShop
{
    namespace MicroScale
    {

    class ICCP: public MicroScaleBase, public SphericalAgglomerateGeometry {

        public:

            static const std::string concrete_name;


            virtual void set_solution(const std::map<VariableNames,SolutionVariable>& sols,const VariableNames& react, const int& index);

            virtual SolutionMap compute_current ( );

            virtual bool has_derivatives(){
                //numerical for the moment
                return false;
            }

            virtual std::string get_name(){
                return concrete_name;
            }

            virtual double aux_volume_fraction(){
                return 0;
            }


            virtual void print_properties();

            virtual void make_thread_safe(ParameterHandler &param, unsigned int thread_index);


        protected:
          
            /*
             * Method for setting up film thickness or porosity depending on parameters
             * called in initialize() or adjust_poly adjust_polydisperse_structure().
             */
            void _initialize_film_porosity();


            typedef FuelCellShop::Layer::MultiScaleCL<deal_II_dimension> CLPropNames;

            ICCP():
            F(Constants::F()),
            pi(Constants::Pi()){
                non_eq_BC = false;
            }

            ICCP(std::string name);


            virtual void set_structure ();
            /*
             * Protected pure virtual member function for declaring parameters, must be implemented
             * by all children. Parameter structure is hierarchical, therefore children
             * should call their parent's declare parameter function from their own declare
             * parameter function.
             */
            virtual void declare_parameters (ParameterHandler &param) const;

            /*
             * Protected pure virtual member function for initializing parameters, must be implemented
             * by all children. Parameter structure is hierarchical, therefore children
             * should call their parent's initialize function from their own initialize function.
             */
            virtual void initialize (ParameterHandler &param);

            virtual boost::shared_ptr<MicroScaleBase> create_replica (){
                return boost::shared_ptr<MicroScaleBase>( new ICCP());
            }
            static ICCP const* PROTOTYPE;

        private:

            double residual(const double & c_inner, const double & c_outer);


            //Stored solutions
            std::vector<SolutionVariable> reactants;
            SolutionVariable proton_pot, electron_pot;

            //Transport factors
            double ActiveArea;      //Scaled to agglomerate surface
            bool non_eq_BC;         //If true then use. See residual()
            double non_eq_BC_coeff;
            double P;               //Pressure
            double H_R_N;            //Henry's number for reactant in Nafion
            double D_R_N;            // diffusion of reactant in bulk nafion

            //Reaction factors
            double molarNumerator;
            VariableNames reactant;
            
            
            //Pointers to assests obtained from the CL
            boost::shared_ptr<FuelCellShop::Kinetics::BaseKinetics> kinetics;
            boost::shared_ptr<FuelCellShop::Material::PolymerElectrolyteBase> electrolyte;
            boost::shared_ptr<FuelCellShop::Material::CatalystBase> catalyst;


  
            //Some constants
            const double F, pi;
        };
    }

}



#endif /* ICCO_H_ */