수치 구현에 따른 빙하 비침투 조건의 차이가 비선형 스토크스 흐름에 미치는 영향

초록

본 연구는 빙하·빙상 모델에서 바닥 슬라이딩을 정확히 표현하기 위해 비선형 스토크스 방정식에 비침투 조건을 적용하는 세 가지 수치 방법(강제 회전, 라그랑주 승수, 압력 기반 근사)을 비교한다. 인공 해법을 이용한 수렴 실험에서는 세 방법 모두 비슷한 수렴률을 보였지만, 실제 지형·고속 슬라이딩 상황에서는 압력 기반 근사법이 강제 회전법에 비해 1 %~5 % 정도의 속도 차이를 나타낸다. 안정성·정확도 측면에서 가장 널리 쓰이는 강제 회전(강한) 방법이 선호된다.

상세 분석

이 논문은 빙하 역학에서 핵심적인 “비침투(impenetrability)” 경계조건을 어떻게 수치적으로 구현하느냐에 따라 시뮬레이션 결과가 얼마나 달라지는지를 체계적으로 조사한다. 비침투 조건은 베이스에서 흐름이 수직으로 침투하지 않도록, 즉 속도 벡터가 베드 표면에 접선 방향만을 갖도록 강제한다. 비선형 스토크스 방정식(전단 점성도가 변하는 비뉴턴 흐름)과 Navier‑slip 경계조건을 동시에 다루는 상황은 수치적으로 매우 까다롭다.

세 가지 구현 방법을 구체적으로 살펴보면, 첫 번째인 강제 회전(strong) 방법은 각 요소(또는 면)마다 로컬 좌표계를 회전시켜 법선 방향을 하나의 축으로 만든 뒤, 그 축에 대한 속도 성분을 0으로 고정한다. 이는 직접적인 강제이므로 수치적으로 가장 정확하지만, 좌표 변환 과정에서 추가적인 연산 비용과 복잡성이 발생한다.

두 번째인 라그랑주 승수(weak) 방법은 베드 면마다 평균 흐름을 0으로 만드는 제약을 라그랑주 승수(λ)로 도입한다. 약한 형태이므로 전역 행렬에 추가적인 자유도를 삽입하게 되며, 3차원 문제에서는 행렬 조건수가 악화돼 수렴이 느려질 위험이 있다. 특히, λ가 물리적 의미를 갖는 압력과 혼동될 수 있어 해석적 해석이 복잡해진다.

세 번째인 압력 기반 근사(approximative) 방법은 기존의 압력 변수(p)를 이용해 비침투 조건을 동시에 만족시키려는 시도다. 압력은 원래 불압축성(∇·u=0) 제약을 위한 라그랑주 승수이지만, 여기서는 베드 법선 응력 σ·n≈p와의 근사 관계를 가정한다. 분석 결과, 이 방법은 불압축성 제약을 완화시키면서 비침투를 “근사적으로” 만족한다는 것이 밝혀졌다. 즉, 압력이 법선 응력에 충분히 근접할 경우에만 정확도가 유지된다.

수치 실험은 두 단계로 진행되었다. 첫 번째는 **제조된 해법(method of manufactured solutions)**을 이용해 인위적인 해를 만들고, 각 방법의 수렴률을 비교했다. 여기서는 모든 방법이 이론적 2차 수렴을 보이며 차이가 거의 없었다. 이는 인위적인 상황이 복잡한 지형·변동 슬라이딩을 충분히 반영하지 못했기 때문이다.

두 번째는 실제와 유사한 2D·3D 테스트 케이스를 설정했다. 여기서는 베드 토폴로지가 급격히 변하고, 슬라이딩 계수 β가 크게 변동하는 영역을 포함한다. 결과는 압력 기반 근사법이 강제 회전법에 비해 평균 속도 오차가 1 %(2D)에서 5 % 이상(3D) 발생함을 보여준다. 특히, 고속 슬라이딩 구역에서 압력과 법선 응력 사이의 차이가 커지면서 오차가 증폭된다.

안정성 측면에서는 약한 라그랑주 승수 방법이 3D에서 행렬 조건수 악화와 수렴 실패를 보였으며, 강제 회전법은 가장 견고하게 동작했다. 따라서 현재 대부분의 빙상 모델이 채택하고 있는 강제 회전(강한) 방법이 정확도·안정성·계산 효율성 모두에서 최적이라고 결론짓는다.

이 연구는 빙하 모델링에서 경계조건 구현이 결과에 미치는 영향을 정량적으로 보여줌으로써, 모델 개발자들에게 구현 선택에 대한 명확한 가이드라인을 제공한다. 특히, 압력 기반 근사법을 사용할 경우 압력‑응력 관계를 사전에 검증하거나, 고속 슬라이딩 구역에서는 보완적인 제약을 추가해야 함을 시사한다.