app_cathode.h

Go to the documentation of this file.

// ----------------------------------------------------------------------------
//
// FCST: Fuel Cell Simulation Toolbox
//
// Copyright (C) 2006-2014 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_cathode.h
// - Description: This class describes diffusion in fuel cell cathodes
//                Ficks, one gas
// - Developers: Marc Secanell,      University of Alberta
//               Valentin N. Zingan, University of Alberta
//               
// ----------------------------------------------------------------------------

#ifndef _FCST_APPLICATION_APP_CATHODE_H_
#define _FCST_APPLICATION_APP_CATHODE_H_

//-- OpenFCST
#include <application_core/optimization_block_matrix_application.h>
#include <application_core/system_management.h>
#include <utils/operating_conditions.h>

#include <materials/PureGas.h>
#include <materials/GasMixture.h>
#include <materials/platinum.h>
#include <materials/nafion.h>
#include <materials/carbon.h>

#include <layers/gas_diffusion_layer.h>
#include <layers/micro_porous_layer.h>
#include <layers/catalyst_layer.h>

#include <equations/ficks_transport_equation.h>
#include <equations/electron_transport_equation.h>
#include <equations/proton_transport_equation.h>
#include <equations/reaction_source_terms.h>
#include <postprocessing/data_out.h>
#include <postprocessing/response_current_density.h>

using namespace dealii;
using namespace FuelCell::ApplicationCore;

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

            
            AppCathode( boost::shared_ptr< FuelCell::ApplicationCore::ApplicationData > data = 
            boost::shared_ptr< FuelCell::ApplicationCore::ApplicationData >() );
            
            ~AppCathode();
            
            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 cell_responses(std::vector<double>&                          dst,
                                        const typename DoFApplication<dim>::CellInfo& cell_info,
                                        const FEVector&);
            void global_responses(std::vector<double>& resp,
                                  const FuelCell::ApplicationCore::FEVector& /*src*/);

            
        protected:


            boost::shared_ptr< FuelCellShop::Geometry::GridBase<dim> > grid;
            

            


            FuelCell::OperatingConditions OC;
            
            FuelCellShop::Material::Oxygen solute;
            FuelCellShop::Material::Hydrogen solute_anode;
            FuelCellShop::Material::Nitrogen solvent;
            


            boost::shared_ptr< FuelCellShop::Layer::GasDiffusionLayer<dim> > CGDL;
            
            boost::shared_ptr< FuelCellShop::Layer::MicroPorousLayer<dim> > CMPL;
            
            boost::shared_ptr< FuelCellShop::Layer::CatalystLayer<dim> > CCL;
            
            bool anode;
            


            FuelCellShop::Equation::FicksTransportEquation<dim>* ficks_transport_equation;
            FuelCellShop::Equation::FicksTransportEquation<dim> ficks_transport_equation_cathode;
            FuelCellShop::Equation::FicksTransportEquation<dim> ficks_transport_equation_anode;
            
            FuelCellShop::Equation::ElectronTransportEquation<dim> electron_transport_equation;
            
            FuelCellShop::Equation::ProtonTransportEquation<dim> proton_transport_equation;
            
            FuelCellShop::Equation::ReactionSourceTerms<dim> reaction_source_terms;
            

            FuelCellShop::PostProcessing::ORRCurrentDensityResponse<dim> ORRCurrent;
            FuelCellShop::PostProcessing::HORCurrentDensityResponse<dim> HORCurrent;
            
            

        private:
            virtual void cell_responses_aux(std::vector<double>&                          dst,
                                            const typename DoFApplication<dim>::CellInfo& cell_info,
                                            const FEVector&);
            
            void set_default_parameters_for_application(ParameterHandler &param)
            {
                param.enter_subsection("System management");
                {
                    param.set("Number of solution variables","3");
                    param.enter_subsection("Solution variables");
                    {
                        param.set("Solution variable 1","oxygen_molar_fraction");
                        param.set("Solution variable 2","protonic_electrical_potential");
                        param.set("Solution variable 3","electronic_electrical_potential");
                    }
                    param.leave_subsection();
                    
                    param.enter_subsection("Equations");
                    {
                        param.set("Equation 1","Ficks Transport Equation - oxygen");
                        param.set("Equation 2","Proton Transport Equation");
                        param.set("Equation 3","Electron Transport Equation");
                    }
                    param.leave_subsection();
                }
                param.leave_subsection();
                param.enter_subsection("Discretization");
                {
                    param.set("Element","FESystem[FE_Q(1)^3]");
                }
                param.leave_subsection();
            }
            
        };
        
    } // Application
    
} // FuelCell

#endif