Facultad de Ingeniería Química

Escuela Profesional De Ingeniería Química

CÁTEDRA : MODELAMIENTO Y SIMULACIÓN DE PROCESOS
CATEDRÁTICO : ING. PASCUAL VISTOR GUEVARA YANQUI

SOLUCIONARIO DEL PRIMER PARCIAL

ESTUDIANTE : MIGUEL ANGEL VASQUEZ CABALLERO

CODIGO : 2016200772B

Huancayo – PERU
2020
1. Resolver la ecuación: 𝑥𝑑𝑦 − 𝑛𝑑𝑦 = 0 y realizar un análisis e interpretación
cuando 𝑛 = 3.

Solución:
𝑥𝑑𝑦 − 𝑛𝑑𝑦 = 0
𝑑𝑦 𝑑𝑥
=3
𝑦 𝑥
Integrando: 𝑥 ≠ 0 & 𝑦 ≠ 0
𝑑𝑦 𝑑𝑥
∫ = 3∫
𝑦 𝑥
ln 𝑦 = 3ln 𝑥 + 𝑐
También se puede poner:
ln 𝑦 = 3ln 𝑥 + ln 𝑐
Aplicando el teorema de logaritmos:
𝑙𝑛𝑥 𝑛 = 𝑛𝑙𝑛𝑥
𝑙𝑛(𝑎. 𝑏) = 𝑙𝑛(𝑎) + 𝑙𝑛(𝑏)
ln(𝑦) = ln(𝑥 3 ) + ln 𝑐
→ 𝑦 ≡ 𝑐 𝑥3
 Análisis: cuando 𝑐 > 0 y 𝑥 > 0, el valor de 𝑦 toma valores
positivos.

𝑎
2. La integral ∫0 𝑥 𝑛−1 𝑒 −𝑥 𝑑𝑥 = Γ(𝑛) se denomina función Gamma de n. evaluar
la integral siguiente:
𝑎
3
∫ 𝑒 −𝑥 𝑑𝑥
0

Solución:
𝑎

∫ 𝑥 𝑛−1 𝑒 −𝑥 𝑑𝑥 = Γ(𝑛) …. (1)

0
Resolviendo la integral:
𝑎
3
∫ 𝑒 −𝑥 𝑑𝑥 … . (2)
0

𝑢 1⁄ 3 = 𝑥
∴ → 𝑥 = 0, 𝑢 = 0 & 𝑥 = 𝑎, 𝑢 = 𝑎3
𝑑𝑥
= 𝑢 1⁄ 3
𝑑𝑢
𝑑𝑥 1 −2
= 𝑢 3
𝑑𝑢 3
1 2
𝑑𝑥 = 𝑢−3 𝑑𝑢
3
Reemplazando en la ecuación(2).
𝑎
3
∫ 𝑒 −𝑥 𝑑𝑥
0
𝑎
1 2
∫ 𝑒 −𝑢 𝑢−3 𝑑𝑢 … . . (3)
3
0

Igualando la ecuación (1) con la ecuación (3):

𝑎 𝑎
1 2
∫ 𝑥 𝑛−1 𝑒 −𝑥 𝑑𝑥 = ∫ 𝑢−3 𝑒 −𝑢 𝑑𝑢
3
0 0

2
→ 𝑛−1= −
3
1
∴ 𝑛=
3
1 1
…….= Γ( )
3 3
3. Escribir las series desarrolladas para la función Bessel, cuando n=0 y cuando
n=3 y demostrar que:
4𝐽0′′′ (𝑥 ) + 3𝐽0′ (𝑥 ) + 𝐽3 (𝑥 ) = 0

Solución:
La función Bessel para números enteros positivos:
𝑛
𝑥 𝑛+2𝑘
𝐽𝑛 (𝑥 ) = ∑(−1)𝑘
2𝑛+2𝑘 . 𝑘! . (𝑛 + 𝑘 )!
𝑘=0

Para n=0
𝑛
𝑥 0+2𝑘
𝐽0 (𝑥 ) = ∑(−1)𝑘
20+2𝑘 . 𝑘! . (0 + 𝑘 )!
𝑘=0
𝑛
(−1)𝑘 𝑥 2𝑘
𝐽0 (𝑥 ) = ∑
(𝑘!)2 22𝑘
𝑘=0

𝑥2 𝑥4 𝑥6 𝑥8 𝑥 10
𝐽0 (𝑥 ) = 1 − + − + −
22 . 1! .1! (2!)2 . 24 (3!)2 . 26 (4!)2 . 28 (5!)2 . 210
2𝑥 4𝑥 3 6𝑥 5 8𝑥 7 10𝑥 9
𝐽0′ (𝑥 ) = − + − + −
22 . 1! .1! (2!)2 . 24 (3!)2 . 26 (4!)2 . 28 (5!)2 . 210
2 12𝑥 2 30𝑥 4 56𝑥 6 90𝑥 8
𝐽0′′ (𝑥 ) = − + − + −
22 . 1! .1! (2!)2 . 24 (3!)2 . 26 (4!)2 . 28 (5!)2 . 210
24𝑥 120𝑥 3 336𝑥 5 720𝑥 7
𝐽0′′′ (𝑥 ) = − + −
(2!)2 . 24 (3!)2 . 26 (4!)2 . 28 (5!)2 . 210
Para n=3
𝛼
𝑥 3+2𝑘
𝐽3 (𝑥 ) = ∑(−1)𝑘
23+2𝑘 . 𝑘! . (3 + 𝑘 )!
𝑘=0

𝑥3 𝑥5 𝑥7
𝐽3 (𝑥 ) = − −
3!. 0!. 23 4!. 1!. 25 5!. 2!. 27

∴
−2𝑥 4𝑥 3 6𝑥 5 8𝑥 7 10𝑥 9
3𝐽0′ (𝑥 ) = 3. ( 2 + − + − ) … (1)
2 . 1! .1! (2!)2 . 24 (3!)2 . 26 (4!)2 . 28 (5!)2 . 210
 24𝑥 120𝑥 3 336𝑥 5 720𝑥 7
4𝐽0′′′ (𝑥 ) = 4. ( − + − ) … . (2)
(2!)2 . 24 (3!)2 . 26 (4!)2 . 28 (5!)2 . 210
𝑥3 𝑥5 𝑥7
𝐽3 (𝑥 ) = − − … . . (3)
3!. 0!. 23 4!. 1!. 25 5!. 2!. 27
 Sumamos todos los valores de 𝑥 de las ecuaciones 1, 2 y 3.

6𝑥 (4)24𝑥
− + =0
22 . 1! .1! (2!)2 . 24

 Sumamos todos los valores que tienen 𝑥 3 de las

ecuaciones 1, 2 y 3.

12𝑥 3 (4)120𝑥 3 𝑥3
− + =0
(2!)2 . 24 (3!)2 . 26 3!. 0!. 23
 Sumamos todos los valores que tienen 𝑥 5 de las
ecuaciones 1, 2 y 3.
(3)6𝑥 5 (4)336𝑥 5 𝑥5
− + − =0
(3!)2 . 26 (4!)2 . 28 4!. 1!. 25
 Comparando sucesivamente los otros valores, se va
haciendo cero; por lo tanto, se determina que:

4𝐽0′′′ (𝑥 ) + 3𝐽0′ (𝑥 ) + 𝐽3 (𝑥 ) = 0
4. Se desea analizar una columna de destilación desde el punto de vista de lo
balances de materia, lo cual se cuenta con el siguiente sistema físico:

Solución:

Aplicando el balance de materia globales en los 3 puntos y las relaciones en

el enunciado. Se sabe que:
𝑙𝑏𝑚𝑜𝑙
𝐴1 = 400
ℎ𝑟
𝐴2 = 1.6𝑥
𝐴3 = 𝑥
𝐴4 = 0.6𝑥
𝐴5 = 1.5𝑦
𝐴6 = 𝑦
𝐴7 = 0.5𝑦
RESOLVIENDO x e y
BALANCE TOTAL
𝐴1 = 𝐴3 + 𝐴6
400 = 𝑥 + 𝑦 (1)
BALANCE PARCIAL
𝑥𝐴1 𝐴1 = 𝑥𝐴3 𝐴3 + 𝑥𝐴6 𝐴6
 400 ∗ 0.35 = 𝑥 (0.75) + 𝑦(0.05)
140 = 0.75𝑥 + 0.05𝑦 (2)

RESOLVIENDO SEL:
0.75 ∗ (𝑥 = 400 − 𝑦) (-)
0.75𝑥 = 140 − 0.05𝑦
0.75𝑥 = 300 − 0.75𝑦
0.75𝑥 = 140 − 0.05𝑦
0 = 160 − 0.70𝑦
160
=𝑦
0.70
𝑙𝑏𝑚𝑜𝑙
𝑦 = 228.5714
ℎ𝑟
REEMPLAZANDO EN (1):
𝑙𝑏𝑚𝑜𝑙
𝑥 = 171.43
ℎ𝑟
ENTONCES:
CORRIENTE FLUJO 𝑥𝐴 𝑥𝐵 𝑥𝐶 FASE
𝐴1 400 0.35 0.35 0.30 VAPOR
𝐴2 274.29 0.75 𝑥𝐵2 𝑥𝐶2 VAPOR
𝐴3 171.43 0.75 𝑥𝐵3 𝑥𝐶3 LIQUIDO
𝐴4 102.86 0.75 𝑥𝐵4 𝑥𝐶4 LIQUIDO
𝐴5 342.86 0.05 𝑥𝐵5 𝑥𝐶5 LIQUIDO
𝐴6 228.57 0.05 𝑥𝐵6 𝑥𝐶6 LIQUIDO
𝐴7 114.29 0.05 𝑥𝐵7 𝑥𝐶7 VAPOR
 APLICAMOS BALANCE PARCIALES CON RESPECTO A 𝑥𝐵 y 𝑥𝐶 DE
TODO EL DESTILADOR:
CONSIDERANDO FLUJOS HOMOGENEOS LAS COMPOSICIONES
SERAN:
𝑥𝐵𝑇 = 𝑥𝐵2 = 𝑥𝐵3 = 𝑥𝐵4
𝑥𝐵𝐹 = 𝑥𝐵5 = 𝑥𝐵6 = 𝑥𝐵7

𝑥𝐶𝑇 = 𝑥𝐶2 = 𝑥𝐶3 = 𝑥𝐶4

𝑥𝐶𝐹 = 𝑥𝐶5 = 𝑥𝐶6 = 𝑥𝐶7

ENTONCES:
0.35 ∗ 400 = 171.43𝑥𝐵𝑇 + 228.57𝑥𝐵𝐹
0.30 ∗ 400 = 171.43𝑥𝐶𝑇 + 228.57𝑥𝐶𝐹
𝑥𝐵𝑇 + 𝑥𝑐𝑇 = 1 − 0.75
𝑥𝐵𝐹 + 𝑥𝑐𝐹 = 1 − 0.05

FORMANDO EL SEL:

𝑥𝐵𝑇 𝑥𝑐𝑇 𝑥𝐵𝐹 𝑥𝑐𝐹

171.43 0 228.57 0 140
0 171.43 0 228.57 120
1 1 0 0 0.25
0 0 1 1 0.95
RESOLVIENDO Y REDONDEANDO A CENTECIMALES:
217
𝑥𝐵𝑇 = = 0.2523 ≅ 0.25
860
1
𝑥𝑐𝑇 = − = −0.0023 ≅ 0
430
2773323
𝑥𝐵𝐹 = = 0.4232 ≅ 0.42
6552340
172570
𝑥𝑐𝐹 = = 0.5267 ≅ 0.53
327617
 ENTONCES COMPLETANDO EL CUADRO QUEDARIA:
CORRIENTE FLUJO 𝑥𝐴 𝑥𝐵 𝑥𝐶 FASE
𝐴1 400 0.35 0.35 0.30 VAPOR
𝐴2 274.29 0.75 0.25 0 VAPOR
𝐴3 171.43 0.75 0.25 0 LIQUIDO
𝐴4 102.86 0.75 0.25 0 LIQUIDO
𝐴5 342.86 0.05 0.42 0.53 LIQUIDO
𝐴6 228.57 0.05 0.42 0.53 LIQUIDO
𝐴7 114.29 0.05 0.42 0.53 VAPOR

5. Un fluido incompresible es bombeado a un tanque cónico invertido. El fluido

sale por un orificio por la parte inferior. Derivar una ecuación que describa el
sistema de vaciado del líquido del depósito y determinar el tiempo necesario
para vaciar el líquido. Analizar e interpretar y discutir los siguientes casos.
a) Flujo laminar
b) Flujo de transición
c) Flujo turbulento
d) Si la salida del orificio se aumenta un tubo de salida de longitud “L”

Solución:
 BALANCE DE MATERIA
Condiciones.
o Generación de masa= 0
o Consumo de masa= 0

Aplicando la expresión integral de la conservación de la masa.

[𝑒𝑛𝑡𝑟𝑎𝑑𝑎 𝑑𝑒 𝑚𝑎𝑠𝑎] = [𝑎𝑐𝑢𝑚𝑢𝑙𝑎𝑐𝑖𝑜𝑛] − [𝑔𝑒𝑛𝑒𝑟𝑎𝑐𝑖𝑜𝑛] + [𝑠𝑎𝑙𝑖𝑑𝑎]
𝜕
∬ 𝜌(𝑣.
̅ 𝑛̅)𝑑𝐴 + ∭ 𝜌𝑑𝑉 = 0
𝑠𝑐 𝜕𝑡 𝑣𝑐
𝜕
∭ 𝜌𝑑𝑉 = 0 , 𝑛𝑜 ℎ𝑎𝑦 𝑎𝑐𝑢𝑚𝑢𝑙𝑎𝑐𝑖𝑜𝑛
𝜕𝑡 𝑣𝑐
Integrando.
𝜌1 𝐴1 𝑉1 = 𝜌2 𝐴2 𝑉2
Fluido incompresible.
𝐴1 𝑉1 = 𝐴2 𝑉2
Despejando.
𝐴2 𝑉2
𝑉1 =
𝐴1

BALANCE DE ENERGÍA
𝛿𝑄 𝛿𝑊𝑠 𝑃 𝜕
− = ∬ (𝑒 + ) 𝜌(𝑣.
̅ 𝑛̅)𝑑𝐴 + ∭ 𝑒𝜌𝑑𝑉
𝛿𝑡 𝑑𝑡 𝑠𝑐 𝜌 𝜕𝑡 𝑣𝑐

𝜕 𝛿𝑄 𝛿𝑊𝑠
∭ 𝑒𝜌𝑑𝑉 = 0 ∧ = . . . . (𝛼 )
𝜕𝑡 𝑣𝑐 𝛿𝑡 𝑑𝑡
𝛿𝑄 𝛿𝑊𝑠 𝑃
− = ∬ (𝑒 + ) 𝜌(𝑣.
̅ 𝑛̅)𝑑𝐴
𝛿𝑡 𝑑𝑡 𝑠𝑐 𝜌
𝑃 𝑣2 𝑃
(𝑒 + ) = 𝑔𝑦 − −𝑢−
𝜌 2 𝜌
 𝑃 𝑣2 2 𝑃2 𝑣1 2 𝑃1
∬ (𝑒 + ) 𝜌(𝑣.
̅ 𝑛̅)𝑑𝐴 = [𝑔𝑦2 + + 𝑢 + ] (𝜌2 𝐴2 𝑉2 ) − [𝑔𝑦1 + + 𝑢 + ] (𝜌1 𝐴1 𝑉1 )
𝑠𝑐 𝜌 2 𝜌 2 𝜌
𝛿𝑄 𝛿𝑊𝑠 𝑣2 2 𝑃2 𝑣1 2 𝑃1
− = [𝑔𝑦2 + + 𝑢 + ] (𝜌2 𝐴2 𝑉2 ) − [𝑔𝑦1 + + 𝑢 + ] (𝜌1 𝐴1 𝑉1 )
𝛿𝑡 𝑑𝑡 2 𝜌 2 𝜌

Asumiendo las siguientes condiciones.

No existe transferencia de calor.
𝛿𝑄
=0
𝛿𝑡
El sistema no realiza trabajo.
𝛿𝑊𝑠
=0
𝑑𝑡

𝑣2 2 𝑃2 𝑣1 2 𝑃1
0 = [𝑔𝑦2 + +𝑢+ ] (𝜌2 𝐴2 𝑉2 ) − [𝑔𝑦1 + +𝑢+ ] (𝜌1 𝐴1 𝑉1 )
2 𝜌 2 𝜌
𝑃𝑖
ℎ𝑖 = 𝑢𝑖 +
𝜌
𝑣2 2 𝑣1 2
0 = [𝑔𝑦2 + + ℎ2 ] − [𝑔𝑦1 + + ℎ1 ]
2 2

Del sistema tenemos.

La entalpia es constante.
𝐴2 𝑉2
La velocidad inicial. 𝑉1 =
𝐴1

Respecto al nivel de referencia. 𝑦2 = 0

𝑣2 2
Por tanto. = 𝑔𝑦1
2

Despejando. 𝑣2 = √2𝑔𝑦1
De acuerdo al teorema de Torricelli para líquidos reales.
𝑣2 = 𝐶𝑑 √2𝑔𝑦1
Tomando un perfil del cono invertido.

𝑅 𝑟
=
𝐻 ℎ
𝑅
𝑟= 𝑥ℎ
𝐻
𝑑ℎ
= 𝑣1
𝑑𝑡
𝑑ℎ 𝐴2
= 𝑣2
𝑑𝑡 𝐴1
→ 𝐴1 = 𝐴 & 𝐴2 = 𝐴0
𝑑ℎ
𝐴 = 𝐶𝑑 √2𝑔ℎ 𝐴0
𝑑𝑡
𝜋(𝑅2 ) 𝑑ℎ ℎ2
( ) ( ) = 𝐶𝑑 √2𝑔 𝐴0 𝑑𝑡
4 √ℎ 𝐻 2
𝐻 𝑡
𝜋 (𝑅 2 )
( ) ∫ ℎ3⁄2𝑑ℎ = 𝐶𝑑 √2𝑔 𝐴0 ∫ 𝑑𝑡
4𝐻 2
0 0
2)
𝜋 (𝑅 2 5⁄2
( ) ( 𝐻 ) = 𝐶𝑑 √2𝑔 𝐴0 (𝑡𝑑 )
4𝐻 2 5

2𝜋𝑅2 √𝐻
𝑡𝑑 =
5𝐶𝑑 √2𝑔 𝐴0
a. Flujo laminar
Se llama flujo laminar o corriente laminar al movimiento de un fluido cuando este
es ordenado, estratificado o suave. En un flujo laminar, el fluido se mueve en
láminas paralelas sin entremezclarse y cada partícula de fluido sigue una
trayectoria suave, llamada línea de corriente
𝑉𝐷
𝑅𝑒 =
𝑣
𝑅𝑒 < 2000
𝑣 = 𝑣𝑖𝑠𝑐𝑜𝑠𝑖𝑑𝑎𝑑 𝑐𝑖𝑛𝑒𝑚𝑎𝑡𝑖𝑐𝑎
 El orificio del tanque y el diámetro del sistema influyen en el sistema como
lo menciona el número de Reynolds.
 Si el orificio es pequeño tendrá una tendencia a ser un flujo laminar.

b. Flujo de transición
Flujo en transición: es aquel en el que hay algunas fluctuaciones intermitentes
del fluido en un flujo laminar, aunque no es suficiente para caracterizar un flujo
turbulento. En la siguiente figura se representa por la imagen.
𝑉𝐷
𝑅𝑒 =
𝑣
2000 < 𝑅𝑒 < 3000
 El orificio del tanque y el diámetro del sistema influyen en el sistema como
lo menciona el número de Reynolds.
 Si el orificio aumenta de diamtro tendrá una tendencia a ser un flujo de
transición hasta cierto punto ya que su rango es muy pequeño.

c. Flujo turbulento
En mecánica de fluidos, se llama flujo turbulento o corriente turbulenta al
movimiento de un fluido que se da en forma caótica, en que las partículas se
mueven desordenadamente y las trayectorias
𝑉𝐷
𝑅𝑒 =
𝑣
 3000 < 𝑅𝑒
 El orificio del tanque y el diámetro del sistema influyen en el sistema como
lo menciona el número de Reynolds.
 Si el orificio aumenta de diámetro tendrá una tendencia a ser un flujo de
turbulento.

d. Si a la salida del orificio se aumenta un tubo de salida de longitud L

Al aumentar un tubo al tanque hace que las velocidades sean estables haciendo
que el coeficiente de descarga sea más adecuado para el sistema.