[논문 리뷰] Origins of Anomalous Transport in Disordered Media: Structural and Dynamic Controls

Yaniv Edery, H. Scher외 2명|arXiv (Cornell University)|2013. 05. 10.

Groundwater flow and contamination studies인용 수 5

한 줄 요약

이 연구는 이방성 투과성 분포를 가진 다공성 매체에서의 이례적 운반 현상의 기원을 로그정규분포된 수압전도도 필드 내에서 유량에 가중치를 부여한 입자 경로에서 기인한다고 규명한다. 연속 시간 랜덤워크(CTRW) 프레임워크에서 전이 시간 ψ(t)의 절단된 멱법칙 분포와 직접적으로 연결된 입자 방문 빈도에 가중치를 부여한 전도도 {wK}의 히스토그램을 통해, 저전도도를 통과하는 시간의 통계적 분포가 멱법칙 행동 ψ(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.

연구 동기 및 목표

• 이종 다공성 매체에서 공간적으로 상관관계가 있는 수압전도도 필드를 가진 비정상적 운반 현상의 물리적 기원을 규명하는 것.
• 구조적 이질성만으로도 이례적 운반 현상을 설명할 수 있는지 여부를 판단하는 것.
• 특히 유량에 가중치를 부여한 전도도 분포 {wK}와 CTRW 프레임워크 내 전이 시간 분포 ψ(t) 사이의 정량적 연관성을 설정하는 것.
• {wK} 매개변수에 기반하여 ψ(t)의 멱법칙 지수 β에 대한 분석적 표현을 유도하여 유량 경로 통계로부터 이례적 운반 현상을 예측할 수 있도록 하는 것.

제안 방법

• 다양한 σ²를 가진 2차원, 로그정규분포, 공간적으로 상관관계가 있는 수압전도도(K) 필드에서 정 steady-state 유량 및 트레이서 입자 추적 시뮬레이션을 수행한다.
• 도메인 출구에서 입자 도착 시간을 기반으로 도착 곡선(BTCs)을 계산하여 이례적 운반 행동를 특성화한다.
• 절단된 멱법칙(TPL) 분포 ψ(t) = n t₁⁻¹ exp(−t/t₂)/(1 + t/t₁)¹⁺ᵝ 를 사용한 연속 시간 랜덤워크(CTRW) 프레임워크를 적용하여 BTCs에 적합하고 β, t₁, t₂ 를 추출한다.
• 주요 유량 경로를 따라 전도도 값의 통계적 분포를 정량화하기 위해 입자 방문 횟수에 가중치를 부여한 히스토그램 {wK} 를 구성한다.
• Darcy의 법칙(t ∝ 1/K)을 통해 {wK} 의 ln(K) 의존성을 전이 시간 t 로 변환함으로써 {wK} 와 ψ(t) 사이의 직접적 등가성을 확립한다.
• CTRW의 멱법칙 영역에서 {wK} 와 TPL 의 로그 도함수를 동일시함으로써 멱법칙 지수 β = (μ − ln K)/σ²² 의 분석적 표현을 유도한다.

실험 결과

연구 질문

• RQ1 이방성 다공성 매체에서의 이례적 운반 현상의 물리적 기원은 K 필드의 구조적 이질성 외에 무엇인가?
• RQ2 특히 입자 경로를 따라 유량에 가중치를 부여한 전도도 값의 통계적 분포가 CTRW 모델 내 전이 시간 분포 ψ(t) 와 어떻게 관련되어 있는가?
• RQ3 절단된 멱법칙(TPL) ψ(t) 의 멱법칙 지수 β 가 주요 유량 경로의 통계적 특성으로 직접 예측될 수 있는가?
• RQ4 선호 경로를 따라 존재하는 저전도도 전이 영역이 도착 곡선에서 관찰되는 긴 尾部 비프ickian 행동을 어떻게 생성하는가?
• RQ5 유량에 가중치를 부여한 전도도 히스토그램 {wK} 의 매개변수와 이례적 운반 지수 β 사이에 정량적이고 분석적인 관계가 존재하는가?

주요 결과

• 이례적 운반을 특징짓는 절단된 멱법칙 분포 ψ(t) 는 주요 입자 경로를 따라 유량에 가중치를 부여한 전도도 히스토그램 {wK} 와 정량적으로 동일하다.
• ψ(t) 의 멱법칙 영역은 높은 입자 유량을 가지는 저전도도 세포( ln K ≲ −1 ) 의 범위와 정확히 일치하며, 이는 희귀하지만 매우 많이 통과되는 저K 영역이 긴 尾部 운반의 주요 원인임을 시사한다.
• 멱법칙 지수 β 는 β = (μ − ln K)/σ²² 로 분석적으로 유도되며, 여기서 μ 와 σ² 는 각각 {wK} 분포의 평균과 분산이다. 이는 유량 경로 통계와 이례적 운반 행동 사이의 직접적 연결 고리를 제공한다.
• 모델은 이례적 운반 현상이 구조적 이질성만으로 설명될 수 없음을 보여주며, 특히 전도도 값의 유량에 가중치를 부여한 분포가 ψ(t) 의 결정적 요소임을 규명한다.
• {wK} 와 ψ(t) 의 등가성은 K 필드의 공간적 상관관계와 복잡한 유량장 역학이 존재하는 상황에서도 정량적으로 정확하다.
• σ² = 5 인 경우 유도된 β ≈ 1.6 이 BTCs에 대한 피팅 결과와 일치하여 {wK} → ψ(t) 매핑의 예측 능력을 확인한다.