이산 확률 모델로 보는 암반 표면 마모 형상

초록

본 논문은 입자 충돌에 기반한 이산 확률 모델을 제시하여, 하천 바닥의 암반 장애물 전면(스토스 면)의 마모 과정을 시뮬레이션한다. 평균 입자 크기, 충돌 방향의 평균값 및 충돌당 제거되는 부피라는 세 가지 제어 변수만으로도 실험실 및 현장 관측과 정량적으로 일치하는 정상 상태 형상이 도출됨을 보인다.

상세 분석

이 연구는 기존의 연속적인 편미분 방정식(PDE) 형태인 Firey 모델과 그 일반화 모델을 이산화하여, 실제 하천에서 발생하는 입자-암반 충돌을 보다 현실적으로 재현한다. 모델은 2차원 다각형 형태의 암석 프로파일 P와 입자 G를 각각 정점과 변으로 표현하고, 충돌 사건을 두 가지 유형(A: 정점‑변 충돌, B: 변‑정점 충돌)으로 구분한다. 사건 발생 확률 p₀은 입자 크기 d와 프로파일 높이 h의 비율에 따라 결정되며, d≫h이면 p₀≈1(즉, Firey 모델에 수렴)이고, d≪h이면 p₀≈0(모래 블라스팅 한계)으로 설정한다.

각 충돌 사건에서는 로그정규 분포를 따르는 부피(면적) Δ가 무작위로 선택되어 프로파일에서 제거된다. 사건 A에서는 충돌 방향 α(평균값)와 그 표준편차 σ₁을 이용해 α가 두 인접 변의 법선 사이에 위치하도록 정점을 선택하고, 해당 정점을 절단해 새로운 정점을 생성한다. 사건 B에서는 변의 길이에 비례한 가중치로 변을 선택하고, 그 변을 평행하게 이동시켜 면적을 감소시킨다. 이러한 절차를 반복함으로써 프로파일은 점진적으로 마모되며, 입자 크기와 충돌 방향에 따라 서로 다른 형태의 정상 상태가 나타난다.

모델의 핵심 파라미터는 α와 p₀이며, σ₁, σ₂, Δ는 시간 스케일을 조정하는 역할만을 한다. 실제 실험에서는 입자 직경이 60–80 mm, 프로파일 높이가 80 mm인 경우 p₀≈0.4가 도출되었으며, α와 프로파일 평균 경사 β 사이에 강한 선형 상관관계가 존재함을 확인했다(α≈58.8°−0.346β−0.0033β²). 이 관계는 p₀ 변화에 크게 영향을 받지 않아 모델의 강인성을 보여준다.

시뮬레이션 결과는 대형 입자(p₀≈0.4)일 때 매끄러운 곡선형 프로파일이, 소형 입자(p₀≈0.03)일 때 뚜렷한 각과 긴 변을 가진 다각형 형태가 형성되는 것을 재현한다. 또한 시간에 따른 마모 속도와 형태 변화가 실험 데이터와 정량적으로 일치한다. 특히, 정상 상태에 도달하는 시간은 Δ값에 비례하여 조정 가능하며, 이는 실제 하천에서의 침식 속도를 추정하는 데 활용될 수 있다.

이 모델은 연속 PDE가 다루기 어려운 유한 입자 충돌과 비균일 충돌 방향을 자연스럽게 포함한다는 점에서 기존 모델보다 물리적 현실성을 크게 향상시킨다. 또한, 입자 크기와 충돌 각도의 직접적인 물리적 의미를 갖는 파라미터를 사용함으로써 현장 조사와 실험 설계에 바로 적용할 수 있는 장점을 제공한다.