고차원 유한차분 압축성 흐름 해석기 검증

초록

본 논문은 6차 고차원 컴팩트 유한차분 스킴을 기반으로 한 압축성 흐름 해석기의 검증·검증을 수행한다. 1차원 소드 충격관, 2차원 충격‑전단층 상호작용, 압축성 채널 흐름, 압축성 난류 경계층, 그리고 충격‑난류 경계층 상호작용 등 다섯 가지 표준 사례에 대해 정확한 충격 포획, 와류 구조 해상도, 1·2차 통계량의 일치성을 확인하였다.

상세 분석

이 연구는 Lele가 제안한 6차 컴팩트 유한차분 스킴을 핵심으로, 고해상도 스펙트럼‑유사 정확도를 유지하면서도 충격 파동을 안정적으로 포획할 수 있는 수치 기법을 구현하였다. 시간 적분은 3차 명시적 Runge‑Kutta 방식을 사용하고, 고차원 비분산 필터와 Kawamura‑Lele 인공 확산 기법을 결합해 고주파 진동을 억제하면서도 물리적 소산을 최소화하였다. 트라이디아고날 시스템을 병렬 Thomas 알고리즘으로 해결함으로써 CPU·GPU 클러스터에서의 확장성을 확보하였다.

검증 사례는 압축성 흐름의 대표적 난제들을 포괄한다. 1차원 소드 충격관에서는 200∼600 격자에서 정확한 라인어 솔루션과 5차 WENO와의 비교를 통해 격자 수가 증가할수록 수치해가 정확히 수렴함을 보였다. 2차원 충격‑전단층 상호작용에서는 동일 격자에서 기존 5차 WENO‑Z 대비 더 선명한 수직 와류와 밀도 변동을 재현했으며, 고해상도 격자와 비교했을 때도 차이가 미미함을 확인하였다.

압축성 채널 흐름 DNS에서는 Mb=1.5, Reb=6000(≈Reτ 218) 조건에서 기존 Morrison·Mo desti 데이터와 평균 속도, Van‑Driest 변환 속도, 온도, rms 변동량이 거의 일치하였다. 이는 근벽층 해상도(Δy⁺≈0.5)와 압축성 처리의 정확성을 입증한다.

압축성 난류 경계층에서는 M∞=2.25, Reθ≈2100 조건에서 1300×250×301 격자를 사용했으며, 비정상적 트리핑 기법을 통해 완전 난류 상태를 유도하였다. Van‑Driest 변환 평균 속도와 밀도 스케일링된 Reynolds 정규응력 프로파일이 Poggie·Pirozzoli·Schlatter 등 기존 DNS와 매우 유사했다. 특히 스트림와이즈 Reynolds 정규응력 피크를 비교적 조밀하지 않은 격자에서도 정확히 포착했다.

마지막으로 24° 압축 램을 통한 충격‑난류 경계층 상호작용 DNS에서는 실험 데이터(Bookie et al.)와 거의 동일한 M∞≈2.9, Reθ≈2400 조건을 재현했으며, 평균 및 변동 통계가 실험·다른 DNS와 일치하였다. 전체적으로 고차원 컴팩트 스킴이 충격 포획과 난류 구조 해상도 사이의 전형적인 트레이드오프를 효과적으로 해소함을 보여준다.