분진 눈사태 흐름의 수학적 모델링과 수치 시뮬레이션

초록

본 논문은 눈사태 중 입자와 공기가 혼합된 에어로졸 영역을 두 개의 가용성 유체로 모델링하고, Fick 법칙에 의한 질량 확산을 포함한 비압축성 연속 방정식과 운동량 방정식을 제시한다. 완전 암시적 유한체적 스킴으로 높은 밀도비를 가진 시스템을 안정적으로 계산하며, 강체 장애물과의 충돌 과정을 상세히 재현한다. 부피분율, 속도, 압력 필드를 통해 실제 눈사태 현상을 정량적으로 분석하고, 실무 적용 가능성을 논의한다.

상세 분석

이 연구는 눈사태를 두 상(고체 입자와 공기)으로 이루어진 가용성 유체 혼합물로 가정하고, 전통적인 다상 흐름 모델 대신 비압축성 연속 방정식과 운동량 방정식을 사용한다는 점에서 혁신적이다. 각각의 상은 부피분율 α와 밀도 ρ로 기술되며, 전체 혼합물의 질량 보존은 ∂(ρ)/∂t+∇·(ρ u)=0 형태로 표현된다. 여기서 ρ는 α₁ρ₁+α₂ρ₂ 로 정의되며, α₁+α₂=1이라는 제약을 갖는다. 입자와 공기 사이의 질량 교환은 Fick 법칙에 기반한 확산 항 D∇α 로 모델링되며, D는 확산계수로서 온도·압력·입자 크기 등에 따라 가변적으로 설정될 수 있다.

운동량 방정식은 혼합물 전체에 적용되며, 압력 p와 점성항 μ∇²u 를 포함한다. 특히, 밀도비가 10³ 이상인 경우에도 수치적 안정성을 유지하기 위해 압력-속도 결합을 투르크-프라스코프 예측‑보정 방식이 아닌 완전 암시적 방식으로 풀었다. 시간 적분은 1차·2차 정확도를 보장하는 백워드·오일러와 크랭크-니콜슨 혼합 스킴을 적용했으며, 비선형 연산자는 뉴턴-라프슨 반복으로 해결한다.

공간 이산화는 제어용적 기반의 유한체적법(FV)으로 수행되며, 각 셀의 면적·부피를 정확히 계산해 보존성을 확보한다. 대류항은 고차 정확도와 비진동성을 동시에 만족하는 TVD(총변동감소) 스키마를 사용하고, 확산항은 중앙차분으로 처리해 수치 확산을 최소화한다. 경계조건은 고정된 강체 장애물에 대해 무슬립(no‑slip)과 비투과성(impermeable) 조건을 적용했으며, 입구·출구는 자유압력 및 지정된 부피분율 프로파일을 부여한다.

수치 실험에서는 (1) 평탄한 경사면에서의 자유 전파, (2) 급격한 경사 전환 구역, (3) 강체 방패와의 충돌 시나리오를 다루었다. 결과는 부피분율이 급격히 변하는 전단면에서 고밀도 입자 구름이 형성되고, 압력 파동이 전파되며, 장애물 뒤쪽에 저압·저속 영역이 발생함을 보여준다. 특히, 기존 단일 상 모델에 비해 입자‑공기 혼합 효과가 정확히 포착되어, 실제 눈사태의 파괴력과 침투 깊이를 정량적으로 예측할 수 있다.

이 모델은 밀도비가 큰 다상 흐름을 다루는 전통적 방법(예: VOF, Level‑Set)보다 구현이 간단하고, 확산 항을 통해 미세 입자 간 혼합을 자연스럽게 반영한다는 장점이 있다. 또한 완전 암시적 스키마 덕분에 시간 스텝 제한이 크게 완화되어, 실시간 혹은 장기 시뮬레이션에 적합하다. 다만, 입자 크기 분포와 비탄성 충돌 효과를 별도 모델링하지 않은 점은 향후 연구 과제로 남는다.