Tag: 압축

벨 노즐의 전산 해석

Computational Analysis of Bell Nozzles

본 보고서는 로켓 엔진에서 널리 사용되는 벨형 노즐의 설계 및 비설계 조건에서의 유동 특성을 전산 유체 역학(CFD)을 통해 분석한 연구 결과를 담고 있습니다. 연구의 핵심은 특성곡선법(MOC)을 이용한 노즐 최적화와 복잡한 충격파 구조의 수치적 규명에 있습니다.

Paper Metadata

• Industry: 항공우주 (Aerospace)
• Material: 티타늄 (Titanium), 공기 (Air), 수소 (Hydrogen)
• Process: 전산 유체 역학 (CFD), 특성곡선법 (MOC)

Keywords

• 벨형 노즐 (Bell type Nozzle)
• 수치 해석 (Numerical Analysis)
• 압축파 (Compressive Waves)
• 충격파 (Shock Wave)
• 특성곡선법 (Method of Characteristics)
• 유동 박리 (Flow Separation)

Executive Summary

Research Architecture

본 연구는 벨형 노즐 내부의 2차원 축대칭 유동을 분석하기 위해 수치적 접근 방식을 채택하였습니다. C 언어로 작성된 특성곡선법(MOC) 및 유선 함수 코드를 사용하여 고효율 노즐 윤곽을 설계하였으며, GAMBIT 2.4.6을 통해 격자를 생성하고 FLUENT 6.3.26을 사용하여 시뮬레이션을 수행하였습니다. 해석은 냉각 유동(300K 공기)과 가열 유동(1000K 수소)의 두 가지 조건에서 설계 압력 및 다양한 비설계 배압(5, 7, 11, 15 bar)을 대상으로 진행되었습니다.

Key Findings

수치 해석 결과, 높은 노즐 압력비(NPR)에서 점성 유동에 의한 비대칭 람다 충격파가 포착되었습니다. 설계 조건에서 마하수 3.0을 목표로 하였으나, 실제 수치 해석에서는 냉각 유동 시 2.932, 가열 유동 시 2.897의 마하수가 도출되어 이론값과 약 2.2~3.4%의 오차를 보였습니다. 비설계 조건에서는 배압이 15 bar에서 5 bar로 감소함에 따라 충격파가 노즐 출구 방향으로 이동하며 강도가 세지는 현상이 관찰되었으며, 1D 이론 대비 약 20%의 위치 편차가 확인되었습니다.

Industrial Applications

본 연구의 결과는 로켓 엔진 노즐 설계 시 비설계 작동 조건에서 발생하는 성능 손실을 예측하고 최소화하는 데 활용될 수 있습니다. 특히 충격파와 경계층의 상호작용으로 발생하는 측면 하중(side loads) 및 유동 박리 지점을 정확히 파악함으로써 노즐의 구조적 안정성과 추진 효율을 높이는 기초 설계 데이터로 사용 가능합니다.

Theoretical Background

특성곡선법 (Method of Characteristics, MOC)

MOC는 초음속 유동 문제를 해결하기 위한 수치적 절차로, 균일하고 파동이 없는 초음속 유동을 생성하는 노즐 윤곽을 설계하는 데 필수적입니다. 이 방법은 유동장의 특성 라인을 따라 미분 방정식을 적분하여 노즐의 확산부 좌표를 계산하며, 등엔트로피, 비점성, 비회전 유동을 가정하여 최적의 팽창 효율을 제공하는 곡선을 정의합니다.

유동 박리 및 충격파 구조 (Flow Separation and Shock Structures)

노즐 내부 배압이 설계치보다 높을 경우 과팽창(Over-expanded) 상태가 되어 유동 박리가 발생합니다. 이때 자유 충격 박리(FSS)와 제한 충격 박리(RSS) 현상이 나타나며, 이는 람다 충격파(Lambda shock)와 마하 디스크(Mach disk)를 형성합니다. 이러한 복잡한 충격파 구조는 노즐 벽면의 압력 분포를 급격히 변화시키며 엔진의 추력 성능에 직접적인 영향을 미칩니다.

Results and Analysis

Experimental Setup

수치 해석을 위해 노즐 벽면 소재는 티타늄으로 설정하였으며, 25,800개의 구조적 격자(Structured grid)를 생성하여 해석의 정밀도를 높였습니다. 난류 모델로는 유동 박리 예측에 탁월한 성능을 보이는 SST k-ω(Shear Stress Transport) 모델을 적용하였습니다. 경계 조건으로 입구는 압력 입구(Pressure Inlet), 출구는 압력 출구(Pressure Outlet)를 사용하였으며, 2차 상향풍 스킴(Second order upwind scheme)을 통해 수렴성을 확보하였습니다.

Visual Data Summary

마하수 등고선 분석 결과, 노즐 목(Throat) 영역에서 반원형의 유동 특성이 관찰되었으며 확산부 중반 이후 유동 변화가 완만해지는 벨 노즐 특유의 설계 특징이 확인되었습니다. 비설계 조건인 15 bar 배압에서는 노즐 내부 깊숙이 충격파가 형성되었으나, 배압이 낮아질수록 충격파가 출구 쪽으로 밀려나며 유동장이 확장되는 시각적 변화가 뚜렷하게 나타났습니다.

Variable Correlation Analysis

배압(Back Pressure)과 충격파 위치 사이에는 강한 상관관계가 존재합니다. 배압이 낮아질수록 노즐 내부의 정압과 외부 압력의 차이가 줄어들어 충격파 발생 지점이 하류로 이동합니다. 또한, 가열 유동(수소)의 경우 냉각 유동(공기)보다 전압 손실이 약 1% 정도 더 크게 발생(5.46% vs 4.49%)하는 것으로 나타났는데, 이는 고온 유동에서의 경계층 효과와 점성 소산이 더 활발하기 때문으로 분석됩니다.

Paper Details

Computational Analysis of Bell Nozzles

1. Overview

• Title: Computational Analysis of Bell Nozzles
• Author: Beena D. Baloni, Sonu P. Kumar, S. A. Channiwala
• Year: 2017
• Journal: Proceedings of the 4th International Conference of Fluid Flow, Heat and Mass Transfer (FFHMT’17)

2. Abstract

벨형 노즐은 로켓 노즐에서 가장 흔히 사용되는 형태입니다. 이 유형의 노즐은 원추형 노즐에 비해 크기와 성능 면에서 상당한 이점을 제공할 뿐만 아니라 환형 노즐에 비해 복잡성을 줄여줍니다. 노즐은 연소실에서 생성된 정체 온도($T_o$)와 정체 압력($P_o$)을 사용하여 연소 가스를 높은 초음속 속도로 가속함으로써 추력을 생성합니다. 노즐 팽창비는 출구 속도에 의해 결정됩니다. 비행 중 제트 유동은 짧은 기간 동안만 주변 유동에 이상적으로 팽창 및 적응하며, 나머지 시간 동안 로켓 엔진은 비설계 조건에서 작동합니다. 본 연구는 GAMBIT 2.4.6 및 FLUENT 6.3.26 소프트웨어를 사용하여 설계 및 비설계 조건에서 벨형 노즐 내의 2D 축대칭 유동 해석을 포함합니다. 해석을 위한 고효율 노즐 윤곽을 정의하기 위해 특성곡선법과 유선 함수를 사용하는 컴퓨터 코드가 개발되었습니다. 시뮬레이션은 냉각 및 가열 유동 조건에 대해 별도로 수행되었습니다. 유동 해석을 위해 SST k-ω 난류 모델이 선택되었습니다. 수렴된 솔루션은 점성 유동에 대해 높은 노즐 압력비(NPR)에서 노즐 내 비대칭 람다 충격파를 포착했습니다. 또한 NPR에 따른 후속 충격파 및 유동 박리를 예측했습니다. 점성 예측에서 마하 스템의 법선 충격파 강도는 일반적으로 NPR이 증가함에 따라 증가합니다. 충격 구조, 충격 위치, 법선 충격 크기, 후속 충격 및 비대칭 람다 충격 측면에서 예측된 시뮬레이션과 분석 결과 사이에 양호한 일치가 관찰되었습니다.

3. Methodology

3.1. 노즐 윤곽 설계: C 언어 프로그램을 통해 MOC와 유선 함수를 결합하여 마하수 3.0에 최적화된 벨 노즐 좌표를 생성함.
3.2. 수치 도메인 및 격자 생성: GAMBIT 2.4.6을 사용하여 25,800개의 사각형 구조적 격자를 생성하고 노즐 목 부근에 격자를 조밀하게 배치함.
3.3. 물리 모델 설정: FLUENT 6.3.26에서 2D 정상 상태 RANS 방정식을 풀었으며, 박리 유동 해석에 적합한 SST k-ω 난류 모델을 적용함.
3.4. 경계 조건 및 수렴: 압력 입구 및 출구 조건을 설정하고, 잔차 기준을 $10^{-4}$로 설정하여 수렴될 때까지 반복 계산을 수행함.

4. Key Results

설계 조건에서 냉각 유동의 출구 마하수는 2.932, 가열 유동은 2.897로 측정되어 이론적 설계치인 3.0과 매우 근접한 결과를 얻었습니다. 비설계 조건(배압 5~15 bar) 분석에서는 배압이 높을수록 노즐 내부에서 충격파가 일찍 발생하며, 이로 인해 유동 박리가 유도됨을 확인하였습니다. 특히 15 bar 조건에서 충격파 이후 마하수가 2.66에서 0.498로 급격히 감소하는 현상이 관찰되었습니다. 수치 해석 결과는 1D 비점성 이론과 비교했을 때 충격파 위치에서 약 20%의 편차를 보였는데, 이는 실제 유동의 점성 효과와 경계층 성장을 반영한 결과입니다.

Figure List

1. Fig. 1: 벨 노즐의 경계 조건이 포함된 유체 도메인
2. Fig. 2: 냉각 유동에 대한 벨 노즐 축 방향 마하수 변화
3. Fig. 3: 냉각 유동의 다양한 비설계 배압에 따른 마하수 등고선
4. Fig. 4: 가열 유동의 다양한 비설계 배압에 따른 마하수 등고선
5. Fig. 5: 15 bar 배압에서 노즐 축 및 벽면의 정압 변화

References

1. G. Hagemann et al., “Advanced Rocket Nozzles,” Journal of Propulsion and Power, 1998.
2. G. P. Sutton, “Rocket Propulsion Elements,” 7th ed, 2001.
3. J. Ostlund, “Flow Process in Rocket Engine Nozzles with Focus on Flow Separation and Side-Loads,” 2002.

Technical Q&A

Q: SST k-ω 난류 모델을 사용한 특별한 이유가 있습니까?

SST k-ω 모델은 표준 k-ε 모델의 원거리 유동 해석 장점과 k-ω 모델의 벽면 근처 해석 장점을 결합한 모델입니다. 특히 로켓 노즐 내부에서 발생하는 강력한 역압력 구배에 의한 유동 박리 현상을 예측하는 데 있어 다른 난류 모델보다 실험 데이터와 더 높은 일치성을 보이기 때문에 선택되었습니다.

Q: 냉각 유동과 가열 유동의 해석 결과에서 나타난 주요 차이점은 무엇입니까?

가열 유동(수소, 1000K)은 냉각 유동(공기, 300K)에 비해 정체 압력 손실이 약 1% 더 높게 나타났습니다(5.46% vs 4.49%). 이는 고온 유동에서 점성 효과가 더 크게 작용하여 경계층 두께에 영향을 미치기 때문이며, 이로 인해 출구 마하수 또한 가열 유동에서 약간 더 낮게 관찰되었습니다.

Q: 비설계 조건에서 발생하는 람다 충격파(Lambda shock)의 구조적 특징은 무엇입니까?

람다 충격파는 입사 충격파, 반사 충격파, 그리고 마하 스템(Mach stem)이 만나는 삼중점(Triple point)을 형성하는 구조를 가집니다. 본 연구에서는 배압 5 bar 조건에서 이러한 분기된 구조를 명확히 포착하였으며, 이는 노즐 벽면의 경계층 박리와 밀접하게 연관되어 비대칭적인 압력 분포를 유발합니다.

Q: 1D 이론적 결과와 CFD 수치 해석 결과 사이에 편차가 발생하는 원인은 무엇입니까?

1D 이론은 유동을 비점성(Inviscid) 및 등엔트로피 과정으로 가정하지만, 실제 CFD 해석에서는 벽면 마찰에 의한 경계층 성장과 점성 소산 효과를 고려합니다. 이러한 점성 효과가 유동의 유효 단면적을 줄이고 에너지 손실을 유발하여 충격파의 위치와 강도에서 약 20%의 차이를 만들게 됩니다.

Q: 연구에서 언급된 ‘측면 하중(Side loads)’은 노즐에 어떤 영향을 미칩니까?

비설계 조건에서 유동 박리가 비대칭적으로 발생하면 노즐 벽면에 불균일한 압력 분포가 형성됩니다. 이는 노즐 구조물에 수직 방향의 힘인 측면 하중을 발생시키며, 심할 경우 노즐의 기계적 파손이나 비행 궤적의 불안정성을 초래할 수 있습니다.

Conclusion

본 연구는 벨형 로켓 노즐의 설계 및 비설계 작동 조건에 대한 2차원 축대칭 수치 해석을 성공적으로 수행하였습니다. MOC를 통해 설계된 노즐은 목표 마하수에 근접한 성능을 보였으며, SST k-ω 모델을 통해 과팽창 조건에서의 복잡한 충격파 구조와 유동 박리 현상을 정밀하게 모사하였습니다.

결론적으로, 노즐 내부의 충격파 위치와 강도는 배압에 따라 민감하게 변하며, 이는 1D 이론만으로는 예측하기 어려운 점성 효과를 포함하고 있습니다. 본 연구 데이터는 향후 다양한 고도에서 작동해야 하는 로켓 엔진의 노즐 최적화 및 구조적 안정성 검토를 위한 중요한 수치적 근거를 제공합니다.

Source Information

Citation: Beena D. Baloni, Sonu P. Kumar, S. A. Channiwala (2017). Computational Analysis of Bell Nozzles. Proceedings of the 4th International Conference of Fluid Flow, Heat and Mass Transfer (FFHMT’17).

DOI/Link: 10.11159/ffhmt17.110

Technical Review Resources for Engineers:

▶ Access the original research paper (PDF)
▶ FLOW-3D 솔루션 팀과 협의하여 기술적 타당성을 검토하시려면..

This material is for informational purposes only. Unauthorized commercial use is prohibited.
Copyright © 2025 STI C&D. All rights reserved.