solid_layer.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//    $Id: solid_layer.h 1354 2013-08-17 00:01:22Z secanell $
//
//    Copyright (C) 2009 by Marc Secanell, University of Alberta
//
//    This file is subject to QPL and may not be  distributed
//    without copyright and license information. Please refer
//    to the file deal.II/doc/license.html for the  text  and
//    further information on this license.
//
//---------------------------------------------------------------------------

#ifndef _FUELCELLSHOP__SOLID_LAYER_H
#define _FUELCELLSHOP__SOLID_LAYER_H

// FCST classes
#include <base_material.h>
#include <base_layer.h>
#include <carbon.h>
#include "fiber_base.h"
#include "carbon_fiber.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 SolidLayer : 
      public  BaseLayer<dim>
      {
      public:
        SolidLayer(std::string name);
        
        ~SolidLayer();
        
        void declare_parameters (ParameterHandler &param) const;

        void set_parameters (const std::vector<std::string>& name_dvar,
                             const std::vector<double>& value_dvar,
                             ParameterHandler &param);
        
        void initialize (ParameterHandler &param);
        
        void set_solid_and_compute (FuelCellShop::Material::FiberBase* ,
                                    double );
        
        virtual void effective_electron_conductivity(double& ) const;
        
        
        
        virtual void effective_electron_conductivity(Tensor<2,dim>& ) const
        {
          const std::type_info& info = typeid(*this);
          deallog << "Pure function " << __FUNCTION__
                  << " called in Class "
                  << info.name() << std::endl;
        }
        virtual void effective_electron_conductivity(const double& /*prop*/, Tensor<2,dim>& /*prop_eff*/) const
        {
          const std::type_info& info = typeid(*this);
          deallog << "Pure function " << __FUNCTION__
                  << " called in Class "
                  << info.name() << std::endl;
        }
        

        virtual void derivative_effective_electron_conductivity(std::vector<Tensor<2,dim> >& ) const
        {
          const std::type_info& info = typeid(*this);
          deallog << "Pure function " << __FUNCTION__
                  << " called in Class "
                  << info.name() << std::endl;
        };
        virtual void effective_thermal_conductivity(double& ) const
        {
          const std::type_info& info = typeid(*this);
          deallog << "Pure function " << __FUNCTION__
                  << " called in Class "
                  << info.name() << std::endl;
        };
        virtual void effective_thermal_conductivity(Tensor<2,dim>& ) const
        {
          const std::type_info& info = typeid(*this);
          deallog << "Pure function " << __FUNCTION__
                  << " called in Class "
                  << info.name() << std::endl;
        };
        virtual void effective_thermal_conductivity(const double&, const double&, double&) const
        {
          const std::type_info& info = typeid(*this);
          deallog << "Pure function " << __FUNCTION__
                  << " called in Class "
                  << info.name() << std::endl;
        };
        virtual void effective_thermal_conductivity(std::vector< Tensor<2,dim> >&) const
        {
          const std::type_info& info = typeid(*this);
          deallog << "Pure function " << __FUNCTION__
                  << " called in Class "
                  << info.name() << std::endl;
        };
        virtual void derivative_effective_thermal_conductivity(std::vector<Tensor<2,dim> >& ) const
        {
          const std::type_info& info = typeid(*this);
          deallog << "Pure function " << __FUNCTION__
                  << " called in Class "
                  << info.name() << std::endl;
        };
        
    const std::type_info& get_base_type() const
    {
        return typeid(SolidLayer<dim>);
    }
    
        virtual void test_class();
        
      protected:
        
        std::vector<unsigned int> boundary_ids;
        
        double temperature;
        
        FuelCellShop::Material::FiberBase* solid;
         
        double electrical_conductivity;
        };
  }
}  // FuelCellShop

#endif // _FUELCELLSHOP__GAS_DIFFUSION_LAYER_H