강우 유발 얕은 산사태 모델링을 위한 분자동역학 접근

초록

**
본 논문은 강우에 의한 수분 침투가 정적 마찰을 감소시켜 얕은 산사태를 유발한다는 가정 하에, 2차원 분자동역학(MD) 기반 알고리즘을 제시한다. 입자 간 상호작용은 Lennard‑Jones 2‑1 포텐셜로 정의하고, 마찰 조건은 Mohr‑Coulomb 파괴 기준을 적용한다. 시뮬레이션 결과는 실제 산사태와 유사한 가속 패턴과 에너지·시간 분포가 파워‑law 형태임을 보여주며, Fukuzono 방법을 통한 파괴 시점 예측 가능성을 제시한다.

**

상세 분석

**
이 연구는 기존 연속체 역학 모델이 갖는 공간 해상도와 복잡도 한계를 극복하기 위해 입자 기반 이산 요소법(DEM)과 유사한 분자동역학(MD) 프레임워크를 도입하였다. 2차원 평면에 N개의 구형 입자를 배치하고, 각 입자는 중력, 정/동 마찰, 점성력, 그리고 인접 입자와의 포텐셜 상호작용을 동시에 받는다. 정적 마찰 계수 μₛ는 시간에 따라 감소하는 함수로 설정했으며, 이는 강우에 의한 토양 포화도가 증가함에 따라 마찰 저항이 약화되는 현상을 수치적으로 모사한다. 마찰 감소 임계값과 입자 속도 임계값 두 가지 조건이 동시에 만족될 때 슬라이딩이 시작되며, 이는 Mohr‑Coulomb 파괴 기준을 입자 수준에서 구현한 형태이다.

포텐셜은 실험 데이터가 부족한 상황을 고려해 Lennard‑Jones 2‑1 형태로 가정했는데, 이는 전통적인 12‑6 포텐셜보다 짧은 거리에서 강한 반발력을 제공하면서도 장거리에서 약한 인력을 유지한다. 이를 통해 입자 간 응집·분산 현상을 자연스럽게 재현한다. 점성 항은 유체 저항을 근사화한 것으로, 점성 계수 η를 조절함으로써 슬라이드 면의 댐핑 효과를 제어한다.

시뮬레이션 결과는 입자 집단이 임계 조건을 초과하면 급격한 가속을 보이며, 이후 속도는 점차 포화 상태에 도달한다. 에너지 방출 이벤트(국부적인 ‘Avalanche’)는 규모와 발생 간격이 Gutenberg‑Richter 법칙과 Omori 법칙에 부합하는 파워‑law 분포를 나타냈다. 이는 모델이 실제 지진‑유사 현상과 동일한 임계 현상을 포착함을 의미한다. 또한, Fukuzono(1985)의 역표면 변위 속도법을 적용해 시뮬레이션 데이터로부터 파괴 시점을 추정했을 때, 실제 파괴 시점과 높은 일치도를 보였다.

이러한 결과는 입자 기반 MD 모델이 복잡한 토양‑수분 상호작용, 마찰 감소 메커니즘, 그리고 비선형 동역학을 통합적으로 다룰 수 있음을 시사한다. 특히 파라미터 공간(μₛ 초기값, 감소율, η, 포텐셜 강도 등)을 광범위하게 탐색했을 때도 전반적인 거동 패턴이 일관되게 유지되는 점은 모델의 강인성을 강조한다. 다만 2차원 구현이라는 제한과, 실제 토양 입자의 비구형성·다공성 등을 고려하지 않은 점은 향후 3차원 확장 및 실험 검증이 필요함을 암시한다.

**