liquid_transport_equation.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_transport_equation.h
//    - Description: Equation class for liquid water transport using Darcy's law
//    - Developers: Jie Zhou and Marc Secanell
//
//---------------------------------------------------------------------------

#ifndef _FCST_FUELCELLSHOP_EQUATION_LIQUID_PRESSURE_EQUATION_H_
#define _FCST_FUELCELLSHOP_EQUATION_LIQUID_PRESSURE_EQUATION_H_

// FCST includes
#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 <microscale/PSD_base.h>
// STD
#include <sstream>
#include <string>

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

            
            LiquidPressureEquation(FuelCell::SystemManagement& system_management,
                                      boost::shared_ptr< FuelCell::ApplicationCore::ApplicationData > data = boost::shared_ptr< FuelCell::ApplicationCore::ApplicationData >());
            
            virtual ~LiquidPressureEquation();
            
            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 assemble_bdry_matrix(FuelCell::ApplicationCore::MatrixVector&                                 bdry_matrices,
                                              const typename FuelCell::ApplicationCore::DoFApplication<dim>::FaceInfo& bdry_info,
                                              FuelCellShop::Layer::BaseLayer<dim>* const              layer);
            
            virtual void assemble_bdry_residual(FuelCell::ApplicationCore::FEVector&                                     bdry_rhs,
                                                const typename FuelCell::ApplicationCore::DoFApplication<dim>::FaceInfo& bdry_info,
                                                FuelCellShop::Layer::BaseLayer<dim>* const              layer);
            
            

            
            virtual void print_equation_info() const;

            
        protected:
            

            
            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_bdry_constant_data(const typename FuelCell::ApplicationCore::DoFApplication<dim>::FaceInfo& bdry_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_assemblers_bdry_variable_data(const typename FuelCell::ApplicationCore::DoFApplication<dim>::FaceInfo& bdry_info,
                                                            FuelCellShop::Layer::BaseLayer<dim>* const layer);
            
            

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


            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 rho_l;
            
            double k_sat;
            
            double k_sat_bdry;
            
            

            
            double p_cell;
            
            std::vector< Tensor<2,dim> > rhok_mu_cell;
            
            std::vector< Tensor<2,dim> > rhok_mu_cell_bdry;
            
            std::vector< Tensor<2,dim> > rho_mu_dk_dp_cell;
            
            std::vector< Tensor<2,dim> > dT_rhok_mu_cell;

            std::vector< std::vector<double> > phi_p_cell;
            
            std::vector< std::vector<double> > phi_p_cell_bdry;
            
            std::vector< std::vector<double> > phi_xwater_cell;
            
            std::vector< std::vector<double> > phi_T_cell;
                   
            std::vector< std::vector< Tensor<1,dim> > > grad_phi_p_cell;
            
            std::vector< std::vector< Tensor<1,dim> > > grad_phi_T_cell;
            
            
            bool cell_residual_counter;

            unsigned int last_iter_cell;
            
            unsigned int last_iter_bdry;
            
            bool breakthrough_phenomenon;
            
            std::vector <bool> reached_breakthrough_pressure;
            
            std::map<unsigned int, double> capillary_pressure_current_flux_map;
            
            double p_channel_GDL;
            
            double p_channel_GDL_second;
            double temperature;
            
        };
    } // Equation    
} // FuelCellShop

#endif