Fenómenos capilares

Meniscos

En las proximidades de la pared de un recipiente, una molécula del líquido (señalada en color rojo) experimenta las siguientes fuerzas:

Supondremos despreciable la fuerza que ejercen sobre la molécula considerada, las moléculas de vapor por encima de la superficie del líquido.

En la figura de la izquierda, se muestran las fuerzas sobre dos moléculas, una que está muy cerca de la pared y otra que está más alejada.

En la figura de la derecha, se muestra la resultante de dichas fuerzas. La superficie es siempre normal a la resultante. Cuando las moléculas están alejadas de la pared, la resultante debida al peso y a las fuerzas de cohesión (las fuerzas de adherencia son despreciables) es vertical hacia abajo, la superficie es entonces, horizontal.

Pueden ocurrir dos casos según sea la intensidad de las fuerzas de cohesión y adherencia.

Recibe el nombre de ángulo θ de contacto, al formado por la tangente a la superficie del menisco en el punto de contacto con la pared. Este ángulo es agudo cuando el líquido moja y es obtuso cuando el líquido no moja.

Fenómenos capilares. Ley de Jurín

Si se coloca un capilar verticalmente en un recipiente de líquido que moje, el líquido asciende por el capilar, hasta alcanzar determinada altura. Si el líquido no moja, el nivel de líquido en el capilar es menor que en el recipiente.

Debido a la curvatura de una superficie se produce una sobrepresión en su interior, que ya hemos estudiado en anteriores páginas.

La superficie del menisco en el capilar se puede considerar como un casquete esférico de radio R. La relación entre el radio del capilar r, el radio del menisco R y el ángulo de contacto θ, es.

r=Rcosθ

Debido a la curvatura de la superficie habrá una sobrepresión hacia el centro del menisco, que de acuerdo con la ley de Laplace (superficie de una cara), valdrá

Δp= 2γ R = 2γ r cosθ

Por efecto de esta sobrepresión, el líquido asciende una altura h.

Δp=ρgh

La altura h a la que asciende el nivel del líquido en el capilar será

h= 2γ rgρ cosθ

Esta expresión es la denominada ley de Jurín:

Otra forma de obtener esta relación, el peso de la columna de líquido se compensa con la componente vertical de la fuerza que ejerce la tensión superficial

( 2πrγ )cosθ=ρ( π r 2 h )g h= 2γcosθ ρrg

La altura a la que se eleva o desciende un líquido en un capilar es directamente proporcional a su tensión superficial y está en razón inversa a la densidad del líquido y del radio del tubo.

El líquido asciende por un tubo capilar

En el apartado anterior, se ha descrito el fenómeno de la elevación h de un líquido en un capilar de radio r conocido, a partir de la cual, se puede obtener la tensión superficial γ del líquido.

En este apartado, se describe el aspecto dinámico de dicha situación física. Introducimos el tubo capilar en posición vertical en un depósito grande de líquido y observaremos como se incrementa la altura de líquido en el capilar a medida que transcurre el tiempo.

Cuando el líquido de viscosidad η, asciende por el tubo capilar de radio r supondremos que fluye en régimen laminar, aplicamos la ley de Poiseuille,

G= π 8 Δp· r 4 η·h

donde G es el gasto o volumen de líquido que fluye en la unidad de tiempo,

G=π r 2 dh dt

dh/dt es el incremento de altura de líquido en el capilar en la unidad de tiempo y h es la altura de la columna de fluido en el tubo capilar en el instante t.

Δp es la diferencia de presión que hace que el líquido ascienda.

Δp= 2γ r ρgh

Cuando esta diferencia de presión se hace cero, Δp=0, se alcanza la situación de equilibrio estudiada en el apartado anterior. El líquido deja de ascender por el tubo capilar. La altura máxima a la que asciende es

h máx = 2γ ρgr

La ecuación diferencial que describe la variación de la altura h de líquido en el tubo capilar en función del tiempo t es,

π r 2 dh dt = π r 4 8ηh ( 2γ r ρgh )

o bien,

dh dt = abh h a= γr 4η b= ρg r 2 8η

Integramos la ecuación

0 h hdh abh = 0 t dt

Después de hacer algunas operaciones obtenemos la ecuación

h+ a b ln abh a =bt

Esta función implícita t=f(h) se puede hacer explícita h=g(t) expresándola en términos de la función W de Lambert

ln( 1 b a h )= b a (h+bt) exp( b a h b 2 a t )= b a h+1

Comparamos la ecuación, ab·x+c=d·x+f con su expresión en términos de la función W de Lambert

x= 1 blna W( b d ( lna ) a cbf/d ) f d

Sustituyendo la variable x por h, obtenemos

h= 1 b a lne W( b a b a (lne)exp( b 2 a t b a ·1 b a ) ) 1 b a h= a b W( exp( b 2 a t1 ) )+ a b

%agua
tensionSup=0.0728; %tensión superficial del agua
densidad=1000; %densidad
viscosidad=0.00105; %viscosidad
radio=0.1/1000; %radio capilar en m
a=tensionSup*radio/(4*viscosidad);
b=densidad*9.8*radio^2/(8*viscosidad);

f=@(t) a*lambertw(0,-exp(-1-b^2*t/a))/b+a/b;
fplot(f,[0,250])
grid on
xlabel('t(s)')
ylabel('h(m)')
title('Fenómeno capilar')

Cuando t tiende a infinito h tiende al cociente a/b.

h máx = a b = 2γ ρgr

>>  hMax=2*tensionSup/(densidad*9.8*radio)
hMax =    0.1486

La altura máxima hmax es independiente de la viscosidad del fluido η, pero el tiempo que tarda en alcanzar aproximadamente esta altura, depende de la viscosidad, una situación análoga a la carga de un condensador a través de una resistencia, o a la velocidad que alcanza una esfera que se deja caer en el seno de un fluido viscoso.

La ecuación que nos da h en función del tiempo t, no es idéntica, a la ecuación que nos da la carga del condensador en función del tiempo.

Actividades:

Se pulsa el botón titulado Nuevo.

El tubo capilar en posición vertical se pone en contacto con el líquido y se observa como asciende el líquido en el tubo capilar.

Los datos de los líquidos examinados se indican en el siguiente cuadro, expresados en el S.I. de unidades.

Líquido Tensión superficial γ Densidad ρ Viscosidad η
Agua 0.0728 1000 0.00105
Aceite 0.0331 900 0.113
Alcohol 0.0228 790 0.00122
Glicerina 0.0594 1260 1.393

Fuente: Manual de Física Elemental. Koshkin N. I., Shirkévich M. G. Edtorial Mir (1975)

El aceite y la glicerina tienen viscosidades η elevadas, por lo que estos líquidos ascienden lentamente en el tubo capilar en comparación con el agua y en el alcohol. El radio r del capilar juega también un papel importante. Aunque un radio mayor disminuye la diferencia de presión Δp, aumenta sin embargo, la rapidez con la que asciende el líquido


Una descripción más completa

Cómo asciende un líquido en un tubo capilar es un fenómeno complejo para el cual existen varias explicaciones, la más sencilla es la descrita en el apartado anterior. Se trata de un problema de masa variable bajo la acción de las fuerzas debidas a la tensión superficial, al peso y a las fuerzas de rozamiento (viscoso) y otras fuerzas.

Se supondrá que el flujo de líquido a través del tubo capilar tiene un perfil de velocidades parabólico como indica la ley de Poiseuille

Se inserta un tubo capilar vertical en un depósito de líquido muy grande. En el instante t=0, h=0, y dh/dt=0. Debido a la componente vertical de la fuerza que ejerce la tensión superficial, el líquido se acelera dentro del tubo capilar. Mientras el fluido asciende, el rozamiento con el tubo y el peso del líquido se oponen y la columna de líquido empieza a decelerarse hasta que finalmente se detiene cuando alcanza la altura máxima h dada por la ley de Jurín para t→∞

h = 2γ rgρ cosθ

Pueden ocurrir dos casos, que la altura de fluido h(t) tienda hacia el valor asintótico h

La ecuación diferencial que describe estos dos comportamientos es

( h+ 37 36 r ) d 2 h d t 2 + 7 6 ( dh dt ) 2 + η ρ ( 8h r 2 + 2 r ) dh dt +gh= 2γ ρr cosθ

En el estado estacionario, se alcanza la altura h

Resolvemos esta ecuación diferencial por procedimientos numéricos con los siguientes datos

Para cinco radios del tubo capilar: r=0.1 0.2, 0.3, 0.4 y 0.5 mm

eta=2.2e-4; %viscosidad
gamma=1.67e-2; %tensión superficial
th=26*pi/180; %ángulo de contacto
rho=710; %densidad
hold on
for k=1:5
    r=0.0001*k; %radio del tubo capilar
    f=@(t,x) [x(2);(2*gamma*cos(th)/(rho*r)-7*x(2)^2/6-
eta*(8*x(1)/r^2+2/r)*x(2)/rho-9.81*x(1))/(x(1)+37*r/36)]; 
    [t,x]=ode45(f,[0,1],[0,0]);
    plot(t,x(:,1))
end
hold off
grid on
legend('i=1','i=2','i=3','i=4','i=5','Location','best')
xlabel('t')
ylabel('h');
title('Capilar')

Observamos los dos tipos de comportamiento, oscilando alrededor de h para r=0.5, 0.4 y 0.3 mm, sin oscilar r=0.2 y 0.1 mm

Una solución analítica que no contempla las oscilaciones alrededor de la altura final h es

h(t)= 2γcosθ ρgr { 1( 1+t 6g 7r +α t 3 )exp( t 6g 7r ) } α 7ηγcosθ 12 ρ 2 r 5 1 15 ( 6g 7r ) 3

En el estado estacionario, t→∞, h(t)→h

La velocidad y la aceleración son, respectivamente,

dh dt = 2γcosθ ρgr ( 6g 7r 3αt+ t 2 α 6g 7r )texp( t 6g 7r ) d 2 h d t 2 = 2γcosθ ρgr { 6g 7r ( 6α+ ( 6g 7r ) 3/2 )t+6α 6g 7r t 2 α 6g 7r t 3 }exp( t 6g 7r )

Representamos la altura h(t), la velocidad dh/dt y la aceleración d2h/dt2, con los datos siguienets

  1. Altura, h(t)

  2. eta=2.2e-4; %viscosidad
    gamma=1.67e-2; %tensión superficial
    th=26*pi/180; %ángulo de contacto
    rho=710; %densidad
    r=0.0005; %radio del tubo capilar
    alfa=7*eta*gamma*cos(th)/(12*rho^2*r^5)-sqrt(6*9.81/(7*r))^3/15;
    h=@(t) 2*gamma*cos(th)*(1-(1+t*sqrt(6*9.81/(7*r))+alfa*t.^3).
    *exp(-t*sqrt(6*9.81/(7*r))))/(rho*9.81*r);
    H=2*gamma*cos(th)/(rho*9.81*r);
    hold on
    line([0,0.15],[H,H],'lineStyle','--')
    fplot(h,[0,0.15])
    hold off
    grid on
    ylim([0,0.009])
    xlabel('t')
    ylabel('h')
    title('Tubo capilar')

  3. Velocidad, dh/dt

  4. eta=2.2e-4; %viscosidad
    gamma=1.67e-2; %tensión superficial
    th=26*pi/180; %ángulo de contacto
    rho=710; %densidad
    r=0.0005; %radio del tubo capilar
    alfa=7*eta*gamma*cos(th)/(12*rho^2*r^5)-sqrt(6*9.81/(7*r))^3/15;
    h=@(t) 2*gamma*cos(th)*(1-(1+t*sqrt(6*9.81/(7*r))+alfa*t.^3).*
    exp(-t*sqrt(6*9.81/(7*r))))/(rho*9.81*r);
    H=2*gamma*cos(th)/(rho*9.81*r);
    hold on
    line([0,0.15],[H,H],'lineStyle','--')
    fplot(h,[0,0.15])
    hold off
    grid on
    ylim([0,0.009])
    xlabel('t')
    ylabel('h')
    title('Tubo capilar')

  5. Aceleración, d2h/dt2

  6. eta=2.2e-4; %viscosidad
    gamma=1.67e-2; %tensión superficial
    th=26*pi/180; %ángulo de contacto
    rho=710; %densidad
    r=0.0005; %radio del tubo capilar
    alfa=7*eta*gamma*cos(th)/(12*rho^2*r^5)-sqrt(6*9.81/(7*r))^3/15;
    v=@(t) 2*gamma*cos(th)*(6*9.81/(7*r)-(6*alfa+(6*9.81/(7*r))^(3/2))*t+
    6*alfa*sqrt(6*9.81/(7*r))*t.^2-6*9.81*alfa*t.^3/(7*r)).
    *exp(-t*sqrt(6*9.81/(7*r)))/(rho*9.81*r);
    fplot(v,[0,0.15])
    grid on
    ylim([-25,150])
    xlabel('t')
    ylabel('d^2h/dt^2')
    title('Tubo capilar')

    La aceleración producida por la componente vertical de la fuerza que origina la tensión superficial, va disminuyendo con el incremento del peso de la columna de fluido y la fuerza de rozamiento

Referencias

Peiris M. G. C. , Tenmakone K.. Rate of rise of a liquid in a capillary tube. Am. J. Phys. 48 (5) May 1980, pp. 415

BoHua Sun. The monotonic rising and oscillating of capillary driven flow in circular cylindrical tubes. AIP Advances. February 2021

BoHua Sun. Singularity-free approximate analytical solution of capillary rise dynamics. SCIENCE CHINA Physics, Mechanics &! Astronomy. August 2018, Vol. 61, No. 8: 084721