response_current_density.h

Go to the documentation of this file.

// ----------------------------------------------------------------------------
//
// FCST: Fuel Cell Simulation Toolbox
//
// Copyright (C) 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: responses.h
// - Description: This is a file containing the declaration of several FCST response evaluators
// - Developers: Marc Secanell Gallart,    University of Alberta
// - $Id: response_current_density.h 2605 2014-08-15 03:36:44Z secanell $
//
// ----------------------------------------------------------------------------

#ifndef _FUELCELLSHOP__RESPONSE_CURRENT_DENSITY_H
#define _FUELCELLSHOP__RESPONSE_CURRENT_DENSITY_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>
#include <exception>      // std::exception

// Include OpenFCST routines:
#include <application_core/fcst_variables.h>
#include <application_core/system_management.h>

#include "equations/equation_base.h"

#include "layers/base_layer.h"
#include "layers/catalyst_layer.h"

#include "postprocessing/base_response.h"

using namespace dealii;

namespace FuelCellShop
{

    namespace PostProcessing
    {

        template <int dim>
        class ORRCurrentDensityResponse : public BaseResponse<dim>
        {
        public:            

            ORRCurrentDensityResponse(const FuelCell::SystemManagement& sm)
            :
            BaseResponse<dim>(sm)
            {}
            
            ~ORRCurrentDensityResponse() {}
            
            void declare_parameters(ParameterHandler& param) const;
            void initialize(ParameterHandler& param);


            void compute_responses(const typename DoFApplication<dim>::CellInfo& info,
                                   FuelCellShop::Layer::BaseLayer<dim>* const layer, 
                                   std::map<FuelCellShop::PostProcessing::ResponsesNames, double>& resp) const;
            void compute_responses(std::vector< FuelCellShop::SolutionVariable > solution_variables,
                                   const typename DoFApplication<dim>::CellInfo& info,
                                   FuelCellShop::Layer::BaseLayer<dim>* const layer, 
                                   std::map<FuelCellShop::PostProcessing::ResponsesNames, double>& resp) const;
            //
            //
        private:
            double S_CL;
            double V_CL;
            FuelCellShop::Equation::VariableInfo xi;
            FuelCellShop::Equation::VariableInfo phiM;
            FuelCellShop::Equation::VariableInfo phiS;
        };
        
        
        //===============================================================================================================
        //===============================================================================================================
        //===============================================================================================================
        //===============================================================================================================
        
        template <int dim>
        class HORCurrentDensityResponse : public BaseResponse<dim>
        {
        public:            

            HORCurrentDensityResponse(const FuelCell::SystemManagement& sm)
            :
            BaseResponse<dim>(sm)
            {}
            
            ~HORCurrentDensityResponse() {}
            
            void declare_parameters(ParameterHandler& param) const;
            void initialize(ParameterHandler& param);


            void compute_responses(const typename DoFApplication<dim>::CellInfo& info,
                                   FuelCellShop::Layer::BaseLayer<dim>* const layer, 
                                   std::map<FuelCellShop::PostProcessing::ResponsesNames, double>& resp) const;
            void compute_responses(std::vector< FuelCellShop::SolutionVariable > solution_variables,
                                   const typename DoFApplication<dim>::CellInfo& info,
                                   FuelCellShop::Layer::BaseLayer<dim>* const layer, 
                                   std::map<FuelCellShop::PostProcessing::ResponsesNames, double>& resp) const;
            //
            //
        private:
            double S_CL;
            double V_CL;
            //FuelCellShop::Equation::VariableInfo x_H2;
            FuelCellShop::Equation::VariableInfo phiM;
            FuelCellShop::Equation::VariableInfo phiS;
        };
        
    }   
}

#endif