vaciado de un tanque

TALLER N. 3 VACIADO DE UN TANQUE
PROFESOR: JOSÉ MIGUEL BERRIO
JULIETH PAOLA SANTANA ALARCÓN
1)

Para un tanque abierto:
BALANCE DE MASABALANCE DE ENERGÍA


Para el caso del tanque abierto, los cálculos son independientes de la presión atmosférica, debido a que tanto el punto 1 y el punto 2 están a la presión atmosférica, por lo tantono hay.
Para un tanque cerrado:
BALANCE DE MASA BALANCE DE ENRGIA



donde: nRT = Pgasinicial * Vgas inicial
En este caso so depende de la Patm,como afecta esta mencionado en los cálculos y será mencionado más adelante.
a)
!Vaciado de un tanque
!j=nRT
!dhdt=ecuacion

double precisionqi,hi,ro,h,dt,dtr,hO,dhdt,f,pa,n,re,j,pg,A,At,Vig,nv,p1
real g,u,t,pi
integer z

open(unit=10,file='resultados de tanque abierto')
open(unit=11,file='resultados de tanque cerrado')
print*,'ingrese el flujo de entrada en [Kg/s]'
read*,fprint*,'ingrese la densidad del liquido en [Kg/m3]'
read*,ro
print*,'ingrese la altura del tanque cilindrico[m]'
read*,hi
print*,'ingrese la altura a la que q se encuentra en el liquido [m]'
read*,hOprint*,'ingrese el diametro del tanque [m]'
read*,dt
print*,'ingrese el diametro de la tuberia [m]'
read*,dtr
print*,'ingrese tamaño de paso'
read*,u

g=9.8
pa=101325
pi=3.141592654A=(pi*dt*dt)/4
At=(pi*dtr*dtr)/4
Vig=(hi-hO)*A
qi=f/ro

h=hO
t=0
print*,'tanque abierto:1, tanque cerrado:2'
read*,z
select case(z)

case(1)
write(10,*)'resultados tanque abierto'write(10,*)'con un tamaño de paso de ',u
write(10,*)' t ',' ',' h'

write(10,*)t,h
n=5
do while (n>=0.00001)
nv=h
dhdt=(qi-At*sqrt(2*g*h))/A
h=h+u*dhdtn=abs(nv-h)
t=t+u
write(10,*)t,h
end do

case(2)
print*,'ingrese la presion del gas'
1 read*,pg
j=pg*Vig
re=(ro*g*hO)+pg
if(re>pa)then
write(11,*)'resultados...