A transient depth-averaged lava flow model with a Herschel-Bulkley rheology accounting for three phases

기술 요약: 허셜-벌클리(Herschel–Bulkley) 유변학을 고려한 세 가지 상(Phase) 기반의 과도적 평균 수심 용암 흐름 모델

문제 제기
용암 흐름은 마그마가 냉각되고 결정화되며 가스가 탈기됨에 따라 발생하는 상당한 유변학적 변화에 의해 복잡한 역학을 나타낸다. 기존의 결정론적 모델들은 종종 용암을 단일 상 유체로 취급하거나 온도, 결정화 및 가스 함량 사이의 명시적인 결합을 무시함으로써 이러한 과정을 단순화한다. 천수(shallow-water, 평균 수심) 모델은 계산 속도와 물리적 완전성 사이의 균형을 제공하지만, 현재의 많은 구현 방식은 용암의 세 가지 상(용융물, 결정, 기포)의 특성을 명시적으로 고려하지 못하거나 결정화의 비평형 역학을 간과한다. 또한, 기포가 유변학에 미치는 영향은 흔히 지나치게 단순화되어 있으며, 점성 가열, 냉각 및 유속 사이의 피드백이 과도적 프레임워크 내에서 항상 완전히 통합되는 것도 아니다. 본 연구는 허셜-벌클리 유변학 하에서 용융물, 결정, 기포 간의 상호작용을 명시적으로 해결하는, 보다 물리적으로 완전하면서도 계산 가능한 1차원 평균 수심 모델의 필요성을 다룬다.

방법론
저자들은 세 가지 상의 현탁액에 대한 질량 및 운동량 보존 법칙에 기초한 과도적 1차원 평균 수심 모델을 개발하였다. 유변학은 항복 응력($\tau_c$), 일관성 계수(겉보기 점도 $\eta$), 그리고 유동 지수($n$)를 포함하는 허셜-벌클리 법칙에 의해 정의된다.

방법론의 핵심은 기존 문헌(구체적으로 Mueller et al., Truby et al., Mader et al.)의 종합을 통해 이러한 유변학적 매개변수에 대한 구성 관계를 도출하는 것이다:

1. 유변학적 폐쇄(Rheological Closures): 겉보기 점도는 입자 농도(Krieger-Dougherty 유형 관계 사용)와 기포 농도(모세관 수에 의존)의 효과를 결합한다. 항복 응력은 결정 농도의 함수로 모델링되며 최대 물리적 한계치를 갖는다. 유동 지수 $n$에 대한 새로운 관계식이 제안되었는데, 이는 입자의 종횡비와 기포 농도 모두에 의해 유도되는 전단 박화(shear-thinning) 거동을 설명하며, 흐름이 뉴턴 유체로 유지되는 임계 고체 부피 분율 임계값을 도입한다.
2. 두 가지 진화 폐쇄 방식: 본 연구는 결정 부피 분율($\Phi_p$)의 진화를 모델링하기 위해 두 가지 구별된 접근 방식을 제안하고 비교한다:
   • 모델 1 (이완 접근법): 결정 분율을 특성 결정 시간($\tau_{crys}$)을 갖는 평형 상태를 향한 이완 항의 영향을 받는 전달량으로 취급한다. 이 방식은 결정화를 명시적인 열 전달 메커즘으로부터 분리함으로써 복잡성을 최소화한다.
   • 모델 2 (열적 접근법): 결정 분율을 규정된 온도-결정도 관계를 통해 용암 온도와 결합한다. 온도의 진화는 복사, 대류, 전도 및 점성 가열의 네 가지 열 전달 메커니즘을 명시적으로 포함하는 평균 수심 내진 방정식에 의해 제어된다.
3. 수치적 구현: 유도된 시스템은 정수압 재구성(hydrostatic reconstruction)과 준암시적 이산화(semi-implicit discretization)를 사용하는 잘 균형 잡힌(well-balanced) 수치 스킴을 사용하여 이산화되며, 이는 강직한 소스 항(stiff source terms)과 비선형 유변학을 처리하기 위함이다.

주요 기여

• 삼상 유변학 프레임워크: 본 논문은 단일 상 근사를 넘어 결정과 기포 모두가 점도, 항복 응력 및 유동 지수에 미치는 영향을 명시적으로 통합하는 일반화된 유변학 모델을 제시한다.
• 이중 폐쇄 전략: 결정 진화를 위한 두 가지 구별된 폐쇄 방식을 도입하여, 단순화된 이완 모델과 복잡한 열역학 모델을 직접 비교할 수 있게 하였다. 이는 계산 효율성과 물리적 충실도 사이의 절충 관계를 강조한다.
• 과도적 역학: 이 프레임워크는 용암 흐름의 시간적, 공간적 변화를 포착하는 흐름 매개변수의 과도적 특성을 포착한다.
• 검증 및 보정: 유동 지수를 포함한 유변학적 매개변수들은 Truby et al. (2015)의 실험 데이터에 맞춰 보정되어, 삼상 현탁액에 대한 모델의 일관성을 확보하였다.

결과
모델들은 수치적으로 테스트되었으며 Mauna Ulu 분출의 채널형 용암 흐름 사례 연구에 적용되었다.

• 모델 간 비교: 실제 용암 흐름 데이터에 대한 모델 1과 모델 2의 비교는 온도의 다매개변수 진화(모델 2)가 이완 기반 폐쇄(모델 1)보다 흐름 매개변수에 대해 더 정확한 예측을 생성함을 시사한다.
• 과도적 거동: 모델은 불규칙하고 가파른 경사면을 따라 이동하는 제한된 용암이 엄격한 정상 상태(steady state)에 도달하지 않음을 성공적으로 예측한다. 대신, 흐름은 일련의 낙하 후 돌출(fill-then-breakout)되는 덩어리들로 특징지어지는 맥동적(pulsatory) 거동을 보인다. 이러한 맥동은 국부적 역학을 지배하며 시스템이 균일한 흐름으로 안정화되는 것을 방해한다.
• 기포의 역할: 본 연구는 탐구된 조건에서 기체 기포가 점도를 약 2배 정도 조절한다는 것을 발견했다. 기포는 전단 박화 거동을 유도하여 느린 흐름은 감속시키고 빠른 흐름은 가속시킨다. 기포를 경성 구체(hard spheres)로 취급한 단순화된 유변학을 사용했을 때는 벌크 점도만이 유일하게 영향을 받는 동적 매개변수였으며, 전단 박화의 미묘한 차이를 포착하지 못했다.
• 정상 해(Stationary Solutions): 모델 2에 대한 이론적 분석 결과, 점성 가열이 냉각을 상쇄할 수 있다면(이는 유속과 수심에 의존하는 균형에 달려 있음) 일정한 용암 플럭스를 유지할 수 있음이 밝혀졌다.

의의
본 논문은 세 가지 상의 유체 특성을 명시적으로 해결하고 열적 및 유변학적 진화 사이의 결합을 통해 용암 흐름 역학에 대한 더 강력한 물리적 설명을 제공한다는 점에서 그 주요한 의의를 주장한다. 흐름의 과도적 특성이 정상 상태보다는 맥동적 역학을 유발한다는 것을 입증함으로써, 본 모델은 불규칙한 지형에서의 흐름 거동을 더 현실적으로 표현한다. 또한, 폐쇄 방식의 비교는 단순화된 이완 모델이 계산 비용은 저렴하지만, 명시적인 열 전달 메커즘을 포함하는 모델 2가 용암 퇴적(emplacement)을 정확하게 예측하는 데 더 우수함을 보여준다. 이 연구는 확률론적 위험 모델과 복잡한 3D 결정론적 모델 사이의 간극을 메우며, 비평형 결정화 및 기포 유도 유변학적 변화와 같은 핵심적인 물리적 과정을 유지하는 평균 수심 도구를 제공한다.