agglomerate_hybrid_1D.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
// C++ Interface: agglomerate_hybrid_1D.h
//
// Description: Used to solve a system of equations representing a
//                      spherical hybrid agglomerate.
//
// Author: Philip Wardlaw 2014,
//         University of Alberta
//
// Copyright: See COPYING file that comes with this distribution
//
//---------------------------------------------------------------------------

#ifndef FUEL_CELL__HYBRID_AGGLOMERATE_1D__H
#define FUEL_CELL__HYBRID_AGGLOMERATE_1D__H


//------------------------------
// Inherits everything necessary from water filled agglomerate
//------------------------------

#include <microscale/agglomerate_water_1D.h>


namespace FuelCellShop
{
namespace MicroScale
{

class HybridAgglomerate : public WaterAgglomerate
{
public:



    static const std::string concrete_name;

    virtual std::string get_name(){
        return concrete_name;
    }

    virtual double aux_volume_fraction(){

        return WaterAgglomerate::aux_volume_fraction()*(1-hybrid_core_volume_fraction);

    }

    virtual void print_properties();


protected:

    void set_structure ( );

    HybridAgglomerate ( int verbose = 1 );
    HybridAgglomerate (std::string concrete_name);


    /*
     * Protected virtual member function for declaring parameters
     */
    virtual void declare_parameters (ParameterHandler &param) const
    {
        NumericalAgglomerateBase::declare_parameters(param);
        water.declare_parameters(param);

        param.enter_subsection(concrete_name);{
            param.declare_entry("Hybrid core fraction", "0.01",
                    Patterns::Double(0.0,0.75),
                    "Fraction of the agglomerate core filled with ionomer.");

        }
        param.leave_subsection();




    }

    /*
     * Protected virtual member function for initializing parameters
     */
    virtual void initialize (ParameterHandler &param)
    {
        NumericalAgglomerateBase::initialize(param);
        water.initialize(param);
        param.enter_subsection(concrete_name);{
            hybrid_core_volume_fraction = param.get_double("Hybrid core fraction");
        }
        param.leave_subsection();


    }

    double compute_thickness_agg ();

    virtual double compute_epsilon_agg ();



    static HybridAgglomerate const* PROTOTYPE;

    virtual boost::shared_ptr<FuelCellShop::MicroScale::MicroScaleBase> create_replica ()
    {
        return boost::shared_ptr<FuelCellShop::MicroScale::MicroScaleBase> (new FuelCellShop::MicroScale::HybridAgglomerate());
    }


    void fsub ( double &, double [], double [], double [] );

    static void fsub_wrapper ( double &, double [], double [], double [] );

    void dfsub ( double &, double [], double [], double [] );

    static void dfsub_wrapper ( double &, double [], double [], double [] );




private:

    double hybrid_interface;
    double hybrid_core_volume_fraction;
};

}//namespace Layer
}//namespace FuelCellShop

#endif