response_water_capillary.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_water_capillary.h
// - Description: This is a file used to output the capillary pressure related responses.
// - Developers: J.Zhou and M. Secanell 2015    University of Alberta
//
// ----------------------------------------------------------------------------

#ifndef _FUELCELLSHOP__RESPONSE_WATER_CAPILLARY_H
#define _FUELCELLSHOP__RESPONSE_WATER_CAPILLARY_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 <materials/PureLiquid.h>
#include <layers/base_layer.h>
#include <layers/catalyst_layer.h>
#include <layers/gas_diffusion_layer.h>

#include <postprocessing/base_response.h>
#include <utils/fcst_units.h>
#include <utils/fcst_constants.h>

using namespace dealii;

namespace FuelCellShop
{
    namespace PostProcessing
    {
        template <int dim>
        class PhaseChangeResponse : public BaseResponse<dim>
        {
        public:            

            PhaseChangeResponse(const FuelCell::SystemManagement& sm)
            :
            BaseResponse<dim>(sm)
            {}
            
            ~PhaseChangeResponse() {}
            
            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:
            FuelCellShop::Equation::VariableInfo t_rev;
            
            FuelCellShop::Equation::VariableInfo p_liquid_water;
            
            FuelCellShop::Equation::VariableInfo x_water;
            
            double evaporation_rate_constant;
            
            double condensation_rate_constant;
        };
        
        template <int dim>
        class BoundaryLiquidResponse : public BaseBoundaryResponse<dim>
        {
        public:            

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


            void compute_responses(const typename DoFApplication<dim>::FaceInfo& bdry_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>::FaceInfo& bdry_info,
                                   std::map<FuelCellShop::PostProcessing::ResponsesNames, double>& resp) const;
            //
            //
        private:
            FuelCellShop::Equation::VariableInfo p_liquid_water;
            FuelCellShop::Equation::VariableInfo t_rev;
        };
        
        template <int dim>
        class BoundaryVaporResponse : public BaseBoundaryResponse<dim>
        {
        public:            

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


            void compute_responses(const typename DoFApplication<dim>::FaceInfo& bdry_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>::FaceInfo& bdry_info,
                                   std::map<FuelCellShop::PostProcessing::ResponsesNames, double>& resp) const;
            //
            //
        private:
            FuelCellShop::Equation::VariableInfo x_water;
            FuelCellShop::Equation::VariableInfo p_liquid_water;
            FuelCellShop::Equation::VariableInfo t_rev;
        };
    }
}
#endif