H2O_N2.hpp

Go to the documentation of this file.

// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
// vi: set et ts=4 sw=4 sts=4:
/*
  Copyright (C) 2009-2013 by Andreas Lauser
  Copyright (C) 2012 by Bernd Flemisch
  Copyright (C) 2010 by Felix Bode

  This file is part of the Open Porous Media project (OPM).

  OPM is free software: you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation, either version 2 of the License, or
  (at your option) any later version.

  OPM is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU General Public License for more details.

  You should have received a copy of the GNU General Public License
  along with OPM.  If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef OPM_BINARY_COEFF_H2O_N2_HPP
#define OPM_BINARY_COEFF_H2O_N2_HPP

#include "HenryIapws.hpp"
#include "FullerMethod.hpp"

#include <opm/material/components/N2.hpp>
#include <opm/material/components/H2O.hpp>

namespace Opm {
namespace BinaryCoeff {

class H2O_N2
{
public:
    template <class Evaluation>
    static Evaluation henry(const Evaluation& temperature)
    {
        const double E = 2388.8777;
        const double F = -14.9593;
        const double G = 42.0179;
        const double H = -29.4396;

        return henryIAPWS(E, F, G, H, temperature);
    }

    template <class Evaluation>
    static Evaluation gasDiffCoeff(const Evaluation& temperature, const Evaluation& pressure)
    {
        typedef Opm::H2O<double> H2O;
        typedef Opm::N2<double> N2;

        // atomic diffusion volumes
        const double SigmaNu[2] = { 13.1 /* H2O */,  18.5 /* N2 */ };
        // molar masses [g/mol]
        const double M[2] = { H2O::molarMass()*1e3, N2::molarMass()*1e3 };

        return fullerMethod(M, SigmaNu, temperature, pressure);
    }

    template <class Evaluation>
    static Evaluation liquidDiffCoeff(const Evaluation& temperature, const Evaluation& /*pressure*/)
    {
        const double Texp = 273.15 + 25; // [K]
        const double Dexp = 2.01e-9; // [m^2/s]

        return Dexp * temperature/Texp;
    }
};

} // namespace BinaryCoeff
} // namespace Opm

#endif