전단 유동에서 상분리를 설명하는 새로운 분리 플럭스 모델

초록

본 논문은 기존 2상 질량 흐름 모델에 ‘분리‑플럭스’를 도입하여, 경사면을 따라 이동하는 토사·물 혼합물에서 고체와 액체가 전방·측면으로 강하게 분리되는 메커니즘을 제시한다. 모델은 압력 구배, 체적분율, 전단률, 흐름 깊이 등 다양한 물리량을 포함하며, 시뮬레이션을 통해 전방의 고체‑풍부 급두와 측면의 레비 형성, 후방의 액체‑풍부 꼬리 흐름을 재현한다.

상세 분석

이 연구는 Pudasaini(2012)의 2상 질량 흐름 방정식을 기반으로, 고체와 액체 사이의 상호작용을 정량화하는 새로운 ‘분리‑플럭스(separation‑flux)’ 항을 추가한다. 기존 모델은 질량 보존식(α_s h, α_f h)과 모멘텀 보존식만을 포함했으며, 두 상이 동일한 속도장을 공유한다고 가정한다. 그러나 실제 토사·물 흐름에서는 고체 입자가 전단에 의해 전방·측면으로 이동하고, 액체는 중앙·후방에 남는 현상이 관찰된다. 이를 설명하기 위해 저자는 고체와 액체의 체적분율(α_s, α_f) 및 그 기울기(∇α)와 압력 구배(∇p) 사이의 비선형 결합을 도입하였다.

분리‑플럭스는 다음과 같은 물리적 요인을 포함한다. 첫째, 압력 구배가 고체‑액체 경계에서 불균형을 일으켜 고체가 높은 압력 영역(전방·측면)으로 이동하도록 한다. 둘째, 체적분율 기울기가 큰 곳에서는 고체 입자가 ‘키네틱 시버링(kinetic sieving)’ 메커니즘에 의해 상향 이동하고, 작은 입자와 액체는 하향 침투한다. 셋째, 전단률(γ̇)과 흐름 깊이(h)의 함수로 정의된 전단‑의존성 마찰계수는 고체‑풍부 구역에서 전단 저항을 증가시켜 전방 급두를 기계적으로 강하게 만든다. 넷째, 유체 점도(η_f)와 가상 질량계수(C) 등은 액체‑풍부 구역에서 저항을 완화시켜 흐름이 뒤쪽으로 늘어지게 한다.

수식적으로는 고체와 액체의 질량 보존식에 각각 추가적인 플럭스 항 J_s 와 J_f 을 삽입한다. J_s = −∇·(D_s ∇α_s) + Π_s ∇p와 같이 확산‑형태와 압력‑구배‑형태를 결합하고, J_f 은 질량 보존을 만족하도록 J_f = −J_s 로 정의한다. 여기서 D_s는 고체 입자의 확산계수이며, Π_s는 압력 구배에 대한 민감도를 나타낸다. 이러한 항들은 모멘텀 방정식의 소스항에 포함되어, 고체‑액체 간의 상대속도 (u_f−u_s) 에 대한 저항을 조절한다.

시뮬레이션 결과는 다음과 같은 핵심 현상을 재현한다. (1) 전방 급두가 고체‑밀도가 높아 기계적 강성을 갖는다. (2) 측면에 레비가 형성되어 고체 입자가 축적되고, 중앙·후방은 점성 액체가 주를 이룬다. (3) 흐름 진행에 따라 고체‑액체 비율이 시간·공간적으로 변하면서, 충격압이 기존 모델 대비 30‑50 % 차이로 크게 변동한다. 이는 재해 위험도 평가와 방지 설계에 중요한 영향을 미친다.

이 모델은 기존의 ‘균일 혼합’ 가정을 탈피하고, 실제 현장에서 관찰되는 복합적인 상분리 현상을 물리적으로 설명한다는 점에서 큰 의의를 가진다. 또한, 압력 구배·체적분율·전단률·점도·밀도비 등 다양한 파라미터를 명시적으로 포함함으로써, 실험·현장 데이터와의 정량적 비교가 가능하도록 설계되었다. 향후 연구에서는 입도분포, 온도·화학 반응, 비정상적인 지형 효과 등을 추가하여 모델의 범용성을 확대할 여지가 있다.