response_capillary_related.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: data_out.h
//    - Description: 
//    - Developers: M. Secanell
//
//---------------------------------------------------------------------------

#ifndef _FUELCELLSHOP__POSTPROCESSING__RESPONSE_CAPILARRY_RELATED_H
#define _FUELCELLSHOP__POSTPROCESSING__RESPONSE_CAPILARRY_RELATED_H

// Include OpenFCST routines:
#include <application_core/optimization_block_matrix_application.h>

// Equation Classes
#include <equations/thermal_transport_equation.h>
#include <equations/reaction_source_terms.h>
#include <equations/proton_transport_equation.h>
#include <equations/lambda_transport_equation.h>
#include <equations/electron_transport_equation.h>
#include <equations/sorption_source_terms.h>
#include <equations/ficks_transport_equation.h>
#include <equations/capillary_transport_equation.h>

#include <layers/porous_layer.h>
#include <postprocessing/data_out.h>

#include <materials/PureGas.h>
#include <application_core/dof_application.h>
#include <application_core/matrix_block.h>
#include <application_core/block_matrix_application.h>

#include <utils/operating_conditions.h>

using namespace dealii;

namespace FuelCellShop
{
    namespace PostProcessing
    {
        
        template <int dim>
        class PermeabilityDataOut 
        : 
        public dealii::DataPostprocessor<dim>
        {
        public:
            


            PermeabilityDataOut(FuelCell::SystemManagement* ,
                                     std::vector< boost::shared_ptr< FuelCellShop::Layer::PorousLayer<dim> > >  ,
                                     FuelCell::OperatingConditions* );
            virtual ~PermeabilityDataOut() {}
            
            virtual std::vector<std::string> get_names() const;
            
            virtual std::vector<DataComponentInterpretation::DataComponentInterpretation> get_data_component_interpretation () const;
            
            virtual UpdateFlags get_needed_update_flags() const;
           
            virtual void compute_derived_quantities_vector (const std::vector< Vector< double > > &uh,
                                                    const std::vector< std::vector< Tensor< 1, dim > > > & /*duh*/,
                                                    const std::vector< std::vector< Tensor< 2, dim > > > & /*dduh*/,
                                                    const std::vector< Point< dim > > & /*normals*/,
                                                    const std::vector< Point<dim> > & /*evaluation_points*/,
                                                    const types::material_id & mat_id,
                                                    std::vector< Vector< double > > &computed_quantities) const;

                                                            
        private:
            std::vector< boost::shared_ptr< FuelCellShop::Layer::PorousLayer<dim> > > porous_layers;
            
            FuelCell::SystemManagement* system_management;

            FuelCell::OperatingConditions* opCond;
            
        };
        
        
        template <int dim>
        class RHSDataOut 
        : 
        public dealii::DataPostprocessor<dim>
        {
        public:
            


            RHSDataOut(FuelCell::SystemManagement* ,
                                     std::vector< boost::shared_ptr< FuelCellShop::Layer::PorousLayer<dim> > >  ,
                                     FuelCell::OperatingConditions* );
            virtual ~RHSDataOut() {}
            
            virtual std::vector<std::string> get_names() const;
            
            virtual std::vector<DataComponentInterpretation::DataComponentInterpretation> get_data_component_interpretation () const;
            
            virtual UpdateFlags get_needed_update_flags() const;
           
            virtual void compute_derived_quantities_vector (const std::vector< Vector< double > > &uh,
                                                    const std::vector< std::vector< Tensor< 1, dim > > > & /*duh*/,
                                                    const std::vector< std::vector< Tensor< 2, dim > > > & /*dduh*/,
                                                    const std::vector< Point< dim > > & /*normals*/,
                                                    const std::vector< Point<dim> > & /*evaluation_points*/,
                                                    const types::material_id & mat_id,
                                                    std::vector< Vector< double > > &computed_quantities) const;

                                                            
        private:
            std::vector< boost::shared_ptr< FuelCellShop::Layer::PorousLayer<dim> > > porous_layers;
            
            FuelCell::SystemManagement* system_management;

            FuelCell::OperatingConditions* opCond;
            
        };
        
        template <int dim>
        class LHS_Pc_DataOut 
        : 
        public dealii::DataPostprocessor<dim>
        {
        public:
            


            LHS_Pc_DataOut(FuelCell::SystemManagement* ,
                                     std::vector< boost::shared_ptr< FuelCellShop::Layer::PorousLayer<dim> > >  ,
                                     FuelCell::OperatingConditions*
                      );
            virtual ~LHS_Pc_DataOut() {}
            
            virtual std::vector<std::string> get_names() const;
            
            virtual std::vector<DataComponentInterpretation::DataComponentInterpretation> get_data_component_interpretation () const;
            
            virtual UpdateFlags get_needed_update_flags() const;
           
            virtual void compute_derived_quantities_vector (const std::vector< Vector< double > > &uh,
                                                    const std::vector< std::vector< Tensor< 1, dim > > > & duh,
                                                    const std::vector< std::vector< Tensor< 2, dim > > > & /*dduh*/,
                                                    const std::vector< Point< dim > > & /*normals*/,
                                                    const std::vector< Point<dim> > & /*evaluation_points*/,
                                                    const types::material_id & mat_id,
                                                    std::vector< Vector< double > > &computed_quantities) const;

                                                            
        private:
            std::vector< boost::shared_ptr< FuelCellShop::Layer::PorousLayer<dim> > > porous_layers;
            
            FuelCell::SystemManagement* system_management;

            FuelCell::OperatingConditions* opCond;
            
        };
        
          template <int dim>
        class ORRSOURCEDataOut 
        : 
        public dealii::DataPostprocessor<dim>
        {
        public:
            


            ORRSOURCEDataOut(FuelCell::SystemManagement* ,
                                     boost::shared_ptr< FuelCellShop::Layer::CatalystLayer<dim> > ,
                                     FuelCell::OperatingConditions* );
            virtual ~ORRSOURCEDataOut() {}
            
            virtual std::vector<std::string> get_names() const;
            
            virtual std::vector<DataComponentInterpretation::DataComponentInterpretation> get_data_component_interpretation () const;
            
            virtual UpdateFlags get_needed_update_flags() const;
           
            virtual void compute_derived_quantities_vector (const std::vector< Vector< double > > &uh,
                                                    const std::vector< std::vector< Tensor< 1, dim > > > & /*duh*/,
                                                    const std::vector< std::vector< Tensor< 2, dim > > > & /*dduh*/,
                                                    const std::vector< Point< dim > > & /*normals*/,
                                                    const std::vector< Point<dim> > & /*evaluation_points*/,
                                                    const types::material_id & mat_id,
                                                    std::vector< Vector< double > > &computed_quantities) const;

                                                            
        private:
            boost::shared_ptr< FuelCellShop::Layer::CatalystLayer<dim> > catalyst_layer;
            
            FuelCell::SystemManagement* system_management;

            FuelCell::OperatingConditions* opCond;
            
        };
        
        template <int dim>
        class RHSData_KelvinOut 
        : 
        public dealii::DataPostprocessor<dim>
        {
        public:
            


            RHSData_KelvinOut(FuelCell::SystemManagement* ,
                                     std::vector< boost::shared_ptr< FuelCellShop::Layer::PorousLayer<dim> > >  ,
                                     FuelCell::OperatingConditions* );
            virtual ~RHSData_KelvinOut() {}
            
            virtual std::vector<std::string> get_names() const;
            
            virtual std::vector<DataComponentInterpretation::DataComponentInterpretation> get_data_component_interpretation () const;
            
            virtual UpdateFlags get_needed_update_flags() const;
           
            virtual void compute_derived_quantities_vector (const std::vector< Vector< double > > &uh,
                                                    const std::vector< std::vector< Tensor< 1, dim > > > & /*duh*/,
                                                    const std::vector< std::vector< Tensor< 2, dim > > > & /*dduh*/,
                                                    const std::vector< Point< dim > > & /*normals*/,
                                                    const std::vector< Point<dim> > & /*evaluation_points*/,
                                                    const types::material_id & mat_id,
                                                    std::vector< Vector< double > > &computed_quantities) const;

                                                            
        private:
            std::vector< boost::shared_ptr< FuelCellShop::Layer::PorousLayer<dim> > > porous_layers;
            
            FuelCell::SystemManagement* system_management;

            FuelCell::OperatingConditions* opCond;
            
        };
        
               
        template <int dim>
        class RHSData_KelvinOut_NoTem 
        : 
        public dealii::DataPostprocessor<dim>
        {
        public:
            


            RHSData_KelvinOut_NoTem(FuelCell::SystemManagement* ,
                                     std::vector< boost::shared_ptr< FuelCellShop::Layer::PorousLayer<dim> > >  ,
                                     FuelCell::OperatingConditions* );
            virtual ~RHSData_KelvinOut_NoTem() {}
            
            virtual std::vector<std::string> get_names() const;
            
            virtual std::vector<DataComponentInterpretation::DataComponentInterpretation> get_data_component_interpretation () const;
            
            virtual UpdateFlags get_needed_update_flags() const;
           
            virtual void compute_derived_quantities_vector (const std::vector< Vector< double > > &uh,
                                                    const std::vector< std::vector< Tensor< 1, dim > > > & /*duh*/,
                                                    const std::vector< std::vector< Tensor< 2, dim > > > & /*dduh*/,
                                                    const std::vector< Point< dim > > & /*normals*/,
                                                    const std::vector< Point<dim> > & /*evaluation_points*/,
                                                    const types::material_id & mat_id,
                                                    std::vector< Vector< double > > &computed_quantities) const;

                                                            
        private:
            std::vector< boost::shared_ptr< FuelCellShop::Layer::PorousLayer<dim> > > porous_layers;
            
            FuelCell::SystemManagement* system_management;

            FuelCell::OperatingConditions* opCond;
            
        };
        
    }
    
}
    


#endif