UNIVERSIDAD NACIONAL DEL CENTRO DEL

PERÚ
FACULTAD DE INGENIERÍA QUÍMICA

EAP INGENIERIA QUIMICA

MODELO MATEMÁTICO PARA DETERMINAR EL

TIEMPO DESCARGA DE UN TANQUE CILÍNDRICO
(6to SMESTRE)

CATEDRATICO Ing. GUEVARA YANQUI, Víctor Pascual

CATEDRA MODELAMIENTO Y SIMULACIÓN DE PROCESOS

SEMESTRE IX
Realizado por:
MENDOZA OSORIO, Michael Yemen IQ

HUANCAYO, PERU

2020
 RESUMEN

En este trabajo se presenta la simulación del tiempo de descarga de un tanque

cilíndrico que contiene agua a una altura inicial (16 cm) del nivel de referencia,
que es la salida del fluido; las otras variable consideradas son el diámetro del
cilindro (D,19.1 cm), el área del tubo de salida (0.50265), , tiempo de descarga (𝑡𝑑,
variable dependiente), la gravedad, el coeficiente de descarga (𝐶𝑑), la altura de
descarga (h, variable independiente) aplicando fenómenos de transporte se

2 ( √ y 0− √ y f ) A 1
obtiene el modelo matemático final t d=
√ 2 g Cd A2
 MARCO TEORICO

FLUIDO

Los fluidos, como todos los materiales, tienen propiedades físicas que permiten
caracterizar y cuantificar su comportamiento, así como distinguirlos de otros.
Algunas de estas propiedades son exclusivas de los fluidos y otras son típicas
de todas las sustancias. Propiedades como la viscosidad, tensión superficial y
presión de vapor solo se pueden definir en los líquidos y gases. Sin embargo la
masa específica, el peso específico y la densidad son atributos de cualquier
materia.[ CITATION DMH80 \l 10250 ]

Propiedades de los fluidos:

 Viscosidad: es una propiedad de los fluidos que se refiera el grado de
fricción interna; se asocia con la resistencia que presentan dos capas
adyacentes moviéndose dentro del fluido. Debido a esta propiedad
parte de la energía cinética del fluido se convierte en energía interna.
 Densidad: es la relación entre la masa y el volumen que ocupa, es
decir la masa de unidad de volumen.
 Volumen específico: es el volumen que ocupa un fluido por unidad de
peso.
 Peso específico: corresponde a la fuerza con que la tierra atrae a una
unidad de volumen.
 Gravedad específica: indica la densidad de un fluido respecto a la
densidad del agua a temperatura estándar. Esta propiedad es
dimensional.

Tipos de Fluidos.
Fluido newtoniano: Un fluido newtoniano es un fluido cuya viscosidad puede
considerarse constante en el tiempo. La curva que muestra la relación entre el
esfuerzo o cizalla contra su tasa de deformación es lineal y pasa por el origen.
El mejor ejemplo de este tipo de fluidos es el agua en contraposición al
pegamento, la miel o los geles que son ejemplos de fluido no newtoniano. Un
 buen número de fluidos comunes se comportan como fluidos newtonianos bajo
condiciones normales de presión y temperatura: el aire, el agua, la gasolina, el
vino y algunos aceites minerales.
Fluido no newtoniano: Es aquél cuya viscosidad varía con la temperatura y
presión, pero no con la variación de la velocidad. Estos fluidos se pueden
caracterizar mejor mediante otras propiedades que tienen que ver con la
relación entre el esfuerzo y los tensores de esfuerzos bajo diferentes
condiciones de flujo, tales como condiciones de esfuerzo cortante oscilatorio.

Ecuaciones matemáticas

Se presentan los modelos matemáticos para determinar los Coeficientes de

descarga, velocidad y contracción, y también para determinar el porcentaje de
error.
Para hallar estos coeficientes se requiere determinar el área del orificio, el área
del chorro contraído, la velocidad real, la velocidad teórica, el caudal real y el
caudal teórico.[ CITATION OLE93 \l 10250 ]
 Teorema de Bernoulli
Una forma especial de la ecuación de Euler derivada por una corriente natural
se llama la ecuación de Bernoulli.
1 1
p1 + ρ gy 1 + ρv 21= p 2 + ρ gy 2 + ρv22 + perdidas
2 2

Dónde:
p= Presión (N/M2)
ρ= Densidad (kg/m3)
g= Constante de
gravitación (m/s2)

 Ecuación de continuidad
y= La elevación (m)

Si dentro de un tubo el caudal o flujo de un fluido que va a una velocidad media

v1   es:
 Q=A1 v 1
Dónde:
A= Área de la sección transversal del tubo
El principio de conservación de la masa en dinámica de fluidos, para flujo en una
dirección es:

ρ1 v 1 A 1= ρ2 v 2 A2 + A ∂ ( ρv )
∂t
En estado estacionario el término de la derivada respecto del tiempo es cero. Un
fluido de densidad constante (como los líquidos) se denomina incompresible.
ρ1 v 1 A 1= ρ2 v 2 A2
Agua.

Las propiedades químicas del agua:

 Cuando las moléculas de agua se atraen unas a otras, se unen. Esta es la
razón del porqué se forma las gotas. Si no fuese por la gravedad de la
Tierra, una gota de agua tendría forma redonda.
 Al agua se le llama el "solvente universal" porque disuelve más substancias
que cualquier otro líquido. Esto significa que a donde vaya el agua, ya sea
a través de la tierra o a través de nuestros cuerpos, lleva consigo valiosos
químicos, minerales y nutrientes.
 El agua pura es neutral.

Propiedades Físicas del Agua:

 El agua es la única substancia natural que se encuentra en sus tres estados
líquidos, sólidos (hielo) y gaseosa (vapor) a las temperaturas encontradas
normalmente en la Tierra. El agua de la Tierra está cambiando
constantemente y siempre está en movimiento.
  El agua se congela a 0 grados Celsius (°C) y hierve a 100 °C (al nivel del
mar). Los puntos de congelamiento y ebullición son la base para medir la
temperatura: 0° En la escala Celsius está el punto de congelamiento del
agua, y 100° es el punto de ebullición del agua. El agua en su forma
sólida, hielo, es menos densa que en su forma líquida, por eso el hielo
flota.
 El agua tiene una tensión superficial muy alta. Esto significa que el agua es
pegajosa y elástica y tiende a unirse en gotas en lugar de separarse en
una capa delgada y fina.

DENSIDAD VISCOSIDAD
kg s
TEMPERATURA °C ρ( 3 ) μ( N . )
m m2
0 1000 1.75∗10−3
5 1000 1.52∗10−3
10 1000 1.30∗10−3
20 998 1.02∗10−3
30 996 8.00∗10− 4
40 992 6.51∗10−4
50 988 5.41∗10−4
60 984 4.60∗10−4
70 978 4.02∗10−4
80 971 3.50∗10−4
90 965 3.11∗10− 4
100 958 2.82∗10−4

VOLUMEN DE CONTROL

Un volumen de control es un espacio arbitrario que se sustituye con el objeto

de estudio. Formado por el espacio de un delimitado por una superficie de
control (SC) cerrada real o virtualmente, donde sus características en general,
será de permanencia en la forma y el tamaño así delimitado.
 PROCEDIMIENTO EXPERIMENTAL

DETERMINACIÓN DEL MODELAMIENTO

Volumen
de control

BALANCE DE MATERIA EN EL TANQUE

Ecuación general de Balance de Materia:

Para las condiciones de problema:

Generación de masa = 0
Consumo de masa = 0
Aplicando la expresión integral de la conservación de la masa:
Rapidez de Rapidez de Rapidez de Rapidez de Rapidez de
( )( )( )(
0= acumulacion − ingreso + la salida + Genración − consumo
de masa de masa de masa de masa de masa )( )
Aplicando la expresión integral de la conservación de la masa:
→ →
∫∫C .S ρ ( v .n) dA+ ∂t∂ ∫∫∫V .C ρdV=0
(1)
∂
∂t
∫∫∫V .C ρdV =0
Considerando que: (2)

→ →

Remplazando (11) en (10):

∫∫C .S ρ ( v .n) dA=0 (3)
 Dónde:
→ → → → → →
∫∫C .S ρ ( v .n) dA=∫∫A 1 ρ ( v .n) dA+∫∫A2 ρ ( v .n) dA=0 (4)

De acuerdo al sentido del fluido:

→ →
∫∫C.S ρ ( v .n) dA=−∫∫A 1 ρ vdA+∫∫A2 ρvdA=0 (5)

Integrando: ρ1 . v 1 . A 1= ρ2 . v 2 . A 2 (6)

Como es un fluido incompresible:

1   2 (7)

Remplazando (7) en (6): v 1 . A 1 =v 2 . A 2 (8)

A2
v 1 =v2 .
Despejando
A1 (9)

BALANCE DE ENERGÍA EN EL TANQUE:

∂Q ∂W p →→ ∂W
− =∫∫C .S e+ ρ ( v . n ) dA + ∂ ∫∫∫C. V eρ dV + μ
( )
∂t ∂ t ρ ∂t ∂t (10)
Considerando que:
∂
∂t
∫∫∫V .C eρdV =0
(Por ser permanente) (11)
∂W μ
=0
∂t (Por no tener fricción) (12)
Remplazando (11) y (12) en (10):

∂Q ∂W p →→
∂t ∂ t ρ ( )
− =∫∫C .S e+ ρ ( v .n ) dA
(13)
Dónde:
 p v2 p
e+ = +gy + μ+
ρ 2 ρ (14)
Remplazando (14) en (13):

∫∫C . S ( e+ ρp ) ρ ( v . n ) dA=∫∫A 1 ( e+ pρ ) ρ ( v . n ) dA+∫∫A 2 ( e+ ρp ) ρ ( v . n ) dA

→ → → → → →

v2 p v 21 p1
∫∫C . S ( e+
p ( → →)
ρ ) 2 [
ρ v . n dA= 2 + gy 2 + μ 2 + 2
ρ2 ] ( ρ 2 v 2 A2 )− [ 2
+ gy 1 + μ1 +
ρ1 ] ( ρ 1 v 1 A1 )
Remplazando en (13):
v 22 p2 v 21 p1
∂Q ∂W
∂t
−
∂t
=
2 [
+ gy 2 + μ2 +
ρ2 ] ( ρ2 v 2 A 2 ) − [ 2
+ gy 1 + μ 1 +
ρ1 ] ( ρ1 v 1 A 1 )

(15)
Asumiendo las siguientes condiciones:
∂Q
=0
No existe transferencia de calor: ∂t
(16)
∂W
=0
El sistema no realiza trabajo: ∂t
(17)
Remplazando (17), (16) en (15):
v 22 p2 v 21 p1
0=
[ 2
+ gy 2 + μ 2 +
ρ2 ] ( ρ2 v 2 A2 )− [ 2
+ gy 1 + μ1 +
ρ1 ] ( ρ1 v 1 A1 )

(18)
Quedando:
v 22 p2 v 21 p1
[ 2
+ gy 2 + μ2 +
ρ2 ] [
( ρ2 v 2 A 2 ) = 2 + gy 1 + μ 1 + ρ ( ρ1 v 1 A 1 )
1
]
(19)

Como el flujo másico es constante: ρ 2 v 2 A 2= ρ 2 v 2 A 2=m

(20)
Remplazando (20) en (19):
v 22 p2 v 21 p1
[ 2
+ gy 2 + μ2 +
ρ2 ][
=
2
+ gy 1 + μ1 +
ρ1 ]
(21)
 pi
hi =μi +
Sabemos que:
ρi (22)
2 2
v2 v1
+ gy 2 + h2 = + gy 1 +h1
Remplazando (22) en (21): 2 2
(23)
Del sistema, tenemos:

La entalpía es constante h1 =h2

(24)

La velocidad inicial v 1 =0

(25)

Respecto al nivel de referencia:

y 2 =0
(26)
Remplazando (24), (25) y (26) en (23):
2
v2
=gy 1
2 (27)

Despejando
v 2 =√ 2gy 1 (28)

HALLANDO EL TIEMPO DE DESCARGA

dy
=−v 1
dt

A 2 v2
v1 =
A1

A2
v1 =√ 2 g y 1 Cd
A1

2 ( √ y 0− √ y f ) A 1
t d=
√ 2 g Cd A2
DATOS EXPERIMENTALES

dato del cilindro

Atura inicial de agua en el tanque (H0) 10 cm
Diámetro del cilindro (D) 7.4 cm
Diámetro de la boquilla de salida (d) 0.5cm
Densidad del Líquido(ρ) 1 g/cm3
Gravedad(g) 981cm cm/s2

altura final del tiempo final de descarga(segundos)

N° agua(centímetros) promedio
corrida 1 corrida 2 corrida 3 de la
corridas
1 10 0 0 0 0
2 9 1,98 1,96 2,01 1,983
3 8 3,77 4,08 4,15 4,07
4 7 5,89 6,21 6,29 6,23
5 6 8,61 8,9 9,02 8,936
6 5 11,31 11,47 11,6 11,65
7 4 14,57 14,42 14,57 14,686

USO DEL PROGRAMA MATLAB PARA LA SIMULACIÓN DE TIEMPO DE

DESCARGA

A continuación, se muestran el modelamiento usando el programa Matlab, así como la

comparación de los datos experimentales registrados y los datos teóricos (obtenidos al
reemplazar los valores de h en el modelo matemático) para los tiempos de descarga
están incluidos en el programa.
clc,clf, clear all
format short
disp('--------------------------------------------------------')
disp(' SIMULACION DEL TIEMPO DE DESCARGA ')
%ingresar datos de la altura establecida
h=[10,9,8,7,6,5,4];
%determinacion de las areas del cilindro
d=0.5;
D=7.4;
Ao=(pi*(d^2))/4;
A=(pi*(D^2))/4;
%datos necesario para la determinacion del tiemo de descarga
H=10;%altura inicial en cm
g=981;%gravedad cm/s^2
Cd=1;%coeficiente de descarga
n=length(h);
disp('--------------------------------------------------------')
disp(' i h tdexperimental td teorico error(%)' )
disp('--------------------------------------------------------')
%INGRESANDO DATOS EXPERIMENTALES
hEXP=[10,9,8,7,6,5,4];
tdEXP=[0,1.98,4.11,6.25,8.96,11.53,14.49];
%HALLANDO EL TIEMPO DE DESCARGA
for i=1:n
td(i)=(2*((H^(1/2))-(h(i)^(1/2)))*A)/(((2*g)^(1/2))*Cd*Ao);
error(i)=abs((td(i)-tdEXP(i))/td(i))*100;
fprintf('%5d %8.1f %12.4f %12.4f %12.4f
\n',i,h(i),tdEXP(i),td(i),error(i));
end
plot(h,td,'*-r',hEXP,tdEXP,'-o'),title('SIMULACION DEL TIEMPO DE
DESCARGA'),xlabel('ALTURA(cm)'),ylabel('TIEMPO(s)'),grid,legend('TEORICO'
,'EXPERIMENTAL')

--------------------------------------------------------
SIMULACION DEL TIEMPO DE DESCARGA
--------------------------------------------------------
i h tdexperimental td teorico error(%)
--------------------------------------------------------
1 10.0 0.0000 0.0000 NaN
2 9.0 1.9800 1.6050 23.3680
3 8.0 4.1100 3.3018 24.4761
4 7.0 6.2500 5.1085 22.3443
5 6.0 8.9600 7.0496 27.0995
6 5.0 11.5300 9.1604 25.8682
7 4.0 14.4900 11.4951 26.0534
 CONCLUSIONES

 Se determinó los tiempos de descarga en un tanque cilíndrico a través de un

orificio de descarga.
El tiempo de descarga teórica es de: 14.49 seg
El tiempo de descarga experimental: 11.49 seg

Se debe a que tomamos el dato teórico de Cd 1

 REFERENCIA BIBLIOGRAFICA

 Chavez, L. M. (s.f.). VALIDACIÓN DE MODELO EXPERIMENTAL EN LA

DESCARGA DE UN TANQUE QUE CONTIENE AGUA. MEXICO.
 D.M.HIMMELBLAU, & K.B.BISCHOFF. (1980). "Análisis y de Simulación de
Procesos. Reverte‚. S.A.
 LEVENSPIEL, O. (1993). Fenomemos de fluidos e intercambiador de calor.
España.
 J.R WELTY, C.E WICKS Y R.E WILSON, ”Fundamentos de la transferencia de
Momento , Calor y Masa”, 1 ra reimpresión, Editorial LIMUSA, MEXICO 1983
ANEXO

LATA