El Principio de Arquímedes

El principio de Arqumedes
El principio de Arqumedes afirma que todo cuerpo sumergido en un fluido experimenta un

empuje vertical y hacia arriba igual al peso de fluido desalojado.
La explicacin del principio de Arqumedes consta de dos partes como se indica en la

figuras:
1. El estudio de las fuerzas sobre una porcin de fluido en equilibrio con el resto del
fluido.
2. La sustitucin de dicha porcin de fluido por un cuerpo slido de la misma forma y
dimensiones.
Porcin de fluido en equilibrio con el resto del fluido.
Consideremos, en primer lugar, las fuerzas sobre una porcin de fluido en

equilibrio con el resto de fluido. La fuerza que ejerce la presin del fluido sobre la
superficie de separacin es igual a pdS, donde p solamente depende de la
profundidad y dS es un elemento de superficie.
Puesto que la porcin de fluido se encuentra en equilibrio, la resultante de las

fuerzas debidas a la presin se debe anular con el peso de dicha porcin de fluido.
A esta resultante la denominamos empuje y su punto de aplicacin es el centro de
masa de la porcin de fluido, denominado centro de empuje.
De este modo, para una porcin de fluido en equilibrio con el resto, se cumple
Empuje=peso=fgV
El peso de la porcin de fluido es igual al producto de la densidad del fluido f por

la aceleracin de la gravedad g y por el volumen de dicha porcin V.
Se sustituye la porcin de fluido por un cuerpo slido de
la misma forma y dimensiones.
Si sustituimos la porcin de fluido por un cuerpo slido de la misma forma y

dimensiones. Las fuerzas debidas a la presin no cambian, por tanto, su resultante
que hemos denominado empuje es la misma y acta en el mismo punto,
denominado centro de empuje.
Lo que cambia es el peso del cuerpo slido y su punto de aplicacin que es el

centro de masa, que puede o no coincidir con el centro de empuje.
Por tanto, sobre el cuerpo actan dos fuerzas: el empuje y el peso del cuerpo, que
no tienen en principio el mismo valor ni estn aplicadas en el mismo punto.
En los casos ms simples, supondremos que el slido y el fluido son homogneos y

por tanto, coinciden el centro de masa del cuerpo con el centro de empuje.
Ejemplo
Supongamos un cuerpo sumergido de densidad rodeado por un fluido de
densidad f. El rea de la base del cuerpo es A y su altura h.
La presin debida al fluido sobre la base superior es p1=fgx, y la presin debida al

fluido en la base inferior es p2=fg(x+h). La presin sobre la superficie lateral es
variable y depende de la altura, est comprendida entre p1 y p2.
Las fuerzas debidas a la presin del fluido sobre la superficie lateral se anulan. Las
otras fuerzas sobre el cuerpo son las siguientes:
Peso del cuerpo, mg
Fuerza debida a la presin sobre la base superior, p1A
Fuerza debida a la presin sobre la base inferior, p2A
En el equilibrio tendremos que
mg+p1A= p2A
mg+fgxA= fg(x+h)A
o bien,
mg=fhAg
Como la presin en la cara inferior del cuerpo p2 es mayor que la presin en la

cara superior p1, la diferencia es fgh. El resultado es una fuerza hacia
arriba fghA sobre el cuerpo debida al fluido que le rodea.
Como vemos, la fuerza de empuje tiene su origen en la diferencia de presin entre

la parte superior y la parte inferior del cuerpo sumergido en el fluido.
Con esta explicacin surge un problema interesante y debatido. Supongamos que

un cuerpo de base plana (cilndrico o en forma de paraleppedo) cuya densidad es
mayor que la del fluido, descansa en el fondo del recipiente.
Si no hay fluido entre el cuerpo y el fondo del recipiente desaparece la fuerza de

empuje?, tal como se muestra en la figura
Si se llena un recipiente con agua y se coloca un cuerpo en el fondo, el cuerpo

quedara en reposo sujeto por su propio peso mg y la fuerza p1A que ejerce la
columna de fluido situada por encima del cuerpo, incluso si la densidad del cuerpo
fuese menor que la del fluido. La experiencia demuestra que el cuerpo flota y llega
a la superficie.
El principio de Arqumedes sigue siendo aplicable en todos los casos y se enuncia

en muchos textos de Fsica del siguiente modo:
Cuando un cuerpo est parcialmente o totalmente sumergido en el fluido que le

rodea, una fuerza de empuje acta sobre el cuerpo. Dicha fuerza tiene direccin
hacia arriba y su magnitud es igual al peso del fluido que ha sido desalojado por el
cuerpo.
Flotacin entre dos lquidos no miscibles

Un cuerpo slido est sumergido en dos lquidos inmiscibles: agua y aceite.
Determinaremos la densidad de dicho cuerpo por dos mtodos distintos:
La ecuacin fundamental de la esttica de fluidos
El aceite que tiene una densidad 0.8 g/cm3 se sita en la parte superior y el agua
que es ms densa 1.0 g/cm3 se sita en la parte inferior del recipiente.
La densidad del bloque es un nmero comprendido entre la densidad del aceite

0.8, y la del agua 1.0. Un cuerpo de esta densidad flota entre los dos lquidos.
Principio de Arqumedes
Conociendo que parte del slido est sumergido en aceite (fluido 1) o en agua
(fluido 2), se determinar la densidad de dicho cuerpo.
El principio de Arqumedes nos dice que si el bloque est en equilibrio, el peso del
bloque debe ser igual al empuje proporcionado por ambos lquidos.
Peso del bloque =empuje del agua + empuje del aceite
Shg=2Sxg+1S(h-x)g
x=121hx=121h
S es el rea de la base del bloque, h su altura, y x es la parte del bloque

sumergida en agua.
En general, un cuerpo que flota en la superficie del agua, est sumergido en dos
fluidos, uno es el agua 2=1000 kg/m3 y el otro es el aire 1=1.29 kg/m3.
Habitualmente, se desprecia la densidad del aire frente a la del slido y la del
fluido 2, aplicndose la frmula aproximada
x2hx2h
Ejemplo
Supongamos un bloque de h=20 cm de altura, observamos que el bloque est

sumergido 13 cm en aceite y 7 cm en agua.
20=0.813+1.07
Despejando en la frmula la densidad del slido, obtenemos el valor de 0.87

g/cm3.
Ecuacin fundamental de la esttica de fluidos

Sean p1 y p2 las presiones sobre la cara superior e inferior del bloque sumergido.
La cara superior est en el aceite a una profundidad y. La presin p1 ser igual a la

atmosfrica p0 ms la correspondiente a la altura y de aceite.
p1=p0+1gy
La cara inferior est en el agua. La presin p2 ser igual a la presin

atmosfrica p0 ms la correspondiente a la altura de aceite (y+h-x) ms la
correspondiente a la altura de la columna de agua (x)
p2=p0+1g(y+h-x)+2gx
La fuerza que ejerce el fluido sobre dichas caras ser el producto de la presin por
el rea de su superficie S.
Como vemos en la figura, para que haya equilibrio se tiene que cumplir que
p1S+mg=p2S
Teniendo en cuenta que m=hS despejamos el valor de x.
x=121hx=121h
Es el mismo resultado que hemos obtenido aplicando el principio de Arqumedes
Ejemplo:
Datos: densidad del agua 1000 kg/m3, densidad del aceite 800 kg/m3, densidad
del mercurio 13550 kg/m3
La cara superior est a 22 cm de la superficie libre, la presin debida al aceite es
p1=0.229.8800=1724.8Pa1724.8135509.8=0.013m=20.65cm en el manmetrop1=0.
229.8800=1724.8Pa1724.8135509.8=0.013m=20.65cm en el manmetro
La cara inferior est a 42 cm de la superficie libre (35 cm de aceite y 7 cm de

agua), la presin debida a ambos lquidos es:
p2=8009.80.35+10009.80.07=3430Pa3430135509.8=0.026m=21.3cm en el manm
etrop2=8009.80.35+10009.80.07=3430Pa3430135509.8=0.026m=21.3cm e
n el manmetro
En el equilibrio se cumple
1724.8S+sS0.29.8=3430S1724.8S+sS0.29.8=3430S
Se obtiene s=870 kg/m3 0.87 g/cm3
El principio de
Arqumedes.Energa
potencial mnima
En este apartado, se estudia el principio de Arqumedes desde el punto de vista energtico.
Supongamos un cuerpo en forma de paraleppedo de altura h, seccin A y de
densidad s. El fluido est contenido en un recipiente de seccin S hasta una
altura b. La densidad del fluido es f> s.
Se libera el cuerpo, oscila hacia arriba y hacia abajo, hasta que alcanza el
equilibrio flotando sobre el lquido sumergido una longitud x. El lquido del
recipiente asciende hasta una altura d. Como la cantidad de lquido no ha
variado Sb=Sd-Ax
d=b+ASxd=b+ASx
Hay que calcular x, de modo que la energa potencial del sistema formado por el
cuerpo y el fluido sea mnima.
Tomamos el fondo del recipiente como nivel de referencia de la energa potencial.
El centro de masa del cuerpo se encuentra a una altura d-x+h/2. Su energa

potencial es Es=(sAh)g(d-x+h/2)
Es=sAhg(b+ASxx+h2)Es=sAhg(b+ASxx+h2)
Para calcular el centro de masas del fluido, consideramos el fluido como una figura
slida de seccin S y altura d a la que le falta una porcin de seccin A y altura x.
El centro de masas de la figura completa, de volumen Sd es d/2

El centro de masas del hueco, de volumen Ax, est a una altura (d-x/2)
yf=(Sd)(d/2)(Ax)(dx/2)SdAx=1Sb{12S(b+ASx)2+Ax22Ax(b+ASx)}yf=(Sd)(
d/2)(Ax)(dx/2)SdAx=1Sb{12S(b+ASx)2+Ax22Ax(b+ASx)}
La energa potencial del fluido es Ef=f(Sb)gyf

La energa potencial total es Ep=Es+Ef
EP=12fgA(1AS)x2sghA(1AS)x+12fSb2g+shA(b+h2)g=12fgA(1AS){x22sfhx
}+cteEP=12fgA(1AS)x2sghA(1AS)x+12fSb2g+shA(b+h2)g=12fgA(1A
S){x22sfhx}+cte
El valor de la constante aditiva cte, depende de la eleccin del nivel de referencia

de la energa potencial.
En la figura, se representa la energa potencial Ep(x) para un cuerpo de

altura h=1.0, densidad s=0.4, parcialmente sumergido en un lquido de
densidad f=1.0.
La funcin presenta un mnimo, que se calcula derivando la energa potencial con

respecto de x e igualando a cero
dEpdx=12fgA(1AS){2x2sfh}=0dEpdx=12fgA(1AS){2x2sfh}=0
En la posicin de equilibrio, el cuerpo se encuentra sumergido
xe=sfhxe=sfh
Energa potencial de un cuerpo que se

mueve en el seno de un fluido
Cuando un globo de helio asciende en el aire actan sobre el globo las siguientes
fuerzas:
El peso del globo Fg=mgj .

El empuje Fe= fVgj, siendo f la densidad del fluido (aire).
La fuerza de rozamiento Fr debida a la resistencia del aire
Dada la fuerza conservativa, la expresin de la energa potencial asociada, se

obtiene integrando
ABFdr=EpAEpB Ep=Ep(y)ABFdr=EpAEpB Ep=Ep(y)
La fuerza conservativa peso Fg=mgj est asociada con la energa

potencial Eg=mgy.
Por la misma razn, la fuerza conservativa empuje Fe= Vg j est asociada
a la energa potencial Ee=-fVgy.
Dada la energa potencial, obtenemos la expresin de la fuerza conservativa,

derivando
Fy=dEp(y)dyFy=dEp(y)dy
La energa potencial asociada con las dos fuerzas conservativas es
Ep=(mg- fVg)y
A medida que el globo asciende en el aire con velocidad constante experimenta

una fuerza de rozamiento Fr debida a la resistencia del aire. La resultante de las
fuerzas que actan sobre el globo debe ser cero.
f Vg- mg-Fr=0
Como fVg> mg a medida que el globo asciende su energa
potencial Ep disminuye.
Empleando el balance de energa obtenemos la misma conclusin
ABFncdr=(Ek+Ep)B(Ek+Ep)AABFncdr=(Ek+Ep)B(Ek+Ep)A
El trabajo de las fuerzas no conservativas Fnc modifica la energa total (cintica

ms potencial) de la partcula. Como el trabajo de la fuerza de rozamiento es
negativo y la energa cintica Ekno cambia (velocidad constante), concluimos que
la energa potencial final EpB es menor que la energa potencia inicial EpA.
Energa potencial de un cuerpo

parcialmente sumergido
En el apartado anterior, estudiamos la energa potencial de un cuerpo totalmente
sumergido en un fluido (un globo de helio en la atmsfera). Ahora vamos a
suponer un bloque cilndrico que se sita sobre la superficie de un fluido (por
ejemplo agua).
Pueden ocurrir dos casos:
Que el bloque se sumerja parcialmente si la densidad del cuerpo slido es

menor que la densidad del fluido, s< f.
Que el cuerpo se sumerja totalmente si s f.
Cuando el cuerpo est parcialmente sumergido, sobre el cuerpo actan dos fuerzas
el peso mg=sShg que es constante y el empuje fSxg que no es constante. Su
resultante es
F=(-sShg+fSxg)j.
Donde S el rea de la base del bloque, h la altura del bloque y x la parte del
bloque que est sumergida en el fluido.
Tenemos una situacin anloga a la de un cuerpo que se coloca sobre un muelle

elstico en posicin vertical. La energa potencial gravitatoria mgy del cuerpo
disminuye, la energa potencial elstica del muelle kx2/2 aumenta, la suma de
ambas alcanza un mnimo en la posicin de equilibrio, cuando se cumple
mg+kx=0, cuando el peso se iguala a la fuerza que ejerce el muelle.
Ep=mgy+12kx2=mgy+12k(y0y)2Fy=dEpdy=mg+k(y0y)=mg+kxEp=mgy+12kx2
=mgy+12k(y0y)2Fy=dEpdy=mg+k(y0y)=mg+kx
El mnimo de Ep se obtiene cuando la derivada de Ep respecto de y es cero, es

decir en la posicin de equilibrio.
La energa potencial del cuerpo parcialmente sumergido ser, de forma anloga
Ep=sShgy+12fSgx2=sShgy+12fSg(y0y)2Fy=dEpdy=sShg+fSg(y0y)=sShg
+fSgxEp=sShgy+12fSgx2=sShgy+12fSg(y0y)2Fy=dEpdy=sShg+fSg
(y0y)=sShg+fSgx
El mnimo de Ep se obtiene cuando la derivada de Ep respecto de y es cero, es

decir, en la posicin de equilibrio, cuando el peso se iguale al empuje. -
sShg+fSxg=0
x=shf=hx=shf=h
El bloque permanece sumergido una longitud x. En esta frmula, se ha

designado como la densidad relativa del slido (respecto del fluido) es decir, la
densidad del slido tomando la densidad del fluido como la unidad.
Fuerzas sobre el bloque
1. Cuando <1 o bien s<f, el cuerpo permanece parcialmente sumergido en

la situacin de equilibrio.
2. Cuando >1 o bien s>f, el peso es siempre mayor que el empuje, la
fuerza neta que acta sobre el bloque es
Fy=-sShg+fShg<0.
No existe por tanto, posicin de equilibrio, el bloque cae hasta que llega al
fondo del recipiente que supondremos muy grande.
3. Cuando =1 o bien s=f, el peso es mayor que el empuje mientras el

bloque est parcialmente sumergido (x<h).
Fy=- Shg+ Sxg<0.
La fuerza neta que acta sobre el bloque cuando est completamente

sumergido (x h) es cero, y cualquier posicin del bloque, completamente
sumergido en el seno del fluido, es de equilibrio.
Curvas de energa potencial
1. La energa potencial correspondiente a la fuerza conservativa peso es
Eg=sShgy
2. La energa potencial correspondiente a la fuerza de empuje tiene dos partes
o Mientras el cuerpo est parcialmente sumergido (x<h)
Ef=12fSgx2=12fSg(y0y)2Ef=12fSgx2=12fSg(y0y)2
Que corresponde al rea del tringulo de la figura de la izquierda.
o Cuando el cuerpo est totalmente sumergido (xh)
Ef=12fSgh2+fSgh(xh)=12fSgh2+fSgh(y0y)Ef=12fSgh2+fSg
h(xh)=12fSgh2+fSgh(y0y)
Que corresponde a la suma del rea de un tringulo de base h, y la

de un rectngulo de base x-h.
3. La energa potencial total es la suma de las dos contribuciones
Ep=Eg+Ef
Cuando la densidad del slido es igual a la del fluido s=f, la energa

potencial total Epes constante e independiente de x (o de y) para x h como
puede comprobarse fcilmente.
Actividades
Se introduce
La densidad del slido relativa al fluido en el control titulado Densidad
relativa.
Se pulsa el botn titulado Nuevo.
El bloque tiene una altura h=1 y una seccin S. Se coloca el bloque justamente
encima de la superficie del fluido. La altura de su centro de masas es y0=1.5
unidades.
Se suelta el bloque, llega hasta la posicin final de equilibrio ye=h, si la

densidad <1, o hasta el fondo del recipiente si la densidad >1.
El programa interactivo no hace ninguna suposicin acerca del modo en el que el

bloque parte de la posicin inicial y llega a la posicin final (no calcula la posicin y
velocidad del cuerpo en cada instante), ya que el objetivo del programa es el de
mostrar los cambios en la energa potencial Ep del cuerpo con la posicin y del c.m.
del mismo.
En la parte derecha, se traza:
la energa potencial debida a la fuerza conservativa peso Eg (en color

negro),
la energa potencial debida al empuje Ef (en color azul)
la suma de ambas contribuciones Ep (en color rojo) en funcin de la
posicin y del c.m. del bloque
La representacin de la energa potencial gravitatoria Eg (en color negro) es una

recta cuyo valor mximo est en la posicin inicial y=1.5 y es cero cuando el
bloque llega al fondo y=0.
La curva de la energa potencial correspondiente al empuje Ef (en color azul) es

algo ms complicada y consta de dos partes: Una parbola mientras el cuerpo est
parcialmente sumergido (x<h) (y>0.5), unida a una lnea recta cuando el cuerpo
est completamente sumergido (x h) (y 0.5). La energa potencial inicial es
cero y se va incrementando a medida que el cuerpo se sumerge en el fluido.
La curva de la energa potencial total Ep (en color rojo) es la suma de las dos
contribuciones, Ep=Eg+Ef
Para trazar estas grficas se ha tomado como unidad de energa, la energa

potencial inicial del bloque sShgy0 con y0=1.5, h=1 y s=, densidad del slido
relativa al fluido f=1. De este modo, la energa potencial inicial del bloque es una
unidad.
Se presentan tres casos:
1. Cuando <1, la energa potencial presenta un mnimo en x=h. En este

caso con x=y0-y, h=1 e y0=1.5, tendremos que la posicin del c.m. en el
equilibrio ser ye=1.5-.
2. Cuando >1, la curva de la energa potencial no tiene mnimo y por tanto,
no hay posicin de equilibrio estable.
3. En el caso lmite en el que =1 la energa potencial para y0.5 es una lnea
recta horizontal y la posicin de equilibrio del c.m. del bloque puede ser
cualquier y0.5.
Medida de la densidad de un
lquido y de un slido
En este pgina, se explica el funcionamiento de un aermetro, dispositivo para medir la
densidad de un lquido, mediante un modelo simple. Se trata de un ejercicio sencillo de
aplicacin del principio de Arqumedes.
Una aplicacin del principio de Arqumedes, la medida de la densidad de un slido se

encuentra en el captulo Unidades y Medidas. En esta pgina, vamos a estudiar otra versin
de dicha experiencia
Medida de la densidad un lquido
Conocida la masa del cuerpo y el volumen de la parte sumergida, determinamos la

densidad del lquido. En esto se basan los aermetros o flotadores de masa
conocida que se sumergen en el lquido de densidad desconocida. Disponen de una
escala graduada, que nos proporcionan mediante lectura directa la densidad del
lquido. La superficie libre del lquido marca el valor de la densidad en la escala del
aermetro.
Dependiendo de la aplicacin concreta los aermetros reciben nombres

especficos: alcohmetros, sacarmetros, etc.
El aermetro de la fotografa es un tubo de vidrio, que dispone de un lastre

formado por bolitas de plomo (en extremo izquierdo). En la parte derecha, dispone
de una escala graduada que mide directamente la densidad del lquido en g/cm3.
Laboratorio de Fsica. Escuela de Ingeniera de Eibar
Para un cuerpo en equilibrio que flota sobre la superficie de un lquido, tenemos
que el peso es igual al empuje
mg=fVg
Nuestro aermetro es un slido de forma cilndrica de 25 cm de altura y densidad

0.5 g/cm3 que se sumerge parcialmente en el lquido cuya densidad se quiere
determinar.
Midiendo en la escala graduada la parte del cilindro que est sumergida,

determinamos la densidad del fluido.
El cuerpo est en equilibrio flotando en el lquido, bajo la accin de dos fuerzas, su

peso y el empuje del fluido.
peso=empuje
sgSh=fgSx
sh=f x
Donde s es la densidad del cuerpo slido, S su seccin, h su altura. f es la

densidad del fluido y x la parte del slido que est sumergido en el lquido.
Ejemplo:
Sea agua el lquido de densidad desconocida
Observamos que el cilindro se sumerge hasta una altura x=12.5 cm

0.525=f 12.5
Se despeja f =1.0 g/cm3
Actividades
Se selecciona el fluido cuya densidad deseamos conocer en la lista de

lquidos: agua, aceite, alcohol, glicerina.
La densidad del slido se ha fijado en s=0.5 g/cm3
La altura del cuerpo slido de forma cilndrica se ha fijado en h=25 cm
Se lee en la escala la longitud x del cuerpo cilndrico que est sumergido y se

calcula densidad del lquido f.
Aceite
Lquidos:
Medida de la densidad de un slido

En la pgina titulada La balanza. Medida de la densidad de un slido se ha
estudiado el procedimiento de medida de la densidad de un cuerpo slido, con una
balanza ordinaria:
1. Se pesa el cuerpo slido y se anota la lectura en gramos m=F de los

indicadores de la balanza cuando su brazo est horizontal en equilibrio.
2. Se sumerge el cuerpo completamente en agua (densidad 1.0 g/cm3) y se
anota la lectura F
La diferencia V=F-F es el volumen del cuerpo expresado en cm3.
La densidad es el cociente =m/V (g/cm3)
En esta seccin, se mide la densidad >1.0 de un cuerpo slido por otro

procedimiento, empleando una balanza electrnica.
1. Llenamos un recipiente parcialmente de agua y lo colocamos sobre una
balanza electrnica, poniendo a cero el indicador del peso.
2. Atamos un hilo a un cuerpo slido, lo sumergimos completamente en agua
pero sin tocar el fondo, anotamos la lectura de la balanza V.
3. Soltamos el cuerpo, dejando de tirar del hilo, de modo que descanse en el
fondo del recipiente, anotamos la lectura de la balanza m.
La densidad del cuerpo es el cociente entre la segunda lectura y la

primera. =m/V
Cuando un cuerpo est suspendido en el fluido completamente sumergido pero sin

tocar el fondo, las fuerzas sobre el cuerpo son:
La tensin del hilo T

El peso mg
El empuje E=fgV, siendo f=1.0 g/cm3 la densidad del lquido (agua)
El cuerpo est en equilibrio, suspendido del hilo T=mg-fgV
El lquido ejerce una fuerza E=fgV sobre el cuerpo cilndrico, por la tercera ley
de Newton, el cuerpo ejerce una fuerza igual y de sentido contrario sobre el lquido
La balanza electrnica mide la fuerza E1=fgV

Cuando el cuerpo descansa en el fondo. La balanza electrnica mide el peso del
cuerpo F2=mg=Vg
La densidad del slido es el cociente F2/F1, si el lquido es agua f=1.0 g/cm3
Actividades
Se selecciona el slido cuya densidad deseamos conocer en la lista de

materiales: hierro, oro, plomo, etc...
La densidad del lquido (agua) se ha fijado en f=1.0 g/cm3
La altura del cuerpo slido de forma cilndrica se ha fijado en h=1.0 cm
Se introduce el slido en el agua, actuando en el control titulado Posicin. Se

simula una cuerda que sostiene el bloque cuya densidad es mayor que la del
fluido. Cuando llega al fondo, se deja de tirar de la cuerda (desaparece), el bloque
descansa en el fondo del recipiente.
Se calcula la densidad del slido a partir de la lectura de la balanza electrnica:
1. Cuando el cuerpo est suspendido del hilo en el lquido, completamente

sumergido
2. Cuando el cuerpo descansa en el fondo del recipiente
Nota: cuando se sumerge el slido, el nivel del lquido en el recipiente debera

cambiar. Esto no afecta a la medida de la densidad, que es independiente de la
altura del nivel siempre que haya suficiente agua para que el cuerpo se encuentre
completamente sumergido.
Ecuacin fundamental de la
esttica de fluidos
En esta pgina, se explica el concepto de presin y se formula la ecuacin fundamental de
la esttica de fluidos
La materia ordinaria se presenta en alguno de los tres estados siguientes: slido,

lquido o gaseoso. Existe un cuarto estado de la materia denominado plasma que
es esencialmente un gas ionizado con igual nmero de cargas positivas que
negativas.
Un slido cristalino es aqul que tiene una estructura peridica y ordenada, como
consecuencia, tiene una forma que no cambia, salvo por la accin de fuerzas
externas. Cuando se aumenta la temperatura, los slidos se funden y cambian al
estado lquido. Las molculas ya no permanecen en posiciones fijas, aunque las
interacciones entre ellas sigue siendo suficientemente grande para que el lquido
pueda cambiar de forma sin cambiar apreciablemente de volumen, adaptndose al
recipiente que lo contiene.
En el estado gaseoso, las molculas estn en continuo movimiento y la interaccin

entre ellas es muy dbil. Las interacciones tienen lugar, cuando las molculas
chocan entre s. Un gas se adapta al recipiente que lo contiene pero trata de
ocupar todo el espacio disponible.
Poner un barco a flote

Un barco est hundido a cierta profundidad, trataremos de ponerlo a flote, inyectando aire
para desalojar el agua que contiene. El barco empieza a ascender cuando el empuje se hace
igual al peso. Se calcular el volumen y la masa de aire a presin atmosfrica que tenemos
que suministrar accionado un compresor. Calcularemos adems, la distancia que recorre el
barco desde el fondo hasta situarse flotando en la superficie del agua.
En el anlisis de esta situacin fsica se pone en relacin
La ecuacin fundamental de la esttica de fluidos
La transformacin isoterma de un gas ideal.
Consideremos un barco que tiene la forma de una caja rectangular (paraleppedo)

de 10 m de altura y de seccin 280 m2.
El barco hundido
El barco empieza a flotar cuando el peso del barco sea igual al empuje. El empuje
es el peso del volumen de agua desalojada
mg=Sxg
Donde es la densidad del agua de mar (tomaremos 1000 kg/m3) aunque que es
algo mayor que la del agua dulce, S el rea de la base del barco y x la altura de
aire en el barco, en color amarillo en la figura.
Ejemplo: para un buque de 863 Tm, el valor de x es aproximadamente 3.08 m.

Lo que corresponde a un volumen de aire de 862 m3 (el rea de la base del barco
es de 280 m2)
La presin del aire ser
p=p0+g(y-h+x)
Donde y es la profundidad, h la altura del barco (10 m). Para una

profundidad y=26 m la presin vale p=287000 Pa y el volumen de aire que hemos
calculado es V=862 m3
El aire que ha de suministrar el compresor es aproximadamente igual al que

ocupara el aire comprimido en el barco a la presin atmosfrica (p0=105 Pa) a la
misma temperatura.
p0V0=pV 100000V0=287000862
Se obtiene un volumen de aire de V0=2475 m3. El tiempo que tarda el compresor

en bombear este volumen depende del caudal. Si el caudal es 1000 m3/min,
entonces tarda 2.5 min.
Calculamos la masa de aire contenida en dicho volumen a la temperatura
ambiente (300 K)
pV=mMRTpV=mMRT
donde M es el peso molecular del aire 28.9 g/mol. Se obtiene 2867 kg, que es una
masa pequea comparada con la del barco 863 Tm=863000 kg.
Ascendiendo hasta la superficie
Cuando el empuje se iguala al peso, el barco inicia su ascenso. Mientras asciende

la presin del aire disminuye y aumenta su volumen, el empuje se hace ms
grande que el peso. Puede ocurrir que el volumen de aire se haga mayor que el del
barco, entonces el barco pierde aire. El empuje tomar su valor mximo y se
mantendr constante hasta que el barco empieza a salir por la superficie del mar.
Los detalles del movimiento del barco, la aceleracin y la velocidad del barco
mientras asciende no son de inters en esta explicacin.
Flotando en la superficie del mar

El barco se para al llegar a la superficie, cuando el peso se vuelve a igualar al
empuje. En ese momento, la diferencia entre el nivel de agua dentro y fuera del
barco ser de z=3.08 m hallada anteriormente. La presin del aire ser
p=p0+gz
p=130184 Pa. A esta presin, el volumen V ocupado por el aire (transformacin

isoterma) ser
287000862=130184V
El resultado es V=1900 m3. Lo que corresponde a una altura de aire en el barco de

6.8 m (el rea de la base del barco es de 280 m2).
El desplazamiento del barco ser la distancia entre su extremo inferior y el fondo,

es decir,
26+6.8-10-3.08=19.72 m
26 m es la profundidad y 10 m es la altura del barco. El resultado es casi 20 m

como vemos en la figura.
Actividades
Se introduce
El peso del barco (en toneladas), en el control titulado Peso del barco.
La profundidad a la que se encuentra hundido el barco, en el control
titulado Profundidad.
El caudal de aire que suministra el compresor, medido en metros cbicos de
aire (a presin atmosfrica) por minuto, en el control titulado Compresor
de aire.
Se pone en marcha el compresor que suministra el aire que desaloja de agua el

barco. El empuje aumenta, hasta que se hace igual al peso, en ese momento el
barco se pone en movimiento ascendente hacia la superficie del mar.
Si el peso o la profundidad es excesiva, puede ocurrir que el barco desaloje todo el

agua, pero el empuje sea insuficiente para igualar al peso, el aire se escapa y el
barco es imposible de ponerlo a flote.
A medida que asciende el barco, la presin disminuye, y el volumen aumenta.

Puede ocurrir que el volumen de aire se haga mayor que el volumen del barco, el
aire se escapa y el empuje se hace constante.
Equilibrio de una varilla

parcialmente sumergida
En esta pgina, se estudia el equilibrio de una varilla parcialmente sumergida que est
sujeta por una de sus extremos.
Consideremos una varilla delgada, de seccin uniforme A y de longitud L y densidad <1

que est sujeta por su extremo P, pero que le permite moverse en el plano vertical. Vamos
a analizar la situacin de equilibrio cuando una parte de la varilla est sumergida en agua
tal como se muestra en la figura.
El punto de sujecin est por encima de la

superficie del agua
Supongamos que el extremo P de la varilla est a una altura y>0 sobre la
superficie del agua y en esa posicin la varilla hace un ngulo con el eje vertical.
Las fuerzas que actan sobre la varilla son:
El peso mg=ALg, acta en el centro de masa de la varilla, en la posicin

cuya abscisa es
xg=L/2sin.
La fuerza de empuje, E=gA(L-y/cos) acta en el centro de la porcin de

varilla sumergida, en la posicin cuya abscisa es
( (
xe= ycos+12 Lycos )) ( )
sin=12 L+ycos sin
La reaccin N se aplica en el punto P.
La barra estar en equilibrio:
Si la resultante de todas las fuerzas es cero,
N+E=mg
Si el momento de las fuerzas respecto del punto P es cero.
Exe-mgxg+N0=0
(
gA Lycos )( )
L+ycos sin2gALL2sin=0
Simplificando, nos queda la ecuacin
( 1y2L21cos2 ) sin=0
Si y>L el primer factor entre parntesis es distinto de cero. La posicin de
equilibrio se obtiene cuando sin=0, es decir, =0. La barra cuelga verticalmente
del punto P, sin estar sumergida
Cuando y<L el primer trmino puede hacerse cero cuando
cos=yL111
y como el coseno tiene que ser menor o igual que la unidad se tiene que cumplir a
la vez que yL1 Cuando no se cumple esta condicin el primer factor no es nulo
y la posicin de equilibrio es =0.
El punto de sujecin est por debajo de la

superficie del agua
Sea y la distancia (positiva) del punto P de sujecin de la varilla al origen O
El peso mg=ALg, acta en el centro de masa de la varilla, en la posicin

cuya abscisa es
xg=L/2sin.
La fuerza de empuje, E=gA(y/cos) acta en el centro de la porcin de

varilla sumergida, en la posicin de abscisa
xe=12ycossin
Para que la varilla est en equilibrio, el momento de estas fuerzas respecto del
punto P es cero.
Exe-mgxg=0
gAycosycossin2gALL2sin=0
Simplificando
( y2L21cos2 ) sin=0
Si y>L el primer factor entre parntesis no puede ser cero. La posicin de

equilibrio se obtiene cuando sen=0, es decir, =0.
Cuando y<L el primer factor puede ser cero cuando
cos=yL11
y como el coseno tiene que ser menor o igual que la unidad se tiene que cumplir a
la vez que yL Cuando no se cumple esta condicin el primer factor no es nulo y
la posicin de equilibrio es =0.
Equilibrio de la varilla
Vamos a representar la varilla a medida que se cambia la posicin y del punto de
sujecin P.
Los ngulos de equilibrio que forma la varilla con la direccin vertical, son
El punto P de sujecin de la varilla est por encima de la superficie agua

(y>0)
y>L1 =0
yL1 cos=yL11
El punto P de sujecin de la varilla est por debajo de la superficie del
agua (y<0)
|y|>L =
|y|L cos=yL1
%L=1, longitud de la varilla
rho=0.3; %densidad
x_area=[0,1,1,0]; %agua
y_area=[0,0,-1,-1];
fill(x_area,y_area,'c')
%calcula el ngulo, que hace la varilla con la vertical
a=sqrt(1-rho);
b=sqrt(rho);
for y=-0.9:0.1:0.9
if (y>0)
if (y<a)
angulo=acos(y/a);
else
angulo=0;
end
else
if abs(y)<b
angulo=acos(y/b);
else
angulo=pi;
end
end
line([0,sin(angulo)],[y,y-cos(angulo)], 'color','r');
end
axis equal
xlabel('x')
ylabel('y')
title('Barra que flota')
Actividades
Se introduce
La densidad de la varilla, en el control titulado Densidad barra

La posicin y del punto P de sujecin de la varilla, en el control
titulado Posicin
La longitud de la varilla se ha fijado en el programa interactivo en el
valor L=1.0
Se pulsa el botn titulado Nuevo
Observamos, la varilla en su posicin de equilibrio. Se muestran las fuerzas, peso

de la varilla y fuerza de empuje sobre la parte sumergida.
En la parte derecha, se representa el ngulo , (en grados) en funcin de la
posicin y del punto de sujecin. Se representa mediante un punto el estado de
equilibrio, el par de valores (y, ).
Fijamos la densidad de la varilla, y vamos cambiando la posicin y del punto P de

sujecin de la varilla.
Ejemplo:
Introducimos los siguientes datos
Densidad =0.3
Posicin y=-0.4 del punto P de sujecin de la varilla.
Como y<0, comprobamos que se cumple la condicin
|y|L 0.4<1.00.3
El ngulo que forma la varilla con el eje vertical se obtiene de la expresin
Oscilaciones de una esfera

que flota en el agua
Consideremos una esfera de radio R que tiene una densidad <1 y que se mantiene
completamente sumergida en agua. Se suelta la esfera y se observa su movimiento
oscilatorio
En esta pgina, vamos a comprobar que su comportamiento difiere del Movimiento

Armnico Simple (M.A.S.).
Supondremos que el agua y el aire son fluidos ideales, que no ejercen fuerzas de
rozamiento sobre la esfera en movimiento.
Para describir el movimiento, situamos el origen del eje X en la superficie del agua
y llamamos x a la posicin del centro de la esfera
Cuando x=-R, la esfera se encuentra completamente sumergida

Cuando x=+R la esfera se encuentra justamente fuera del agua
Cuando la esfera se encuentra parcialmente sumergida las fuerzas que actan son:
El peso mg
El empuje E
Para una esfera de densidad relativa al agua (cuya densidad es la unidad) la

masa es
m= 4 3 R 3
El empuje es el peso en agua del volumen de la parte sumergida. Calculamos el

volumen de la parte de la esfera sumergida en agua. Este volumen V es la suma
(integral) de los elementos diferenciales de volumen de radio y y de altura dx, uno
de los cuales se muestra en la figura.
V= x R y 2 dx = x R ( R 2 x 2 )dx = R 3 3 ( 23 x R + x 3 R 3 )
Cuando x=-R obtenemos el volumen de la esfera 4R3/3
La ecuacin del movimiento es
m d 2 x d t 2 =Emg d 2 x d t 2 =g( 1 4 ( 23 x R x 3 R 3 )1 )
Para calcular la posicin x del centro de la esfera en funcin del tiempo t,
resolvemos esta ecuacin diferencial por procedimientos numricos con las
siguientes condiciones iniciales, en el instante t=0, x=-R, dx/dt=0. Cuando la
esfera se encuentra completamente sumergida x=-R se suelta (su velocidad inicial
es cero)
Transformamos la ecuacin diferencial de segundo orden, en la ecuacin

diferencial de primer orden
d 2 x d t 2 = dv dt = dv dx dx dt =v dv dx v dv dx =g( 1 4 ( 23 x R x 3 R 3 )
1 )
Integramos esta ecuacin entre x=-R donde la velocidad de la esfera es nula

(posicin inicial) y la posicin xR, donde la velocidad es v.
v 2 =2gR( 1 x R + 1 16 ( 13+8 x R 6 x 2 R 2 + x 4 R 4 ) )
Cuando x=R la esfera se encuentra completamente fuera del agua y su velocidad

es
v 1 2 =2gR( 12 )
a partir de este momento se eleva en el aire hasta una altura mxima y vuelve a
caer en el agua, volviendo a sumergirse hasta alcanzar la posicin inicial de
partida, repitindose de nuevo la oscilacin.
La esfera sale completamente fuera del agua en el caso de que su densidad

sea <0.5. Si la densidad es mayor que esta cantidad, la esfera permanece
parcialmente sumergida durante su movimiento oscilatorio.
Cuando la densidad relativa de la esfera es <0.5.
A partir de la posicin x=R, la esfera sale completamente fuera del agua y

sobre la esfera solamente acta la fuerza constante peso. Las ecuaciones
del movimiento son
v= v 1 +(g)(t t 1 ) x=R+ v 1 (t t 1 )+ 1 2 (g)( t t 1 ) 2
Siendo t1 el instante en el que el centro de la esfera pasa por la

posicin x=R.
Alternativamente, aplicando el principio de conservacin de la energa
1 2 m v 2 +mgx= 1 2 m v 1 2 +mgR v 2 = v 1 22g(xR)
El centro de la esfera se eleva hasta la altura mxima cuando v=0.
x mx =R( 1 )
Por ejemplo:
o Densidad =0.4
o El radio de la esfera R=1
La esfera parte de la posicin x=-1 con velocidad v=0 y llega a la

posicin x=+1 con una velocidad v1.
v 1 2 =29.8( 120.4 0.4
Oscilaciones de una boya

que flota en el agua
Tenemos un boya de forma cilndrica flotando en el mar. Se deja caer un objeto sobre la
boya (por ejemplo, una persona que salta encima). La boya empieza a oscilar.
Determinaremos el periodo de la oscilacin y la ecuacin del M.A.S.
Supongamos una boya de forma cilndrica o paralepipdica de densidad s menor que la del
agua, de seccin S y altura h.
Situacin de equilibrio
En el equilibrio, la boya estar sumergida una altura h1 dada por el principio de
Arqumedes:
peso=empuje
sghS=fgh1S , es decir,
sh=fh1
Supongamos que colocamos un bloque de masa m sobre la boya (por ejemplo,

una persona que salta sobre la boya).
La nueva posicin de equilibrio h2 se deduce del principio de Arqumedes
mg+sghS=fgh2S
Oscilaciones
Al colocar el bloque sobre la boya y soltarlo el sistema bloque-boya comienza a

oscilar. Hallaremos el periodo de las oscilaciones
Calculamos la fuerza neta que acta cuando la boya se ha desplazado x de la
posicin de equilibrio. Como vemos en la figura, si el desplazamiento x es hacia
arriba, la resultante es hacia abajo. La fuerza es de signo contrario al
desplazamiento.
F=empuje-peso=fgS(h2-x)g-(mg+sghS)= -fSxg
La fuerza es proporcional al desplazamiento y de signo contrario a ste. El

sistema describe un M.A.S. cuya frecuencia y periodo hallamos a partir de la
segunda ley de Newton
(m+shS)a=-fSxg
o bien, expresado en forma de ecuacin diferencial del MAS
d2xdt2+fSgm+shSx=0d2xdt2+fSgm+shSx=0
El periodo es, por tanto,
P=2m+shSfSgP=2m+shSfSg
La ecuacin del MAS, solucin de la ecuacin diferencial es
x=Asin(t+ )
v=A cos( t+ )
Las condiciones iniciales determinan la amplitud A y la fase inicial .
El bloque se suelta cuando la boya se ha sumergido h1, al poner el bloque, en la

nueva posicin de equilibrio la boya se sumerge h2.
Luego, en el instante t=0, el desplazamiento de la boya respecto a su nueva

posicin de equilibrio es x=h2-h1 y su velocidad v=0.
h2-h1=Asin
0=A cos
La fase inicial es =/2 y la amplitud A= h2-h1
La ecuacin del M.A.S. es finalmente,
x=(h2-h1)sin( t+ /2)=(h2-h1)cos( t)
Ejemplo:
Sea, s=600 kg/m3, S=0.5 m2, y m=100 kg. La altura de la boya es de 1.0 m, y
est fijada en el programa interactivos
La boya se sumerge en el agua hasta la altura

6001=1000h1, es decir, h1=0.6 m 60 cm
Situamos el bloque de color negro sobre la boya. El sistema boya-bloque empieza

a oscilar.
La nueva posicin de equilibrio h2 se calcula aplicando de nuevo el principio de

Arqumedes
100+60010.5=10000.5h2, es decir, h2=0.8 m 80 cm
Medimos esta altura, parando el movimiento, cuando el sistema oscilante pasa por
la posicin de equilibrio (usar los botones || y >|)
La amplitud de la oscilacin es 0.8-0.6 =0.2 m 20 cm tal como podemos ver en

la representacin grfica posicin-tiempo.
El periodo de las oscilaciones es
P=2100+60010.510000.59.8=1.8sP=2100+60010.5100
00.59.8=1.8s
Medimos el periodo sobre el eje horizontal en la representacin grfica de x en

funcin de t.
Cuando situamos el bloque sobre la boya, el centro de masas deja de estar en el

centro de la boya. La nueva posicin del c.m. relativo al centro de la boya se
calcula mediante la siguiente frmula
xcm=m0.5m+(s1S)xcm=m0.5m+(s1S)
En este ejemplo, xcm=0.125 m, o 12.5 cm por encima del centro de la boya.
El c.m. del sistema boya-bloque oscilar alrededor de la posicin 0.8-0.5-

0.125=0.175 m 17.5 cm por debajo de la superficie del agua.
Actividades
Se introduce
La densidad de la boya s en kg/m3, en el control titulado Densidad

El rea de la base de la boya S en m2, en el control titulado rea.
La altura h=1 m de la boya est fijada en el programa.
La masa m del bloque en kg se introduce en el control titulado Masa
bloque.
Se pulsa el botn titulado Nuevo y se ve la boya flotando sobre en el agua

Se pulsa el botn titulado y el bloque de color negro se coloca encima de la
boya que empieza a oscilar. Supondremos que no hay rozamiento y la amplitud de
la oscilacin se mantiene constante. Se utiliza los botones pausa || y paso a
paso >| para medir la amplitud de la oscilacin y el periodo
Si los datos introducidos hacen que el sistema bloque-boya quede completamente

sumergido durante la oscilacin, se pide al usuario para que cambie los datos.
El ludin o diablillo de
Descartes
En esta pgina, se describe el comportamiento complejo, del denominado ludin o diablillo
de Descartes en base a leyes simples:
La ley de Boyle de los gases perfectos
Esta pgina es una continuacin de la pgina titulada la flotacin de un barco, que estudia
las fuerzas sobre un recipiente que contiene una burbuja de aire en su interior, y de la
pgina titulada las oscilaciones de una boya, que estudia el MAS de una boya que flota en
el mar y sobre la cual salta un marinero en el instante inicial.
El diablillo o ludin es un tubo de ensayo de longitud L, de dimetro interior d1 y

de dimetro exterior d2 que se llena parcialmente de agua hasta una altura L-l0.
Luego, se invierte en un recipiente grande que contiene agua, que est cerrado y
conectado por una parte, a un manmetro para medir la presin del aire encerrado
en su interior y por otra, a una jeringa que nos permite variar la presin P del aire
del recipiente.
Esttica del diablillo
En este apartado, estudiaremos el equilibrio del tubo de ensayo (diablillo) con la
burbuja de aire dentro
Control de la presin
Para variar la presin del aire contenido en el recipiente se emplea una jeringa. Si
el volumen total de aire del recipiente y de la jeringa es V0. Al disminuir en x el
volumen de la jeringa la presin aumenta de P0 a P. Si suponemos
una transformacin isotrmica, tendremos que:
P(V0-x)=P0V0
Si x<<V0 hacemos la siguiente aproximacin
P=P01x/V0=P0(1xV0)1P0(1+xV0)P=P01x/V0=P0(1xV0)1P0(1+xV0)
El incremento de presin P=P-P0 es proporcional a la variacin del volumen x de

la jeringa o bien, al desplazamiento de su mbolo.
Peso del tubo de ensayo
Despus de haber invertido el tubo de ensayo dentro del agua del recipiente, la
longitud de la burbuja de aire que se ha formado en el interior del tubo invertido
es l=x+z
La base del tubo se encuentra a una altura x por encima de la superficie del agua
del recipiente y el agua en el tubo se encuentra a una altura z por debajo de dicha
superficie tal como se indica en la primera figura.
Sobre el tubo actan dos fuerzas: el peso y el empuje
El peso
El peso del tubo es igual a la densidad del vidrio multiplicado por el volumen
del tubo, un recipiente en forma de capa cilndrica de dimetro interior d1 y
exterior d2.
V=4(d22d21)LV=4(d22d12)L
despreciamos el volumen de la base del tubo.
El rea interior A del tubo de vidrio es
A=4d21A=4d12
El peso del tubo de vidrio es el producto de la densidad del vidrio v, por el

volumen Vde la capa cilndrica de vidrio y por la aceleracin de la
gravedad g,
vVg
El empuje
El empuje no admite una expresin nica y es distinta
o Cuando el tubo est parcialmente sumergido x>0

o Cuando el tubo est completamente sumergido x0
Analizaremos cada uno de los casos por separado
En el primer caso, el empuje se compone de la suma de dos trminos:
o El empuje del tubo de vidrio que est parcialmente sumergido una

longitud L-x.
gV(1-x/L)
o El empuje de la burbuja de aire que ha desalojado el agua contenida
en el volumen cilndrico de rea A (el rea de la seccin trasversal
del tubo) y de altura z.
gAz
La resultante es la diferencia entre el empuje y el peso
F=gAz+gV(1xL)vgVF=gAz+gV(1xL)vgV (1)
Siendo =1000 kg/m3 la densidad del agua, g=9.8 m/s2 la aceleracin de la

gravedad, y v la densidad del vidrio, aproximadamente 2300 kg/m3.
Trasformacin isoterma
La ecuacin fundamental de la esttica de fluidos nos permite calcular la

presin del aire en la burbuja, que es la presin existente a una
profundidad z por debajo de la superficie del agua contenida en el recipiente
P+gz
Siendo P la presin del aire en el recipiente.
Si suponemos que en todo momento la temperatura permanece constante,

aplicamos la ecuacin de los gases ideales a la burbuja de aire:
o En la situacin inicial, con el tubo en la posicin derecha, tenemos un

volumen de aire Al0 a la presin atmosfrica P0
o En la situacin final, con el tubo en la posicin invertida dentro del
recipiente, tenemos una burbuja de aire cuyo volumen es A(x+z) a
la presin P+gz.
Se cumplir que
P0l0=(P+gz)(x+z) (2)
Situacin de equilibrio
El tubo permanece en equilibrio cuando el peso se igual al empuje o bien,
cuando F=0.
Conocida la presin P determinamos x y z resolviendo un sistema de dos

ecuaciones (1) y (2) con dos incgnitas.
De la segunda ecuacin despejamos x
x=P0l0P+gzzx=P0l0P+gzz (3)
y la sustituimos en la primera, para despejar z, quedando una ecuacin de
segundo grado, az2+bz+c=0, con
a=2g(VL+A)b=P(VL+A)(v)Vgc=VLP0l0(v)VPa=2g(VL+A)b=P(VL
+A)(v)Vgc=VLP0l0(v)VP
Se calcula la raz positiva de la ecuacin de segundo grado
z=b+b24ac2az=b+b24ac2a
Una vez que se ha calculado z, se determina mediante la ecuacin (3), la posicin

de equilibriox de la parte superior del tubo de ensayo.
Presin crtica.
Si se incrementa la presin hasta un valor lmite P, el tubo de ensayo se va

sumergiendo en el agua hasta que la posicin de la parte superior del tubo est en
el origen xe=0. Si incrementamos un poco ms la presin el tubo se hunde
completamente.
La presin crtica P se determina poniendo x=0, en la ecuacin (1) con F=0

(situacin de equilibrio) y despejando P en la ecuacin (2).
z=VA(v1) P=P0l0zgzz=VA(v1) P=P0l0zgz
El tamao z de la burbuja es independiente de su tamao l0inicial, solamente

depende de la geometra del tubo de ensayo, de la densidad del vidrio v y del
agua
Actividades
Se introduce
El tipo de tubo que se va a usar en la experiencia, en el control

titulado Tubos
Se pulsa en el botn titulado Nuevo
Se modifica la altura l0 inicial de la burbuja de aire en el tubo, moviendo con

el puntero del ratn el crculo de color rojo, situada al lado de la imagen
amplificada del tubo.
Se modifica la presin P=P-P0 en mm de mercurio, del aire contenido en el
recipiente, en el control titulado Presin.
Observaremos el tubo en la posicin de equilibrio estable F=0, si es que existe. En

caso contrario, un mensaje nos lo indica.
Al lado del tubo, se muestran dos vectores que representan el peso y al empuje, y
en la parte superior del recipiente, se muestra el valor de la fuerza F.
Datos de los tubos utilizados, tomados del artculo citado en las referencias
Dimetro Dimetro
Tubo Altura L (cm) Densidad v(g/cm3)
interior d1(cm) exterior d2(cm)
1 9.8 1.38 1.6 2.35
2 16.0 1.64 1.64 2.29
3 17.1 4.44 4.97 2.16
4 20.1 1.74 2.00 2.17
La presin atmosfrica es P0=101300 Pa. La densidad del mercurio (Hg) es 13.55

g/cm3
Ejemplo de presin crtica
Si elegimos Tubo 2.
Con los datos de la tabla anterior calculamos el rea de la seccin trasversal y el

volumen de vidrio del tubo de ensayo A=1.61 cm2, V=8.10 cm3.
Calculamos la presin crtica para los estados iniciales siguientes:
Altura inicial l0 de la burbuja Presin Diferencia de presin P=P-

de aire en cm crtica P (Pa) P0en mm de Hg
4 61632 -299
6 92766 -64
8 123901 170
Para comenzar una nueva experiencia se pulsa en el botn titulado Nuevo.

Cuando el tubo est completamente sumergido x0
Empuje
El empuje se compone de la suma de dos trminos:
o El empuje del tubo de vidrio que ahora est completamente

sumergido
gV
o El empuje de la burbuja de aire que ha desalojado el agua contenida

en el volumen cilndrico de rea A (el rea de la seccin trasversal
del tubo) y de altura z.
gAz
La resultante de las fuerzas que actan sobre el tubo es la diferencia entre

el empuje y el peso
F=gAz+gV- vgV
Trasformacin isoterma
La presin del aire contenido en la burbuja es, para x<0
P+g(z+|x|)=P+g(z-x)
Suponiendo que el aire experimenta una transformacin isotrmica entre el

estado inicial y el final, tendremos
P0l0= (P+g(z-x))z
Dinmica del diablillo

De nuevo, escribimos la fuerza que acta sobre el tubo de ensayo
F={gAz+gV(1x/L)vgV x>0gAz+gVvgV
x0F={gAz+gV(1x/L)vgV x>0gAz+gVvgV x0
Relacionamos z y x mediante la ley de Boyle
P0l0={(P+gz)(z+x) x>0(P+g(zx))z x0P0l0={(P+gz)(z+x)

x>0(P+g(zx))z x0
Tenemos que despejar z de las ecuaciones de segundo grado
gz2+(P+gx)z+PxP0l0=0 x>0gz2+(Pgx)zP0l0=0
x0z=b+b24ac2aa=g b={P+gxPgx c={PxP0l0 x>0P0l0
x0gz2+(P+gx)z+PxP0l0=0 x>0gz2+(Pgx)zP0l0=0
x0z=b+b24ac2aa=g b={P+gxPgx c={PxP0l0 x>0P0l0
x0
La fuerza F es ahora, una funcin de la posicin x de la base del tubo de ensayo

invertido
Para x>0
F=gV(1xL)vgV+A2(P+gx+(gxP)2+4gP0l0)F=
gV(1xL)vgV+A2(P+gx+(gxP)2+4gP0l0)
La fuerza F depende de la posicin x, se trata de una fuerza conservativa cuya

energa potencial se calcula del siguiente modo
xexFdx=Ep(xe)Ep(x)xexFdx=Ep(xe)Ep(x)
Donde xe es la posicin de equilibrio al cual le vamos a asignar una energa

potencial cero. Establecemos por tanto, en esta posicin el nivel cero de energa
potencial Ep(xe)=0.
Para x>0 tenemos que hacer la integral
E1(x)=F(x)dx+C1E1(x)=g2(A2+VL)x2+(AP2+V(v))xA2f(x)+C1E1(x)=F(x)dx
+C1E1(x)=g2(A2+VL)x2+(AP2+V(v))xA2f(x)+C1
donde C1 es una constante de integracin que se determina a partir de la condicin

de que E1(xe)=0, f(x) es una funcin que se define ms abajo.
Para x0
La fuerza es
F=(v)gV+A2(Pgx+(gxP)2+4gP0l0)F=(v)
gV+A2(Pgx+(gxP)2+4gP0l0)
Para calcular la energa potencial, tenemos que hacer la integral
E2(x)=F(x)dx+C2E2(x)=Ag4x2+(AP2+V(v))xA2f(x)+C2E2(x)=F(x)dx+C2
E2(x)=Ag4x2+(AP2+V(v))xA2f(x)+C2
La fuerza F es una funcin discontinua en x=0, pero la energa potencial Ep(x) es

una funcin continua de modo que, en x=0 se tiene que cumplir que E1(0)=E2(0),
esta condicin determina el valor de la constante C2 de integracin.
La funcin f(x) es una integral que no es inmediata
f(x)=(gxP)2+4gP0l0dxf(x)=(gxP)2+4gP0l0dx
Esta integral se puede escribir de la forma
4P0l0t2+1dt t=(gxP4gP0l0)4P0l0t2+1dt
t=(gxP4gP0l0)
Se descompone en la suma de dos integrales
t2+1dt=t2+1t2+1dt=ttt2+1dt+dtt2+1=tt2+1
+ln(t+t2+1)2t2+1dt=t2+1t2+1dt=ttt2+1dt+dtt2+1=tt2+1+ln(t
+t2+1)2
El primer trmino se integra por partes, y el segundo es una integral inmediata. El

resultado final es
f(x)=(x2P2g)(gxP)2+4gP0l0+2P0l0ln(gxP+(gxP)2
+4gP0l0)f(x)=(x2P2g)(gxP)2+4gP0l0+2P0l0ln(
gxP+(gxP)2+4gP0l0)
Una de las posibles representaciones grficas de Ep(x) se muestra en la figura.

Para presiones bajas del aire contenido en el recipiente, la funcin potencial
tiene un mnimo local en xe>0, que es una posicin de equilibrio estable, y
un mximo local en una posicin xm<0 (equilibrio inestable)
Cuando la presin P se hace igual a la crtica, el mximo y el mnimo
coinciden en x=0. Para esta presin y valores superiores de la misma el
tubo de ensayo se hunde en el recipiente.
Si en la posicin de equilibrio, xe le damos al tubo una velocidad inicial v0, la

energa total del mvil es la energa cintica ya que la energa potencial es nula.
E=12mv20=12(vV)v20E=12mv02=12(vV)v02
La energa E total es constante en todos los puntos de la trayectoria. Si E es menor

que el mximo local de la energa potencial, el tubo que sale de la posicin de
equilibrio xe oscila entre dos posiciones, x1 y x2 determinadas por las races de la
ecuacin trascendente Ep(x)-E=0, abscisas de los puntos de interseccin de la
curva de la energa potencial y la recta horizontal Ep(x)=E, tal como se indica en al
figura.
Cuando la funcin energa potencial Ep(x) no tiene mnimo local o posicin de

equilibrio estable, el tubo sale de la posicin x=0 con velocidad v0 y se hunde en el
recipiente hasta que llega al fondo.
Se simula el movimiento del diablillo de Descartes resolviendo numricamente, la

ecuacin diferencial del movimiento, con las condiciones iniciales especificadas.
Actividades
Se introduce
El tipo de tubo de ensayo que se va a usar en la experiencia,

seleccionndolo en el control de seleccin titulado Tubos

Se modifica la altura l0 inicial de la burbuja de aire en el tubo, moviendo con
el puntero del ratn el crculo de color rojo, situada al lado de la imagen
amplificada del tubo.
Se modifica la presin P=P-P0 en mm de mercurio, del aire contenido en el
recipiente, actuando en la barra de desplazamiento, o introduciendo su
valor en el control titulado Presin.
Se introduce el valor de la velocidad inicial v0 del tubo en el control
titulado Velocidaden la posicin inicial xe de equilibrio estable, si existe. En
caso contrario, la posicin inicial es x=0,
Se pulsa en el botn titulado
Observamos el movimiento del tubo, oscilatorio, si la energa total E del tubo es

menor que la del mximo local de la energa potencial. Se representa a la
izquierda, la curva Ep(x), la energa total E, y los puntos x1 y x2 de retorno.
En el caso de que E sea mayor que el mximo local, el tubo se hunde en el

recipiente y llega hasta el fondo. Lo mismo ocurre si la curva de energa potencial
no tiene mnimo ni mximo local, es decir, la presin del aire contenido en el
recipiente es superior a la crtica.
Si no hay posicin de equilibrio estable, modificamos la presin o bien,

modificamos el tamao inicial l0 de la burbuja de aire actuando con el puntero del
ratn en el crculo de color rojo.
Para comenzar una nueva experiencia se pulsa en el botn titulado Nuevo. Se

utiliza los botones pausa || y paso a paso >| para medir la amplitud de la
oscilacin y el periodo
Densidad de un fluido
La densidad de una sustancia se define como el cociente de su masa entre el
volumen que ocupa.
=mV=mV
La unidad de medida en el S.I. de Unidades es kg/m3, tambin se utiliza

frecuentemente la unidad g/cm3
Densidad de slidos y lquidos a (20C)
Densidad Densidad
Sustancia Sustancia
(g/cm3) (g/cm3)
Acero 7.7-7.9 Oro 19.31

Aluminio 2.7 Plata 10.5
Cinc 7.15 Platino 21.46
Cobre 8.93 Plomo 11.35
Cromo 7.15 Silicio 2.3
Estao 7.29 Sodio 0.975
Hierro 7.88 Titanio 4.5
Magnesio 1,76 Vanadio 6.02
Nquel 8.9 Volframio 19.34
Densidad Densidad
Sustancia Sustancia
(g/cm3) (g/cm3)
Aceite 0.8-0.9 Bromo 3.12
Acido sulfrico 1.83 Gasolina 0.68-0.72
Agua 1.0 Glicerina 1.26
Agua de mar 1.01-1.03 Mercurio 13.55
Alcohol etlico 0.79 Tolueno 0.866
Fuente: Manual de Fsica Elemental. Koshkin N. I., Shirkvich M. G.. Edtorial Mir
(1975) (pgs. 36-37).
Concepto de presin
Se define presin como el cociente entre la componente normal de la fuerza sobre

una superficie y el rea de dicha superficie.
p=FnSp=FnS
La unidad de medida recibe el nombre de pascal (Pa).
La fuerza que ejerce un fluido en equilibrio sobre un cuerpo sumergido en

cualquier punto es perpendicular a la superficie del cuerpo. La presin es una
magnitud escalar y es una caracterstica del punto del fluido en equilibrio, que
depender nicamente de sus coordenadas como veremos en la siguiente pgina.
En la figura, se muestran las fuerzas que ejerce un fluido en equilibrio sobre las
paredes del recipiente y sobre un cuerpo sumergido. En todos los casos, la fuerza
es perpendicular a la superficie, su magnitud y el punto de aplicacin se calculan a
partir la ecuacin fundamental de la esttica de fluidos.
Variacin de la presin con la profundidad

Consideremos una porcin de fluido en equilibrio de altura dy y de seccin S,
situada a una distancia y del fondo del recipiente que se toma como origen.
Las fuerzas que mantienen en equilibrio a dicha porcin de fluido son las
siguientes:
El peso, que es igual al producto de la densidad del fluido, por su volumen y

por la intensidad de la gravedad, (Sdy)g.
La fuerza que ejerce el fluido sobre su cara inferior, pS
La fuerza que ejerce el fluido sobre su cara superior, (p+dp)S
La condicin de equilibrio establece que
( Sdy)g+pS=(p+dp)S
dp=-gdy
Integrando esta ecuacin entre los lmites que se indican en la figura
pApBdp=yAyBgdy pBpA=gyAgyBpApBdp=yAyBgdy
pBpA=gyAgyB
Situamos el punto B en la superficie y el punto A a una profundidad h. Si p0 es la

presin en la superficie del fluido (la presin atmosfrica), la presin p a la
profundidad h es
p=p0+ gh
Medida de la presin. Manmetro

Para medir la presin empleamos un dispositivo denominado manmetro. Como A
y B estn a la misma altura la presin en A y en B debe ser la misma. Por una
rama la presin en B es debida al gas encerrado en el recipiente. Por la otra rama
la presin en A es debida a la presin atmosfrica ms la presin debida a la
diferencia de alturas del lquido manomtrico.
p=p0+gh
Experiencia de Torricelli
Para medir la presin atmosfrica, Torricelli emple un tubo largo, cerrado por uno
de sus extremos, lo llen de mercurio y le dio la vuelta sobre una vasija de
mercurio. El mercurio descendi hasta una alturah=0.76 m al nivel del mar. Dado
que el extremo cerrado del tubo se encuentra casi al vaco p=0 y sabiendo la
densidad del mercurio es 13.55 g/cm3 13550 kg/m3 el valor de la presin
atmosfrica es
p=gh=135509.810.76=101023 Pap=gh=135509.810.76=101023 Pa
Actividades
En esta actividad comprobamos la ecuacin fundamental de la esttica de fluidos,
es decir, que la presin vara linealmente con la profundidad.
Se conecta un tubo por un extremo a un manmetro y por el otro a un elemento o

cpsula de presin consistente en un cilindro de metal con un diafragma de goma,
dispuesto para medir la presin hidrosttica. El elemento de presin se introduce
en el fluido a una profundidad h. En la prctica real, el elemento de presin se
puede girar a fin de demostrar que la presin solamente depende de la posicin,
pero es independiente de la direccin en la que se mide.
Seleccionamos uno de los fluidos cuyas densidades se recogen en la tabla y a

continuacin, se pulsa en el botn titulado Nuevo.
Densidad
Sustancia
(kg/m3)
Agua 1000
Aceite 900
Alcohol 790
Glicerina 1260
Mercurio 13550
La ltima sustancia es el lquido manomtrico, el mercurio.
Establecemos la profundidad deseada en el control titulado Profundidad.

Medimos con el manmetro la presin, o tambin en la grfica de la derecha,
donde se representa la profundidad en el eje vertical y la presin en el eje
horizontal.
Ejemplo:
Supongamos que el fluido es agua. Bajamos la cpsula de presin hasta una

profundidad de 60 cm. La presin debida a la altura de fluido es
p=10009.80.6=5880 Pa
El manmetro marca 2.2 cm por ambas ramas, que corresponde a una presin de
p=135509.820.022=5843
Como el manmetro est abierto por el otro extremo, no nos mide la presin total
(atmosfrica ms la altura de fluido) sino solamente la presin debida al fluido.
Como vemos en la grfica de la derecha a la profundidad de 60 cm le

corresponden algo menos de 106 000 Pa, que corresponden a la presin
atmosfrica (aproximadamente 100 000 Pa) ms la presin debida a la altura de la
columna de fluido (6000 Pa).
La grfica de la derecha est trazada de forma no usual, ya que la presin

(variable dependiente) debera estar en el eje vertical y la altura (variable
independiente) en el eje horizontal. La grfica por tanto, nos muestra la
dependencia lineal de la presin p con la profundidad h.
p=p0+gh
Dinmica de fluidos
El movimiento de un fluido real es muy complejo. Para simplificar su descripcin
consideraremos el comportamiento de un fluido ideal cuyas caractersticas son las
siguientes:
1. Fluido no viscoso. Se desprecia la friccin interna entre las distintas partes del
fluido
2. Flujo estacionario. La velocidad del fluido en un punto es constante con el tiempo
3. Fluido incompresible. La densidad del fluido permanece constante con el tiempo
4. Flujo irrotacional. No presenta torbellinos, es decir, no hay momento angular del
fluido respecto de cualquier punto.
Ecuacin de la continuidad
Consideremos una porcin de fluido en color amarillo en la figura, el instante
inicial t y en el instante t+t.
En un intervalo de tiempo t la seccin S1 que limita a la porcin de fluido en la

tubera inferior se mueve hacia la derecha x1=v1t. La masa de fluido desplazada
hacia la derecha es m1=S1x1=S1v1t.
Anlogamente, la seccin S2 que limita a la porcin de fluido considerada en la

tubera superior se mueve hacia la derecha x2=v2t. en el intervalo de tiempo
t. La masa de fluido desplazada es m2=S2v2 t. Debido a que el flujo es
estacionario la masa que atraviesa la seccin S1 en el tiempo t, tiene que ser
igual a la masa que atraviesa la seccin S2 en el mismo intervalo de tiempo. Luego
v1S1=v2S2
Esta relacin se denomina ecuacin de continuidad.
En la figura, el radio del primer tramo de la tubera es el doble que la del segundo
tramo, luego la velocidad del fluido en el segundo tramo es cuatro veces mayor
que en el primero.
Ejemplo:
Cuando se abre poco a poco un grifo, se forma un pequeo chorro de agua, un hilo
cuyo radio va disminuyendo con la distancia al grifo y que al final, se rompe
formando gotas.
La ecuacin de continuidad nos proporciona la forma de la superficie del chorrito

de agua que cae del grifo, tal como apreciamos en la figura.
La seccin trasversal del chorro de agua cuando sale del grifo es S0 y la velocidad
del agua es v0. Debido a la accin de la gravedad la velocidad v del agua se
incrementa. A una distancia h del grifo la velocidad es
v2=v20+2ghv2=v02+2gh
Aplicando la ecuacin de continuidad
S0v0=Svr20v0=r2vS0v0=Svr02v0=r2v
Despejamos el radio r del hilo de agua en funcin de la distancia h al grifo.
r=r0v20v20+2gh4r=r0v02v02+2gh4
Ecuacin de Bernoulli
Evaluemos los cambios energticos que ocurren en la porcin de fluido sealada
en color amarillo, cuando se desplaza a lo largo de la tubera. En la figura, se
seala la situacin inicial y se compara la situacin final despus de un tiempo t.
Durante dicho intervalo de tiempo, la cara posterior S2 se ha desplazado v2 t y la
cara anterior S1 del elemento de fluido se ha desplazado v1t hacia la derecha.
El elemento de masa m se puede expresar como m=S2v2t= S1v1t=V
Comparando la situacin inicial en el instante t y la situacin final en el

instante t+t.Observamos que el elemento m incrementa su altura, desde la
altura y1 a la altura y2
La variacin de energa potencial es
Ep=mgy2-mgy1=V(y2-y1)g
El elemento m cambia su velocidad de v1 a v2,
La variacin de energa cintica es
Ek=12mv2212mv21=12V(v22v21)Ek=12mv2212mv12=12
V(v22v12)
El resto del fluido ejerce fuerzas debidas a la presin sobre la porcin de fluido
considerado, sobre su cara anterior y sobre su cara posterior F1=p1S1 y F2=p2S2.
La fuerza F1 se desplaza x1=v1t. La fuerza y el desplazamiento son del mismo

signo
La fuerza F2 se desplaza x2=v2 t. La fuerza y el desplazamiento son de signos

contrarios.
El trabajo de las fuerzas exteriores es Wext=F1 x1- F2 x2=(p1-p2) V

El teorema del trabajo-energa nos dice que el trabajo de las fuerzas exteriores
que actan sobre un sistema de partculas modifica la energa del sistema de
partculas, es decir, la suma de las variaciones de la energa cintica y la energa
potencial del sistema de partculas
Wext=Ef-Ei=(Ek+Ep)f-(Ek+Ep)i=Ek+Ep
Simplificando el trmino V y reordenando los trminos obtenemos la ecuacin de

Bernoulli
p1+gy1+12v21=p2+gy2+12v22p1+gy1+12v12=p2+gy2+12v22
Efecto Venturi
Cuando el desnivel es cero, la tubera es horizontal. Tenemos entonces, el

denominado tubo de Venturi, cuya aplicacin prctica es la medida de la velocidad
del fluido en una tubera. El manmetro mide la diferencia de presin entre las dos
ramas de la tubera.
La ecuacin de continuidad se escribe
v1S1=v2S2
Que nos dice que la velocidad del fluido en el tramo de la tubera que tiene menor
seccin es mayor que la velocidad del fluido en el tramo que tiene mayor seccin.
Si S1>S2, se concluye que v1<v2.
La ecuacin de Bernoulli con y1=y2
p1+12v21=p2+12v22p1+12v12=p2+12v22
Como la velocidad en el tramo de menor seccin es mayor, la presin en dicho
tramo es menor.
Si v1<v2 se concluye que p1>p2. El lquido manomtrico desciende por el lado

izquierdo y asciende por el derecho
Obtenemos las velocidades v1 y v2 en cada tramo de la tubera a partir de la

lectura de la diferencia de presin p1-p2 en el manmetro.
v2=S12(p1p2)(S21S22)v2=S12(p1p2)(S12S22)
Ejemplo:
Determinar la velocidad del agua en ambos tramos de la tubera, sabiendo que:
Radio del tramo izquierdo de la tubera, 20 cm.

Radio del tramo derecho de la tubera, 5 cm.
Medida de la diferencia de presin, 1275 Pa.
Los datos son:
S1=(0.2)2 m2, S2=(0.05)2 m2, =1000 kg/m3, y p1-p2=1275 Pa.
Empleando la expresin anterior, obtenemos el valor de v2=1.6 m/s.

Calculamos v1 a partir de la ecuacin de continuidad (v1S1=v2S2)
obteniendo v1=0.1 m/s 10 cm/s.
Introducimos en el programa interactivo los siguientes datos y comprobamos el

resultado:
Radio del tramo izquierdo de la tubera, 20 cm.

Velocidad del fluido en el tramo izquierdo, 10 cm/s
Diferencia de alturas entre los dos tramos, 0
Actividades
Se introduce
El radio del tramo izquierdo de la tubera, en el control titulado Radio.

El radio del tramo derecho est fijado en 5 cm.
El valor de la velocidad del tramo izquierdo, en el control
titulado Velocidad.
El desnivel, (un nmero positivo, nulo o negativo) o diferencia de alturas
entre los dos tramos, en el control titulado Desnivel.

El valor de la velocidad en el tramo derecho se obtiene aplicando la ecuacin de
continuidad. Si el radio del tramo izquierdo es el doble que el radio del tramo
derecho, la velocidad en el tramo derecho es cuatro veces mayor que en el
izquierdo. Es decir, mientras que la seccin anterior S1 del elemento de fluido se
desplaza 10 cm, la seccin posterior S2 se desplaza 40.
A continuacin, nos fijaremos en los cambios energticos.
A medida que el elemento de fluido (coloreado de amarillo) se mueve hacia la

derecha su energa cambia. En la parte inferior izquierda, se muestra la variacin
de energa cintica, de energa potencial y el trabajo de las fuerzas exteriores (que
ejerce el resto del fluido sobre el elemento de fluido considerado). Las fuerzas
exteriores se sealan mediante flechas. Comprobamos que la suma de las
variaciones de energa cintica y potencial nos da el trabajo de las fuerzas
exteriores.
Vaciado de un depsito
abierto
Se calcula y se representa la altura h de un fluido en un depsito en funcin del tiempo,
cuando se abre un orificio en el fondo.
Teorema de Torricelli
Un depsito cilndrico, de seccin S1 tiene un orificio muy pequeo en el fondo de
seccin S2 mucho ms pequea que S1. Aplicamos el teorema de Bernoulli a los
puntos (1) y (2) situados en la superficie libre del fluido y en el centro del orificio
inferior.
p1+gy1+12v21=p2+gy2+12v22p1+gy1+12v12=p2+gy2+12v22
suponiendo que la velocidad del fluido en la seccin mayor S1es despreciable v1 0

comparada con la velocidad del fluido v2en la seccin menor S2.
Por otra parte, el elemento de fluido delimitado por las secciones S1 y S2 est en
contacto con el aire a la misma presin. Luego, p1=p2=p0.
La diferencia de alturas es y1-y2=h. Siendo h la altura de la columna de fluido
Con estos datos la ecuacin de Bernoulli se escribe
gh=12v22 v2=2ghgh=12v22 v2=2gh
En la deduccin del teorema de Torricelli hemos supuesto que la velocidad del
fluido en la seccin mayor S1 es despreciable v1 0 comparada con la velocidad del
fluido v2 en la seccin menor S2.
Supondremos ahora, que v1 no es despreciable frente a v2.
La ecuacin de continuidad se escribe
v1S1=v2S2
y la ecuacin de Bernoulli
gh+12v21=12v22gh+12v12=12v22
De estas dos ecuaciones obtenemos v1 y v2
v2=S12ghS21S22v2=S12ghS12S22
Si S1>>S2 obtenemos el resultado de Torricelli
El volumen de fluido que sale del depsito en la unidad de tiempo es S2v2 y en el

tiempo dtser S2v2dt . Como consecuencia disminuir la altura h del depsito
-S1dh= S2v2dt
Si la altura inicial del depsito en el instante t=0 es H. Integrando esta ecuacin

diferencial, obtenemos la expresin de la altura h en funcin del tiempo.
Hhdhh=S22gS21S220tdt2H2h=S22gS21S22tHhdhh=S
22gS12S220tdt2H2h=S22gS12S22t
Cuando h=0, despejamos el tiempo t que tarda el depsito en vaciarse por

completo.
t=(S21S221)2Hg t=(S12S221)2Hg
Si S1>>S2, se puede despreciar la unidad
t=S1S22Hgt=S1S22Hg
Ejemplo.
Radio del depsito 10 cm, luego, S1= (0.1)2 m2

Radio del orificio 0.8 cm, luego, S2=(0.008)2 m2
Altura inicial 45 cm, H=0.45 m
Sustituyendo estos datos en la frmula del tiempo obtenemos t=47.34 s, que es el

tiempo que tarda en vaciarse completamente el depsito. Si aplicamos la
aproximacin S1>>S2, obtenemos prcticamente el mismo tiempo t=47.35 s.
S1=pi*0.1^2; %seccin del depsito

S2=pi*0.008^2; %seccin del orificio
H=0.45; %altura inicial del fluido
tMax=sqrt(2*H*(S1^2/S2^2-1)/9.8); %tiempo que tarda en salir el fluido
t=linspace(0,tMax,100);
h=(sqrt(H)-S2*sqrt(2*9.8/(S1^2-S2^2))*t/2).^2;
plot(t,h)
grid on
xlabel('t(s)')
ylabel('h(m)')
title('Vaciado de un depsito')
tMax = 47.3499
Actividades
Se introduce:
el radio del depsito R1, en el control titulado Radio depsito

el radio del orificio R2 situado en el fondo, en el control titulado Radio del
orificio
la altura inicial H de agua en el depsito, se ha fijado en 45 cm.
El fluido comienza a salir por el orificio, a la vez que se representa grficamente la

altura de la columna de fluido en funcin del tiempo.
cerrado
En la pgina anterior, hemos estudiado el vaciado de un depsito, suponiendo que est
abierto por arriba.
Vamos a estudiar en esta pgina, el vaciado de un depsito de agua que est cerrado por la
parte superior mediante una tapa hermtica y que contiene aire en su interior a una presin
inicial dada.
Este ejemplo nos va a servir de introduccin al estudio del cohete impulsado por agua, un
problema interdisciplinar en el que intervienen, tres partes de la Fsica: Fluidos, Dinmica y
Termodinmica.
A medida que se vaca el depsito, el volumen de aire aumenta y la presin disminuye.

Supondremos que esta disminucin de presin se realiza a temperatura constante, es decir,
se trata de un proceso isotrmico.
En la figura, se muestra un depsito que tiene una altura H y una seccin S1, la
seccin del orificio de salida en el fondo del depsito es S2, la altura inicial de agua
es h0 y la presin del aire en su interior p0.
Se abre el orificio de salida del agua y se mide la altura h de la columna de agua

en funcin del tiempo t.
Aplicamos el teorema de Bernoulli comparando dos puntos del fluido. El punto 1 en

la interfase aire-agua y el punto 2 en el orificio de salida.
Sea p1 la presin del aire en el interior del depsito, v1 la velocidad del agua en el
punto 1 y hla altura de agua en el depsito en el instante t. La presin p2 en el
orificio de salida es la atmosfrica pat y la velocidad del fluido es v2.
Consideramos los puntos (1) y (2) situados en la superficie libre del fluido y en el
centro del orificio inferior. Las ecuaciones que describen el comportamiento de
este sistema fsico son:
1. Ecuacin de continuidad
S1v1=S2v2
2. Ecuacin de Bernoulli
p1+12v21+gh=pat+12v22p1+12v12+gh=pat+12v22
3. Expansin isotrmica del gas
p0S1(H-h0)=p1S1(H-h)
Altura del fluido en equilibrio

La consecuencia ms importante de estas ecuaciones es que el agua deja de salir
por el orificio cuando v2 y por tanto v1 sean nulos.
La presin del aire en el interior del depsito ser algo menor que la presin
atmosfrica. La diferencia ser la presin correspondiente a la columna de agua de
altura h.
Las ecuaciones de Bernoulli con v1=0 y v2=0, y la transformacin isoterma
p1+gh=pat
p0 (H-h0)=p1(H-h)
Obtenemos la ecuacin de segundo grado en h
gh2(gH+pat)hH(p0pat)+p0h0=0gh2(gH+pat)hH(p0pat)+p0h0=0
La raz real positiva se obtiene con p0h0>H(p0-pat)
Los valores de las races no dependen del rea de la seccin del depsito S1, ni del
orificio S2.
Ejemplo:
Sea el radio del depsito, r1=10 cm

El radio del orificio, r2=0.8 cm
La altura del depsito, H=50 cm
La altura inicial de agua en el depsito, h0=40 cm
La presin inicial de aire en el depsito, p0=4 atm
Tomando como presin atmosfrica pat=101293 Pa, la densidad del agua =1000
kg/m3, resolvemos la ecuacin de segundo grado en h y calculamos la altura final
del agua en el depsito.
Las dos races son h1=0.096 m=9.6 cm, y h2=10.74 m que es mayor que H=0.5
m. Cuando la altura de agua en el depsito alcanza 9.6 cm deja de salir por el
orificio. Calculamos la presin final del aire en el depsito
p1=101293-10009.80.09=100411 Pa
un poco menos que la presin atmosfrica
Si cambiamos la presin inicial de aire en el depsitio a p0=6 atm
obtenemos dos races h1=-0.0946 y h2=10.9306. El depsito se vaca

completamente. No se cumple que p0h0>H(p0-pat), 60.4<0.5(6-1)

H=0.5; %altura del depsito
h0=0.4; %altura inicial de agua en el depsito
pAtm=101293; %presin atmosfrica
p0=4*pAtm; %presin inicial del aire en el depsito
rho=1000; %densidad del agua kg/m3
%equilibrio
h1=((rho*9.8*H+pAtm)-
sqrt((rho*9.8*H+pAtm)^2-4*rho*9.8*(p0*h0-H*(p0-pAtm))))/(2*rho*9.8)
h2=((rho*9.8*H+pAtm)+
sqrt((rho*9.8*H+pAtm)^2-4*rho*9.8*(p0*h0-H*(p0-pAtm))))/(2*rho*9.8)
h1 = 0.0962
h2 = 10.7398
Variacin de la altura de agua en el depsito con el
tiempo
Despejamos v1 en el sistema de tres ecuaciones
v1=p0(Hh0)Hh+ghpat12(S21S221) v1=p0(Hh0
)Hh+ghpat12(S12S221)
Para hallar como cambia la altura h del agua en el depsito con el tiempo,
tenemos en cuenta que,
v1=dhdtv1=dhdt
Integramos
0tdt=12(S21S221)h0hdhp0(Hh0)Hh+ghpat
t=12(S21S221)Hp0(Hh0) Hh0/Hh/Hdx11x+gH2p0(Hh0)xpatHp0(Hh0)
t=ch0/Hh/Hdx11x+axb0tdt=12(S
12S221)h0hdhp0(Hh0)Hh+ghpatt=12(S12S221)Hp0(Hh0)Hh0/H
h/Hdx11x+gH2p0(Hh0)xpatHp0(Hh0)t=ch0/Hh/Hdx11x+axb
Consideremos un depsito cerrado de 50 cm de altura y 10 cm de radio, con un

orificio de 8 cm de radio. La altura inicial de agua en el depsito es de 40 cm. La
presin inicial del aire en su interior es de 4 atm que vamos a ir cambiando.
Calculamos el tiempo que tarda en alcanzarse la altura de equilibrio

hFin=0; %altura final del agua en el depsito
%equilibrio
if p0*h0>H*(p0-pAtm)
hFin=((rho*9.8*H+pAtm)-
sqrt((rho*9.8*H+pAtm)^2-4*rho*9.8*(p0*h0-H*(p0-pAtm))))/(2*rho*9.8);
end
c=-sqrt((S1^2/S2^2-1)*H*rho/(2*p0*(H-h0)))*H;
a=rho*9.8*H^2/(p0*(H-h0));
b=pAtm*H/(p0*(H-h0));
f=@(x) 1./sqrt(1./(1-x)+a*x-b);
t=c*integral(f,h0/H,hFin/H);
fprintf('La altura final es %1.3f y el tiempo %1.2f\n',hFin,t);
La altura final es 0.096 y el tiempo 6.50
Si cambiamos la presin inicial de aire en el depsitio a p0=6 atm
La altura final es 0.000 y el tiempo 4.77
Como vemos en el cdigo, se integra numricamente para obtener el tiempo t
Finalmente, modificamos el script para representar la altura del fluido en funcin

del tiempo

hFin=0;
%equilibrio
hFin=((rho*9.8*H+pAtm)-sqrt((rho*9.8*H+pAtm)^2-
4*rho*9.8*(p0*h0-H*(p0-pAtm))))/(2*rho*9.8);
end
c=-sqrt((S1^2/S2^2-1)*H*rho/(2*p0*(H-h0)))*H;
a=rho*9.8*H^2/(p0*(H-h0));
f=@(x) 1./sqrt(1./(1-x)+a*x-b);
t=zeros(50,1);
h=hFin:(h0-hFin)/50:h0;
i=1;
for i=1:length(h)
t(i)=c*integral(f,h0/H,h(i)/H);
end
plot(t,h)
grid on
xlabel('t(s)');
ylabel('x(m)');
Ecuacin diferencial
Resolvemos numricamente la ecuacin diferencial de primer orden
dhdt=p0(Hh0)Hh+ghpat12(S21S221) dhdt=p0
(Hh0)Hh+ghpat12(S12S221)
o de forma equivalente
dxdt=c11x+axb c=2p0(Hh0)(S21S221)H3
a=2gHp0(Hh0) b=patHp0(Hh0)dxdt=c11x+axb c=2p0(Hh0)(S12S221)H3
a=2gHp0(Hh0) b=patHp0(Hh0)
donde x=h/H, en el instante inicial t=0, x=h0/H
El proceso de integracin cesa cuando la altura final de fluido es cero o bien, si la

presin no es suficiente la latura de equilibrio h1, una de las races de la ecuacin
de segundo grado en h.
Observaremos en el cdigo MATLAB que ya no utilizamos la rutina habitual ode45,

que aqu no produce buenos resultados, sino ode15s. Bsquese en Internet el
trmino MATLAB ode, la lista de funciones para resolver ecuaciones diferenciales.

%equilibrio
hFin=((rho*9.8*H+pAtm)-sqrt((rho*9.8*H+pAtm)^2
-4*rho*9.8*(p0*h0-H*(p0-pAtm))))/(2*rho*9.8);
end
c=sqrt(2*p0*(H-h0)/((S1^2/S2^2-1)*H*rho))/H;
a=rho*9.8*H^2/(p0*(H-h0));
f=@(t,x) -c*sqrt(1/(1-x)+a*x-b);
tspan=[0 10];
opts=odeset('events',@(t,x) vaciado_ode45(t,x,hFin/H));
[t,x]=ode15s(f,tspan,h0/H,opts);
plot(t,real(x*H));
grid on
xlabel('t(s)')
ylabel('h(cm)')
La funcin vaciado_ode45, termina el proceso de integracin cuando se alcanza la
altura final de equilibrio x=hfin/H
function [value,isterminal,direction]=vaciado_ode45(~,x,xFin)
value=x-xFin;
isterminal=1;
direction=0;
end
Cuando se eleva la presin inicial del aire en el depsito a 5 atm, el agua sale
completamente del depsito, tal como se muestra en la figura
Se pueden cambiar varios parmetos: altura inicial de agua en el depsito, presin
inicial del aire, la seccin del orificio y ver su efecto, mediante el script de MATLAB
o en la simulacin ms abajo.
Forma alternativa de resolver la ecuacin diferencial
Definimos la funcin v_deposito que calcula el cuadrado de la velocidad de salida

del agua del depsito por el orificio inferior
function v2 = v_deposito(t,x,p)
r1=p(1);
r2=p(2);
h0=p(3)/100;
p0=p(4);
pAtm=101293; %en Pa
H=0.5; %cm, altura depsito
den=(r1^4/r2^4-1)*500; %1000 es la densidad del agua
v2=(p0*pAtm*(H-h0)/(H-x)+1000*9.8*x-pAtm)/den;
end
Definimos la funcin events_1que detiene el proceso de integracin cuando:
la velocidad de salida del agua por el orificio inferior sea cero.

cuando la altura del agua en el depsito sea cero
function [detect,stopin,direction] = events_1(t,x,p)

detect = [v_deposito(t,x,p) x];
stopin = [1 1];
direction = [0 0];
end
Creamos el script
r1=10; %radio del depsito en cm

r2=0.8; %radio del orificio en cm
h0=40; %altura inicial de agua en el depsito en cm
p0=4; %presin inicial del aire en el depsito en atm
p=[r1 r2 h0 p0];
f=@(t,x,p) -sqrt(v_deposito(t,x,p));
options = odeset('Events',@events_1);
[t,x,te,xe,ie] = ode15s(f, [0 20],h0/100,options,p);
plot(t,x)
ylim([0 0.5]);
grid on
xlabel('t')
ylabel('x');
title('Vaciado de un depsito cerrado')
La variable te nos proporciona el tiempo hasta que se detiene el proceso de
integracin y xe el valor de x en dicho instante. El tiempo de vaciado del depsito
y la altura del agua en el depsito son, respectivamente
te = 6.4578
xe= 0.0962
Cambiamos la presin inicial p0=5 atm del aire en el depsito cerrado. El agua sale
completamente del depsito
te = 8.1176
xe= 1.0e-06 * 0.2018 - 0.0000i
Obtenemos dos grficas idnticas a las anteriores
Actividades
El programa interactivo estudia el comportamiento de un depsito de 50 cm de
altura cerrado con un orificio en la parte inferior.
Se introduce:
el radio del depsito r1, en el control titulado Radio depsito

el radio del orificio r2 situado en el fondo, en el control titulado Radio
orificio
la presin inicial p0 de aire (en atm) en el depsito cerrado, en el control
titulado Presin aire
la altura inicial de agua en el depsito h0 en el control titulado Altura del
agua.

Se representa en la parte derecha, la altura h de agua en el depsito en funcin
del tiempo t.
El caso ms importante, ocurre cuando el agua deja de salir por el orificio, se

cumple que la presin del aire en el interior del depsito se hace menor que la
presin atmosfrica, es decir, la diferencia de presin del aire en el interior y en el
exterior del depsito se hace igual a la presin que ejerce la columna de agua de
altura h.
p1-pat=gh
Como hemos visto, esta altura no depende de los radios del depsito r1 ni del
orifico r2.
En la parte superior derecha, se muestra los valores de la presin atmosfrica 101

293 Pa y la presin p1, a medida que sale el agua por el orificio inferior.
Cohete propulsado por agua

Se continua en esta pgina, el estudio de un recipiente cerrado lleno parcialmente de
aguaque contiene aire en su interior a una presin elevada. Cuando se abre el orificio en la
parte inferior del recipiente, el agua expulsada ejerce una fuerza sobre el recipiente similar
al empuje que experimenta un cohete al expulsar el combustible quemado por sus toberas.
El anlisis del sistema fsico tiene las siguientes partes:
Llenado de aire del recipiente mediante una bomba de bicicleta o similar

Apertura del orificio en la parte inferior del recipiente y expulsin del agua, que ya
hemos estudiado en la pgina anterior.
Empuje que experimenta el recipiente al expulsar el agua
Ecuaciones del movimiento.
El cohete consta de un recipiente de forma cilndrica de 10 cm de radio y 50 cm de

altura.
El radio del orificio situado en la parte inferior se puede modificar entre 1/2 y 1/10
del radio del recipiente. Por ejemplo, al elegir , el radio del orificio es 10/4=2.5
cm
Otro dato es la proporcin de agua en el recipiente. Por ejemplo, una proporcin

del 70% equivale a una altura de agua de 0.750=35 cm.
El cohete puede transportar una carga que es la suma de la carga til ms la masa
de las paredes del recipiente.
Se introduce aire comprimido en el cohete con una bomba de volumen Vb=5 litros.
Llenado de aire
Antes de accionar la bomba tenemos n0 moles de aire en el recipiente a la presin

atmosfrica y a la temperatura ambiente T.
patS1(H-h0)=n0RT
Cada vez que accionamos la bomba de volumen Vb, introducimos en el

recipiente n moles de aire a la misma temperatura T.
patVb=nRT
Si accionamos la bomba N veces, tendremos que la presin p0 del aire contenido

en en el recipiente es
p0S1(H-h0)=(n0+nN)RT
El manmetro marcar una presin final p0 dada por la expresin
p0=pat ( )
1+NVbS1(Hh0)
Ejemplo:
Supongamos que el tanto por ciento de agua en el recipiente es del 70%, la altura
inicial de agua es h0=0.7H=0.750=35 cm.
Sabiendo que el volumen de la bomba Vb=5 litros y el recipiente tiene un

radio r1=10 cm. Si accionamos la bomba N=4 veces, la presin del aire en el
recipiente cerrado ser de p0=5.24 atm que es lo que marca el manmetro.
Empuje que experimenta el cohete

El recipiente experimenta un empuje que es el producto de la velocidad de salida
del agua ve(medida en el sistema de referencia del cohete) por la masa de agua
expulsada en la unidad de tiempo dM/dt. La velocidad de salida del agua es v2 y el
volumen de agua expulsada en la unidad de tiempo (gasto) es S2v2.
| |
E= vedMdt =v2S2v2=S2v22
Como hemos visto en la pgina anterior el las ecuaciones que describen este
sistema son:
1. La ecuacin de Bernoulli,
p1+12v21+gh=pat+12v22
2. La ecuacin de continuidad,
S1v1=S2v2
3. Expansin isotrmica del gas
p0S1(H-h0)=p1S1(H-h)
que nos permiten obtener la expresin de v1 v2 en funcin de la altura h de agua

en el recipiente.
v22=p0(Hh0)Hh+ghpat12 ( )
1S22S21
Variacin de la altura del agua en el recipiente con el

tiempo
A partir de la ecuacin de continuidad, obtenemos la variacin de la altura h del
agua en recipiente en funcin del tiempo t.
v1=dhdt dhdt=S2S1v2
Ecuaciones del movimiento

El movimiento del cohete se divide en dos etapas
1. Mientras sale agua por el orificio
La masa del recipiente no es constante, sino disminuye con el tiempo. La

masa del recipiente es la suma de la carga til, de la masa de las paredes
del recipiente y del agua que contiene en el instante t.
m=mu+S1h
La ecuacin del movimiento vertical de un cohete, es la de una partcula de

masa m bajo la accin de dos fuerzas el empuje y el peso.
ma=E-mg
En forma de ecuacin diferencial
d2xdt2=Emg
Tenemos que resolver un sistema de dos ecuaciones diferenciales

simultneas:
o Una ecuacin diferencial de primer orden, que nos calcula la
variacin de h con el tiempo.
o La ecuacin del movimiento. El empuje E y la masa m del cohete son
funciones de h (altura de agua en el recipiente).
( )
v22=f(h)=p0(Hh0)Hh+ghpat12 1S22S21 dhdt=S2S1f(h)d2xdt2=S2S1h+muf(h)g
En el programa interactivo, se ha resuelto el sistema de dos ecuaciones

diferenciales por el mtodo de Runge-Kutta, sin realizar ninguna
aproximacin. Lo que nos permite incluso examinar el caso de que la
presin del aire en el interior del recipiente no sea suficiente para expulsar
todo el agua y se alcance una altura del fluido en equilibrio, tal como vimos
en la pgina anterior.
2. Cuando se ha agotado el agua
Una vez que se ha agotado el agua del depsito, el aire en el interior del
depsito tiene una presin p mayor que la presin atmosfrica, pero
supondremos despreciable el impulso adicional proporcionado por la salida
del aire por el orificio inferior hasta que se igualan las presiones en el
interior y exterior del recipiente. Sobre el cohete acta solamente el peso,
por lo que el movimiento es uniformemente acelerado
( )
a=gv=v0g(tt0)x=x0+v0(tt0)12g tt0 2
donde x0, y v0 son la posicin y la velocidad del mvil en el instante t0 en el

que se ha agotado el combustible, en este caso, agua.
El rozamiento del aire

Al moverse un cuerpo en el aire con velocidad v, experimenta una fuerza de
rozamiento, que es proporcional al cuadrado de la velocidad
Esta fuerza de rozamiento no es importante durante la fase de lanzamiento que

dura poco tiempo y durante la cual predomina la fuerza de empuje, pero puede ser
importante en la fase de vuelo libre desde que se agota el combustible hasta que
alcanza la mxima altura.
La fuerza de rozamiento no se ha tenido en cuenta en la simulacin del cohete
propulsado por agua.
Resultados
En este apartado vamos a examinar la dinmica del cohete de 10 cm de radio, que
transporta una carga til de 5 kg, el orificio de salida del agua es de 10/6 cm y la
proporcin de agua en el depsito es del 70%. Mientras el agua sale por el orificio,
hay que resolver numricamente un sistema de dos ecuaciones diferenciales, o tes
ecuaciones diferenciales de primer orden. Representamos la velocidad del cohete y
su altura sobre el suelo, para dos presiones iniciales de aire distintas: 3 atm y 5
atm. En el primer caso, la presin no es suficiente para que salga todo el agua y
en el segundo, sale completamente expulsada del depsito del cohete.

S2=pi*(0.1/6)^2; %seccin del orificio
h0=0.7*0.5; %altura inicial de agua en el depsito (70%)
mu=5; %carga til que transporta
c=2*p0*(H-h0)/((1-S2^2/S1^2)*rho);
a=rho*9.8/(p0*(H-h0));
b=pAtm/(p0*(H-h0));
f=@(h) c*(1./(H-h)+a*h-b);
if(rho*S2*f(h0)/(rho*S1*h0+mu)<9.8)
disp('No despega, el empuje es menor que el peso')
return;
end

hFin=((rho*9.8*H+pAtm)-sqrt((rho*9.8*H+pAtm)^2-
4*rho*9.8*(p0*h0-H*(p0-pAtm))))/(2*rho*9.8);
end
%x(1) es la altura de agua en el depsito, x(2) la posicin del cohete y
%x(3)su velocidad
x0=[h0,0,0];
fg=@(t,x) [-S2*sqrt(f(x(1)))/S1; x(3); rho*S2*f(x(1))/(rho*S1*x(1)+mu)-
9.8];
tspan=[0 100];
opts=odeset('events',@(t,x) vaciado_cohete(t,x,hFin));
[t,x]=ode15s(fg,tspan,x0,opts);
subplot(2,1,1)
plot(t,x(:,2))
grid on
xlabel('t(s)')
ylabel('x(m)')
subplot(2,1,2)
plot(t,x(:,3))
grid on
xlabel('t(s)')
ylabel('v(m/s)')
function [value,isterminal,direction]=vaciado_cohete(~,x,hFin)
value=[x(1)-hFin, x(2), x(3)];
isterminal=[1,1,0];
direction=0;
end
Con una presin de 3 atm, el cohete acelera, el empuje es mayor que el peso, la
velocidad alcanza un mximo cuando el empuje se iguala al peso y disminuye
cuando el empuje es menor que el peso. El cohete regresa al suelo antes de que
termine de salir el agua por el orificio.
Con una presin de 5 atm, el cohete acelera, el empuje es siempre mayor que el
peso, la velocidad se incrementa
Con el siguiente script investigamos la proporcin ptima de agua en el depsito
para una presin de 5 atm, una carga til de 5 kg y un orificio de radio 10/6 cm.
Vemos que la proporcin ptima es del 46%

S2=pi*(0.1/6)^2; %seccin del orificio
h0=0.7*0.5; %altura inicial de agua en el depsito (70%)
mu=5; %carga til que transporta
i=0;
vMax=zeros(0,81);
hAgua=zeros(0,81);
for percent=0.1:0.01:0.9
i=i+1;
h0=percent*H;
c=2*p0*(H-h0)/((1-S2^2/S1^2)*rho);
a=rho*9.8/(p0*(H-h0));
b=pAtm/(p0*(H-h0));
f=@(h) c*(1./(H-h)+a*h-b);
if(rho*S2*f(h0)/(rho*S1*h0+mu)<9.8)
disp('No despega, el empuje es menor que el peso')
hAgua(i)=percent*100;
vMax(i)=0;
continue;
end
hFin=((rho*9.8*H+pAtm)-sqrt((rho*9.8*H+pAtm)^2
-4*rho*9.8*(p0*h0-H*(p0-pAtm))))/(2*rho*9.8);
end
x0=[h0,0,0];
fg=@(t,x) [-S2*sqrt(f(x(1)))/S1; x(3); rho*S2*f(x(1))
/(rho*S1*x(1)+mu)-9.8];
tspan=[0 100];
opts=odeset('events',@(t,x) vaciado_cohete(t,x,hFin));
[t,x]=ode15s(fg,tspan,x0,opts);
hAgua(i)=percent*100;
vMax(i)=x(end,3);
end
plot(hAgua,vMax);
grid on
xlabel('f(%)')
ylabel('v_m(m/2)')
title('Cohete de agua')
>> [xmax nmax]=max(vMax);
>> xmax
xmax = 16.2649
>> hAgua(nmax)
ans = 46
Con un orificio un poco mayor de radio 10/4 de cm, mejoramos la velocidad final
del cohete con una proporcin del 50% de agua en el depsito
>> [xmax nmax]=max(vMax);

>> vMax(nmax)
ans = 18.4223
>> hAgua(nmax)
ans = 50
Probaremos el efecto de los distintos parmetros, modificando sus valores
Actividades
Se introduce
La proporcin de agua en el recipiente, en el control titulado % de agua

El radio del orificio de salida del agua, eligiendo una fraccin del radio del
depsito, en el control titulado Radio orificio.
La carga que transporta el cohete (que incluye la carga til y las paredes
del recipiente) introduciendo un valor en el control titulado Carga.

A continuacin, se pulsa varias veces en el botn titulado Aire, para introducir aire
en el interior del recipiente. En el manmetro situado a la derecha del cohete
vemos como va aumentando la presin.
Se pulsa el botn titulado .
Veremos que el agua es expulsada por el orificio inferior, el aire se expande en el

interior del recipiente disminuyendo la presin.
Dos flechas indican las magnitudes relativas del peso y del empuje durante la fase
de lanzamiento. Cuando se se termina de salir el agua, el empuje se hace cero y el
cohete se mueve con movimiento uniformemente acelerado hasta que alcanza la
altura mxima.
En la parte izquierda, un punto de color rojo indica la posicin (altura) del cohete.
Los datos del tiempo, velocidad y posicin que aparecen en color rojo,
corresponden al instante en el que deja de salir agua por el orificio y por tanto, el
empuje es nulo. El cohete alcanza la mxima velocidad, contina su movimiento
ascendente hasta que su velocidad se hace cero, alcanzando en ese instante la
mxima altura.
Oscilaciones de un lquido
en un tubo en forma de U.
Deducimos la frmula del periodo de las oscilaciones de un lquido en un tubo en forma de
U, a partir de las ecuaciones de la dinmica para dos casos:.
El tubo en forma de U tiene los extremos abiertos

El tubo en forma de U est cerrado por uno de sus extremos
que nos permite disear una experiencia de medida del ndice adiabtico del aire
El tubo en forma de U tiene los extremos

abiertos.
Consideremos un lquido contenido en un tubo en forma de U de seccin
uniforme S. La longitud del lquido es L y su volumen SL. En la situacin de
equilibrio la altura del lquido en ambas ramas es la misma. Supongamos que por
algn procedimiento, se desnivela el lquido en las dos ramas.
Si consideramos el lquido como un fluido ideal, libre de rozamiento con las
paredes del recipiente que lo contiene, el lquido empezar a oscilar con un periodo
que depende nicamente de la longitud L de de la columna de lquido.
Cuando el lquido se desplaza x hacia la derecha, el desnivel entre las dos ramas
del tubo en forma de U es 2x, tal como se aprecia en la figura. La fuerza F que se
opone al movimiento de todo el fluido, es el peso de la columna de fluido de
seccin S y altura 2x.
La masa de todo el fluido es SL

La fuerza F tiene por mdulo, S(2x)g, es de sentido contrario al
desplazamiento x.
La ecuacin del movimiento se escribe (masa por aceleracin igual a fuerza)
SLd2xdt2=S(2x)g d2xdt2+2gLx=0SLd2xdt2=S(2x)g
d2xdt2+2gLx=0
Se trata de la ecuacin diferencial de un MAS cuya frecuencia angular

es 2=2g/L o cuyo periodo es
P=2L2gP=2L2g
La solucin de la ecuacin diferencial es
x=Asin(t+).
La velocidad de la columna de lquido en el instante t es
v=dx/dt= Acos(t+).
donde A y se determinan a partir de las condiciones iniciales.
Para conseguir un desnivel inicial entre las dos ramas del tubo, se pone un corcho
en el extremo de una de las ramas, se extrae algo de aire, el lquido en esta rama
se eleva h0 y desciende la misma longitud en la otra rama. Se descorcha la rama y
se empieza a contar el tiempo, en el instante inicial t=0, tenemos x=h0 y v=0.
Con estas condiciones iniciales =/2, x=h0cos(t), v=-h0sin(t)
Balance energtico
Para efectuar el balance energtico, comparamos la situacin inicial en el
instante t=0 y la situacin final en el instante t, tal como se muestra en la figura
En el instante t, cuando la columna de lquido de longitud L y seccin S lleva una

velocidad v, la energa cintica es
Ek=12(SL)v2=12(SL)h202sin2(t)=Sgh20sin2(t)Ek=12(SL)v2=12(SL)h02
2sin2(t)=Sgh02sin2(t)
Comparamos la energa potencial en la situacin de equilibrio con la energa

potencial cuando el desnivel del lquido en las dos ramas es 2x.
El lquido contenido en la rama horizontal no cambia su energa potencial. Los

centros de masas de las columnas verticales de fluido se han marcado mediante
un punto de color rojo, cuya altura es la mitad de la longitud de dicha columna.
En la situacin de equilibrio, la energa potencial, referida al suelo, vale
Sdgd/2+Sdgd/2=Sgd2
En el instante t, el desnivel del lquido en las dos ramas es 2x, la energa potencial
referida al suelo es.
S(d-x)g(d-x)/2+ S(d+x)g(d+x)/2=Sg(d2+x2)
La energa potencial con respecto a la situacin de equilibrio estable, es la

diferencia entre ambas cantidades.
Ep=Sgx2=Sgh20cos2(t)Ep=Sgx2=Sgh02cos2(t)
La suma de la energa potencial y de la energa cintica es constante e igual a la

energa potencial en la situacin inicial, t=0, x=h0, v=0
E=Ek+Ep=Sgh20E=Ek+Ep=Sgh02
El tubo en forma de U est cerrado por uno

de sus extremos
Consideremos un tubo en forma de U con uno de sus extremos taponado con un
corcho. Inicialmente, la altura del aire encerrado entre la superficie del lquido y el
tapn es H y est a la presin atmosfrica p0.
Cuando el fluido se desplaza x de la posicin de equilibrio, hacia la derecha, actan

tres fuerzas sobre todo el lquido (SL) contenido en el tubo en forma de U.
El peso de la columna de lquido de altura 2x, opuesta al desplzamiento del

lquido
La fuerza que ejerce la presin del aire comprimido sobre la superficie libre
del lquido en la rama derecha, opuesta al movimiento del fluido
La fuerza que ejerce la presin de aire exterior sobre la superficie libre del
lquido en la rama izquierda, favoreciendo el movimiento.
La fuerza neta que se opone al movimiento del todo el lquido del tubo es
-S(2x)g-pS+p0S
En el apartado anterior, al estar abiertas las dos ramas del tubo, p coincide con p0,
por lo que las fuerzas que ejerca el aire sobre las dos superficies libres del lquido
se cancelaba.
Si suponemos que las oscilaciones son rpidas, el aire encerrado entre el lquido y
el tapn, experimenta una transformacin aproximadamente adiabtica.
p0(SH)=p(S(Hx))p0(SH)=p(S(Hx))
Supondremos tambin que la amplitud de las oscilaciones son pequeas respecto

de la longitud de la columna de aire, H, por lo que x<<H.
p=p0(1xH)p0(1+xH)p=p0(1xH)p0(1+xH)
La ecuacin del movimiento del lquido contenido en el tubo en forma de U ser
(SL)d2xdt2=S(2x)gp0SxHd2xdt2+(2gL+p0HL)x=0(SL)d2xdt2=S(2x)gp0Sx
Hd2xdt2+(2gL+p0HL)x=0
El lquido contenido en el tubo, describe un M.A.S. cuyo periodo es
P=2L2g+p0H P=2L2g+p0H
Si elegimos como condiciones iniciales t=0, x=0, y v=v0. La solucin de la

ecuacin diferencial es
x=(v0/)sin(t)
v=dx/dt=v0cos(t)
Balance energtico
Para efectuar el balance energtico, comparamos la situacin inicial en el
instante t=0 y la situacin final en el instante t, tal como se muestra en la figura
En la situacin inicial t=0, la energa es cintica

Ek=12(SL)v20Ek=12(SL)v02
La energa potencial respecto de la situacin de equilibrio estable es cero
Ep=0
En el instante t, cuando la columna de lquido de longitud L y seccin S lleva una

velocidad v, la energa cintica es
Ek=12(SL)v20cos2(t)Ek=12(SL)v02cos2(t)
La expresin de la energa potencial con respecto a la situacin de equilibrio

estable es la misma que la que se dedujo en el apartado anterior.
Ep=Sgx2=Sgv202sin2(t)Ep=Sgx2=Sgv022sin2(t)
La suma de la energa cintica y de la energa potencial no es una cantidad

constante, ya que hemos de considerar el trabajo realizado por las fuerzas
exteriores (el aire que rodea a la columna de lquido). El aire encerrado en la rama
derecha se comprime adiabticamente y el aire exterior se expande a presin
constante durante media oscilacin y durante la otra media, el aire encerrado se
expande y el aire de la atmsfera que rodea al dispositivo se comprime.
La variacin de energa, es la energa final menos la inicial
E=EfEi=12SLv20cos2(t)+Sgv202sin2(t)12SLv20==Sv20(12L+g2)sin2(t)
=p0S2Hv202sin2(t)E=EfEi=12SLv02cos2(t)+Sgv022sin2(t)12SLv0
2==Sv02(12L+g2)sin2(t)=p0S2Hv022sin2(t)
La diferencia entre la energa final Ef y la energa inicial Ei, es el trabajo de las

fuerzas exteriores Wext, como vamos a comprobar ms abajo.
Trabajo realizado por las fuerzas exteriores
Consideremos la situacin que se muestra en la figura.

La atmsfera que rodea al dispositivo y que acta sobre la superficie libre de la
rama izquierda, realiza un trabajo que es el producto de la presin
constante p0 por el incremento de volumen (Sx).
W1=p0(Sx)
El aire encerrado en la rama derecha se comprime adiabticamente, por lo que el

trabajo realizado es
W2=0xpSdx0xp0(1+xH)Sdx=p0Sxp0S2Hx2W2=0xpSdx0xp0(1+x
H)Sdx=p0Sxp0S2Hx2
Se ha tenido en cuenta, la relacin existente entre p y p0 en la aproximacin de

que x<<H, deducida al principio de este apartado
La suma de ambas contribuciones es el trabajo realizado por las fuerzas exteriores
Wext=W1+W2=p0S2Hx2=p0S2Hv202sin2(t)Wext=W1+W2=p0S2Hx2=
p0S2Hv022sin2(t)
El trabajo de las fuerzas exteriores modifica la energa del sistema de partculas,

es decir, es igual a la diferencia entre la energa final y la inicial del sistema. Ya
empleamos esta frmula para deducir la ecuacin de Bernoulli.
Como hemos comprobado, el trabajo realizado por el aire que rodea a la columna
de lquido es igual a la diferencia entre la energa final Ef y la energa inicial Ei.
Medida del ndice adiabtico del aire

Si llamamos P1 al periodo de las oscilaciones del lquido con las dos ramas abiertas
y P2 el periodo de las oscilaciones del mismo lquido con una rama taponada en
uno de sus extremos. La relacin entre los cuadrados de los periodos es
P21P22=2g+p0H2gP12P22=2g+p0H2g
Despejando el ndice adiabtico de aire
=2gHp0(P21P221)=2gHp0(P12P221)
Si el lquido contenido en el tubo en forma de U es mercurio, la presin

atmosfrica es p0=gH0, siendo H0=0.76 m, la altura de la columna de mercurio
en el experimento de Torricelli. La frmula final que nos permite calcular el ndice
adiabtico es
=2HH0(P21P221)=2HH0(P12P221)
En la experiencia de laboratorio, descrita en el artculo citado en las referencias,
los datos son los siguientes:
Se emplea mercurio como lquido, que se introduce un un tubo de Pyrex de

19 mm.
La longitud de la columna de lquido es de L=75 cm.
La columna de aire entre la superficie libre y el tapn que cierra uno de los
extremos del tubo en U tiene una longitud H=84 cm.
La presin atmosfrica p0 corresponde a la altura de una columna de
mercurio de H0=75.4 cm.
El valor obtenido para el ndice adiabtico de aire fue =1.38
Actividades
La longitud de la parte horizontal del tubo en U es de 20 cm
Establecemos la altura inicial de fluido en el tubo en U, en el control Altura fluido
Activamos el botn de radio titulado extremos abiertos
Medimos el periodo P1 de las oscilaciones de la columna de lquido con los

extremos del tubo en forma de U abiertos.
Activamos el botn de radio titulado extremo cerrado. Observamos en la rama

derecha una columna de aire que se va a comprimir y expandir adiabticamente.
Medimos el periodo P2 de las oscilaciones. El dato que necesitamos es la presin

atmosfrica que es la presin que ejerce una columna de mercurio de H0=0.76 m
En la parte derecha, se representa el MAS que describe la superficie libre del

lquido, la altura x del lquido en funcin del tiempo t, medida respecto de la
situacin de equilibrio estable, es decir, cuando el lquido en ambas ramas est a
la misma altura. Arbitrariamente, se ha asignado una amplitud de 0.05=5 cm,
para que se cumpla en todos los casos, que el desplazamiento x es siempre
pequeo respecto de la altura H de la columna de aire.
En la parte inferior derecha, una barra horizontal nos muestra las energas que
intervienen.
Cuando los extremos del tubo estn abiertos, la suma de la energa cintica y
potencial se mantiene constante. La energa cintica se transforma en potencial y
viceversa, tal como ocurre en un sistema formado por una partcula unida a un
muelle elstico, o en un pndulo simple.
Cuando un extremo del tubo est cerrado, la suma de la energa cintica y
potencial no es constante, ya que hay que tener en cuenta el trabajo de las
fuerzas exteriores (del aire), este trabajo es igual a la diferencia entre la energa
final y la energa inicial.
Ejemplo:
Establecemos la altura inicial de fluido en 0.40 m.
Medimos el periodo de las oscilaciones con los extremos abiertos del tubo en forma
de U, P1=1.42 s.
Activamos el botn de radio titulado extremo cerrado, que equivale a colocar un

tapn en el extremo derecho del tubo en U. La longitud de la columna de aire
encerrado en la rama derecha es H=0.6
Medimos el periodo de las oscilaciones P2=1.03
Con estos dos datos obtenemos el ndice adiabtico del aire, aplicando la frmula
deducida en el apartado anterior
=20.60.76(1.42211.0321)=1.42=20.60.76(1.42121.0321)=1.42
Comprobacin:
La longitud de la columna de mercurio es L=20.4+0.2=1.0 m, ya que la longitud

de la parte horizontal del tubo en U mide 20 cm.
El periodo de las oscilaciones del lquido con los extremos abiertos es:
P1=21.029.8=1.42sP1=21.029.8=1.42s
La longitud de la columna de aire encerrado en la rama derecha del tubo

es H=1.0-0.4=0.6 m
Sabiendo que el ndice adiabtico del aire es =1.4, y que la presin atmosfrica
equivale a la presin que ejerce una columna de mercurio de H0=0.76 m de altura,
se simplifica la densidad del mercurio en la frmula que calcula el periodo P2.
P2=21.029.8+1.49.80.760.6=1.03sP2=21.029.8+1.49.80.7
60.6=1.03s
Los periodos calculados coinciden con los tiempos P1 y P2 medidos

El Principio de Arquímedes

Cargado por

Información del documentohacer clic para expandir la información del documento

Copyright:

Formatos disponibles

El Principio de Arquímedes

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

El Principio de Arquímedes

Cargado por

Copyright:

Formatos disponibles

El principio de Arqumedes

El principio de Arqumedes afirma que todo cuerpo sumergido en un fluido experimenta un

La explicacin del principio de Arqumedes consta de dos partes como se indica en la

Porcin de fluido en equilibrio con el resto del fluido.

Consideremos, en primer lugar, las fuerzas sobre una porcin de fluido en

Puesto que la porcin de fluido se encuentra en equilibrio, la resultante de las

El peso de la porcin de fluido es igual al producto de la densidad del fluido f por

Si sustituimos la porcin de fluido por un cuerpo slido de la misma forma y

Lo que cambia es el peso del cuerpo slido y su punto de aplicacin que es el

En los casos ms simples, supondremos que el slido y el fluido son homogneos y

La presin debida al fluido sobre la base superior es p1=fgx, y la presin debida al

En el equilibrio tendremos que

Como la presin en la cara inferior del cuerpo p2 es mayor que la presin en la

Como vemos, la fuerza de empuje tiene su origen en la diferencia de presin entre

Con esta explicacin surge un problema interesante y debatido. Supongamos que

Si no hay fluido entre el cuerpo y el fondo del recipiente desaparece la fuerza de

Si se llena un recipiente con agua y se coloca un cuerpo en el fondo, el cuerpo

El principio de Arqumedes sigue siendo aplicable en todos los casos y se enuncia

Cuando un cuerpo est parcialmente o totalmente sumergido en el fluido que le

Flotacin entre dos lquidos no miscibles

La densidad del bloque es un nmero comprendido entre la densidad del aceite

Peso del bloque =empuje del agua + empuje del aceite

S es el rea de la base del bloque, h su altura, y x es la parte del bloque

Supongamos un bloque de h=20 cm de altura, observamos que el bloque est

Despejando en la frmula la densidad del slido, obtenemos el valor de 0.87

Ecuacin fundamental de la esttica de fluidos

La cara superior est en el aceite a una profundidad y. La presin p1 ser igual a la

La cara inferior est en el agua. La presin p2 ser igual a la presin

Teniendo en cuenta que m=hS despejamos el valor de x.

Es el mismo resultado que hemos obtenido aplicando el principio de Arqumedes

La cara superior est a 22 cm de la superficie libre, la presin debida al aceite es

La cara inferior est a 42 cm de la superficie libre (35 cm de aceite y 7 cm de

Se obtiene s=870 kg/m3 0.87 g/cm3

Tomamos el fondo del recipiente como nivel de referencia de la energa potencial.

El centro de masa del cuerpo se encuentra a una altura d-x+h/2. Su energa

El centro de masas de la figura completa, de volumen Sd es d/2

La energa potencial del fluido es Ef=f(Sb)gyf

El valor de la constante aditiva cte, depende de la eleccin del nivel de referencia

En la figura, se representa la energa potencial Ep(x) para un cuerpo de

La funcin presenta un mnimo, que se calcula derivando la energa potencial con

En la posicin de equilibrio, el cuerpo se encuentra sumergido

Energa potencial de un cuerpo que se

El peso del globo Fg=mgj .

Dada la fuerza conservativa, la expresin de la energa potencial asociada, se

ABFdr=EpAEpB Ep=Ep(y)ABFdr=EpAEpB Ep=Ep(y)

La fuerza conservativa peso Fg=mgj est asociada con la energa

Dada la energa potencial, obtenemos la expresin de la fuerza conservativa,

La energa potencial asociada con las dos fuerzas conservativas es

A medida que el globo asciende en el aire con velocidad constante experimenta

Empleando el balance de energa obtenemos la misma conclusin

El trabajo de las fuerzas no conservativas Fnc modifica la energa total (cintica

Energa potencial de un cuerpo

Pueden ocurrir dos casos:

Que el bloque se sumerja parcialmente si la densidad del cuerpo slido es

Tenemos una situacin anloga a la de un cuerpo que se coloca sobre un muelle

El mnimo de Ep se obtiene cuando la derivada de Ep respecto de y es cero, es

El mnimo de Ep se obtiene cuando la derivada de Ep respecto de y es cero, es

El bloque permanece sumergido una longitud x. En esta frmula, se ha