[论文解读] Origins of Anomalous Transport in Disordered Media: Structural and Dynamic Controls

Yaniv Edery, H. Scher+2 more|arXiv (Cornell University)|May 10, 2013

Groundwater flow and contamination studies被引用 5

一句话总结

本研究将无序多孔介质中异常输运的起因归因于对数正态分布的水力传导率场中通量加权的粒子路径。通过在连续时间随机游走（CTRW）框架下，将传导率的粒子访问加权直方图 {wK} 与过渡时间 ψ(t) 的截断幂律分布直接关联，作者表明 ψ(t) 的幂律行为源于沿优选流路的低传导率传输时间的统计分布，并从 {wK} 参数中推导出幂律指数 β 的定量表达式。

ABSTRACT

We quantitatively identify the origin of anomalous transport in a representative model of a heterogeneous system---tracer migration in the complex flow patterns of a lognormally distributed hydraulic conductivity ($K$) field. The transport, determined by a particle tracking technique, is characterized by breakthrough curves; the ensemble averaged curves document anomalous transport in this system, which is entirely accounted for by a truncated power-law distribution of local transition times $ψ(t)$ within the framework of a continuous time random walk. Unique to this study is the linking of $ψ(t)$ directly to the system heterogeneity. We assess the statistics of the dominant preferred pathways by forming a particle-visitation weighted histogram $\{wK\}$. Converting the ln($K$) dependence of $\{wK\}$ into time yields the equivalence of $\{wK\}$ and $ψ(t)$, and shows the part of $\{wK\}$ that forms the power-law of $ψ(t)$, which is the origin of anomalous transport. We also derive an expression defining the power law exponent in terms of the $\{wK\}$ parameters. This equivalence is a remarkable result, particularly given the correlated $K$-field, the complexity of the flow field and the statistics of the particle transitions.

研究动机与目标

• 确定具有空间相关水力传导率场的异质多孔介质中异常（非Fickian）输运的物理起源。
• 确定结构异质性本身是否足以解释异常输运，还是动态路径统计是必不可少的。
• 在连续时间随机游走（CTRW）框架下，建立粒子路径统计——特别是通量加权传导率分布 {wK}——与过渡时间分布 ψ(t) 之间的定量联系。
• 从 {wK} 参数出发，推导出 ψ(t) 幂律指数 β 的解析表达式，从而实现从流路统计预测异常输运。

提出的方法

• 在二维、对数正态分布、具有不同 σ² 的空间相关水力传导率（K）场中，模拟稳态流和示踪粒子追踪。
• 通过计算域出口处粒子到达时间的突破曲线（BTCs），表征异常输运行为。
• 应用具有截断幂律（TPL）分布 ψ(t) = n t₁⁻¹ exp(−t/t₂)/(1 + t/t₁)¹⁺ᵝ 的连续时间随机游走（CTRW）框架，拟合 BTCs 并提取 β、t₁ 和 t₂。
• 构建粒子访问加权直方图 {wK}，以量化沿主导流路的 K 值统计分布。
• 通过达西定律（t ∝ 1/K）将 {wK} 的 ln(K) 依赖关系转换为传输时间 t，从而在 {wK} 与 ψ(t) 之间建立直接等价关系。
• 通过将 {wK} 的对数导数与 TPL 在幂律区域的对数导数相等，推导出幂律指数 β = (μ − ln K)/σ²² 的解析表达式。

实验结果

研究问题

• RQ1 在 K 场的结构异质性之外，无序多孔介质中异常输运的物理起源是什么？
• RQ2 动态路径统计——特别是粒子路径上传导率值的通量加权分布——与 CTRW 模型中的过渡时间分布 ψ(t) 有何关系？
• RQ3 是否能直接从主导流路的统计特性预测 TPL ψ(t) 的幂律指数 β？
• RQ4 沿优选路径的低传导率传输区的存在在多大程度上导致了突破曲线中观察到的长尾、非Fickian行为？
• RQ5 通量加权传导率直方图 {wK} 的参数与异常输运指数 β 之间是否存在定量、解析的关系？

主要发现

• 表征异常输运的截断幂律分布 ψ(t) 的过渡时间，与主导粒子路径上水力传导率的通量加权直方图 {wK} 在定量上等价。
• ψ(t) 的幂律区域精确对应于低传导率单元（ln K ≲ −1）且粒子通量高的范围，表明这些稀少但高度穿越的低 K 区域是长尾输运的主要来源。
• 幂律指数 β 被解析推导为 β = (μ − ln K)/σ²²，其中 μ 和 σ² 分别为 {wK} 分布的均值和方差，从而在路径统计与异常输运行为之间建立了直接联系。
• 模型表明，异常输运无法仅由结构异质性解释；粒子路径的动态统计——特别是 K 值的通量加权分布——是决定 ψ(t) 的根本因素。
• 即使在 K 场存在空间相关性及复杂流动场动力学的情况下，{wK} 与 ψ(t) 之间的等价关系在定量上仍完全成立。
• 当 σ² = 5 时，推导出的 β ≈ 1.6 与通过拟合 BTCs 获得的值一致，证实了 {wK} → ψ(t) 映射的预测能力。