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:
- La componente vertical de la fuerza que ejerce la presión hidrostática p sobre la superficie mojada del cuerpo.
- La componente vertical de la fuerza que ejerce la tensión superficial a lo largo de la línea de contacto líquido-aire. La fuerza es el producto de la tensión superficial γ por la longitud de la línea de contacto líquido-cuerpo
De acuerdo con Keller
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
-
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
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
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
El peso del fluido, mfg
La fuerza que ejerce la tensión superficial en el punto de contacto C, γΔx que hace un ángulo α con la vertical. Esta fuerza es igual y de sentido contrario a la que ejerce el fluido sobre el cuerpo
La fuerza que ejerce la tensión superficial del resto del fluido sobre dicha porción en en el punto A, γΔx. Supondremos que A estará suficientemente alejado del cuerpo, para que la superficie libre del fluido en A sea casi horizontal
La fuerza vertical que ejerce la presión hidrostática pDC en la cara inferior DC (no interesa)
La fuerza horizontal que ejerce la presión hidrostática pAB en la cara AD de área h·Δx. Como vemos a la derecha de la figura, esta presión se incrementa linealmente con la profundidad, un ejemplo que hemos estudiado en la página titulada Ecuación fundamental de la estática de fluidos en el apartado 'Fuerza sobre el muro de un embalse'. La fuerza es
La resultante de las fuerzas debidas a la presión está aplicada a la profundidad 2h/3
Equilibrio en la dirección horizontal
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
Ejemplos sencillos

Un clindro de madera de pequeñas dimensiones flota en agua
- Densidad de la madera, ρs=510 kg/m3
- Densidad del agua a 20 °C, ρf=998 kg/m3
- Tensión superficial del agua, γ=0.0728 N/m
- Radio del cilindro, R 9 mm
- Altura del cilindro, L=10 mm
Utilizando solamente el principio de Arquímedes y despreciando el efecto de la tensión superficial
Incluimos los efectos de la tensión superficial, sabiendo que el ángulo de contacto es θc
El ángulo de contacto es θc=0°
Se recubre la madera con una capa de pintura repelente al agua, el ángulo de contacto cambia a θc=165°

El peso es igual al empuje
La fuerza que ejerce la tensión superficial es T=γ·2πR. En el equilibrio

La longitud capilar
Que es de la mismo orden de magnitud que el radio R=9 mm, por lo que la tensión superficial juega un papel relevante. Comprobamos que el número de Bond es pequeño

Flotación de un largo cilindro
El área sombreada de color gris vale

Haciendo el cambio de variable
Se obtiene
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
La fuerza de empuje es proporcional a S1
La componente vertical de la fuerza que ejerce la tensión superficial es proporcional a
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
Definiendo magnitudes adimensionales
La densidad relativa ρs/ρf 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 ρs/ρf 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
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
Multiplicamos ambos miembros por dy/dx e integramos
Las condiciones de contorno son: cuando x→∞, y→0, dy/dx→0. El resultado es c=1
Por otra parte, en la posición
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
En términos de magnitudes adimensionales
Expresamos la ecuación diferencial
en términos de magnitudes adimensionales X=x/lc e Y=y/lc. Despejamos dY/dX e integramos
Para la primera integral hacemos el cambio de variable
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)
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)
La integración de la ecuación diferencial es
La constante C se calcula, sabiendo que para x=r·sinφ entonces y=-h o bien, para X=R·sinφ entonces, Y=-H
El resultado final, es una ecuación implícita, con Y<0
Los datos de este ejemplo son
- Radio adimensional, R=0.5
- Densidad relativa, D=ρs/ρf=4
- Angulo de contacto, θc=2π/3, (120°)
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

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
La fuerza de empuje es proporcional a V1
La componente vertical de la fuerza que ejerce la tensión superficial es proporcional a
2π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
En términos de magnitudes adimensionaless
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
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