UNIVERSIDAD NACIONAL DE INGENIERÍA

FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL

MÉTODOS NUMÉRICOS PARA INGENIERÍA QUÍMICA (PI-524/A)
SEMESTRE 2020 - I

3° Práctica Calificada

DATOS
APELLIDOS Y NOMBRES: Rey de Castro Gutierrez, Daniel Alejandro CÓDIGO: 20180111E

INDICACIONES
1. Grabar su práctica con sus apellidos
2. La solución del problema se debe de colocar en orden
3. Detallar claramente el procedimiento seguido

1. Para m = 3 gr en un caso de transferencia de calor se tiene la siguiente ecuación

m ⋅ Cp ⋅ ((T - 20)) ＝ Q - Qp

Para los datos de la tabla adjunta se desea encontrar los parámetros de los siguientes
modelos propuestos:

A) Qp ＝ a + a ⋅ T + a ⋅ T 2 B) Qp ＝ b + 2.31 ⋅ T + b ⋅ T 2
0 1 2 0 2

2. Sea una sustancia pura cuya presión de vapor Pv puede calcularse según la expresión
Frost Kalkwarf Thodos (F K T )

B D ⋅ Pv
ln ((Pv)) ＝ A - ―+ C ⋅ ln ((T)) + ―― . . . . . . . . . . . . . . (1)
T T2

A) Utilizando los criterios de regresión determinar A, B, C y D para los

datos mostrados en la tabla adjunta donde la constante de equilibrio ( Ks )
depende de la temperatura (T en °C)

Donde:

Los datos de la tabla se tomaron a:

P ＝ 0.95 ⋅ Atm
t

Pv
K ＝ ――
s P
t

B) Considerando los datos de la tabla indicar si la ecuación (2) es mejor a

la ecuación ( (1)

B
ln ((Pv)) ＝ A - ―― . . . . . . . . . . . . . . (2)
T+C

APELLIDOS Y NOMBRES: Rey de Castro Gutierrez, Daniel Alejandro 3 Practica Calificada

 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
MÉTODOS NUMÉRICOS PARA INGENIERÍA QUÍMICA (PI-524/A)
SEMESTRE 2020 - I

3. Determinar
W＝a + b + c + d + e
Donde:
a, b, c, d, e son números enteros

Si se sabe que al aplicar el método de aproximación polinomial simple con los 4 primeros
puntos de la tabla adjunta se encontró un polinomio de interpolación P(x) de tercer
grado .
Para los datos ya mencionados se cumple la siguiente ecuación matricial:

Además se cumple también:

⎛5 b⎞
f ⎜―― ⎟ ＝ 10
⎝ 2 ⎠

4. mol
La tabla adjunta muestra los valores de la concentración (en ―― ) de una sustancia a
L
diferentes valores de presión P y temperatura T. Utilizando Lagrange con 3 puntos

Kg
A) Encuentre el valor de la concentración a la presión de 3.8 ⋅ ――y T = 172 K
cm 2

Kg mol
B) Trabajando a 4 ―― hallar a que temperatura se cumple: C ＝ 0.117 ⋅ ――
2
cm L

NOTA
En todos los ejercicios comente sus resultados

Solución
APELLIDOS Y NOMBRES: Rey de Castro Gutierrez, Daniel Alejandro 3 Practica Calificada
 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
MÉTODOS NUMÉRICOS PARA INGENIERÍA QUÍMICA (PI-524/A)
SEMESTRE 2020 - I

Solución

Pregunta 1

Se tiene la siguiente ecuación para un caso de transferencia de calor:

m ⋅ Cp ⋅ ((T - 20)) ＝ Q - Qp m≔3

Despejando: Qp ＝ Q - m ⋅ Cp ⋅ ((T - 20))

Se tienen los siguientes datos de un caso de transferencia de calor:

T Q Cp ⎡ 142.08 ⎤
⎢ 201.25 ⎥
―――――→ ⎢ 231.17 ⎥
25 150 0.528 Qp ≔ Q - m ⋅ Cp ⋅ ((T - 20)) = ⎢ ⎥
⎢ 270.74 ⎥
30 220 0.625 ⎢ 322.425 ⎥
⎢⎣ 342.54 ⎥⎦
35 266 0.774
40 320 0.821
Se coloca vectorización para que se multiplique
45 390 0.901
elemento a elemento.
50 432 0.994

Entonces, se desea plantear los siguientes modelos:

a) Qp ＝ a + a ⋅ T + a ⋅ T 2 Se identifica que es un caso de regresión

0 1 2
polinomial de segundo orden respecto a T.

Variable dependiente: Qp
Variable independiente: T

Aplicando el método:
⎡ 0 -152.51714286 ⎤
⎢ 1 13.82806071 ⎥ ⟨⟨1⟩⟩
R_P ((T , Qp , 2)) = ⎢ ⎥ a ≔ R_P ((T , Qp , 2))
⎢ 2 -0.07729643 ⎥
⎣ “ r” 0.99572413 ⎦

cal cal
a = -152.51714286 cal a = 13.82806071 ―― a = -0.07729643 ――
0 1 C° 2
C° 2

Qpa ((t)) ≔ a + a ⋅ t + a ⋅ t 2
0 1 2

APELLIDOS Y NOMBRES: Rey de Castro Gutierrez, Daniel Alejandro 3 Practica Calificada

 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
MÉTODOS NUMÉRICOS PARA INGENIERÍA QUÍMICA (PI-524/A)
SEMESTRE 2020 - I

b) Qp ＝ b + 2.31 ⋅ T + b ⋅ T 2 Variable dependiente: Qp

0 2
Variable independiente: T

Linealizando: Qp - 2.31 ⋅ T ＝ b + b ⋅ T 2
0 2

Se identifica un modelo de regresión lineal: x ≔ T2 y ≔ Qp - 2.31 ⋅ T

Aplicando el método:

⎡ “ao =” 54.38198122 ⎤ ⟨⟨1⟩⟩

R_L ((x , y)) = ⎢ “a1 =” 0.07483528 ⎥ c ≔ R_L ((x , y))
⎢ ⎥
⎣ “r =” 0.97278967 ⎦

Se tiene: Qp - 2.31 ⋅ T ＝ c + c ⋅ T 2
0 1

Qp ＝ c + 2.31 ⋅ T + c ⋅ T 2
0 1

Entonces:
cal
b ≔ c = 54.38198122 cal b ≔ c = 0.07483528 ――
0 0 2 1
C° 2

Qpb ((t)) ≔ b + 2.31 ⋅ t + b ⋅ t 2

0 2

Graficando ambos modelos: t ≔ 22 , 22.001 ‥ 52

385

360

335

310

285

260
Qp
235

210
Qpa ((t))
185
Qpb ((t))
160

135

110

22 25 28 31 34 37 40 43 46 49 52
T

APELLIDOS Y NOMBRES: Rey de Castro Gutierrez, Daniel Alejandro 3 Practica Calificada

 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
MÉTODOS NUMÉRICOS PARA INGENIERÍA QUÍMICA (PI-524/A)
SEMESTRE 2020 - I

Comentarios

Ÿ Se hallaron ambos modelos para representar el calor perdido por una masa
de 3 g en un cierto caso de transferencia de calor, el primero según una
regresión polinomial de segundo orden y el segundo mediante una ecuación
de segundo grado cuyo coeficiente de primer grado fue dado, ambas
respecto a la temperatura como variable independiente.

Ÿ Se observa que ambos modelos representan bien el conjunto de datos

obtenidos, siendo la gran mayoría de puntos cercanos a las gráficas de
dichos modelos. Además, esto se confirma mediante el coeficiente de
regresión, siendo de 0.995 para el primer caso y de 0.973 para el segundo.

Ÿ Si bien ambos son buenos modelos, se puede afirmar matemáticamente

que el primer modelo se ajusta mejor a los datos que el segundo, teniendo
un coeficiente de regresión muy cercano a 1 en comparación al del
segundo.

Pregunta 2
APELLIDOS Y NOMBRES: Rey de Castro Gutierrez, Daniel Alejandro 3 Practica Calificada
 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
MÉTODOS NUMÉRICOS PARA INGENIERÍA QUÍMICA (PI-524/A)
SEMESTRE 2020 - I

Pregunta 2

Se tienen los datos experimentales para el equilibrio líquido - vapor de una determinada
sustancia:

Ks T Se sabe: PT ≔ 0.95 atm

Pv
Además, Ks ＝ ――
5.2 300 PT

5.8 400
⎡ 4.94 ⎤
6.4 500 ⎢ 5.51 ⎥
⎢ ⎥
6.8 600 Pv ≔ Ks ⋅ PT = ⎢ 6.08 ⎥
7.5 700 ⎢ 6.46 ⎥
⎢⎣ 7.125 ⎥⎦

a) Entonces, para hallar los parámetros de la ecuación de Frost-Kalkward-Thodos (FKT) se

realizará una regresión empleando los datos de presión de vapor (Pv) y temperatura (T).

1 Pv
ln ((Pv)) ＝ A - B ⋅ ―+ C ⋅ ln ((T)) + D ⋅ ―― ((1))
T T2

Puesto que se tienen 4 parámetros, se tomará un dato para reducir el número de

ecuaciones:

Pv0 ≔ 6.08 T0 ≔ 500 Pv T

1 Pv0
ln ⎛⎝Pv0⎞⎠ ＝ A - B ⋅ ― + C ⋅ ln ⎛⎝T0⎞⎠ + D ⋅ ―― 4.94 300
T0 T0 2
5.51 400

((2)) 6.46 600

7.125 700

Restando (2) de (1):

⎛ 1 1⎞ ⎛ Pv Pv0 ⎞
ln ((Pv)) - ln ⎛⎝Pv0⎞⎠ ＝ B ⋅ ⎜― - ―⎟ + C ⋅ ⎛⎝ln ((T)) - ln ⎛⎝T0⎞⎠⎞⎠ + D ⋅ ⎜―― - ―― ⎟
⎝ T0 T ⎠
2
⎜⎝ T T0 2 ⎟⎠
⎛ Pv Pv0 ⎞
⎜―― 2
- ―― ⎟
ln ((Pv)) - ln ⎛⎝Pv0⎞⎠ ⎛⎝ln ((T)) - ln ⎛⎝T0⎞⎠⎞⎠ ⎜⎝ T T0 2 ⎟⎠
―――――― ＝ B + C ⋅ ―――――― + D ⋅ ――――
⎛ 1 1⎞ ⎛ 1 1⎞ ⎛ 1 1⎞
⎜― - ―⎟ ⎜― - ―⎟ ⎜― - ―⎟
⎝ T0 T ⎠ ⎝ T0 T ⎠ ⎝ T0 T ⎠

Entonces, se presenta un modelo de regresión multivariable: y ＝ a + a ⋅ x1 + a ⋅ x2

0 1 2

APELLIDOS Y NOMBRES: Rey de Castro Gutierrez, Daniel Alejandro 3 Practica Calificada

 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
MÉTODOS NUMÉRICOS PARA INGENIERÍA QUÍMICA (PI-524/A)
SEMESTRE 2020 - I

⎛ Pv Pv0 ⎞
ln ((Pv)) - ln ⎛⎝Pv0⎞⎠ ⎛⎝ln ((T)) - ln ⎛⎝T0⎞⎠⎞⎠
⎜―― - ―― ⎟
y ≔ ―――――― x1 ≔ ―――――― 2
T0 2 ⎟⎠
⎛ 1 1⎞ ⎛ 1 1⎞ ⎜⎝ T
⎜― - ―⎟ ⎜― - ―⎟ x2 ≔ ――――
⎛ 1 1⎞
⎝ T0 T ⎠ ⎝ T0 T ⎠ ⎜― - ―⎟
⎝ T0 T ⎠

Aplicando el método:

⟨ ⟩ ⟨ ⟩ ⟨ ⟩
D⟨0⟩ ≔ y D⟨1⟩ ≔ x1 D⟨2⟩ ≔ x2

⎡ “ ao = ” 1300.73923542 ⎤
⎢ “ a1 =” -0.5835627 ⎥ ⟨⟨1⟩⟩
RegMult ((D)) = ⎢ ⎥ a ≔ RegMult ((D))
⎢ “ a2 =” 40856.68602364 ⎥
⎣ “ r =” 0.92273933 ⎦

⎛ Pv Pv0 ⎞
⎜―― 2
- ―― ⎟
ln ((Pv)) - ln ⎛⎝Pv0⎞⎠ ⎛⎝ln ((T)) - ln ⎛⎝T0⎞⎠⎞⎠ ⎜⎝ T T0 2 ⎟⎠
―――――― ＝ a + a ⋅ ―――――― + a ⋅ ――――
⎛ 1 1⎞ 0 1 ⎛ 1 1⎞ 2 ⎛ 1 1⎞
⎜― - ―⎟ ⎜― - ―⎟ ⎜― - ―⎟
⎝ T0 T ⎠ ⎝ T0 T ⎠ ⎝ T0 T ⎠

⎛ 1 1⎞ ⎛ Pv Pv0 ⎞
ln ((Pv)) - ln ⎛⎝Pv0⎞⎠ ＝ a ⋅ ⎜― - ―⎟ + a ⋅ ⎛⎝ln ((T)) - ln ⎛⎝T0⎞⎠⎞⎠ + a ⋅ ⎜―― - ―― ⎟
⎝ T0 T ⎠
2
0 1 2
⎜⎝ T T0 2 ⎟⎠

1 Pv0 1 Pv
ln ((Pv)) ＝ ln ⎛⎝Pv0⎞⎠ + a ⋅ ― - a ⋅ ln ⎛⎝T0⎞⎠ - a ⋅ ―― - a ⋅ ―+ a ⋅ ln ((T)) + a ⋅ ――
0 T 2 0 T
0
1 2
T0 1 2
T2

1 Pv0
A ≔ ln ⎛⎝Pv0⎞⎠ + a ⋅ ― - a ⋅ ln ⎛⎝T0⎞⎠ - a ⋅ ―― = 7.03946207
0 T
0
1 2
T0 2
B ≔ a = 1300.73923542 °C
0

C ≔ a = -0.5835627
1

D ≔ a = 40856.68602364 (°C)^2 / atm

2

b) Se planteará el modelo dado por la siguiente ecuación:

B Pv T
ln ((Pv)) ＝ A - ――
T+C
ln ((Pv)) ⋅ ((T + C)) ＝ A ⋅ ((T + C)) - B 4.94 300
ln ((Pv)) ⋅ T + ln ((Pv)) ⋅ C ＝ A ⋅ T + A ⋅ C - B 5.51 400
ln ((Pv)) ⋅ T ＝ ((A ⋅ C - B)) + A ⋅ T + ((-C)) ⋅ ln ((Pv))
6.08 500
6.46 600
7.125 700

Se propone el modelo de regresión multivariable: y ＝ a + a ⋅ x1 + a ⋅ x2

APELLIDOS Y NOMBRES: Rey de Castro Gutierrez, Daniel Alejandro 0 1 2 3 Practica Calificada
 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
MÉTODOS NUMÉRICOS PARA INGENIERÍA QUÍMICA (PI-524/A)
SEMESTRE 2020 - I

Se propone el modelo de regresión multivariable: y ＝ a + a ⋅ x1 + a ⋅ x2

0 1 2
―――→
y ≔ ln ((Pv)) ⋅ T x1 ≔ T x2 ≔ ln ((Pv))

Aplicando el método:

⟨ ⟩ ⟨ ⟩ ⟨ ⟩
D⟨0⟩ ≔ y D⟨1⟩ ≔ x1 D⟨2⟩ ≔ x2

⎡ “ ao = ” 298.38605155 ⎤
⎢ “ a1 =” 2.5598809 ⎥ ⟨⟨1⟩⟩
RegMult ((D)) = ⎢ ⎥ a ≔ RegMult ((D))
⎢ “ a2 =” -372.93863684 ⎥
⎣ “ r =” 0.99926304 ⎦

ln ((Pv)) ⋅ T ＝ a + a ⋅ T + a ⋅ ln ((Pv))
0 1 2

ln ((Pv)) ⋅ T - a ⋅ ln ((Pv)) ＝ a + a ⋅ T
2 0 1

ln ((Pv)) ⋅ ⎛T - a ⎞ ＝ a + a ⋅ T
⎝ 2⎠ 0 1

a a
0 0
―+ T ―+ T - a + a
a +a ⋅T a a 2 2 a +a ⋅a
0 1 1 1 0 1 2
ln ((Pv)) ＝ ―――＝ a ⋅ ――― ＝ a ⋅ ―――――＝ a + ――――
T-a 1 T-a 1 T-a 1 T-a
2 2 2 2

A ≔ a = 2.5598809
1

B ≔ -⎛a + a ⋅ a ⎞ = 656.29244195 °C
⎝ 0 1 2⎠

C ≔ -a = 372.93863684 °C
2

Graficando con Desmos puesto que Mathcad no permite graficar ecuaciones implícitas:

Pv ((atm ))

T ((°C ))

Donde
APELLIDOS Y NOMBRES: Rey de Castro Gutierrez, Daniel Alejandro 3 Practica Calificada
 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
MÉTODOS NUMÉRICOS PARA INGENIERÍA QUÍMICA (PI-524/A)
SEMESTRE 2020 - I

Donde
1 Pv
---- ln ((Pv)) ＝ A - B ⋅ ―+ C ⋅ ln ((T)) + D ⋅ ――
T T2
B
---- ln ((Pv)) ＝ A - ――
T+C

Comentarios

Ÿ Se hallaron los parámetros para ambos modelos, tanto la ecuación de FKT como el
segundo modelo proporcionado. En ambos se obtuvo una buena correlación, dado el
alto coeficiente de regresión (>0.9 en ambos casos).

Ÿ En el caso del primer modelo, se tuvo que fijar un dato para poder utilizar el programa
de regresión multivariable sin problemas, de tal forma que se redujo el número de
variables independientes de 3 a 2. Dicho problema se hubiese evitado de tener un
programa que trabaje con tres variables independientes.

Ÿ Se observa que la ecuación de FKT presenta otra rama al ser graficada, implicando
que para mismos valores de temperatura existen diferentes valores de presión. Puesto
que esto no tiene sentido físico, se tomará únicamente la rama que se adecúe a los
datos experimentales obtenidos. Por cuestiones técnicas no se puede atenuar ni
eliminar la rama no utilizada.

Ÿ Además, se observa que la gráfica de la ecuación de FKT presenta una singularidad en

T = 500 K. Ello se debe a que para poder hallar los coeficientes con un programa de
dos variables, se tuvo que descartar el dato correspondiente a dicha temperatura.

Ÿ Se puede observar gráficamente que ambos modelos representan adecuadamente los

datos. Sin embargo, por el coeficiente de regresión se puede afirmar que el segundo
modelo es mejor que la ecuación de FKT, teniendo un coeficiente de regresión de
0.999 en comparación al 0.923.

Pregunta 3
APELLIDOS Y NOMBRES: Rey de Castro Gutierrez, Daniel Alejandro 3 Practica Calificada
 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
MÉTODOS NUMÉRICOS PARA INGENIERÍA QUÍMICA (PI-524/A)
SEMESTRE 2020 - I

Pregunta 3

De la ecuación matricial resuelta brindada, se deduce que los términos son:

A -1 ⋅ y ＝ a

Entonces, los coeficientes del polinomio de interpolación para los cuatro primeros puntos
de la tabla son:
La ecuación del polinomio de tercer grado:
⎡ -5 ⎤
⎢ -12 ⎥
a≔⎢ ⎥ y ≔ a + a ⋅ x + a ⋅ x2 + a ⋅ x3
⎢ -2 ⎥
0 1 2 3

⎣ 1 ⎦
y ≔ -5 - 12 ⋅ x - 2 ⋅ x 2 + x 3

Puesto que el valor de y para cada valor de x, para los cuatro primeros puntos, debe ser
exactamente el que se muestra en la tabla, reemplazando:

Para x ＝ -3 a y ＝ -14 ⎡ 1 ⎤
⎢ ⎥
2 3 solve ‾‾
13 - 5 ⎥
Se tiene: -14 ＝ -5 - 12 ⋅ ((-3 a)) - 2 ⋅ ((-3 a)) + ((-3 a)) ―― ⎢
→ ―――
⎢ 6 ⎥
⎢ ⎥
⎢ - ‾‾
13 - 5 ⎥
――――
Puesto que a es entero, a≔1 ⎢⎣ 6 ⎥⎦

-9 e
Para x＝a y ＝ ――
2
Sea e＝t para evitar problemas al calcular:
-9 ⋅ t 2 3
Se tiene: ――＝ -5 - 12 ⋅ ((a)) - 2 ⋅ ((a)) + ((a))
2
-2 ⎛
t ≔ ―― ⋅ ⎝-5 - 12 ⋅ ((a)) - 2 ⋅ ((a)) 2 + ((a)) 3 ⎞⎠ = 4 e≔4
9

5⋅b
Para x ＝ ―― y ＝ 10
2
⎡ 2 ⎤
2 ⎢ 3 ⎥
⎛5⋅b⎞ ⎛5⋅b⎞ ⎛ 5 ⋅ b ⎞ solve ⎢ 1i ⋅ ‾‾3 -3 ⎥
Se tiene: 10 ＝ -5 - 12 ⋅ ⎜―― ⎟ - 2 ⋅ ⎜―― ⎟ + ⎜―― ⎟ ―― → ――――
⎝ 2 ⎠ ⎝ 2 ⎠ ⎝ 2 ⎠ ⎢ 5 ⎥
⎢ ⎛ ⎞ ⎥
⎢ -⎝1i ⋅ ‾‾
3⎠ - 3 ⎥
―――――
Puesto que b es entero y real, b≔2 ⎢⎣ 5 ⎥⎦

Para x ＝ -b y＝3 d

Se tiene: 3 d ＝ -5 - 12 ⋅ ((-b)) - 2 ⋅ ((-b)) 2 + ((-b)) 3

1
d ≔ ―⋅ ⎛⎝-5 - 12 ⋅ ((-b)) - 2 ⋅ ((-b)) 2 + ((-b)) 3 ⎞⎠ = 1 d≔1
3

Para x＝2 b y＝4 c - 1

APELLIDOS Y NOMBRES: Rey de Castro Gutierrez, Daniel Alejandro 3 Practica Calificada
 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
MÉTODOS NUMÉRICOS PARA INGENIERÍA QUÍMICA (PI-524/A)
SEMESTRE 2020 - I

Para x＝2 b y＝4 c - 1

2 3
Se tiene: 4 c - 1 ＝ -5 - 12 ⋅ ((2 b)) - 2 ⋅ ((2 b)) + ((2 b))
1 ⎛ 2 3⎞
c ≔ ―⋅ ⎝-5 + 1 - 12 ⋅ ((2 b)) - 2 ⋅ ((2 b)) + ((2 b)) ⎠ = -5 c ≔ -5
4

Entonces, finalmente W≔a+b+c+d+e

W=3

Comentarios

Ÿ Contrastando con la ecuación matricial para la interpolación por

aproximación polinomial simple de tercer grado, se deduce que los valores
brindados son los correspondientes a los coeficientes del polinomio de
interpolación.

Ÿ Puesto que una interpolación cumple con que se conservan todos los
valores utilizados, en este caso los cuatro primeros, se puede reemplazar en
la ecuación igualando los datos mostrados en la tabla.

Ÿ Además, se utiliza la condición dada de que el valor de y en 5b/2 es igual a

10, el cual, coincidentemente, es igual al brindado en la tabla. No obstante,
esto es una mera coincidencia y no es implicado por la teoría, puesto que
para la interpolación dicho valor no se tomó en cuenta.

Ÿ Reemplazando para cada valor de la tabla, se calcularon los valores de a, b,

c, d y e, para al final dar la respuesta W = a + b+ c + d + e deseada de 3.

Pregunta 4
APELLIDOS Y NOMBRES: Rey de Castro Gutierrez, Daniel Alejandro 3 Practica Calificada
 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
MÉTODOS NUMÉRICOS PARA INGENIERÍA QUÍMICA (PI-524/A)
SEMESTRE 2020 - I

Pregunta 4

Se tienen los siguientes datos de concentración de la sustancia a diferentes condiciones de

presión y temperatura:

P T150 T170 T190 T200

1 0.078 0.006 0.005 0.005

3 0.143 0.053 0.047 0.045
5 0.167 0.147 0.132 0.125
6 0.156 0.138 0.129 0.117
8 0.143 0.119 0.112 0.108

a) Primero se hallará, por interpolación de Lagrange con tres puntos, la concentración a

diferentes temperaturas, considerando una presión de 3.8 kg/cm2.

Los 3 valores más cercanos: Pint ≔ 3.8 ⎡ 2.8 ⎤

⎢ 0.8 ⎥ 3
――→ ⎢
||P - Pint|| = 1.2 ⎥ 5
⎢ ⎥
⎢ 2.2 ⎥ 6
⎢⎣ 4.2 ⎥⎦

P T150 T170 T190 T200 Para Pint ≔ 3.8 siendo CT

la concentración a una temperatura T para
una presión de 3.8 kg/cm2:
3 0.143 0.053 0.047 0.045
5 0.167 0.147 0.132 0.125 C150 ≔ L ⎛⎝Pint , P , T150⎞⎠ = 0.15996

6 0.156 0.138 0.129 0.117 C170 ≔ L ⎛⎝Pint , P , T170⎞⎠ = 0.10852

C190 ≔ L ⎛⎝Pint , P , T190⎞⎠ = 0.09556

C200 ≔ L ⎛⎝Pint , P , T200⎞⎠ = 0.09236

Luego, para la presión fija de 3.8 kg/cm2 se tiene:

T C Se desea hallar la concentración para Tint ≔ 172 K

Los tres valores más cercanos:

150 C150 T C
⎡ 22 ⎤
150
170 C170 ――→ ⎢ 2 ⎥
||T - Tint|| = ⎢ ⎥ 170
190 C190 ⎢ 18 ⎥ 190 150 C150
⎣ 28 ⎦
200 C200 170 C170
190 C190

Entonces, aplicando Lagrange con tres puntos:

APELLIDOS Y NOMBRES: Rey de Castro Gutierrez, Daniel Alejandro 3 Practica Calificada
 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
MÉTODOS NUMÉRICOS PARA INGENIERÍA QUÍMICA (PI-524/A)
SEMESTRE 2020 - I

Entonces, aplicando Lagrange con tres puntos:

C ≔ L ⎛⎝Tint , T , C⎞⎠ = 0.1054924 mol/L

b) Primero se hallará, por interpolación de Lagrange con tres puntos, la concentración a

diferentes temperaturas, considerando una presión de 4 kg/cm2.

P T150 T170 T190 T200

1 0.078 0.006 0.005 0.005

3 0.143 0.053 0.047 0.045
5 0.167 0.147 0.132 0.125
6 0.156 0.138 0.129 0.117
8 0.143 0.119 0.112 0.108

Los 3 valores más cercanos: Pint ≔ 4 ⎡3⎤

⎢1⎥ 3
――→ ⎢ ⎥
||P - Pint|| = 1 5
⎢ ⎥
⎢2⎥ 6
⎢⎣ 4 ⎥⎦

P T150 T170 T190 T200 Para Pint ≔ 4 siendo CT

la concentración a una temperatura T para
una presión de 3.8 kg/cm2:
3 0.143 0.053 0.047 0.045
5 0.167 0.147 0.132 0.125 C150 ≔ L ⎛⎝Pint , P , T150⎞⎠ = 0.16266667

6 0.156 0.138 0.129 0.117 C170 ≔ L ⎛⎝Pint , P , T170⎞⎠ = 0.11866667

C190 ≔ L ⎛⎝Pint , P , T190⎞⎠ = 0.10466667

C200 ≔ L ⎛⎝Pint , P , T200⎞⎠ = 0.101

Luego, para la presión fija de 4 kg/cm2 se tiene:

mol
C T Se desea hallar la temperatura para Cint ≔ 0.117 ――
Los tres valores más cercanos: L

C150 150
⎡ 0.04566667 ⎤
C170 170 ―――→ ⎢ 0.00166667 ⎥ C T
||C - Cint|| = ⎢ ⎥ 0.1186666
C190 190 ⎢ 0.01233333 ⎥ 0.1046666
⎣ 0.016 ⎦ C170 170
0.101
C200 200
C190 190
C200 200

Entonces, aplicando Lagrange con tres puntos:

APELLIDOS Y NOMBRES: Rey de Castro Gutierrez, Daniel Alejandro 3 Practica Calificada
 UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
MÉTODOS NUMÉRICOS PARA INGENIERÍA QUÍMICA (PI-524/A)
SEMESTRE 2020 - I

Entonces, aplicando Lagrange con tres puntos:

T ≔ L ⎛⎝Cint , C , T⎞⎠ = 170.86988483 K

Comentarios

Ÿ Haciendo una doble interpolación en cada caso, se pudo hallar los valores
de concentración y temperatura solicitados empleando el cuadro de doble
entrada.

Ÿ Se utilizó el método de interpolación de Lagrange trabajando con 3 puntos,

por lo que cada aplicación del método requirió determinar cuáles serían
estos tres puntos más cercanos. Ello se determinó calculando el valor
absoluto de la diferencia entre cada elemento con el valor de interpolación
deseado para P, T y C.

Ÿ En el caso del item b), para hallar la temperatura en función de la

concentración se tuvo que realizar el inverso de lo hecho en el item a) tras
calcular por interpolación la concentración a distintas temperaturas con una
presión fija: en lugar de tomar como variable dependiente y como valor de
interpolación a la temperatura, se tomó la concentración solicitada.

Ÿ En todos los casos se verifica que los datos obtenidos por interpolación
sigan la tendencia de los datos en el cuadro original, así como que los
valores sean relativamente cercanos a aquellos de los valores dados que
sean cercanos al valor de interpolación.

APELLIDOS Y NOMBRES: Rey de Castro Gutierrez, Daniel Alejandro 3 Practica Calificada