design_MPL.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_MPL.h
//    - Description: Class used to represent a design MPL where effective properties are computed based on the porosity etc.
//    - Developers: Madhur Bhaiya
//    - $Id: design_MPL.h 2605 2014-08-15 03:36:44Z secanell $ 
//
//---------------------------------------------------------------------------

#ifndef _FUELCELLSHOP__DESIGN_MPL_H
#define _FUELCELLSHOP__DESIGN_MPL_H

// FCST classes
#include "base_material.h"
#include "micro_porous_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 DesignMPL : 
        public  MicroPorousLayer<dim>
        {
        public:


            static const std::string concrete_name;
            


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

                
            
            


            virtual void effective_gas_diffusivity(const double& property,
                                                   const double& saturation,
                                                   double& effective_property) const;                        
            virtual void effective_gas_diffusivity(const double& property,
                                                   const double& saturation,
                                                   Tensor<2,dim>& effective_property) 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(double& ) const;
            
            virtual void effective_thermal_conductivity(Tensor<2,dim>& ) const;
            
            virtual void 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;
            
            
            virtual void effective_transport_property_solid(const double& property, 
                                                    double& effective_property) const;
                                                    
            virtual void effective_transport_property_solid(const Tensor<2,dim>& property, 
                                                    Tensor<2,dim>& effective_property) const;
            

        protected:


            void declare_parameters (const std::string& mpl_section_name, 
                                     ParameterHandler &param) const;
                                     
                                     
            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);
                FcstUtilities::log << "Pure function " << __FUNCTION__
                << " called in Class "
                << info.name() << std::endl;
            }
                                     

            virtual boost::shared_ptr<FuelCellShop::Layer::MicroPorousLayer<dim> > create_replica (const std::string &name)
            {
                return boost::shared_ptr<FuelCellShop::Layer::MicroPorousLayer<dim> > (new FuelCellShop::Layer::DesignMPL<dim> (name));
            }    
            
            static DesignMPL<dim> const* PROTOTYPE;
            
            
            


            //--------------------------
            bool anisotropy;
            double porosity;
            double solid_phase;
            std::string method_eff_property_pores;
            std::string method_eff_property_fibres;
            std::string method_eff_thermal_conductivity;
            double electrical_conductivity;
            double thermal_conductivity;
            Tensor<2,dim> matrix_electrical_conductivity;
            Tensor<2,dim> matrix_thermal_conductivity;
            
            std::string method_rel_liquid_permeability;
            double s_irr;
            
            //----------------------
            std::vector<double> porosity_th;
            std::vector<double> porosity_mu;
            std::vector<double> porosity_gamma;
            std::vector<double> D_O2;
            std::vector<double> D_wv;
            std::vector<double> fibre_th;
            std::vector<double> fibre_mu;
            
            std::vector<double> abs_permeability;
            
            std::string method_capillary_function;
            double compaction_pressure;
            double PTFE_loading;
            double kumbur_factor;
            
        };
    }
}  // FuelCellShop

#endif // _FUELCELLSHOP__DESIGN_MPL_H