비구조 격자를 위한 2차 유한체적 스킴의 개선형 리미터 연구

초록

본 논문은 비구조 격자 위에서 2차 정확도를 갖는 유한체적법에 적용되는 세 가지 리미터(Venkatakrishnan, Wang‑수정형, Nishikawa R³)를 비교한다. 전이음속 NACA 0012 에어포일을 세 가지 공격각에서 RANS‑SA‑neg 모델로 시뮬레이션하고, 각 리미터의 수렴성, 인공 점성, 충격파 포획 능력을 평가한다. 결과는 적절한 파라미터 범위 내에서 세 리미터가 모두 실험 데이터와 일치하지만, 인공 점성 및 수렴 속도에서 차이를 보임을 보여준다.

상세 분석

이 연구는 비구조 격자에서 2차 정확도의 MUSCL‑형 선형 재구성을 구현할 때 필수적인 리미터 함수의 성능을 체계적으로 검증한다. 첫 번째로 Venkatakrishnan 리미터는 기존에 널리 사용되던 형태로, 셀 내부와 이웃 셀 간의 최소·최대값을 이용해 제한계수를 계산한다. 이때 ε² 항을 도입해 거의 일정한 영역에서 수치적 불안정을 방지한다. 그러나 ε²가 셀 크기에 비례하도록 정의되면, 격자 비균일성이 큰 경우 ε² 값이 급격히 변동해 제한계수가 과도하게 억제되거나 반대로 완화되는 현상이 발생한다. 이를 보완하기 위해 Wang은 전역 최소·최대값을 이용한 ε² 정의를 제안했으며, 이는 격자 비균일성에 대한 민감도를 크게 낮춘다.

Nishikawa가 최근 제시한 R³ 리미터는 p=3인 Rᵖ 계열 중 하나로, Δ⁺와 Δ⁻의 절대값을 3제곱 형태로 결합하고, ε³ 항을 통해 작은 변화에도 부드러운 전이를 보장한다. 이론적으로는 3차 정확도를 유지할 수 있으나, 현재 논문에서 사용된 전체 스키마는 2차에 머물기 때문에 R³ 리미터는 추가적인 인공 점성을 제공한다. 실제 시뮬레이션에서는 R³ 리미터가 충격파 근처에서 가장 낮은 인공 점성을 보이며, 결과적으로 압력·온도 급변 구간을 더 선명하게 포착한다.

세 리미터 모두 동일한 재구성 프레임워크(Barth‑Jespersen 기반)와 Roe 플럭스 차분을 사용했으며, 시간 적분은 암시적 Euler와 GMRES(m=200)으로 수행했다. 수렴성 측면에서, 제한계수를 0으로 설정하면 1차 스킴으로 전환되어 잔차가 10⁻⁶ 수준까지 감소하지만 수렴 속도가 매우 느리다. 반면 제한을 완전히 해제하면(ψ=1) 충격파 근처에서 비물리적 진동이 급증해 발산한다. 특히 공격각 2.03°(구성 3)에서는 모든 리미터가 수렴 정체 현상을 보였으며, 이는 높은 CFL(≤10⁴)과 제한된 2차 재구성의 조합이 비선형 잔차를 충분히 억제하지 못하기 때문이다.

파라미터 민감도 실험에서는 εV∈