app_diffusion.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//
//    FCST: Fuel Cell Simulation Toolbox
//
//    Copyright (C) 2013 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: app_diffusion.cc
//    - Description:
//    - Developers: M. Sabharwal, M. Secanell
//    - $Id: app_diffusion.h 2605 2014-08-15 03:36:44Z secanell $
//
//---------------------------------------------------------------------------

#ifndef _FUELCELL__APP_DIFFUSION__H
#define _FUELCELL__APP_DIFFUSION__H

// Include deal.II classes
#include "base/parameter_handler.h"
#include "base/function_lib.h"
#include "base/function.h"
#include "base/quadrature_lib.h"

#include "lac/block_vector.h"
#include "lac/full_matrix.h"
#include "lac/solver_cg.h"
#include "lac/solver_gmres.h"
#include "lac/precondition.h"
#include "lac/precondition_block.h"
#include "lac/block_matrix_array.h"
#include "lac/sparse_ilu.h"
#include "lac/sparse_direct.h"

#include "grid/grid_generator.h"
#include "grid/tria_accessor.h"
#include "grid/tria_iterator.h"
#include "grid/tria_boundary_lib.h"

#include "fe/fe_values.h"

#include "numerics/vector_tools.h"
#include "numerics/matrix_tools.h"

#include "boost/shared_ptr.hpp"

// Include appframe classes
#include "block_matrix_application.h"

// Include FuelCell classes
#include "fcst_constants.h"
#include "application_core/optimization_block_matrix_application.h"
#include "solver_utils.h"

#include "operating_conditions.h"
#include "geometries.h"

#include "PureGas.h"
#include "GasMixture.h"
#include "platinum.h"
#include "nafion.h"
#include "carbon.h"

#include "gas_diffusion_layer.h"


#include "new_ficks_transport_equation.h"


#include "linear_solvers.h"

#include "data_out.h"


//Include STL
#include "fstream"
#include "iostream"
#include "sstream"
#include <algorithm>

// Use namespace of deal.II
using namespace dealii;

namespace FuelCell
{

    namespace Application
    {
        //---------------------------------------------------------------------------
        //---------------------------------------------------------------------------
        //---------------------------------------------------------------------------
        template<int dim>
        class AppDiffusion : public OptimizationBlockMatrixApplication<dim>
        {
        public:
            

            
            AppDiffusion( boost::shared_ptr< FuelCell::ApplicationCore::ApplicationData > data = 
            boost::shared_ptr< FuelCell::ApplicationCore::ApplicationData >() );
            
            ~AppDiffusion();
            
            virtual void declare_parameters(ParameterHandler& param);
            
            virtual void initialize(ParameterHandler& param);
            
              virtual void initialize_solution (FEVector& initial_guess,
                                                std::shared_ptr<Function<dim> > initial_function = std::shared_ptr<Function<dim> >());

            

            
            virtual void cell_matrix(MatrixVector&                                 cell_matrices,
                                     const typename DoFApplication<dim>::CellInfo& cell_info);
            
            virtual void cell_residual(FuelCell::ApplicationCore::FEVector&          cell_res,
                                       const typename DoFApplication<dim>::CellInfo& cell_info);
            
            

            
            virtual void dirichlet_bc(std::map<unsigned int, double>& boundary_values) const;

            virtual double evaluate (const FuelCell::ApplicationCore::FEVectors& src);
            
            virtual void data_out(const std::string&         filename,
                                  const FEVectors& src);
            
            

            
            virtual void bdry_responses(std::vector<double>&                                                     dst,
                                        const typename FuelCell::ApplicationCore::DoFApplication<dim>::FaceInfo& bdry_info,
                                        const FuelCell::ApplicationCore::FEVector& src);
            
        protected:            


            


            FuelCell::OperatingConditions OC;
            
            boost::shared_ptr<FuelCellShop::Material::PureGas> solute;
            
            
            
            boost::shared_ptr<FuelCellShop::Material::PureGas> solvent;
            
                        
            


            boost::shared_ptr< FuelCellShop::Layer::GasDiffusionLayer<dim> > CGDL;
            
            


            FuelCellShop::Equation::NewFicksTransportEquation<dim> ficks_transport_equation;
            
        private:
        };
        
    } // Application
    
} // FuelCell

#endif