agglomerate_ionomer_1D.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
// C++ Interface: agglomerate_ionomer_1D.h
//
// Description: Used to solve a system of equations representing a
//                      spherical ionomer-filled agglomerate.
//
// Authors: Jason Boisvert <jjb701@mail.usask.ca>, (C) 2011
//                      University of Saskatchewan,
//          Peter Dobson <pdobson@ualberta.ca>,
//          Philip Wardlaw,
//          Michael Moore,
//                  University of Alberta.

// Copyright: See COPYING file that comes with this distribution
//
//---------------------------------------------------------------------------

#ifndef FUEL_CELL__IONOMER_AGGLOMERATE_NUMERICAL_1D__H
#define FUEL_CELL__IONOMER_AGGLOMERATE_NUMERICAL_1D__H


//------------------------------
// FUEL CELL DECLARATIONS
//-----------------------------
#include <microscale/numerical_agglomerate_base.h>
#include <contribs/DAE_wrapper.h>




namespace FuelCellShop
{
namespace MicroScale
{

class IonomerAgglomerate : public NumericalAgglomerateBase, public FuelCell::ApplicationCore::DAEWrapper
{
public:


    static const std::string concrete_name;

    virtual SolutionMap compute_current ( );

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

    virtual double aux_volume_fraction(){
        return 0;
    }

protected:

    void setup_DAE_solver ();

    int cont_tolerance (double start_tol, double end_tol);

    /*
     * Virtual function for returning film thickness in nano meters.
     *
     */
    virtual double get_film_thickness(){
        return delta_agg*1e9;
    }

    /*
     * Virtual function for returning agglomerate radius in nano meters.
     *
     */
    virtual double get_radius(){
        return r_agg*1e9;
    }

    virtual void set_structure ();


    static IonomerAgglomerate const* PROTOTYPE;

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

    IonomerAgglomerate ( int verbose = 1 );
    IonomerAgglomerate(std::string name);

    /*
     * Protected virtual member function for declaring parameters
     */
    virtual void declare_parameters (ParameterHandler &param) const
    {
        NumericalAgglomerateBase::declare_parameters(param);
        param.enter_subsection(concrete_name);{
            param.declare_entry("Proton Conductivity factor", "1.0",
                    Patterns::Double());
            param.declare_entry("Non Equilibrium BC Rate constant", "1.0", Patterns::Double());
            param.declare_entry("Use non equilibrium BC", "false",
                    Patterns::Bool());
        }
        param.leave_subsection();




    }

    /*
     * Protected virtual member function for initializing parameters
     */
    virtual void initialize (ParameterHandler &param)
    {
        NumericalAgglomerateBase::initialize(param);
        param.enter_subsection(concrete_name);{
            cond_factor = param.get_double("Proton Conductivity factor");
            k_O2 = param.get_double("Non Equilibrium BC Rate constant");
            non_equil_bc = param.get_bool("Use non equilibrium BC");
        }
        param.leave_subsection();


    }

private:



    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 [] );

    void gsub ( int &, double [], double & );

    static void gsub_wrapper ( int &, double [], double & );

    void dgsub ( int &, double [], double [] );

    static void dgsub_wrapper ( int &, double [], double [] );

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

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


    double cond_factor;

    double k_O2;

    bool non_equil_bc;

    // Transport Parameters
    double sigma_p;
    double lambda; //used for continuity
    bool uFD; //use finte differences 
    double endTol; //set tolerance for final solution
        double prev_effectiveness;

};

}
}

#endif