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
//
// ----------------------------------------------------------------------------

#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; }
      
      extern inline double K_SI() { return 1.3806488e-23; }
      
      extern inline double N_A() { return 6.022140857e23; }

      // --- 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;  }
      
      // --- GROUP 6 ---
      // NOTE: Constants taken from "Toward a Unified Theory of Isotropic Molecular Transport Phenomena" 
      //       by Kerkhof and Geboers, Table C2. If you wish to implement other \Omega^{*(l,s)} values
      //       see this paper for constants
      
      extern inline double A11_visc_R1() { return 1.340794; }

      extern inline double B11_visc_R1() { return 0.326244; }

      extern inline double C11_visc_R1() { return 1.546648; }

      extern inline double D11_visc_R1() { return 2.768179; }
      
      extern inline double A11_visc_R2() { return 1.066993; }

      extern inline double B11_visc_R2() { return 0.157384; }

      extern inline double C11_visc_R2() { return 0.424013; }

      extern inline double D11_visc_R2() { return 0.698873; }
      
      extern inline double A22_visc_R1() { return 26.425725; }

      extern inline double B22_visc_R1() { return 0.045563; }

      extern inline double C22_visc_R1() { return -25.232304; }

      extern inline double D22_visc_R1() { return 0.016075; }
      
      extern inline double A22_visc_R2() { return 1.151508; }

      extern inline double B22_visc_R2() { return 0.145812; }

      extern inline double C22_visc_R2() { return 0.437374; }

      extern inline double D22_visc_R2() { return 0.670219; }
      
}

#endif