fcst_constants.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: fcst_constants.h
// - Description: This namespace contains all necessary FCST constants
// - Developers: Valentin N. Zingan,    University of Alberta
//               Marc Secanell Gallart, University of Alberta
// - Id: $Id: fcst_constants.h 2605 2014-08-15 03:36:44Z secanell $
//
// ----------------------------------------------------------------------------

#ifndef _FCST_CONSTANTS_H_
#define _FCST_CONSTANTS_H_

#include <deal.II/base/tensor.h>
#include <deal.II/base/symmetric_tensor.h>

const unsigned int dim = deal_II_dimension;

namespace Constants
{
      // --- GROUP 1 ---

      extern inline double R() { return 8.314462176; }

      extern inline double F() { return 9.648533992e4; }

      extern inline double Pi() { return 3.141592654; }

      extern inline double E0() { return 8.854187818e-12; }

      extern inline double K() { return 8.617332478e-5; }

      // --- GROUP 2 ---

      extern inline const dealii::Tensor<1,dim> gravity_acceleration()
      {
             dealii::Tensor<1,dim> result;
             result[dim-1] = -9.81;
             return result;
      }

      extern inline const dealii::SymmetricTensor<2,dim> unit_tensor()
      {
             dealii::SymmetricTensor<2,dim> result;

             for(unsigned int i = 0; i < dim; ++i)
                    for(unsigned int j = 0; j < dim; ++j)
                           if( i == j )
                                  result[i][j] = 1.0;

             return result;
      }

      // --- GROUP 3 ---

      extern inline double A_vk() { return 1.16145; }

      extern inline double B_vk() { return 0.14874; }

      extern inline double C_vk() { return 0.52487; }

      extern inline double D_vk() { return 0.77320; }

      extern inline double E_vk() { return 2.16178; }

      extern inline double F_vk() { return 2.43787; }

      // --- GROUP 4 ---

      extern inline double A_diff() { return 1.06036; }

      extern inline double B_diff() { return 0.15610; }

      extern inline double C_diff() { return 0.19300; }

      extern inline double D_diff() { return 0.47635; }

      extern inline double E_diff() { return 1.03587; }

      extern inline double F_diff() { return 1.52996; }

      extern inline double G_diff() { return 1.76474; }

      extern inline double H_diff() { return 3.89411; }

      // --- GROUP 5 ---

      extern inline double b_0() { return -2.1794;    }

      extern inline double b_1() { return 2.9530e-2;  }

      extern inline double b_2() { return -9.1837e-5; }

      extern inline double b_3() { return 1.4454e-7;  }
}

#endif