liquid_water_source_terms.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//
//    FCST: Fuel Cell Simulation Toolbox
//
//    Copyright (C) 2015 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: liquid_water_source_terms.h
//    - Description: Equation class for liquid water source terms
//    - Developers: Jie Zhou, Marc Secanell
//
//---------------------------------------------------------------------------

#ifndef _FCST_FUELCELLSHOP_EQUATION_LIQUID_SOURCE_TERMS_EQUATION_H_
#define _FCST_FUELCELLSHOP_EQUATION_LIQUID_SOURCE_TERMS_EQUATION_H_

// FCST includes
#include "utils/fcst_constants.h"
#include "utils/fcst_utilities.h"
#include "utils/fcst_units.h"
#include "materials/PureLiquid.h"

#include "equations/equation_base.h"

#include "layers/gas_diffusion_layer.h"
#include "layers/micro_porous_layer.h"
#include "layers/catalyst_layer.h"
#include "PSD_base.h"
// STD
#include <sstream>
#include <string>

namespace FuelCellShop
{
    namespace Equation
    {
        template<int dim>
        class LiquidSourceEquation : public EquationBase<dim>
        {
        public:
            

            
            LiquidSourceEquation(FuelCell::SystemManagement& system_management,boost::shared_ptr< FuelCell::ApplicationCore::ApplicationData > data = 
            boost::shared_ptr< FuelCell::ApplicationCore::ApplicationData >());
            
            virtual ~LiquidSourceEquation();
            
            virtual void declare_parameters(ParameterHandler& param) const;
            
            virtual void initialize(ParameterHandler& param);
            
            


            virtual void assemble_cell_matrix(FuelCell::ApplicationCore::MatrixVector&                                 cell_matrices,
                                              const typename FuelCell::ApplicationCore::DoFApplication<dim>::CellInfo& cell_info,
                                              FuelCellShop::Layer::BaseLayer<dim>* const              layer);
            
            virtual void assemble_cell_residual(FuelCell::ApplicationCore::FEVector&                                     cell_rhs,
                                                const typename FuelCell::ApplicationCore::DoFApplication<dim>::CellInfo& cell_info,
                                                FuelCellShop::Layer::BaseLayer<dim>* const              layer);
            
            

            
            virtual void print_equation_info() const;
            
            double get_kevap() const
            {
                return k_e;
            }
            
            double get_kcond() const
            {
                return k_c;
            }
            
            void set_relative_humility(double const x_wv) 
            {
                this->Relative_humility = x_wv;
                Relative_humility_is_set = true;
            }
                        
            virtual void adjust_internal_cell_couplings(std::vector< couplings_map >& ) const;

            
        protected:
            inline void get_temp_related_parameters(double Temp,double c_pressure,double &psat,double &dpsat_dT,
                                                    double &h_lv,double &dhlv_dT,double &Kelvin_factor,
                                                    double &psat_kelvin, double &dpsat_kelvin_dpc, double &dpsat_kelvin_dT,
                                                    FuelCellShop::Material::WaterVapor water) const
            {
                psat = water.get_water_vapor_saturation_pressure(Temp);  // saturation pressure [Pa]
                dpsat_dT = water.get_Dwater_vapor_saturation_pressure_Dtemperature(Temp); // derivative of saturation pressure w.r.t temperature
                h_lv = FuelCellShop::Material::LiquidWater::latentVap_heat(Temp); // latent heat of vaporization [J/mol]
                dhlv_dT = FuelCellShop::Material::LiquidWater::deriv_latentVap_heat(Temp); // derivative of latent heat w.r.t temperature 
                Kelvin_factor = std::exp(c_pressure *  M_water * Units::convert(1, Units::PER_C_UNIT3, Units::PER_UNIT2) / (Constants::R() * rho_l * Temp )  );
                
                psat_kelvin = psat * Kelvin_factor;
                dpsat_kelvin_dpc = psat * (M_water * Units::convert(1, Units::PER_C_UNIT3, Units::PER_UNIT2) / (Constants::R() * rho_l * Temp )) * Kelvin_factor;
                dpsat_kelvin_dT = dpsat_dT * Kelvin_factor 
                    + psat * (-1.0 * (c_pressure *  M_water * Units::convert(1, Units::PER_C_UNIT3, Units::PER_UNIT2) / (Constants::R() * rho_l)  ) / std::pow(Temp,2.0)) * Kelvin_factor;
            }
            

            
            virtual void make_assemblers_generic_constant_data();
            
            virtual void make_assemblers_cell_constant_data(const typename FuelCell::ApplicationCore::DoFApplication<dim>::CellInfo& cell_info);
            
            virtual void make_assemblers_cell_variable_data(const typename FuelCell::ApplicationCore::DoFApplication<dim>::CellInfo& cell_info,
                                                            FuelCellShop::Layer::BaseLayer<dim>* const layer);
            
            


            virtual void make_boundary_types(){};
            
            virtual void make_output_types(){};
            


            std::map<unsigned int, double> capillary_pressure_current_flux_map;
            
            VariableInfo p_liquid_water;
            
            VariableInfo x_water;
            
            VariableInfo t_rev;
            
            unsigned int fickswater_blockindex_xwater;
            
            unsigned int fickswater_blockindex_pliquid;
            
            unsigned int fickswater_blockindex_trev;
            
            unsigned int thermal_blockindex_xwater;
            
            unsigned int thermal_blockindex_pliquid;
            
            unsigned int thermal_blockindex_trev;
            
            double M_water;
            
            double k_e;
            double k_c;
            double rho_l;
            
            

            
            double p_cell;
            
            std::vector<double> area_lv_cell;
            
            std::vector<double> darea_lv_dp_cell;

            std::vector< std::vector<double> > phi_p_cell;
            
            std::vector< std::vector<double> > phi_xwater_cell;
            
            std::vector< std::vector<double> > phi_T_cell;
            
            
            bool cell_residual_counter;

            unsigned int last_iter_cell;
            
            double T_constant;
            
            double Relative_humility;
            
            bool Relative_humility_is_set;
            
        };
    } // Equation    
} // FuelCellShop

#endif