system_management.h

Go to the documentation of this file.

// ----------------------------------------------------------------------------
//
// FCST: Fuel Cell Simulation Toolbox
//
// Copyright (C) 2006-2013 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: system_management.h
// - Description: This class manages systems of equations
// - Developers: Valentin N. Zingan,    University of Alberta
//                       Marc Secanell Gallart, University of Alberta
//                       Madhur Bhaiya, University of Alberta
// - Id: $Id: system_management.h 1449 2013-08-30 21:18:58Z madhur $
//
// ----------------------------------------------------------------------------

#ifndef _FCST_FUELCELL_SYSTEM_MANAGEMENT_H_
#define _FCST_FUELCELL_SYSTEM_MANAGEMENT_H_

#include "optimization_block_matrix_application.h"

using namespace dealii;
using namespace AppFrame;

typedef std::map<   std::string , std::map<std::string, DoFTools::Coupling>   > couplings_map;

enum VariableNames
{
  nothing = 0,

  oxygen_molar_fraction,
  water_molar_fraction,
  hydrogen_molar_fraction,
  nitrogen_molar_fraction,
  protonic_electrical_potential,
  electronic_electrical_potential,
  membrane_water_content,
  temperature_of_solid_phase,
  temperature_of_fluid_phase,
  temperature_of_REV,

  oxygen_concentration,
  water_concentration,
  hydrogen_concentration,
  nitrogen_concentration,
  proton_concentration,

  total_pressure
};

namespace FuelCell
{

class SystemManagement : public Subscriptor
{
public:



  SystemManagement();

  SystemManagement(MeshWorker::BlockInfo&        block_info,
                   Table<2, DoFTools::Coupling>& cell_couplings,
                   Table<2, DoFTools::Coupling>& flux_couplings);

 ~SystemManagement();

  void initialize(MeshWorker::BlockInfo&        rblock_info,
                  Table<2, DoFTools::Coupling>& rcell_couplings,
                  Table<2, DoFTools::Coupling>& rflux_couplings)
  {
    block_info     = &rblock_info;
    cell_couplings = &rcell_couplings;
    flux_couplings = &rflux_couplings;
  }

  void declare_parameters(ParameterHandler& param) const;

  void initialize(ParameterHandler& param);




  const bool solution_in_userlist(const std::string& name) const;

  const unsigned int solution_name_to_index(const std::string& name) const;

  const unsigned int equation_name_to_index(const std::string& name) const;

  const unsigned int matrix_block_index(const std::string& equation_name,
                                        const std::string& solution_name) const;

  void make_cell_couplings(const std::vector<couplings_map>& src);

  void make_flux_couplings(const std::vector<couplings_map>& src);




  const std::vector<std::string>& get_all_solution_names() const
  {
    return all_solution_names;
  }

  const std::vector<std::string>& get_all_equation_names() const
  {
    return all_equation_names;
  }

  const std::vector<std::string>& get_solution_names() const
  {
    return solution_names;
  }

  const std::vector<std::string>& get_equation_names() const
  {
    return equation_names;
  }

  const unsigned int& get_number_of_solution_names() const
  {
    return n_solution_names;
  }

  void print_system_info() const;




  DeclException2(VariableNotFoundInFCSTVariables,
                 std::string,
                 std::string,
                 << "A "  << arg1 << " with name \"" << arg2 << "\" is not stored in available FCST solution variables");

  DeclException2(EquationNotFoundInFCSTEquations,
                 std::string,
                 std::string,
                 << "An " << arg1 << " with name \"" << arg2 << "\" is not stored in available FCST equations");

  DeclException2(VariableNotFoundInUserVariables,
                 std::string,
                 std::string,
                 << "A "  << arg1 << " with name \"" << arg2 << "\" is not stored in user defined solution variables");

  DeclException2(EquationNotFoundInUserEquations,
                 std::string,
                 std::string,
                 << "An " << arg1 << " with name \"" << arg2 << "\" is not stored in user defined equations");

  DeclException5(SystemMatrixBlockDoesNotExist,
                 std::string,
                 unsigned int,
                 unsigned int,
                 std::string,
                 std::string,
                 << "A "  << arg1 << " (" << arg2 << "," << arg3 << ") " << "does not exist, because the variable \"" << arg4 << "\" is not coupled with the equation \"" << arg5 << "\"");


  // DATA //



  MeshWorker::BlockInfo* block_info;

  Table<2, DoFTools::Coupling>* cell_couplings;

  Table<2, DoFTools::Coupling>* flux_couplings;


protected:



  void set_all_solution_names();

  void set_all_equation_names();

  void check_solution_names() const;

  void check_equation_names() const;


  // DATA //



  std::vector<std::string> all_solution_names;

  std::vector<std::string> all_equation_names;

  std::vector<std::string> solution_names;

  std::vector<std::string> equation_names;

  unsigned int n_solution_names;


};

} // FuelCell

#endif