diffusivityconstanttau.hh 2.54 KB
 Holger Class committed Oct 24, 2013 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 // -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- // vi: set et ts=4 sw=4 sts=4: /***************************************************************************** * See the file COPYING for full copying permissions. * * * * This program 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. * * * * This program 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 this program. If not, see . * *****************************************************************************/ /*! * \file * * \brief Relation for the saturation-dependent effective diffusion coefficient */ #ifndef DIFFUSIVITY_CONSTANT_TAU_HH #define DIFFUSIVITY_CONSTANT_TAU_HH  Bernd Flemisch committed Oct 25, 2013 27 #include  Thomas Fetzer committed Mar 23, 2015 28 #include  Holger Class committed Oct 24, 2013 29 30 31 32 33 34 35 36 37 38  namespace Dumux { /*! * \ingroup fluidmatrixinteractionslaws * * \brief Relation for the saturation-dependent effective diffusion coefficient * * The material law is: * \f[  Thomas Fetzer committed Mar 23, 2015 39 40  * D_\text{eff,pm} = \phi * S_w * \tau * D * \f]  Holger Class committed Oct 24, 2013 41  *  Thomas Fetzer committed Mar 23, 2015 42  * with a constant tau.  Holger Class committed Oct 24, 2013 43  */  Thomas Fetzer committed Mar 23, 2015 44 template  Holger Class committed Oct 24, 2013 45 46 47 48 class DiffusivityConstantTau { public: /*!  Thomas Fetzer committed Mar 23, 2015 49 50  * \brief Returns the effective diffusion coefficient \f$[m^2/s]\f$ based * on a constant tortuosity value  Holger Class committed Oct 24, 2013 51 52 53 54 55 56 57 58 59 60  * * \param porosity The porosity * \param saturation The saturation of the phase * \param diffCoeff The diffusion coefficient of the phase */ static Scalar effectiveDiffusivity(const Scalar porosity, const Scalar saturation, const Scalar diffCoeff) {  Christoph Grueninger committed May 16, 2014 61  Scalar tau = GET_RUNTIME_PARAM(TypeTag, Scalar, tau);  Holger Class committed Oct 24, 2013 62 63 64 65 66 67  return porosity * saturation * tau * diffCoeff; } }; } #endif