Anomalous dispersion in correlated porous media: a coupled continuous time random walk approach

Comolli, Alessandro; Dentz, Marco

doi:10.1140/epjb/e2017-80370-6

Cited by 30 publications

(30 citation statements)

References 87 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Meyer and Tchelepi (2010) have shown that the evolution of Lagrangian velocities may also be modeled as a Markov process in time whose evolution follows a nonlinear Langevin equation. The CTRW and TDRW frameworks model non-Fickian solute transport through broad distributions of independent advective transition times over given length scales intrinsic to the medium structure (Berkowitz & Scher, 1997;Comolli & Dentz, 2017;Painter & Cvetkovic, 2005). In this sense, these approaches assume that particle velocities form Markov processes when sampled equidistantly along trajectories (Benke & Painter, 2003;Cvetkovic et al, 1996;Fiori et al, 2006;Le Borgne et al, 2007, 2008a, this means at time points that are spaced by the advective transition time over a given distance (Cvetkovic et al, 1996;Le Borgne et al, 2007).…”

Section: 1029/2018wr023810mentioning

confidence: 99%

Upscaling and Prediction of Lagrangian Velocity Dynamics in Heterogeneous Porous Media

Hakoun

Comolli

Dentz

2019

Water Resources Research

Self Cite

View full text Add to dashboard Cite

The understanding of the dynamics of Lagrangian velocities is key for the understanding and upscaling of solute transport in heterogeneous porous media. The prediction of large‐scale particle motion in a stochastic framework implies identifying the relation between the Lagrangian velocity statistics and the statistical characteristics of the Eulerian flow field and the hydraulic medium properties. In this paper, we approach both challenges from numerical and theoretical points of view. Direct numerical simulations of Darcy‐scale flow and particle motion give detailed information on the evolution of the statistics of particle velocities both as a function of travel time and distance along streamlines. Both statistics evolve from a given initial distribution to different steady‐state distributions, which are related to the Eulerian velocity probability density function. Furthermore, we find that Lagrangian velocities measured isochronally as a function of travel time show intermittency dominated by low velocities, which is removed when measured equidistantly as a function of travel distance. This observation gives insight into the stochastic dynamics of the particle velocity series. As the equidistant particle velocities show a regular random pattern that fluctuates on a characteristic length scale, it is represented by two stationary Markov processes, which are parametrized by the distribution of flow velocities and a correlation distance. The velocity Markov models capture the evolution of the Lagrangian velocity statistics in terms of the Eulerian flow properties and a characteristics length scale and shed light on the role of the initial conditions and flow statistics on large‐scale particle motion.

show abstract

Section: 1029/2018wr023810mentioning

confidence: 99%

Upscaling and Prediction of Lagrangian Velocity Dynamics in Heterogeneous Porous Media

Hakoun

Comolli

Dentz

2019

Water Resources Research

Self Cite

View full text Add to dashboard Cite

show abstract

“…Under ergodic conditions, this means for a sufficiently large injection volume and flow domain, the steady space Lagrangian PDF

p_{ℓ} false(v false) = \lim_{ℓ \to \infty} {\overset{p}{true}}_{ℓ} false(v, ℓ false)

and the Eulerian velocity PDFs are related through flux‐weighting as shown in (Comolli & Dentz, ; Dentz et al, ; Kang et al, ):

p_{ℓ} false(v false) = \frac{v p_{e} false(v false)}{⟨ v_{e} ⟩} .

…”

Section: Flow and Transport Behaviormentioning

confidence: 99%

Linking Structural and Transport Properties in Three‐Dimensional Fracture Networks

et al. 2019

View full text Add to dashboard Cite

We investigate large‐scale particle motion and solute breakthrough in sparse three‐dimensional discrete fracture networks characterized by power law distributed fracture lengths. The three networks we consider have the same fracture intensity values but exhibit different percolation densities, geometric properties, and topological structures. We considered two different average transport models to predict solute breakthrough, a streamtube model and a Bernoulli continuous time random walk model, both of which provide insights into the flow fields within the networks. The streamtube model provides acceptable predictions at short distances in two of the networks but fails in all cases to predict breakthrough times at the outlet plane, which indicates that particle motion in such fracture networks cannot be characterized by a constant velocity between the inlet and control plane at which the breakthrough curve is detected. Rather, the structure of the network requires that frequent velocity transitions be made as particles move through the system. Despite the relatively broad distribution of fracture radii and relatively small number of independent velocity transitions, the continuous time random walk approach conditioned on the initial velocity distribution provides reasonable predictions for the breakthrough curves at different distances from the inlet. The application of these averaged transport models provides a richer understanding of the link from the fracture network structure to flow and transport properties.

show abstract

“…Continuous time random walk (CTRW) and time domain random walk models (Berkowitz et al, 2006;Noetinger et al, 2016) provide natural frameworks to upscale transport in media with spatially variable flow properties (Berkowitz & Scher, 1997;Comolli & Dentz, 2017;Dentz et al, 2016;Painter & Cvetkovic, 2005;Puyguiraud et al, 2019b). In these approaches, a solute plume is conceptualized as an assembly of idealized solute particles who transition through time and space by sampling the local flow velocities.…”

Section: Introductionmentioning

confidence: 99%

Characterizing the Influence of Fracture Density on Network Scale Transport

et al. 2020

Self Cite

View full text Add to dashboard Cite

The topology of natural fracture networks is inherently linked to the structure of the fluid velocity field and transport therein. Here we study the impact of network density on flow and transport behaviors. We stochastically generate fracture networks of varying density and simulate flow and transport with a discrete fracture network model, which fully resolves network topology at the fracture scale. We study conservative solute trajectories with Lagrangian particle tracking and find that as fracture density decreases, solute channelization to large local fractures increases, thereby reducing plume spreading. Furthermore, in sparse networks mean particle travel distance increases and local network features, such as velocity zones where flow is counter to the primary pressure gradient, become increasingly important for transport. As the network density increases, network statistics homogenize and such local features have a reduced impact. We quantify local topological influence on transport behavior with an effective tortuosity parameter, which measures the ratio of total advective distance to linear distance at the fracture scale; large tortuosity values are correlated to slow-velocity regions. These large tortuosity, slow-velocity regions delay downstream transport and enhance tailing on particle breakthrough curves. Finally, we predict transport with an upscaled, Bernoulli spatial Markov random walk model and parameterize local topological influences with a novel tortuosity parameter. Bernoulli model predictions improve when sampling from a tortuosity distribution, as opposed to a fixed value as has previously been done, suggesting that local network topological features must be carefully considered in upscaled modeling efforts of fracture network systems. and topology, and the corresponding flow field. The flow field within an individual fracture is typically highly correlated, commonly causing solute velocity to display persistent, low variability behavior over the in-fracture scale; consequently, the greatest Lagrangian accelerations occur at fracture intersections . As the fracture density increases, solute encounters more intersections on average, and the velocity correlation scale decreases. Furthermore, strong preferential flow paths form within interconnected networks of large fractures and channel a significant portion of mass, enabling solute to persist at high velocities for distances greater than the single fracture scale Sherman et al., 2019). This channelization becomes enhanced in sparse networks, where particles encounter fewer intersections, enabling them to persist on single fractures for longer distances. Resolving all these intranetwork features in 3-D DFN models is still computationally costly, and so upscaled modeling approaches, which account for network variability through effective parameter schemes, while maintaining a parsimonious framework, present an attractive alternative. However, how to properly parameterize network properties, such as velocity correlation and geometry, and inco...

show abstract

Anomalous dispersion in correlated porous media: a coupled continuous time random walk approach

Cited by 30 publications

References 87 publications

Upscaling and Prediction of Lagrangian Velocity Dynamics in Heterogeneous Porous Media

Upscaling and Prediction of Lagrangian Velocity Dynamics in Heterogeneous Porous Media

Linking Structural and Transport Properties in Three‐Dimensional Fracture Networks

Characterizing the Influence of Fracture Density on Network Scale Transport

Contact Info

Product

Resources

About