UNIVERSIDAD AUTÓNOMA DE CHIHUAHUA

Ingenierı́a Aeroespacial

Propulsión
Dr. Hermes Moreno Alvarez

Parcial 3

Corrales Acosta Evelyn Elena 286060

Soto Delgado Mario Alan 281836
Gil Beltran Cesar Omar 281889

Grupo: 9AE1

30 de Noviembre de 2017
1. Tobera de dos posiciones
1.1. Planteamiento del problema
El empuje más alto a una presión externa dada pa se obtiene cuando la tobera
está dimensionada para tener una presión de plano de salida pe que coincida con
la presión externa. Para los cohetes que deben escalar a través de los motores de
etapa completa), esta coincidencia no se puede lograr con cada presión, y se debe ha-
cer algún compromiso. Desde el más simple hasta el más complejo, estos podrı́an ser:

1.Una tobera de campana de geometrı́a fija (el enfoque más común). La tobera
está dimensionada para coincidir a alguna presión intermedia, y está sobre expan-
dida en el suelo. La elección precisa de la relación de área o relación de presión se
puede optimizar a través de simulaciones de trayectoria, y es común encontrar que
este óptimo implicarı́a la separación de flujo en el suelo. Como esto no es deseable
para un motor robusto, la elección es diseñar para una separación incipiente al mo-
mento del lanzamiento.

2.Una tobera de dos posiciones, como la que se muestra en el boceto: la parte ini-
cial del ascenso se realiza con una tobera de baja expansión relativamente baja, que
podrı́a diseñarse como en (a), pero en algún punto una tobera, que podrı́a diseñarse
como en (a), pero en algún punto se coloca una extensión de tobera en su lugar,
para proporcionar empuje adicional a mayores altitudes. El área de salida para esta
extensión y el punto de transición son parámetros de diseño a ser investigados en
este HWK.

3.Una tobera de goma que se extenderá gradualmente para que siempre se ajuste
a presión durante el ascenso. Esto obviamente no es práctico, pero el concepto es
útil como un ideal de rendimiento.

(a) Considere un cohete cuya presión en la cámara es de 70 atm., La temperatura

de su cámara es de 3300K, el gas tiene una relación de calor especı́fica.
La campana interna de la tobera de dos posiciones está diseñada, como se indicó
anteriormente, para pe1 /pe0 ∼= 0,4. Calcule el área de salida Ae1 , velocidad de salida
ue1 , empuje de vacı́o, FV 1 = ṁue1 + pe1 Ae1 y empuje de tierra, F0 = FV 1 − pa Ae1

(b) Sea t1 el ”tiempo de transición”, en el cual la campana exterior se coloca

en su lugar y tb . el tiempo de agotamiento de la etapa. La masa inicial es m0 , la
masa en la transición es m1 , y la masa de que se consume es mb0 . Muestre que las
velocidades obtenidas en el ascenso vertical antes y después de la transición están
dadas por
  Z t1   Z tb
FV 1 m0 pa (t) FV 2 m1 pa (t)
v1 = In −gt1 −Ae1 dtvb −v1 = In −g(tb −t1 )−Ae2 dt
ṁ m1 0 m(t) ṁ mb t1 m(t)

donde m(t) = m0 − ṁt Suponiendo ahora que el área de salida Ae2 , se ha seleccionado
(y por lo tanto la nueva velocidad de escape ue2 , y la presión de salida pe2 , también),

1
muestre por diferenciación que la velocidad final vb , se maximiza si la transición
ocurre cuando la presión atmosférica satisface.
FV 1 − FV 2
pa (t1 ) =
Ae2 − Ae2
y que si se cumple, el empuje no tiene discontinuidad cuando se agrega la extensión.
(c) Suponga que la nueva configuración de tobera también comienza a estar cerca
de la separación, es decir, pe2 , = 0,4pa (t). Esta condición, junto con la condición
de optimalidad anterior, determina pe2 oAe2 . Realice el cálculo para el caso bajo
consideración (iteración u otro método será necesario). Muestre en forma de tabla
los valores de empuje en el suelo, a la altura donde se produce la transición, y en
vacı́o, para la campana de geometrı́a fija, la campana de dos posiciones y la campana
de ”motor de goma”. Usando la aproximación pa (atm) = e−z/6,8km , trace los tres
perfiles de empuje contra z (hasta aproximadamente 60 km).

a) √
∗ RTc
c = (1)
Γ(γ)
 γ+1
 2(γ−1)
√ 2
γ (2)
γ+1
2,2
p 2 2·0,2
Γ 1,2 = 0,6485
2,2

J
8,314
R= mol · K = 461,88 J ≈ 461,9
J
kg mol · K mol · K
0,018
mol
Sustituyendo los valores para c∗ y ṁ obtenemos:
r
J
461,88 ∗ 3300K
mol · K m
c∗ = = 1903,77
0,6485 s

N
Pc At 70 atm ∗ 1,013 × 105 2 ∗ 0,1 m2 kg kg
ṁ = ∗ = m = 372,47 ≈ 372,5
c m s s
1903,77
s
Según el inner bell designado para la tobera de dos posciones:
Pe1
= 0,4
Pa0

2
Por lo tanto, el numero de Mach a la salida esta dado por:
 
70 1,2 − 1 2
= 1+ Me1
0,4 2
  0,2
70 1,2 0,2 2
=1+ M
0,4 2 e1
2 (2,365 − 1)2
Me1 =
0,2
s
(2,365 − 1)2
Me1 =
0,2
∴ Me1 = 3,694 ≈ 3,695
Ası́ pues, la razón de área se obtiene a partir de la siguiente ecuación:
γ+1
! 2(γ−1)
Ae1 1 1 + γ−1
2
M 2
e1
= γ+1 (4)
At Me1 2

Sustituyendo en la siguiente ecuación:

1,2+1
! 2(1,2−1)
1,2−1
Ae1 1 1+ 2
· 3,69462
= 1,2+1
At 3,6946 2

Ae1
= 18,233 ≈ 18,234
At
Ae1 = 18,234(0,1 m2 )
∴ Ae1 = 1,8234 m2
Recordando la ecuación del empuje de vacı́o
Fv1 = ṁue1 + Pe1 Ae1 (5)
Podemos obtener la velocidad de salida despejando la de la ecuación anterior:
v
u
u γ
"   γ−1 #
P e1
γ
ue1 = t2 RTc 1 − (6)
γ−1 Pc
v
u "   0,2 #
u 1,2 0,4 1,2
ue1 = t2 ∗ ∗ 461,9 ∗ 3300K 1 −
0,2 70
m
∴ ue1 = 3249,193 ≈ 3249,2
s
En la superficie
F0 = Fv1 − Pa0 Ae1 (7)

F0 = 1,284 × 106 − 1,01 × 105 ∗ 1,8234 = 1,0995 × 106 N

Fv1 = ṁue1 +Pe1 Aw1 = 372,5∗3421+0,4∗1,013×105 ∗1,8234 = 1,2842×106 N

3
b) Desarrollando
dv
m = F − mg = Fv − Pa Ae − mg (8)
dt
dv Fv Ae Pa
= −g− (9)
dt m m
m = m0 − ṁt
dm = −ṁdt (10)

Despejando dv de la ecuación 9 y sustituyendo dt de la ecuación 10 reexpre-

samos:
Fv dm Ae Pa
dv = − − gdt − dt (11)
ṁ m m
Representando la integral de la ecuación anterior nos da:
Z Z Z Z
fv dm A e Pa
dv = − gdt − dt
ṁ m m

Integrando para: t = 0 y t = t1 ; asi como Fv = Fv1 y Ae = Ae1 :

Z Z Z t=t1 Z t=t1
fv1 dm Pa
dv = − gdt − Ae1 dt
ṁ m t=0 t=0 m
  Z t1
fv1 m0 Pa (t)
v1 = In − gdt1 − Ae1 dt (12)
ṁ m1 0 m(t)
Ahora tomando los valores t = t1 y t = tb asi como Fv = Fv1 y Ae = Ae2 se
vuelve a hacer la integración:
Z Z Z t=tb Z t=tb
fv2 dm Pa
dv = − gdt − Ae2 dt
ṁ m t=1 t=1 m
  Z tb
Fv2 m1 Pa (t)
vb − v1 = In − g(tb − t1 ) − Ae2 dt (13)
ṁ mb t1 m(t)

Reemplazando la ecuación 12 en la ecuación 13v1 pasa sumando:

   Z tb 
Fv2 m1 Pa (t)
vb = In − g(tb − t1 ) − Ae2 dt + v1
ṁ mb t1 m(t)

Agrupando
    Z t1 Z tb
Fv1 m0 Fv2 m1 Pa (t) Pa (t)
vb = In + In − gtb − Ae1 dt − Ae2 dt
ṁ m1 ṁ mb 0 m(t) t1 m(t)
(14)
Sabiendo:
m1 = m0 − ṁt1 (15)

4
 dvb dm1
Para optimizar se establece = 0 y con la relación = −ṁ
dt1 dt1
Fv1 1 Fv2 1 Pa (t1 ) Pa (t1 )
− (−m) + (−m) − Ae1 + Ae2 =0
m m1 m m1 m1 m1
Fv2 − Fv1
Pa (t1 ) = para optimizar t1
Ae2 − Ae1
Si t1 es seleccionado, el empuje es:

Justo antes de la transición

Fv2 − Fv1
F1 (t1 − ε) = Fv1 − Ae1
Ae2 − Ae1
Fv1 Ae2 − Fv2 Ae1
F1 (t1 − ε) = (16)
Ae2 − Ae1
Justo después de la transición
Fv2 − Fv1
F2 (t1 + ε) = Fv2 − Ae2
Ae2 − Ae1
Fv1 Ae2 − Fv2 Ae1
F2 (t1 − ε) = (17)
Ae2 − Ae1
∴ F (t1 − ε) = F2 (t1 + ε) no hay descontinuidad en el empuje

c) Si determinamos ahora Pe2 = 0,4Pa (t1 ) (correspondiente a la separación en la

tobera) tenemos que:
Fv2 − Fv1
Pe2 = 0,4
Ae2 − Ae1
Teniendo: Fv2 = ṁue2 + pe2 Ae2

(ṁue2 + Pe2 Ae2 ) − Fv1

Pe2 = 0,4 (18)
Ae2 − Ae1
Aquı́, ṁ, Fv1 , y Ae1 ya las conocemos ( sacadas en el inciso a) ), y las otras
cantidades (ue2 , Pe2 , Ae2 ) dependen únicamente de un solo parámetro que
viene siendo Me2 . Por lo tanto la ecuación de arriba determina Me2 , y todas
las demás cantidades.

Considera la ecuación 18 y resuelve mediante prueba y error lo siguiente:

a) Supone un valor para Me2

 γ+1
γ − 1 2 2(γ − 1)

At  1 + 2 Me2 
b) Calcula Ae2 =
Me2
 γ+1 
2

5
 Tc √
c) Te2 = y ue2 = Me2 γRTe2
γ−1 2
1+ Me
2
  γ
Te2 γ − 1
d) Pe2 = Pc
Tc
e) Calcular el valor para el lado derecho de la ecuación 18 y comparela con
el valor de Pe2 en el lado izquierdo.

Los resultados calculados mediante matlab se muestran en las figuras.

La solución viene siendo Me2 = 4,815
Ae2
Encontrando: = 90,13
At
Ae2 = 9,013 m2
m
ue2 = 3, 575
s
Pe2 = 0,05236 atm
Fv2 = ṁue2 + Pe2 Ae2 = 1,379 × 106 N
Fv2 − Fv1
Pa (t1 ) = = 13,83 km
Ae2− Ae1 
1
z1 = 6,8 In = 13,83 km
0,1308

El empuje de transición de la tobera 1 a la 2 es:

F (t1 ) = Fv1 − Pa (t1 )Ae1 = 1,2570 × 106 N

Para la expansión ideal en la tobera, Pe = Pa (z) entonces:

v
u
u

  γ − 1
u γ Pa (z) γ 
ue = ut2 RTc 1 −

γ−1 Pc


F = ṁue
Para z = 0, z = z1 yz → ∞(Pa = 0) y para los 3 tipos de toberas la tabla a
continuación recolecta los valores de empuje

Empuje (MN)
Tobera z=0 z = z1 = 13,83 km z→∞
Geometrı́a Fija 1.0995 1.2570 1.2842
2 Posiciones 1.0995 1.2570 1.3795
Expansión ideal 1.1348 1.2836 1.5931

Tabla 1: Cálculos de empuje para las tres toberas.

Nota

6
El empuje incrementa con la altitud en todos los casos

La tobera de dos posiciones es equivalente al la tobera fija en z = 1, pero

claramente superior a la condición de vacı́o.

La expansión ideal de la tobera supera en rendimiento en todas las altitudes.

Observa como la curva de expasión ideal conecta con la tobera 1 en z ≈
1
6,23 km; que es donde éstas interceptan (6,8In = 6,23). También co-
  0,4
1
necta a la tobera 2 en z = 6,8In = 20,1 km donde igualmente
0,05236
interceptan.