El principio de Arquímedes y el efecto de la tensión superficial

Sea un objeto de masa m parcialmente sumergido en un líquido de densidad ρf. La superficie del fluido se deforma estableciéndose una línea de contacto entre el cuerpo y el líquido a una profundidad h por debajo de la superficie libre sin deformar. La recta tangente a la superficie del fluido en los puntos de contacto hace un ángulo α con la vertical

Hay dos fuerzas que se oponen al peso del objeto:

De acuerdo con Keller

  1. La componente vertical de la fuerza debida a la presión es equivalente al peso del líquido desplazado por la porción del cuerpo sumergido y la porción de altura h delimitada por un contorno rojo por encima de la línea de contacto. El volumen desplazado corresponde al sombreado de color rosa de la figura

  2. La componente vertical de la fuerza debido a la tensión superficial es igual al volumen sombreado de color amarillo (menisco).

La tensión superficial contribuye a la flotación de un cuerpo siempre que su longitud característica R sea mucho menor que la longitud capilar lc definida

l c = γ ρ f g

donde γ es la tensión superficial, ρf es la densidad del líquido y g=9.8 m/s2 la aceleración de la gravedad.

Alternativamente, la tensión superfical juega un papel relevante, siempre que el número adimensional de Bond, Bo sea pequeño

Bo= ρ f g R 2 γ = ( R l c ) 2

Equilibrio de una porción de fluido

Consideremos la porción de fluido ACDA de espesor Δx que rodea al cuerpo parcialmente sumergido de altura h.

Las fuerzas sobre dicha porción de fluido son

Equilibrio en la dirección horizontal

γsinα·Δx+ ρ f g h 2 2 ΔxγΔx=0 h= 2γ( 1sinα ) ρ f g = l c 2( 1sinα )

De esta manera relacionamos dos magnitudes importantes h y el ángulo α

Angulo de contacto

El ángulo de contacto θc entre una superficie sólida plana y la superficie del líquido se muestra a la izquierda en la figura. Es el ángulo entre dicho plano y la recta tangente a la superficie del líquido en el punto de contacto

A la derecha de la figura, se traza la tangente a la superficie del sólido parcialmente sumergido en el punto de contacto (sólido-líquido). Se traza la tangente a la superficie del líquido en dicho punto (hace un ángulo α con la vertical). El ángulo de contacto θc entre las dos rectas se muestra en la figura

La relación entre águlos θc, α y φ es

θ c ( π 2 α )+ φ =π θ c +α+ φ = 3π 2

Ejemplos sencillos

Un clindro de madera de pequeñas dimensiones flota en agua

  1. Utilizando solamente el principio de Arquímedes y despreciando el efecto de la tensión superficial

  2. El peso es igual al empuje

    ρ s π R 2 L= ρ f π R 2 h h= ρ s ρ f L h= 510 998 10=5.1mm

  3. Incluimos los efectos de la tensión superficial, sabiendo que el ángulo de contacto es θc

  4. La fuerza que ejerce la tensión superficial es T=γ·2πR. En el equilibrio

    mg+Tcos θ c =E ρ s π R 2 Lg+2πRγcos θ c = ρ f π R 2 hg h= ρ s ρ f L+ 2γcos θ c ρ f Rg

     

Flotación de un largo cilindro

El área sombreada de color gris vale

y 0 r 2xdy=2 y 0 r r 2 y 2 dy

Haciendo el cambio de variable

y=rsinθ,dy=rcosθ·dθ

Se obtiene

2 r 2 θ 0 π/2 cos 2 θ·dθ=2 r 2 θ 0 π/2 1+cos2θ 2 dθ= r 2 ( θ+ sin2θ 2 | θ 0 π/2 )= r 2 2 ( π2 θ 0 sin2 θ 0 )

Fuerzas sobre el cilindro

Consideremos un cilindro de radio r y densidad ρs de longitud infinita flotando en la superficie de un fluido de densidad ρf. La geometría plana nos facilita los cálculos

El área S1 sombreada de color rosa en la figura es

S 1 = r 2 2 ( π2( π 2 φ )sin( 2( π 2 φ ) ) )+2rsinφ·h S 1 = r 2 φ r 2 sinφcosφ+2rhsinφ

La fuerza de empuje es proporcional a S1

F B = ρ f g S 1 = ρ f g( r 2 φ r 2 sinφ·cosφ+2hrsinφ )

La componente vertical de la fuerza que ejerce la tensión superficial es proporcional a

F S =2γcosα

El peso del cilindro se equilibra con la fuerza de empuje FB y la componente vertical de la fuerza que ejerce la tensión superficial FS

F B + F S = ρ s π r 2 g ρ f g r 2 ( φsinφ·cosφ+2 h r sinφ )+2γcosα= ρ s π r 2 g

Definiendo magnitudes adimensionales

R= r l c ,H= h l c , l c = γ ρ f g π ρ s ρ f =φsinφ·cosφ+2 H R sinφ+ 2 R 2 cosα { θ c +α+φ= 3π 2 H= 2( 1sinα )

La densidad relativa ρsf es una función del ángulo α

Sea R=0.5 y el ángulo de contacto θc=2π/3

R=0.5; % radio adimensional
th_c=2*pi/3; %ángulo de contacto
H=@(x) sqrt(2*(1-sin(x)));
phi=@(x) 3*pi/2-th_c-x; %x es el ángulo alfa
f=@(x) (phi(x)-sin(phi(x)).*cos(phi(x))+2*H(x).*sin(phi(x))/R+2*cos(x)/R^2)/pi;
fplot(f,[0,pi/2])
grid on
set(gca,'XTick',0:pi/12:pi/2)
set(gca,'XTickLabel',{'0','\pi/12','\pi/6','\pi/4','\pi/3','5\pi/12','\pi/2'})
xlabel('\alpha')
ylim([0,5])
ylabel('\rho_s/\rho_f')
title('Densidad')

La densidad relativa ρsf presenta un máximo para un cierto valor del ángulo α. Por encima de este valor máximo, el cilindro no puede mantenerse en equilibrio y por tanto, se hundirá

Perfil

En la figura, se muestra el perfil y=y(x) de la superficie de agua que rodea el cilindro

De la ecuación de Young-Laplace

ρgy= γ R c , R c = ( 1+ ( dy dx ) 2 ) 3 2 d 2 y d x 2

Donde Rc es el radio de curvatura que es negativo para esta función, cuya pendiente disminuye con x. Una función similar es y=-e-x, la derivada primera, dy/dx=e-x>0, la derivada segunda, d2y/dx2=-e-x<0

Tenemos una situación análoga a la estudiada en la página titulada Soluciones simples de la ecuación de Young-Laplace

l c 2 y= d 2 y d x 2 ( 1+ ( dy dx ) 2 ) 3 2

Multiplicamos ambos miembros por dy/dx e integramos

l c 2 y dy dx = d 2 y d x 2 dy dx ( 1+ ( dy dx ) 2 ) 3 2 l c 2 y·dy= d 2 y d x 2 dy dx ( 1+ ( dy dx ) 2 ) 3 2 dx l c 2 y·dy= 1 2 l c 2 y 2 u ( 1+ u 2 ) 3/2 ·du= 1 1+ u 2 ,u(x)= dy dx ,du= d 2 y d x 2 dx 1 2 l c 2 y 2 = 1 1+ ( dy dx ) 2 +c

Las condiciones de contorno son: cuando x→∞, y→0, dy/dx→0. El resultado es c=1

Por otra parte, en la posición

x=rsinφ,y=h, dy dx =tan( π 2 α )=cotα

Estableceremos la relación entre h y el ángulo α que ya habíamos obtenido, al principio de la página, a partir del equilibrio horizontal de dicha porción de líquido

1 2 l c 2 h 2 =1 1 1+ ( cosα sinα ) 2 1 2 l c 2 h 2 =1sinα h= l c 2( 1sinα )

En términos de magnitudes adimensionales

H= h l c = 2( 1sinα )

Expresamos la ecuación diferencial

1 2 l c 2 y 2 =1 1 1+ ( dy dx ) 2

en términos de magnitudes adimensionales X=x/lc e Y=y/lc. Despejamos dY/dX e integramos

1 2 Y 2 =1 1 1+ ( dY dX ) 2 dY dX = Y 4 Y 2 2 Y 2 2 Y 2 Y 4 Y 2 dY= dX +C

Para la primera integral hacemos el cambio de variable

Y=2sint,dY=2cost·dt 24 sin 2 t 2sint 44 sin 2 t 2cost·dt= dt sint 2 sint·dt =ln( tan t 2 )+2cost

La primera integral, está resuelta en la página titulada Integrales. La segunda es inmediata

>> syms x;
>> int(1/sin(x))
ans =log(tan(x/2))

Para expresar los resultados de la primera integral en términos de Y, relacionamos tan(t) y tan(t/2)

tant= 2tan t 2 1 tan 2 t 2 tant· tan 2 t 2 +2tan t 2 tant=0 tan t 2 = 2± 4+4 tan 2 t 2tant = 1± 1 cost tant = cost±1 sint = 1 sin 2 t ±1 sint tan t 2 = 1 Y 2 4 ±1 Y 2 = 4 Y 2 ±2 Y

Como Y es negativo descartamos la solución positiva de la ecuación de segundo grado. El resultado de la integral en Y (primer miembro)

ln( 2 4 Y 2 Y )+2 1 Y 2 4 =ln( 2 4 Y 2 Y )+ 4 Y 2

La integración de la ecuación diferencial es

ln( 2+ 4 Y 2 Y )+ 4 Y 2 =X+C

La constante C se calcula, sabiendo que para x=r·sinφ entonces y=-h o bien, para X=R·sinφ entonces, Y=-H

ln( 2+ 4 Y 2 H )+ 4 H 2 =Rsinφ+C

El resultado final, es una ecuación implícita, con Y<0

X=ln( 2+ 4 Y 2 Y )+ 4 Y 2 ln( 2+ 4 H 2 H ) 4 H 2 +Rsinφ

Los datos de este ejemplo son

Representamos el perfil de la superficie del líquido en las proximidades del cilindro

R=0.5; % radio adimensional
th_c=2*pi/3; %ángulo de contacto
D=2; %cociente entre densidades rho_s/rho_f
h=@(x) sqrt(2*(1-sin(x))); %altura
phi=@(x) 3*pi/2-th_c-x; %x es el ángulo alfa
f=@(x) phi(x)-sin(phi(x))*cos(phi(x))+2*h(x)*sin(phi(x))/R+2*cos(x)/R^2-pi*D;
alfa=fzero(f,pi/6);
phi_0=phi(alfa);
H=h(alfa);
c=-log((2+sqrt(4-H^2))/H)-sqrt(4-H^2)+R*sin(phi_0);
y=linspace(-H,-0.01,50);
x=@(y) log(-(2+sqrt(4-y.^2))./y)+sqrt(4-y.^2)+c;
hold on
plot(x(y), y,'b')
plot(-x(y), y,'b')
fp=fplot(@(x) R*sin(x), @(x) -H+R*cos(phi_0)+R*cos(x), [0,2*pi],'k');
fill(fp.XData, fp.YData,[0.7,0.7,0.7])
plot(0,-H+R*cos(phi_0),'ro','markersize',3,'markerfacecolor','r')
plot(R*sin(phi_0),-H,'bo','markersize',3,'markerfacecolor','b')
quiver(R*sin(phi_0),-H, sin(alfa),cos(alfa),'r')
hold off
grid on
axis equal
xlim([0,4])
xlabel('X')
ylabel('Y')
title('Perfil')

Una flecha de color rojo, señala la dirección de la fuerza que ejerce la tensión superficial, tangente a la superficie del líquido en los puntos de contacto (Rsinφ, -H) con el cilindro (en color azul). El ángulo α es el que hace la tangente con la vertical.

Los resultados son: α=66°, φ=84°, H=-0.42

>> alfa*180/pi
ans =   65.9442
>> phi_0*180/pi
ans =   84.0558
>> H
H =    0.4168

Mostramos la parte derecha de la figura

Observamos que la deformación de superficie libre del líquido es notable cerca del cilindro y es casi inapreciable, a una distancia de más de ocho radios R

Flotación de una esfera

Volumen de un casquete esférico, comprendido entre z0 y r

V= z 0 r π x 2 dz= z 0 r π( r 2 z 2 )dz= π r 3 3 ( 23 z 0 r + z 0 3 r 3 )

Fuerzas sobre la esfera

Consideremos una esfera de radio r y densidad ρs flotando en la superficie de un fluido de densidad ρf.

El volumen V1 sombreado de color rosa en la figura es la suma de dos volúmenes: un casquete esférico con z0=rcosφ y un cilindro de radio rsinφ y altura h

V 1 = π r 3 3 ( 23 rcosφ r + ( rcosφ ) 3 r 3 )+π ( rsinφ ) 2 h V 1 = π r 3 3 ( 23cosφ+ cos 3 φ+3 h r sin 2 φ )

La fuerza de empuje es proporcional a V1

F B = ρ f g V 1 = ρ f g π r 3 3 ( 23cosφ+ cos 3 φ+3 h r sin 2 φ )

La componente vertical de la fuerza que ejerce la tensión superficial es proporcional a

F S =γ( 2πrsinφ )cosα

r·sinφ es la longitud de la línea de contacto

El peso de la esfera se equilibra con la fuerza de empuje FB y la componente vertical de la fuerza que ejerce la tensión superficial FS

F B + F S = ρ s π r 2 g ρ f g π r 3 3 ( 23cosφ+ cos 3 φ+3 h r sin 2 φ )+γ( 2πrsinφ )cosα= ρ s 4 3 π r 3 g

En términos de magnitudes adimensionaless

R= r l c ,H= h l c , l c = γ ρ f g ρ f g R 2 3 γ ρ f g ( 23cosφ+ cos 3 φ+3 H R sin 2 φ )+γ( 2sinφ )cosα= ρ s 4 3 R 2 γ ρ f g g ρ s ρ f = 1 4 ( 23cosφ+ cos 3 φ+3 H R sin 2 φ )+ 3 2 R 2 sinφcosα

En este caso, no disponemos una relación entre H y α por lo que hemos de resolver la ecuación diferencial que determina la forma de la superficie del líquido en las proximidades de la esfera

Perfil

En este caso, la superficie de líquido que rodea la esfera es una superficie de revolución generada por la función z(x) al girar alrededor del eje Z, semejante a la de la figura.

Tal como se ha descrito en la página titulada Presión producida por la curvatura de una superficie tiene dos radios de curvatura

De la ecuación de Young-Laplace, similar al cilindro

ρ f gz=γ( 1 r 1 + 1 r 2 ) r 1 = ( 1+ ( dz dx ) 2 ) 3 2 | d 2 z d x 2 | , r 2 = x 1+ ( dz dx ) 2 | dz dx |

Queda pendiente la solución de la ecuación diferencial con las condiciones de contorno especificadas, por procedimientos numéricos, para determinar la función z(x) que al girar alrededor del eje Z genera la superficie de revolución que describe la superficie del líquido en las proximidades de la esfera

Referencias

David Naylor, Scott S.H. Tsai. Archimedes’ principle with surface tension effects in undergraduate fluid mechanics. International Journal of Mechanical Engineering Education, SAGE, DOI: 10.1177/03064190211055431

Dominic Vella, Duck-Gyu Lee, Ho-Young Kim. The Load Supported by Small Floating Objects. Langmuir 2006, 22, 5979-5981