🔬 condensed matter

Investigating frictional instability due to pressurization in granular media: insights from coupled computational fluid dynamics discrete element method

본 연구는 3D 결합 CFD-DEM 시뮬레이션을 활용하여, 입상 매질 내 유체 유도 불안정성이 단순히 간극 수압에 의해서가 아니라 배수 조건, 입자 골격(granular fabric)의 진화, 그리고 팽창성(dilatancy) 중심의 강도 증가에서 간극 수압 중심의 약화로의 전이를 포함하는 복잡한 수리-역학적 피드백에 의해 지배된다는 것을 입증한다.

원저자: Bimal Chhushyabaga (Department of Civil and Environmental Engineering, University of Houston), Behrooz Ferdowsi (Department of Civil and Environmental Engineering, University of Houston)

게시일 2026-06-05

📖 4 분 읽기☕ 가벼운 읽기

CC BY 4.0

원저자: Bimal Chhushyabaga (Department of Civil and Environmental Engineering, University of Houston), Behrooz Ferdowsi (Department of Civil and Environmental Engineering, University of Houston)

원본 논문은 CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) 라이선스로 제공됩니다. ✨ 이것은 아래 논문에 대한 AI 생성 설명입니다. 저자가 작성하거나 승인한 것이 아닙니다. 기술적 정확성을 위해서는 원본 논문을 참조하세요. 전체 면책 조항 읽기

모래 더미를 상상해 보세요. 모래 놀이터의 모래 같지만, 이 모래는 물에 흠뻑 젖어 있고 두 개의 무거운 판 사이에 꽉 끼여 있습니다. 이것이 물의 압력이 내부에서 높아질 때 젖은 알갱이 물질(토양, 모래, 또는 암석 가루와 같은)이 어떻게 갑자기 파괴되거나 미끄러지는지를 연구하기 위한 실험 설정입니다.

연구진은 초고성능 컴퓨터 시뮬레이션을 사용하여 이 현상이 입자 하나하나 단위로 어떻게 발생하는지 관찰했습니다. 이것은 모든 모래 알갱이와 물방울 하나하나를 개별적으로 추적하는 고해상도 영화와 같습니다.

연구 결과는 다음과 같으며, 쉬운 비유를 통해 설명합니다.

1. 설정: "샌드위치"와 "압력 밥솥"

빵이 무거운 판이고 속재료가 젖은 모래인 샌드위치를 상상해 보세요.

응력(Stress): 위쪽 판은 아래로 누르고 있고(중력), 옆으로 미끄러지려고 합니다(전단 응력).
트리거(Trigger): 그들은 샌드위치 속재료 내부에서 풍선을 부풀리는 것처럼 아래쪽에서 수압을 주입하기 시작했습니다.
목표: 그들은 어느 시점에서 수압이 모래의 결속력을 잃게 만들어 통제 불능의 미끄러짐을 유발하는지 알고 싶었습니다.

2. 두 가지 시나리오: "밀폐된 상자" vs "열린 창문"

연구진은 물이 어떻게 행동하는지에 따라 결과가 극적으로 달라지는 두 가지 다른 방식을 테스트했습니다.

"밀폐된 상자" (비배수 상태, Undrained): 샌드위치의 윗부분이 꽉 막혀 있다고 상상해 보세요. 물을 주입하면 압력이 갈 곳이 없습니다. 압력이 빠르게 쌓이면서 알갱이들을 밀어내고 알게 약간 떠 있게 만듭니다.
- 결과: 이는 갑작스럽고 격렬한 붕괴를 일으켰습니다. 모래는 순식간에 힘을 잃었고, 위쪽 판은 빠르게 미끄러졌습니다. 마치 갑작스러운 지진과 같았습니다.
"열린 창문" (배수 상태, Drained): 샌드위치 윗부분에 작은 구멍이나 물을 빼낼 수 있는 스펀지가 있다고 상상해 보세요. 물을 주입할 때 일부가 밖으로 빠져나갑니다.
- 결과: 압력이 높게 쌓이지 않았습니다. 모래는 한꺼번에 무너지지 않았습니다. 대신, 갑작스러운 추락보다는 언덕을 천천히 미끄러져 내려가는 것처럼 느리고 꾸준하게 약해졌습니다.

3. 비밀 재료: 단순히 수압 때문만이 아니다

수압만이 미끄러짐을 유발한다는 것은 흔한 오해입니다. 논문은 아니라고 말합니다.

모래 알갱이들을 원형으로 손을 잡고 있는 사람들의 팀("힘의 사슬", force chain)이라고 생각해 보세요.

팀워크: 사람들이 손을 꽉 잡고 있는 동안 팀은 강합니다.
물: 수압은 사람들을 서로 밀어내려 합니다.
실제 원인: 미끄러짐은 물이 강하게 밀 때 그리고 사람들이 이미 서로의 손을 놓았을 때만 발생합니다.
- 만약 팀이 이미 약하다면(모래가 이미 미끄러지기 직전이라면), 약간의 수압만으로도 충돌이 일 발생합니다.
- 만약 팀이 강하다면(모래가 미끄러지기에서 멀리 떨어져 있다면), 그 결속을 깨기 위해 많은 수압이 필요합니다.

연구진은 이 "손을 놓는 현상"이 알갱이들이 회전하고 재배열되면서 구조를 유지하는 강한 사슬을 끊기 때문에 발생한다는 것을 발견했습니다.

4. 알갱이들의 "교통 체증"

모래 내부에서 움직임은 균일하지 않습니다.

팽창과 압축(Dilation and Compaction): 어떤 부분의 모래는 사람들이 퍼져 나가려는 것처럼 팽창(dilate)하려고 했고, 다른 부분은 더 꽉 조여졌습니다(compacted).
투수성 변화(Permeability Shift): 이러한 재배열은 물이 모래를 통해 얼마나 쉽게 흐를 수 있는지를 변화시켰습니다. 이것은 마치 교통 체증이 갑자기 한 차선에서 풀리는 동시에 다른 차선은 막히는 것과 같습니다. 물은 새로운 경로를 찾아 이동했으며, 이는 압력이 가장 높은 곳을 변화시켰습니다.

5. "크리프(Creep)" vs "충돌(Crash)"

연구진은 모래가 얼마나 빨리 움직이는지 관찰했습니다.

느린 크리프(Slow Creep): 처음에는 모래가 매우 느리게, 거의 망설이는 것처럼 움직였습니다. 이 기간 동안 모래는 알갱이들이 더 단단하고 안정적인 위치를 찾아가면서(마치 사람들이 발 디딜 곳을 찾는 것처럼) 일시적으로 오히려 더 강해졌습니다.
충돌(The Crash): 결국 수압이 이 안정성을 압도했습니다. 알갱이들은 발 디딜 곳을 찾지 못하고 격렬하게 회전하기 시작했습니다. "힘의 사슬"(손잡기)이 무너졌고, 마찰력이 떨어지면서 빠른 미끄러짐이 일어났습니다.

6. "점성 수" (속도 제한)

연구진은 자신들의 모래를 걸쭉한 수프나 밀도가 높은 현탁액(예: 스무디)과 비교하려 했습니다. 그들은 모래가 유체처럼 행동하는지 보기 위해 "점성 수(viscous number)"라는 수학적 개념을 사용했습니다.

발견: 모래는 표준적인 "수프" 규칙에 딱 들어맞지 않았습니다. 모래는 너무 빽빽했고 알갱이들은 너무 고체적이었습니다. 그것은 매끄러운 액체의 흐름이라기보다, 자동차(알갱이)들이 가끔 서로 부딪히고 재배열되는 느린 교통 체증처럼 행동했습니다. 두꺼운 유체를 위한 표준 규칙들이 여기서 일어나는 일을 완벽하게 예측하지는 못했습니다.

요약

논문은 젖은 모래가 실패할 때, 단순히 수압이 밀어내는 것만이 아니라고 결론짓습니다. 그것은 복잡한 춤입니다:

수압은 알갱이들을 밀어내려 합니다.
알갱이들은 버티기 위해 스스로를 재배열하려고 노력합니다.
만약 물이 빠져나갈 수 없다면(밀폐된 상자), 압력이 승리하고 알갱이들은 한꺼번에 결속력을 잃습니다.
만약 물이 빠져나갈 수 있다면(열린 창문), 압력이 해소되어 실패가 천천히 일어납니다.
결정적으로, 실패는 알갱이들이 무게를 지탱할 수 없을 정도로 내부의 "손잡기"(접촉 네트워크)가 충분히 약화되었을 때만 발생합니다.

기술 요약: 과압에 의한 입상 매질의 마찰 불안정성 조사

문제 정의
유체 가압은 단층 구구(fault gouge), 강우 영향을 받는 사면, 주입으로 인해 교란된 저류층을 포함한 포화 지질 재료에서 변형과 불안정성을 유발하는 주요 트리거이다. 거시적 파괴 기준(예: Mohr–Coulomb)은 간극 수압에 의한 유효 수직 응력의 감소를 설명하지만, 압력 확산, 배수, 그리고 접촉 네트워크의 퇴화 사이의 연결 고리인 입자 규모의 피드백—즉, 미세한 압력 변화가 어떻게 동적 슬립(dynamic slip)으로 전이되는지를 결정하는 메커니즘—을 설명하는 데는 한계가 있다. 특히, 규정된 전단 응력 하에서 아임계 응력을 받는 입상 층이 어떻게 재활성화되는지, 그리고 배수 조건이 이 과정을 어떻게 조절하는지에 대한 메커니즘은 입자 수준에서 해결하기 어려운 과제로 남아 있다.

연구 방법론
본 연구는 구속된 유체 포화 입상 전단층 내에서 압력 유도 재활성화를 조사하기 위해 3차원 양방향 결합 전산유체역학-이산요소법(CFD–DEM) 시뮬레이션을 채택하였다.

수치 프레임워크: 유체 상(phase)은 체적 평균 질량 및 운동량 보존 방정식을 사용하여 고정된 오일러 격자(Eulerian grid) 상에서 계산되며, 개별 입자는 뉴턴의 법칙을 따르는 라그랑주 입자로 추적된다. 상 간 운동량 교환은 압력 구배 힘과 Syamlal–O'Brien 저해 항력(hindered drag)을 통해 모델링된다. 입자 간 상호작용에는 Coulomb 마찰을 포함한 연성 구체 Hertzian 접촉 법칙이 사용된다.
하중 프로토콜: 연구는 두 단계 접근 방식을 활용한다. 첫째, 서보 제어된 수직 응력 하에서의 변형률 제어 시뮬레이션을 통해 기준 마찰 강도와 Mohr–Coulomb 파괴 포락선을 설정한다. 둘째, 응력 제어 시뮬레이션에서는 기저 간극 수압을 점진적으로 증가시키면서 아임계 전단 응력( $\tau = 0.46\tau_{ss}$ ~ $0.90\tau_{ss}$ )을 가한다.
경계 조건: 시뮬레이션은 비배수(undrained) 조건(불투과성 상부 경계로 인해 과잉 간극 압력이 유지됨)과 배수(drained) 조건(상부에 대기압이 존재하여 수직 압력 소산이 허용됨)을 비교한다.
진단 지표: 분석은 거시적 지표(겉보기 마찰, 슬립 속도, 과잉 간극 압력)와 미시적 진단 지표(입자 회전, 힘 사슬(force-chain) 재편성, 배위수, 공극률 분포) 및 2차 구조 이방성 지표를 통합한다.

주요 결과

배수의 역할: 불안정성은 간극 압력 단독에 의해 결정되는 것이 아니라, 유효 응력 및 배수와의 결합된 진화에 의해 결정된다.
- 비배수: 과잉 간극 압력이 유지되어 급격한 유효 응력 감소, 급격한 마찰 약화 및 동적 슬립 가속화를 초래한다.
- 배수: 수직 압력 구배가 유지되어 과잉 압력 축적을 억제한다. 유체 탈출이 간극 압력 유도가 유효 응력 손실로 전환되는 것을 완충하기 때문에, 약화 현상이 더 부드럽고 점진적으로 나타난다.
응력 근접도와 구조 퇴화: 파괴를 유발하는 데 필요한 압력 증가는 초기 응력 상태와 전단 강도 사이의 관계에 따라 달라진다.
- 파괴에 가까운 층( $\tau \approx 0.9\tau_{ss}$ )은 완만한 압력 증가 후 파괴되며, 구조적 퇴화가 덜 광범위하게 일어난다.
- 파괴에서 먼 층( $\tau \approx 0.46\tau_{ss}$ )은 더 큰 압력 축적이 필요하며, 힘 사슬 붕괴 및 배위수 감소를 포함한 더 심각한 접촉 네트워크 퇴화를 겪는다.
미시적 진화: 파괴는 분산된 변형에서 국부화된 슬립으로의 전이로 특징지어진다.
- 힘 사슬(Force Chains): 거시적 약화는 강한 힘의 골격(strong-force backbone)의 붕괴 및 재편성과 일치한다. 파괴 후, 낮은 전단 응력 사례는 방향성 힘 사슬 조직의 상당한 상실을 보이는 반면, 높은 전단 응력 사례는 더 일관된 이방성 구조를 유지한다.
- 공극률 및 투과율: 파괴는 팽창(dilation)과 압축(compaction)의 교대 띠를 포함한다. 이러한 체적 변화는 공극률 유래 유효 투과율을 재편성하며, 이는 유압 경로가 변형 중에 동적으로 진화함을 나타낸다.
유변학적 스케일링: 마찰 응답은 낮은 점성수( $I_v$ ) 영역( $10^{-8} \lesssim I_v \lesssim 10^{-5}$ )에서 점성수( $I_v$ )와 부분적으로 조직화된다. 그러나 데이터는 고유한 국부 유변학 법칙에 수렴하지 않으며, 이는 단순한 국부 $\mu(I_v)$ 법칙보다는 간헐적인 재배열 및 비국부적 응력 전달이 지배하는 크리프(creeping) 영역의 특성을 반영한다.

의의 및 주장
본 논문은 배수 제어된 수리-역학적 피드백과 구조 퇴화가 어떻게 간극 압력 유도를 마찰 불안정성으로 전환하는지에 대한 입자 해상도의 통찰을 제공한다고 주장한다. 주요 기여는 다음과 같다:

불안정성으로의 전이가 압력 생성, 유압 확산, 그리고 변형 유도 공극 부피 변화 사이의 경쟁에 의해 조절됨을 입증함.
배수 조건이 파괴 경로를 근본적으로 변화시킨다는 점을 보여줌: 비배수 시스템은 압력 유지에 따른 급격한 약화를 선호하는 반면, 배수 시스템은 압력 구배를 통한 점진적 파괴를 보인다.
접촉 네트워크의 퇴화 정도가 초기 응력 상태에 의존함을 밝힘: 파괴에 대한 근접도가 입자 골격이 부과된 전단 하중을 더 이상 지탱할 수 없게 되기 전까지 필요한 구조 재편성의 범위를 결정한다.
거시적 약화가 단순히 벌크 간극 압력이나 접촉 손실의 함수가 아니라, 방향성 하중 지지 구조의 상실인 접촉 네트워크 퇴화의 시스템 규모 발현임을 명확히 함.

저자들은 이러한 결과가 단층 구구, 강우 유발 사면, 주입 교란 시스템에서 발생하는 압력 유도 불안정성을 해석하기 위한 메커니즘적 기초를 제공하며, 배수가 불량한 층이 광범위한 압력 유지 맥락 내에서 국부적 구조 퇴화가 발생할 때 급격한 약화에 취약하다는 점을 강조한다고 결론짓는다.