Abstract
Fractional-derivative models (FDMs) are promising tools for characterizing non-Fickian transport in natural geological media. Hydrologic applications of FDMs, however, have been limited in the last two decades, due to the lack of feasible models and solvers to quantify multi-dimensional anomalous diffusion for pollutants in bounded aquifers. This study develops and applies FDM tools to capture vector fractional dispersion for both conservative and reactive pollutants in fractional Brownian motion (fBm) random fields with bounded domains. A d-dimensional anisotropic fBm field for hydraulic conductivity (K) is first generated numerically. A particle-tracking based, fully Lagrangian solver is then developed to approximate particle dynamics in the fBm K fields under various boundary conditions, where the governing equation is the vector FDM subordination to regional flow. Numerical experiments show that the Lagrangian solver can combine nonlocal anomalous transport and local aquifer properties to quantify pollutant transport in bounded aquifers. Application analyses further reveal that the K correlation can significantly enhance the spreading of conservative pollutant particles, and increase the reaction rate by enhancing the mobility and mixing of reactant particles undergoing bimolecular reactions.
Extension of the Lagrangian solver is also discussed, including modeling transient flow, generalizing boundary conditions, and capturing complex chemical reactions. This study therefore provides the hydrologic community an efficient Lagrangian solver to model reactive anomalous transport in bounded anisotropic aquifers with any dimension, size, and boundary conditions.
This paper is dedicated to the memory of the late Professor Wen Chen
Appendix A. Lagrangian solution of the multi-scaling fADE (2.8)
The M-fADE (2.8) can define three types of vector diffusion, to capture complex multi-dimensional transport in geological media. The three types of vector diffusion exhibit different relationships between the mixing measure M and the scaling matrix H∗.
Type-1 vector diffusion: Eigenvectors in H∗ differ from the jumping angles in M(dθ). In some heterogeneous porous media, groundwater flow and contaminant transport do not need to follow the same directions. For example, water can flow along depositional stratigraphic (strike or dip) directions which can be defined by the scaling directions (i.e., eigenvectors) in H∗, while the contaminants move along preferential pathways formed by interconnected high-permeability deposits which can be defined by the angles dθ in the mixing measure M. This leads to model (2.8) with (usually orthogonal) eigenvectors in H∗ and an arbitrary mixing measure M(dθ); see [51] for examples.
Type-2 vector diffusion: Eigenvectors in H∗ are consistent with the jumping angles in M(dθ). For example, in discrete fracture networks (DFNs), water and contaminant move along the DFN’s orientations [38], which can be captured by model (2.8) with H∗ containing nonorthogonal eigenvectors and M(dθ) concentrated on eigenvectors. A specific example of type-2 vector diffusion is the classical 2-d fADE:
where Lt denotes the temporal operator, and the space indexes 1 < αx < 2 and 1 < αy < 2. Here the Riemann-Liouville fractional derivative is used for the space fractional derivative, to have correct solutions in bounded domains (while the Caputo fractional derivative leads to negative solutions at the boundary) [3]. Note that Professor W. Chen had also developed the Eulerian approaches to solve the vector fADE (6.1). For example, Chen et al. [9] developed the meshless Kansa method to solve anomalous diffusion in irregular flow fields, Chen and Wang [11] applied the singular boundary method to solve transient diffusion, and Chen et al. [12] also solved anomalous diffusion described by a structural derivative model. Hence, this study extends Professor Chen’s work using the Lagrangian solver.
Type-3 vector diffusion: Super-diffusion with a streamline-dependent mixing measure M. The directions of pollutant dispersion may not be fixed, but rather change in space due to the local variation of preferential pathways (note this assumption is widely used in hydrological models [21]), which can be captured by model (2.8) with a streamline-dependent M.
The Lagrangian solver for the M-fADE (2.8) is similar to that for the S-fADE (2.1), except for step 2 (random spatial jumps). Here the random displacement can be calculated using the following multi-scaling compound Poisson process for the i-th jump [51]:
where Z⃗(τ) is the random walker’s location at the operational time τ, [τ/dM] is the number of jumps by time τ with the step size dM, and Ri is the random length of the i-th jump along the independent jumping direction θi:
where R1/αk denotes the random displacement along the k-th eigenvector of H∗; Θ = ∂D/∂xk, –∂D/∂xk, or 0 if ∂D/∂xk > 0, < 0, and 0, respectively; and
Step 2 expressed by Eq. (6.2) and Eq. (6.3) can be modified for the three types of anomalous transport discussed above to achieve the best computational efficiency in the Lagrangian solver. For type-1 transport, components of the jump vector described by (6.3) are dependent; while for type-2 transport, the jump along each eigenvector is independent to the others, resulting in a simplified Lagrangian solver. For type-3 transport, the random displacement (6.3) can be projected to the adjacent streamline.
Acknowledgements
The work was partially supported by the National Natural Science Foundation of China (under grants 41931292, 41330632, and 11972148) and Guangdong Provincial Key Laboratory of Soil and Groundwater Pollution Control. This paper does not necessary reflect the view of the funding agencies.
References
[1] E.E. Adams, L.W. Gelhar, Field study of dispersion in a heterogeneous aquifer 2. spatial moment analysis. Water Resour. Res. 28, No 12 (1992), 3293–3307.10.1029/92WR01757Suche in Google Scholar
[2] B. Baeumer, Y. Zhang, R. Schumer, Incorporating super-diffusion due to sub-grid heterogeneity to capture non-Fickian transport. Ground Water53, No 5 (2015), 699–708.10.1111/gwat.12267Suche in Google Scholar PubMed
[3] B. Baeumer, M. Kovács, M.M. Meerschaert, H. Sankaranarayanan, Boundary conditions for fractional diffusion. J. Comput. Appl. Math. 336 (2018), 408–424.10.1016/j.cam.2017.12.053Suche in Google Scholar
[4] B. Baeumer, M. Kovács, H. Sankaranarayanan, Fractional partial differential equations with boundary conditions. J. Differ. Equ. 264, No 2 (2018), 1377–1410.10.1016/j.jde.2017.09.040Suche in Google Scholar
[5] D.A. Benson, The fractional advection-dispersion equation: Development and application. Ph.D. Dissertation, Univ. of Nev., Reno (1998).Suche in Google Scholar
[6] D.A. Benson, M.M. Meerschaert, B. Baeumer, H.P. Scheffler, Aquifer operator scaling and the effect on solute mixing and dispersion. Water Resour. Res. 42 (2006), W01415; 10.1029/2004WR003755.Suche in Google Scholar
[7] J.M. Boggs, S.C. Young, L.M. Beard, Field study of dispersion in a heterogeneous aquifer 1. Overview and site description. Water Resour. Res. 28, No 12 (1992), 3281–3291.10.1029/92WR01756Suche in Google Scholar
[8] A. Chechkin, V.Y. Gonchar, J. Klafter, R. Metzler, L.V. Tanatarov, Lévy flights in a steep potential well. J. Stat. Phys. 115, No 516 (2004), 1505–1535.10.1023/B:JOSS.0000028067.63365.04Suche in Google Scholar
[9] W. Chen, L. Ye, H.G. Sun, Fractional diffusion equations by the Kansa method. Comput. Math. Appl. 59, No. 5 (2010), 1614–1620.10.1016/j.camwa.2009.08.004Suche in Google Scholar
[10] W. Chen, G.F. Pang, A new definition of fractional Laplacian with application to modeling three-dimensional nonlocal heat conduction. J. Comput. Phys. 309, No 1 (2016); 10.1016/j.jcp.2016.01.003.Suche in Google Scholar
[11] W. Chen, F. Wang, Singular boundary method using time-dependent fundamental solution for transient diffusion problems. Eng. Anal. Bound. Elem. 68, No 7 (2016), 115–123.10.1016/j.enganabound.2016.04.004Suche in Google Scholar
[12] W. Chen, Y. Liang, X. Hei, Structural derivative based on inverse Mittag-Leffler function for modeling ultraslow diffusion. Fract. Calc. Appl. Anal. 19, No 5 (2016), 1250–1261; 10.1515/fca-2016-0064; https://www.degruyter.com/view/j/fca.2016.19.issue-5/issue-files/fca.2016.19.issue-5.xml.Suche in Google Scholar
[13] V.T. Chow, D.R. Maidment, L.W. Mays, Applied Hydrology. McGraw-Hill Publishing Companies (2013).Suche in Google Scholar
[14] O.A. Cirpka, A.J. Valocchi, Debates—Stochastic subsurface hydrology from theory to practice: Does stochastic subsurface hydrology help solving practical problems of contaminant hydrogeology? Water Resour. Res. 52 (2016); 10.1002/2016WR019087.Suche in Google Scholar
[15] C.J.M. Cremer, I. Neuweiler, M. Bechtold, J. Vanderborght, Solute transport in heterogeneous soil with time-dependent boundary conditions. Vadose Zone J. 15, No 6 (2016); 10.2136/vzj2015.11.0144.Suche in Google Scholar
[16] F. Delay, P. Ackerer, C. Danquigny, Simulating solute transport in porous or fractured formations using random walk particle tracking. Vadose Zone J. 4, No 2 (2005), 360–379.10.2136/vzj2004.0125Suche in Google Scholar
[17] M. Dentz, T. Le Borgne, A. Englert, B. Bijeljic, Mixing, spreading and reaction in heterogeneous media: A brief review. J. Contam. Hydrol. 120–121 (2011), 1–17.10.1016/j.jconhyd.2010.05.002Suche in Google Scholar PubMed
[18] G.E. Fogg, Groundwater flow and sand body interconnectedness in a thick multiple-aquifer system. Water Resour. Res. 22 (1986), 679–694.10.1029/WR022i005p00679Suche in Google Scholar
[19] G.E. Fogg, Y. Zhang, Debates-stochastic subsurface hydrology from theory to practice: A geologic perspective. Water Resour. Res. 53, No 12 (2016), 9235–9245.10.1002/2016WR019699Suche in Google Scholar
[20] R. Gorenflo, F. Mainardi, A. Vivoli, Continuous-time random walk and parametric subordination in fractional diffusion. Chaos Soliton. Fract. 34, No 1 (2007), 87–103.10.1016/j.chaos.2007.01.052Suche in Google Scholar
[21] C.T. Green, Y. Zhang, B.C. Jurgens, J.J. Starn, M.K. Landon, Accuracy of travel time distribution (TTD) models as affected by TTD complexity, observation errors, and model and tracer selection. Water Resour. Res. 50, No 7 (2014), 6191–6213.10.1002/2014WR015625Suche in Google Scholar
[22] I. Gupta, A.M. Wilson, B.J. Rostron, Groundwater age, brine migration, and large-sale solute transport in the Alberta Basin. Canada. Geofluids15 (2015), 608–620.10.1111/gfl.12131Suche in Google Scholar
[23] A.W. Harbaugh, MODFLOW-2005, The U.S. Geological Survey Modular Ground-Water Model - the Ground-Water Flow Process. U.S. Geological Survey Techniques and Methods 6-A16, Reston, Virginia (2005).10.3133/tm6A16Suche in Google Scholar
[24] K. Kang, S. Redner, Fluctuation-dominated kinetics in diffusion-controlled reactions. Phys. Rev. E32, No 7 (1985), 435–447.10.1103/PhysRevA.32.435Suche in Google Scholar
[25] M. Karamouz, R. Kerachian, B. Zahraie, Monthly water resources and irrigation planning: case study of conjunctive use of surface and groundwater resources. J. Irrig. Drain Eng. 130, No 5 (2004), 391–402.10.1061/(ASCE)0733-9437(2004)130:5(391)Suche in Google Scholar
[26] J.F. Kelly, H. Sankaranarayanan, M.M. Meerschaert, Boundary conditions for two-sided fractional diffusion. J. Comput. Phy. 376 (2019), 1089–1107.10.1016/j.jcp.2018.10.010Suche in Google Scholar
[27] J.C. Koch, R.C. Toohey, D.M. Reeves, Tracer-based evidence of heterogeneity in subsurface flow and storage within a boreal hillslope. Hydrol. Porcess. 31, No. 13 (2017), 2453–2463.10.1002/hyp.11205Suche in Google Scholar
[28] E.M. LaBolle, G.E. Fogg, A.F.B. Tompson, Random-walk simulation of transport in heterogeneous porous media: Local mass-conservation problem and implementation methods. Water Resour. Res. 32, No 3 (1996), 583–593.10.1029/95WR03528Suche in Google Scholar
[29] E.M. LaBolle, J. Quastel, G.E. Fogg, J. Gravner, Diffusion processes in composite porous media and their numerical integration by random walks: Generalized stochastic differential equations with discontinuous coefficients. Water Resour. Res. 36 (2000), 651–662.10.1029/1999WR900224Suche in Google Scholar
[30] E.M. LaBolle, G.E. Fogg, Role of molecular diffusion in contaminant migration and recovery in an alluvial aquifer system. Transp. Porous Media42 (2001), 155–179.10.1007/978-94-017-1278-1_8Suche in Google Scholar
[31] E.M. LaBolle, RWHet: Random Walk Particle Model for Simulating Transport in Heterogeneous Permeable Media, Version 3.2, User’s Manual and Program Documentation. Univ. of Calif., Davis (2006).Suche in Google Scholar
[32] B.Q. Lu, Y. Zhang, H.G. Sun, C.M. Zheng, Lagrangian simulation of multi-step and rate-limited chemical reactions in multi-dimensional porous media. Water Sci. Eng. 11, No 2 (2018), 101–113.10.1016/j.wse.2018.07.006Suche in Google Scholar
[33] B. Mandelbrot, J.W. van Ness, Fractional Brownian motions, fractional noises and applications. SIAM Rev. 10, No 4 (1968), 422–437.10.1137/1010093Suche in Google Scholar
[34] M. Marseguerra, A. Zoia, Monte Carlo investigation of anomalous transport in presence of a discontinuity and of an advection field. Physica A377, No 2 (2007), 448–464.10.1016/j.physa.2006.11.083Suche in Google Scholar
[35] M.M. Meerschaert, D.A. Benson, B. Baeumer, B., Operator Lévy motion and multiscaling anomalous diffusion. Phys. Rev. E63 (2001), ID 021112; 10.1103/PhysRevE.63.021112.Suche in Google Scholar PubMed
[36] F. Obiri-Nyarko, S.J. Grajales-Mesa, G. Malina, An overview of permeable reactive barriers for in situ sustainable groundwater remediation. Chemosphere111 (2014), 243–259.10.1016/j.chemosphere.2014.03.112Suche in Google Scholar PubMed
[37] L.J. Perez, J.J. Hidalgo, M. Dentz, Upscaling of mixing-limited bimolecular chemical reactions in Poiseuille flow. Water Resour. Res. 55, No 1 (2019), 249–269.10.1029/2018WR022730Suche in Google Scholar
[38] D.M. Reeves, D.A. Benson, M.M. Meerschaert, H.P. Scheffler, Transport of conservative solutes in simulated fracture networks: 2. Ensemble solute transport and the correspondence to operator-stable limit distribution. Water Resour. Res. 44 (2008), W05410; 10.1029/2008WR006858.Suche in Google Scholar
[39] K.R. Rehfeldt, J.M. Boggs, L.W. Gelhar, Field study of dispersion in a heterogeneous aquifer 3. Geostatistical analysis of hydraulic conductivity. Water Resour. Res. 28, No 12 (1992), 3309–3324.10.1029/92WR01758Suche in Google Scholar
[40] Salamon, P., D. Fernàndez-Garcia, J. J. Gómez-Hernández, Modeling mass transfer processes using random walk particle tracking. Water Resour. Res. 42 (2006), W11417; 10.1029/2006WR004927Suche in Google Scholar
[41] X. Sanchez-Vila, D. Fernàndez-Garcia, A. Guadagnini, Interpretation of column experiments of transport of solutes undergoing an irreversible bimolecular reaction using a continuum approximation. Water Resour. Res. 46 (2010), W12510; 10.1029/2010WR009539.Suche in Google Scholar
[42] W. Shao, Z. Yang, J. Ni, Y. Su, W. Nie, X. Ma, Comparison of single- and dual-permeability models in simulating the unsaturated hydro-mechanical behavior in a rainfall-triggered landslide. Landslides15, No 12 (2018), 2449–2464.10.1007/s10346-018-1059-0Suche in Google Scholar
[43] E.R. Siirila, R.M. Maxwell, Evaluating effective reaction rates of kinetically driven solutes in large-scale statistically anisotropic media: human health risk implications. Water Resour. Res. 48 (2012), W04527.10.1029/2011WR011516Suche in Google Scholar
[44] R. Sinha, M. Israil, D.C. Singhal, A hydrogeophysical model of the relationship between geoelectric and hydraulic parameters of anisotropic aquifers. Hydrogeol. J. 17 No 5 (2009), ID 495; 10.1007/s10040-008-0424-9.Suche in Google Scholar
[45] A.Y. Sun, R.W. Ritzi, D.W. Sims, Characterization and modeling of spatial variability in a complex alluvial aquifer: implications on solute transport. Water Resour. Res. 44 (2008), W04402.10.1029/2007WR006119Suche in Google Scholar
[46] H.G. Sun, A.L. Chang, Y. Zhang, W. Chen, A review on variable-order fractional differential equations: Mathematical foundations, physical models, numerical methods and applications. Frac. Calc. Appl. Anal. 22, No 1 (2019), 27–59; 10.1515/fca-2019-0003; https://www.degruyter.com/view/j/fca.2019.22.issue-1/issue-files/fca.2019.22.issue-1.xml.Suche in Google Scholar
[47] J.J.A. van Kooten, Groundwater contaminant transport including adsorption and first order decay. Stoch. Hydrol. Hydraul. 8, No 3 (1994), 185–205.10.1007/BF01587234Suche in Google Scholar
[48] E. Vidal-Henriquez, V. Zykov, E. Bodenschatz, A. Gholami, Convective instability and boundary driven oscillations in a reaction-diffusion-advection model. Chaos27 (2017), ID 103110; 10.1063/1.4986153.Suche in Google Scholar
[49] T. Vogel, H.H. Gerke, R. Zhang, M.Th. Van Genuchen, Modeling flow and transport in a two-dimensional dual-permeability system with spatially variable hydraulic properties. J. Hydro. 238, No 1-2 (2000), 78–89.10.1016/S0022-1694(00)00327-9Suche in Google Scholar
[50] Y. Wang, S.A. Bradford, J. Simunek, Estimation and upscaling of dual-permeability model parameters for the transport of E. coli D21g in soils with preferential flow. J. Contam. Hydro. 159 (2014), 57–66.10.1016/j.jconhyd.2014.01.009Suche in Google Scholar PubMed
[51] Y. Zhang, D.A. Benson, M.M. Meerschaert, E.M. LaBolle, H.P. Scheffler, Random walk approximation of fractional-order multiscaling anomalous diffusion. Phys. Rev. E74 (2006), ID 026706; 10.1103/PhysRevE.74.026706.Suche in Google Scholar PubMed
[52] Y. Zhang, E.M. LaBolle, K. Pohlmann, Monte Carlo approximation of anomalous diffusion in macroscopic heterogeneous media. Water Resour. Res. 45 (2009), W10417; 10.1029/2008WR007448.Suche in Google Scholar
[53] Y. Zhang, B. Baeumer, D.M. Reeves, A tempered multiscaling stable model to simulate transport in regional-scale fractured media. Geophy. Res. Lett. 37 (2010), L11405; 10.1029/2010GL043609.Suche in Google Scholar
[54] Y. Zhang, C. Papelis, Particle-tracking simulation of fractional diffusion-reaction processes. Phys. Rev. E84 (2011), 066704; 10.1103/PhysRevE.84.066704.Suche in Google Scholar PubMed
[55] Y. Zhang, C.T. Green, G.E. Fogg, Subdiffusive transport in alluvial settings: The influence of medium heterogeneity. Adv. Water Resour. 54 (2013), 78–99.10.1016/j.advwatres.2013.01.004Suche in Google Scholar
[56] Y. Zhang, M.M. Meerschaert, B. Baeumer, E.M. LaBolle, Modeling mixed retention and early arrivals in multidimensional heterogeneous media using an explicit Lagrangian scheme. Water Resour. Res. 51 (2015), 6311–6337; 10.1002/2015WR016902.Suche in Google Scholar
[57] Y. Zhang, C.T. Green, E.M. LaBolle, R.M. Neupauer, H.G. Sun, Bounded fractional diffusion in geological media: Definition and Lagrangian approximation. Water Resour. Res. 52, No 11 (2016), 8561–8577.10.1002/2016WR019178Suche in Google Scholar
[58] Y. Zhang, H.G. Sun, H.H. Stowell, M. Zayernouri, S.E. Hansen, A review of applications of fractional calculus in Earth system dynamics. Chaos Soliton. Frac. 102 (2017), 29–46.10.1016/j.chaos.2017.03.051Suche in Google Scholar
[59] Y. Zhang, B. Baeumer, L. Chen, D.M. Reeves, H.S. Sun, A fully subordinated linear flow model for hillslope subsurface stormflow. Water Resour. Res. 53 (2017), 3491–3504; 10.1002/2016WR020192.Suche in Google Scholar
[60] Y. Zhang, M.M. Meerschaert, Particle tracking solutions of vector fractional differential equations: A review. In: Handbook of Fractional Calculus with Applications, Volume 3: Numerical Methods (2019), 275–285.10.1515/9783110571684-010Suche in Google Scholar
[61] C. Zheng, M. Bianchi, S.M. Gorelick, Lessons learned from 25 years of research at the MADE site. Ground Water49, No 5 (2010), 649–662.10.1111/j.1745-6584.2010.00753.xSuche in Google Scholar PubMed
© 2019 Diogenes Co., Sofia
Artikel in diesem Heft
- Frontmatter
- Editorial Note
- FCAA special issue – In memory of late professor Wen Chen (FCAA–Volume 22–6–2019)
- Survey Paper
- State-of-art survey of fractional order modeling and estimation methods for lithium-ion batteries
- An investigation on continuous time random walk model for bedload transport
- Tutorial Survey
- Porous functions
- Research Paper
- A time-space Hausdorff derivative model for anomalous transport in porous media
- High-order algorithms for riesz derivative and their applications (IV)
- Mass-conserving tempered fractional diffusion in a bounded interval
- Dispersion analysis for wave equations with fractional Laplacian loss operators
- Lagrangian solver for vector fractional diffusion in bounded anisotropic aquifers: Development and application
- Some further results of the laplace transform for variable–order fractional difference equations
- Robust stability analysis of LTI systems with fractional degree generalized frequency variables
- Discovering a universal variable-order fractional model for turbulent Couette flow using a physics-informed neural network
Artikel in diesem Heft
- Frontmatter
- Editorial Note
- FCAA special issue – In memory of late professor Wen Chen (FCAA–Volume 22–6–2019)
- Survey Paper
- State-of-art survey of fractional order modeling and estimation methods for lithium-ion batteries
- An investigation on continuous time random walk model for bedload transport
- Tutorial Survey
- Porous functions
- Research Paper
- A time-space Hausdorff derivative model for anomalous transport in porous media
- High-order algorithms for riesz derivative and their applications (IV)
- Mass-conserving tempered fractional diffusion in a bounded interval
- Dispersion analysis for wave equations with fractional Laplacian loss operators
- Lagrangian solver for vector fractional diffusion in bounded anisotropic aquifers: Development and application
- Some further results of the laplace transform for variable–order fractional difference equations
- Robust stability analysis of LTI systems with fractional degree generalized frequency variables
- Discovering a universal variable-order fractional model for turbulent Couette flow using a physics-informed neural network