design_fibrous_GDL.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: design_fibrous_GDL.h
//    - Description: Class used to represent a fibrous GDL where effective properties are computed based on the porosity etc.
//    - Developers: M. Secanell and Madhur Bhaiya
//    - $Id: design_fibrous_GDL.h 1460 2013-08-30 22:28:49Z madhur $ 
//
//---------------------------------------------------------------------------

#ifndef _FUELCELLSHOP__DESIGN_FIBROUS_GDL_H
#define _FUELCELLSHOP__DESIGN_FIBROUS_GDL_H

// FCST classes
#include <base_material.h>
#include <gas_diffusion_layer.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 STL
#include<cmath>
#include<iostream>

using namespace dealii;

namespace FuelCellShop
{
    namespace Layer
    {
        template <int dim>
        class DesignFibrousGDL : 
        public  GasDiffusionLayer<dim>
        {
        public:

            static const std::string concrete_name;
            


            DesignFibrousGDL();
            
            DesignFibrousGDL(const std::string& name);
            
            ~DesignFibrousGDL()
            {}
            
            void declare_parameters (ParameterHandler &param) const
            {
                declare_parameters(this->name, param);
            };
            
            virtual void set_parameters (const std::vector<std::string>& name_dvar,
                                         const std::vector<double>& value_dvar,
                                         ParameterHandler &param) const {}
                                 
            virtual void initialize (ParameterHandler &param);
            


            virtual void effective_gas_diffusivity(const double& prop, double& prop_eff) const;
            virtual void effective_gas_diffusivity(const double& prop, Tensor<2,dim>& prop_eff) const;
            
            virtual void effective_gas_diffusivity(Table< 2, Tensor<2,dim> > &D_eff) const;
            
            virtual void effective_electron_conductivity(double&) const;
            
            virtual void effective_electron_conductivity(Tensor<2,dim>&) const;
            

        private:


            virtual void declare_parameters (const std::string& name, 
                                             ParameterHandler &param) const;
                                             
            virtual void set_parameters (const std::vector<std::string>& name_dvar,
                                         const std::vector<double>& value_dvar,
                                         const std::string& name,
                                         ParameterHandler &param) const 
            {
                const std::type_info& info = typeid(*this);
                deallog << "Pure function " << __FUNCTION__
                << " called in Class "
                << info.name() << std::endl;                            
                                             
            }



            virtual boost::shared_ptr<FuelCellShop::Layer::GasDiffusionLayer<dim> > create_replica (const std::string &name)
            {
                return boost::shared_ptr<FuelCellShop::Layer::GasDiffusionLayer<dim> > (new FuelCellShop::Layer::DesignFibrousGDL<dim> (name));
            }
            

            //--------------------------            
            double porosity;
            double solid_phase;
            std::string method_eff_property_pores;
            std::string method_eff_property_solid;
            std::vector<double> sigma_e;
            
            //----------------------
            bool anisotropy;
            std::vector<double> porosity_th;
            std::vector<double> porosity_mu;
            std::vector<double> porosity_gamma;
            std::vector<double> solid_th;
            std::vector<double> solid_mu;
            


            static DesignFibrousGDL<dim> const* PROTOTYPE;
        };
    }
}

#endif