Cálculo Termodinámico de Un Explosivo
Cálculo Termodinámico de Un Explosivo
Cálculo Termodinámico de Un Explosivo
Curso 09/10
%peso
AGUA 13,230
FUEL OIL 5,110
NITRATO
80,190
AMONICO
Monoleato de
1,470
Sorbitano
C H O N PM %peso %PM
0 2 1 0
AGUA 18 13,23 7,34
0,000 14,676 7,338 0,000
11 24 0 0
FUEL OIL 156,310 5,11 0,33
3,590 7,833 0,000 0,000
0 4 3 2
NIT. AMON 80 80,19 10,01
0,000 40,030 30,022 20,015
24 44 6 0
MonDSorbtn 428,6 1,47 0,03
0,822 1,507 0,205 0,000
∑ 4,412 64,046 37,566 20,015 998,373 100
Sumamos los moles de C, H, O y N de cada sustancia y hallamos así la fórmula global que es:
C4,412H64,064O37,566N20,015
1
Energía y Entalpía de formación
Miramos en tablas la energía (Ef ) y entalpía (Hf) de formación de cada sustancia que
necesitaremos más adelante. Y hallamos el %ΔEf o energía interna, EO:
NITRATO
-1057 80,19 -847,608 -1091
AMONICO
Monoleato de
43,64 1,47 0,642 9,1
Sorbitano
∑ -5250,880 100 -1377,724 -5414,180
Balance de oxígeno
BO = −5,25803556 %
En nuestro caso, como la fórmula sólo tiene C, H, O y N los únicos elementos que consumen
oxígeno son el C y el H, el primero consume 2 moles de O2 por mol y el segundo 0,5 moles de
O2 por mol. Así hallamos el oxígeno necesario para oxidar los elementos.
Tenemos un balance de oxígeno negativo, quiere decir que por cada 100 gramos de explosivo
tenemos un defecto de oxígeno de 5,258 gramos.
2
Temperatura de explosión
Normalmente la temperatura de explosión está entre 2000 y 5000K. Como no sabemos cuál es
la del explosivo muestra, vamos a ir probando.
Empezamos suponiendo una temperatura de 2000K. Calculamos los moles de productos que
se forman a dicha temperatura. Para ello resolvemos el siguiente sistema de ecuaciones:
Supuesta una temperatura de 2000K, miramos en las tablas los valores de las constantes de
equilibrio K1 y K2
3
N n6
20,015 10,0074414
n2 0,71912306
n1 3,17955897
n2 0,71912306
n3 30,5487953
n4 1,526228
n5 0,521
n6 10,0074414
Calor de explosión y calor sensible de los productos (seguimos con una temperatura de
2000K)
Para hallar el calor de explosión y el sensible de los productos debemos averiguar primero la
entalpía en función de la temperatura, , y la energía de formación a 298K, de los
productos mirando en las tablas.
2000
298
ΔEf n H-H298
CO2 -94,05 3,17955897 21,857
CO -26,76 0,71912306 13,561
H20 -57,5 30,4875761 17,37
H2 0 1,53526112 12,651
C 0 0,51333216 8,442
N2 0 10,0074414 13,418
ng=45,9289606 ∑(H-H298)*ni
∑ n*ΔEf 298 -2071,31688 766,852832
Q1 693,593175 kcal/kg
4
El calor sensible de los productos es inferior al calor de explosión, esto quiere decir que
debemos aumentar la temperatura, ¿pero cuánto? Una primera estimación es aumentar la
temperatura supuesta inicialmente unos 3 grados por Kcal/kg de diferencia de calores.
N n6
20,015 10,0074414
n1 3,202182
n2 1,21964508
n3 29,941808
n4 2,0810292
n5 -0,010
Ha salido n5 negativa, quiere decir que no se forma grafito, C, entonces debemos resolver el
siguiente sistema:
n1 3,1944892877
n2 1,055343469
n3 29,95931361
n4 2,063523598
5
2300
ΔEf 298 n H-H298
CO2 -94,05 3,19448929 26,212
CO -26,76 1,05534347 16,175
H20 -57,5 29,9593136 21,1
H2 0 2,0635236 15,146
N2 0 10,0074414 16,015
46,2702984 ∑ (H-H298)*ni
∑ n*ΔEf 298 -2051,34324 924,468953
Q2 673,619538 kcal/kg
Ahora el calor sensible sale mayor que el de explosión, por lo que para hallar la nueva
temperatura de estimación interpolaremos linealmente a partir de las dos temperaturas halladas
anteriormente.
Tenemos así:
− 611 ,526874
T = 2300 + ( 2000 − 2300 ) 673 ,619538 = 2144 ,535956 K
740 ,406909 − 611 ,526874
Pero para esta temperatura no tenemos los valores de las constantes de equilibrio en las tablas
por lo que tendremos que interpolar los valores a partir de los que tenemos para las
temperaturas 2100K (T1) y 2200K (T2), intervalo en el cual se encuentra la temperatura de
nuestros cálculos. Tendremos que interpolar para hallar el valor de K1 y de K2. En ambos casos
aplicaremos la siguiente fórmula:
6
N n6
20,015 10,0074414
n1 3,34223117
n2 0,97636723
n3 29,9049875
n4 2,11784969
n5 0,093
Para hallar las entalpías interpolaremos al igual que hicimos con las constantes de equilibrio.
2144,535956
298
ΔEf N H-H298
CO2 -94,05 3,34223117 23,9496621
CO -26,76 0,97636723 14,8179082
H20 -57,5 29,9049875 19,1566995
H2 0 2,11784969 13,8455392
C 0 0,09341579 9,2922075
N2 0 10,0075 14,6656814
ng=46,4423514 ∑ (H-H298)*ni
∑ n*ΔEf 298 -2060,00121 844,351501
Q2 682,277507 kcal/kg
El calor sensible sale menor que el de explosión por lo que la temperatura real de explosión
será mayor que la estimada. Para estimar una nueva temperatura volvemos a interpolar, esta
vez entre 2157,326877K y 2300K.
− 740 ,406909
T = 2300 + ( 2144 ,535956 − 2300 ) 682 ,277507 = 2164 ,01379 K
673 ,951405 − 740 ,406909
Sale una temperatura 20K mayor que la anterior, es un error demasiado grande por lo que
seguiremos iterando, hasta que el error sea menor de 10K entre la temperatura nueva y la
antigua.
7
K1 5,2046 4,876 5,08634933
K2 67,67 45,61 59,731443
N n6
20,015 10,0074414
n1 3,3312894
n2 1,00727651
n3 29,8959618
n4 2,12687542
n5 0,073
2164,01379
Q 681,556586 kcal/kg
Qs 682,549519 kcal/kg
− 673 ,951205
T = 2144 ,535956 + ( 2164 ,01379 − 2144 ,535956 ) 681 ,556586 = 2161 ,76444 K
682 ,549519 − 673 ,951205
8
Es un error aceptable, podríamos parar aquí las iteraciones, sin embargo continuamos con una
iteración más para comprobar que efectivamente con más iteraciones el error disminuye,
aproximándose más a cero.
N n6
20,015 10,0074414
n1 3,33352095
n2 1,00394173
n3 29,8948335
n4 2,12800374
n5 0,075
2161,76444
Q 681,612346 kcal/kg
Qs 681,688647 kcal/kg
Ahora vemos claramente que el calor sensible y el de explosión se aproximan bastante.
− 682 ,549519
T = 2164 ,01379 + ( 2161 ,76444 − 2164 ,01379 ) 681 ,612346 = 2161 ,56509 K
681 ,688647 − 682 ,549519
9
ε = Tanterior − Tnueva = 2161 ,76444 − 2161 ,01379 = 0,19935386
Paramos las iteraciones. En los cálculos siguientes redondearemos resultados, así que para
facilitar el cálculo tomaré como ng = 46,44 moles.
Hallamos la masa molecular media de los productos gaseosos, en g/mol. Se puede hallar
de dos maneras, bien dividiendo los gramos de productos gaseosos entre los moles
totales o bien restando al kg inicial de explosivo los gramos de productos no gaseosos y
dividiendo por los moles totales. En ambos casos salen valores parecidos.
∑n
gases
i ⋅ PM i
1000 − PM C ⋅ nC
M = =
ng ng
1000 − 12 ⋅ 0,075
M = ≈ 21,5 g
46 ,44 mol
1 1
φ = ng ⋅ M 2
⋅Q 2
1 1
φ = 46,44 ⋅ 21,5 2 ⋅ 2850 2 = 11503,138
DC J = A ⋅ φ 2 ⋅ ( 1 + B ⋅ ρ O )
1
• Velocidad de detonación:
A = 22,328
B = 1,3
1,47 ⋅ ρ O
• Densidad de detonación: ρ CJ =
1 + 0,05625 ⋅ ρ O
1,47 ⋅1,09
ρCJ = = 1,509734409 ≈ 1,51 g / cm 3
1 + 0,05625 ⋅1,09
10
ΓCJ + 1
• Coeficiente adiabático: ρ CJ = ⋅ ρO
ΓCJ
1 + 0,05625 ⋅ ρ O 1 + 0,05625 ⋅ 1,09
ΓCJ = = = 2,59688
0,47 − 0,05625 ⋅ ρ O 0,47 − 0,05625 ⋅ 1,09
2,59688 +1
ρCJ = ⋅1,09 = 1,5097344 ≈ 1,51
2,59688
Q 2850
PRP = = = 0,732647 ≈ 73,26 %
Q ANFO 3890
Q ρ 2850 1,09
PRV = = = 0,9982 ≈ 99 ,82 %
Q ANFO ρ ANFO 3890 0,8
RESULTADOS OBTENIDOS
Fórmula C4,412H64,064O37,566N20,015
BO (%) -5,258
T (K) 2160
Q (kJ/kg) 2850
F (kJ/kg) 834
11