Constants.h

#ifndef __CONSTANTS_CONSTANTS_H__
#define __CONSTANTS_CONSTANTS_H__

#include <boost/math/constants/constants.hpp>
#include <cmath>
#include <limits>
#include <map>

#include "../Vector/Vector.hpp"
#include "./Helper.h"

namespace MDSimulation 
{
    namespace Constants
    {
        const double PHI = 1.6180339887;

        const double AVOGADRO = 6.02214129e23;

        const double GAS_CONSTANT = 8.3144621;
        
        const double BOLTZMANN_CONSTANT = 1.3806488e-23;

        const int DIMENSIONS = 3;

        const double NEVER = std::numeric_limits<double>::infinity();

        const double PI = boost::math::constants::pi<double>();

        const int NULL_PARTICLE = -1;

        const double P0 = 1.0e5;                   // atmospheric pressure

        const double KB = 1.380662e-23;            // Boltzmann's const

        const double GAMMA = 5.0 / 3.0;

        const double ExcludedVolume = 0.00005105;  // Van Der Waal's constant

        const double FusionBarrier = 4.5e7 * KB * 3.0 / 2.0;
        const double Plank = 6.626068e-34;

        const IntVector X(1, 0, 0);
        const IntVector Y(0, 1, 0);
        const IntVector Z(0, 0, 1);

        const IntVector STD_BASIS[3] =
        {
            X,
            Y,
            Z
        };

        const IntVector NEG_STD_BASIS[3] =
        {
            IntVector(-X),
            IntVector(-Y),
            IntVector(-Z)
        };

        const int LIQUID_TYPES = 5;

        const std::string LIQUID_NAMES[LIQUID_TYPES] =
        {
            "lithium",
            "mercury",
            "water",
            "diesel",
            "ig-4"
        };

        const std::map<std::string, int> LIQUID_FROM_NAME
            = invert_array_mapping<std::string>(LIQUID_NAMES, LIQUID_TYPES);

        const double Viscosity[LIQUID_TYPES] =
        {
            1.245e-6, // Lithium
            0.114e-6, // Mercury
            1.004e-6, // Water
            2.000e-6, // Diesel
            118.0e-6  // IG-4 from MultiTherm
        };

        const double SurfaceTension[LIQUID_TYPES] =
        {
            0.396,   // Lithium
            0.485,   // Mercury
            0.0728,  // Water
            23.8e-3, // Diesel
            23.8e-3  // IG-4 from MultiTherm (estimated from diesel)
        };

        const double Density[LIQUID_TYPES] =
        {
            516,    // Lithium
            5430,   // Mercury
            999.07, // Water
            832.00, // Diesel
            863.0   // IG-4 from MultiTherm
        };

        const double SpeedOfSound[LIQUID_TYPES] =
        {
            4490, // Lithium
            1450, // Mercury
            1497, // Water
            1250, // Diesel
            1300  // IG-4 from MultiTherm (estimate)
        };
    }
}

#endif /* __CONSTANTS_H__ */