로켓 노즐 초음속 유동과 배기에서 발생하는 팽창파 해석

MAE_290C_Final.pdf

2.66MB

1. 수치 해석적 정밀도 (Numerical Methods)

이 연구의 핵심은 맥코맥 방법(MacCormack Method)을 사용하여 시공간에서 모두 2차 정확도(2nd-order accuracy)를 확보하려 했다는 점이다.

• 예측자-교정자 기법: 유동 방정식을 풀 때, 먼저 전방 차분(Forward Difference)으로 다음 단계의 값을 예측(Predictor)한 뒤, 후방 차분(Backward Difference)으로 그 값을 수정(Corrector)한다. 이를 통해 수치적 불안정성을 줄이고 정확도를 높였다.
• 물리량 산출: 유체의 밀도($\rho$), x·y축 속도($u, v$), 온도($T$), 압력($p$)을 동시에 계산하며, 실제 기체에 가깝게 계산하기 위해 서덜랜드 법칙(Sutherland’s Law)으로 점성 계수를, 이상 기체 상태 방정식으로 압력을 구했다.

2. 격자 변환과 블록 구조 (Grid & Mesh)

단순한 사각형 계산 영역을 복잡한 노즐 모양에 맞추기 위해 두 가지 고급 기술이 사용되었다.

• 격자 변환 (Grid Transformation): 수학적인 매핑(Mapping)을 통해 직교 좌표계(x, y)를 노즐의 곡선 형태를 따라가는 물리적 좌표계($\xi, \eta$)로 변환했다. 이를 통해 곡면 형태의 노즐 벽면에서도 정확한 계산이 가능해진다.
• 블록 구조와 고스트 레이어 (Block Structure & Ghost Layer):
  • 전체 영역을 4개의 블록(노즐 내부와 외부 배기 영역들)으로 나누어 계산했다.
  • 각 블록이 만나는 경계에서는 '고스트 레이어(Ghost Layer)'라는 가상의 격자 층을 두어, 블록 간에 데이터가 끊김 없이 흐르도록(연속성 유지) 설계했다.

3. 결과 분석 및 물리적 현상

실험 결과에서 가장 흥미로운 부분은 노즐 입구 압력 변화에 따른 유동의 변화이다.

• 부족 팽창(Under-expanded) 유동: 입구 압력이 높아질수록 노즐 출구 압력이 대기압보다 커진다. 이때 가스가 밖으로 나오면서 갑자기 사방으로 퍼지는데, 이때 발생하는 팽창파(Expansion Waves)를 시뮬레이션이 성공적으로 포착했다.
• 마하 수(Mach Number)의 증가: 노즐을 통과하며 유동이 가속되어 마하 1.5 이상의 초음속에 도달하며, 배기 영역에서 압력이 떨어짐에 따라 속도가 더욱 빨라지는 현상을 확인했다.
• 이론과의 비교 (Prandtl-Meyer 함수): 시뮬레이션에서 나타난 유동의 굴절 각도(Deflection Angle)를 이론적인 계산값과 비교했을 때, 약 10도~22도 사이의 분포를 보이며 이론적 예측 범위 내에 있음을 검증했다.

4. 이 연구의 한계점 (Critical Limitations)

논문은 스스로 몇 가지 한계점을 솔직하게 밝히고 있다.

• 벽면 오차: 실제보다 노즐 벽면 근처의 경계층(Boundary Layer) 두께를 두껍게 예측하는 경향이 있어, 마하 수 계산에서 약간의 오차가 발생했다.
• 과팽창(Over-expanded) 유동 해결 불가: 압력이 너무 낮아 유동이 벽면에서 떨어져 나가는 복잡한 현상은 현재의 수치 모델로는 계산이 불가능하다고 명시했다.

---

5. 기준 입구 조건 (Stagnation Conditions)

연구에서 설정을 위해 사용된 기준 압력과 온도는 다음과 같다.

• 정체 압력 ($P_0$): $1.2 \times 10^6 \, \text{Pa}$ (약 12기압)
• 정체 온도 ($T_0$): $2,000 \, \text{K}$
• 작동 유체: 공기 (비열비 $\gamma = 1.4$, 기체 상수 $R = 287$)

이 조건에서 계산된 노즐목(Throat, $M=1$)에서의 물리량은 다음과 같다.

• 속도 ($u_{CH}$): $818.38 \, \text{m/s}$
• 온도 ($T_{CH}$): $1,667 \, \text{K}$
• 압력 ($P_{CH}$): $6.340 \times 10^5 \, \text{Pa}$

---

6. 입구 압력 변화에 따른 팽창파 분석

논문에서 가장 중요한 데이터는 입구 압력($P$)을 조절했을 때, 노즐 출구 외부에서 발생하는 유동의 굴절 각도($\theta$)이다. 이론값(Prandtl-Meyer 함수)과 시뮬레이션 관측값을 비교한 수치는 다음과 같다.

입구 압력 조건	구분	출구 마하 수 ($M_e$)	하류 마하 수 ($M_2$)	수치적 예측 각도	실제 관측 각도
$0.75 P_0$	중심선 기준	1.71	1.99	$8.09^\circ$	$10.01^\circ$
$1.0 P_0$	중심선 기준	1.76	2.26	$11.94^\circ$	$16.65^\circ$
$1.7 P_0$	중심선 기준	1.85	2.69	$21.36^\circ$	$22.53^\circ$

• 해석: 입구 압력이 높아질수록($0.75 P_0 \rightarrow 1.7 P_0$) 노즐 밖으로 분사되는 가스의 속도(하류 마하 수)가 빨라지며, 유동이 꺾이는 각도 또한 약 $10^\circ$에서 $22.5^\circ$까지 크게 증가하는 것을 수치로 확인할 수 있다.

---

7. 위치별 유동 분포 데이터 (Figure 11)

중심선(Centerline)을 따라 측정한 거리($x$)별 데이터 변화는 다음과 같다.

• 압력 분포 ($P/P_{atm}$):
  • 노즐 입구($x=0$)에서 약 5배 높았던 압력이 노즐을 통과하며 급격히 하강한다.
  • 노즐 출구($x \approx 3.5 \times 10^{-5}$)를 지나면서 대기압($1.0$) 수준으로 수렴하기 위해 매끄럽게 감소한다.
• 마하 수 분포 ($M$):
  • 입구에서 $M=1$로 시작한 유동이 노즐 출구 근처에서 약 1.7 ~ 1.8에 도달한다.
  • 배기 영역(노즐 밖)으로 나가면서 팽창파의 영향으로 가속되어 최종적으로는 마하 2.0 ~ 2.7 사이의 초음속에 도달한다.

---

8. 오차 및 시뮬레이션 격자 데이터

• 격자 해상도: 각 블록당 $40 \times 61$ 개의 격자점을 사용하여 총 4개의 블록으로 시뮬레이션을 수행했다.
• 검증 결과: Olson의 기존 코드(FVM 방식)와 비교했을 때, 본 DNS 코드는 압력을 다소 높게 예측하는 경향을 보였으나, 팽창파의 각도 변화는 이론적 범위 내에 있어 물리적 타당성을 갖춘 것으로 판단되었다.

이처럼 구체적인 수치들은 이 시뮬레이션이 단순히 그림을 그리는 수준을 넘어, 실제 초음속 엔진 설계에서 발생하는 팽창 현상을 수치적으로 예측할 수 있음을 보여준다. <끝>