(한동안 토익이랑 오픽하고 한참을 방치 했더니 기억이 살짝 가물가물하지만, 기억을 더듬어내면서 복습할 겸 적는다.)
로켓을 설계하기 위해서는 연료(Fuel+Oxid)의 연소속도, 성능 등 '연료가 가진 기본적인 특성'을 알아야 한다. 하지만 연료의 특성은 Fuel과 Oxid의 종류, 조성비, 챔버 내 압력 등에 따라 변한다. 물론 실험으로 알아낼 수도 있고 정확한 로켓 설계를 위해서는 실험이 필요하지만, 기존 데이터가 아에 없거나 효율을 가늠하기 위한 기준 성능 데이터가 필요한 경우에는 이론으로 풀어낸 연료의 특성 값을 찾아야 한다.
뭐 당연히 KNSB같은 경우는 아마추어 로켓티어들이 많이 쓰는 연료이므로 이미 실험이나 이론 계산으로 도출된 연료의 특성을 구글링으로 쉽게 찾을 수 있다. 하지만 나중에 다른 연료를 쓰거나 특별한 경우가 있을 수 있으니 이러한 툴을 사용하는게 좋다고 생각해, 이 방법으로 진행을 했다.
내가 쓸 툴은 NASA CEA(https://cearun.grc.nasa.gov/) 여기에서 CEA는 Chemical Equilibrium Analysis인데, 말 그대로 이 프로그램은 화학적 조성을 계산할 수 있는 웹프로그램이다. 연소 유동장 내에서는 일반적으로 푸는 생성물 조성과 다르고, 각종 변수들이 존재하기 때문에 깁스 에너지 최소화법으로 이를 풀어내야 하는데, 이 NASA CEA가 해당 툴이다. 하지만 이 프로그램의 문제는 내가 원하는 화학 물질이 데이터베이스에 없다면, 조합을 할 수 없고 계산도 할 수 없다는 단점이 있다.
그래서 나는 NASA CEA 프로그램을 오프라인버전(cmd)으로 만든 새로운 프로그램(http://akrmys.com/public/cea/cea_install_cmd.html.en)을 찾았다. 이 프로그램은 기존에 없는 데이터베이스를 직접 넣어 계산할 수 있는 장점이 있다. 그래서 나는 KNSB니까 Fuel로 쓸 Sorbitol과 Oxid로 쓸 KNO3의 물질의 특성을 각각 데이터베이스에 입력하여 계산하였다.
<데이터베이스 입력값>
[Sorbitol]
C6 H14 O6
h : -1353.7 kJ/mol
t(K) : 298.15
[KNO3]
K1 N1 O3
h: -494.5 kJ/mol
t(K) : 298.15
사용법은 (추후 글 작성 예정~~)을 참고.
<설계값>
챔버 압력 : 35bar
챔버 압력 / 출구압력(대기압 1bar) : 35
산화제와 연료의 비율 : 65/35

띄리리리리릭...
계산을 하면....
짜잔~~~~!!!


결과값 1에 있는 것이 중요한데,
rho : 1.841 g/cm^3
T_0 : 1590.32 K
M : 39.720 kg/kmol
Cp : 2.2615 kJ/kg*K (M을 곱해서 변환, 89.8268kJ/kmol*K)
Cv : 81.5128 kJ/kmol*K (C_p - R)
k : 1.102(Cp/Cv)
이와 같이 연료의 특성을 알 수 있다.
또한 이 연료의 특성을 이용하여
Ae/At 6.3746
C* 910.5 m/s
CF 1.5386
C 1400.8 m/s
Isp 142.793 s (Isp = C/g)
와 같은 로켓에 쓰이는 수식값들도 얻을 수 있다.
물론 이전 글을
Ae/At
노즐목에서 마하수 1이 되는 값을 찾는 것이 가장 이상적인 노즐목 면적이라고 할 수 있다.
노즐의 내부를 단열, 가역과정으로 가정(등엔트로피)해서 에너지 보존 법칙, 질량보존법칙, 열역학 제2법칙 등으로 이거 저거 해보면 (https://youtu.be/wzK0_AO8cLQ 블로그에 수식 적는게 너무 힘들어서 잘 설명하신 분 링크 남깁니다. 이해점 ㅠㅠㅠ)
(공식 첨부 ㄱㄱ)
라는 최적노즐확장비 공식이 도출된다.
여기에서 위에서 구한 연료의 특성을 대입하여 역수를 취하면, 우리가 원하는 Ae/At 값이 나온다.
Isp(비추력)
추진제 단위중량을 사용하여 얻을 수 있는 총임펄스
이 문장을 식으로 나타내면
(공식)
추력이 시간에 대해 일정하면
(공식)
이렇게 나타낼 수 있다.
여기에서 Ve를 에너지 보존법칙을 이용해 비열비 형태로 나타내면
(ve 공식) 이고,
대입하면
(Isp 공식)
이 나온다.
여기에서 위에서 구한 연료의 특성을 대입하면, 우리가 원하는 Isp 값이 나온다.
C(유효배기속도)
노즐 출구에서 배기가스 속도를 균일하다고 가정할 때 이 배기가스의 속도를 이야기 한다.
(C 공식)
위 식에서 Isp = (공식)이라고 했으므로
C= Isp * g_0
여기에서 위에서 구한 Isp값을 대입하면, 우리가 원하는 C 값이 나온다.
C*과 CF도 중요하지만 추후 기초 계산에는 사용되지 않기 때문에 그냥 패스하도록 하겠다.
자세한 증명 과정과 C*, CF도 알고 싶다면(http://www.nakka-rocketry.net/th_nozz.html 와 http://www.nakka-rocketry.net/th_thrst.html 를 참고!)
'고체로켓' 카테고리의 다른 글
| [고체로켓 프로젝트 SESANG] #1 로켓의 설계 요소 정하기 (1) | 2022.11.29 |
|---|


