diff --git a/doc/handbook/dumux-handbook.bib b/doc/handbook/dumux-handbook.bib index 7455bca77d8db3bfa35d243999bb7ec87c64d069..4e8ca0dbc545855749061c91c7a2fd0c833749fc 100644 --- a/doc/handbook/dumux-handbook.bib +++ b/doc/handbook/dumux-handbook.bib @@ -1929,3 +1929,27 @@ year={1999} title = {Accelerating the Uzawa Algorithm}, journal = {{SIAM} Journal on Scientific Computing} } + +@article{Alanis2000, + author = {Elvio Alanis and Graciela Romero and Carlos C. Martinez}, + title = {{Interferometric measurement of diffusion coefficients through a scanning laser beam}}, + volume = {39}, + journal = {Optical Engineering}, + number = {3}, + publisher = {SPIE}, + pages = {744 -- 750}, + year = {2000}, + doi = {10.1117/1.602422}, + url = {https://doi.org/10.1117/1.602422} +} + +@article{Rard1979, + author = {Rard, Joseph A. and Miller, Donald G.}, + title = {{The mutual diffusion coefficients of Na2SO4−H2O and MgSO4−H2O at 25°C from Rayleigh interferometry}}, + volume = {8}, + journal = {Journal of Solution Chemistry}, + number = {10}, + year = {1979}, + doi = {10.1007/BF00648779}, + url = {https://doi.org/10.1007/BF00648779} +} diff --git a/dumux/material/fluidsystems/brine.hh b/dumux/material/fluidsystems/brine.hh index 98a55c712d7b01fec67d83023561cfa5a20c6140..7429f028da2b716165b75f59967b31896be08af5 100644 --- a/dumux/material/fluidsystems/brine.hh +++ b/dumux/material/fluidsystems/brine.hh @@ -440,6 +440,10 @@ public: * \param phaseIdx Index of the fluid phase * \param compIIdx Index of the component i * \param compJIdx Index of the component j + * + * The implemented value for NaCl is for a molar concentration of 2.5984 mol/l and a temperature of 25°C, see Rard and Miller, 1979 \cite Rard1979. + * Dependent on the salt concentration the coefficient can vary between 1.47e-9 m^2/s and 1.6e-9 m^2/s, see Rard and Miller, 1979. + * It also depends on temperature, values for different temparatures can e.g. found here: Alanis et al., 2000 \cite Alanis2000. */ template <class FluidState> static Scalar binaryDiffusionCoefficient(const FluidState& fluidState, @@ -454,11 +458,9 @@ public: using std::swap; swap(compIIdx, compJIdx); } - //! \todo TODO implement binary coefficients - // http://webserver.dmt.upm.es/~isidoro/dat1/Mass%20diffusivity%20data.pdf - // The link above was given as a reference in brine_air fluid system. + if (compJIdx == NaClIdx) - return 0.12e-9; + return 1.54e-9; else DUNE_THROW(Dune::NotImplemented, "Binary diffusion coefficient of components " << compIIdx << " and " << compJIdx