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 1460 2013-08-30 22:28:49Z madhur $ 
//
//---------------------------------------------------------------------------


#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 set_parameters (const std::vector<std::string>& name_dvar, 
                                 const std::vector<double>& value_dvar, 
                                 ParameterHandler& param);
            
            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 oxygen_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;

            
        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_sorption = method;
            }
            
        };
    }
}
#endif