data_out.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
//    - $Id: data_out.h 2605 2014-08-15 03:36:44Z secanell $
//
//---------------------------------------------------------------------------

#ifndef _FUELCELLSHOP__POSTPROCESSING__DATAOUT_H
#define _FUELCELLSHOP__POSTPROCESSING__DATAOUT_H

// Include deal.II classes
#include <deal.II/numerics/data_out.h>

//Include STL
#include <boost/shared_ptr.hpp>
#include <mutex>
#include <typeinfo>

// Include OpenFCST routines:

#include <layers/porous_layer.h>
#include <layers/catalyst_layer.h>
#include <layers/gas_diffusion_layer.h>
#include <layers/micro_porous_layer.h>

#include <materials/PureGas.h>

#include <utils/operating_conditions.h>

using namespace dealii;

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


            ORRCurrentDensityDataOut(FuelCell::SystemManagement* ,
                                     boost::shared_ptr< FuelCellShop::Layer::CatalystLayer<dim> > ,
                                     FuelCell::OperatingConditions* );
            virtual ~ORRCurrentDensityDataOut() {}
            
            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;
           
            void set_oxygen_density_name(std::string name)
            {
                oxygen_density_name = name;
            }
            
            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;
            
            mutable std::mutex catalyst_layer_mutex;
            
            std::string oxygen_density_name;
        };
        
        //-------------------------------------------------------------------------------------------
        //-------------------------------------------------------------------------------------------
        //-------------------------------------------------------------------------------------------
        
        template <int dim>
        class HORCurrentDensityDataOut 
        : 
        public dealii::DataPostprocessor<dim>
        {
        public:
            


            HORCurrentDensityDataOut(FuelCell::SystemManagement* ,
                                     boost::shared_ptr< FuelCellShop::Layer::CatalystLayer<dim> > ,
                                     FuelCell::OperatingConditions* );
            virtual ~HORCurrentDensityDataOut() {}
            
            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;
            
            void set_hydrogen_density_name(std::string name)
            {
                hydrogen_density_name = name;
            }
            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;
            
            mutable std::mutex catalyst_layer_mutex;
            
            std::string hydrogen_density_name;
        };
       
        //-------------------------------------------------------------------------------------------
        //-------------------------------------------------------------------------------------------
        //-------------------------------------------------------------------------------------------
        
        template <int dim>
        class RelativeHumidityDataOut 
        : 
        public dealii::DataPostprocessorScalar<dim>
        {
        public:
            


            RelativeHumidityDataOut(FuelCell::SystemManagement* ,
                                    std::vector< boost::shared_ptr< FuelCellShop::Layer::PorousLayer<dim> > > ,
                                    FuelCell::OperatingConditions* );
            virtual ~RelativeHumidityDataOut() {}
           
                                                            
            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:
            
            FuelCellShop::Material::WaterVapor water;
            
            FuelCell::SystemManagement* system_management;
            
            std::vector< boost::shared_ptr< FuelCellShop::Layer::PorousLayer<dim> > > porous_layers;
            
            FuelCell::OperatingConditions* opCond;
            
        };
        
        //-------------------------------------------------------------------------------------------
        //-------------------------------------------------------------------------------------------
        //-------------------------------------------------------------------------------------------
        
        template <int dim>
        class CapillaryPressureDataOut 
        : 
        public dealii::DataPostprocessorScalar<dim>
        {
        public:
            


            CapillaryPressureDataOut(FuelCell::SystemManagement* ,
                                     std::vector< boost::shared_ptr< FuelCellShop::Layer::PorousLayer<dim> > > pls,
                                     FuelCell::OperatingConditions* );
            virtual ~CapillaryPressureDataOut() {}
            
            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:
            FuelCell::SystemManagement* system_management;
                                    
            std::vector< boost::shared_ptr< FuelCellShop::Layer::PorousLayer<dim> > > porous_layers;
            
            FuelCell::OperatingConditions* opCond;

        };
        
        template <int dim>
        class SaturationDataOut 
        : 
        public dealii::DataPostprocessorScalar<dim>
        {
        public:
            


            SaturationDataOut(FuelCell::SystemManagement* ,
                                     std::vector< boost::shared_ptr< FuelCellShop::Layer::PorousLayer<dim> > > pls,
                                     FuelCell::OperatingConditions* );
            virtual ~SaturationDataOut() {}
            
            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:
            FuelCell::SystemManagement* system_management;
                                    
            std::vector< boost::shared_ptr< FuelCellShop::Layer::PorousLayer<dim> > > porous_layers;
            
            FuelCell::OperatingConditions* opCond;

        };        
    }
    
}
    


#endif