Màster Universitari en Enginyeria Aeronàutica

Combustió i Propulsió de Coets

Problemas entregables

Estudiantes:
Eduardo Simón Gil de Muro
Mª Magdalena Lirola Salas
Vı́ctor Thomas Ruiz
Zoe Higueras Arena
Índice
1. Problemas de cohetes de combustible sólido 3
1.1. Estimación de parámetros operacionales . . . . . . . . . . . . . . . . . . . . . 3
1.1.1. Regresión logarı́tmica de las propiedades en función de Pc . . . . . . . 3
1.1.2. Variación del área de combustión y relación de Klemmung en función
del tiempo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.1.3. Evolución de la presión en la cámara de combustión . . . . . . . . . . . 4
1.1.4. Estudio y gráfica de la velocidad de regresión (Vr ), el flujo másico del
combustible quemado (ṁp ) y la masa total consumida . . . . . . . . . . 5
1.1.5. Evolución de la relación de la presión de salida y la de cámara de
combustión . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.1.6. Evolución de la velocidad a la sección de salida de la tobera (ue ) . . . . 8
1.1.7. Evolución del empuje que proporciona el motor . . . . . . . . . . . . . 9
1.2. Estudio del efecto de aire y humedad en la combustión de PBAN 1516 . . . . 10
1.2.1. Estudio y gráfica de la temperatura de combustión (Tc) en función del
porcentaje de humedad (χH2 o ) . . . . . . . . . . . . . . . . . . . . . . . 11
1.2.2. Estudio y gráfica de la temperatura de combustión (Tc ) en función del
porcentaje de aire (χair ) . . . . . . . . . . . . . . . . . . . . . . . . . . 13
1.2.3. Conclusiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.3. Estudio comparativo del efecto de metales ligeros en la combustión PBAN y AP 15
1.3.1. Estudio de la densidad en función del % de metal ligero . . . . . . . . . 16
1.3.2. Estudio de la temperatura de combustión en condiciones estándar (1
bar) en función del % de metal ligero . . . . . . . . . . . . . . . . . . . 18
1.3.3. Estudio del impulso especı́fico a nivel del mar en función del % de metal
ligero . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
1.4. Estudio de la geometrı́a del grano de un cohete con alma troncocónica . . . . . 24
1.4.1. Determina el área de combustión, la masa de combustión, el volumen
del alma central y el ratio entre el volumen del grano y el volumen total
de la carcasa en el instante inicial. . . . . . . . . . . . . . . . . . . . . . 25

1
 1.4.2. Determina el área de combustión, la masa de combustión, el volumen
del alma central y el ratio entre el volumen del grano y el volumen total
de la carcasa en el instante en el que el diámetro interior en la parte
inferior del alma central logra el diámetro exterior del grano. . . . . . . 27
1.4.3. Estudia la evolución de la relación de Klemmung (KN ) en función de
la masa de combustible quemado . . . . . . . . . . . . . . . . . . . . . 28
1.4.4. Estudia la evolución de la presión de la cámara de combustión en fun-
ción de la relación de Klemmung. . . . . . . . . . . . . . . . . . . . . . 28
1.4.5. Estudia el instante final de la combustión . . . . . . . . . . . . . . . . . 29

2. ANEXO: Código de MATLAB empleado 31

2
 Report combustible sólido

1. Problemas de cohetes de combustible sólido

1.1. Estimación de parámetros operacionales

Se pretende obtener los parámetros operacionales de un motor cohete de combustible sóli-

do que consta únicamente de una etapa, conocido el combustible seleccionado, ası́ como la
geometrı́a del grano.

1.1.1. Regresión logarı́tmica de las propiedades en función de Pc

Para realizar el estudio se realiza un barrido en la presión de la cámara de combustión, Pc,

desde 0,1 M P a hasta 50 M P a, utilizando un incremento de 0,1 M P a.

(a) Temperatura de combustión (b) Relación de calores

(c) Masa molar media (d) Constante de los gases

Figura 1: Regresión obtenida para las distintas propiedades en función de ln(Pc)

De esta forma, se obtienen las siguientes ecuaciones para las respectivas regresiones:

3
 Report combustible sólido

Tc = 3075,5 + 58,624 · ln(Pc )

γ = 1,1528 + 0,0028 · ln(Pc )

(1)
M = 26,326 + 0,1810 · ln(Pc )

R = 315,86 + 2,1326 · ln(Pc )

1.1.2. Variación del área de combustión y relación de Klemmung en función del

tiempo

Utilizando el programa RPA, se obtienen las siguientes funciones del área de combustión y
relación de Klemmung en función del tiempo:

(a) Área de combustión (b) Relación de Klemmung

Figura 2: Variación de los parámetros dados en función del tiempo

Para la función de la relación de Klemmung, se obtiene un mı́nimo de KN,min = 480 y un

máximo de KN,max = 1440

1.1.3. Evolución de la presión en la cámara de combustión

Para un único instante dado, la presión de cámara del equilibrio se define de la siguiente
manera.
1
Pc(i+1) = (bST D · ρ · c ∗ (Pc(n) ) · KN ) 1−n (2)

Donde b, ρ y KN son datos para el instante dado, y c∗ se obtiene utilizando una presión de
cámara inicial (Pc(n) ). De esta forma se tiene un sistema iterativo, el cuál se detendrá cuando
la diferencia entre Pc(n+1) y Pc(n) sea de un orden inferior a las milésimas, como se puede

4
 Report combustible sólido

apreciar en el código de Matlab incluido en el anexo.

Figura 3: Evolución de Pc en el tiempo

Además, se obtienen los siguientes valores mı́nimo y máximo, ası́ como final:

Pceq,min = 12,619 M P a

Pceq,max = 90,812 M P a (3)

Pceq,f in = 43,757 M P a

Se puede observar que los datos no son del todo exactos con los dados en el enunciado del
problema (se tiene un error en torno a un 5 %). Este error, puesto que es bastante bajo, puede
ser debido a una diferencia en el barrido utilizado en el tiempo o a aproximaciones realizadas,
sin embargo, se puede concretar que el procedimiento está bien realizado.

1.1.4. Estudio y gráfica de la velocidad de regresión (Vr ), el flujo másico del

combustible quemado (ṁp ) y la masa total consumida

Para este estudio, de la misma manera que en el apartado anterior, se ha elaborado un código
Matlab que se puede encontrar en el Anexo.

Primeramente, se va a estudiar la velocidad de regresión (Vr ). Para ello, será necesaria la

5
 Report combustible sólido

siguiente expresión:
Vr = b(Tm )Pcn (4)

Por tanto, teniendo en cuenta esta ecuación, la gráfica obtenida de la velocidad de regresión
es la siguiente (Figura 4):

Figura 4: Velocidad de regresión en función del tiempo

Seguidamente, se procede al análidid del flujo másico del combustible quemado (ṁp ). Para
poder calcularlo, se va a utilizar la siguiente expresión:

ṁ = AC · ρp · vr (5)

Como se puede observar, presenta la velocidad de regresión máxima, con un valor de 0,0222
m/s en un t/tc de 0,33. Además, presenta su velocidad de regresión mı́nima de 0,0091 para
un t/tc de 0. Finalmente, para t/tc igual a 1 se obtiene una velocidad de regresión de 0,016.

Ası́ que, la gráfica obtenida del flujo másico del combustible quemado es la siguiente (Figura
5):

6
 Report combustible sólido

Figura 5: Flujo másico del combustible quemado en función del tiempo

Finalmente, se puede calcular el valor de la masa total consumida gracias a la función trapz
de Matlab. Por lo que se obtiene el siguiente resultado:

Mtotal = 216144, 17 kg (6)

1.1.5. Evolución de la relación de la presión de salida y la de cámara de com-

bustión

Conocida la relación de expansión, ϵ se puede obtener la relación de presiones utilizando la

siguiente expresión.
v
 1 
 γ−1  γ1 u "   γ−1 #
1 γ+1 Pe t γ + 1 1 − Pe
u γ
= (7)
ϵ 2 Pc γ−1 Pc

Puesto que γ es una variable que depende del instante en el que se estudie, la relación de
presiones también lo será. Con ello se obtiene una evolución como la siguiente:

7
 Report combustible sólido

Figura 6: Evolución de la relación de presiones, Pe /Pc

Como se puede apreciar en esta gráfica (Figura 6), la relación de presión de cámara y de
salida presenta un pico negativo en t/tc = 0,33. Anteriormente a este punto, es decir, de 0 a
0,33 la relación presenta un decrecimiento rápido. Aunque la variación de relación de presión
es minúscula, ya que sólo varı́a 3, 566 · 10−5 . Después, del valor de 0,33 hasta 1, la recta va
creciendo paulatinamente.

1.1.6. Evolución de la velocidad a la sección de salida de la tobera (ue )

En primer lugar, para poder realizar el estudio de la evolución de la ue , es necesario presentar

la expresión con la que se calcula:
v
u
u γ
"   γ−1 #
P e
γ
ue = t2 RTc 1 − (8)
γ−1 Pc

Por ende, aplicando esta expresión, se puede elaborar la siguiente gráfica (Figura 7):

8
 Report combustible sólido

Figura 7: Velocidad a la sección de salida de la tobera en función del tiempo

Como se puede observar, la ue presenta la forma opuesta a la gràfica de la relación de presión

de cámara y de salida (Figura 6) justo en el mismo valor de t/tc . Ası́ que, la velocidad de la
salida de la tobera presenta un pico positivo en t/tc = 0,33 teniendo un crecimiento rápido.
Aunque, del mismo modo que la relación de presiones de cámara y salida, la variación es
mı́nima, ya que sólo se observa una variación de 7,44 m/s de diferencia. Seguidamente, desde
t/tc = 0,33 hasta el final, la ue posee un decrecimiento no tan acentuado como el crecimiento
de 0 a 0,33.

1.1.7. Evolución del empuje que proporciona el motor

El empuje del motor se puede obtener mediante la siguiente expresión:

F = ṁ · ue + Ae · (Pe − Pa ) (9)

Donde ue es la velocidad de salida obteinda en el apartado anterior, y la Pa la presión

atmosférica a nivel del mar, Pa = 101325 M P a

9
 Report combustible sólido

Figura 8: Evolución del empuje proporcionado del motor

Además se obtienen los siguientes datos de empuje inicial, final y máximo:

Fini = 474 kN

Ff in = 1952 kN (10)

Fmax = 4182 kN

Se puede observar que los datos difieren un poco de los dados en el enunciado. Este hecho es
probablemente debido al solver utilizado, que pueda estar generando un error suficiente en los
cálculos que tras varias operaciones aumente lo suficiente, sin embargo, se puede comprobar
que los resultados obtenidos se mantienen dentro del mismo orden de magnitud.

1.2. Estudio del efecto de aire y humedad en la combustión de

PBAN 1516

La presencia de impurezas y sustancias en la matriz de los combustibles sólidos pueden

modificar significativamente las caracterı́sticas de la combustión de los combustibles sólidos,
ası́ como en la vida útil del grano.

Cuando se habla de impurezas se refiere a la presencia de humedad, aire (debido al grado

de compactación del grano) o la presencia de otras sustancias como metales, óxidos, aceites,

10
 Report combustible sólido

grasas, etc.

En esta ocasión, se pretende estudiar el efecto de la presencia de humedad y aire en la matriz

del combustible sólido PBAN 1516 (15 % de aluminio, 16 % de PBAN y 69 % de perclorato
amónico).

Primeramente, se tienen que estipular las hipótesis con las que se va a trabajar para realizar
el estudio. Dichas hipótesis son las siguientes:

La combustión da lugar en condiciones estándar: 1 bar y 21 ºC.

La humedad absorbida por el grano será en fase lı́quida y estará en el rango 0 % ÷ 5 %

w/w.

El aire atrapado en el interior del grano se encontrará en una proporción de 0 % ÷ 1 %

w/w.

La proporción entre los tres componentes principales se mantendrá constante:

16
χP BAN = (1 − χH2 O − χair ) [ %] (11)
100

15
χAl = (1 − χH2 O − χair ) [ %] (12)
100
69
χAp = (1 − χH2 O − χair ) [ %] (13)
100

Para realizar este estudio, se ha utilizado el programa RPA para obtener los datos y paráme-
tros necesarios para el análisis del estudio de aire y humedad en la combustión de PBAN
1516.

1.2.1. Estudio y gráfica de la temperatura de combustión (Tc) en función del

porcentaje de humedad (χH2 o )

Para este caso, al traterse del estudio en función del porcentaje de humedad, se van a compren-
der los valores entre 0 % y 5 % para cumplir el segundo punto de las hipótesis manifestadas
anteriormente.

Primeramente, se van a presentar los valores de las tres ecuaciones presentadas anteriormente
(Equaciones 11, 12 y 13) en la siguiente tabla (Tabla 1) teniendo en cuenta el porcentaje de

11
 Report combustible sólido

humedad. Por lo tanto,

Tabla 1: Valores de la proporción entre los tres componentes principales por cada
porcentaje de humedad

Porcentaje de H2 O ( %) χAl χP BAN χAp

0 0,1500 0,1600 0,6900
1 0,1485 0,1584 0,6831
2 0,1470 0,1568 0,6762
3 0,1455 0,1552 0,6693
4 0,1440 0,1536 0,6624
5 0,1425 0,1520 0,6555

Seguidamente, con los datos obtenidos con el RPA, es posible graficar la temperatura de
combustión en función del porcentaje de humedad. Por ende, ası́ como se puede observar en
la siguiente figura (Figura 9), la gráfica queda de la siguiente manera:

Figura 9: Temperatura de combustión en función del porcentaje de humedad

Gracias a la recta de esta gráfica (Figura 9), se puede determinar la ecuación de la recta:

Tc = −1929, 7 · χH2 O + 2936, 9 (14)

12
 Report combustible sólido

Además, por otra parte, se puede obtener el valor de la temperatura de combustión por cada
porcentaje de humedad. Estos resultados se presentan en la siguiente tabla (Tabla 2):

Tabla 2: Temperatura de combustión por cada porcentaje de humedad

Porcentaje de H2 O ( %) Tc (K)
0 2935,2956
1 2917,9564
2 2899,5806
3 2880,2376
4 2859,9855
5 2838,8692

1.2.2. Estudio y gráfica de la temperatura de combustión (Tc ) en función del

porcentaje de aire (χair )

Para el estudio de la temperatura de combustión en funcion del porcentaje de aire, se va a

seguir el mismo procedimiento que para el estudio anterior del porcentaje de humedad.

Sin embargo, en este caso, al tratarse de porcentaje en humedad, los datos estudiados se van
a comprender entre 0 % y 1 %.

Por consiguiente, se pueden presentar los valores de la proporción entre los tres componentes
principales por cada porcentaje de aire:
Tabla 3: Valores de la proporción entre los tres componentes principales por cada
porcentaje de aire

Porcentaje de aire ( %) χAl χP BAN χAp

0,0 0,15000 0,16000 0,69000
0,2 0,14970 0,15968 0,68862
0,4 0,14940 0,15936 0,68724
0,6 0,14910 0,15904 0,68586
0,8 0,14880 0,15872 0,68448
1,0 0,14850 0,15840 0,68310

Como en el análisis anterior, gracias a los datos cobtenidos con el RPA, se puede trazar la

13
 Report combustible sólido

Figura 10: Temperatura de combustión en función del porcentaje de aire

gráfica de la temperatura de combustión en función del porcentaje de aire. Por lo tanto, la

gráfica queda de la siguiente forma (Figura):

Gracias a la recta de esta gráfica (Figura 10), se puede determinar la ecuación de la recta:

Tc = 83, 37 · χair + 2935, 30 (15)

Como en el caso anterior del porcentaje de la humedad, gracias a esta gráfica, es posible
determinar los valores de las temperaturas por cada porcentaje de aire. Estos resultados se
presentan en la siguiente tabla (Tabla 4):

Tabla 4: Temperatura de combustión por cada porcentaje de aire

Porcentaje de aire ( %) Tc (K)

0,0 2935,2956
0,2 2935,4657
0,4 2935,6341
0,6 2935,8009
0,8 2935,9659
1,0 2936,1293

14
 Report combustible sólido

1.2.3. Conclusiones

Primeramente, como se puede observar en la gráfica de la Figura 9 y en su respectiva recta

(Ecuación 14), la temperatura de combustión en función del porcentaje de humedad presenta
una pendiente negativa; es decir, es una recta decreciente. En cambio, en la gráfica de la
Figura 10 y en su respectiva recta (Equación 15), la temperatura de combustión en función
del porcentaje de aire presenta una pendiente positiva; es decir, es una recta creciente.

Además, también se puede observar como tanto para el estudio de la humedad como para el
del aire, ambos empiezan en un 0 % con una temperatura de combustión de 2935,30 K. En
cambio, es interesante remarcar que, en el caso de 1 % de humedad (el segundo valor estudiado
del rango de porcentajes), la temperatura de combustión tiene un valor de 2917,96 K. Mientras
que para el 1 % de aire (último valor estudiado del rango de porcentajes), la temperatura de
combustión es de 2936,13 K. De modo que, se puede concluir que en el caso del porcentaje
de la humedad, la temperatura difiere mucho más que en el caso del porcentaje del aire.
De hecho, en cuanto al porcentaje de aire, la temperatura sólo varı́a 0,83 K (aumentando
su valor) del 0 % al 1 %. Por el contrario, en el estudio de porcentaje de humedad, su valor
disminuye 17,34 K en tan sólo 1 % (de 0 % a 1 %).

1.3. Estudio comparativo del efecto de metales ligeros en la com-

bustión PBAN y AP

Los combustibles de tipo compuesto están constituidos mayoritariamente por cristales oxi-
dantes - como el perclorato de amoniaco N H4 ClO4 - en una matriz, que a su vez, hace de
combustible (blinder) - como el derivado del polibutadieno, el PBAN -. A esta mezcla se le
añaden aditivos minoritarios como los metales ligeros, necesarios para completar la reacción
REDOX exotérmica.

En este estudio se analizará el efecto de variar la proporción másica de metales ligeros en un

grano formado por una proporción constante del 16 % de PBAN y perclorato de amoniaco
(AP). Los metales ligeros que se consideran son:
kg
Magnesio (Mg), con una densidad de ρM g = 1738 m3

kg
Aluminio (Al), con una densidad de ρAl = 2700 m3

15
 Report combustible sólido

kg
Berilio (Be), con una densidad de ρBe = 1848 m3

kg
Titanio (Ti), con una densidad de ρT i = 4540 m3

kg
Zirconio (Zr), con una densidad de ρZr = 6511 m3

Se ha utilizado el programa RPA para extraer todos los datos necesarios para realizar los
cálculos y operaciones deseadas.

1.3.1. Estudio de la densidad en función del % de metal ligero

En este apartado se llevará a cabo un estudio de la variación de la densidad del combustible

en función del porcentaje de metal ligero que se le añada.

Para ello se usará la siguiente fórmula:

xmet xP BAN xAP −1

ρ=( + + ) (16)
ρmet ρP BAN ρAP

De esta ecuación se conocen todos los valores, la mayorı́a son constantes, excepto xmet que
será el que se irá variando de 0 % a 84 %, ya que este último es el máximo dado que xP BAN
permanecerá constante a 16 % y la suma de ambos no puede superar el 100 % lógicamente.
Por otro lado, xAP será la resta entre xmet y xP BAN .

Una vez el programa ha proporcionado todos los valores necesarios para resolver la ecuación,
se ha procedido a resolver-la y graficar los resultados con Matlab.

Los resultados obtenidos se muestran representados a continuación en la siguiente gráfica

(Figura 11):

16
 Report combustible sólido

Figura 11: Densidad del combustible en función del % de metal ligero

En la Figura 11 se puede ver como en función de si el metal añadido tiene una densidad inferior
o superior al perclorato de amoniaco, la densidad del combustible aumenta o disminuye. En
los casos del Magnesio y el Berilio, cuya densidad es inferior a la del perclorato de amoniaco
a medida que se va aumentando el % de estos en la mezcla, la densidad del combustible
disminuye. Mientras que en los casos del Aluminio, Titanio y Zirconio pasa lo contrario, ya
que la densidad de estos es mayor a la del perclorato de amoniaco.

Cuando el % de metal ligero añadido es del 19 % la densidad del combustible obtenida es la

siguiente.

Tabla 5: Densidad del combustible con un 19 % de metal ligero en la mezcla

Metal ligero ρ
Magnesio 1629,7406
Aluminio 1740,2082
Berilio 1647,2094
Titanio 1831,0869
Zirconio 1874,5728

17
 Report combustible sólido

1.3.2. Estudio de la temperatura de combustión en condiciones estándar (1 bar)

en función del % de metal ligero

En este apartado se llevará a cabo el estudio de la Temperatura de combustión en condiciones

estándard en función del % de metal ligero. En este caso se ha establecido un rango de análisis
de 0 a 55 % de metal ligero, con variaciones de 5 %.

Los resultados obtenidos se muestran en la tabla siguiente (Tabla 6):

Tabla 6: Temperatura de combustión en función del % de metal ligero

Porcentaje de metal ligero ( %) Magnesio Aluminio Berilio Titanio Zirconio

0 2647,2774 2647,2774 2647,2774 2647,2774 2647,2774
5 2667,0833 2759,4632 2867,1947 2682,1615 2684,1658
10 2710,6753 2858,8926 3017,0092 2718,205 2720,8826
15 2729,973 2935,2956 2980,2855 2752,6718 2757,7401
20 2725,975 2948,1863 2562,3345 2760,9262 2794,9824
25 2694,3553 2843,3197 2537,3594 2783,3495 2832,7952
30 2614,7994 2547,4894 2544,6747 2691,2943 2871,318
35 2382,8215 2315,3363 2549,8657 2588,6264 2859,9257

Como en el apartado anterior, mediante un programa de Matlab se han graficado los resulta-
dos para verlos de forma más visual y para, además, poder apreciar mejor la evolución y el
comportamiento de la temperatura de combustión por cada variación de porcentaje de cada
metal. Ası́ que, en la siguiente gráfica se pueden apreciar estos valores (Figura 12).

18
 Report combustible sólido

Figura 12: Temperatura de combustión del combustible en función del % de metal ligero

Se ha calculado también la temperatura para un porcentaje de metal ligero del 19 % y los

resultados son los siguientes:

Tabla 7: Temperatura de combustión del combustible con un 19 % de metal ligero en la

mezcla

Metal ligero Tc (K)

Magnesio 2728,7259
Aluminio 2954,1455
Berilio 2705,3982
Titanio 2755,5678
Zirconio 2787,4938

En cuanto a las temperaturas máximas, como se puede apreciar en la gráfica de la Figura

12, serı́an las siguientes:

19
 Report combustible sólido

Tabla 8: Temperatura de combustión máxima y % de metal ligero al que pertenece

Metal ligero Tcmax % metal

Magnesio 2729,973 15
Aluminio 2954,1455 19
Berilio 3017,0092 10
Titanio 2783,3495 25
Zirconio 2871,318 30

Por otra parte, también es necesario calcular la presencia de HCl (χHCl ). Por lo tanto, se ha
seguido el mismo procedimiento que para analizar la temperatura de combustión. Es decir,
se ha variado el porcentaje entre el AP y el metal para observar la presencia de HCl por cada
variación. Ası́ que, con este método, se han obtenido los siguientes resultados gracias al RPA
(Tabla 9):

Tabla 9: Presencia del HCl (χHCl ) en función del % de metal ligero

Porcentaje del metal liguero ( %) Magnesio Aluminio Berilio Titanio Zirconio

0 0,2375539 0,2375539 0,2375539 0,2375539 0,2375539
5 0,1627465 0,2065619 0,2011984 0,223974 0,223163
10 0,1405014 0,1759676 0,1659169 0,2099896 0,2086675
15 0,1217233 0,1423029 0,0619571 0,1958208 0,194089
19 0,1082925 0,1122324 0,013121 0,1855288 0,1794508
20 0,1050905 0,1046115 0,0066656 0,1826851 0,1794508
25 0,0897049 0,0691168 0,0052255 0,1676035 0,1647793
30 0,0748125 0,0415253 0,0046977 0,1452191 0,1501032
35 0,0594334 0,0227266 0,0041373 0,102277 0,1372839

Como en el caso de la temperatura de combustión, todos estos valores de la presencia de

HCl, se han expresado en una gráfica (Figura 13) para poder observar el comportamiento
por cada porcentaje de cada metal.

20
 Report combustible sólido

Figura 13: Presencia de HCl (χHCl ) en función del % de metal ligero

Primeramente, en todos los metales ligeros, la presencia de HCl disminuye.

Como se puede observar, en el caso del Berilio, la presencia del HCl sufre un decrecimiento
muy rápido a medida que su porcentaje aumenta. Por el contrario, en los otros metales ligeros,
decrece más paulatinamente. Sobretodo en el caso del Titanio y el Zirconio.

En el caso del Manganesio, sufre una caı́da acentuada en los primeros porcentajes. Sin em-
bargo, despues va decreciendo más suavemente.

1.3.3. Estudio del impulso especı́fico a nivel del mar en función del % de metal
ligero

Para este apartado se considera un cohete que trabaja a una presión en la cámara de com-
bustión de 2 MPa y que tiene una relación de expansión de 15.

En este apartado se repetirá el proceso del apartado anterior, pero teniendo en cuenta el
valor del impulso especı́fico en lugar de la temperatura de combustión.

Ası́ pues, los resultados son los siguientes:

21
 Report combustible sólido

Tabla 10: Impulso especı́fico en función del % de metal ligero

Porcentaje de metal ligero ( %) Magnesio Aluminio Berilio Titanio Zirconio

0 156,19 2156,19 156,19 156,19 156,19
5 160,27 162,47 179,97 157,04 155,43
10 163,98 167,82 189,19 157,61 154,31
15 167,65 172,16 190,81 157,73 153,37
20 170,73 174,56 181,28 157,92 160,35
25 170,73 166,65 174,86 156,11 163,17
30 162,15 161,56 168,93 152,93 161,73
35 156,02 156,73 162,36 150,16 160,81

Para este caso también se ha generado una gráfica para coparar los resultados.

Figura 14: Impulso especı́fico del combustible en función del % de metal ligero

Otra vez se han extraido los valores para un porcentaje de metal ligero del 19 %.

22
 Report combustible sólido

Tabla 11: Temperatura de combustión del combustible con un 19 % de metal ligero en la

mezcla

Metal ligero Tc (K)

Magnesio 170,27
Aluminio 174,64
Berilio 185,85
Titanio 157,27
Zirconio 158,88

Como se puede observar en los datos mostrados el impulso especı́fico máximo serı́a:

Tabla 12: Temperatura de combustión máxima y % de metal ligero al que pertenece

Metal ligero Iespmax % metal

Magnesio 170,73 20
Aluminio 174,64 19
Berilio 190,81 15
Titanio 157,92 20
Zirconio 163,17 25

En este caso podemos ver como al añadir metales ligeros de menor densidad, como el Mag-
nesio, Berilio y Aluminio, el impulso especı́fico aumenta notablemente, sobre todo en el caso
del Berilio, que añadiendo un pequeño porcentaje a la mezcla ya se experimenta un aumento
considerable.

Finalmente la conclusión que se extrae viendo las 3 gráficas obtenidas en que a mayor densidad
del grano, como en el caso del Titanio y el Zirconio, menor incremento del impulso especı́fico.

En cuanto a la relación del impulso especı́fico y la temperatura de combustión, en las gráficas

podemos ver cierto parecido, y que para una temperatura de combustión más elevada se
obtiene un impulso especı́fico más elevado también.

23
 Report combustible sólido

1.4. Estudio de la geometrı́a del grano de un cohete con alma

troncocónica

Se pretende diseñar un cohete de combustible sólido con un grano de geometrı́a cilı́ndrica y

alma central troncocónica. Los datos que se muestran a continuación son conocidos en relación
con la geometrı́a de la carcasa, geometrı́a de la tobera, geometrı́a del grano y propiedades
del combustible.

Figura 15: Esquema de la configuración del motor cohete formado por la carcasa cilı́ndrica,
el grano con alma troncocónica, el ignitor en la culata y la tobera convergente-divergente

24
 Report combustible sólido

En primer lugar, se deben estipular las hipótesis con las que se va a trabajar. Dichas hipótesis
son las siguientes:

En cada instante, la cámara de combustión logra la presión de equilibrio instantánea-

mente.

El frente de llama se propaga perpendicularmente a la superficie en toda la superficie

de combustión y a la misma velocidad de regresión

La carcasa, el grano y la tobera se consideran infinitamente rı́gidas.

1.4.1. Determina el área de combustión, la masa de combustión, el volumen del

alma central y el ratio entre el volumen del grano y el volumen total de
la carcasa en el instante inicial.

Para calcular el área del alma se calculará en alma del tronco de un cono cortado, para ello
se usara el siguiente esquema:

Figura 16: Esquema del grano y nomenclatura usada.

En este caso el Rint es el radio interior de la culata, 5 mm. Conociendo este dato y beta
podemos calcular L′t y Ht′ .
5
L′t = = 143, 18 mm (17)
tg(2)

25
 Report combustible sólido

5
Ht′ = = 143, 27 mm (18)
sin(2)

Una vez concoidos estos datos se puede proceder a calcular Ri2 y Ht , que son los que realmente
interesan.
Ri2 = tg(2) · (LD + L′t ) = 25, 95 mm (19)
LD + L′t
Ht = = 743, 63 mm (20)
cos(2)

Con todos estos datos calculados, se obtienen faáilmente el área de combustión, el volumen
del alma, volumen del cilindro, volumen del grano y volumen de la carcasa.

Ac = π((Rint )2 ‘ + (Ri2 )2 + (Ht − Ht′ )(Rint + Ri2 )) = 0, 0607 m2 (21)

Lp · π
Valma = · ((Rint )2 ‘ + (Ri2 )2 + Rint · Ri2 ) = 5, 203 · 10−4 m3 (22)
3
Vcilindro = π · (Rext )2 · Lp = 47, 12 · 10−4 m3 (23)

Vcarcasa = π · (Rext )2 · Lc = 48, 69 · 10−4 m3 (24)

Vgrano = Vcilindro − Valma = 41, 92 · 10−4 m3 (25)

Sabiendo el valor de todos los volúmenes y que la densidad del grano es 1750 kg/m3 es fácil
obtener la masa y el ratio entre el volumen del grano y la carcasa.

m = ρ · Vgrano = 7, 336 kg (26)

Vgrano
r= = 0, 86 (27)
vcarcasa

Ası́ pues, finalmente los resultados obtenidos se muestran en la siguiente tabla resumen.

Tabla 13: Resultados en el instante inicial

Área de combustión 0,06057 m2

Masa de combustión 7,336 kg
Volumen del alma central 5,203 · 10−4 m3
Ratio 0,86

26
 Report combustible sólido

1.4.2. Determina el área de combustión, la masa de combustión, el volumen del

alma central y el ratio entre el volumen del grano y el volumen total de
la carcasa en el instante en el que el diámetro interior en la parte inferior
del alma central logra el diámetro exterior del grano.

Para este caso se seguira el mismo esquema que en el apartado anterior pero teniendo en
cuente que ahora Rint es desconocido y Ri2 pasa a ser Rext que es de 50 mm.

50
Lt = = 1431, 81 mmL′t = Lt − LD = 831, 81 mm (28)
tg(2)

LT L′T
Ht = = 1432, 685 mmHt′ = = 832, 319 mm (29)
cos(2) cos(2)
Rint = tg(2) · L′t = 29, 047 mm (30)

Ac = π((Rint )2 ‘ + (Ri2 )2 + (Ht − Ht′ )(Rint + Ri2 )) = 0, 1596 m2 (31)

Lp · π
Valma = · ((Rint )2 ‘ + (Ri2 )2 + Rint · Ri2 ) = 3, 0135 · 10−3 m3 (32)
3
Vgrano = Vcilindro − Valma = 17, 02 · 10−4 m3 (33)

m = ρ · Vgrano = 2, 97 kg (34)
Vgrano
r= = 0, 3495 (35)
vcarcasa

El volumen del cilindro y el volumen de la carcasa seguirán siendo los mismos que en el
apartado anterior, ya que sus dimensiones no varı́an.

Los resultados para este apartado son los siguientes:

Tabla 14: Resultados en el instante en el que el diámetro interior en la parte inferior del
alma central logra el diámetro exterior del grano.

Área de combustión 0,1596 m2

Masa de combustión 2,97 kg
Volumen del alma central 3,0135 · 10−3 m3
Ratio 0,3495

27
 Report combustible sólido

1.4.3. Estudia la evolución de la relación de Klemmung (KN ) en función de la

masa de combustible quemado

En el instante inicial, la relación de Klemmung se expresa como el área de combustión inicial,

entre el áera de la gola, lo que de un resultado de:

Ac,0 0, 06057
KN,0 = = −3 = 301, 25 (36)
Ag π(( 16·10
2
)2

Calcular la relación de Klemmung en el instante final es un poco más complejo, para ello se
usará la presión mı́nima de combustión.

1 (Pc,min )1−n
Pc,min = (b · ρp · c∗ · KN,f ) 1−n KN,f = = 186, 18 (37)
b · ρ p · c∗

Podemos calcular la relación de Klemmung máxima sabiendo que esta se producirá cuando el
área de combustión sea máxima, y que esto ocurre en el instante en que en el que el diámetro
interior en la parte inferior del alma central logra el diámetro exterior del grano.

Ac,max 0, 1596
KN,max = = −3 = 793, 78 (38)
Ag π(( 16·10
2
)2

Es decir, la relación de Klemmung tiene forma de pico, parte de un valor de 301,25, aumenta
hasta que el área de combustión es la máxima y toma un valor de 793,78 y después de eso
disminuye hasta que se acaba el combustible y toma un valor de 186,18.

1.4.4. Estudia la evolución de la presión de la cámara de combustión en función

de la relación de Klemmung.

Se conoce el valor mı́nimo de la presión de la cámara de combustión, que es de 1,7 MPa.

Conociendo el valor de la relación de Klemmung máxima, se puede saber también el valor
máximo de la presión en la cámara de equilibrio.

1
Pc,max = (b · ρp · c∗ · KN,max ) 1−n = 19, 058 M P a (39)

Con la misma ecuación se calculará la presión de equilibrio, pero usando la relación de

28
 Report combustible sólido

Klemmung inicial en lugar de la máxima.

1
Peq,0 = (b · ρp · c∗ · KN,0 ) 1−n = 3, 791 M P a (40)

Como era de esperar, la presión en la cámara de combustión aumentara si la relación de

Klemmung aumenta.

1.4.5. Estudia el instante final de la combustión

En primer lugar se calculará el área de combustión final, usando la relación de Klemmung

final.
Ac,f
KN,f = Ac,f = KN,f · Ag = 0, 037 m2 (41)
Ag

Sabiendo el valor del área de combustión, se puede saber también la longitud final del com-
bustible, que serı́a de 0,1245 m. Y con ese dato se calcula también el área del alma que es de
8,96 · 10− 4 m3 .

Para calcular el volumen del grano final y ası́ obtener la masa, al volumen del cilindro se le
restará el volumen del alma y el volumen consumido.

Vgrano,f = Vcilindro − Valma − Vconsumido = 81, 82 · 10−6 m3 (42)

Conociendo el volumen final del grano se procede a calcular la masa de este.

mf = ρp · Vgrano,f = 0,143 kg (43)

mf
= 0, 0195 (44)
m0

Los resultados obtenidos para el instante final son los siguientes:

29
 Report combustible sólido

Tabla 15: Resultados en el instante final

Área de combustión final 0,037 m2

Masa de combustible restante 0,143 kg
Volumen del alma central 8,96 · 10−4 m3
Relación entre la massa final e inicial 0.0195

30
 Report combustible sólido

2. ANEXO: Código de MATLAB empleado

En este anexo se muestra el código Matlab utilizado para la resolución de los apartados 2-7
del perimer problema.

1 clc ;
2 clear all ;
3 close all ;
4

5 % % Datos del problema

6

7 rho = 1717; %Densidad de grano ( Kg / m3 )

8 bstd = 7e -6; %Coef . de temperatura ( est á ndar )
9 tao = 0.15; %Exponente de temperatura (1/ K )
10 n = 0.44; %Í ndice de combusti ó n ()
11 tcomb = 240; %Tiempo de combusti ó n ( s )
12 Ag = 0.025; %Á rea de garganta ( m2 )
13 epsi = 50; %Relaci ó n de expansi ó n ()
14 Pa = 101325; %Presi ó n atmosf é rica ( Pa )
15

16

17 % % Resoluci ó n del problema

18

19 %Inicializaci ó n de variables
20 Pc = 0.1 e6 ;
21 i = 1;
22 syms PePc
23 for ttc = 0:0.01:1
24 %Definici ó n de Ac y Kn ( obtenidos con excel )
25 if ( ttc <0.33)
26 Ac = 72* ttc + 12; %( m2 )
27 Kn = 2880* ttc + 480; %()
28 else
29 Ac = -18*( ttc -1/3) + 36; %( m2 )
30 Kn = -720* ttc + 1680; %()
31 end
32 e = 1;
33 while (e >0.001)

31
 Report combustible sólido

34 %Par á metros de los gases

35 Tc = 58.624* log ( Pc ) + 3075.5;
36 gamma = 0.0028* log ( Pc ) + 1.1528;
37 R = -2.1326* log ( Pc ) + 315.86;
38 %Velocidad caracter ı́ stica
39 c = sqrt ( gamma * R * Tc ) /( gamma *(2/( gamma +1) ) ^(( gamma +1) /(2*( gamma -1)
)));
40 %Presi ó n de c á mara de equilibrio
41 Peq = ( bstd * rho * c * Kn ) ^(1/(1 - n ) ) ;
42 %Diferencia entre Pc obtenida e inicial
43 e = abs ( Peq - Pc ) ;
44 Pc = Peq ;
45 end
46 %Presi ó n de c á mara de equilibrio ( Pa )
47 PC (i ,1) = Pc ;
48 %Velocidad de regresi ó n ( m / s )
49 vr (i ,1) = bstd * Pc ^ n ;
50 %Flujo m á sico ( Kg / s )
51 m (i ,1) = Ac * rho * vr (i ,1) ;
52 %Relaci ó n entre Pe y Pc
53 PEPC (i ,1) = vpasolve ((1/ epsi ) ==(( gamma +1) /2) ^(1/( gamma -1) ) * PePc ^(1/
gamma ) * sqrt ( (( gamma +1) /( gamma -1) ) *(1 - PePc ^(( gamma -1) / gamma ) ) ) , PePc
,0.5) ;
54 %Velocidad de salida de los gases
55 ue (i ,1) = sqrt (2*( gamma /( gamma -1) ) * R * Tc *(1 - PEPC (i ,1) ^(( gamma -1) / gamma
)));
56 %Empuje proporcionado por el motor
57 F (i ,1) = m (i ,1) * ue (i ,1) + epsi * Ag *( PEPC (i ,1) * Pc - Pa ) ;
58 i = i +1;
59 end
60

61 ttc = (0:0.01:1) ’;
62

63 % % Apartado 3
64 figure
65 plot ( ttc , PC , ’ LineWidth ’ ,1) ;
66 grid on
67 grid minor

32
 Report combustible sólido

68 xlabel ( ’t / tc ’) ;
69 ylabel ( ’ Pc [ MPa ] ’) ;
70 title ( ’ Presi ó n de Combusti ó n en Equilibrio en funci ó n del tiempo ’) ;
71 set ( findall ( gcf , ’ - property ’ , ’ FontSize ’) , ’ FontSize ’ ,18) ;
72

73 % % Apartado 4
74

75 figure
76 plot ( ttc , vr , ’ LineWidth ’ ,1) ;
77 grid on
78 grid minor
79 xlabel ( ’t / tc ’) ;
80 ylabel ( ’ Vr [ m / s ] ’) ;
81 title ( ’ Velocidad de Regresi ó n en funci ó n del tiempo ’) ;
82 set ( findall ( gcf , ’ - property ’ , ’ FontSize ’) , ’ FontSize ’ ,18) ;
83

84 figure
85 plot ( ttc ,m , ’ LineWidth ’ ,1) ;
86 grid on
87 grid minor
88 xlabel ( ’t / tc ’) ;
89 ylabel ( ’ Flujo m á sico [ kg / s ] ’) ;
90 title ( ’ Flujo M á sico en funci ó n del tiempo ’) ;
91 set ( findall ( gcf , ’ - property ’ , ’ FontSize ’) , ’ FontSize ’ ,18) ;
92

93 M_total = trapz (( ttc * tcomb ) ,m ) ;

94

95 % % Apartado 5
96

97 figure
98 plot ( ttc , PEPC , ’ LineWidth ’ ,1) ;
99 grid on
100 grid minor
101 xlabel ( ’t / tc ’) ;
102 ylabel ( ’ Pe / Pc ’) ;
103 title ( ’ Relaci ó n entre Presi ó n de C á mara y de Salida en funci ó n del tiempo
’) ;
104 set ( findall ( gcf , ’ - property ’ , ’ FontSize ’) , ’ FontSize ’ ,18) ;

33
 Report combustible sólido

105

106 % % Apartado 6
107

108 figure
109 plot ( ttc , ue , ’ LineWidth ’ ,1) ;
110 grid on
111 grid minor
112 xlabel ( ’t / tc ’) ;
113 ylabel ( ’ u_ { e } [ m / s ] ’) ;
114 title ( ’ Velocidad a la Salida de la tobera en funci ó n del tiempo ’) ;
115 set ( findall ( gcf , ’ - property ’ , ’ FontSize ’) , ’ FontSize ’ ,18) ;
116

117 % % Apartado 7
118

119 figure
120 plot ( ttc ,F , ’ LineWidth ’ ,1) ;
121 grid on
122 grid minor
123 xlabel ( ’t / tc ’) ;
124 ylabel ( ’F ’) ;
125 title ( ’ Empuje en funci ó n del tiempo ’) ;
126 set ( findall ( gcf , ’ - property ’ , ’ FontSize ’) , ’ FontSize ’ ,18) ;

34