response_ohmic_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_ohmic_heat.h
// - Description: This is header file for electronic and protonic ohmic heat response evaluator classes.
// - Developers: Madhur Bhaiya
// - $Id: 
//
// ----------------------------------------------------------------------------

#ifndef _FUELCELLSHOP__RESPONSE_OHMIC_HEAT_H
#define _FUELCELLSHOP__RESPONSE_OHMIC_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/thermal_transport_equation.h"

#include "layers/base_layer.h"
#include "layers/catalyst_layer.h"
#include "layers/gas_diffusion_layer.h"
#include "layers/micro_porous_layer.h"
#include "layers/membrane_layer.h"

#include "postprocessing/base_response.h"

using namespace dealii;

namespace FuelCellShop
{

    namespace PostProcessing
    {

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

            ElectronOhmicHeatResponse(const FuelCell::SystemManagement& sm,
                                      const FuelCellShop::Equation::ThermalTransportEquation<dim>* tte)
            :
            BaseResponse<dim>(sm),
            thermal_equation(tte),
            factor_GDL(0),
            factor_MPL(0),
            factor_CL(0)
            {}
            
            ~ElectronOhmicHeatResponse() {}
            
            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("ElectronOhmicHeatResponse::compute_responses(solution_variables, info, layer, respMap) not implemented");
            }
            //
        private:
            const FuelCellShop::Equation::ThermalTransportEquation<dim>* thermal_equation;
            
            FuelCellShop::Equation::VariableInfo phiS;
            
            unsigned int factor_GDL;
            
            unsigned int factor_MPL;
            
            unsigned int factor_CL;
        };
        
        
        //===============================================================================================================
        //===============================================================================================================
        //===============================================================================================================
        //===============================================================================================================
        
        
        template <int dim>
        class ProtonOhmicHeatResponse : public BaseResponse<dim>
        {
        public:            

            ProtonOhmicHeatResponse(const FuelCell::SystemManagement& sm,
                                    const FuelCellShop::Equation::ThermalTransportEquation<dim>* tte)
            :
            BaseResponse<dim>(sm),
            thermal_equation(tte),
            factor_CL(0),
            factor_ML(0)
            {}
            
            ~ProtonOhmicHeatResponse() {}
            
            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("ProtonOhmicHeatResponse::compute_responses(solution_variables, info, layer, respMap) not implemented");
            }
            //
        private:
            const FuelCellShop::Equation::ThermalTransportEquation<dim>* thermal_equation;
            
            FuelCellShop::Equation::VariableInfo phiM;
            
            FuelCellShop::Equation::VariableInfo lambda;
            
            FuelCellShop::Equation::VariableInfo tRev;
            
            unsigned int factor_CL;
            
            unsigned int factor_ML;
        };
    }   
}

#endif