data_out_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: data_out_capillary.h
//    - Description: This file is used to output all the liquid water related variables for example, permeability, interfacial surface area.
//    - Developers: J. Zhou and M. Secanell 2015
//
//---------------------------------------------------------------------------

#ifndef _FUELCELLSHOP__POSTPROCESSING__DATA_OUT_CAPILARRY_H
#define _FUELCELLSHOP__POSTPROCESSING__DATA_OUT_CAPILARRY_H

// deal.II parent class:
#include <deal.II/numerics/data_out.h>

// OpenFCST example class:
#include <postprocessing/data_out.h>

// Other required OpenFCST classes:
#include <utils/operating_conditions.h>
#include <layers/porous_layer.h>
#include <layers/catalyst_layer.h>

using namespace dealii;

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


            LiquidVelocityPermeabilityDataOut(FuelCell::SystemManagement* ,
                                     std::vector< boost::shared_ptr< FuelCellShop::Layer::PorousLayer<dim> > >  ,
                                     FuelCell::OperatingConditions*
                      );
            virtual ~LiquidVelocityPermeabilityDataOut() {}
            
            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 PhaseChangeDataOut 
        : 
        public dealii::DataPostprocessor<dim>
        {
        public:
            


            PhaseChangeDataOut(FuelCell::SystemManagement* ,
                              std::vector< boost::shared_ptr< FuelCellShop::Layer::PorousLayer<dim> > >  ,
                              FuelCell::OperatingConditions* ,
                              double,
                              double);
            virtual ~PhaseChangeDataOut() {}
            
            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;
            
            double kc;
            
            double ke;
            
        };
        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