TAREA 4 CASO PRACTICO 1

GRUPO: 2404
MATERIA: METODOS NUMERICOS
NOMBRE: ALBINO SOSA ALFREDO

Caso práctico

Planteamiento del problema

Corpogas, una empresa dedicada a venta de tanques de gas para uso particular dentro de la
Ciudad de México, desea conocer el volumen que deben tener los tanques de gas que distribuye.
Para ello considera que la presión a la que estará sometido el gas es de 200 kPa. Para tener esta
información firma un convenio con el Laboratorio de Caracterización de Mezclas del Instituto
Mexicano del Petróleo (IMP) ubicado en eje central, Ciudad de México.

Francisco, quien es ingeniero químico y responsable del laboratorio, realiza el experimento. El

experimento consiste en someter el gas a diferentes presiones y caracterizar el volumen de la
muestra. A partir del experimento obtiene los siguientes datos de presión y volumen:

A=3

Presión (kPa) 336 294.4 266.4 260.8 260.5 249.6 193.6 165.6

Volumen(cm3) 0.5+ 2 +A 3 +A 4 +A 6+A 8 +A 10 +A 11+A

A

Un detalle técnico al que se enfrenta Francisco, es que no cuenta con un dispositivo dentro de su
laboratorio que le permita fijar la presión del gas a 200 kPa por lo cual solo cuenta con los datos
de sus mediciones.

Francisco decide acudir al laboratorio de análisis de datos del instituto, donde la doctora Ximena
Ruíz, quien es especialista en análisis numérico, le dice que pondrá a alguien de sus becarios a
analizar los datos. Tú acabas de entrar a realizar tu servicio social y la doctora te asigna el análisis
de los datos de Francisco y te da las siguientes indicaciones para determinar el volumen del gas
cuando está sometido a una presión de 200 kPa.

a) Identificar las variables independiente y dependiente que describen el experimento de

Francisco.

b) Determinar dos modelos matemáticos diferentes que relacionen a la variable dependiente con
la variable independiente.
c) Diseñe un diagrama de bloques que permita relacionar las variables con cada uno de los
métodos que utilizó para obtener los modelos matemáticos.

Resultado de
presión 200 kPa
Polinomio

d) Generar en Excel un programa que te permita obtener la información asociada a cada uno de
los modelos matemáticos.

e) Con base en estos dos modelos matemáticos y tu programa en Excel, determinar el volumen de
la muestra cuando se encuentra sometida a una presión de 200 kPa.
Polinomio de Lagrange tomando 8 puntos

i Presión Volumen
0 336 3.5
1 294.4 5
2 266.4 6
3 260.8 7
4 260.5 9
5 249.6 11
6 193.6 13
7 165.6 14

7
(200 − 𝑥𝑗 )
𝐿(0) = ∏ = −0.0073
𝑖=0
(336 − 𝑥𝑗 )
𝑗≠0

7
(200 − 𝑥𝑗 )
𝐿(1) = ∏ = 1.7599
𝑖=1
(294.4 − 𝑥𝑗 )
𝑗≠1

7
(200 − 𝑥𝑗 )
𝐿(2) = ∏ = −64.9652
𝑖=2
(266.4 − 𝑥𝑗 )
𝑗≠2

7
(200 − 𝑥𝑗 )
𝐿(3) = ∏ = 1851.6709
𝑖=3
(260.8 − 𝑥𝑗 )
𝑗≠3

7
(200 − 𝑥𝑗 )
𝐿(4) = ∏ = −1805.3603
𝑖=4
(260.5 − 𝑥𝑗 )
𝑗≠4

7
(200 − 𝑥𝑗 )
𝐿(5) = ∏ = 18.4868
𝑖=5
(249.6 − 𝑥𝑗 )
𝑗≠5

7
(200 − 𝑥𝑗 )
𝐿(6) = ∏ = 0.7263
𝑖=6
(193.6 − 𝑥𝑗 )
𝑗≠6
 7
(200 − 𝑥𝑗 )
𝐿(7) = ∏ = −0.4234
𝑖=7
(165.6 − 𝑥𝑗 )
𝑗≠7

𝑓(200) = −0.0073(3.5) + 1.7599(5) − 64.9552(6) + 1851.6709(7) − 1805.3603(9)

+ 18.4868(11) + 0.7263(13) − 0.4234(14) = −3460.63335

Polinomio de Lagrange tomando 3 puntos

i Presión Volumen
0 249.6 11
1 193.6 13
2 165.6 14

2
(200 − 𝑥𝑗 )
𝐿(0) = ∏ = 0.0468
𝑖=0
(249.6 − 𝑥𝑗 )
𝑗≠0

2
(200 − 𝑥𝑗 )
𝐿(1) = ∏ = 1.0881
𝑖=1
(193.6 − 𝑥𝑗 )
𝑗≠1

2
(200 − 𝑥𝑗 )
𝐿(2) = ∏ = −0.1349
𝑖=2
(165.6 − 𝑥𝑗 )
𝑗≠2

𝑓(200) = 0.0468(11) + 1.0881(13) − 0.1349(14) = 12.7715

 Polinomio de Newton tomando 8 puntos

i Presión Volumen 1 2 3 4 5 6 7
0 336 3.5
1 294.4 5 -0.03605769
2 266.4 6 -0.03571429 -4.934E-06
-5.66042E-
260.8 7
3 -0.17857143 0.0042517 05
-
260.5 9
4 -6.66666667 1.09967716 0.032313435 0.00042724
- 3.91055E-
249.6 11
5 -0.18348624 -0.5788554 0.099912652 0.00295148 05
- -4.4065E- 5.84064E-
193.6 13
6 -0.03571429 0.00220885 -0.00858105 0.0014903 05 07
-2.32755E- -8.9893E- 1.56765E- -4.6383E- 6.14963E-
165.6 14
7 -0.03571429 0 05 05 05 07 09

a0 3.5
a1 -0.036057692
a2 -4.934E-06
a3 -5.66042E-05
a4 0.000427243
a5 3.91055E-05
a6 5.84064E-07
a7 6.14963E-09

𝑓(200) = 3.5 − 0.036(200 − 336) − 4.93𝑥10−6 (200 − 336)(200 − 294.4)

− 5.66𝑥10−5 (200 − 336)(200 − 294.4)(200 − 266.4)
+ 0.0004(200 − 336)(200 − 294.4)(200 − 266.4)(200 − 260.8)
+ 3.91𝑥10−5 (200 − 336)(200 − 294.4)(200 − 266.4)(200 − 260.8)(200
− 260.5)
+ 5.84𝑥10−7 (200 − 336)(200 − 294.4)(200 − 266.4)(200 − 260.8)(200
− 260.5)(200 − 249.6)
+ 5.84𝑥10−9 (200 − 336)(200 − 294.4)(200 − 266.4)(200 − 260.8)(200
− 260.5)(200 − 249.6)(200 − 193.6) = −3449.51
 Polinomio de Newton tomando 3 puntos

i Presión Volumen 1 2
0 249.6 11
1 193.6 13 -0.03571429
2 165.6 14 -0.03571429 0

a0 11
a1 -0.03571429
a2 0

𝑓 (200) = 11 − 0.0357 (200 − 249.6) + 0(200 − 249.6)(200 − 193.6) = 12.77072

d) De acuerdo con los valores encontrados en la literatura, el volumen del gas al someterse a una
presión de 200 kPa es de 12 cm3. Tomando como referencia este valor del volumen para el gas,
compara con los resultados que obtuviste con tus modelos matemáticos y determina cuál de ellos
es más conveniente utilizar.

Comparando los 2 modelos matemáticos para determinar cuál de ellos es más conveniente
utilizar, en ambos casos al interpolar en 8 puntos el volumen tiene un valor negativo por lo que se
utilizo 3 puntos en los cuales sí se logró obtener un volumen positivo con un error en ambos
modelos no tan significante, por lo que ambos son una buena opción, sin embargo por ser
prácticos el método de newton es más rápido si interpolamos por puntos no tan grandes.

A Francisco le interesa además caracterizar el dispositivo que utilizó para realizar el experimento,
en específico conocer el trabajo realizado durante el experimento.

Para determinar el trabajo realizado durante el experimento, la doctora Ximena te menciona que
la definición de trabajo termodinámico W está dada por la expresión:

𝑊 = − ∫ 𝑃𝑑𝑉

Donde la presión P está expresada como una función del volumen. Además, te menciona que el
número de moles de gas en el experimento fue de n= 125 moles. Por otro lado, te recomienda tratar
el gas como un gas ideal sometido a un proceso isotérmico. Al tomar en cuenta estas
consideraciones y mediante un desarrollo algebraico, se determina que el trabajo lo puedes calcular
con la expresión:
Donde R es la constante universal de los gases (R= 8.31 J/mol K). El experimento fue realizado a
temperatura ambiente (T=20° C= 293 K).

A partir de la definición de trabajo termodinámico, la doctora Ximena te asigna las siguientes

indicaciones:
e) Mediante tres métodos de integración diferentes y utilizando los datos que midió Francisco,
determina el trabajo realizado durante el experimento.

𝑛𝑅𝑇
𝑊 = −∫ 𝑑𝑣
𝑉

𝑅 = 8.31

𝑇 = 293

𝑛 = 125

304353.75
𝑊 = −∫ 𝑑𝑣
𝑉

f) Diseñe un diagrama de flujo que permita relacionar las variables con cada uno de los
métodos de integración que propone y que describa el algoritmo para determinar la solución.
g) Realiza un programa en Excel que te permita extraer la información necesaria para modelar
cada método de integración que utilizaste.
 Regla del trapecio múltiple

(4 segmentos)
14
304353.75
−∫ 𝑑𝑣
3.5 𝑉

i v f(v)
0 3.5 -86958
1 6.125 -49690
2 8.75 -34783
3 11.375 -26756
4 14 -21739

𝑓(𝑥0 ) + 2 ∑𝑛−1
𝑖=1 𝑓(𝑥𝑖 ) + 𝑓(𝑥𝑛 )
𝐼 = (𝑏 − 𝑎)
2𝑛

−86958 + 2(−49690 − 34783 − 26756) − 21739

𝐼 = (14 − 3.5)
2(4)

−86958 − 222458 − 21739

𝐼 = (10.5) = −434640.9375
8

Regla de Simpson 1/3 múltiple

(Aplicando la regla 3 veces)

i v f(v)

0 3.5 -86958

1 6.5 -46823

2 8 -38044

3 9.5 -32037

4 11 -27668

5 12.5 -24348

6 14 -21739
 𝑓(𝑥0 ) + 4 ∑𝑛−1 𝑛−2
𝑖=1,3,5 𝑓(𝑥𝑖 ) + 2 ∑𝑖=2,4,6 𝑓(𝑥𝑖 ) + 𝑓𝑥𝑛
I= (𝑏 − 𝑎)
3𝑛

−86958 + 4(−46823 − 32037 − 24348) + 2(−38044 − 27668) − 21739

I= (14 − 3.5)
3(6)

−86958 − 412832 − 131424 − 21739

I= (10.5) = −380889.25
18

Regla de Simpson 3/8

i v f(v)

0 3.5 -86958

1 7 -43479

2 10.5 -28986

3 14 -21739

𝑓 (𝑥0 ) + 3𝑓 (𝑥1 ) + 3𝑓 (𝑥2 ) + 𝑓(𝑥3 )

𝐼 = (𝑏 − 𝑎)
8

−86958 + 3(−43479) + 3(−28986) − 21739

𝐼 = (14 − 3.5)
8
−86958 − 130437 − 86958 − 21739
𝐼 = (10.5) = −427995.75
8