전천후 로우형 스킴을 이용한 대규모 와류 시뮬레이션 개선 방안

초록

본 논문은 전통적인 로우(Roe) 충격 포획 스킴이 대규모 와류 시뮬레이션(LES)에서 과도한 수치 점성으로 k‑5/3 스펙트럼을 재현하지 못함을 분석한다. 로우 스킴을 기본 업윈(dissipation)·압력·속도 차이에 의한 다섯 부분(ξ, pUΔ, uUΔ, ppΔ, upΔ)으로 분해하고, 균질 감쇠 와류(HDT) 실험을 통해 각 부분의 LES 기여를 평가한다. 결과는 ξ와 upΔ가 주된 점성원이며, pUΔ는 안정성 확보에만 필요하고, 나머지는 LES에 거의 영향을 주지 않음을 보여준다. 이를 토대로 ξ의 계수를 0.5로 감소시키고 upΔ를 제거한 전천후 LES‑Roe 스킴을 제안, 저해상도 격자에서도 k‑5/3 스펙트럼과 와류 구조를 성공적으로 재현한다.

상세 분석

본 연구는 로우 스킴의 수치 점성을 체계적으로 해부하기 위해 스킴을 다섯 개의 물리적 의미를 가진 항으로 분해하였다. 첫 번째 항 ξ는 전통적인 업윈(upwind) 점성으로, 흐름의 안정성을 보장하지만 과도한 점성으로 인해 고주파 영역에서 에너지 스펙트럼이 -5 정도로 급격히 감소한다. 두 번째 항 pUΔ는 압력 차이에 의해 인터페이스 속도를 보정하는 항으로, 실제 물리적 점성은 거의 없으며, 저마하 흐름에서 체커보드 현상을 억제하는 역할만 수행한다. 세 번째 항 uUΔ는 속도 차이에 의한 보정으로, 실험 결과 거의 무시할 수 있는 점성을 제공한다. 네 번째 항 ppΔ는 압력 차이에 의한 압력 보정으로, 미세한 음의 점성을 나타내지만 LES에 큰 영향을 주지 않는다. 마지막 항 upΔ는 속도 차이에 의해 압력을 보정하는 항으로, 로우 스킴 특유의 비물리적 행동(non‑physical behaviour) 문제를 일으키는 주된 원인이다. 이 항은 O(c) 규모의 큰 점성을 포함해 고주파 에너지 소멸을 과도하게 촉진한다.

균질 감쇠 와류 실험에서는 9가지 케이스를 설계해 각 항의 독립적·조합적 영향을 검증하였다. 기본 중심 차분(Cen)만 사용한 경우와 pUΔ, uUΔ, ppΔ만 포함한 경우는 모두 고주파 영역에서 에너지 누적이 발생해 수치 발산을 초래한다. 반면 ξ만 사용한 경우와 upΔ만 사용한 경우는 각각 큰 점성과 비물리적 압력 변동을 보여주었으며, 특히 ξ는 점성만으로도 안정성을 유지하지만 스펙트럼이 크게 왜곡된다. upΔ는 ξ와 유사한 수준의 점성을 제공하지만, 비물리적 압력 스케일링을 야기한다.

이러한 분석을 바탕으로 저자들은 두 가지 핵심 개선 전략을 도출한다. 첫째, ξ의 계수를 0.5로 감소시켜 기본 업윈 점성을 절반으로 줄임으로써 고주파 영역에서의 과도한 점성을 완화한다. 둘째, upΔ 항을 완전히 제거해 비물리적 압력 스케일링을 억제한다. pUΔ는 여전히 유지해 체커보드 현상을 방지한다. 최종 제안된 전천후 LES‑Roe 스킴은 식 (29)와 (30)으로 표현되며, 저해상도(32³) 격자에서도 k‑5/3 스펙트럼을 넓은 파수대에 걸쳐 재현하고, 와류 구조와 스큐니스 텐서 통계도 기존 Smagorinsky 모델을 적용한 중심 차분(Cen‑SMA)과 거의 동등한 수준을 보인다.