response_reaction_heat.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: response_reaction_heat.h
// - Description: This is header file for response class computing heat generated due to electrochemical reactions.
// - Developers: Madhur Bhaiya
// - $Id: 
//
// ----------------------------------------------------------------------------

#ifndef _FUELCELLSHOP__RESPONSE_REACTION_HEAT_H
#define _FUELCELLSHOP__RESPONSE_REACTION_HEAT_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 "equations/reaction_source_terms.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 ORRReactionHeatResponse : public BaseResponse<dim>
        {
        public:            

            ORRReactionHeatResponse(const FuelCell::SystemManagement& sm,
                                    const FuelCellShop::Equation::ReactionSourceTerms<dim>* rst)
            :
            BaseResponse<dim>(sm),
            reaction_source(rst),
            reaction_heat(NULL)
            {}
            
            ~ORRReactionHeatResponse()
            {
                delete reaction_heat;
            }
            
            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>& respMap) 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>& respMap) const
            {
                throw std::runtime_error("ORRReactionHeatResponse::compute_responses(solution_variables, info, layer, respMap) not implemented");
            }
            //
        private:
            const FuelCellShop::Equation::ReactionSourceTerms<dim>* reaction_source;
            
            FuelCellShop::Kinetics::BaseKinetics* kinetics;
            
            FuelCellShop::Equation::ReactionHeat* reaction_heat;
            
            FuelCellShop::Equation::VariableInfo xOxygen;
            
            FuelCellShop::Equation::VariableInfo phiS;
            
            FuelCellShop::Equation::VariableInfo phiM;
            
            FuelCellShop::Equation::VariableInfo tRev;
        };
        
        
        //===============================================================================================================
        //===============================================================================================================
        //===============================================================================================================
        //===============================================================================================================
        
        
        template <int dim>
        class HORReactionHeatResponse : public BaseResponse<dim>
        {
        public:            

            HORReactionHeatResponse(const FuelCell::SystemManagement& sm,
                                    const FuelCellShop::Equation::ReactionSourceTerms<dim>* rst)
            :
            BaseResponse<dim>(sm),
            reaction_source(rst),
            reaction_heat(NULL)
            {}
            
            ~HORReactionHeatResponse()
            {
                delete reaction_heat;
            }
            
            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>& respMap) 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>& respMap) const
            {
                throw std::runtime_error("HORReactionHeatResponse::compute_responses(solution_variables, info, layer, respMap) not implemented");
            }
            //
        private:
            const FuelCellShop::Equation::ReactionSourceTerms<dim>* reaction_source;
            
            FuelCellShop::Kinetics::BaseKinetics* kinetics;
            
            FuelCellShop::Equation::ReactionHeat* reaction_heat;
            
            FuelCellShop::Equation::VariableInfo xHydrogen;
            
            FuelCellShop::Equation::VariableInfo xWater;
            
            FuelCellShop::Equation::VariableInfo phiS;
            
            FuelCellShop::Equation::VariableInfo phiM;
            
            FuelCellShop::Equation::VariableInfo tRev;
        };
    }   
}

#endif