Hong-Guang Sun | USGS Science

USGS2018

Comparison of time nonlocal transport models for characterizing non-Fickian transport: From mathematical interpretation to laboratory application

Non-Fickian diffusion has been increasingly documented in hydrology and modeled by promising time nonlocal transport models. While previous studies showed that most of the time nonlocal models are identical with correlated parameters, fundamental challenges remain in real-world applications regarding model selection and parameter definition. This study compared three popular time nonlocal transport models, including the multi-rate mass transfer (MRMT) model, the continuous time random walk (CTRW) framework, and the tempered time fractional advection–dispersion equation (tt-fADE), by focusing on their physical interpretation and feasibility in capturing non-Fickian transport. Mathematical comparison showed that these models have both related parameters defining the memory function and other basic-transport parameters (i.e., velocity v and dispersion coefficient D ) with different hydrogeologic interpretations. Laboratory column transport experiments and field tracer tests were then conducted, providing data for model applicability evaluation. Laboratory and field experiments exhibited breakthrough curves with non-Fickian characteristics, which were better represented by the tt-fADE and CTRW models than the traditional advection–dispersion equation. The best-fit velocity and dispersion coefficient, however, differ significantly between the tt-fADE and CTRW. Fitting exercises further revealed that the observed late-time breakthrough curves were heavier than the MRMT solutions with no more than two mass-exchange rates and lighter than the MRMT solutions with power-law distributed mass-exchange rates. Therefore, the time nonlocal models, where some parameters are correlated and exchangeable and the others have different values, differ mainly in their quantification of pre-asymptotic transport dynamics. In all models tested above, the tt-fADE model is attractive, considering its small fitting error and the reasonable velocity close to the measured flow rate.

Water↗

USGS2016

Bounded fractional diffusion in geological media: Definition and Lagrangian approximation

Spatiotemporal Fractional-Derivative Models (FDMs) have been increasingly used to simulate non-Fickian diffusion, but methods have not been available to define boundary conditions for FDMs in bounded domains. This study defines boundary conditions and then develops a Lagrangian solver to approximate bounded, one-dimensional fractional diffusion. Both the zero-value and non-zero-value Dirichlet, Neumann, and mixed Robin boundary conditions are defined, where the sign of Riemann-Liouville fractional derivative (capturing non-zero-value spatial-nonlocal boundary conditions with directional super-diffusion) remains consistent with the sign of the fractional-diffusive flux term in the FDMs. New Lagrangian schemes are then proposed to track solute particles moving in bounded domains, where the solutions are checked against analytical or Eularian solutions available for simplified FDMs. Numerical experiments show that the particle-tracking algorithm for non-Fickian diffusion differs from Fickian diffusion in relocating the particle position around the reflective boundary, likely due to the non-local and non-symmetric fractional diffusion. For a non-zero-value Neumann or Robin boundary, a source cell with a reflective face can be applied to define the release rate of random-walking particles at the specified flux boundary. Mathematical definitions of physically meaningful nonlocal boundaries combined with bounded Lagrangian solvers in this study may provide the only viable techniques at present to quantify the impact of boundaries on anomalous diffusion, expanding the applicability of FDMs from infinite do mains to those with any size and boundary conditions.

Water Resources Research↗