동적 서브스케일을 이용한 항별 변분 다중스케일 방법의 난류 항공역학 적용

초록

본 논문은 등차 순서의 속도‑압력 보간을 허용하면서도 대규모 3차원 외부 유동에서 안정적인 수치를 제공하는, 동적 서브스케일을 포함한 항별 변분 다중스케일(VMS) 안정화 방식을 제안한다. 압력‑보정 분할 단계에 직접 삽입된 최소한의 안정화 항들로 구성된 이 방법은 Ahmed body와 실제 Formula 1 차량에 대해 3 ~ 40 백만 요소의 비구조화 테트라hedral 메쉬에서 검증되었으며, 드래그·리프트 예측과 점별 속·압력 스펙트럼을 통해 관성 구간의 에너지 전달과 비해석적 소산 제어가 가능함을 보였다.

상세 분석

본 연구는 변분 다중스케일(VMS) 이론을 기반으로, 동적 서브스케일(dynamic subscale)과 직교 투영(orthogonal projection)을 결합한 항별(term‑by‑term) 안정화 프레임워크를 제시한다. 기존 VMS‑LES는 잔차 기반(residual‑based) 서브스케일을 사용해 전체 연산자를 수정하는 반면, 저자들은 각 연산자 블록(질량, 대류, 확산, 압력‑속도 결합)에 필요한 최소한의 항만을 삽입함으로써 ‘구조적 명료성’을 확보하였다. 이는 특히 압력‑보정(pressure‑correction) 분할 스킴에 적합하도록 설계되었으며, 각 단계에서 추가적인 속도‑압력 결합을 도입하지 않아 기존의 투영 기반 분할 알고리즘과 원활히 통합된다.

동적 서브스케일은 시간에 따라 진화하는 미세 스케일 변수로, BDF2 2차 정확도와 함께 해석적으로 해석된 스케일 분리 메커니즘을 제공한다. 직교 투영을 적용함으로써 해결 가능한(RESOLVED) 스케일과 서브스케일 사이의 에너지 교환이 명확히 정의되고, 비잔류 항이 아닌 ‘항별’ 형태로 소산(dissipation)이 발생한다. 이는 전통적인 인공 점성(artificial viscosity)이나 모델 기반 난류 모델에 비해 물리적 일관성을 유지하면서도 수치적 안정성을 확보한다는 장점을 가진다.

또한, 저자들은 선택적으로 압력 서브스케일에 grad‑div 항을 추가하여 비선형 수렴성을 강화하였다. 이 항은 기본 안정화 스킴에 필수적이지 않지만, 고속·고레놀즈 흐름에서 발생할 수 있는 압력 진동을 억제한다.

검증 사례로는 Ahmed body(슬랜트 각도 0°, 15°, 30°)와 실제 Formula 1 차량을 사용하였다. Ahmed body에서는 7.68 × 10⁵ Re에서 3 ~ 40 백만 요소 메쉬를 이용해 드래그 계수와 분리점 위치, 그리고 와류 구조를 정확히 재현하였다. 특히, 점별 속도·압력 스펙트럼을 통해 해결 가능한 주파수 대역에서 Kolmogorov −5/3 기울기를 보이며, 미해결 고주파 대역에서는 동적 서브스케일에 의해 적절히 소산이 이루어짐을 확인했다.

Formula 1 사례에서는 자유류속 56 m/s(≈10⁶ Re) 조건에서 전체 차량(전·후 윙, 바디, 휠 등)을 포함한 복합형 3차원 흐름을 시뮬레이션하였다. 5 ~ 10 백만 요소 메쉬에서도 압력‑보정 스킴이 수렴했으며, 전체 드래그·리프트와 휠 주변의 복잡한 와류 조직을 실험 데이터와 비교해 높은 일치도를 보였다. 특히, 대규모 병렬 실행 시 스케일링 효율이 80 % 이상 유지되어 실용적인 산업 적용 가능성을 시사한다.

요약하면, 이 논문은 (1) 항별 직교 VMS 안정화가 압력‑보정 분할 스킴에 자연스럽게 삽입될 수 있음을, (2) 동적 서브스케일이 물리적 소산 메커니즘을 제공하면서도 수치적 안정성을 보장함을, (3) 복잡한 외부 유동(Ahmed body, Formula 1)에서 실험적 정밀도와 병렬 효율성을 동시에 달성함을 입증하였다. 이러한 접근은 고레놀즈, 복합형 차량 흐름을 다루는 산업 CFD에서 기존 LES·RANS·Hybrid 모델을 대체하거나 보완할 수 있는 강력한 대안으로 기대된다.