[論文レビュー] 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.

研究の動機と目的

• 空間相関のある透水係数場を有する不均質な多孔質媒体における異常（非Fick的）輸送の物理的起源を特定すること。
• 構造的不均質性のみが異常輸送を説明できるのか、あるいは動的パスウェイ統計が不可欠であるかを特定すること。
• 特に、フラックス重み付き透水係数分布 {wK} と CTRW フレームワークにおける遷移時間分布 ψ(t) の間の定量的関係を確立すること。
• ψ(t) の累乗則指数 β について、{wK} パラメータを用いた解析的表現を導出し、流れ経路統計から異常輸送を予測可能とすること。

提案手法

• 異なる σ² を有する2次元、対数正規分布、空間相関のある透水係数（K）場における定常流とトレーサー粒子追跡のシミュレーション。
• ドメイン出口における粒子到着時刻からブレークスルー曲線（BTCs）を算出し、異常輸送行動を特徴づける。
• 切り捨てられた累乗則（TPL）分布 ψ(t) = n t₁⁻¹ exp(−t/t₂)/(1 + t/t₁)¹⁺ᵝ を用いた連続時間ランダムウォーク（CTRW）フレームワークを適用し、BTCs にフィットさせて β、t₁、t₂ を抽出する。
• 主要な流れ経路に沿った K 値の統計的分布を定量化するため、粒子通過回数重み付きヒストグラム {wK} を構築する。
• Darcy の法則（t ∝ 1/K）を用いて、{wK} の ln(K) 依存性を透過時間 t に変換し、{wK} と ψ(t) の直接的な同等性を確立する。
• ψ(t) の累乗則領域における {wK} の対数微分と TPL の対数微分を等置することで、累乗則指数 β = (μ − ln K)/σ²² の解析的表現を導出する。

実験結果

リサーチクエスチョン

• RQ1 不規則な多孔質媒体における異常輸送の物理的起源は、K フィールドの構造的不均質性を越えて何に起因するか？
• RQ2 動的パスウェイ統計、特に粒子パスに沿ったフラックス重み付き透水係数分布 {wK} は、CTRW モデルにおける遷移時間分布 ψ(t) とどのように関係するか？
• RQ3 切り捨てられた累乗則（TPL）ψ(t) の累乗則指数 β は、主な流れ経路の統計的特性から直接予測可能か？
• RQ4 好ましい経路に沿った低透水係数透過ゾーンの存在が、ブレークスルー曲線で観察される長尾的・非Fick的挙動をどれほど生成するか？
• RQ5 フラックス重み付き透水係数ヒストグラム {wK} のパラメータと異常輸送指数 β の間には、定量的・解析的関係が存在するか？

主な発見

• 異常輸送を特徴づける切り捨てられた累乗則分布 ψ(t) は、主な粒子パスに沿ったフラックス重み付きヒストグラム {wK} と定量的に同等である。
• ψ(t) の累乗則領域は、低透水係数セル（ln K ≲ −1）に正確に対応しており、これらの希少だが高頻度通過の低Kゾーンが長尾的輸送の主な要因であることを示している。
• 累乗則指数 β は解析的に β = (μ − ln K)/σ²² として導出され、{wK} の平均 μ と分散 σ² を用いて、パス統計と異常輸送挙動の直接的リンクを提供している。
• モデルは、異常輸送が構造的不均質性のみでは説明できないことを示しており、特に透水係数値のフラックス重み付き分布という、粒子パスの動的統計が ψ(t) の根本的決定要因である。
• 空間相関のある K フィールドおよび複雑な流れ場ダイナミクスが存在する状況でも、{wK} と ψ(t) の同等性は定量的に正確である。
• σ² = 5 の場合、導出された β ≈ 1.6 は BTCs のフィッティングから得られた値と一致し、{wK} → ψ(t) のマッピングの予測能力を確認している。