Fast Prediction of Aero-Thermal Performance During Hypersonic Vehicle Design

기술 요약: 극초음속 차량 설계 시 공력-열 성능의 신속한 예측

문제 정의
극초음속 차량의 설계에는 개발 초기 단계에서 공력 및 표면 열유속(heat flux)을 신속하게 추정하는 과정이 필요하다. 고충실도 전산유체역학(CFD) 시뮬레이션은 상세한 정확도를 제공하지만, 예비 설계 단계에서 요구되는 광범위한 설계 반복 및 궤적 최적화에 적용하기에는 계산 비용이 너무 높다. 접선 쐐기(tangent-wedge), 수정 뉴턴 이론(modified Newtonian), 프란틀-마이어(Prandtl-Meyer)와 같은 국소 표면 경사 기법(local surface inclination techniques)을 사용하는 전통적인 공학 수준의 방법들은 속도는 빠르지만 점성 효과를 모델링하는 데 있어 충실도가 부족한 경우가 많다. 특히, 표준적인 접근 방식은 경계층이 완전히 층류이거나 완전히 난류라고 가정하는 경향이 있다. 완전 난류 흐름 가정은 보수적인 방법이지만, 이는 종종 피부 마찰 항력과 표면 열유속을 과대 예측하는 결과를 초러온다. 이는 열방호 시스템(TPS)의 과다 설계, 차량 질량 증가, 연료 소비 증가로 이어진다. 또한, 현재의 빠른 공학 도구 중 국소 경사 프레임워크 내에서 경계층 천이 효과를 통합하여 구현한 사례는 드물며, 천이 유동이 공력-열 하중에 미치는 영향은 매우 크다.

방법론
본 논문은 고충실도 CFD와 고전적 해석 방법 사이의 간극을 메우기 위해 설계된 경량 공력-열역학 분석 도구인 HipeX를 제시한다. HipeX는 무관한 유동 특성을 추정하기 위해 국소 표면 경사법에 기반한 하이브리드 공학 수준 접근 방식을 사용하며, 점성 효과 및 천이를 다루기 위해 혼합 경계층 모델(blended boundary-layer model)을 채택한다.

• 비점성 유동 모델링: 솔버는 비정형 표면 메쉬(STL 형식)를 처리하며, 입사 속도 벡터에 대한 국 local 표면 경사도에 따라 적절한 유동 모델을 선택한다:

  • 수정 뉴턴 이론(Modified Newtonian Theory): 박리 충격파(blatazed shocks, 둔체) 및 큰 편향각에 적용된다.
  • 경사 충격파 관계식(Oblique Shock Relations): 부착 충격파 및 압축 영역에 사용된다.
  • 프란틀-마이어 팽창(Prandtl-Meyer Expansion): 유동 팽창 영역에 적용되며, 수치적 안정성을 보장하기 위해 이론적 최대 편향각의 50%로 적용을 제한하는 안전장치를 둔다.
  • Gaubeaud 공식: Prandtl-Meyer 해법이 실패할 경우 사용하는 경험적 대체 모델이다.
  • 대기 물성은 비행 고도에 따른 표준 모델을 통해 결정된다.

• 점성 및 열 모델링:

  • 압축성: 층류 및 난류 영역 모두에 적용 가능한 에커트 참조 온도법(Eckert's reference temperature method)을 사용하여 보정한다.
  • 천이 모델링: 천이의 시작을 결정하기 위해 마하 수의 함수로서 경험적인 임계 레이놀즈 수($Re_{cr}$)를 계산한다. 모델은 0(층류)에서 1(난류) 사이의 간섭 계수(intermittency factor, $\gamma_{it}$)를 사용하여 정의된 천이 영역 내에서 층류 및 난류의 피부 마찰 계수를 혼합한다.
  • 열유속 계산: 레이놀즈 유추(Reynolds analogy)와 스탠튼 수(Stanton number)를 사용하여 피부 마찰 계수로부터 표면 열유속을 도출한다. 모델은 대류 열전달과 복사 열전달(방사율 0.4 가정)을 모두 고려하여, 젖은 면적(wetted area)에 대해 적분함으로써 순 열유속과 총 열하중을 계산한다.

• 워크플로우: 이 도구는 3D 형상과 비행 조건(마하 수, 고도, 받음각, 요각)을 입력값으로 받는다. 유동 특성을 계산하고, 점성 보정을 적용하며, 힘을 적분하여 공력 계수를 결정한 후, 결과를 VTK 형식으로 내보낸다. 코드는 궤적 분석을 위한 병렬 실행을 지원한다.

주요 기여

1. 혼합 천이 모델: 주요 기여는 빠른 국소 경사 프레임워크 내에 혼합 경계층 모델을 통합한 것이다. 이를 통해 솔버는 이진적인 층류/난류 가정에 의존하는 대신, 국소 레이놀즈 수와 마하 수를 기반으로 층류에서 난류로의 변화를 추정할 수 있다.
2. 계산 효율성: HipeX는 신속한 실행이 가능하도록 설계되어, 표준 워크스테이션에서 약 1분 만에 전체 궤적(약 300개 지점)에 대한 성능을 계산할 수 있다.
3. 벤치마크 데이터 검증: 본 도구는 다양한 마하 수(M=2 ~ M=7.5) 및 비행 궤적에 대해 HEXAFLY-INT 형상(마하 7~8 순항을 위한 혼합 리프팅 바디 설계)에 대한 CFD 데이터 및 실험 참조 자료를 통해 검증되었다.

결과
검증은 다양한 마하 수(M=2 ~ M=7.5)와 비행 궤적에 걸쳐 HEXAFLY-INT 형상을 대상으로 수행되었다.

• 공력 계수: 고마하수(M=6 및 M=7)에서 HipeX의 양력($C_L$), 항력($C_D$), 압력 중심($X_{cp}$) 예측은 CFD 데이터와 밀접하게 일치한다. 혼합 천이 해법은 순수 층류 또는 완전 난류 가정보다 일반적으로 더 정확한 결과를 제공한다. 낮은 마하수(M=4 및 M=2)에서는 절대적인 $C_L$ 및 $C_D$ 값에서 차이가 커지지만, 양항비($L/D$)는 상당히 정확하게 유지된다. 이 도구는 유동이 강한 극초음속 조건에서 벗어날 때 국소 경사법의 한계를 정확히 식별한다.
• 열 하중: 유동이 주로 층류인 경우(고고도), 혼합 해법은 완전 층류 사례와 유사한 결과를 낸다. 그러나 천이 영역에서는 혼합 모델이 층류와 난류 극한 사이의 중간 열유속 수준을 포착한다.
• TPS 영향: 연구 결과, 완전 난류 흐름을 가정하면 특히 전연(leading edge) 근처에서 가열을 크게 과대 예측한다는 것이 입증되었다. 궤적의 상당 부분에서 경계층의 층류 또는 부분 난류 특성을 포착함으로써, HipeX는 더 낮고 현실적인 열유속 추정치를 제공한다. 이는 안전 마진을 저해하지 않으면서도 TPS 질량을 줄일 수 있는 가능성을 시사한다.
• 천이 개시: BLITZ 코드와의 비교 결과, 천이 위치에 대해 유사한 경향을 보였으나, HipeX는 난류의 조기 발생를 예측하는 경향이 있다. 현재 모델은 특정 팽창 영역에서 재층류화(relaminarization) 예측을 방지하는 데 어려움이 있다.

의의 및 주장
본 논문은 HipeX를 초기 단계의 공력-열 평가 및 설계 최적화를 위한 견고한 토대로 규정한다. HipeX의 의의는 설계의 보수성과 임무 효율성 사이의 "균형 잡힌 절충안"을 제공한다는 데 있다. 완전 난류 흐름이라는 지나치게 보수적인 가정에서 벗어남으로써, 이 도구는 열방호 시스템의 크기를 더 정확하게 산정할 수 있게 하며, 잠재적으로 차량 질량과 연료 소비를 줄일 수 있게 한다.

저자들은 HFE (HipeX)가 충격파 상호작용, 박리 구역, 강력한 3차원 효과를 예측할 수 없는 국소 표면 경사법 고유의 한계를 가지고 있음을 인정한다. 또한, 표면 거칠기 및 자유류 교란에 따른 의존성 때문에 실제 비행 조건에서의 정확한 천이 개시 결정은 여전히 도전적인 과제이다. 그럼에도 불구하고, HipeX는 고전적 해석 방법과 고충실도 CFD 사이의 간극을 성공적으로 메움으로써, 예비 설계에 필수적인 전체 궤적 및 대규모 기하학적 구성의 신속한 평가를 가능하게 한다.