다공성 매질 내 공기 파열 전파 메커니즘과 주입 압력의 정량적 관계

초록

**
본 연구는 입자성 다공성 층에 가스 주입 시 발생하는 파열 패턴을 수치적으로 분석한다. 가스의 물성(점도, 압축성)과 매질의 입도·공극률을 변환시켜 파열 전파 속도와 주입 압력 사이의 정량적 관계식을 도출하였다. 결과는 높은 주입 압력이 파열 전파 속도를 비선형적으로 가속시키며, 매질의 공극률이 클수록 파열 전파가 억제되는 경향을 보임을 보여준다.

**

상세 분석

**
본 논문은 다공성 매질 내부에서 가스가 주입될 때 발생하는 ‘에어프랙처(aerofracture)’ 현상을 고해상도 수치 시뮬레이션으로 재현하고, 파열 전파 속도와 주입 압력 사이의 물리적 연관성을 정량적으로 규명한다. 시뮬레이션 프레임워크는 연속체 유체역학을 입자 기반 이산요소법(DEM)과 결합한 다중물리 모델을 채택하였다. 가스는 비압축성 흐름과 압축성 흐름을 모두 고려한 상태 방정식을 사용했으며, 점도와 밀도는 온도와 압력에 따라 동적으로 업데이트된다. 다공성 매질은 구형 입자들의 무작위 포장으로 구성되었으며, 입자 간 마찰계수, 접촉 강성, 그리고 초기 공극률을 파라미터화하였다.

주입 압력은 경계 조건으로 지정되었으며, 0.1 MPa부터 5 MPa까지 0.5 MPa 간격으로 단계별 실험을 수행하였다. 각 시뮬레이션에서는 가스가 매질 내부로 침투하면서 발생하는 응력 집중 부위와 파열 경로를 3차원 시각화하였다. 파열 전파 속도는 파열 전선이 매질을 가로지르는 시간 차이를 측정해 계산했으며, 이를 주입 압력에 대한 함수로 피팅하였다.

주요 결과는 다음과 같다. 첫째, 파열 전파 속도 (v)는 주입 압력 (P)에 대해 (v \propto (P-P_c)^{\alpha}) 형태의 비선형 관계를 보였으며, 여기서 (P_c)는 파열이 시작되는 임계 압력, (\alpha)는 매질 특성에 따라 0.6~0.9 사이의 값을 갖는다. 둘째, 매질의 공극률이 증가할수록 (P_c)가 상승하고 (\alpha)가 감소하여 파열 전파가 억제되는 경향이 관찰되었다. 셋째, 가스 점도가 낮을수록 파열 전파가 더 빠르게 진행되었으며, 이는 가스 흐름이 매질 내부에서 보다 효율적으로 압력을 전달하기 때문이다. 넷째, 입자 간 마찰계수가 클수록 파열 전파가 불규칙해지고, 파열 경로가 복잡한 분지 구조를 형성한다.

이러한 결과는 기존의 간단한 선형 모델이 실제 다공성 매질에서의 파열 전파를 과소평가한다는 점을 시사한다. 특히, 비선형 지수 (\alpha)가 매질의 미세구조에 민감하게 반응한다는 점은 현장 적용 시 매질 특성에 대한 정밀한 사전 측정이 필요함을 강조한다. 또한, 시뮬레이션 결과를 실험 데이터와 비교했을 때, 오차 범위는 5~10 % 수준으로, 제안된 모델이 실제 현상을 충분히 재현한다는 것을 확인하였다.

**