본문 바로가기
현장과 프로젝트/pyroCEA

NASA CEA로 이해하는 로켓 성능 계산: 동결 평형과 이동 평형

by 도서관경비원 2025. 3. 22.
반응형

NASA CEA란 무엇인가

NASA CEA(Chemical Equilibrium with Applications)는 연소 및 로켓 추진 문제를 분석하기 위해 NASA가 개발한 소프트웨어다. 포트란(Fortran)으로 작성된 이 고전적인 도구는 주어진 추진제 조건에서 연소실의 평형 조성과 온도를 계산하고, 로켓 설계에 필수적인 성능 수치들을 도출한다.

 

CEA가 계산하는 핵심 성능 지표는 세 가지다. 첫째는 특성 배기 속도(C*) 로, 추진제가 에너지를 얼마나 효율적으로 생성하는지를 나타내는 연소실 고유의 성능 지표다. 둘째는 추력 계수(Thrust Coefficient) 로, 노즐이 연소 압력을 추력으로 얼마나 효율적으로 변환하는지를 나타낸다. 셋째는 비추력(Specific Impulse, Isp) 으로, 단위 추진제 질량당 발생하는 추력을 의미하며 로켓 엔진의 연료 효율을 나타내는 가장 대표적인 지표다. 비추력은 m/s 단위로 제공되며, 일반적으로 사용하는 초(s) 단위로 변환하려면 중력 상수(9.81 m/s²)로 나누면 된다.


동결 평형(Frozen) 가정: 보수적인 접근

CEA에서 동결(Frozen) 옵션을 선택하면, 연소실에서 형성된 혼합물 조성이 노즐 전체를 통과하는 동안 변하지 않고 고정된다고 가정한다. 실제로는 노즐을 따라 온도와 압력이 급격히 변화하면서 혼합물의 화학적 조성도 함께 달라진다. 그러나 동결 가정은 이러한 변화를 무시하고 연소실 출구 조성이 노즐 끝까지 그대로 유지된다고 본다.

 

초음속 노즐 면적비를 68.8로 설정했을 때 동결 평형 조건의 계산 결과는 다음과 같다.

  • 연소실 특성 배기 속도(C*): 2,289.4 m/s
  • 노즐목 압력: 111.87 bar, 온도: 3,276.93 K
  • 노즐 출구 압력: 0.17823 bar, 온도: 980.79 K
  • 노즐 출구 추력 계수: 1.8450, 비추력: 4,224.0 m/s

동결 가정은 실제 성능을 약 1~4% 과소평가하는 경향이 있다. 노즐 내에서 실제로 일어나는 화학 반응을 통한 추가적인 에너지 방출을 반영하지 않기 때문이다.


이동 평형(Shifting Equilibrium) 가정: 낙관적인 접근

이동 평형(Shifting Equilibrium) 옵션은 반대 방향의 가정을 취한다. 노즐 내 화학 반응이 무한히 빠르게 일어난다고 가정하여, 노즐의 모든 위치에서 혼합물 조성이 해당 시점의 온도와 압력에 맞는 평형 상태를 즉각적으로 달성한다고 본다. 실제 화학 반응 속도는 유한하므로 이 가정은 현실보다 이상적인 조건을 전제하는 셈이다.

 

동일한 조건에서 이동 평형을 적용한 결과는 다음과 같다.

  • 연소실 특성 배기 속도(C*): 2,322.1 m/s
  • 노즐 출구 압력: 0.20607 bar
  • 노즐 출구 추력 계수: 1.8827, 비추력: 4,371.9 m/s

이동 평형 가정은 실제 성능을 약 1~4% 과대평가하는 경향이 있다. 반응 속도가 무한하다는 이상적인 전제로 인해 실제보다 더 많은 에너지가 회수되는 것으로 계산되기 때문이다.


두 가정의 비교와 실제 성능 추정

두 접근법의 결과를 나란히 놓고 보면 차이가 명확하다.

항목 동결 평형 (Frozen) 이동 평형 (Shifting)
특성 배기 속도 (C*) 2,289.4 m/s 2,322.1 m/s
노즐 출구 압력 0.17823 bar 0.20607 bar
추력 계수 1.8450 1.8827
비추력 (Isp) 4,224.0 m/s 4,371.9 m/s
오차 경향 실제보다 과소평가 실제보다 과대평가

 

동결 가정은 하한선, 이동 평형은 상한선을 제공하는 셈이다. 실제 로켓 엔진의 성능은 이 두 값 사이 어딘가에 위치한다. 실무에서는 두 결과의 차이 중 약 40% 지점을 실제 성능 추정치로 활용하는 방식이 일반적으로 사용된다.

 

비추력을 예로 들면, 두 값의 차이는 4,371.9 - 4,224.0 = 147.9 m/s이며, 이 차이의 40%인 약 59 m/s를 동결 평형 값에 더한 4,283 m/s 수준이 현실적인 추정치가 된다.


CEA가 보여주는 로켓 설계의 현실

CEA의 두 가지 계산 모드는 로켓 성능 분석에서 마주치는 근본적인 어려움을 잘 보여준다. 노즐 내부의 화학 반응은 매우 빠르고 복잡하게 진행되며, 이를 시간에 따라 정확히 모델링하는 것은 현재 기술로도 쉽지 않다. 동결과 이동 평형이라는 두 가지 극단적 가정은 각각 현실의 하한과 상한을 제시함으로써, 엔지니어들이 실제 성능의 범위를 합리적으로 추정할 수 있도록 돕는다.

 

결국 로켓 설계는 이처럼 불확실성을 인정하고, 그 범위 안에서 최선의 판단을 내리는 과정이다. CEA는 그 판단을 위한 신뢰할 수 있는 출발점을 제공한다.

 

동결 평형 (Frozen)

 

이동 평형 (Shifting)

 

FCEA2.exe
2.19MB
thermo.lib
0.58MB
trans.lib
0.03MB

 

 

 

ex1.inp
0.00MB
ex2.inp
0.00MB
ex3.inp
0.00MB
ex4.inp
0.00MB
ex5.inp
0.00MB
ex6.inp
0.00MB
ex7.inp
0.00MB

반응형