Stability-Based Optimization of Planar Quiet Wind Tunnel Nozzles: Two-Dimensional Design and Three-Dimensional Assessment

기술 요약: 평면형 정숙 초음속 풍동 노즐의 안정성 기반 최적화

문제 정의
정숙 초음속 풍동의 설계는 고속 흐름에서의 자연 경계층(BL) 천이를 조사하는 데 매우 중요하다. 기존 시설들은 지상-비행 간 외삽을 저해하는 높은 자유류 교란을 겪기 때문이다. 축대칭 노즐은 흐름 대칭성이 본질적으로 2차원(2D)인 반면, 평면형 노즐은 측벽(sidewalls)을 가지고 있어 본질적인 3차원(3D) 효과를 유발한다. 표준 설계 방법론은 종종 2차원 가정(스팬 방향 균일 흐름)에 의존하며, 수직 팽창에 따른 횡방향 팽창 부재로 인한 이방성 팽창을 간과한다. 이러한 단순화는 스팬 방향 압력 구배와 이차 유동 구조를 생성하며, 이는 2차원 모델에는 존재하지 않는 크로스플로우(crossflow)와 같은 불안정성 메커니즘을 유발할 수 있다. 본 연구에서 다루는 핵심 문제는 3차원 측벽 효과가 존재하는 평면형 노즐에 2차원 안정성 기반 최적화 전략을 적용했을 때 그 견고성(robustness)이 유지되는지 여부이다.

방법론
본 연구는 베이지안 최적화(BO) 루프 내에 통합된 모듈형 안정성 기반 최적화 프레임워크를 채택한다. 과정은 세 가지 주요 단계로 구성된다:

1. 비점성 윤곽 생성: HYPNOZE 툴킷은 특성법(Method of Characteristics, MoC)을 사용하여 기하학적 파라미터(목 곡률, 개구각, 직선 구간 길이) 및 운전 조건에 따른 노즐 윤곽을 생성한다.
2. 점성 경계층 분석: DEKAF 솔버는 층류 경계층 프로파일과 변위 두께를 계산하여 유효 형상에 대한 점성 보정을 제공한다.
3. 안정성 평가: VESTA 툴킷은 국소 선형 안정성 이론(LST)을 적용하여 교란의 증폭을 평가한다. 분석은 1차 모드(first-mode), 2차 모드(second-mode), 고틀러(Görtler) 불안정성을 대상으로 한다.

최적화는 약한 증폭에 대해 둔감한 제곱합의 제곱근(root-sum-square) 공식의 한계를 피하기 위해 N-인자들의 선형 합($N_{tot} = N_{1st} + N_{2nd} + N_{G\ddot{o}}$)으로 정의된 비용 함수를 최소화한다. 설계 변수는 노즐 개구각($\omega$), 직선 구간 길이($X_{BC}$), 그리고 매개변수화된 벽면 온도 분포($T_w(x)$)이다.

견고성을 평가하기 위해, 최적화된 2차원 설계들을 측벽을 명시적으로 포함하는 완전한 3차원 전산유체역학(CFD) 시뮬레이션(CFD++ 사용)을 통해 재평가한다. 이러한 3차원 기저 유동에 대한 안정성 분석에는 측벽을 따라 발생하는 이차 유동을 포착하기 위해 LST와 선형 포물선 안정성 방정식(LPSE)이 모두 사용된다.

주요 기여

• 2차원 최적화 프레임워크: 본 연구는 기하학적 제어와 능동적 열 관리를 결합한 2차원 안정성 기반 최적화가 평면형 노즐 상단 벽의 1차 모드 및 고틀러 불안정성을 성공적으로 억제하여, 아임계(subcritical) 증폭 수준($N_{combined} < 7.5$)을 달고 정숙 흐름 영역을 확장할 수 있음을 입증한다.
• 3차원 견고성 평가: 본 논문은 2차원 설계를 3차원 환경에서 체계적으로 평가하여, 상단 벽의 안정성 예측은 정확하게 유지되는 반면, 측벽 영역은 이차 유동에 의해 구동되는 뚜렷한 불안정성 메커니즘을 보임을 밝혀냈다.
• 3차원 불안정성에 대한 열 제어: 2차원 메커니즘(특히 목 하류의 급격한 온도 강하)을 위해 최적화된 벽면 온도 분포가 측벽을 따른 크로스플로우형 증폭 또한 감쇄시킨다는 중요한 발견은, 열 제어와 3차원 불안정성 메커니즘 사이의 결합을 시사한다.
• 안정성 이론의 유효 범위: LST와 LPSE의 비교는 3차원 평면 유동에서 선형 안정성 이론의 유효 영역을 구분한다. LPSE 통합은 강한 비평행성(non-parallelism)과 급격한 층류 방향 변화(일반적으로 변곡점 하류)가 발생하는 영역에서 무너지며, 이는 더 높은 충실도의 접근 방식이 필요한 지점을 나타낸다.

결과

• 최적화 결과: 최적의 구성(IT048)은 얕은 개구각($\approx 2^\circ$), 짧은 직선 구간, 그리고 특정 열 프로파일을 특징으로 한다: 1차 모드 파동을 억제하기 위해 목 근처에서는 높은 온도($350$K)를 유지하고, 이후 불안정성을 더욱 감쇄시키기 위해$\approx 292$ K로 급격히 떨어뜨린 뒤 출구 쪽으로 갈수록 약간 상승하는 형태이다. 이 설계는 조기 천이를 예측하는 참조 케이스와 달리 상단 벽을 따라 층류를 유지한다.
• 상단 벽 성능: 3차원 CFD 시뮬레이션에서 최적화된 설계의 상단 벽 안정성 특성(N-인자 진화)은 2차원 예측과 일치하며, 상대 오차는 5% 미만이다. 따라서 2차원 프레임워크는 상단 벽에 대한 신뢰할 수 있는 1차 근사치 역할을 한다.
• 측벽 성능: 측벽 경계층은 코너 와류와 횡방향 속도 성분을 포함한 강한 3차원 효과를 보인다. 이는 2차원 모델이 포착하지 못하는 크로스플로우형 불안정성을 유발한다. 최적화된 기하 구조와 열 제어가 참조 케이스에 비해 증폭 수준을 줄여주기는 하지만, N-인자 임계값($N_{crit} \approx 7.5$)은 측벽을 따라 약 $x/L \approx 0.5$ 지점에서 여전히 도달한다. 이는 최적화된 설계에서도 측벽 유도 천이가 의도된 균일 핵 영역보다 상류에서 발생할 수 있음을 시사한다.
• 메커니즘 상호작용: 최적화된 열 프로파일은 1차 모드와 고틀러 증폭 영역의 공간적 중첩을 효과적으로 분리한다. 또한, 2차원 물리학으로부터 도출된 열 제어 전략은 측벽을 따른 크로스플ow 증폭을 감소시킴으로써 이차적인 이점을 제공한다.

의의 및 주장
본 논문은 안정성 기반 최적화가 평면형 정숙 노즐의 상단 벽 천이를 제어하기 위한 효과적이고 물리적으로 해석 가능한 도구라고 주장한다. 연구는 2차원 최적화가 상단 벽 설계에는 충분하지만, 측벽 유도 크로스플로우 불안정성을 포착하지 못하기 때문에 완전한 노즐 평가에는 불충분하다는 점을 확립한다.

저자들은 평면형 구성을 위한 신뢰할 수 있는 노즐 설계에는 3차원 효과에 대한 명시적인 고려가 필요하다고 결론짓는다. 현재의 2차원 최적화가 견고한 기준(baseline)을 제공하기는 하지만, 측벽 천리 위험을 고려할 때 향후 설계 프레임워크에 면 방향(plane-marching) PSE나 스팬 방향 BiGlobal 안정성 이론과 같은 고충실도 접근 방식을 통합하여 다방향 불안정성 역학을 완전히 해결해야 한다. 본 연구는 벽면 온도 제어가 2차원 및 3차원 불안정성 메커니즘 모두에 영향을 미치는 다재다능한 설계 레버임을 강조하며, 측벽 천이를 완화하기 위한 새로운 경로를 제시한다.