Facultad de Ingeniería

Escuela de Ingeniería de Procesos Industriales

MÉTODOS DE RESOLUCIÓN DE
ECUACIONES NO LINEALES

París Paredes Inostroza

Claudia Aravena
Abril, 2021
Solución de ecuaciones no lineales
• En una ecuación no lineal la variable que se necesita determinar (x) se encuentra en términos
como: 𝑒 𝛽𝑥 , 𝑙𝑛 𝛽 + 𝑥
• No es posible despejar la variable que se desea determinar.
• La solución de estas ecuaciones es un proceso iterativo: se necesita de una valor inicial de la
variable que se calcula.

Ejemplos:
• Cálculo de la velocidad (v) con la que circula un fluido a través de una tubería (BEM):
∆𝑝 ∆𝑣 2 𝐿 𝑣2
− + + 4𝑓 = 0, no puede despejarse v ya que el factor de fricción (f) es una
𝜌 2 𝐷 2
función no lineal de v.

• Cálculo de las temperaturas de salida de un intercambiador de calor de área conocida.

Métodos de solución de ecuaciones no lineales
• Método de sustitución sucesiva
• Método de interpolación lineal
• Método de Newton-Raphson

• Los 3 métodos se explicarán usando el mismo ejemplo.

• Ejemplo: El factor de fricción f para el régimen turbulento de un fluido incompresible que

circula a través de una tubería puede calcularse desde la ecuación no lineal de Colebrook:
𝜀
1 2,51
= −0,86 × 𝐿𝑛 𝐷
+ (1)
𝑓 3,7 𝑁𝑅𝑒× 𝑓

• En donde e/D es la rugosidad relativa del material de la tubería, D el diámetro de la tubería, y

Re en número de Reynolds:
𝑣𝐷𝜌
N𝑅𝑒 = (2)
𝜇
Método de sustitución sucesiva
• Este método requiere que la ecuación que contiene la incógnita (x) que se desea calcular que se
escriba de la siguiente forma:
𝑥=𝑔 𝑥 (3)
• En donde g(x) una función de x.
• La solución corresponde al punto donde se intersecta la recta y=x y la curva x=g(x).
• La fórmula iterativa es:
𝑥𝑛+1 = 𝑔 𝑥𝑛 (4)
• Se necesita una estimación inicial de x, x0, para iniciar el cálculo.
• Una condición suficiente para la convergencia es:
𝑑𝑔
<1 (5)
𝑑𝑥
(f0,g(f0))
• La solución al problema se obtiene cuando:
𝑥𝑛+1 − 𝑥𝑛 ≤ 𝑇𝑜𝑙𝑒𝑟𝑎𝑛𝑐𝑖𝑎 (6) (f1,f1)

• En donde Tolerancia es un valor pequeño (ej. 0.1; 0.0005)

Dependerá de las magnitudes con las que se trabaje.
(f1,g(f1))

• La Figura 1 muestra el procedimiento del método de (f2,f2)

sustitución sucesiva para el ejemplo.
• Reordenando la ecuación (1) para obtener la fórmula iterativa
en (3) se tiene:
1
𝑓𝑛+1 = 𝑔 𝑓𝑛 = 2 (7)
𝜀/𝐷 2,51
0,86×𝐿𝑛 + f0, valor para iniciar la búsqueda
3,7 𝑁𝑅𝑒× 𝑓𝑛
Figura 1. Método sustitución sucesiva
 • En el archivo napellido_T2.xlsx, Hoja1, implemente la solución del ejemplo usando el método de
sustitución sucesiva como se muestra en la Figura 2.
• Para verificar que el punto f0 cumple con la condición de convergencia derive la ecuación (5) o
una aproximación de diferencias finitas:
𝑑𝑔 𝑔 𝑓0 +∆𝑓 −𝑔 𝑓0
ቚ ≅ (8)
𝑑𝑓 𝑓 ∆𝑓
0
• Verifique que obtiene una solución con distintos valores de e/D, NRe.
Como son pasos iterativos, es
necesario fijar algunas celdas
(datos principalmente) para
Datos del problema esto usar tecla F4. En el
ejemplo la celda B4 está fijada
Valor de f para iniciar y queda con los signos $$
el procedimiento
iterativo
dg/df <0, por lo que f0
seleccionado puede ser usado
Con los datos del
para iniciar el proceso iterativo
problema calcule
dg/df usando Ec. (7)

Prueba lógica:
Si abs(fn –fn-1)<Tolerancia entonces
En la 2° iteración y “mostrar valor de fn” (solución)
posteriores calcula Sino “” (mostrar celda en blanco)
usando el valor en
la iteración
anterior Cálculo g(f) (que es igual a fn)
con datos del problema y valor
de fn-1
Figura 2. Solución del ejemplo 1 usando el método sustitución sucesiva en planilla de cálculo.
Método de interpolación lineal

La ecuación no lineal debe expresarse como:

ℎ 𝑥 =0 (9)
• Paso 1. Se necesita conocer 2 valores de x, x1, x2, (x1<x2) entre los cuales se encuentre la raíz
o cero de la Ec. (8). Esta condición se cumple cuando:
ℎ 𝑥1 × ℎ 𝑥2 < 0 (10)
• Paso 2. Si 𝑥2 − 𝑥1 ≤ 𝑇𝑜𝑙𝑒𝑟𝑎𝑛𝑐𝑖𝑎, entonces se ha encontrado la raíz de la Ec. (9).
• Paso 3. Se debe reducir el tamaño del intervalo que contiene la raíz. Para ello se busca el
valor de x3 (en el intervalo [x1 x2]) usando una interpolación lineal.
• La recta (y=mx+n) que une los puntos (x1 , h(x1)) y (x2 , h(x2)), tiene pendiente (m) e
intersección (n) iguales a:
ℎ 𝑥1 −ℎ 𝑥2
𝑚= ; 𝑛 = ℎ 𝑥1 − 𝑚 × 𝑥1 (11)
𝑥1 −𝑥2
• El valor de x3 se obtiene cuando y=0:
𝑛
𝑦 = 0 = 𝑚 × 𝑥3 + 𝑛 → 𝑥3 = − (12)
𝑚
• Paso 4. Se debe determinar en que intervalo se encuentra la raíz de la Ec. (9), [x1, x3] o [x3, x2].
• Si ℎ 𝑥1 × ℎ 𝑥3 < 0 entonces 𝑥2 = 𝑥3 y ℎ 𝑥2 = ℎ 𝑥3 , y se repite el
procedimiento desde Paso 2.
• Si ℎ 𝑥3 × ℎ 𝑥2 < 0 entonces 𝑥1 = 𝑥3 y ℎ 𝑥1 = ℎ 𝑥3 , y se repite el
procedimiento desde Paso 2.
 (x1, h(x1))

x3, 2° iteración

x3, 1° iteración

(x2, h(x2)), 1° iteración

(x2, h(x2)), 3° iteración(x2, h(x2)), 2° iteración

Figura 3. Solución del ejemplo 1 usando el método interpolación lineal.
 • En el archivo napellido_T2.xlsx, Hoja2, implemente la solución del ejemplo usando el método de
interpolación lineal como se muestra en la Figura 4.
• Verifique que obtiene una solución con distintos valores de e/D, NRe.

Datos del problema

Calcula
pendiente e
intercepto de la
Busca recta
intervalo
en el que
se
encuentra
la solución

Figura 4. Solución del ejemplo 1 usando el método interpolación lineal en

planilla de cálculo. Prueba lógica:
Si abs(f1 –f2)<Tolerancia entonces
“mostrar valor de fn” (solución)
Sino “” (mostrar celda en blanco)
En la segunda iteración para f1 y f2 debe usar =si para determinar
a cual de los valores remplazará el f3, paso 4 del método
Método de Newton-Raphson
• La ecuación no lineal debe expresarse como en Ec. (9).
• La función h(x) se expande en series de Taylor en torno a un valor estimado de la raíz, x1:
ℎ′′ 𝑥
ℎ 𝑥 = ℎ 𝑥1 + ℎ′ 𝑥1 𝑥 − 𝑥1 + 1
𝑥 − 𝑥1 2 + ⋯ (13)
2
• Se retienen los 2 primeros términos (lado derecho), lo que equivale a linealizar la función h(x).
• Como se desea calcular un cero o raíz de la función h, se tiene,
ℎ 𝑥 = ℎ 𝑥1 + ℎ′ 𝑥1 𝑥 − 𝑥1 = 0 (14)
• Despejando x, se obtiene una fórmula general de iteración:
ℎ 𝑥𝑛
𝑥𝑛+1 = 𝑥𝑛 − (15)
ℎ′ 𝑥𝑛
• Este método al igual que los anteriores converge a la mínimos locales (ver Figura 5).
(a) (b)
h(x)
h(x)

Solución h(x)=0

x1 x1 x1 x2 x3

Figura 5. (a) Convergencia, (b) funcionamiento del método Newton-Raphson.

• En el archivo napellido_T2.xlsx, Hoja3, implemente la solución del ejemplo usando el método
de Newton-Raphson como se muestra en la Figura 6.
• Use una aproximación de diferencias finitas (Ec. 8) para calcular la derivada (h’(x)) en Ec. (15).
• Verifique que obtiene una solución con distintos valores de e/D, NRe.

Datos del problema

Para el cálculo de h’(f)
Valor inicial

Prueba lógica:
Si abs(fn –fn+1)<Tolerancia
entonces “mostrar valor de
fn” (solución)
Sino “” (mostrar celda en
blanco)
Figura 6. Solución del ejemplo 1 usando el método Newton-Raphson en planilla de cálculo.
• En el ejemplo estudiado el valor del factor f depende del valor del número de Reynolds (Re), y la rugosidad
relativa de la tubería (e/D).
• Si se necesita calcular f con distintos valores de Re y e/D, se pueden cambiar los valores en las celdas
respectivas en la sección de datos.
• Una forma de generar resultados para distintas combinaciones de datos es usando tablas de doble entrada
(Tabla 1).
Tabla 1. Valores de f calculados para distinto valores de Re (filas) y e/D
• Especifique los valores de Re, y e/D para los (columnas) (antes de ingreso de datos)
cuales desea calcular f como se muestra en la
Tabla 1.
• En la celda del vértice superior de la tabla
(O26 en ejemplo) ingrese “+nombre de la
Valores de e/D
celda que contiene la fórmula”. Esta celda es
la que contiene el valor de f luego de 10 Celda que contiene la
fórmula que depende
iteraciones (por ejemplo +H12) de e/D y Re
• Seleccione el rango de celdas de la tabla
(O26:T33 en el caso del ejemplo).
• En la barra de herramientas busque Valores de Re
Datos→Análisis de hipótesis →Tabla de
datos
• En el menú desplegable ingrese la celda que
contiene el valor de e/D en la sección datos
en Celda de entrada (fila)
• Seleccione la celda que contiene el valor de
Re en la sección datos para Celda de
entrada (columna).
• Para verificar que obtiene resultados
correctos asegúrese que la tabla contiene los
datos con los que realizó el cálculo. Valor de f calculado para e/D=1.4×10-4, Re=50000
Taller 2

Un intercambiador de calor que opera con flujos en contracorriente, las corrientes de

fluido caliente y fluido frío entran con temperaturas de 480°F (T1) y 120 °F (t1),
respectivamente.
Implemente en una planilla los cálculos para determinar el valor de las temperaturas de
salida de ambos fluidos (T2, t2) y el calor intercambiado Q entre las corrientes. Suponga
que el intercambiador se encuentra perfectamente aislado.

t2= ??? °F t1= 120 °F

T1= 480 °F T2= ??? °F

Figura 7. Diagrama de intercambiador de calor por contracorriente

Ecuaciones:

Q = WCpc (T1 − T2 ) = wCp f (t2 − t1 )

Del balance de energía térmica se tiene que: Q es el calor, W y w son flujos másicos,
Cp calor específico,T y t son las temperaturas de los fluidos

UA(T1 − t 2 ) − (T2 − t1 ) Ecuación de diseño del

Q= intercambiador de calor.
 (T1 − t 2 )
Ln  
(
 2 1 
T − t )
Temperatura promedio logarítmico.

Tml =
(T1 − t 2 ) − (T2 − t1 )
 (T − t )
Ln  1 2  Q = U * A * Tml
 (T2 − t1 )
Datos del problema:

Intercambiador de calor (A=662 pies2; U=55 BTU/(h pie2 °F)) que opera con
flujos en contracorriente, las corrientes de fluido caliente (W=43800 lb/h;
Cpc=0.508 BTU/(lb °F)) y fluido frío (w=1490000 lb/h; Cpf= 0.63 BTU/(lb °F))
entran con temperaturas de 480°F (T1) y 120 °F (t1), respectivamente.
a) Implemente en una planilla los cálculos para determinar el valor de las
temperaturas de salida de ambos fluidos (T2, t2) y el calor intercambiado Q
entre las corrientes. Suponga que el intercambiador se encuentra
perfectamente aislado. La planilla debe calcular la solución para distintos
valores de los flujos másicos de las corrientes, coeficiente de transferencia
de calor, área de transferencia, y propiedades de los fluidos.
b) Grafique el procedimiento de cálculo correspondiente a las 3 primeras
iteraciones*.
c) Utilice una tabla de doble entrada para determinar el efecto del área y el
coeficiente global de transferencia de calor.
Facultad de Ingeniería
Escuela de Ingeniería de Procesos Industriales

MÉTODOS DE RESOLUCIÓN DE
ECUACIONES NO LINEALES

París Paredes Inostroza

Claudia Aravena
Abril, 2021