Problemas de cohetes de

combustible líquido
220305-Combustió i propulsió de coets
Màster Universitari en Enginyeria Aeronàutica

Autores:
Jordi Ferrís Roig
Mireia Garcia Granell
Alberto Simón Félix
Raimon Fernández Temprado

Fecha:
09/06/2022

Escola Superior d’Enginyeries Industrial, Aeroespacial i

Audiovisual de Terrassa.
Universitat Politècnica de Catalunya
Problemas de combustible líquido Combustió i propulsió de coets

Índice
1. Problema térmico tobera ............................................................................................... 2
2. Problema .......................................................................................................................... 4
3. Problema ........................................................................................................................ 16

1
Problemas de combustible líquido Combustió i propulsió de coets

1. Problema térmico tobera

Calcula la densidad de flujo térmico a través de las paredes de una tobera, dado el
siguiente registro de temperaturas:

t [s] T [ºC]
0,0 -23,3
0,2 754,4
0,4 885,6
0,6 972,8
0,8 1027,2
1,0 1061,1
1,2 1088,9
1,4 1098,9
1,6 1116,1
1,8 1126,7
2,0 1158,3
2,2 1172,2
2,4 1178,9
2,6 1172,2
2,8 1147,8
3,0 1109,4
3,2 1047,2
3,4 1000,0
3,6 959,4
3,8 898,9
4,0 839,4

Datos:

- Conductividad térmica  = 0,073 cal/cm s ºC = 30,5432 · 10-3 KJ/m s ºC

- Difusividad térmica  = 0,086 cm2/s = 8,6 · 10-8 m2/s

Para calcular la densidad de flujo térmico a través de las paredes de una tobera entre dos
intervalos de tiempo determinado y teniendo un registro temporal de temperatura se puede
hacer utilizando la Ecuación (1):
𝑡𝑛
𝜆 𝑑Φ(𝜏) 𝑑𝜏 (1)
𝑞(𝑡𝑛 ) = ∫
√𝜋𝛼 0 𝑑𝜏 √𝑡 − 𝜏

La cual se puede resolver numéricamente si se desarrolla de la siguiente forma:

𝑛
𝜆 1 1 𝑇𝑖 − 𝑇𝑖−1
𝑞(𝑡𝑛 ) ≈ ∑ 𝑚𝑖 · 2 · [(𝑡𝑛 − 𝑡𝑖−1 )2 − (𝑡𝑛 − 𝑡𝑖 )2 ] ; (𝑚𝑖 = )
√𝜋𝛼 𝑖=1 𝑡𝑖 − 𝑡𝑖−1

Para agilizar el cálculo, se utiliza un doble bucle de Matlab y posteriormente, se suma

todos los valores de q y se multiplican por el intervalo de tiempo (0,2s), para obtener la
densidad de flujo térmico Q.

2
Problemas de combustible líquido Combustió i propulsió de coets

El doble bucle utilizado es el siguiente, siendo “t” el vector de tiempos, “T” el vector de
temperaturas, “sum_general” el sumatorio grande de la ecuación y “sum_t” otro
sumatorio para llevar el calculo de la parte entre corchetes. Los demás valores son los
mismos que en la ecuación original.

for i=2:size(t,2)

m(i-1)=(T(i)-T(i-1))/(t(i)-t(i-1));

for j=2:i
sum_t(j-1)=m(j-1)*2*((t(i)-t(j-1))^0.5-(t(i)-t(j))^0.5);
sum_general(i-1)=sum_general(i-1)+sum_t(j-1);
end

q(i)=(lambda/((pi*alfa)^0.5))*sum_general(i-1);

end

Q = sum(q)* 0.2;

Finalmente se obtiene un valor de densidad de flujo térmico de:

𝑄 = 23.9976,47 𝐾𝐽/𝑚2

3
Problemas de combustible líquido Combustió i propulsió de coets

2. Problema

Se quiere diseñar un motor cohete que funcione con biérgol (H2O2 – MMH) y utilice
un sistema de alimentación “pressure feed”.

Teniendo en cuenta lo siguiente:

- El reductor se utiliza para refrigerar la tobera y la cámara de combustión.

- Se utiliza un único sistema de presurización para presurizar los 2 ergoles.

T (kN) Tc (ºC) pc (MPa) r  pa (Pa) M (kg/mol)

30 2986 7 3,46 1,14 101300 0,02461

Oxidante Reductor
Inyectores P circuito  Inyectores P circuito 
Nº P (bar) Cd (bar) (kg/m ) Nº P (bar) Cd
3 (bar) (kg/m3)
50 6 0,4 3 1141 15 2,5 0,6 8 920

1) Propón un esquema del sistema de alimentación que abastezca desde el sistema

de presurización hasta la cámara de combustión.

Figura 1. Esquema del sistema de alimentación “pressure feed”.

4
Problemas de combustible líquido Combustió i propulsió de coets

2) Calcula:

a. El caudal másico total y el de cada ergol.

Para resolver este primer apartado, se parte de la Ecuación (2) que

relaciona el empuje con el caudal másico:

𝑇 = 𝑚̇𝑝 · 𝑣𝑒 + 𝐴𝑒 (𝑃𝑒 − 𝑃𝑎 ) (2)

Donde T es el empuje del cohete, 𝑚̇𝑝 es el caudal másico de combustible

total, 𝑣𝑒 es la velocidad de salida de los gases de la tobera, 𝐴𝑒 es el área de
salida de la tobera y 𝑃𝑒 y 𝑃𝑎 son las presiones de salida de la tobera y
ambiente, respectivamente.

Además, se considera que la tobera está adaptada (𝑃𝑒 = 𝑃𝑎 ) de modo que

el segundo término de la ecuación anterior se anula:

𝑇 = 𝑚̇𝑝 · 𝑣𝑒

La variable que se quiere despejar es el caudal másico y el empuje es un

dato del enunciado, de modo que es necesario obtener previamente el valor
de la velocidad de salida de los gases. Este se calcula aplicando la
Ecuación (3):

𝛾−1
𝛾 𝑅 𝑇𝐶 𝑝𝑒 𝛾
𝑣𝑒 = √2 [1 − ( ) ] (3)
𝛾−1 𝑀 𝑝𝑐

Sustituyendo por los valores correspondientes y teniendo en cuenta que la

presión de salida de los gases es igual a la presión ambiente al estar la
tobera adaptada:

1,14−1
1,14 8,3145 · (2986 + 273) 101300 1,14
𝑣𝑒 = √2 [1 − ( ) ]
1,14 − 1 0,02461 7 · 106

𝑣𝑒 = 2696,76 𝑚/𝑠

Una vez conocido el valor de la velocidad de salida de los gases de la

tobera, ya que se puede calcular el caudal másico total haciendo uso del
empuje:

𝑇 30 · 103
𝑇 = 𝑚̇𝑝 · 𝑣𝑒 → 𝑚̇𝑝 = =
𝑣𝑒 2696,76

𝑚̇𝑝 = 11,12 𝑘𝑔/𝑠

5
Problemas de combustible líquido Combustió i propulsió de coets

Por último, se conoce del enunciado el ratio de mezcla oxidante-reductor

(r), de modo que se puede aplicar el siguiente sistema de ecuaciones para
obtener el caudal másico del oxidante y del reductor:

𝑚̇𝑝 = 𝑚̇𝑜 + 𝑚̇𝑟 (4.1)

{ 𝑚̇𝑜
𝑟= = 3,46 (4.2)
𝑚̇𝑟

Haciendo los cambios necesarios y sustituyendo por los valores

correspondientes se llega a:

𝑚̇𝑝 11,12
𝑚̇𝑟 = =
1+𝑟 4,46

𝑚̇𝑟 = 2,49 𝑘𝑔/𝑠

𝑚̇𝑜 = 𝑚̇𝑝 − 𝑚̇𝑟 = 11,12 − 2,49

𝑚̇𝑜 = 8,63 𝑘𝑔/𝑠

b. El diámetro de los inyectores.

Para calcular el diámetro de los inyectores es necesario conocer el flujo

másico que debe entregar cada inyector. Del enunciado se conoce que hay
50 inyectores para el oxidante y 15 inyectores para el reductor, de modo
que el caudal másico que debe entregar cada uno se obtiene dividiendo el
caudal másico de cada ergol por el número de inyectores que tiene cada
uno.

𝑚̇𝑜 8,63
𝑚̇𝑜,𝑖 = = = 0,173 𝑘𝑔/𝑠
𝑁𝑖,𝑜 50

𝑚̇𝑟 2,49
𝑚̇𝑟,𝑖 = = = 0,166 𝑘𝑔/𝑠
𝑁𝑖,𝑟 15

Una vez conocido el caudal másico que entrega cada inyector, se puede
calcular el diámetro de cada uno aplicando la Ecuación (5) y teniendo en
cuenta el área es circular:

(5)
𝑚̇ = 𝐶𝐷 𝐴√2𝜌∆𝑃

𝐷2
𝐴=𝜋
4

Combinando ambas ecuaciones, se puede despejar la variable del

diámetro:

6
Problemas de combustible líquido Combustió i propulsió de coets

𝐷2 4 · 𝑚̇
𝑚̇ = 𝐶𝐷 𝜋 √2𝜌∆𝑃 → 𝐷 = √
4 𝐶𝐷 𝜋√2𝜌∆𝑃

Finalmente, se sustituye por los valores correspondientes de cada ergol;

4 · 0,173
𝐷𝑖,𝑜 = √ = 0,003858 𝑚
0,4 · 𝜋 · √2 · 1141 · 6 · 105

𝐷𝑖,𝑜 = 3,858 𝑚

4 · 0,166
𝐷𝑖,𝑟 = √ = 0,004052 𝑚
0,6 · 𝜋 · √2 · 920 · 2,5 · 105

𝐷𝑖,𝑟 = 4,052 𝑚

c. El volumen de los tanques y cuál será su geometría si el tiempo de

operación es de 65 segundos y no se puede superar un diámetro de
0,90 metros. En caso de que alguno de los tanques sea cilíndrico,
dimensionar también los fondos tipo Klöpper (DIN 28011).

El primer paso que se debe realizar para obtener el volumen de los tanques,
tanto del oxidante como del reductor, es calcular la masa total que estos
deben almacenar. Para ello se multiplica el caudal másico por el tiempo de
operación:

𝑚𝑜 = 𝑚̇𝑜 · 𝑡 = 8,63 · 65 = 560,95 𝑘𝑔

𝑚𝑟 = 𝑚̇𝑟 · 𝑡 = 2,49 · 65 = 161,85 𝑘𝑔

Con todo ello, se puede obtener el volumen de los tanques aplicando la

Ecuación (6) que relaciona la densidad, con la masa y el volumen:

𝑚 𝑚
𝜌= →𝑉= (6)
𝑉 𝜌

Sustituyendo por los valores correspondientes:

𝑚𝑜 560,95
𝑉𝑜 = =
𝜌𝑜 1141

𝑉𝑜 = 0,49 𝑚3

𝑚𝑟 161,85
𝑉𝑟 = =
𝜌𝑟 920

𝑉𝑟 = 0,18 𝑚3

7
Problemas de combustible líquido Combustió i propulsió de coets

En segundo lugar, se desea conocer la geometría de los tanques. Para ello

se ha de tener en cuenta que siempre es preferible un tanque esférico
debido a que puede soportar el doble de tensiones con el mismo espesor
que un tanque cilíndrico, con lo cual se reduce el peso. Sin embargo, está
la limitación del diámetro máximo, de modo que hay que ver si es posible
o no.

El volumen de la esfera viene dado por la Ecuación (7), donde se despeja

el valor del radio y se calcula el diámetro necesario para almacenar el
volumen de los tanques calculado anteriormente.

4 3 3 3·𝑉
(7)
𝑉𝑒𝑠𝑓𝑒𝑟𝑎 = 𝜋𝑟 → 𝑟 = √
3 4·𝜋

3 3 · 0,49
𝑟𝑜 = √ = 0,489 𝑚
4·𝜋

𝐷𝑜 = 2 · 𝑟𝑜 = 2 · 0,489 = 0,978 𝑚 > 0,9 𝑚

3 3 · 0,18
𝑟𝑟 = √ = 0,35 𝑚
4·𝜋

𝐷𝑟 = 2 · 𝑟𝑜 = 2 · 0,35 = 0,7 𝑚 < 0,9 𝑚

Con todo ello, se puede concluir que el tanque para almacenar el reductor
será un tanque esférico de 0,7 metros de diámetro, mientras que el tanque
del reductor no puede ser cilíndrico porque el diámetro necesario supera
el máximo permitido.

El espesor mínimo de las paredes del depósito esférico destinado a

almacenar el reductor viene dado por la Ecuación (8):

𝑃 · 2𝑅𝑖𝑛𝑡 (8)
𝑒𝑒 ≥
4 · 0,9 · 𝜎 · 𝑧 − 𝑃

Donde P es la presión requerida dentro del tanque, 2𝑅𝑖𝑛𝑡 es el diámetro

interior, 𝜎 es la tensión nominal del límite elástico y 𝑧 es el factor de
soldadura (0,7).

La presión requerida dentro del tanque será la presión de la cámara de

combustión más la suma de las pérdidas de presión al circuito e inyector:

𝑃 = 𝑃𝑐 + ∆𝑃 = 7 + 0,3 + 0,6 = 7,9 𝑀𝑃𝑎

Por otro lado, la tensión nominal del límite elástico depende del material.
En este caso se ha optado por un Aluminio 7039-T64 por tener la mejor

8
Problemas de combustible líquido Combustió i propulsió de coets

relación masa-precio respecto los materiales vistos en clase. Este aluminio

tiene una tensión nominal del límite elástico de 380 MPa.

Sustituyendo los valores correspondientes en la Ecuación (8):

7,9 · 106 · 0,7

𝑒≥ = 0,00582 𝑚 = 5,82 𝑚𝑚
4 · 0,9 · 380 · 106 · 0,7 − 7,9 · 106

Con todo ello, el diámetro exterior se puede obtener a partir del diámetro
interior y el espesor de las paredes:

𝐷𝑒𝑥𝑡 = 𝐷𝑖𝑛𝑡 + 2𝑒

𝐷𝑖𝑛𝑡 = 0,7 − 2 · 0,00582 = 0,7116 𝑚

En la Tabla 1 se recogen las dimensiones del tanque esférico necesario

para almacenar el reductor.

Tanque esférico reductor

Diámetro interior 0,7 m
Diámetro exterior 0,7116 m
Espesor paredes 5,82 mm
Tabla 1. Dimensiones del tanque esférico.

Por otro lado, el tanque el oxidante será cilíndrico con el fondo toroesférico
tipo Köppler. Dado que el diámetro máximo es de 0,9 metros, se fija este
valor como diámetro exterior (𝐷𝑒𝑥𝑡 = 0,9 𝑚)

El espesor mínimo de las paredes de la parte central cilíndrica del depósito

se obtiene aplicando la Ecuación (9):

𝑃 · 𝐷𝑒𝑥𝑡
𝑒≥ (9)
2 · 0,9 · 𝜎 · 𝑧 − 𝑃

Sustituyendo los valores correspondientes:

7,9 · 106 · 0,9

𝑒≥ = 0,0151 𝑚 = 15,1 𝑚𝑚
2 · 0,9 · 380 · 106 · 0,7 − 7,9 · 106

Con todo ello, el diámetro interior se puede obtener a partir del diámetro
exterior y el espesor de las paredes:

𝐷𝑖𝑛𝑡 = 𝐷𝑒𝑥𝑡 − 2𝑒

𝐷𝑖𝑛𝑡 = 0,9 − 2 · 0,0151 = 0,8698 𝑚

Una vez conocidas las dimensiones de la parte central cilíndrica del

tanque, se calcula el espesor de la parte inferior del tanque, el cual viene
dado por la Ecuación (10):

9
Problemas de combustible líquido Combustió i propulsió de coets

𝑒 = max (𝑒𝑠 , 𝑒𝑦 , 𝑒𝑏 ) (10)

Cuyas componentes se calculan a continuación. De la primera componente

son conocidos todos los valores:

𝑃 · 𝐷𝑒𝑥𝑡
𝑒𝑠 = (11)
2 · 0,9 · 𝜎 · 𝑧 − 0,5 · 𝑃

7,9 · 106 · 0,9

𝑒𝑠 = = 0,01497 𝑚
2 · 0,9 · 380 · 106 · 0,7 − 0,5 · 7,9 · 106

De la segunda componente 𝛽 es el coeficiente de corrección, el cual se

obtiene del gráfico de la Figura 2, y 𝑆𝑅𝑖𝑛𝑡 es del radio interior de la esfera.

3 1
𝛽𝑃 (4 𝑆𝑅𝑖𝑛𝑡 + 𝐷𝑖𝑛𝑡 ) (12)
𝑒𝑦 = 5
0,9 · 𝜎

Figura 2. Coeficiente de corrección.

Para obtener el valor de 𝛽 se debe calcular la Ecuación (13):

3 1 𝐷𝑖𝑛𝑡 𝑃 (13)
( + )
4 5 𝑆𝑅𝑖𝑛𝑡 𝜎𝑐á𝑙𝑐𝑢𝑙𝑜

3 1 0,8698 7,9 · 106

( + ) = 0,01961 ≈ 0,02
4 5 0,9 380 · 106

Haciendo uso del gráfico de la Figura 2, se obtiene un valor de 𝛽 = 0,75

Con todo ello, ya se pueden sustituir los valores correspondientes en la

Ecuación (12):

10
Problemas de combustible líquido Combustió i propulsió de coets

3 1
0,75 · 7,9 · 106 (4 0,9 + 0,8698)
𝑒𝑦 = 5 = 0,01471 𝑚
0,9 · 380 · 106

Por último, de la tercera componente son conocidos todos los valores:

1
1,5
0.825
3 1 𝑃 𝐷𝑖𝑛𝑡
𝑒𝑏 = ( 𝑆𝑅𝑖𝑛𝑡 + 𝐷𝑖𝑛𝑡 ) [ 𝜎 ( ) ] (14)
4 5 111 𝑟2
1,5

Donde:

1 1 1
𝑟2 = = =
10 · 𝑆𝑅𝑖𝑛𝑡 10 · 0,9 9
1
0,825 1,5
3 1 7,9 · 106 0,8698
𝑒𝑏 = ( 0,9 + 0,8698) [ ( ) ] = 0,01129 𝑚
4 5 380 · 106 1
111 9
1,5

Finalmente, el espesor de la parte inferior del tanque es:

𝑒 = max(0,01487; 0,01471; ,0,01129) = 0,01487 𝑚 = 14,97 𝑚𝑚

Llegados a este punto únicamente queda por calcular las alturas del tramo
cilíndrico (ℎ1) y del tramo esférico (ℎ2), las cuales se obtiene aplicando
las Ecuaciones (15) y (16) respectivamente:

ℎ1 = 3,5 · 𝑒 (15)

ℎ1 = 3,5 · 0,01497 = 0,0524 𝑚 = 52,4 𝑚𝑚

ℎ2 = 0,1935 · 𝑆𝑅𝑖𝑛𝑡 − 0,455 · 𝑒 (16)

ℎ2 = 0,1935 · 0,9 − 0,455 · 0,01497 = 0,1673 𝑚 = 167,3 𝑚𝑚

De esta forma se obtiene un volumen interno del fondo toroesférico de:

1
𝑉𝑓𝑜𝑛𝑑𝑜 = (𝐷 − 2𝑒)3 (17)
10 𝑒𝑥𝑡

1
𝑉𝑓𝑜𝑛𝑑𝑜 = (0,9 − 2 · 0,01497)3 = 0,0659 𝑚3
10

11
Problemas de combustible líquido Combustió i propulsió de coets

Este volumen se tiene que multiplicar por 2 ya que hay dos fondos
toroesférico:

𝑉𝑓𝑜𝑛𝑑𝑜 𝑡𝑜𝑡𝑎𝑙 = 2 · 𝑉𝑖𝑛𝑡 = 2 · 0,0659 = 0,132 𝑚3

Finalmente, sabiendo que el volumen necesario para almacenar la masa de

oxidante es 𝑉𝑜 = 0,49 m3 (calculado anteriormente) y el volumen de los
fondos toroesféricos, se puede calcular el volumen de la parte cilíndrica
central:

𝑉𝑐𝑖𝑙𝑖𝑛𝑑𝑟𝑜 = 𝑉𝑜 − 𝑉𝑓𝑜𝑛𝑑𝑜 𝑡𝑜𝑡𝑎𝑙

𝑉𝑐𝑖𝑙𝑖𝑛𝑑𝑟𝑜 = 0,49 − 0,132 = 0,358 𝑚3

Con todo ello, conociendo el volumen de la parte cilíndrica, se puede

obtener la altura de esta parte:

𝐷𝑖𝑛𝑡 2 𝑉𝑐𝑖𝑙𝑖𝑛𝑑𝑟𝑜 · 4
𝑉𝑐𝑖𝑙𝑖𝑛𝑑𝑟𝑜 = ℎ𝑐𝑖𝑙𝑖𝑛𝑑𝑟𝑜 𝜋 → ℎ𝑐𝑖𝑙𝑖𝑛𝑑𝑟𝑜 =
4 𝜋 · 𝐷𝑖𝑛𝑡 2

0,358 · 4
ℎ𝑐𝑖𝑙𝑖𝑛𝑑𝑟𝑜 = = 0,602 𝑚
𝜋 · 0,86982

En la Tabla 2 se recogen las dimensiones del tanque cilíndrico necesario

para almacenar el oxidante.

Tanque cilíndrico con fondo toroesférico oxidante

Diámetro exterior 0,9 m
Diámetro interior 0,8698 m
Espesor paredes cilindro 15,1 mm
Espesor paredes fondo 14,97 mm
Altura cilindro fondo 52,395 mm
Altura esfera fondo 167,3 mm
Altura cilindro central 0,602 m
Tabla 2. Dimensiones del tanque cilíndrico con fondo toroesférico.

d. Si la presurización del sistema durante de fase de ensayos en el banco

de pruebas se hace con depósitos de aire de 50 litros a 200 bar (a 15
ºC) i la temperatura ambiente se considera constante a 25 ºC.
¿Cuántos depósitos se necesitarán para cada ensayo? Se puede asumir
que el sistema está bastante aislado del exterior para considerar que
el proceso se realiza de forma isentrópica.

Para obtener el número de depósitos necesarios para el ensayo, en primer

lugar, de debe calcular la masa de aire necesaria, tanto para el tanque del
oxidante como del reductor, aplicando la Ecuación (18):

12
Problemas de combustible líquido Combustió i propulsió de coets

𝑃𝑝 · 𝑉𝑝 𝛾
𝑚𝑜 = ( ) (18)
𝑅 · 𝑇0 𝑃𝑓
1−𝑃
𝑜

Donde 𝑃𝑝 y 𝑉𝑝 son la presión y el volumen de los tanques. La presión de

cada tanque es la presión de la cámara de combustión más la suma de las
pérdidas de presión en los inyectores y los circuitos, mientras que los
volúmenes son datos conocidos del apartado c):

𝑃𝑝,𝑜 = 𝑃 + 𝛥𝑃𝑖,𝑜 + 𝛥𝑃𝑡,𝑜 = 7 + 0,6 + 0,3 = 7,9 𝑀𝑃𝑎

𝑃𝑝,𝑟 = 𝑃𝑐 + 𝛥𝑃𝑖,𝑟 + 𝛥𝑃𝑡,𝑟 = 7 + 0,8 + 0,25 = 8,05 𝑀𝑃𝑎

𝑉𝑝,𝑜 = 𝑉𝑡,𝑜 = 0,49 𝑚3

𝑉𝑡,𝑟 = 0,18 𝑚3

Además, se asume que 𝑃𝑓 = 𝑃𝑝

Por otro lado, es necesario calcular la 𝑃𝑜 ya que la temperatura ambiente

no es la misma a la cual se conoce la presión de los depósitos.

𝑃 200 · 298
= 𝑐𝑡𝑒 → 𝑃𝑜 = = 206,94 𝑏𝑎𝑟
𝑇 288

Con todo ello ya se pueden obtener las masas de reductor y oxidante

sabiendo que:

7,9 · 106 · 0,49 1,4

𝑚𝑜,𝑜 = ( ) = 102,49 𝑘𝑔
287 · 298 7,9 · 106
1−
206,94 · 105

8,05 · 106 · 0,18 1.4

𝑚𝑜,𝑟 = ( ) = 38,82 𝑘𝑔
287 · 298 8,05 · 106
1−
206,94 · 105

Por tanto, la masa total del gas presurizado es:

𝑚𝑜 = 𝑚𝑜,𝑜 + 𝑚𝑜,𝑟 = 102,49 + 38,82 = 141,31 𝑘𝑔

En segundo lugar, de debe calcular la masa de gas que hay en nada depósito:

𝑃𝑉 200 · 105 · 50 · 10−3

𝑚𝑑𝑒𝑝ó𝑠𝑖𝑡𝑜 = = = 12,01 𝑘𝑔
𝑅𝑇 287 · (15 + 273)

13
Problemas de combustible líquido Combustió i propulsió de coets

Finalmente, para obtener el número de depósitos necesarios basta con dividir

la masa de aire necesario para presurizarlos entre la masa de aire que contiene
cada deposito:

𝑚𝑜 141,31
𝑁𝑑𝑒𝑝ó𝑠𝑖𝑡𝑜𝑠 = = = 11,76 → 12
𝑚𝑑𝑒𝑝ó𝑠𝑖𝑡𝑜 12,01

𝑁𝑑𝑒𝑝ó𝑠𝑖𝑡𝑜𝑠 = 12

3) Dimensiona la tobera, teniendo en cuenta que las pérdidas por divergencia

admisibles sean del 2% del empuje que daba la tobera.

En el último apartado del problema se quiere dimensionar la tobera, es decir, se deben

obtener los valores de las áreas de la garganta y la sección de salida, así como su
longitud.

En este caso, dado que las perdidas por divergencia admisibles son del 2%, se sabe
que  = 1-0,2 = 0, 98 y partir de este valor se puede calcular el ángulo de expansión
aplicando la Ecuación (19):

𝜆 = 0,5 · (1 + cos(𝛼) (19)

𝛼 = acos(2𝜆 − 1) = acos(2 · 0,98 − 1) = 16,26 º

Por otro lado, se puede calcular el área de la garganta de la tobera aplicando la

Ecuación (20) que lo relaciona con el coeficiente de empuje, aunque primero hay que
calcular dicho coeficiente:

𝛾+1 𝛾−1
𝑇 2𝛾 2 2 𝛾−1 𝑃𝑒 𝛾
(20)
𝑐𝐹 = =√ ( ) (1 − ( ) )
𝑃𝑐 · 𝐴𝑔 𝛾−1 𝛾+1 𝑃𝑐

1,14+1 1,14−1
2 · 1,142 2 1,15−1 0,1013 1,14
𝑐𝐹 = √ ( ) (1 − ( ) ) = 1,636
1,14 − 1 1,14 + 1 7

Despejando el área de la garganta de la Ecuación (20) y suponiendo que es circular:

𝐹 30 · 106
𝐴𝑔 = = = 2,6196 · 10−3 𝑚2
𝑃𝑐 · 𝑐𝐹 7 · 106 · 1,636

4 · 𝐴𝑔
𝐷𝑔 = √ = 0,0578 𝑚
𝜋

14
Problemas de combustible líquido Combustió i propulsió de coets

𝐷𝑔 = 0,0578 𝑚

Por otro lado, calculando la relación de expansión mediante la Ecuación (21), se

puede obtener el área de salida:
1 1
2 1−𝛾 𝑃𝑐 𝛾
𝐴𝑒 (𝛾 + 1) (𝑃 )
𝑒 (21)
𝜀= =
𝐴𝑓 𝛾−1
√𝛾 + 1 (1 − ( ) 𝑃𝑒 𝛾
)
𝛾−1 𝑃𝑐

1 1
2 1,14−1 7 · 106 1,14
(1,14 + 1) ( 6 )
0,1013 · 10
𝜀= = 10,17
1,14−1
√1,14 + 1 (1 − (0,1013)
1,14
1,14 − 1 7 )

𝐴𝑒 = 𝐴𝑓 · 𝜀 = 2,6196 · 10−3 · 10,1748 = 0,02665 𝑚2

Del mismo modo que el área de la garganta, suponiendo que la sección es circular:

4 · 𝐴𝑒 4 · 0,02665
𝐷𝑒 = √ =√ = 0,1842 𝑚
𝜋 𝜋

𝐷𝑒 = 0,1842 𝑚
Con todo ello, mediante la Ecuación (22), se puede obtener la longitud de la tobera
cónica en función del radio de la garganta:

𝑅𝑔 (√𝜀 − 1) + 𝑅(sec 𝛼 − 1) (22)

𝐿𝑡𝑜𝑏𝑒𝑟𝑎 𝑐ó𝑛𝑖𝑐𝑎 (15º) =
tan 𝛼

1
0,0289(√10,1748 − 1) + 1,5 · 0,0289 ( − 1)
cos(15)
𝐿𝑡𝑜𝑏𝑒𝑟𝑎 𝑐ó𝑛𝑖𝑐𝑎 (15º) =
tan(15)

𝐿𝑡𝑜𝑏𝑒𝑟𝑎 𝑐ó𝑛𝑖𝑐𝑎 (15º) = 0,242 𝑚

Finalmente, en la Tabla 3 se resumen las dimensiones de la tobera calculadas.

Dimensiones de la tobera
Diámetro de la garganta 0,0578 m
Diámetro de salida 0,1842 m
Longitud de la tobera 0,223 m
Ángulo de expansión 16.26º
Tabla 3. Dimensiones de la tobera.

15
Problemas de combustible líquido Combustió i propulsió de coets

3. Problema

Se quiere diseñar un motor cohete que funcione con biérgol (H2 – O2) y utilice un
sistema de alimentación “pump feed” con generador de gases.

Teniendo en cuenta lo siguiente:

T (kN) Tc (ºC) pc (MPa) r  pa (Pa) M (kg/mol)

2 3352,32 7 6 1,1402 101300 0,0134604

Oxidante Reductor
Inyectores P circuito  Inyectores P circuito 
Nº  (mm) Cd (bar) (kg/m3) Nº  (mm) Cd (bar) (kg/m3)
500 7 0,4 30 1141 150 9 0,6 80 70,8

- La tobera está adaptada al nivel del mar.

- El reductor se utiliza para refrigerar la tobera y la cámara de combustión.
- Se utiliza un generador de gases alimentado por los mismos ergoles para accionar
la turbina con una relación de mezcla r = 5. Potencia específica disponible a la
turbina (9475 kJ/kg).
- Inyectores coaxiales
- La pc y la p del generador de gases coinciden con los de la cámara principal.
- La p hasta el generador de cases se considera negligible.

1) Propón un esquema del sistema de alimentación que abastezca desde el sistema

de presurización hasta la cámara de combustión.

Figura 3. Esquema del sistema de alimentación “pump feed”.

16
Problemas de combustible líquido Combustió i propulsió de coets

2) Calcula:

a. El caudal másico total y el de cada ergol.

Para el cálculo del caudal másico total se utiliza la Ecuación (2) que
relaciona el empuje con el caudal másico:

𝑇 = 𝑚̇𝑝 · 𝑣𝑒 + 𝐴𝑒 (𝑃𝑒 − 𝑃𝑎 )

Dónde T es el empuje y es dato del enunciado, 𝑣𝑒 es la velocidad de salida

de los gases de la tobera, 𝑝𝑒 es la presión a la salida y 𝑝𝑎 es la presión
ambiente. Dado que la tobera se encuentra adaptada 𝑃𝑎 = 𝑃𝑒 , siendo 𝑣𝑒 la
única incógnita.

La velocidad de salida de los gases se calcula aplicando la Ecuación (3):

𝛾−1
𝛾 𝑅 𝑇𝐶 𝑝𝑒 𝛾
𝑣𝑒 = √2 [1 − ( ) ]
𝛾−1 𝑀 𝑝𝑐

Sustituyendo los datos del enunciado queda un valor de 𝑣𝑒 :

1,1402−1
1,1402 8,3145 · (3352,32 + 273) 101300 1,1402
𝑣𝑒 = √2 [1 − ( ) ]
1,1402 − 1 0,014604 7 · 106

𝑣𝑒 = 3845,35 𝑚/𝑠

Pudiendo ya hallar el caudal másico total:

𝑇 2 · 106
𝑇 = 𝑚̇𝑝 · 𝑣𝑒 → 𝑚̇𝑝 = =
𝑣𝑒 3845,35

𝑚̇𝑝 = 520,108 𝑘𝑔/𝑠

Para el cálculo del caudal másico de cada ergol se tiene en cuenta la

relación oxidante reductor.

𝑚̇𝑜
𝑟=
𝑚̇𝑟

Además, sabiendo que:

𝑚̇𝑝 = 𝑚̇𝑜 + 𝑚̇𝑟

Combinando ambas ecuaciones se llega a los siguientes resultados:

17
Problemas de combustible líquido Combustió i propulsió de coets

𝑚̇𝑝 520,108
𝑚̇𝑟 = =
1+𝑟 7

𝑚̇𝑟 = 74,301 𝑘𝑔/𝑠

𝑚̇𝑜 = 𝑚̇𝑝 − 𝑚̇𝑟 = 520,108 − 74,301

𝑚̇𝑜 = 8445,807 𝑘𝑔/𝑠

b. Las pérdidas de carga en los inyectores (en bares).

Sabiendo que corresponden 500 inyectores para el oxidante y 150 para el

reductor y que el flujo másico de cada ergol se tiene que dividir entre los
inyectores se llega al resultado de:

445,807
𝑚̇𝑜𝑖 = = 0,8916 𝑘𝑔/𝑠
500

74,301
𝑚̇𝑟𝑖 = = 0,4953 𝑘𝑔/𝑠
150

Una vez conocido el caudal másico que entrega cada inyector, se puede
calcular el diámetro de cada uno aplicando la Ecuación (5), teniendo en
cuenta el área es circular:

𝑚̇ = 𝐶𝐷 𝐴√2𝜌∆𝑃

𝐷2
𝐴=𝜋
4

Combinando ambas ecuaciones, se puede despejar la variable de la caída

de presiones para el oxidante y el reductor:

𝐷𝑜𝑖 2
𝑚̇𝑜𝑖 = 𝐶𝐷 𝑜 · 𝜋 √2 · 𝜌𝑜 · ∆𝑃𝑖𝑛𝑦𝑜𝑥
4
2
1 4 · 𝑚̇𝑜𝑖
∆𝑃𝑖𝑛𝑦𝑜𝑥 = ( )
2 · 𝜌𝑜 𝐶𝐷 𝑜 · 𝜋 · 𝐷𝑜𝑖 2

2
1 4 · 0,8916
∆𝑃𝑖𝑛𝑦𝑜𝑥 = ( )
2 · 1141 0,4 · 𝜋 · (7 · 10−3)2

∆𝑃𝑖𝑛𝑦𝑜𝑥 = 14,7 𝑏𝑎𝑟

𝐷𝑟𝑖 2
𝑚̇𝑟𝑖 = 𝐶𝐷 𝑟 · 𝜋 √2 · 𝜌𝑟 · ∆𝑃𝑖𝑛𝑦𝑟𝑒𝑑
4

18
Problemas de combustible líquido Combustió i propulsió de coets

2
1 4 · 𝑚̇𝑟𝑖
∆𝑃𝑖𝑛𝑦𝑟𝑒𝑑 = ( )
2 · 𝜌𝑟 𝐶𝐷 𝑟 · 𝜋 · 𝐷𝑟𝑖 2

2
1 4 · 0,4953
∆𝑃𝑖𝑛𝑦𝑜𝑥 = ( )
2 · 70,8 0,6 · 𝜋 · (9 · 10−3 )2

∆𝑃𝑖𝑛𝑦𝑟𝑒𝑑 = 11,89 𝑏𝑎𝑟

c. El caudal másico total y el de cada ergol del generador de gases.

Para el cálculo del caudal másico total y de cada ergol del generador de
gases se tendrá en cuenta que el caudal que mueve el generador será tanto
como el oxidante 𝑚̇𝑜𝑥 , como el del reductor 𝑚̇𝑟𝑒𝑑 . El generador lo realiza
a una presión correspondiente a la suma de presiones de la cámara, 𝑃𝑐 , las
pérdidas del circuito, ∆𝑃𝑐𝑖𝑟𝑐𝑢𝑖𝑡𝑜 , y las pérdidas de los inyectores calculadas
en el apartado b). Por lo tanto, se utiliza la Ecuación (23):

𝑃𝑐 + ∆𝑃𝑐𝑖𝑟𝑐 𝑜𝑥 + ∆𝑃𝑖𝑛𝑗 𝑜𝑥 𝑃𝑐 + ∆𝑃𝑐𝑖𝑟𝑐 𝑟𝑒𝑑 + ∆𝑃𝑖𝑛𝑗 𝑟𝑒𝑑

𝑃𝑔 𝑚̇𝑔 = 𝑚̇𝑜𝑥 ( ) + 𝑚̇𝑟𝑒𝑑 ( ) (23)
𝜌𝑜𝑥 𝜌𝑟𝑒𝑑

Donde 𝑚̇𝑜𝑥 y 𝑚̇𝑟𝑒𝑑 se pueden expresar de la siguiente forma:

𝑚̇𝑜𝑥 = 𝑚̇𝑜𝑥 𝑔 + 𝑚̇𝑜𝑥 𝑐

𝑚̇𝑟𝑒𝑑 = 𝑚̇𝑟𝑒𝑑 𝑔 + 𝑚̇𝑟𝑒𝑑 𝑐

Y sabiendo que:

𝑚̇𝑔
𝑚̇𝑜𝑥 𝑔 = 1
𝑚̇𝑔 = 𝑚̇𝑜𝑥 𝑔 + 𝑚̇𝑟𝑒𝑑 𝑔 1+𝑟
𝑚̇𝑜𝑥 𝑔 → 𝑔
𝑟𝑔 = 𝑚̇ 𝑚̇𝑔
{ 𝑟𝑒𝑑 𝑔 } 𝑚̇𝑟𝑒𝑑 𝑔
= 1
1+𝑟
𝑔

Por lo tanto, sustituyendo en la primera expresión queda:

𝑚̇𝑔 𝑃𝑐 + ∆𝑃𝑐𝑖𝑟𝑐 𝑜𝑥 + ∆𝑃𝑖𝑛𝑗

𝑜𝑥
𝑃𝑔 𝑚̇𝑔 = (𝑚̇𝑜𝑥 𝑐 + )( )
1 𝜌 𝑜𝑥
1+𝑟
𝑔

𝑚̇𝑔 𝑃𝑐 + ∆𝑃𝑐𝑖𝑟𝑐 𝑟𝑒𝑑 + ∆𝑃𝑖𝑛𝑗 𝑟𝑒𝑑

+ (𝑚̇𝑟𝑒𝑑 𝑐 + )( )
1 𝜌 𝑟𝑒𝑑
1+𝑟
𝑔

19
Problemas de combustible líquido Combustió i propulsió de coets

Donde resolviendo se llega a los siguientes resultados:

𝑚̇𝑔 = 2,277 𝑘𝑔/𝑠

𝑚̇𝑟𝑒𝑑 = 0,379 𝑘𝑔/𝑠

𝑚̇𝑜𝑥 = 1,898 𝑘𝑔/𝑠

d. Diámetro y longitud de la cámara de combustión principal.

Teniendo en cuenta que la cámara de combustión presenta una forma

cilíndrica se debe seleccionar un valor de la longitud característica, 𝐿∗ , de
entre el intervalo para la combustión del biergol (𝐻2 + 𝑂2 ). Este intervalo
está comprendido entre 76 cm y 102 cm, y se ha seleccionado el valor
medio de 89 cm.

Una vez establecido este parámetro se desea conocer el diámetro de la

cámara de combustión principal. Para ello, mediante la Ecuación (24) se
puede hallar dicho valor ya que el resto de los parámetros, o bien han sido
calculados con anterioridad, o son datos del enunciado:

𝛾+1 𝛾−1
𝐹 2𝛾 2 2 𝛾−1 𝑝𝑒 𝛾 𝑝𝑒 − 𝑝𝑎
𝐶𝐹 = = √ ( ) [1 − ( ) ]+𝜀( ) (24)
𝑝𝑐 · 𝐴𝑔 𝛾−1 𝛾+1 𝑝𝑐 𝑝𝑐

Sabiendo que 𝐴𝑔 tiene un valor de 0,17465 𝑚2 y, por lo tanto, el diámetro

es 𝑑𝑔 = 47,149 𝑐𝑚. Posteriormente, se puede hallar también el volumen
de la cámara de combustión:

𝑉𝑐 = 𝐿∗ · 𝐴𝑔 = 155440,28 𝑐𝑚3

Con este valor ya se puede calcular el área de la cámara cilíndrica y por

tanto el diámetro de la cámara cilíndrica. Para ello primero se tiene que
realizar una aproximación de la longitud de la cámara, 𝐿𝑐 :
2 (𝑑
𝐿𝑐 ≈ 𝑒 0,029·𝑙𝑛 𝑔)+0,47·𝑙𝑛(𝑑𝑔)+1,94 = 65,47 𝑐𝑚

𝐴𝑐 = 2374,09𝑐𝑚2 → 𝐷𝑐′ = 54,979 𝑐𝑚

Correspondiendo este diámetro a aproximaciones y valores de referencia

experimentales. Para conocer el diámetro interno de la cámara de
combustión se itera varias veces la Ecuación (25) suponiendo un ángulo
de 45º:

24 ′
𝑑𝑔3 + · 𝑉 · 𝑡𝑎𝑛𝜃
𝐷𝑐 = √ 𝜋 𝑐 (25)
𝐷𝑐′ + 6𝐿𝑐 · 𝑡𝑎𝑛𝜃

20
Problemas de combustible líquido Combustió i propulsió de coets

Tras las iteraciones se llega al siguiente resultado:

𝐷𝑐 = 47,15𝑐𝑚

Siendo el volumen de la cámara de combustión correspondiente:

𝑉𝑐 = 114318,81 𝑐𝑚3

21