HUMIDIFICACIÓN Final
HUMIDIFICACIÓN Final
HUMIDIFICACIÓN Final
EXPERIMENTAL DE LOS
COEFICIENTESINDIVIDUALES
DE TRANSPORTE EN
INTERACCIÓN AIRE-AGUA
Autores: Iliyan Blas Zahonero Sánchez y Álvaro Gento Núñez
Asignatura: Experimentación en ingeniería química III Gr. A
(34765)
Grupo: L5
Pareja: P8
Experimentación en Ingeniería Química III Interacción Aire-Agua
ÍNDICE
1. INTRODUCCIÓN ...................................................................................... 2
2. OBJETIVO ................................................................................................. 5
7. ÁPENDICE .............................................................................................. 27
1
Experimentación en Ingeniería Química III Interacción Aire-Agua
1. INTRODUCCIÓN
La humidificación es una operación que consiste en aumentar la cantidad de vapor
presente en una corriente gaseosa. El vapor puede aumentar pasando el gas a través de un
líquido que se evapora en el gas. Esta transferencia hacia el interior de la corriente gaseosa
tiene lugar por difusión y en la interfase hay, simultáneamente, transferencia de calor y
de materia.
• Una corriente de agua caliente se pone en contacto con una de aire seco (o con
bajo contenido en humedad), normalmente aire atmosférico.
• Parte del agua se evapora, enfriándose así la interfase.
• El seno del líquido cede entonces calor a la interfase, y por lo tanto se enfría.
• A su vez, el agua evaporada en la interfase se transfiere al aire, por lo que se
humidifica.
𝑝𝑤 𝑀𝑤 𝑝𝑤
𝑋= ∙ = 0.622 ∙ (Ec.1)
𝑃 − 𝑝𝑤 𝑀𝑎 𝑃 − 𝑝𝑤
2
Experimentación en Ingeniería Química III Interacción Aire-Agua
𝑝𝑤 (Ec.3)
𝑀𝑀 = 𝑀𝑎 − ∙ (𝑀𝑎 − 𝑀𝑤 )
𝑃
Donde ca es el calor específico del aire y cw es el calor específico del agua (ca =1
kJ/K·Kg) y cw =1,86 kJ/K·kg, respectivamente).
• Entalpía del aire húmedo: es igual a la suma de la entalpía del aire seco
(ia) y la del vapor de agua que contiene (iw)
𝑖𝐺 = 𝑖𝑎 + 𝑖𝑤 ∙ 𝑋
(Ec.5)
(Ec.6)
𝑖𝑎 = 𝑐𝑎 ∙ (𝑡 − 0)
(Ec.7)
𝑖𝑤 = 𝑟0 + 𝑐𝑤 ∙ (𝑡 − 0)
3
Experimentación en Ingeniería Química III Interacción Aire-Agua
𝑖𝐺2 (Ec.10)
𝑑𝑖𝐺
𝑁𝑖 = ∫
𝑖𝐺𝑖 − 𝑖𝐺
𝑖𝐺1
4
Experimentación en Ingeniería Química III Interacción Aire-Agua
2. OBJETIVO
Compresión del funcionamiento de una columna de enfriamiento y análisis del
impacto de los caudales de aire y agua que tienen sobre la transferencia de materia y
energía en la operación, además de la familiarización con las ecuaciones de diseño
pertinentes para el cálculo de los coeficientes de transferencia ya mencionados.
3. DISPOSITIVO EXPERIMENTAL
Se compone por una columna adiabática de paredes mojadas donde el agua y el
aire circulan en contracorriente.
5
Experimentación en Ingeniería Química III Interacción Aire-Agua
4. PROCEDIMIENTO EXPERIMENTAL
4.1 Experimentos con caudal de agua constante y caudal de aire variable
Se fija un caudal de agua y se comprueba que el agua moja las paredes de la columna de humidificación. A continuación, se efectuarán 12
experimentos variando el caudal de aire y en cada uno de ellos se esperará hasta alcanzar el estado estacionario.
4
Experimentación en Ingeniería Química III Interacción Aire-Agua
4 9,5 0,8 100 3,38E-03 28,4 17,4 27,2 80,9 25,5 27,7
7 9,5 0,8 160 5,52E-03 28,3 17,6 27,3 82,2 26,1 27,5
8 9,5 0,8 180 6,19E-03 28,2 16,9 27,4 82,3 26,1 27,5
5
Experimentación en Ingeniería Química III Interacción Aire-Agua
10 9,5 0,8 220 7,67E-03 28,2 17,3 27,4 82,5 26,3 27,6
11 9,5 0,8 240 8,49E-03 28,2 17,4 27,4 82,7 26,3 27,5
12 9,5 0,8 260 9,08E-03 28,3 17,3 27,3 82,7 26,5 27,4
Para el cálculo de caudal de aire que entra al sistema, QG, se precisa de una ecuación de recalibrado anexada al aparato:
también, para deducir la humedad, 1 y 2, a través de la humedad relativa leída en las sondas, 1sonda y 2sonda, se recurre a
unas expresiones de calibrado adjuntadas al dispositivo experimental.
6
Experimentación en Ingeniería Química III Interacción Aire-Agua
Una vez conocidos los caudales, composiciones y temperaturas de entradas y salidas, se puede comprobar el balance de
energía aplicado a la columna de interacción:
Donde M´ es el caudal másico de aire seco, iG1 e iG2 las entalpías del aire a la entrada y a la salida respectivamente, Lm el
caudal másico medio de agua, cL el calor específico medio del agua y T2 y T1 las temperaturas de entrada y salida del agua.
Con las lecturas de los manómetros y mediante la ecuación de calibrado (Ec.11), se calcula el caudal volumétrico de aire a
la entrada (QG). También se debe calcular las humedades a la salida y a la entrada de la columna con sus correspondientes calibrados:
(Ec.12) y (Ec.13).
7
Experimentación en Ingeniería Química III Interacción Aire-Agua
2 24,7 19 3,87 3,57E-03 27,6 17,2 26,8 77,7 13,88 73,01 24,4 27,5
3 22,1 18 3,66 3,57E-03 27,7 17,2 26,8 78,2 13,88 73,64 24,6 27,6
4 19,4 16 3,44 3,57E-03 27,9 17,2 26,9 77,1 13,88 72,26 25,0 27,6
5 16,6 14 3,19 3,57E-03 27,9 16,6 27,1 78,1 13,37 73,51 25,1 27,6
6 14,4 12 2,98 3,57E-03 28,0 16,6 27,2 78,6 13,37 74,14 25,3 27,7
7 11,5 10 2,67 3,57E-03 28,0 17,1 27,2 79,9 13,80 75,78 25,4 27,7
8 8,5 8 2,30 3,57E-03 28,1 16,6 27,4 81,4 13,37 77,69 25,7 27,7
9 6,5 7 2,02 3,57E-03 28,2 16,5 27,5 82,7 13,29 79,35 25,8 27,7
10 3,6 6 1,51 3,57E-03 28,2 16,5 27,6 84,2 13,29 81,29 26,1 27,5
11 2,2 5 1,18 3,57E-03 28,2 16,5 27,4 85,7 13,29 83,25 26,3 27,5
12 0,9 3 0,76 3,57E-03 28,4 17,2 27,8 87,0 13,88 84,97 26,6 27,5
8
Experimentación en Ingeniería Química III Interacción Aire-Agua
A partir de la información experimental, haciendo uso de las ecuaciones necesarias, se determina la variación de energía que experimentan
el aire y el agua a lo alto de la columna (Tabla 4), con el fin de comprobar el cumplimiento o no del balance de energía. Se considerarán válidos
aquellos experimentos que cumplan el balance de energía con un error inferior al 10%.
2 1,07E-03 3,105E-03 1,589E-02 35,525 67,352 28,866 1,218 1,31E-03 1,30E-03 4,15E-02 3,57E-03 4,62E-02 10,21
3 1,02E-03 3,123E-03 1,603E-02 35,672 67,708 28,866 1,216 1,24E-03 1,23E-03 3,95E-02 3,57E-03 4,47E-02 11,67
4 9,55E-04 3,160E-03 1,582E-02 35,968 67,272 28,866 1,212 1,16E-03 1,15E-03 3,61E-02 3,57E-03 3,87E-02 6,77
5 8,86E-04 3,043E-03 1,630E-02 35,669 68,691 28,867 1,209 1,07E-03 1,07E-03 3,52E-02 3,57E-03 3,73E-02 5,40
6 8,27E-04 3,061E-03 1,654E-02 35,815 69,412 28,867 1,205 9,97E-04 9,9E-04 3,34E-02 3,57E-03 3,58E-02 6,68
7 7,41E-04 3,159E-03 1,692E-02 36,065 70,370 28,866 1,202 8,91E-04 8,88E-04 3,05E-02 3,57E-03 3,43E-02 11,14
8 6,39E-04 3,079E-03 1,756E-02 35,962 72,231 28,866 1,199 7,66E-04 7,64E-04 2,77E-02 3,57E-03 2,98E-02 7,06
9 5,60E-04 3,077E-03 1,806E-02 36,058 73,601 28,866 1,197 6,71E-04 6,69E-04 2,51E-02 3,57E-03 2,83E-02 11,36
10 4,19E-04 3,077E-03 1,863E-02 36,058 75,152 28,866 1,195 5,01E-04 4,99E-04 1,95E-02 3,57E-03 2,09E-02 6,49
11 3,29E-04 3,077E-03 1,886E-02 36,058 75,542 28,866 1,194 3,92E-04 3,91E-04 1,54E-02 3,57E-03 1,79E-02 13,65
12 2,12E-04 3,254E-03 1,973E-02 36,710 78,178 28,865 1,190 2,52E-04 2,51E-04 1,04E-02 3,57E-03 1,34E-02 22,36
9
Experimentación en Ingeniería Química III Interacción Aire-Agua
Con los experimentos que cumplan adecuadamente el balance de energía, se procede a la determinación de la pendiente de la
línea de unión (-hL/k’) y el número de unidades de transferencia (Ni) mediante la construcción analítica de la curva de Mickley (Tabla
5). El procedimiento para la construcción de la curva de Mickley se describe detalladamente en el Apéndice II.
Se comprobará asimismo que la aproximación de la media logarítmica resulta adecuada para el cálculo del número de unidades
de transferencia (Niml).
2 1,23E-03 3,57E-03 26,405 68,750 82,647 33,078 14,939 22,820 1,40 1,30 7,42
3 1,15E-03 3,57E-03 26,354 70,143 83,572 34,175 16,301 24,145 1,30 1,25 3,32
4 1,07E-03 3,57E-03 26,408 70,457 83,203 34,788 14,512 23,191 1,42 1,24 12,71
5 9,94E-04 3,57E-03 26,472 71,177 84,017 35,362 14,606 23,474 1,43 1,33 7,26
6 8,88E-04 3,57E-03 26,616 71,582 83,798 35,517 13,428 22,710 1,51 1,38 8,44
7 7,64E-04 3,57E-03 26,555 72,607 83,836 36,645 11,605 21,777 1,67 1,35 18,86
8 6,69E-04 3,57E-03 26,541 72,968 83,390 36,911 9,789 20,434 1,84 1,37 25,61
9 4,99E-04 3,57E-03 26,553 72,968 83,390 36,911 8,238 19,118 2,04 1,38 32,31
10
Experimentación en Ingeniería Química III Interacción Aire-Agua
Conocidos los valores del número de unidades de transferencia, la determinación del coeficiente de transferencia de materia (k’) y
de los coeficientes de transferencia de energía (hL y hC) es inmediata (Tabla 6).
2 1,23E-03 3,57E-03 26,40 1,30 0,571 15,07 1,0058 1,0298 1,0178 0,581
3 1,15E-03 3,57E-03 26,35 1,25 0,554 14,60 1,0059 1,0294 1,0177 0,564
4 1,07E-03 3,57E-03 26,41 1,24 0,517 13,65 1,0057 1,0303 1,0180 0,526
5 9,94E-04 3,57E-03 26,47 1,33 0,450 11,92 1,0057 1,0308 1,0182 0,459
6 8,88E-04 3,57E-03 26,62 1,38 0,386 10,28 1,0059 1,0315 1,0187 0,393
7 7,64E-04 3,57E-03 26,55 1,35 0,340 9,03 1,0057 1,0327 1,0192 0,347
8 6,69E-04 3,57E-03 26,54 1,37 0,294 7,81 1,0057 1,0336 1,0197 0,300
9 4,99E-04 3,57E-03 26,55 1,38 0,217 5,76 1,0057 1,0346 1,0202 0,221
Para comprobar la dependencia del coeficiente de transferencia de materia y de los coeficientes de transferencia de energía con el
caudal de aire que circula por la columna, se representará gráficamente log (coef) = f (log M’), siendo coef k’, h L o hC. A partir de estas
11
Experimentación en Ingeniería Química III Interacción Aire-Agua
representaciones se propondrán, para el caudal de líquido constante, tres ecuaciones de la forma k’ = (M’), hC = (M’), hL = (M’) y se realizará
un análisis crítico de los resultados obtenidos.
k' = f(M')
k' (kg/s·m2·X)
y = 1058,8x1,1198
R² = 0,994
0,200
4,50E-04
M' (kg a.s/s)
Gráfica 1 Evolución del coeficiente de transferencia de materia en función del caudal de aire.
12
Experimentación en Ingeniería Química III Interacción Aire-Agua
hL = f(M')
hL (kg/s·m2·ºC)
y = 23911x1,0976
R² = 0,9943
5,00
4,50E-04
M' (kg a.s/s)
Gráfica 2 Evolución del coeficiente de transferencia de calor desde el agua en función del caudal de aire.
13
Experimentación en Ingeniería Química III Interacción Aire-Agua
hc = f(M')
hc (kg/s·m2·ºC)
y = 1057,2x1,117
R² = 0,994
0,20
4,50E-04
M' (kg a.s/s)
Gráfica 3 Evolución del coeficiente de transferencia de calor desde el aire en función del caudal de aire.
14
Experimentación en Ingeniería Química III Interacción Aire-Agua
2 9,5 0,70 2,43 1,97E-03 28,1 16,7 26,6 79,8 13,46 75,65 24,1 27,4
3 9,5 0,70 2,43 2,61E-03 28,4 16,9 27,0 80,4 13,63 76,41 24,9 27,6
4 9,5 0,80 2,43 3,38E-03 28,4 17,4 27,2 80,9 14,06 77,05 25,5 27,7
5 9,5 0,80 2,43 4,23E-03 28,5 17,0 27,3 81,2 13,71 77,43 25,9 27,7
6 9,5 0,80 2,43 4,86E-03 28,4 17,0 27,3 81,6 13,71 77,94 25,9 27,5
7 9,5 0,80 2,43 5,52E-03 28,3 17,6 27,3 82,2 14,23 78,71 26,1 27,5
8 9,5 0,80 2,43 6,19E-03 28,2 16,9 27,4 82,3 13,63 78,84 26,1 27,5
15
Experimentación en Ingeniería Química III Interacción Aire-Agua
10 9,5 0,80 2,43 7,67E-03 28,2 17,3 27,4 82,5 13,97 79,10 26,3 27,6
11 9,5 0,80 2,43 8,49E-03 28,2 17,4 27,4 82,7 14,06 79,35 26,3 27,5
12 9,5 0,80 2,43 9,08E-03 28,3 17,3 27,3 82,7 13,97 79,35 26,5 27,4
2 1,07E-03 3,105E-03 1,589E-02 35,525 67,352 28,866 1,218 8,02E-04 8,04E-04 2,572E-02 1,975E-03 2,724E-02 -5,90
3 1,02E-03 3,123E-03 1,603E-02 35,672 67,708 28,866 1,216 8,01E-04 8,03E-04 2,676E-02 2,612E-03 2,948E-02 -10,16
4 9,55E-04 3,160E-03 1,582E-02 35,968 67,272 28,866 1,212 8,01E-04 8,03E-04 2,744E-02 3,379E-03 3,107E-02 -13,22
5 8,86E-04 3,043E-03 1,630E-02 35,669 68,691 28,867 1,209 8,01E-04 8,03E-04 2,796E-02 4,229E-03 3,182E-02 -13,83
6 8,27E-04 3,061E-03 1,654E-02 35,815 69,412 28,867 1,205 8,01E-04 8,03E-04 2,833E-02 4,862E-03 3,252E-02 -14,80
7 7,41E-04 3,159E-03 1,692E-02 36,065 70,370 28,866 1,202 8,01E-04 8,04E-04 2,857E-02 5,522E-03 3,232E-02 -13,10
8 6,39E-04 3,079E-03 1,756E-02 35,962 72,231 28,866 1,199 8,01E-04 8,04E-04 2,935E-02 6,190E-03 3,623E-02 -23,42
9 5,60E-04 3,077E-03 1,806E-02 36,058 73,601 28,866 1,197 8,01E-04 8,04E-04 2,968E-02 7,042E-03 4,415E-02 -48,77
10 4,19E-04 3,077E-03 1,863E-02 36,058 75,152 28,866 1,195 8,01E-04 8,04E-04 2,931E-02 7,671E-03 4,169E-02 -42,22
16
Experimentación en Ingeniería Química III Interacción Aire-Agua
12 2,12E-04 3,254E-03 1,973E-02 36,710 78,178 28,865 1,190 8,01E-04 8,04E-04 2,900E-02 9,078E-03 3,415E-02 -17,76
2 8,03E-04 2,61E-03 26,77 69,94 82,10 33,41 12,25 21,09 1,58 1,21 23,07
3 8,03E-04 3,38E-03 26,71 72,04 105,43 35,26 34,49 34,87 0,98 0,77 21,49
4 8,03E-04 4,23E-03 26,65 73,43 82,68 36,70 11,15 21,44 1,62 1,14 29,86
5 8,03E-04 4,86E-03 26,59 73,40 83,43 36,82 11,60 21,83 1,61 1,49 7,65
6 8,04E-04 5,52E-03 27,20 74,24 83,63 37,51 11,34 21,87 1,63 1,47 9,38
7 8,04E-04 9,08E-03 26,91 75,60 116,46 39,02 43,80 41,36 0,87 0,79 9,68
17
Experimentación en Ingeniería Química III Interacción Aire-Agua
2 8,03E-04 2,61E-03 26,77 1,21 0,398 10,65 1,0058 1,0298 1,0178 0,405
3 8,03E-04 3,38E-03 26,71 0,77 0,629 16,79 1,0059 1,0294 1,0177 0,640
4 8,03E-04 4,23E-03 26,65 1,14 0,424 11,31 1,0057 1,0303 1,0180 0,432
5 8,03E-04 4,86E-03 26,59 1,49 0,324 8,62 1,0057 1,0308 1,0182 0,330
6 8,04E-04 5,52E-03 27,20 1,47 0,328 8,93 1,0059 1,0315 1,0187 0,334
7 8,04E-04 9,08E-03 26,91 0,79 0,614 16,51 1,0057 1,0327 1,0192 0,626
18
Experimentación en Ingeniería Química III Interacción Aire-Agua
5.2.3 Determinación de la influencia de los caudales de agua y aire en los coeficientes de transporte
Se propondrán tres ecuaciones que expresen la dependencia de los coeficientes de transporte con los caudales de agua y aire, en la forma
k’ = (Lh, M’), hC = ((Lh, M’), hL = ((Lh, M’).
19
Experimentación en Ingeniería Química III Interacción Aire-Agua
k' = f(Lm)
0,300
0,001 0,010
Lm (kg/s)
Gráfica 4 Evolución del coeficiente de transferencia de materia en función del caudal de agua.
20
Experimentación en Ingeniería Química III Interacción Aire-Agua
hL = f(Lm)
hL (kg/s·m 2·ºC)
y = 65,245x0,296
R² = 0,9748
8,00
0,001 0,010
Lm (kg/s)
Gráfica 5 Evolución del coeficiente de transferencia de calor desde el agua en función del caudal de agua.
21
Experimentación en Ingeniería Química III Interacción Aire-Agua
hC = f(Lm)
hC (kg/s·m2·ºC)
y = 2,2566x0,2777
R² = 0,9543
0,300
0,001 0,010
Lm (kg/s)
Gráfica 6 Evolución del coeficiente de transferencia de calor desde el aire en función del caudal de agua.
22
Experimentación en Ingeniería Química III Interacción Aire-Agua
Los datos aportados son tales como los que deben aparecer en la Tabla 1 y 2 (hm,
Ppo, N, L1, … ) tanto en la experiencia de caudal de agua constante como en el de caudal
de aire constante.
El error permitido entre los balances de energía del agua y el aire debía rondar el
10% de error (relativo) y se precisaba de 7 valores para continuar con el cálculo de los
coeficientes de transporte. Sin embargo, los resultados que aportaban los datos tomados
en el laboratorio no cumplían las premisas ya mencionadas.
22
Experimentación en Ingeniería Química III Interacción Aire-Agua
Tabla 11 Comparación entre los resultados obtenidos con los datos tomados en el laboratorio y
los proporcionados (Caudal de agua constante).
DATOS LABORATORIO DATOS ANEXADOS
M'·(iG2-iG1) Lm·cL·(T2-T1) eT (%) M'·(iG2-iG1) Lm·cL·(T2-T1) eT (%)
(kJ/s) (kJ/s) (kJ/s) (kJ/s)
3,94E-02 4,36E-02 -10,52 4,44E-02 4,77E-02 7,00
3,65E-02 4,18E-02 -14,76 4,15E-02 4,62E-02 10,21
3,35E-02 4,01E-02 -19,56 3,95E-02 4,47E-02 11,67
3,08E-02 3,49E-02 -13,20 3,61E-02 3,87E-02 6,77
2,54E-02 3,73E-02 -47,12 3,52E-02 3,73E-02 5,40
2,52E-02 4,01E-02 -59,07 3,34E-02 3,58E-02 6,68
2,34E-02 4,51E-02 -92,46 3,05E-02 3,43E-02 11,14
2,22E-02 2,61E-02 -17,75 2,77E-02 2,98E-02 7,06
1,99E-02 2,44E-02 -22,46 2,51E-02 2,83E-02 11,36
1,76E-02 2,09E-02 -18,63 1,95E-02 2,09E-02 6,49
1,43E-02 1,57E-02 -10,02 1,54E-02 1,79E-02 13,65
1,03E-02 1,22E-02 -18,17 1,04E-02 1,34E-02 22,36
Tabla 12 Comparación entre los resultados obtenidos con los datos tomados en el laboratorio y
los proporcionados (Caudal de aire constante).
DATOS LABORATORIOS DATOS ANEXADOS
M'·(iG2-iG1) Lm·cL·(T2-T1) eT (%) M'·(iG2-iG1) Lm·cL·(T2-T1) eT (%)
(kJ/s) (kJ/s) (kJ/s) (kJ/s)
2,76E-02 1,92E-02 30,59 2,44E-02 2,57E-02 -5,65
23
Experimentación en Ingeniería Química III Interacción Aire-Agua
Una posible causa de que los datos no cumplan con el error mínimo se encuentra
en la variación de entalpía del aire, y a su vez a la mínima variación temperatura que
experimenta el aire a lo largo de la columna; ya que a medida que se disminuye el caudal
de aire, la variación de la temperatura debería reducirse gradualmente pues existe una
menor trasferencia de calor desde la fase líquida (teniendo en cuenta también el agua
evaporada que se transfiere al aire) y esto no ocurre de este modo.
Una vez seleccionados los datos que sí cumplan la igualdad de la entalpía del aire
y la temperatura del agua (balances de energía), se procede al cálculo de los coeficientes
de transporte:
24
Experimentación en Ingeniería Química III Interacción Aire-Agua
Se realizan varias experiencias variando tanto el caudal del agua como del aire
para estudiar cual de las dos opciones proporcionan un mejor enfriamiento. Para analizar
como el caudal de ambos, aire y agua, afecta a la operación se debe prestar atención a la
evolución de los coeficientes de transporte.
25
Experimentación en Ingeniería Química III Interacción Aire-Agua
26
Experimentación en Ingeniería Química III Interacción Aire-Agua
7. ÁPENDICE
• Los datos de entalpía y temperatura del sistema aire-agua se ajustan al
polinomio:
𝑖𝐺𝑖 = 7.806793 + 1.976076 ∙ 𝑡𝑖 − 2.093384 ∙ 10−4 ∙ 𝑡𝑖2 + 1.172473 ∙ 10−3 ∙ 𝑡𝑖3 (Ec.21)
Donde la entalpía (iGi) viene dada en kJ/kg a.s y la temperatura (ti) en ⁰C.
ℎ (Ec.22)
𝑡𝑗+1 = 𝑡𝑗 + ∙ (𝐾1 + 2𝐾2 + 2𝐾3 + 𝐾4 )
6
𝑗
𝐾1 = 𝑓(𝑖𝐺 , 𝑡𝑗 )
𝑗 1 1
𝐾2 = 𝑓(𝑖𝐺 + ∙ ℎ, 𝑡𝑗 + ∙ ℎ ∙ 𝐾1 )
2 2
𝑗 1 1
𝐾3 = 𝑓(𝑖𝐺 + ∙ ℎ, 𝑡𝑗 + ∙ ℎ ∙ 𝐾2 )
2 2
𝑗 1
𝐾4 = 𝑓(𝑖𝐺 + ∙ ℎ, 𝑡𝑗 + ℎ ∙ 𝐾3 )
2
𝑖𝐺2 −𝑖𝐺1
Si se selecciona h como se puede resolver la ecuación (8) paso a paso.
10
27