해양파와 해빙 상호작용을 위한 이산입자 모델 연구

초록

본 논문은 이산입자법(DEM)을 이용해 파동에 노출된 해빙의 기계적 거동을 모사한다. 구형 입자를 촘촘히 배열하고 접점에 인장·압축·모멘트를 모두 전달할 수 있는 결합을 부여해 파동에 의한 응력 분포와 미세 균열, 파편 크기 변화를 추적한다. 파고가 클수록 파단 결합 비율이 증가하고, 파편 평균 크기는 감소한다는 결과를 제시한다.

상세 분석

이 연구는 해양파와 해빙 사이의 복합적인 물리 현상을 정량적으로 해석하기 위해 이산입자법(Discrete Element Method, DEM)을 채택한 점이 가장 큰 특징이다. 기존의 연속체 기반 유한요소법(FEM)이나 스칼라 파동 이론은 파동이 입체적인 입자 구조에 미치는 미세한 응력 집중과 결합 파단 메커니즘을 포착하기 어려웠다. 본 논문은 입자를 구형으로 가정하고, 입자 간 접점에 ‘기계적 결합(bond)’을 설정함으로써 인장·압축·전단력뿐 아니라 모멘트까지 전달할 수 있게 설계하였다. 이러한 결합은 선형 탄성 거동을 따르며, 특정 임계 응력 혹은 모멘트에 도달하면 파괴되도록 구현되었다. 따라서 파동이 해빙 표면에 도달하면 입자 간 결합이 순차적으로 파단되며, 미세 균열이 전파되고 최종적으로 파편화가 일어나는 과정을 자연스럽게 재현한다.

모델 설정에서는 입자 직경, 결합 강도, 파동의 진폭·주기 등 주요 파라미터를 실험 데이터와 비교해 보정하였다. 파동은 1차원 사인파 형태로 입구에서 가해지며, 입자 군집 전체에 전파되는 과정에서 각 입자에 작용하는 힘과 모멘트를 실시간으로 기록한다. 결과 분석에서는 ‘파단 결합 비율(φ)’과 ‘평균 파편 크기(l)’를 핵심 지표로 삼았다. φ는 파동 진폭이 증가함에 따라 비선형적으로 상승했으며, 특히 진폭이 임계값을 초과할 때 급격히 증가하는 경향을 보였다. 반면 l은 진폭이 커질수록 감소했으며, 이는 파동 에너지가 결합 파단에 더 많이 사용되어 작은 파편이 다수 생성된다는 물리적 해석을 가능하게 한다.

이러한 결과는 해빙 파편화 메커니즘을 미시 수준에서 설명하는 데 기여한다. 특히, 파동 주기가 길어질수록 입자 간 결합에 전달되는 응력이 더 오래 지속되어 파단 가능성이 높아진다는 점이 시뮬레이션을 통해 확인되었다. 또한, 결합 파괴가 국소적으로 집중되는 영역과 전역적으로 퍼지는 양상을 시각화함으로써, 실제 해양 환경에서 발생할 수 있는 ‘파동‑해빙 상호작용에 의한 균열 전파 경로’를 예측하는 데 활용 가능하다.

하지만 모델에는 몇 가지 한계점도 존재한다. 첫째, 입자를 구형으로 단순화함으로써 실제 해빙의 비정형 결정 구조와 미세결정립의 이방성 효과를 완전히 반영하지 못한다. 둘째, 결합 파괴 기준을 단순 임계 응력/모멘트로 설정했기 때문에 온도·염도·염분 등에 의한 물성 변화가 고려되지 않았다. 셋째, 파동을 1차원 사인파로 가정했기 때문에 복합 파동 스펙트럼이나 비선형 파동 현상을 다루기 어렵다. 향후 연구에서는 비구형 입자, 온도 의존성 결합 모델, 다중 주파수 파동 입력 등을 도입해 모델의 현실성을 높일 필요가 있다.

전반적으로 본 논문은 DEM을 해양파‑해빙 상호작용 연구에 적용한 선구적인 시도이며, 파동 진폭·주기에 따른 결합 파단 비율과 파편 크기 변화를 정량적으로 제시함으로써 해빙 파편화 예측 모델링에 새로운 방향을 제시한다.