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 2605 2014-08-15 03:36:44Z secanell $ 
//
//---------------------------------------------------------------------------

#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();
            
                                 
            virtual void initialize (ParameterHandler &param);
            

            
            virtual void effective_gas_diffusivity(const double&, const double&, double&) const;
            virtual void effective_gas_diffusivity(const double&, const double&, Tensor<2,dim>&) const;
            
            virtual void effective_gas_diffusivity(std::vector< Tensor<2,dim> >& ) const;
            virtual void derivative_effective_gas_diffusivity(std::map< VariableNames, std::vector< Tensor<2,dim> > >&) 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;
            
            virtual void effective_thermal_conductivity(std::vector< Tensor<2,dim> >&) const;
            virtual void derivative_effective_thermal_conductivity(std::vector< Tensor<2,dim> >&) const;
            
            virtual void liquid_permeablity(std::vector< Tensor<2,dim> >& ) const;
            virtual void derivative_liquid_permeablity(std::map< VariableNames, std::vector< Tensor<2,dim> > >& ) const;
            
            virtual void pcapillary(std::vector<double>&) const;
            virtual void dpcapillary_dsat(std::vector<double> &) const;
            virtual void derivative_dpcapillary_dsat(std::map< VariableNames, std::vector<double> > &) const;
            
            virtual void interfacial_surface_area(std::vector<double>&) const;
            virtual void derivative_interfacial_surface_area(std::map< VariableNames, std::vector<double> >&) const;
            

        private:

            DesignFibrousGDL(const std::string& name);
                        

            virtual void declare_parameters (const std::string& name, 
                                             ParameterHandler &param) const;



            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::string method_eff_property_thermal;
            
            std::string method_rel_liquid_permeability;
            double s_irr;
            
            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;
            
            std::vector<double> sigma_e;
            
            std::vector<double> k_thermal;
            
            std::vector<double> abs_permeability;
            
            std::string method_capillary_function;
            double compaction_pressure;
            double PTFE_loading;
            double kumbur_factor;
            
            


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

#endif