UNIVERSIDAD NACIONAL DEL CENTRO DEL PERÚ

FACULTAD DE INGENIERÍA QUÍMICA

CARRERA PROFESIONAL DE INGENIERÍA QUÍMICA

DEPARTAMENTO ACADÉMICO DE INGENIERÍA

DETERMINCAIÓN DEL MODELO

MATEMÁTICO PARA CALCULAR EL
TIEMPO DE DESCARGA DE UN TANQUE
CILÍDRICO

PRESENTADO A:

Ms. GUEVARA YANQUI, Pascual Víctor

REALIZADO POR:

Arlette Canut N.

ALUMNO DEL IX SEMESTRE

HUANCAYO, 19-09-2017
RESUMEN

Para el presente informe, se empezó por determinar el modelo matemático para

determinar el tiempo de descarga de un tanque cilíndrico luego se procedió a
diseñar el equipo (tanque cilíndrico recto), y luego hizo las corridas. Nuestro
módulo tiene un diseño para ser utilizado también con recirculación, sus
dimensiones son altura de pruebas experimentales 0.145 m con un diámetro de
0.068m. El diámetro de la tubería de descarga es de 0.008m. A continuación, se
evaluó los valores de los datos obtenidos y se llevó a análisis para ver si cumplían
o no con los parámetros establecidos. Los análisis y los cálculos se trabajaron
comparando los datos experimentales con los teóricos, se utilizó un coeficiente de
descarga de 0.99985 y el diámetro de la tubería de descarga es 0.14 cm.
INDICE

RESUMEN..................................................................................................................... 2
INDICE........................................................................................................................... 3
INTRODUCCIÓN........................................................................................................... 4
OBJETIVOS................................................................................................................... 5
OBJETIVO PRINCIPAL..............................................................................................5
OBJETIVO ESPECÍFICOS.........................................................................................5
MARCO TEÓRICO........................................................................................................ 6
FLUIDOS.................................................................................................................... 6
TIPOS DE FLUIDOS...............................................................................................6
OTRAS DEFINICIONES.............................................................................................7
MODELAMIENTO TEÓRICO.....................................................................................7
PLANTEO DE ECUACIONES:................................................................................7
PARTE EXPERIMENTAL............................................................................................10
MATERIALES Y HERRAMIENTAS:.........................................................................10
METODOLOGÍA.......................................................................................................10
PROCEDIMIENTO...................................................................................................10
CÁLCULOS..............................................................................................................11
CONCLUSIONES........................................................................................................12
BIBLIOGRAFÍA............................................................................................................13
INTRODUCCIÓN

La mecánica del medio continuo tiene como finalidad estudiar los esfuerzos que se
manifiestan en el interior de los sólidos, líquido y gases. La mecánica de medio
continuo se originó con los estudios de Galileo y sus discípulos. Galileo planteó y
resolví ó los primeros problemas de resistencia de materiales en su libro “Discorsi e
dimostrazioni matematiche intorno a due nuo scienze”, impreso en 1638. Las dos
nuevas ciencias de ese entonces eran precisamente la mecánica de los sólidos
deformables y la cinemática de los proyectiles Benedetto Castelli y Evangelista
Torricelli, por su parte, se ocuparon del movimiento de los fluidos. Torricelli, en la obra
De “motu gravium naturaliter descendentium et proiectorum”, publicada en 1644, pudo,
con intuición realmente genial, deducir la ley de descarga de un líquido, a través de un
orificio practicado en un depósito, a partir de la ley de caída de los sólidos. Fue así
como se sentaron, casi simultáneamente, las bases de la mecánica del medio continuo
relacionada con sus dos objetivos principales: el sólido deformable y el fluido en
movimiento.

Dentro de la mecánica de fluidos tenemos el vaciado de tanques el cual es realizado

en las labores diarias de una industria, fábricas, empresas mineras, etc., de ahí
nuestro interés de determinar el tiempo de vaciado de un tanque. La forma de los
tanques pueden variar, entre ellos tenemos: cilindros, cubos, piramidales, cónicos,
trapezoidales, etc.

Para poder determinar el tiempo de vaciado se utilizará diferentes principios como son:
balance de masa, balance de movimiento y balance de energía. También
necesitaremos de ecuaciones auxiliares para hallar los coeficientes de velocidad, de
contracción y de descarga entre otros.
OBJETIVOS

OBJETIVO PRINCIPAL
 Determinar el modelo matemático para el cálculo del tiempo de descarga de un
tanque cilíndrico.

OBJETIVO ESPECÍFICOS
 Plantear las ecuaciones matemáticas que rigen este fenómeno y obtener la
relación finalizada
 Aprender cómo afectan algunas variables al fenómeno de descarga.
 Construir un módulo de vaciado de tanques para la parte experimental.
 Conocer las desviaciones experimentales respecto a los datos teóricos.
MARCO TEÓRICO

FLUIDOS
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.

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.
Es importante clasificar los fluidos no newtonianos en independientes del tiempo o
dependientes del tiempo. Algunos ejemplos de fluidos independientes del tiempo
son: el plasma sanguíneo, polietileno fundido, látex, almibares, adhesivos,
malazas y tintas. Los fluidos que dependen del tiempo son difíciles de analizar
porque su viscosidad aparente varía con el tiempo. Ejemplos de ellos son
petróleos crudos a temperaturas bajas, tinta para impresoras, nylon, ciertas
gelatinas, mezclas de harina y varias solucione de polímeros.

OTRAS DEFINICIONES
 Un orificio se define como una abertura por lo general redonda por la
cual fluye el fluido.
 La velocidad real de salida del chorro es menor que la teórica, pues
en la salida se presentan pérdidas por fricción. La relación entre las
velocidades real y teórica se llama coeficiente de velocidad (Cv).
 El caudal real es menor que el caudal teórico
 El chorro a la salida del orificio se contrae y en esta sección el chorro
se llama vena contraída. La relación entre el área de la sección
contraída y el área del orificio se llama coeficiente de contracción (Cc).
 Se define el coeficiente de descarga Cd como la relación entre el
caudal real y el caudal teórico.

MODELO DE TORRICELLI

Aplicando el teorema de Bernoulli en los

puntos 1 y 2, del diagrama ilustrado,
podemos escribir la siguiente expresión:

1 2 1 2
P1 + gρ. h1 + ρ .u 1 =P2+ gρ . h2 + ρ .u2
2 2

Donde  es la densidad del fluido, P1 y P2

son las presión de los puntos 1 y 2
respectivamente. De igual modo u1 y u2
designan las velocidades del fluido en los
puntos 1 y 2 receptivamente. La presión
en la interfase aire – agua superior (punto
Ilustración 1: MODELO DE
2) es la presión atmosférica (Patm = P2). BERNOULLI
BALANCE MACROSCÓPICO DE ENERGÍA MECÁNICA

El balance macroscópico de energía mecánica en estado no estacionario para flujo

isotérmico:

Ilustración 2: Volumen de control con flujo unidimensional a través de las fronteras

Para un sistema en estado estacionario y sin pérdidas debidas a la fricción:

NÚMERO DE REYNOLDS (Re)

FACTOR DE FRICCIÓN (f)

Flujo laminar (Re < 2,1x103)

Flujo Turbulento (2,1x103 < Re <105)

MODELAMIENTO TEÓRICO

MODELO FÍSICO

MODELO
MATEMÁTICO

SIMULACIÓN

SOLUCIÓN DEL
MODELO
EXPERIMENTACIÓN
COMPARACIONEXPE
MODELO FINAL

REVISIÓN COMPARACIÓN

MODELO FINAL

APLICACIONES
PARA EL CONTROL
Y OPERACIONES

Ilustración 3: DIAGRAMA DE FLUJO PARA EL MODELAMIENTO DE UN FENÓMENO

PRIMER MODELO MATEMÁTICO

PLANTEO DE ECUACIONES:

Por la ecuación de Bernoulli:

2
v2
z 1=z 2+ … … ...(1)
2g
Despreciando las pérdidas:
 H=z 2−z 1

v 2=√ 2 g(z 1−z 2 )=√ 2 gH … … … … .(2)

Considerando las pérdidas, la velocidad real de descarga a través de

la boquilla es:

v=C v √ 2 gH

C v : Coeficiente de velocidad.

Área de vena contractada:

Ac =C c A 0

Ac : Área de la vena contractada.

C c : Coeficiente de contracción.

A0 : Área de la conexión de salida.

Caudal descargado:
V
Q= = A c v =(C ¿ ¿ c A 0) ( C v √2 gH ) … … … .(3)¿
t
V : Volumen descargado.

Q=(C ¿ ¿ c C v ) A0 √ 2 gH ¿ … … … ..(4)

Por teoría

Q=C d A 0 √ 2 gH … … … …(5)

C d : Coeficiente de descarga.

Ahora determinamos el caudal como ecuación diferencial.

dV
=C d A 0 √ 2 gH … … … …(6)
dt
Diferenciando:
dV
=C d A 0 √ 2 gH
dt

dV =C d A 0 √2 gH dt … … …(7)

El volumen descargado también se puede expresar como:

dV =− AdH
A : Sección transversal del tanque.
A que es supuesta constante.
Reemplazando el diferencial de volumen en la ecuación se tiene:

− AdH =C d A 0 √ 2 gH dt

A dH
dt = … … ..(8)
A0 C d √2 gH

Donde:
2
πD
A=
4
D : Diámetro del tanque

Integrando entre los tiempos 0 y t y los niveles H 0 y H f , para áreas

transversales constantes se tiene:
A0Cd dH
dt = … … ..(9)
A √2 gH
4 A0Cd dH
dt=
πD
2
√2 gH
H =H f
4 A 0 C d t =t dH
πD 2 ∫ dt= ∫ √ 2 gH … … …..(10)
t =0 H =H 0

Resolviendo y despejando t .

π D ( √ H 0− √ H f )
2
t= … … .(11)
√ 8 A 0 Cd √ g

H 0 : Altura inicial que cubre el agua.

H f : Altura final transcurrido un tiempo t .

g : Gravedad.

Como la conexión de salida es circular entonces:

π d2
A0 =
4
Reemplazando en la ecuación de tiempo de descarga se tiene el
modelo matemático que predice el tiempo de descarga a diferentes
alturas.

2 D ( √ H 0− √ H f ) π . D ( √ H 0− √ H f )
2 2
t= = … … …(12)
d Cd √2 g
2
√ 8 . Cd . A 0 √ g
 D

H0

Ilustración 4: DISEÑO DEL TANQUE ESTUDIADO

DETERMINAMOS EL COEFICIENTE DESCARGA

De donde definimos:
1
n= ……………(13)
m

Para la recta realizamos un ajuste lineal:

y=mx+b …………..(14)

log Q=m. logh+b ………..(15)

Despejando Q:
Q=10 √h ………………….(16)
bn

Sabemos:
Q=Cd . A 2 . √n 2 gh ……………(17)
Evaluamos el caudal (Q):

Ideal.

Q=Cd I . A 2 . √n 2 gh ………….(18)
Dónde:

CdI = coeficiente de descarga ideal = 1

Remplazando (18) en (17):

b
10
Cd I = …………..(19)
A2 . √ 2 gh
n

Experimental.
Qe =Cd e . A 2 . √ 2 g h e …………….(20)
n

Dónde:
Cde = coeficiente de descarga experimental

Igualando (63) en la ec (60):

Qe =10b . √n he =Cd e . A 2 . √n 2 g h e
(65)

Despejando:
b
10
Cd e = ……………….
A2 . √ 2 g
n

(21)

Coeficientes de descarga (cd) en un orificio de pared delgada. [4]

Un orificio se define como una abertura por lo general redonda por la
cual fluye el fluido. La velocidad real de salida del chorro es menor
que la teórica, pues en la salida se presentan pérdidas por fricción.

Se define el coeficiente de descarga Cd como la relación entre el

caudal real y el caudal teórico.

a. Coeficiente de descarga (Cd).

CD = QR /QT …………………(22)
Dónde:
CD: Coeficiente de descarga
QR: Caudal real (m3).
QT: Caudal teórico (m3).

b. Caudal Real (Qr).

QR = V / t…………………(23)
Dónde:
QR: Caudal real (m3/s).
V: volumen (m3).
t: tiempo (s).
 c. Caudal Teórico (Qt)

QT = Ao ( 2g * h )½………….(24)
Dónde:
QT: Caudal teórico (m3/s).
A: Área (m2).
g: Gravedad (m/s2).
h: Altura piezométrica (m).

A. Diámetro de salida del fluido. [4]

La descarga del fluido para un cilindro atmosférico se desarrolla en
flujo constante, tomadas las consideraciones y ecuaciones
anteriores, el modelo matemático se presentara de la siguiente
manera:

dρv
ℜ=
μ
(70)

Reemplazando (28) en (67), se tiene:

dρ √ 2 gh ……………(25)
ℜ=
μ

ℜ. μ ……………….(26)
d=
ρ √ 2 gh

Para Re < 2100 por ser flujo laminar

SEGUNDO MODELO MATEMÁTICO

DETERMINANDO EL MODELO EMPÍRICO: “METODO π DE BUCKINGHAN”

Como t=f (h, D, d, e, ρ, u, g)

- Lista de variables influyentes y sus dimensiones:

N Variable Símbolo Dimensión

°
1 Tiempo T θ
2 Altura H L
3 Diámetro del tanque D L
4 Diámetro del orificio D L
5 Espesor del material E L
 6 Densidad del fluido Ρ M-1L-3
7 Viscosidad del fluido U ML-1 θ-1
8 Gravedad g L θ-2

- Dimensiones fundamentales: M; L; θ.
- Número de parámetros: 8 – 3 = 5
-Número de variables: 8
-Número de dimensiones elementales: 3
- Como tenemos 5 grupos adimensionales lo calculamos así:

[ θ ] =f {[ L ] , [ L ] , [ L ] , [ L ] , [ M L−3 ] , [ M L−1 θ−1 ] , [ L θ−2 ] }

a b c d e f g

- Como la ecuación debe ser dimensionalmente homogénea entonces:

[M] 0=e+f
[ θ] 1 = -f -2g
[ L] 0 = a + b + c + d – 3e –f -2g

- Nos produce un sistema de ecuaciones de 3 ecuaciones con 7

incógnitas, por lo que se escogen 4 variables (Que queramos que
se repitan en los diferentes grupos adimensionales). Y se pone en
función de los demás. Para nuestro caso escogeremos: el diámetro
(b), la viscosidad (f), la gravedad (g), la altura(a).
e = -f
e = 1- 2g
a = e –b-c-d +1
b = e – a- c- d +1

- Reemplazando en la relación :

[ θ ] =f {[ L ] }
1−e
e −e
, [ L ] , [ L] , [ M L ] , [M L θ ] , [ L θ ]
e−b−c−d +1 e−a−c−d+1 c d 2
,[ L]
−3 −1 −1 −2

- Agrupando las potencias obtenemos:

( )
2 2
u∗t ρ ∗h ∗g D d e
2
=f 2
, , , Relación Empírica
ρ∗h u h h h
PARTE EXPERIMENTAL

MATERIALES Y HERRAMIENTAS:
 Recipiente de Plástico y botella
 Bomba eléctrica.
 Manguera de nivel de 1 cm de diámetro.
 Silicona en barra
 Flexómetro
 Cronómetro

METODOLOGÍA
La metodología empleada es experimental.

PROCEDIMIENTO
1. Construcción del módulo de descarga de un tanque cilíndrico.
2. Se llena el tanque con agua.
3. Se deja correr el agua y se controla el tiempo cada 5 segundos y se mide la
altura respectiva.
4. Se anota los datos obtenidos.
5. Se procede al tratamiento de datos mediante la ecuación (12), del
modelamiento teórico.

CÁLCULOS
DATOS:
Diámetro del tanque=14cm
Diámetro de la tubería de descarga=0.8
Altura del tanque=14cm TABLA N°1:Datos experimentales

G=981cm2/s t h
Mediante la ecuación (12) 0 14
5 11
Datos experimentales: 10 9.4
15 7.8
20 6.2
25 4.8
30 3.7
35 2.5
40 1.3
45 0.5
50 0
Datos teóricos: TABLA N°2: Datos teóricos
T h
0 14
9.07864960
9 11
14.4332077
1 9.4
20.2664641
7 7.8
26.7356933
2 6.2
33.1242162
6 4.8
38.8348218
7 3.7
GRÁFICA: 46.1484392
4 2.5
55.5673691
7 1.3
64.8176678
7 0.5
79.9213911 0

DESCARGA DE TANQUES t vs h
16

14
f(x) = 0.00323543123543124 x² − 0.433407925407925 x + 13.5678321678322
12 R²
f(x)==0.997418905400747
− 0.271636363636364 x + 12.3545454545455
R² = 0.970568590958899
10
h(cm)

8 experimental
Linear (experimental)
6
Polynomial (exper-
imental)
4
teorico
2

0
0 10 20 30 40 50 60 70 80 90

t(s)
Se considero el coeficiente de descarga ideal:
Cd=1

DETERMINAMOS EL DIÁMETRO DE LA TUBERÍA DE DESCARGA PARA UN

FLUJO LAMINAR, MEDIANTE EL NUMERO DE REYNOLDS
TABLA N° 3: Datos Bibliográficos.

DATOS DE REFERENCIA TEORICA

Gravedad g 981 cm/s2
Densidad  0.999 1 g/cm3
Viscosidad µ 0.011 098 g/cm.s

ρ∙ v ∙ d
ℜ= <2 100
μ
TABLA N°4: Calculo del número de Reynolds para un diámetro de la tubería de
descarga=0.8 cm
  d 0.8  
      
h (cm) V (cm3) v (cm/s) Re
14 153.938 165.735 11725.578
13 142.942 159.706 11299.050
12 131.947 153.441 10855.776
11 120.951 146.908 10393.614
10 109.956 140.071 9909.922
9 98.960 132.883 9401.378
8 87.965 125.284 8863.704
7 76.969 117.192 8291.236
6 65.973 108.499 7676.193
5 54.978 99.045 7007.373
4 43.982 88.589 6267.585
3 32.987 76.720 5427.888
CORRECCIÓN DEL DIÁMETRO DE LA TUBERÍA DE DESCARGA
TABLA N°5: Calculo del número de Reynolds para un diámetro de la tubería de
descarga=0.14 cm

  D 0.14
       
h (cm) V (cm3) v (cm/s) Re
14 153.938 165.735 2051.976
13 142.942 159.706 1977.334
12 131.947 153.441 1899.761
11 120.951 146.908 1818.882
 10 109.956 140.071 1734.236
9 98.960 132.883 1645.241
8 87.965 125.284 1551.148
7 76.969 117.192 1450.966
6 65.973 108.499 1343.334
5 54.978 99.045 1226.290
4 43.982 88.589 1096.827
3 32.987 76.720 949.880

CALCULAMOS EL COEFICIENTE DE DESCARGA REAL

TABLA N°6: CALCULOS DE VELOCIDAD DE DESCARGA Y CAUDAL TEÓRICOS

h (m) t (s) v (m/s) V (m3) Caudal (m3/s) Ln h Ln Q

0.140 316.06655 1.65735 0.00216 0.000002551 -1.96611 -12.87891
0.130 293.49037 1.59706 0.00200 0.000002458 -2.04022 -12.91597
0.120 270.91418 1.53441 0.00185 0.000002362 -2.12026 -12.95599
0.110 248.33800 1.46908 0.00169 0.000002261 -2.20727 -12.99949
0.100 225.76182 1.40071 0.00154 0.000002156 -2.30259 -13.04715
0.090 203.18564 1.32883 0.00139 0.000002046 -2.40795 -13.09983
0.080 180.60946 1.25284 0.00123 0.000001929 -2.52573 -13.15872
0.070 158.03327 1.17192 0.00108 0.000001804 -2.65926 -13.22549
0.060 135.45709 1.08499 0.00092 0.000001670 -2.81341 -13.30256
0.050 112.88091 0.99045 0.00077 0.000001525 -2.99573 -13.39372
0.040 90.30473 0.88589 0.00062 0.000001364 -3.21888 -13.50529
0.030 67.72855 0.76720 0.00046 0.000001181 -3.50656 -13.64913

Ln(h)vsLn(Q)
Series2 Linear (Series2)
-12.60000
-3.20000 -3.00000 -2.80000 -2.60000 -2.40000 -2.20000 -2.00000 -1.80000
-12.70000
-12.80000
-12.90000
f(x) = 0.499999999999974 x − 11.8958558329238
-13.00000
LnQ

R² = 1
-13.10000
-13.20000
-13.30000
-13.40000
-13.50000
Lnh
b -11.896
A0 1.5394E-06
Cd 0.99985584

Mi coeficiente de descarga es 0.999855

CONCLUSIONES

Se determinó el modelo matemático para calcular el tiempo de descarga teórico dado

en la ecuación (12):

2 D 2 ( √ H 0− √ H f )
t= … … …(12)
d2 C d √ 2 g

Se puede decir que el tiempo depende de la altura en la que se encuentra el nivel del
líquido del tanque.
Se construyó un módulo para realizar los experimentos y compara los datos obtenidos
con los datos teóricos.
Se muestra una desviación luego del segundo 15 respecto a los datos teóricos, se
puede observar en la ilustración un aproximado a los datos.
Logramos calcular el diámetro corregido el cual es 0.14 cm.
Se calculó también el nuevo coeficiente de descarga que es 0.99985
BIBLIOGRAFÍA

[1.] Arlette Canut N. y otros. Descarga de tanques. M_exico D.F., 12 de marzo de

2008.

[2.] José Cruz Toledo Matus. MODELO MATEMATICO ARA EL DRENADO DE UN

TANQUE ATMOSFERICO. Artículo de Divulgación C. y Tecnológica.

[3.] Méndez Chávez, L. y otros. VALIDACIÓN DE MODELO EXPERIMENTAL EN

LA DESCARGA DE UN TANQUE QUE CONTIENE AGUA. México.

[4.] Fernández Larrañaga Bonifacio (1999). INTRODUCCION A LA MECANICA DE

FLUIDOS. México: Alfa omega Grupo Editorial.

[5.] Njock J, Mechanics of the slow draining of a large tank under gravity. Am. J.
Phys. 71 (11) November 2003, pp. 1204-1207

[6.] https://es.scribd.com/doc/92135563/DETERMINACION-DEL-VACIADO-
DE-UN-TANQUE-ANALISIS-Y-SIMULACION-doc