solid_layer.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//
//    FCST: Fuel Cell Simulation Toolbox
//
//    Copyright (C) 2011-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: solid_layer.h
//    - Description: Solid layer class
//    - Developers: M. Secanell and Jie Zhou
//    - $Id: solid_layer.h 2605 2014-08-15 03:36:44Z secanell $
//
//---------------------------------------------------------------------------

#ifndef _FUELCELLSHOP__SOLID_LAYER_H
#define _FUELCELLSHOP__SOLID_LAYER_H

// FCST classes
#include <materials/base_material.h>
#include <layers/base_layer.h>
#include <materials/carbon.h>
#include <materials/fiber_base.h>
#include <materials/carbon_fiber.h>
#include <materials/PureSolid.h>

// Include deal.II classes
#include <deal.II/base/parameter_handler.h>
#include <deal.II/base/point.h>
#include <deal.II/base/function.h>
#include <deal.II/lac/vector.h>
#include <deal.II/fe/fe_values.h>

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

using namespace dealii;

namespace FuelCellShop
{
    namespace Layer
    {
        template <int dim>
        class SolidLayer : 
        public  BaseLayer<dim>
        {
        public:
            


             virtual void effective_electron_conductivity(Tensor<2,dim>& ) const;
            
             virtual void effective_electron_conductivity(std::vector<Tensor<2,dim>>& ) const;
            
            virtual void derivative_effective_electron_conductivity(std::map<VariableNames, Tensor<2,dim>>& ) const;
             
             virtual void derivative_effective_electron_conductivity(std::vector< std::map<VariableNames, Tensor<2,dim> > >& ) const;
            
             virtual void effective_thermal_conductivity(Tensor<2,dim>& ) const;
            
              virtual void effective_thermal_conductivity(std::vector< Tensor<2,dim> >&) const;
            
            virtual void derivative_effective_thermal_conductivity(std::map<VariableNames, Tensor<2,dim> >& ) const;
            
             virtual void derivative_effective_thermal_conductivity(std::vector< std::map<VariableNames, Tensor<2,dim> > >& ) const;
            
            inline void set_solution (std::map<VariableNames, SolutionVariable>& sol)
            {
                
                for (std::map<VariableNames, SolutionVariable>::iterator iterator = sol.begin();
                     iterator != sol.end();
                     iterator++)
                     {
                         Assert( iterator->second.get_variablename() == iterator->first, ExcMessage("Wrong solution variable passed in PureSolid::set_solution.") );
                     }
                     
                     this->Solutions = sol;
            }
            


            const std::type_info& get_base_type() const
            {
                return typeid(SolidLayer<dim>);
            }
            
//             /**
//              * Class test
//              */
//             virtual void test_class();

        protected:
                     


            SolidLayer(const std::string &section_layer_name, ParameterHandler& param);
            
            ~SolidLayer();
            
            void declare_parameters (const std::string& solidlayer_section_name,
                                     ParameterHandler &param) const;
                                     
            void initialize (ParameterHandler &param) ;
            
            

            
            std::vector<unsigned int> boundary_ids;
            
            boost::shared_ptr<FuelCellShop::Material::PureSolid>  solid;
            
            std::string concrete_solid_name;

            SolutionVariable T_vector;
            
            std::map <VariableNames, SolutionVariable> Solutions;
            
            /*
             * Define an iterator for the mapping from VariableNames to SolutionVariable
             */
            std::map <VariableNames, SolutionVariable>::iterator solution_iterator;
            
            /*
             * Define an VariableNames of Enumerator of VariableName_of_REV
             */            
            VariableNames temperature_enum;
            
        };
    }
}  // FuelCellShop

#endif // _FUELCELLSHOP__SOLID_LAYER_H