공동 격자 기반 Navier‑Stokes 유한체적 해법의 절단오차 정밀 추정법

초록

본 논문은 공동(colocated) 격자에서 비압축성 Navier‑Stokes 방정식을 유한체적(FV) 방식으로 풀 때 발생하는 절단오차(truncation error)를, 미세 격자의 해를 동일한 차분 스키마를 갖는 더 거친 격자에 제한(restrict)하고, 그 격자에서의 이산 연산자를 적용해 추정하는 방법을 제시한다. 핵심은 질량 플럭스의 압력 부분을 선형화된 모멘텀 방정식의 계수와 무관하게 만드는 새로운 모멘텀 보간 기법이며, 이를 통해 기존에 존재하던 격자‑의존적 오류 지시자와 달리 실제 절단오차를 정량적으로 얻을 수 있다. 2차 정확도를 갖는 차분식에 대해 수치 실험으로 검증하였다.

상세 분석

이 논문은 유한체적(FV) 방법이 구조화된 격자에서 2차 이상 정확도를 달성함에도 불구하고, 실제 시뮬레이션에서는 격자 비대칭, 스큐, 비균일성 등에 의해 지역적인 절단오차가 크게 변동한다는 점을 지적한다. 기존의 절단오차 지시자는 미분 연산자의 근사식만을 기반으로 하여 격자 형태에 대한 정보를 충분히 반영하지 못한다. 저자들은 “제한‑재구성”(restriction‑reconstruction) 아이디어를 차용해, 미세 격자 h에서 얻은 수치해 φ_h를 격자 연산 I_{h}^{2h}를 통해 거친 격자 2h로 제한하고, 그 격자에서 동일한 차분 연산자 N_{2h}를 적용한다. 이때 b_{2h}=∫{ΔΩ_P}b dΩ/ΔΩ_P 를 우변으로 사용하면, τ{2h}^h = b_{2h} – N_{2h}(I_{h}^{2h}φ_h) 라는 “상대 절단오차”를 직접 계산할 수 있다.

핵심적인 기술적 난관은 질량 플럭스(mass flux)의 압력 기여를 기존의 ‘momentum interpolation’(MIM) 방식에서 발생하는 압력‑계수 연관성을 제거하는 것이다. 저자들은 압력 부분을 전적으로 격자 기하와 제한 연산에만 의존하도록 재구성함으로써, N_{2h}와 N_h가 동일한 차분 스키마를 유지하도록 보장한다. 이는 공동 격자에서 흔히 발생하는 ‘pressure‑velocity decoupling’ 문제를 회피하면서도, 압력-속도 연계성을 정확히 유지한다는 점에서 혁신적이다.

이론 전개에서는 τ_h와 ε_h(해와 정확해의 차이) 사이의 관계를 Taylor 전개를 통해 τ_h = O(h^p) , ε_h = O(h^p) 로 보이며, 격자 정제 시 τ와 ε가 동일한 형태를 유지함을 증명한다. 이후 τ_{2h}^h 를 τ_h 로 환산하기 위해 τ_h ≈ (2^p/(2^p–1)) τ_{2h}^h 라는 스케일링 식(식 3.12)을 도출한다. 이는 격자 정제 비율 r에 대해 일반화될 수 있어, 임의의 r·h 격자에서도 동일한 추정식이 적용 가능함을 의미한다.

수치 실험에서는 2차 정확도를 갖는 차분식과 다양한 격자 비대칭·스큐 사례를 대상으로, 추정된 절단오차와 실제 절단오차(정밀 해와의 차이)를 비교하였다. 결과는 평균 절대 오차가 10% 이하이며, 특히 급격히 변하는 경계층 근처에서도 안정적인 추정이 가능함을 보여준다. 또한, 추정 오차는 제한 연산 I_{h}^{2h}의 차수가 차분식 차수(p)보다 높을 경우 더욱 감소한다는 실험적 근거를 제시한다.

이 방법의 장점은 (1) 추가적인 미세 격자 해를 요구하지 않으며, (2) 기존 FV 코드에 제한 연산과 τ_{2h}^h 계산 루틴만 삽입하면 바로 적용 가능하고, (3) 다중 격자(Multigrid) 프레임워크와 자연스럽게 연계될 수 있다는 점이다. 한계점으로는 (가) 압력‑속도 연계가 강한 비정상 흐름(예: 고레시스 수, 강한 회전)에서 제한 연산이 충분히 정확하지 않을 수 있으며, (나) 3차원 확장 시 격자 제한·보간 연산의 구현 복잡도가 급증한다는 점을 들 수 있다. 향후 연구에서는 고차 차분식, 비정형 격자, 그리고 동적 격자 적응(AMR) 환경에서의 적용 가능성을 탐색할 필요가 있다.