/home/hauberg/Dokumenter/Capture/humim-tracker-0.1/src/OpenTissue/OpenTissue/dynamics/sph/sph_idealgas.h

Go to the documentation of this file.

#ifndef OPENTISSUE_DYNAMICS_SPH_SPH_IDEAL_GAS_H
#define OPENTISSUE_DYNAMICS_SPH_SPH_IDEAL_GAS_H
//
// OpenTissue Template Library
// - A generic toolbox for physics-based modeling and simulation.
// Copyright (C) 2008 Department of Computer Science, University of Copenhagen.
//
// OTTL is licensed under zlib: http://opensource.org/licenses/zlib-license.php
//
#include <OpenTissue/configuration.h>

#include <cassert>

namespace OpenTissue
{
  namespace sph
  {

    template<typename real_type>
    class IdealGas
    {
    public:

      real_type mol(real_type P,real_type V,real_type T) const
      {
        assert(V>0);
        assert(T>0);
        real_type R = static_cast<real_type>(8.314510);   
        return (P*V) / (R*T);
      }

      real_type mol(real_type V,real_type molar_volume) const
      {
        assert(V>0 && molar_volume>0);
        return V/molar_volume;
      }

      real_type molekyles(real_type P,real_type V,real_type T) const
      {
        assert(V>0);
        assert(T>0);
        real_type Na = static_cast<real_type>(6.0221367e23);  
        return mol(P,V,T)*Na;
      }

      real_type mass(real_type mol, real_type molar_mass) const
      {
        assert(mol>0 && molar_mass>0);
        return mol*molar_mass;
      }

      real_type celsius(real_type K) const
      {
        return K - static_cast<real_type>(273.15);
      }

      real_type kelvin(real_type C) const
      {
        return C + static_cast<real_type>(273.15);
      }

      real_type bar (real_type Pa) const
      {
        return Pa/static_cast<real_type>(1e5);
      }

      real_type Pa (real_type bar) const
      {
        return bar*static_cast<real_type>(1e5);
      }

      real_type kPa(real_type Pa) const
      {
        return Pa/static_cast<real_type>(1e3);
      }

      real_type molar_mass(real_type rho, real_type V, real_type mol) const
      {
        assert(V>0 && rho>0 && mol>0);
        return (rho*V)/mol;
      }

    };

  } // namespace sph
} // namespace OpenTissue

// OPENTISSUE_DYNAMICS_SPH_SPH_IDEAL_GAS_H
#endif