INSTITUTO POLITÉCNICO NACIONAL

ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E NDUSTRIAS EXTRACTIVAS

SEGUNDO EXAMEN PARCIAL SIMULACIÓN DE PROCESOS
GARCÍA PINEDA RONALDO – OLEA DÍAZ GEORGINA LIZET MV92

Problema 1
Realizar el análisis de la transferencia de calor en estado estacionario utilizando la ecuación
de Laplace, en la sección transversal de una placa rectangular de 25x25 cm. Considerando
4 y 9 nodos.
Área: (0.25m) (0.25m) = 0.0625 m2
Propiedades termo-físicas del Aluminio 6061
ρ: 2680 Kg/m3 k: 229 W/m K Cp: 1007.46 J/Kg K
Ecuación de Laplace Ecuación de transferencia de calor
∂2 T ∂2 T (Conducción) por diferencias finitas
+ =0 derivada de la ecuación de Laplace.
∂x 2 ∂y 2
a) 4 nodos: Ti+1,j + Ti-1,j + Ti,j+1 + Ti,j-1 = 0

Nodo T11 Nodo T22

T21 + 620 + T12 + 560 – 80 + T12 + 100 + T21 –
4T11 = 0 4T22 = 0
-4T12 + T21 + T12 = -1180 -4T22 + T12 + T21 = -180
Nodo T12 Nodo T21
T22 + 620 + 100+ T11 – 4T12 = 0 80 + T11 + T22 + 560 – 4T21 = 0
-4T22 + T22 + T12 = -720 4T21 + T11 + T22 = -640
Matriz
-4T11 + T21 + T12 + 0T22 = -1180
T11 + 0T21 - 4T12 + T22 = -720
0T11 + T21 + T12 - 4T22 = -640
T11 - 4T21 + 0T12 + T22 = -640
T11=465 °C T12=350 °C
T21=330 °C T22=215 °C
 INSTITUTO POLITÉCNICO NACIONAL
ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E NDUSTRIAS EXTRACTIVAS
SEGUNDO EXAMEN PARCIAL SIMULACIÓN DE PROCESOS
GARCÍA PINEDA RONALDO – OLEA DÍAZ GEORGINA LIZET MV92

b) 9 Nodos

Nodo T11 Nodo T12 Nodo T13

T21 + 620 + T12 + 560 - 4T11 = 0 T22 + 620 + T13 + T11 - 4T12 = 0 T23 + 620 + 100 + T12 - 4T13 = 0

-4T11 + T21 + T12 = -1120 T11 - 4T12 + T13 + T22 = -620 T12 - 4T13 + T23 = -720

Nodo T21 Nodo T22 Nodo T23

T31 + T11 + T22 + 560 -4T21 = 0 T32 + T12 + T23 + T21 - 4T22 = 0 T33 + T13 + 100 + T22 - 4T23 = 0

T11 -4T21 + T21 + T22 = -560 T12 + T21 - 4T22 + T23 + T32 = 0 T13 + T22 - 4T23 + T33 = -100

Nodo T31 Nodo T32 Nodo T33

80 + T21 + T32 + 560 -4T31 = 0 80 + T22 + T33 + T31 - 4T32 = 0 80 + T23 + 100 + T32 - 4T33 = 0

T21 -4T31 + T32 = -640 T22 + T31 - 4T32 + T33 = -80 T23 + T32 - 4T33 = -180

Matriz

-4T11 + T21 + 0T31 + T12 + 0T22 + 0T32 + 0T13 + 0T23 + 0T33 = -1180

T11 - 4T21 + T31 + 0T12 + T22 + 0T32 + 0T13 + 0T23 + 0T33 = -560

0T11 + T21 - 4T31 + 0T12 + 0T22 + T32 + 0T13 + 0T23 + 0T33 = -640

T11 + 0T21 + 0T31 - 4T12 + 0T22 + 0T32 + T13 + 0T23 + 0T33 = -620

0T11 + T21 + 0T31 + T12 - 4T22 + T32 + 0T13 + T23 + 0T33 = 0

0T11 + 0T21 + T31 + 0T12 + T22 - 4T32 + 0T13 + 0T23 + T33 = -80

0T11 + 0T21 + 0T31 + T12 + 0T22 + 0T32 - 4T13 + T23 + 0T33 = -720

0T11 + 0T21 + 0T31 + 0T12 + T22 + 0T32 + T13 - 4T23 + T33 = -100

0T11 + 0T21 + 0T31 + 0T12 + 0T22 + T32 + 0T13 + T23 - 4T33 = -180

T11= 488.88 °C T12= 358.37 °C T13= 324.6 °C

T21= 417.18 °C T22= 300.22 °C T23= 220.03 °C
T31= 317.61 °C T32= 213.29 °C T33= 153.33 °C
 INSTITUTO POLITÉCNICO NACIONAL
ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E NDUSTRIAS EXTRACTIVAS
SEGUNDO EXAMEN PARCIAL SIMULACIÓN DE PROCESOS
GARCÍA PINEDA RONALDO – OLEA DÍAZ GEORGINA LIZET MV92

Comparativa.
Discretizar el espacio en un mayor número de nodos permite obtener temperaturas más
aproximadas a los valores reales, además de dar una idea más precisa de cómo fluye el calor
en la placa. Por lo que se recomienda discretizar el espacio en un mayor número de nodos
sin caer en excesos que pudieran ralentizar el software.
Problema 2
Realizar el modelado de una placa inclinada y abierta a la superficie, donde el espesor de la
capa de 2 cm de espesor, 0.4 cm de longitud y presenta una inclinación de 60º respecto a
la vertical. Determinar los perfiles de velocidad.
p=1050 kg/m3
µ= 22.5 mPa*s
T= 25 °C
𝜌𝑔𝑐𝑜𝑠𝛽 𝑥 2 𝑉2 (2𝜇) 1 𝑥 2 2𝜇
V2=( 𝛿 2 )(1 − (𝛿) ) → ∗ = (1 − (𝛿) ) 𝑥𝛿 2 → 𝑉2 ∗ 𝜌𝑔𝑐𝑜𝑠𝛽 ∗ (𝛿 2 − 𝑥 2 )
2𝜇 𝜌𝑔𝑐𝑜𝑠𝛽 𝛿2

𝐾𝑔
2𝜇 2(0.0225 2 ∗ 𝑠)
𝐾𝑛 = = 𝑚𝑠 = 8.7374𝑥10−6 𝑚 ∗ 𝑠
𝜌𝑔𝑐𝑜𝑠𝛽 (1050 𝐾𝑔)(9.8 𝑚 )(cos 60°)
𝑚3 𝑠2
Matrices de rigidez.

Elemento 1 (Nodo 1 y 2)

1 K1 (V21 – V22) = (δ2) 𝐾1 −𝐾1 𝑉1 0.0004

[ ] ∗ [ 22] = [ ]
−𝐾1 𝐾1 𝑉3 0.000393
2 -K1 (V21 – V22) = (δ2 – Xi2)

Elemento 2 (Nodo 2 y 3)

2 K2 (V22 – V23) = (δ2-X12) 𝐾2 −𝐾2 𝑉2 0.000393

[ ] ∗ [ 2 3] = [ ]
−𝐾2 𝐾2 𝑉3 0.000375
3 -K2 (V22 – V23) = (δ2 – X22)

Elemento 3 (Nodo 3 y 4)

3 K3 (V23 – V24) = (δ2-X32) 𝐾3 −𝐾3 𝑉3 0.000375

[ ] ∗ [ 2 4] = [ ]
−𝐾3 𝐾3 𝑉3 0.000343
4 -K3 (V23 – V24) = (δ2 – X32)

Elemento 4 (Nodo 4 y 5)

4 K4 (V24 – V25) = (δ2-X32) 𝐾4 −𝐾4 𝑉4 0.000343

[ ] ∗ [ 2 5] = [ ]
−𝐾4 𝐾4 𝑉3 0.000300
5 -K4 (V24 – V25) = (δ2 – X42)

Elemento 5 (Nodo 5 y 6)
 INSTITUTO POLITÉCNICO NACIONAL
ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E NDUSTRIAS EXTRACTIVAS
SEGUNDO EXAMEN PARCIAL SIMULACIÓN DE PROCESOS
GARCÍA PINEDA RONALDO – OLEA DÍAZ GEORGINA LIZET MV92

5 K5 (V25 – V26) = (δ2-X42) 𝐾5 −𝐾5 𝑉5 0.000300

[ ] ∗ [ 2 6] = [ ]
−𝐾5 𝐾5 𝑉3 0.000243
6 -K5 (V25 – V26) = (δ2 – X42)

Elemento 6 (Nodo 6 y 7)

6 K6 (V26 – V27) = (δ2-X52) 𝐾6 −𝐾6 𝑉5 0.000243

[ ] ∗ [ 2 6] = [ ]
−𝐾6 𝐾6 𝑉3 0.000175
7 -K6 (V26 – V27) = (δ2 – X62)

Elemento 7 (Nodo 7 y 8)

7 K7 (V27 – V28) = (δ2-X62) 𝐾7 −𝐾7 𝑉6 0.000175

[ ] ∗ [ 2 7] = [ ]
−𝐾7 𝐾7 𝑉3 0.000937
8 -K7 (V27 – V28) = (δ2 – X72)

Elemento 8 (Nodo 8 y 9)

8 K8 (V28 – V29) = (δ2-X72) 𝐾 −𝐾8 𝑉8 0.0000937

[ 8 ] ∗ [ 29] = [ ]
−𝐾8 𝐾8 𝑉3 0
9 -K8 (V28 – V29) = (δ2 – X82)

Matriz acoplada de rigidez.

K1 -K1 0 0 0 0 0 0 0 0.0004

-K1 K1+k2 -K2 0 0 0 0 0 0 0.000343

0 -K2 K2+K3 -K3 0 0 0 0 0 0.000375

0 0 +K3 K3+K4 -K4 0 0 0 0 0.000343

0 0 0 -K4 K4+K5 -K5 0 0 0 0.0003000

0 0 0 0 -K5 K5+K6 -K6 0 0 0.000243

0 0 0 0 0 -K6 K6+K7 -K7 0 0.000175

0 0 0 0 0 0 -K8 K7+K8 -K8 0.0000937

0 0 0 0 0 0 0 0 K8 0
 INSTITUTO POLITÉCNICO NACIONAL
ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E NDUSTRIAS EXTRACTIVAS
SEGUNDO EXAMEN PARCIAL SIMULACIÓN DE PROCESOS
GARCÍA PINEDA RONALDO – OLEA DÍAZ GEORGINA LIZET MV92

Por cuestiones de practicidad la Kn se consideró como entero. Sin embargo, tenía valores de 10-6,
por lo que las velocidades obtenidas se multiplicaron por este mismo valor, para obtener el valor
real.

V21= V22= V23= V24= V25= V26= V27= V28= V29 = 1.3075x1012*(10-6) = 1.30756x106 m/s

Sin embargo, los valores obtenidos no tienen lógica puesto que se están alcanzando velocidades
muy altas, por lo que se debe calcular el valor del Reynold medio para determinar si este flujo sí es
laminar.

Calcular Vz media para calcular Re

2 Kg m
δ2 ρgcosβ (0.02m) (1030 m3 )(9.8 s 2 )cos(60°) m
Vz = = 2 = 30.52
3μ Kgm s
3(0.0225 s )

Calculando el Reynolds medio

m Kg
4δV2 ρ 4(0.02m)(30 s )(1050 m3 )
Re = = = 114545.45
μ Kgm2
0.022 s

El Re es mayor al valor del de un flujo laminar, así que este modelo se debe analizar como flujo
turbulento. Por lo que se concluye que los perfiles de velocidad calculados no son correctos. Y se
debe simular en algún software especializado.