Computational Fluid Dynamics Comparisons forthe Virginia Tech Stability Wind Tunnel

기술 요약: 버지니아 공과대학교 안정성 풍동을 위한 전산유체역학 비교 연구

문제 정의
본 연구는 대규모 아음속 풍동에 대한 전산유체역학(CFD) 모델의 검증 및 타당성 확인(V&V)을 다룬다. 구체적으로, 본 연구는 비어 있는 상태의 버지니아 공과대학교 안정성 풍동(VTSWT)에 초점을 맞춘다. 주요 목적은 기하학적 충실도(geometric fidelity)의 필요 수준, 비균일 유입 경계 조건의 영향, 그리고 실제 시설 내 유동장을 예측하는 데 있어 다양한 난류 모델의 효용성을 결정하는 것이다. 본 연구는 9개 국제 팀의 계산 결과와 광범위한 실험 데이터를 비교함으로써 아음속 풍동 시설 모델링을 위한 최적의 실무 지침(best practices)을 수립하고자 한다.

방법론
본 연구는 레이놀즈 수 $1.7, 2.8, 3.9 \times 10^6 \text{ m}^{-1}$에서 VTSWT에 대한 CFD 시뮬레이션과 실험 측정값 사이의 체계적인 비교를 수행한다.

• 기하 구조: 두 가지 뚜렷한 기하 구조가 사용되었다: "설계된 형상(as-designed)"(이상적인 정방형 시험부)과 수축부 및 시험부의 매끄러운 광학 스캔으로부터 도출된 "제작된 형상(as-built)"이다. 디퓨저는 모든 계산에서 제외되었으며, 디퓨저가 시험부 상류에 미치는 영향이 미미한 것으로 밝혀짐에 따라 시험부를 하류 방향으로 연장하여 유출 경계로 활용하였다.
• 격자: 두 기하 구조 모두에 대해 체계적으로 세분화된 8개의 구조적 격자군(structured grids)이 생성되었다. 격자 크기는 약 124만 개의 셀(가장 거친 격자)에서 6억 3,800만 개의 셀(가장 미세한 격자)까지 다양했다. 모든 격자는 벽면 $y^+$ 값을 1 미만으로 유지하였으며 경계층 에지를 충분히 해상하였다.
• 참여자 및 모델: 9개의 계산 팀이 참여하였으며, 5가지의 서로 다른 CFD 코드(Fluent, STAR-CCM+, USM3D, FUN3D, Kestrel)와 6가지의 서로 다른 난류 모델을 활용하였다. 여기에는 Spalart–Almaras(SA) 변형 모델(QCR 포함), Menter SST, Langtry–Menter 천이 모델, Realizable $k$-$\epsilon$, 그리고 Reynolds Stress 모델(RSM)이 포함된다.
• 경계 조건: 시뮬레이션에는 균일 유입 조건과 수축부 상류의 실험 측정치로부터 도출된 비균일 유입 조건이 모두 적용되었다. 계산과 실험 간의 정확한 무차원화를 보장하기 위해 참조 조건(정압, 속도, 밀도)을 정밀하게 일치시켰다.

주요 결과

• 격자 민감도: 대부분의 난류 모델에서, Level 6(약 1억 6,000만 개 셀)보다 미세한 격자에서는 질량 유량 및 벽면 압력 계수에 대한 격자 민감도가 미미하게 나타났다.
• 기하학적 충실도 (설계된 형상 vs. 제작된 형상):
  • 압력 계수: 제작된 형상을 사용했을 때 터널 벽면의 0.61 m 패널 치수와 상관관계가 있는 국부적인 벽면 압력 계수의 변화가 나타났다. 그러나 시험부 내의 전체적인 압력 구배는 설계된 형상에 의해 예측된 것과 유사하게 유지되었다.
  • 경계층: 제작된 형상이 경계층 프로파일 및 코너 유동에 미치는 영향은 미미한 것으로 밝혀졌다.
  • 불일치: 계산 결과가 수축부에서 압력 계수를 약간 과대 예측하였다. 이러한 편향 오차(bias error)의 원인은 현재 조사 중이다.
• 난류 모델링:
  • 경계층: SST 모델과 Reynolds Stress 모델은 Spalart–Almaras 모델이 경계 에지를 부드럽게 만드는 경향이 있는 것과 달리, 더 날카로운 경계층 에지를 포착하였다.
  • 코너 유동: 비선형 모델(QCR 변형 모델)과 Reynolds Stress 모델은 코너에서의 이차 유동을 성공적으로 포착하였으나, 선형 에디 점성(linear eddy viscosity) 모델은 이를 포착하지 못했다. 그러나 이차 유동 모델링의 포함 여부는 비균일 유입 조건을 사용하는 것보다 실험적 코너 유동을 일치시키는 데 있어 덜 결정적인 요소였다.
• 유입 조건: 실험적으로 측정된 비균일 유입 경계 조건을 사용하는 것이 기하학적 충실도(설계된 형상 vs. 제작된 형상)나 난류 모델 변형을 선택하는 것보다 코너 유동 및 횡단 유동(crossflow) 스트림라인 패턴에 훨씬 더 큰 영향을 미쳤다. 특히, 실험 데이터에서 관찰되는 코너 영역의 저운동량 돌출부(low-momentum bulge)를 포착하기 위해서는 비균일 유입이 필수적이었다.
• 질량 유량: 서로 다른 난류 모델은 경계층 성장 예측의 차이로 인해 약간씩 다른 질량 유량을 예측하였다. Spalart–Almaras 모델은 일반적으로 SST 또는 천이 모델보다 약간 높은 질량 유량을 예측하였다.

의의 및 결론
본 논문은 아음속 풍동의 CFD 모델링을 위한 몇 가지 최적의 실무 지침을 제시한다:

1. 참조 조건: 특정 시뮬레이션 또는 실험에서 직접 추출한 참조 조건을 사용하여 모든 물리량을 무차원 형태로 보고하는 것이 중요하다. 참조 압력의 작은 오차는 동압(dynamic pressure)과 같은 유도된 물리량에서 상당한 오차를 초래할 수 있기 때문이다.
2. 기하 구조: 제작된 형상이 국부적인 표면 변화(패널 불일치 및 굽힘)를 포착하기는 하지만, 전체적인 유동 특성(압력 구배, 경계층)은 설계된 형상에 의해 잘 예측된다. 제작된 형상의 격자를 생성하는 노력은 매우 까다로우며, 항상 전역 유동 예측 정확도 향상에 비례하는 이득을 주는 것은 아니다.
3. 유입: 코너 유동의 특정 특징을 포착하는 데 있어서는 기하학적 충실도나 고급 난류 모델링보다 정확한 비균일 유입 모델링이 더 중요하다.
4. 모델링 범위: 시험부를 충분히 하류로 연장한다면, 비어 있는 상태의 풍동 시뮬레이션에서 디퓨저와 그 내부의 와류 발생기(vortex generators)를 생략해도 시험부 유동에 영향을 주지 않는다.
5. 향후 과제: 현재 생략되어 날카로운 모서리로 모델링된 코너 가스킷(corner gaskets)은 코너 유동에 상당한 영향을 미칠 수 있으며, 비구조적 격자(unstructured grids)를 이용한 조사가 필요하다.

본 연구는 고급 난류 모델과 제작된 형상이 가치 있는 국부적 세부 사항을 제공하기는 하지만, 비균일 유입의 정확한 표현이 비어 있는 풍동 시뮬레이션에서 실험적 코너 유동 데이터와 일치시키는 데 있어 가장 지배적인 요인이라는 결론을 내린다.