8 de septiembre de 2023

Instituto Tecnológico de Durango

Ingeniería química

Herramientas computacionales para ingenieros químicos

UNIDAD II

POLYMATH

Alumna: García Rodríguez Yeimi Alejandra

No. de control: 20041167

Docente: Pope Solís José Domingo

Grupo: 7WP
 Índice
Ejercicio 1. Volumen molar y factor de compresibilidad usando la ecuación de Van Der Waals .......2
Ejercicio 2. Cálculo de concentraciones ............................................................................................6
Ejercicio 3. Balances de materia en estado estacionario en un tren de separación. ............................8
Ejercicio 4. Matrices y sistemas de ecuaciones lineales ..................................................................10
Ejercicio 5. Tres ecuaciones diferenciales ordinarias simultáneas y dos ecuaciones algebraicas
explícitas………………………………………………………………………………………………………………………………………12
Ejercicio 6. Sistema amortiguado resorte-masa ...............................................................................15
Ejercicio 7. Ajuste de los datos de calibración de un caudalímetro a una recta ...............................18
Ejercicio 8. Ajuste de datos no lineales a una recta .........................................................................21
Ejercicio 9. Ajuste de datos no lineales a una recta .........................................................................24
Ejercicio 10. Ajuste de datos no lineales a una recta .......................................................................27

1
Ejercicio 1. Volumen molar y factor de compresibilidad usando la
ecuación de Van Der Waals
Conceptos demostrados
Usar la ecuación de Van Der Waals para calcular el volumen molar y el factor de
compresibilidad de un gas.

Método numérico utilizado

Ecuación algebraica no lineal.

Problema
La ley de los gases ideales puede representar la relación presión-volumen-temperatura (PVT)
de los gases sólo a presiones bajas (cerca de la atmosférica). Para presiones más altas, se
deben utilizar ecuaciones de estado más complejas. El cálculo del volumen molar y el factor
de compresibilidad utilizando ecuaciones de estado complejas típicamente requiere una
solución numérica cuando se especifican la presión y la temperatura.
La ecuación de Van Der Waals está dada por:
𝑎
(𝑃 + ) (𝑉 − 𝑏) = 𝑅𝑇
𝑉2
Donde:
27 𝑅 2 𝑇𝑐2
𝑎= ( )
64 𝑃𝑐

𝑅𝑇𝐶
𝑏=
8𝑃𝐶

Las variables son definidas por:

P= Presión en atm

V= Volumen molar en L/g·mol

T= Temperatura en K
R= Constante de los gases (R = 0.08206 atm·L / g-mol·K)

Tc= Temperatura crítica (405.5 K para amoniaco)

Pc= Presión crítica (111.3 atm para amoniaco)

2
Presión reducida está dada por:
𝑃
𝑃𝑟 =
𝑃𝑐

Y el factor de compresibilidad está dado por:

𝑃𝑉
𝑍=
𝑅𝑇
A) Calcular el volumen molar y el factor de compresibilidad de los gases amonia a una
presión de 56 atm y a una temperatura de 450 K usando la ecuación de Van Der Waals.

Solución
La ecuación (2-1) no se puede reorganizar en una forma en la que V pueda expresarse
explícitamente, expresado como una función de T y P. Sin embargo, se puede resolver
fácilmente numéricamente utilizando técnicas para ecuaciones no lineales.
Para resolverlo, la ecuación de Van Der Waals debe de estar reescrita de la siguiente manera:
𝑎
𝑓(𝑉) = (𝑃 + ) (𝑉 − 𝑏) − 𝑅𝑇
𝑉2
Donde la solución se obtiene cuando la función es cercana a cero, 𝑓(𝑉) = 0.

Para resolverte este ejercicio debemos de tomar en cuenta:

1. Tomar en cuenta los siguientes datos y ecuaciones

Line Ecuación
1 f(V)=(P+a/(V^2))*(V-b)-
R*T
2 P=56
3 R=0.08206
4 T=450
5 Tc=405.5
6 Pc=111.3
7 Pr=P/Pc
8 a=27*(R^2*Tc^2/Pc)/64
9 b=R*Tc/(8*Pc)
10 Z=P*V/(R*T)
11 V(min)=0.4
12 V(max)=1
2. Agregar la función de V. Para la parte de este problema, el volumen es calculado a
partir del gas ideal.

3
 3. Para al final agregar las demás ecuaciones y datos numéricos.

Resultados

4
5
Ejercicio 2. Cálculo de concentraciones
Conceptos demostrados
Usar las ecuaciones de velocidad de reacción para determinar concentraciones y volumen.

Método número utilizado

Sistema de ecuaciones algebraicas no lineales.

Problema
Considerando las siguientes ecuaciones:
2
𝑘 ∗ 𝐶𝐴1 = 𝑣 ∗ (𝐶𝐴0 − 𝐶𝐴1 )/𝑉
2
𝑘 ∗ 𝐶𝐴2 = 𝑣 ∗ (𝐶𝐴1 − 𝐶𝐴2 )/𝑉

Donde:

k = 0.075
v = 30
CA0 = 1.6

CAQ = 0.2 CA0

Solución
Para resolver este sistema de ecuaciones, cada ecuación algebraica no lineal debe reescribirse
en la forma f(x) = (una expresión que debe tener el valor de cero en la solución). Las formas
apropiadas para estas ecuaciones son:
2
𝑓(𝐶𝐴1 ) = 𝑘 ∗ 𝐶𝐴1 − 𝑣 ∗ (𝐶𝐴0 − 𝐶𝐴1 )/𝑉
2
𝑓(𝑉) = 𝑘 ∗ 𝐶𝐴2 − 𝑣 ∗ (𝐶𝐴1 − 2)/𝑉

Dentro de del programa de Polymath los datos deben ser escritos de manera correcta,
respetando las variables y sus valores.

6
Por último, tener la ventana principal con las ecuaciones y datos correspondientes.

Resultados

7
Ejercicio 3. Balances de materia en estado estacionario en un tren de
separación.
Conceptos demostrados
Balances de materia en un proceso en estado estacionario sin reciclaje.

Método numérico utilizado

Solución de ecuaciones lineales simultáneas.

Problema
El paraxileno, el estireno, el tolueno y el benceno deben separarse con la matriz de columnas
de destilación como se muestra en la siguiente imagen.

Calcular los caudales molares de D1, D2, B1 y B2.

Solución
Determinar el rendimiento de los balances de materiales sobre componentes individuales.
𝑃𝑎𝑟𝑎𝑥𝑖𝑙𝑒𝑛𝑜 = 0.07𝐷1 + 0.18𝐵1 + 0.15𝐷2 + 0.24𝐵2 = 0.15 ∗ 70

𝐸𝑠𝑡𝑖𝑟𝑒𝑛𝑜 = 0.04𝐷1 + 0.24𝐵1 + 0.10𝐷2 + 0.65𝐵2 = 0.25 ∗ 70

𝑇𝑜𝑙𝑢𝑒𝑛𝑜 = 0.54𝐷1 + 0.42𝐵1 + 0.54𝐷2 + 0.10𝐵2 = 0.40 ∗ 70

𝐵𝑒𝑛𝑐𝑒𝑛𝑜 = 0.35𝐷1 + 0.16 𝐵1 + 0.21𝐷2 + 0.01𝐵2 = 0.20 ∗ 70

8
Los coeficientes y las constantes de estas ecuaciones se pueden introducir directamente en
POLYMATH Linear Equation Solver en forma matricial de la siguiente manera:

Resultados

9
Ejercicio 4. Matrices y sistemas de ecuaciones lineales
Conceptos demostrados
Obtener los valores x1, x2, x3.

Método numérico utilizado

Sistemas de ecuaciones lineales.

Problema
Con los datos del diagrama siguiente (donde los porcentajes están dados en peso), encuentre
posibles valores de las corrientes M1, M2, M3 y M4.

Solución
Mediante balances de materia por componente y global, se tiene:

Solo se tienen tres ecuaciones linealmente independientes, pues la ecuación del balance
global de materia es la suma de las otras tres. Por el sistema homogéneo es consistente, y
como el rango de la matriz coeficiente es menor que el número de incógnitas, el sistema tiene
un número infinito de soluciones. Fijando una base de cálculo, por ejemplo M4 = 100 kg, se
obtiene el sistema.

10
Resultados

11
Ejercicio 5. Tres ecuaciones diferenciales ordinarias simultáneas y
dos ecuaciones algebraicas explícitas
Conceptos demostrados
Calcular los valores de las variables.

Método numérico utilizado

Runge-Kutta45.

Problema
Considerar las siguientes ecuaciones:
𝑑(𝐴)
= −𝐾1 𝐴
𝑑(𝑡)
𝑑(𝐵)
= 𝐾1 𝐴 + 𝐾2 𝐵
𝑑(𝑡)
𝑑(𝐶)
= 𝐾2 𝐵
𝑑(𝑡)

Que se integra desde t = 0 hasta t = 3, donde:

A (0) = 1

B (0) = 0
C (0) = 0 (valores iniciales)
Siendo los valores constantes: k1=1 y k2 = 2.

Solución
Introducir cada ecuación diferencial con su correspondiente valor inicial.

12
Después agregar los valores de las constantes.

Asigne un valor a k2 de manera similar. Antes de que se pueda llevar a cabo la integración,
también se deben especificar los valores inicial y final de la variable independiente. A
continuación, se muestra la forma final del sistema de ecuaciones con los valores iniciales y
finales apropiados, listo para ser resuelto.

13
Resultados

14
Ejercicio 6. Sistema amortiguado resorte-masa
Conceptos demostrados
Resolver a partir de ecuaciones diferencial.

Método numérico utilizado

Runge-Kutta45.

Problema
La ecuación diferencial ordinaria siguiente describe el movimiento de un sistema
amortiguado resorte-masa.
𝑑2𝑥 𝑑𝑥 𝑑𝑥
𝑚2 2 +𝑎| | + 𝑏𝑥 3 = 0
𝑑𝑡 𝑥𝑡 𝑑𝑡
Donde x = desplazamiento a partir de la posición de equilibrio, t = tiempo, m = 1 kg masa, y
a = 5 N/(m/s)2. El término de amortiguamiento es no lineal y representa el amortiguamiento
del aire.
El resorte es un resorte cúbico y también es no lineal con b = 5 N/m3. Las condiciones
iniciales son:
𝑑𝑥
Velocidad inicial 𝑑𝑡
= 0.5 𝑚/𝑠

Desplazamiento x = 1

El caso siguiente es una ecuación lineal similar:

𝑑2𝑥 𝑑𝑥
𝑚 2
+2 + 5𝑥 = 0
𝑑𝑡 𝑑𝑡
Resolver esta ecuación con algún método numérico para el periodo de tiempo 0 ≤ t ≤ 8 s.

Solución
Convertir la ecuación diferencial anterior en una ecuación diferencial ordinaria.
𝑑 𝑑𝑥 𝑑𝑥
( )+2 + 5𝑥 = 0
𝑑𝑡 𝑑𝑡 𝑑𝑡
𝑥 = 𝑥1
𝑑𝑥1
= 𝑥2
𝑑𝑡
𝑑𝑥2
+ 2𝑥2 + 5𝑥1 = 0
𝑑𝑡

15
 𝑑𝑥2
= −2𝑥2 − 5𝑥1
𝑑𝑡
𝑑𝑥1
𝑡=0 = 0.5
𝑑𝑡
𝑡=0 𝑥=1

𝑥2 = 0.5 𝑥1 = 1

Añadir las ecuaciones al programa.

De la misma forma agregar las ecuaciones, valores y límites.

16
Resultados

17
Ejercicio 7. Ajuste de los datos de calibración de un caudalímetro a
una recta
Conceptos demostrados
Obtener los datos de las variables.

Método numérico utilizado

Regresión lineal

Problema
Se obtuvieron los datos siguientes para la calibración de un rotámetro (velocidad de flujo
contra lectura del rotámetro):

Lectura del rotámetro

Velocidad de flujo V(L/min)
R
20 10
52.1 30
84.6 50
118.3 70
151.0 80

Solución
1. Agregar los datos de las velocidades de flujo en las celdas.

2. Agregar los datos de la lectura del rotámetro.

18
 3. Definir cuál es la variable dependiente e independiente

Resultados

19
20
Ejercicio 8. Ajuste de datos no lineales a una recta
Conceptos demostrados
Obtener los datos de las variables.

Método numérico utilizado

Regresión no lineal.

Problema
Se mide la velocidad de flujo de la masa m(g/s) como función de la temperatura T°C.

T 10 20 40 80
m 14.76 20.14 27.73 38.47

Hay motivos para creer que m varía linealmente con la raíz cuadrada de T:

m = a T ^ (1/2) + b
Use una gráfica lineal para verificar esta fórmula y determinar los valores de a y b.

Solución
Si la fórmula es correcta, la gráfica de m contra T ^ (1/2) sería lineal, con pendiente = a e
intersección = b. A continuación, se agrega la fila de T ^ (1/2) a la tabla de datos:

Y se grafica m contra T ^ (1/2).

21
Agregar la ecuación correspondiente.

Resultados

22
23
Ejercicio 9. Ajuste de datos no lineales a una recta
Conceptos demostrados
Obtener los datos de las variables a, b y c.

Método numérico utilizado

Regresión no lineal.

Problema
Se cree que la lectura de un instrumento de proceso, Z(volts), se relaciona con la velocidad
de flujo de la corriente de proceso F(L/s) y la presión P(kPa) por la siguiente expresión:
𝑍 = 𝑎𝑉 𝑏 𝑃𝑐

Se obtuvieron datos del proceso en dos conjuntos de corridas —en una se mantuvo V
constante y en la otra P constante—. Los datos obtenidos fueron:

Punto 1 2 3 4 5 6 7
V (Ls) 0.65 1.02 1.75 3.43 1.02 1.02 1.02
P(kPa) 11.2 11.2 11.2 11.2 9.1 7.6 5.4
Z (volts) 2.27 2.58 3.72 5.21 3.50 4.19 5.89

Solución
1. Se agregan los datos de la tabla anterior dentro de Polymath.

2. Linealizar la ecuación tomando en cuenta:

𝑦 = 𝑎0 + 𝑎1 𝑥1 + 𝑎2 𝑥2

𝑙𝑛𝑍 = 𝑙𝑛𝑎 + 𝑏𝑙𝑛𝑉 + 𝑐𝐿𝑛𝑃

3. Añadir X1 y X2 y Y en relación con su logaritmo natural correspondiente.

24
 4. Dentro del apartado regresión no lineal agregar la expresión 𝑍 = 𝑎𝑉 𝑏 𝑃𝑐 .

Resultados

25
26
Ejercicio 10. Ajuste de datos no lineales a una recta
Conceptos demostrados
Obtener los datos de las variables a, b y c.

Método numérico utilizado

Regresión polinomial.

Problema
Un tanque de almacenamiento se carga con una solución que contiene desechos peligrosos y
ésta se somete a un tratamiento químico para descomponer los desechos en productos
inofensivos. Se ha visto que la concentración del desecho que se descompone. C, varía con
el tiempo de acuerdo con la fórmula:

𝐶 = 1/(𝑎 + 𝑏𝑡)

Cuando ha transcurrido el tiempo suficiente para que la concentración descienda a 0.01 g/L,
el contenido del tanque se descarga en un río que pasa junto a la planta.
Se obtienen los siguientes datos para C y t:

t (h) 1.0 2.0 3.0 4.0 5.0

C (g/L) 1.43 1.02 0.73 0.53 0.38

Solución
1. Agregar los datos a polymath.

2. Linealizar la formula
1
= 𝑎 + 𝑏𝑡
𝐶
3. Con base a lo anterior determinar el valor de y con base a 1/C.

27
 4. Siendo x el valor de t.

5. Establecer 3er grado polinómico.

Resultados

28
29