Premium
A direct simulation demonstrating the role of spacial heterogeneity in determining anomalous diffusive transport
Author(s) -
Voller Vaughan R.
Publication year - 2015
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1002/2014wr016082
Subject(s) - fractal , sierpinski carpet , permeability (electromagnetism) , anomalous diffusion , power law , exponent , fractal dimension , sierpinski triangle , porous medium , porosity , diffusion process , infiltration (hvac) , mechanics , diffusion , geometry , statistical physics , geology , mathematics , physics , mathematical analysis , geotechnical engineering , chemistry , thermodynamics , computer science , biochemistry , knowledge management , innovation diffusion , statistics , linguistics , philosophy , membrane
Typically the spreading length‐scale in a diffusion process increases in time as ℓ ∼ t n , n = 1 2. In the presence of spatial heterogeneities, however, so‐called anomalous diffusion can occur, here the time exponent n ≠ 1 2. The objective of this paper is to present a numerical simulation that directly demonstrates the link between spacial heterogeneity and anomalous diffusion. The simulation is of the infiltration of a fluid into a horizontal unit area Hele‐Shaw cell in which the permeability at specified locations, laid out as a Sierpinski fractal carpet pattern, can differ from the nominal value K = 1. When there is no permeability contrast, the fluid infiltration has a diffusion‐like behavior, i.e., the filled plan‐form area increases in time as F = A t n , n = 1 2. When there is a permeability contrast, however, although the evolution of infiltration still follows a power law, anomalous behavior is observed; subdiffusive ( n < 1 2) when K < 1 and superdiffusive ( n > 1 2) when K > 1. These anomalous behaviors persist, even when the permeability contrast is only imposed over the largest sized fractal pattern element. But, if the pattern over which the permeability contrast is imposed has a subdomain length scale (obtained by filling in the largest sized element in the Sierpinski carpet with the smaller sized elements), normal diffusion n = 1 2behavior is recovered. In the special case that the imposed permeability in the pattern elements is K ≪ 1 , an approximate model, directly relating the coefficient and exponent in the infiltration power law to the porosity and fractal dimension of the carpet pattern, is derived and validated.