Universidad San Carlos de Guatemala

Centro Universitario de Occidente

División Ciencias de la Ingeniería
Matemática Aplicada 3
Ing. Luis Hernández

Proyecto De Aplicación “Métodos Numéricos En Ecuación De Van Der

Waals”

Pablo Andrés Morales Vidal 201530903

Quetzaltenango 8 de noviembre del 2,020
 Introducción

En el campo de la ingeniería y en cualquiera de sus ramas actualmente existen diferentes cantidades

de problemas que requieren soluciones o una aproximación a una solución real, es por eso que los
métodos numéricos han tomado muchísima importancia en el campo laboral debido a la factibilidad
de estos, ya que por medio de una gran cantidad de cálculos algebraicos permiten obtener datos que
se aproximen de manera muy precisa a la solución de estos problemas ya que estos se ayudan
conociendo el porcentaje de error aproximado en las iteraciones del cálculo, el uso de plataformas
que desarrollen la aplicación de estos tipo de cálculos ha aumentado exponencialmente con el paso
de los años ya que estos ayudan a facilitar el cálculo de los métodos numéricos, estos programas
pueden ser desarrollados de manera independiente o ser adquiridos, su desempeño hacen que sea una
herramienta confiable para múltiples aplicaciones.
 Objetivos

General
Realizar la aplicación de métodos numéricos en ejercicios de termodinámica.
Específicos

• Utilizar métodos numéricos abiertos, Newton-Raphson para calcular el volumen de

un gas aplicado a la ecuación de Van der Waals para gases ideales
• Determinar el volumen de un gas ideal aplicando los métodos numéricos cerrados,
Bisección en la ecuación de Van der Waals para gases ideales
• Comparar los resultados de los volúmenes calculados por los dos métodos Newton-
Raphson y Bisección con el resultado del volumen calculado con la ecuación que
proporciona el libro de termodinámica y conocer la exactitud de las diferentes
soluciones.
• Utilizar hojas de calculo de Excel para desarrollar el calculo de los dos métodos
numéricos utilizados.
Descripción de la investigación
Para el desarrollo del proyecto utilizamos un problema de gases ideales proporcionado por el libro de
termodinámica de Yunes Cengel, séptima edición, el cual debe de ser resuelto con la ayuda de una de
las primeras ecuaciones de estado, ecuación Van der Waals para gases ideales, esta ecuación al ser
operada de manera algebraica da un polinomio de tercer grado en términos de la variable que en este
caso es el volumen del gas y la mayor parte de los valores de las variables de la ecuación pueden ser
proporcionados por el problema o se pueden averiguar por medio de tablas que tiene el libro.
Se opto por aplicar el método abierto de Newton-Raphson por la velocidad de convergencia del
método y la exactitud al momento de dar una respuesta ya que traza una tangente del punto (𝑥𝑖 , 𝑓(𝑥𝑖 ))
de la curva de la función siendo esta una aproximación mejorada de la raíz.
El segundo método para solucionar el problema fue el método cerrado de Bisección ya que este
método aprovecha la característica de localizar un intervalo en el que la función cambie de signo.
Entonces, la localización del cambio de signo se logra con más exactitud al dividir el intervalo en
varios subintervalos. Se investiga cada uno de estos subintervalos para encontrar el cambio de signo.
El proceso se repite y la aproximación a la raíz mejora cada vez más.

El ejercicio dice lo siguiente.

Un dispositivo de cilindro-embolo contiene 100 gramos de monóxido de carbono. Inicialmente el
monóxido de carbono está a 1,000 Kpa y 200 ⁰C. Determine el volumen final del monóxido de
carbono, considerándolo
a. Como un gas ideal
b. La ecuación de Van der Waals
 Marco Teórico
Métodos Numéricos
Los métodos numéricos constituyen técnicas mediante las cuales es posible formular problemas
matemáticos de tal forma que puedan resolverse utilizando operaciones aritméticas. Siempre van a
requerir de múltiples cálculos aritméticos, por medio de diferentes tipos de Software ayuda a resolver
problemas no lineales de formas complejas y reduciendo tiempos emitiendo soluciones aproximadas.
“Chapra, S., Canale, R. Métodos numéricos para ingenieros. 5ta edición”

Tipos de Métodos numéricos

Métodos Cerrados
Son métodos que aprovechan el hecho de que una función cambia de signo en la vecindad de una
raíz. Estos métodos necesitan de dos valores iniciales para poder “encerrar” o estar a ambos lados de
la raíz, estos métodos emplean diferentes estrategias para reducir sistemáticamente el tamaño del
intervalo y así converger a la respuesta correcta.
“Chapra, S., Canale, R. Métodos numéricos para ingenieros. 5ta edición”

Método de Bisección
Conocido también como de corte binario, de partición de intervalos o de Bolzano, es un tipo de
búsqueda incremental en el que el intervalo se divide siempre a la mitad. Si la función cambia de
signo sobre un intervalo, se evalúa el valor de la función en el punto medio. La posición de la raíz se
determina situándola en el punto medio del subintervalo, dentro del cual ocurre un cambio de signo.
El proceso se repite hasta obtener una mejor aproximación.

“Chapra, S., Canale, R. Métodos numéricos para ingenieros. 5ta edición”

Métodos Abiertos
Se basan en fórmulas que requieren únicamente de un valor de inicio “x” o que empiecen con un par
de ellos, pero que no necesariamente encierren la raíz. Estos métodos divergen en algunas ocasiones
o se alejan de la raíz verdadera a medida que se avanza en el cálculo, tienden a converger con mayor
rapidez que los métodos cerrados
“Chapra, S., Canale, R. Métodos numéricos para ingenieros. 5ta edición”
Método de Newton-Raphson
Método abierto que se deduce a partir de la interpretación geométrica de un método alternativo basado
en la serie de Taylor. Se entiende que la primera derivada de “x” es equivalente a la pendiente
𝑓(𝑥𝑖 ) − 0
𝑓´(𝑥𝑖 ) =
𝑥𝑖 − 𝑥𝑖+1
Que se arregla para obtener la fórmula de Newton-Raphson
𝒇(𝒙𝒊 )
𝒙𝒊+𝟏 = 𝒙𝒊 −
𝒇´(𝒙𝒊 )
Aunque en general el método de Newton-Raphson es muy eficiente, hay situaciones donde se
comporta de manera deficiente
“Chapra, S., Canale, R. Métodos numéricos para ingenieros. 5ta edición”

Ecuaciones de estado

La ecuación de estado de gas ideal es muy simple, pero su ámbito de aplicabilidad es limitado. Es
deseable tener ecuaciones de estado que representen con precisión y sin limitaciones el
comportamiento P-v-T de las sustancias en una región más grande.
“Yunes, C., Michael, B. Termodinámica. 7ma edición”

Ecuación de estado de Van der Waals

La ecuación de estado de Van der Waals se propuso en 1873, tiene dos constantes que se
determinan del comportamiento de una sustancia en el punto crítico. Esta dada por

𝑎
(𝑃 + ) (𝑣 − 𝑏) = 𝑅𝑇
𝑣2

las fuerzas de atracción intermoleculares y el volumen que ocupan las moléculas por sí mismas. El
término a/v2 toma en cuenta las fuerzas intermoleculares y b el volumen que ocupan las moléculas
de gas. La determinación de las dos constantes que aparecen en esta ecuación se basa en la
observación de que la isoterma crítica en un diagrama P-v tiene un punto de inflexión horizontal en
el punto crítico

27𝑅 2 𝑇𝑐𝑟
2 𝑅𝑇𝑐𝑟
𝑎= 𝑦 𝑏=
64𝑃𝑐𝑟 8𝑃𝑐𝑟

“Yunes, C., Michael, B. Termodinámica. 7ma edición”

 Cálculos Realizados

Solución del problema con ecuaciones de estado para gas ideales

Datos
Gas = monóxido de carbono
P = 1000 Kpa
T= 200 ⁰C
m = 100g
a. Como un gas ideal

𝑃𝑉 = 𝑅𝑇, despejando para la variable del volumen tenemos que

𝑅𝑇
𝑉= para conocer el valor de la variable R debemos recurrir a las tablas del libro de
𝑃
termodinámica (Tabla A-1)
𝑚3
El valor de la variable R para monóxido de carbono es de 0.2968 𝐾𝑝𝑎 ∗
𝑘𝑔 ∗°𝐾
Convertimos los 200 ⁰C a ⁰K, ya que la constante R la temperatura es en ⁰K
°𝐾 = 2000 C + 273° = 473°K
Sustituimos valores en la ecuación para calcular el volumen

𝑚3
(0.2968𝐾𝑝𝑎∗ 𝑘 ∗ °𝐾)(473°𝐾)
𝑔 𝒎𝟑
𝑉= = 𝟎. 𝟏𝟒𝟎𝟑𝟖𝟔𝟒
1000 𝐾𝑝𝑎 𝑲𝒈

b. Ecuación de Van der Waals

Para calcular los coeficientes “a” y “b” necesitamos encontrar el valor de las variables

𝑃𝑐𝑟 𝑦 𝑇𝑐𝑟 y sustituirlos en las ecuaciones de los coeficientes.

Los valores de las variables Pcr = 3.50 Mpa y para Tcr = 133 ⁰K obtenidos por las tablas del
libro de termodinámica (Tabla A-1)

La presión critica del gas debe convertirse de Mpa a Kpa

1,000 𝐾𝑝𝑎
3.50 𝑀𝑝𝑎 ∗ = 3500 𝐾𝑝𝑎
1 𝑀𝑝𝑎

𝑚3
27(0.2968𝐾𝑝𝑎∗ )² (1330 𝐾)²
27𝑅2 𝑇𝑐𝑟
2 𝑘𝑔 ∗ °𝐾 𝑲𝒑𝒂
𝑎= = = 𝟎. 𝟏𝟖𝟕𝟖𝟐𝟐𝟏𝒎𝟔 ∗
64𝑃𝑐𝑟 64(3500 𝐾𝑝𝑎) 𝑲𝒈

𝑚3
(0.2968𝐾𝑝𝑎 ∗ 𝑘 ∗ °𝐾 ) (133 °𝐾)
𝑅𝑇𝑐𝑟 𝑔 𝒎𝟑
𝑏= = = 𝟎. 𝟎𝟎𝟏𝟒𝟎𝟗𝟖
8𝑃𝑐𝑟 8(3500 𝐾𝑝𝑎) 𝑲𝒈
 Ahora que ya se tienen calculados los valores de los coeficientes “a” y “b” se sustituyen en

la ecuación de Van der Waals para calcular el volumen del gas.

𝑎
(𝑃 + ) (𝑉 − 𝑏) = 𝑅𝑇
𝑉2

𝐾𝑝𝑎
0.1878221 𝑚6 ∗ 𝐾
𝑔 𝑚3 𝑚³
(1,000 𝐾𝑝𝑎 + ) (𝑉 − 0.0014098 ) = (0.2968 𝐾𝑝𝑎 ∗ )(473 °𝐾)
𝑉2 𝐾𝑔 𝐾𝑔 ∗ °𝐾

𝒎³
𝑽 = 𝟎. 𝟏𝟒𝟎𝟒𝟕𝟐𝟓 , solución con un error de 0.05639568 o del 5.639569%
𝑲𝒈

Solución para métodos numéricos

Primero comenzamos simplificando la formula para tener un polinomio de tercer grado en función
del volumen del gas.
𝑎 𝑎 𝑎𝑏
(𝑃 + ) (𝑉 − 𝑏) = 𝑅𝑇 → 𝑃𝑉 − 𝑃𝑏 + + − 𝑅𝑇 = 0
𝑉2 𝑉 𝑉2

Multiplicamos la ecuación resultante por V², el resultado de la ecuación multiplicada por V² la

dividimos por 1/P.

1
𝑃𝑉 3 − 𝑃𝑏𝑉 2 + 𝑎𝑉 + 𝑎𝑏 − 𝑅𝑇𝑉 2 = 0, 𝑑𝑖𝑣𝑖𝑑𝑖𝑚𝑜𝑠 𝑝𝑜𝑟
𝑃

𝑎𝑉 𝑎𝑏 𝑅𝑇𝑉²
𝑉 3 − 𝑏𝑉 2 + + −
𝑃 𝑃 𝑃
Factorizamos V² de la ecuación y ordenamos para obtener la ecuación que será utilizada en el
método de Newton-Raphson y Müller

𝑷𝒃 + 𝑹𝑻 𝟐 𝒂 𝒂𝒃
𝑽𝟑 − ( )𝑽 + ( )𝑽 − =𝟎
𝑷 𝑷 𝑷
O bien sustituyendo valores tenemos la siguiente ecuación

𝑽𝟑 − 𝟎. 𝟏𝟒𝟏𝟕𝟗𝟔𝟐𝑽𝟐 + 𝟎. 𝟎𝟎𝟎𝟏𝟖𝟖𝑽 − 𝟎. 𝟎𝟎𝟎𝟎𝟎𝟎𝟓𝟓𝟗𝟖𝟏 = 𝟎

Solución del problema utilizando el método de Newton-Raphson
Tenemos la función

𝑓(𝑉) = 𝑉 3 − 0.1417962𝑉 2 + 0.000188𝑉 − 0.00000055981

𝑓´(𝑉) = 3𝑉 2 − 0.2835924𝑉 + 0.000188

Implementamos el método de Newton-Raphson de la siguiente manera
𝒇(𝒙𝒊 )
𝒙𝒊+𝟏 = 𝒙𝒊 −
𝒇´(𝒙𝒊 )

𝑉³ − 0.1417962𝑉 2 + 0.000188𝑉 − 0.00000055981

𝑉𝑖+1 = 𝑉𝑖 −
3𝑉 2 − 0.2835924𝑉 + 0.000188
Primera Iteración

𝑉𝑖 = 0.1404725
(0.1404725)³ − 0.1417962(0.1404725)2 + 0.000188(0.1404725) − 0.00000055981
𝑉𝑖+1 = (0.1404725) −
3(0.1404725)2 − 0.2835924(0.1404725) + 0.000188

𝑽𝒏+𝟏 = 𝟎. 𝟎𝟖𝟒𝟎𝟕𝟕 ← 𝑹𝒂𝒊𝒛 𝑴𝒆𝒋𝒐𝒓𝒂𝒅𝒂

Error de la primera iteración

𝐸𝑟𝑟𝑜𝑟 = |𝑉𝑖+1 − 𝑉𝑖 |
𝐸𝑟𝑟𝑜𝑟 = |0.08407682 − 0.1404725| = 𝟎. 𝟎𝟓𝟔𝟑𝟗𝟓𝟔𝟖
Seguimos calculando ya que el Error debe ser
𝐸𝑟𝑟𝑜𝑟 ≤ 1 ∗ 10−6
Segunda Iteración

𝑉𝑖 = 0.084077
(0.08407682)³ − 0.1417962(0.08407682)2 + 0.000188(0.08407682) − 0.00000055981
𝑉𝑖+1 = (0.084077) −
3(0.08407682)2 − 0.2835924(0.08407682) + 0.000188

𝑽𝒏+𝟏 = 𝟎. 𝟎𝟒𝟖𝟒𝟒𝟒𝟗𝟗 ← 𝑹𝒂𝒊𝒛 𝑴𝒆𝒋𝒐𝒓𝒂𝒅𝒂

Error de la primera iteración

𝐸𝑟𝑟𝑜𝑟 = |𝑉𝑖+1 − 𝑉𝑖 |
𝐸𝑟𝑟𝑜𝑟 = |0.04844499 − 0.08407682| = 𝟎. 𝟎𝟑𝟓𝟔𝟑𝟏𝟖𝟑
Seguimos calculando ya que el Error debe ser
𝐸𝑟𝑟𝑜𝑟 ≤ 1 ∗ 10−6
Tercera Iteración

𝑉𝑖 = 0.0484499
(0.0484499)³ − 0.1417962(0.0484499)2 + 0.000188(0.0484499) − 0.00000055981
𝑉𝑖+1 = (0.0484499) −
3(0.0484499)2 − 0.2835924(0.0484499) + 0.000188

𝑽𝒏+𝟏 = 𝟎. 𝟎𝟐𝟔𝟕𝟒𝟑𝟐𝟗 ← 𝑹𝒂𝒊𝒛 𝑴𝒆𝒋𝒐𝒓𝒂𝒅𝒂

Error de la primera iteración

𝐸𝑟𝑟𝑜𝑟 = |𝑉𝑖+1 − 𝑉𝑖 |
𝐸𝑟𝑟𝑜𝑟 = |0.02674329 − 0.04844499| = 𝟎. 𝟎𝟐𝟏𝟕𝟎𝟏𝟕
Seguimos calculando ya que el Error debe ser
𝐸𝑟𝑟𝑜𝑟 ≤ 1 ∗ 10−6
Seguimos calculando hasta que se cumpla la condición del error llegando hasta la undécima
iteración del problema, ya que en esta iteración se cumple la condición
VALORES INICIALES.

Cálculos desarrollados por una hoja de Excel programada para desarrollar el método
El volumen del monóxido de carbono calculado con el método de Newton-Raphson es de
0.00142615 m³/Kg con un error de 0%
Solución del problema utilizando el método de Bisección
Partimos de la misma ecuación

𝑓(𝑉) = 𝑉 3 − 0.1417962𝑉 2 + 0.000188𝑉 − 0.00000055981

Para utilizar este método de la bisección realizaremos la siguiente tabla

Implementamos el método de Bisección

𝑋𝑎 +𝑋𝑏
𝑋𝑟 = ; 𝑓(𝑋𝑟 ); 𝑓(𝑋𝑎 ); 𝑓(𝑋𝑏 ) 𝑦 𝑓 (𝑋𝑎 ) ∗ 𝑓 (𝑋𝑟 )
2

Como
𝑉𝑎 + 𝑉𝑏
𝑉𝑟 = ; 𝑓(𝑉𝑟 ); 𝑓(𝑉𝑎 ); 𝑓(𝑉𝑏 ) 𝑦 𝑓(𝑉𝑎 ) ∗ 𝑓 (𝑉𝑟 )
2

Primera Iteración
Tomamos como valores iniciales 𝑉𝑎 = 0 𝑦 𝑉𝑏 = 0.5
0 + 0.5
𝑉𝑟 = = 𝟎. 𝟐𝟓
2
𝑓(𝑉𝑟 ) = (0.25)3 − 0.1417962(0.25)2 + 0.000188(0.25) − 0.00000055981 = 0.024533703

𝑓(𝑉𝑎 ) = (0)3 − 0.1417962(0)2 + 0.000188(0) − 0.00000055981 = −0.00000055981

𝑓(𝑉𝑏 ) = (0.5)3 − 0.1417962(0.5)2 + 0.000188(0.5) − 0.00000055981 = 0.16054249

𝑓(𝑉𝑎 ) ∗ 𝑓(𝑉𝑟 ) = (−0.00000055981)(0.0245337) = −𝟏. 𝟒𝑬 − 𝟖

𝑉𝑟 (𝑎𝑐𝑡𝑢𝑎𝑙) − 𝑉𝑟 (𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟)
𝐸𝑝 = | | ∗ 100
𝑉𝑟 (𝑎𝑐𝑡𝑢𝑎𝑙)

0.25 − 0
𝐸𝑝 = | | 100 = 𝟏𝟎𝟎%
0.25
Seguimos calculando ya que el Error debe ser
𝐸𝑟𝑟𝑜𝑟 ≤ 1 ∗ 10−6
Ya que el valor de f(Va)* f(Vr) no es mayor que 0 nuestros valores iniciales para la siguiente
iteración son los siguientes, el nuevo valor de Vr se debe sustituir en Vb.
Segunda Iteración

Tomamos como valores iniciales 𝑉𝑎 = 0 𝑦 𝑉𝑏 = 0.25

0 + 0.25
𝑉𝑟 = = 𝟎. 𝟏𝟐𝟓
2
𝑓(𝑉𝑟 ) = (0.125)3 − 0.1417962(0.125)2 + 0.000188(0.125) − 0.00000055981 = 0.004191631

𝑓(𝑉𝑎 ) = (0)3 − 0.1417962(0)2 + 0.000188(0) − 0.00000055981 = −0.00000055981

𝑓(𝑉𝑏 ) = (0.25)3 − 0.1417962(0.25)2 + 0.000188(0.25) − 0.00000055981 = 0.024533703
𝑓(𝑉𝑎 ) ∗ 𝑓(𝑉𝑟 ) = (−0.00000055981)(0.0245337) = −𝟐𝑬 − 𝟎𝟗
𝑉𝑟 (𝑎𝑐𝑡𝑢𝑎𝑙) − 𝑉𝑟 (𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟)
𝐸𝑝 = | | ∗ 100
𝑉𝑟 (𝑎𝑐𝑡𝑢𝑎𝑙)

0.125 − 0.25
𝐸𝑝 = | | 100 = 𝟏𝟐. 𝟓%
0.125
Ya que el valor de f(Va)* f(Vr) no es mayor que 0 nuestros valores iniciales para la siguiente
iteración son los siguientes, el nuevo valor de Vr se debe sustituir en Vb.

Seguimos calculando ya que el Error debe ser

𝐸𝑟𝑟𝑜𝑟 ≤ 1 ∗ 10−6
Tercera Iteración

Tomamos como valores iniciales 𝑉𝑎 = 0 𝑦 𝑉𝑏 = 0.125

0 + 0.125
𝑉𝑟 = = 𝟎. 𝟎𝟔𝟐𝟓
2
𝑓(𝑉𝑟 ) = (0.0625)3 − 0.1417962(0.0625)2 + 0.000188(0.0625) − 0.00000055981 = 0.000809222

𝑓(𝑉𝑎 ) = (0)3 − 0.1417962(0)2 + 0.000188(0) − 0.00000055981 = −0.00000055981

𝑓(𝑉𝑏 ) = (0.25)3 − 0.1417962(0.25)2 + 0.000188(0.25) − 0.00000055981 = 0.004191631

𝑓(𝑉𝑎 ) ∗ 𝑓(𝑉𝑟 ) = (−0.00000055981)(0.0245337) = 𝟒. 𝟓𝟑𝑬 − 𝟏𝟎

𝑉𝑟 (𝑎𝑐𝑡𝑢𝑎𝑙) − 𝑉𝑟 (𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟)
𝐸𝑝 = | | ∗ 100
𝑉𝑟 (𝑎𝑐𝑡𝑢𝑎𝑙)

0.0625 − 0.125
𝐸𝑝 = | | 100 = 𝟔. 𝟐𝟓%
0.0625
Seguimos calculando ya que el Error debe ser
𝐸𝑟𝑟𝑜𝑟 ≤ 1 ∗ 10−6
Seguimos calculando hasta que se cumpla la condición del error llegando hasta la decimotercera
iteración del problema, ya que en esta iteración se cumple la condición del error
VALORES INICIALES.

Cálculos desarrollados por una hoja de Excel programada para desarrollar el método
El volumen del monóxido de carbono calculado con el método de Bisección es de 0.00140381
m³/Kg con un error de 6.1035E-05%
 Conclusiones
Como pudimos apreciar en la solución del problema de gas ideal utilizando diferentes métodos para
resolverlo concluimos lo siguiente:
1. La solución por medio de la ecuación de Van der Vaals tiene un error del 5.639569% es un
error aceptable ya que el libro de termodinámica explica la inexactitud de la ecuación, la
razón por la cual se utilizó este tema para la aplicación de métodos numéricos.
2. El numero de iteraciones entre los métodos utilizados afecta al momento de encontrar una
solución en el problema ya que el método cerrado necesito de trece iteraciones y el abierto
once.
3. El resultado del método de Newton-Raphson es más exacto comparándolo con el resultado
del método de Bisección ya que a la hora de que se encontró una “raíz” en los dos métodos
el error del método de Newton fue del 0% y el error del método de Bisección cuando
encontramos la “raíz” es del 6.1035E-05%, notándose esto en la solucion de los métodos ya
que al momento de hacer la diferencia entre las dos, tenemos una diferencia del 0.00002234
m³/Kg
4. Los métodos numéricos abiertos convergen con menor cantidad de iteraciones que los
métodos cerrados y tienden a ser mas exactos al momento de encontrar una solución o raíz.
5. El uso de hojas de calculo o diferentes softwares ayudan a reducir el tiempo de realizar las
múltiples operaciones aritméticas que se deben hacer para encontrar una solución al
problema, llevando de la mano la exactitud en los cálculos que estos ofrecen.
 Bibliografía

Chapra, S., Canale, R. (2007). Métodos numéricos para ingenieros. 5ta edición
Pps. 3, 120, 124, 142, 148 y 149
Yunes, C., Michael, B. (2012). Termodinámica. 7ma edición
Pps. 144, 145, 160 y 908