theta_scheme.h

Go to the documentation of this file.

//---------------------------------------------------------------------------
//
//    FCST: Fuel Cell Simulation Toolbox
//
//    Copyright (C) 2009-2015 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: theta_scheme.h
//    - Description: Time loop with theta scheme transient discretization.
//    - Developers: A. Kosakian, University of Alberta
//
//---------------------------------------------------------------------------

#ifndef __theta_scheme_h
#define __theta_scheme_h

#include <solvers/transient_base.h>

namespace FuelCell
{
namespace ApplicationCore
{
    class ThetaScheme : public TransientBase
    {
    public:

        ThetaScheme (ApplicationBase& app);

        virtual void declare_parameters (ParameterHandler& param);

        virtual void initialize (ParameterHandler& param);

        virtual void solve(FuelCell::ApplicationCore::FEVector& u,
                           const FuelCell::ApplicationCore::FEVectors& in_vectors);
        
        static const FuelCell::ApplicationCore::Event new_assembly;
        
    private:
        double theta;
        
        void solve_operations(FuelCell::ApplicationCore::FEVector& u, const FuelCell::ApplicationCore::FEVectors& in_vectors, bool print_solution);    
        
        /*
         * Variable that contains current time (i.e. time that has passed from the initial time layer till the current layer).
         */
        double current_time;
        
        float current_layer;
        
        /*
         * The order of the given time discretization scheme.
         */
        unsigned short int scheme_order;
    };
}
}

#endif