인터페이스 점프 조건을 위한 보정함수 기반 고정밀 포아송 해법

초록

본 논문은 고정계수 포아송 방정식에서 인터페이스에 따른 불연속성을 처리하기 위해, 기존 고스트 플루이드 방식의 보정항을 함수 형태로 확장한 “보정함수 방법(CFM)”을 제안한다. 보정함수는 좁은 밴드 영역에서 정의되며, 자체가 포아송 형태의 PDE와 경계조건을 만족하도록 설계된다. 2차와 4차 정확도를 달성한 구현을 통해, 기존 방법보다 높은 차수의 정확도와 대칭성·밴드 구조를 유지한 선형 시스템을 얻을 수 있음을 보인다. 다양한 인터페이스 형상에 대한 수치 실험으로 방법의 견고함과 수렴성을 검증하였다.

상세 분석

이 연구는 인터페이스를 가로지르는 격자 스템플이 포아송 연산자를 직접 적용할 수 없게 되는 전통적인 문제점을 고스트 플루이드 메소드(GFM)의 아이디어를 확장함으로써 해결한다. GFM에서는 “보정항”을 미리 계산해 우변에 삽입하지만, 보정항을 정확히 구하는 것이 1차 혹은 2차 정확도에 머무는 한계가 있다. 저자들은 보정항을 점이 아닌 함수, 즉 보정함수 D(x)로 정의하고, 이 함수가 인터페이스 주변 좁은 영역 ΩΓ에서 자체적인 포아송 방정식 ∇²D = f⁺−f⁻와 인터페이스 상의 점프 조건 D=a, ∂ₙD=b를 만족하도록 한다. 이렇게 하면 보정함수는 해석적으로는 Cauchy 문제이지만, 실제 수치 해석에서는 고주파 성분이 차단된 데이터와 제한된 거리 내에서만 해를 구하므로 ill‑posed 특성이 억제된다.

보정함수의 수치 구현은 두 단계로 구성된다. 첫째, 인터페이스를 레벨셋 함수로 암시적으로 표현하고, Gradient‑Augmented Level Set(GA‑LS) 기법을 이용해 고정밀 곡률·법선 정보를 얻는다. 둘째, ΩΓ 내부에서 D(x)를 bicubic 보간 기반의 근사함수로 전개하고, 최소제곱법을 적용해 PDE와 경계조건을 만족하도록 계수를 결정한다. 이 접근법은 스템플이 인터페이스를 가로지르는 경우에도 동일한 9‑point 2차 중앙 차분 스키마를 그대로 사용할 수 있게 하며, 선형 시스템 자체는 기존 포아송 해법과 동일한 대칭·밴드 구조를 유지한다.

정밀도 측면에서 저자들은 2차와 4차 스킴을 각각 구현하고, 수치 실험을 통해 L₂와 최대오차에서 기대한 차수 수렴을 확인한다. 특히 4차 스킴은 bicubic 보간과 최소제곱 해석을 결합함으로써, 인터페이스 근처에서의 급격한 변화에도 불구하고 전역적으로 4차 정확도를 유지한다. 또한, 복잡한 곡선형 인터페이스, 다중 인터페이스, 그리고 인터페이스가 격자와 거의 일치하거나 거의 평행한 경우에도 안정적인 결과를 보여, 방법의 일반성과 견고함을 입증한다.

이와 같은 설계는 기존 Immersed Interface Method(IIM)이나 고차 GFM이 요구하는 복잡한 스템플 수정·비대칭 행렬을 피하면서도, 높은 차수의 정확도와 sharp interface capture를 동시에 달성한다는 점에서 큰 의미가 있다. 또한, 보정함수 PDE를 다양한 차수와 차원(3D)으로 확장할 수 있는 구조적 장점을 제공하므로, 향후 복합 물리 현상(예: 열‑전기‑유체 결합)이나 비정형 격자와의 결합에도 적용 가능할 것으로 기대된다.