저감소수도 수치 스키마 비교 연구

초록

본 논문은 고전적 HLL 스키마와 대비하여 HLLC, TV 플럭스 분할, 저감소수도 중앙-업윈드(LDCU), 그리고 특성분해 기반 중앙-업윈드(LCDCU) 네 가지 저감소수도 방법을 1차부터 5차까지의 정확도로 확장하고, 1·2차원 유체역학 방정식에 적용해 정확도·견고성·계산 효율성을 종합적으로 평가한다.

상세 분석

논문은 먼저 HLL 스키마의 구조적 한계, 즉 접촉파와 전단파가 두 파동 모델에 포함되지 않아 과도한 확산이 발생한다는 점을 지적하고, 이를 보완하는 HLLC, TV 분할, LDCU, LCDCU 네 가지 저감소수도 스키마를 체계적으로 소개한다. HLLC는 중간 접촉파를 명시적으로 재구성함으로써 접촉면을 선명하게 포착하면서도 HLL의 견고함을 유지한다. TV 분할은 플럭스를 충격, 접촉, 팽창 파동으로 분해해 각각에 맞는 전파 속도를 적용해 확산을 최소화하고, 특히 정지 접촉면을 정확히 보존한다. LDCU는 중앙-업윈드 프레임워크에 서브셀 해상도 기반 투영 단계와 안티디퓨전 항을 도입해 접촉·전단 파동을 날카롭게 재현하면서 비진동성을 확보한다. LCDCU는 확산 항에 지역 특성 분해를 적용해 파동별 확산을 조절, 특히 고속 충격 영역에서 과도한 확산을 억제한다.

다음으로 저자들은 이들 스키마를 1차~5차 정확도로 확장한다. 2차는 제한자를 이용한 선형 재구성, 3차·5차는 A‑WENO 기반의 고차 차분 보정항(두 번째 및 네 번째 미분)과 비선형 가중치를 적용한다. 재구성은 특성 변수에 적용해 진동을 억제하고, θ 파라미터를 통해 확산 강도를 조절한다.

수치 실험에서는 1차원 소음 파동, Sod 충격, Lax 충격, 복합 파동, 그리고 2차원 복셀, 이중 충격, 켈리-헬름홀츠 흐름, 복합 난류와 같은 표준 벤치마크를 사용한다. 결과는 HLLC와 TV 분할이 접촉·전단 파동을 가장 선명히 포착하지만, 강한 충격 앞에서는 약간의 비진동성 문제를 보인다. LDCU와 LCDCU는 충격 전후 모두 안정적이며, 특히 LCDCU는 고차 정확도에서 가장 낮은 전반적 확산을 보인다. 계산 비용 측면에서는 HLLC와 TV 분할이 가장 저렴하고, LCDCU가 약간 높은 연산량을 요구한다(특성 분해와 고차 보정 때문). 전반적으로 네 스키마는 모두 기존 HLL 대비 확연히 낮은 수치 확산을 보이며, 선택은 문제의 파동 특성(접촉·전단 강조 vs. 강충격 강조)과 허용 가능한 계산 비용에 따라 달라진다.