경계 최적화 업윈드 SBP 연산자: 비등거리 격자와 고차 정확도

초록

본 논문은 경계 근처에 비등거리 격자를 배치하여 차수 2~9까지의 업윈드(상류) 유한차분 연산자를 설계한다. 대각 노름(SBP) 프레임워크 내에서 경계 폐쇄식을 최적화함으로써 전통적인 등거리 격자 기반 SBP 연산자에 비해 경계 정확도를 크게 향상시키고, 인공 점성 없이도 선형 안정성을 확보한다. SAT와 프로젝션 방식 두 가지 경계 조건 적용법을 모두 지원하며, 1차원 선형 하이퍼볼릭 문제와 2차원 압축성 유체 방정식에 대한 수치 실험을 통해 정확도·안정성·계산 효율성을 검증한다.

상세 분석

이 연구는 대각 노름을 갖는 SBP(합산-부분) 연산자를 업윈드 형태로 확장하고, 특히 경계 영역에서 비등거리 격자를 도입함으로써 기존 중앙 차분 기반 SBP 연산자의 경계 정확도 저하 문제를 해결한다. 전통적인 대각 노름 SBP 연산자는 등거리 격자에서 내부는 p차 정확도를 유지하지만, 경계에서는 최소 p/2(짝수 p) 혹은 (p‑1)/2(홀수 p) 차수만 보장한다. 이는 전체 수렴률을 제한하고, 고차 정확도를 요구하는 장기 시뮬레이션에서 비효율을 초래한다. 저자들은 정의 2.4·2.5에 따라 경계 오차 벡터 e(q)를 ℓ₂‑노름으로 측정하고, 이 오차의 선도항을 최소화하는 최적화 문제를 설정한다.

구성 절차는 다음과 같다. (1) 내부 스텐실은 기존 업윈드 SBP 연산자와 동일하게 p차 정확도를 갖는 중앙 좁은 스텐실을 기반으로, 앞·뒤 차분 연산자 Δ±와 적절한 가중 계수 αp를 이용해 한두 격자를 제거해 오프‑센터링을 만든다. (2) 경계 폐쇄식은 최소 s = p(짝수) 혹은 s = p‑1(홀수) 개의 경계 격점을 사용하고, 이 격점들의 위치 d₁, d₂를 자유 변수로 둔다. (3) 노름 행렬 H는 대각 원소 h₁…h_s 로 파라미터화하고, 양의 정부호 조건을 유지한다. (4) Q⁺ 행렬을 대칭·반대칭 제약(Q⁺+Q⁺ᵀ=0, Q⁺+Q⁻ᵀ=2S≤0)과 경계 정확도 조건(e(q)=0, q≤p/2 등)을 동시에 만족하도록 설계한다.

p>3인 경우 자유 파라미터가 남아 있어, (i) S가 음의 반정치(semi‑definite)임을 보장하고, (ii) ∥e(p/2+1)∥₂와 ∥e(p/2+2)∥₂의 가중합을 최소화하는 다목적 최적화를 수행한다. 이때 d₁, d₂는