Journal
INTERNATIONAL JOURNAL OF HEAT AND MASS TRANSFER
Volume 133, Issue -, Pages 606-618Publisher
PERGAMON-ELSEVIER SCIENCE LTD
DOI: 10.1016/j.ijheatmasstransfer.2018.12.141
Keywords
Simulation; Fracture; Self-affine; Multiscale roughness; Non-Fickian transport; Mobile-immobile model
Categories
Funding
- National Natural Science Foundation of China [41602239, 41877171]
- Natural Science Foundation of Jiangsu Province [SBK20160861]
- Fundamental Research Funds of the Central Universities [2019816714]
- International Postdoctoral Exchange Fellowship Program [20150048]
Ask authors/readers for more resources
In this article, the influence of multiscale roughness on transport of a conservative solute through a self-affine fracture was investigated. The fracture roughness was decomposed into two different scales (i.e., a small-scale stationary secondary roughness superimposed on a large-scale non-stationary primary roughness) by a wavelet analysis technique. The fluid flow in the single fracture was characterized by Forchheimer's law and exhibited nonlinear flow features such as eddies and tortuous streamlines. The results indicated that the small-scale secondary roughness was primarily responsible for the nonlinear flow features. Numerical simulations of conservative solute transport showed non-Fickian transport characteristics (i.e., early arrivals and long tails) in breakthrough curves (BTCs) and in residence time distributions (RTDs) with and without consideration of the secondary roughness. Analysis of multiscale BTCs and RTDs showed that the small-scale secondary roughness played a significant role in enhancing the non-Fickian transport characteristics. Removing small-scale secondary roughness delayed the arrival time and shortened the tail. The peak concentrations in BTCs decreased as the secondary roughness was removed, implying that the secondary roughness could also enhance the solute dilution. Fitting the one-dimensional (1D) Fickian advection-dispersion equation (ADE) to the numerical BTCs resulted in considerable errors that decreased with the small-scale secondary roughness being removed. The 1D mobile-immobile model (MIM) provided a better fit to the numerical BTCs and inclusion of the small-scale secondary roughness in numerical simulations resulted in a decreasing MIM mobile domain fraction and an increasing mass exchange rate between immobile and mobile domains. (C) 2018 Elsevier Ltd. All rights reserved.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available