Exploring Particle Clogging Mechanisms through Deterministic Pore Network Modeling

기술 요약: 결정론적 기공 네트워크 모델링을 통한 입자 폐쇄 메커니즘 탐구

문제 정의
다공성 매질 내 입자 폐쇄(Clogging)는 투과율과 유체 흐름을 변화시키는 중요한 과정으로, 지하수 정화, 여과 시스템 설계에서부터 지질 에너지 운영에 이르기까지 다양한 응용 분야에 영향을 미친다. 직접 수치 시뮬레이션(DNS) 및 격자 볼츠만 방법(LBM)은 고해상도의 통찰을 제공하지만, 대표 체적 요소(REV)를 나타내기에는 계산적으로 제한된 작은 영역에 국한된다. 기공 네트워크 모델링(PNM)은 중간 규모(meso-scale) 분석을 위한 계산 효율적인 대안을 제공하지만, 종종 기공 목(pore throat)이 언제 막히는지를 결정하기 위해 확률 기반의 기준에 의존한다. 이러한 확률적 접근 방식은 예측에 불확실성을 초래한다. 또한, 기존 모델들은 폐쇄가 시작된 이후 유체 플럭스(flux)와 입자 플럭스의 분리, 특히 부분적으로 막힌 영역에서의 유체 흐름과 입자 수송 사이의 역동적인 상호작용을 정확하게 포착하는 데 어려움을 겪는 경우가 많다.

방법론
본 연구는 Raoof 등(2013)이 개발한 종합적인 PNM 도구인 PoreFlow 프레임워크를 활용하여 입자 폐쇄 역학을 시뮬레이션한다. 방법론은 세 가지 핵심 구성 요소로 구축된다:

1. 결정론적 폐쇄 기준: 무작위 확률 분포에 의존하는 대신, 본 모델은 Wyss 등(2006)이 제안하고 Dersoir(2015)가 개선한 임계 입자 수 $N^*$에 기반한 결정론적 기준을 채택한다. $N^*$는 특정 기공 목을 막는 데 필요한 입자 수를 정의한다. 모델은 다음 요소들을 바탕으로 각 기공 목에 대해 $N^*$를 동적으로 계산한다:

   • 부착 잠재력(Attachment potential): DLVO 이론을 통해 추정되며 이온 강도(예: NaCl 농도)의 영향을 받는다.
   • 기하학적 요인: 구체적으로 기공 목 너비($W$)와 입자 직경($D$)의 비율을 따르며, 거듭제곱 법칙 관계를 가진다.
   • 유량: 유량이 임계값($Q_c$) 이상일 때 $N^*$에 영향을 미치는 문턱 효과가 포함된다.
   • 모델은 특정 기공 목을 통과하는 누적 입자 수가 $N^*$를 초과하면 해당 기공 목이 "폐쇄 상태(clogged state)"로 전이된다고 가정한다.

2. 유체-입자 디커플링(자기 여과): 폐쇄 발생 시, 모델은 유체 흐름과 입자 플럭스를 분리한다. 유체는 폐쇄된 기공 목을 통해 (비록 감소된 전도도로 인해) 계속 흐르지만, 입자 플럭스는 차단되어 하류의 유체 내 입자 농도를 사실상 zero로 만든다. 이는 기공 목 내에서 성장하는 "필터 케이크(filter cake)" 또는 플러그(plug)로 모델링된다. 이 플러그의 길이($L_{cl}$)는 잔류 입자의 부피를 기반으로 계산되며, 기공 목의 수리학적 전도도는 직렬 저항 모델(개방된 구간과 폐쇄된 구간의 결합)을 사용하여 업데이트된다.

3. 시뮬레이션 설정: $1000 \times 1000 \times 2000$ $\mu$m 도메인 내에 3D 기공 네트워크(1656개 기공, 3082개 기공 목)를 생성하였다. 시뮬레이션은 단상 유동 조건 하에서 일정한 입구 압력(2.3 kPa)과 일정한 입구 입자 농도를 적용하여 수행되었다. 본 연구는 입자 크기(1–5 $\mu$m)와 폐쇄 전도도 비율(폐쇄 상태와 개방 상태의 전도도 비)을 체계적으로 변화시켜 민감도 분석을 수행하였다.

주요 기여

• 결정론적 프레임워크: 본 연구는 확률적 모델의 불확실성을 줄이면서 물리적 해석 가능성을 유지하는 투명하고 모듈화된 결정론적 폐쇄 기준을 확립한다.
• 동적 전도도 업데이트: 폐쇄된 기공 목이 여전히 부분적인 투과성을 유지하도록 하여, 입자 수송이 중단되더라도 유체 흐름은 지속되는 "필터 케이크" 효과를 시뮬레이션할 수 있는 메커니즘을 구현한다.
• 민감도 분석: 입자 크기와 폐쇄된 기공 목의 전도도 비율이 폐쇄의 시간적 진화와 결과적인 네트워크 투과율 변화에 어떻게 영향을 미치는지에 대한 체계적인 조사를 제공한다.

결과

• 투과율 진화: 모델은 실험적 관찰과 일치하는 지수 함수적 투과율 감소를 성공적으로 포착한다. 폐쇄는 좁은 기공 목에서 시작되어 시간이 지남에 따라 넓은 기공 목으로 전파된다.
• 입자 크기 효과:
  • 큰 입자: 폐쇄 과정을 가속화하여 급격한 전도도 손실과 네트워크 연결성이 크게 저해되는 "완전 폐쇄" 상태를 유발한다. 이러한 시나리오는 입구 근처에 집중된 폐쇄(얕은 평균 폐쇄 깊이)를 초래하며, 흐름 재분배에 따라 속도 불균질성(변동 계수, $CV$)이 감소한다.
  • 작은 입자: 전도도가 점진적으로 감소하는 "부분 폐쇄" 상태를 유발한다. 이러한 시나리오는 더 오랜 기간 높은 네트워크 연결성을 유지하며, 우선적 흐름 경로(preferential flow paths)의 형성으로 인해 지속적인 속도 불균질성($CV$) 증가를 보인다.
• 전도도 비율의 영향: 높은 전도도 비율(폐쇄된 목이 여전히 상당한 흐름을 허용함을 의미)은 폐쇄 진행을 지연시킨다. 이 효과는 빠른 폐쇄 가능성(예: 큰 입자)이 있는 시나리오에서 더 두드러지게 나타나며, 이는 입자 크기와 폐쇄된 플러그의 투과성 사이의 임계한 상호작용을 시사한다.
• 공간적 분포: 폐쇄가 진행됨에 따라 평균 폐쇄 깊이는 감소하며(입구 쪽으로 이동), 폐쇄 깊이의 분산이 감소하여 폐쇄 이벤트가 입구 지점에 집중됨을 나타낸다.

의의 및 한계점
본 논문은 이 PNM 프레임워크가 지하수 정화, 여과 설계, 지질 에너지 분야에서 광범위하게 적용 가능한, 다공성 매질 내 폐쇄 거동을 예측하기 위한 강력한 도구임을 주장한다. 결정론적 접근 방식은 미세 규모 현상과 그것이 거시적 특성에 미치는 집단적 영향에 대한 해석 가능성을 높인다.

그러나 저자들은 다음과 같은 특정 한계점을 인정한다:

• 현재 모델은 완전한 폐쇄 이벤트가 발생하기 전에 기공 목 벽면에 입자가 부착되는 현상을 고려하지 않으므로, 총 입자 보유량을 과소평가할 가능성이 있다.
• 현재 모델은 단분산 입자와 단상 유동에 국한되어 있어, 이분산 입자 혼합물이나 다상 유동(예: 불포화대) 조건을 처리하는 능력이 부족하다.
• 본 연구는 향-실제 세계의 복잡성을 더 잘 반영하기 위해 입자 부착 역학과 다분산 시스템을 위한 확률적 시나리오를 포함하는 향후 개선 방향을 제시한다.