VOLNA 전천후 쓰나미 시뮬레이션을 위한 차세대 수치 모델

초록

본 논문은 쓰나미의 발생, 전파, 해안 침수 전 과정을 통합적으로 모사할 수 있는 VOLNA 코드를 소개한다. 비정형 삼각형 메쉬 기반의 유한체적 스킴을 적용해 복잡한 지형에서도 높은 정확도를 확보했으며, 습·건 전이 영역을 물리적으로 처리하는 새로운 알고리즘을 제시한다. 분석 해와 실험 데이터를 통한 검증 결과, 기존 반경험적 방법에 비해 침수 예측 정확도가 크게 향상되었음을 확인하였다.

상세 분석

VOLNA 코드는 쓰나미 모델링에 필요한 전 단계—지진·해저 변위에 의한 파동 발생, 개방해역에서의 비선형 전파, 그리고 해안가에서의 급격한 런업·런다운—를 하나의 프레임워크 안에서 구현한다는 점에서 혁신적이다. 핵심 수치 기법은 비정형 삼각형 메쉬 위에 구축된 2차원 비보존형 슈뢰딩거 방정식(Shallow Water Equations, SWE)의 유한체적(FV) 스킴이다. 삼각형 메쉬는 복잡한 해안선, 수중 지형, 그리고 인공 구조물까지 자유롭게 재현할 수 있어, 기존 정사각형 격자 기반 모델이 갖는 지형 해상도 제한을 극복한다.

시간 적분은 2차 정확도를 유지하는 강인한 TVD(총변동 감소) 룰 기반의 러너-쿠타(RK) 방법을 사용한다. 특히, 경계조건 처리에서 고차 정확도를 유지하기 위해 반사·투과 경계와 개방 경계를 각각 특성 기반으로 구현하였다. 이러한 접근은 파동 에너지의 인위적 손실을 최소화하고, 장거리 전파 시 발생하는 수치 확산을 억제한다.

가장 두드러진 기술적 기여는 습·건 전이(wet/dry) 영역의 물리적 처리이다. 기존 다수의 쓰나미 모델은 물이 거의 없는 건조 구역을 단순히 고정된 최소 수심으로 대체하거나, 인위적인 마찰 항을 삽입하는 반경험적 방법을 사용한다. VOLNA는 ‘수위‑깊이’ 변수의 비선형 결합을 통해 건조 구역을 동적으로 활성화·비활성화하며, 물이 다시 채워지는 과정에서도 질량 보존과 운동량 보존을 동시에 만족한다. 이를 위해 ‘수위 재설정(Height Re‑initialization)’과 ‘마스크 함수(Mask Function)’를 결합한 알고리즘을 도입했으며, 수치적 안정성을 확보하기 위해 CFL 조건을 동적으로 조정한다.

검증 단계에서는 해석적 해(예: 선형 파동 이론에 기반한 해양 전파 해)와 실험실 파동 탱크 실험 데이터를 활용했다. 특히, 복합 지형을 포함한 실험 사례에서 VOLNA는 침수 최대 높이와 침수 거리에서 평균 오차가 5% 이하로, 기존 상용 모델 대비 현저히 낮은 오차를 기록했다. 또한, 실제 2004년 인도양 쓰나미와 2011년 동일본 대지진에 대한 재현 실험에서도 관측된 해일 높이와 도착 시간을 정확히 재현했으며, 이는 모델의 실용성을 강하게 뒷받침한다.

한계점으로는 2차원 수심-수위 모델에 기반하기 때문에 수직 구조(예: 급격한 수심 변화가 있는 해저 절벽)와 비선형 파동 붕괴 현상을 완전하게 포착하기 어렵다는 점을 언급한다. 향후 연구에서는 3차원 비보존형 방정식과 고차 정확도 스키마를 결합하거나, GPU 기반 병렬 연산을 도입해 실시간 운영 능력을 강화하는 방향을 제시한다.