nafion.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: nafion.h
//    - Description:  Class to represent and return bulk properties for Nafion (Polymer Electrolyte Membrane)
//    - Developers: M. Secanell (2009-13) and Madhur Bhaiya (2012-13)
//    - Id: $Id: nafion.h 2605 2014-08-15 03:36:44Z secanell $ 
//
//---------------------------------------------------------------------------


#ifndef _FUELCELLSHOP__NAFION_H
#define _FUELCELLSHOP__NAFION_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 "polymer_electrolyte_material_base.h"
#include "PureGas.h"

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

class NafionTest;

namespace FuelCellShop
{
    namespace Material
    {
        class Nafion :
        public PolymerElectrolyteBase
        {
        public:
            
            static const std::string concrete_name;
            

            
            friend class ::NafionTest;
            


            Nafion(const bool);
            
            Nafion(std::string name = "Nafion");
            
            ~Nafion();
            
            virtual void declare_parameters(ParameterHandler &param) const;
            
            virtual void initialize (ParameterHandler &param);
            
            

            
            virtual void sorption_isotherm(std::vector<double>&) const;
            virtual void sorption_isotherm_derivative(std::map < VariableNames, std::vector<double> >&) const;
            
            virtual void proton_conductivity(double&) const;
            virtual void proton_conductivity(std::vector<double>&) const;
            virtual void proton_conductivity_derivative(std::map< VariableNames, std::vector<double> >&) const;
            
            virtual void water_diffusivity(double&) const;
            virtual void water_diffusivity(std::vector<double>&) const;
            virtual void water_diffusivity_derivative(std::map< VariableNames, std::vector<double> >&) const;
            
            virtual void electroosmotic_drag(std::vector<double>&) const;
            virtual void electroosmotic_drag_derivative(std::map< VariableNames, std::vector<double> >&) const;
            
            virtual void thermoosmotic_coeff(std::vector<double>&) const;
            virtual void thermoosmotic_coeff_derivative(std::map< VariableNames, std::vector<double> >&) const;
            
            virtual void oxygen_diffusivity(double&) const;

            virtual void hydrogen_diffusivity(double&) const;
            virtual void oxygen_diffusivity(std::vector<double>&) const;
            virtual void oxygen_diffusivity_derivative(std::map< VariableNames, std::vector<double> >&) const;
            
            virtual void proton_diffusivity(double&) const;

            virtual void sorption_enthalpy(std::vector<double>& ) const;
            virtual void sorption_enthalpy_derivative(std::map< VariableNames, std::vector<double> >&) const;

            

            
            virtual double get_Hlambda(const double&) const;
            
            virtual double get_dHlambda_dT(const double&) const;
            
            virtual double get_d2Hlambda_dT2(const double&) const;
            
            
        private:
            


            virtual boost::shared_ptr<FuelCellShop::Material::PolymerElectrolyteBase > create_replica ()
            {
                return boost::shared_ptr<FuelCellShop::Material::PolymerElectrolyteBase > (new FuelCellShop::Material::Nafion ());
            }
            static Nafion const* PROTOTYPE;
            
            inline void modify_parameters(const unsigned int& index, const std::string& method)
            {
                Assert(index>=1 && index<=6, ExcMessage("Wrong index input in Nafion::modify_parameters."));
                
                if (index == 1)
                    method_conductivity = method;
                else if (index == 2)
                    method_diffusivity = method;
                else if (index == 3)
                    method_electroosmotic_drag = method;
                else if (index == 4)
                    method_thermoosmosis = method;
                else if (index == 5)
                    method_sorption = method;
                else if (index == 6)
                    method_enthalpy_sorption = method;
            }
        };
    }
}
#endif