Fuerza sobre un dieléctrico (II)

Potencial constante

Consideremos un condensador plano-paralelo cuyas armaduras son perpendiculares a la superficie de un dieléctrico líquido de constante dieléctrica k. Cuando introducimos el condensador en el recipiente que contiene el dieléctrico su superficie libre se eleva entre las placas del condensador, tal como se muestra en la figura.

Sean a y b las dimensiones de la placa rectangular. Si el dieléctrico está introducido una longitud x en el condensador su capacidad es la suma de las capacidades de dos condensadores uno de longitud x con dieléctrico y otro de longitud a-x en el vacío.

C= k ε bx d + ε 0 b(ax) d

Si la diferencia de potencial en el condensador V permanece constante, entonces la fuerza vertical que se ejerce sobre la superficie libre del líquido dieléctrico, perpendicularmente a las líneas del campo es

F x = ( U x ) V = 1 2 V 2 C x = ε 0 (k1)b 2d V 2

Esta es la fuerza que hace que el líquido ascienda

En el equilibrio la fuerza que ejerce el campo eléctrico sobre el dieléctrico Fx se iguala al peso de la columna de líquido mg=ρ·g·bhd.

h= ε 0 (k1) 2ρg ( V d ) 2

La presión p que hace la columna de líquido de altura h es

p= F x bd = ε 0 (k1) 2 d 2 V 2 = 1 2 ε 0 (k1) E 2

Donde E=V/d es la intensidad del campo eléctrico entre las placas del condensador vacío.

La presión, en la superficie de separación entre dos medios dieléctricos es proporcional al cuadrado de la intensidad del campo eléctrico E.

Balance energético

La capacidad del condensador con dieléctrico aumenta, para que la diferencia de potencial entre las placas se mantenga constante en V, la batería tiene que suministrar carga al condensador. La batería suministra energía para incrementar la energía electrostática del condensador y a la vez, elevar el centro de masas del líquido dieléctrico en h/2.

En la situación inicial la energía asociada al campo eléctrico del condensador vacío es

C 0 = ε 0 ab d U 0 = 1 2 C 0 V 2

En la situación de equilibrio, el condensador tiene en su interior una porción h de dieléctrico líquido. La energía asociada al campo eléctrico del condensador es

C= k ε bh d + ε 0 b(ah) d U ε = 1 2 C V 2 U ε = 1 2 ε 0 (k1)b d h V 2 + 1 2 ε 0 ab d V 2 = 1 4 ε 0 2 (k1) 2 b ρg d 3 V 4 + 1 2 ε 0 ab d V 2 = 1 4 ε 0 2 (k1) 2 b ρg d 3 V 4 + U 0

La porción h de dieléctrico que ha ascendido ha adquirido una energía potencial gravitatoria, su centro de masas se encuentra a h/2 de altura.

U g =ρ(bdh)g h 2 = 1 2 ρgbd h 2 = 1 2 ρgbd ε 0 2 (k1) 2 4 ρ 2 g 2 ( V d ) 4 = 1 8 ε 0 2 (k1) 2 b ρg d 3 V 4

La energía suministrada por la batería será la diferencia entre la energía final Ug+Uε y la energía inicial U0.

W= U g + U ε U 0 = 3 8 ε 0 2 (k1) 2 bd ρg ( V d ) 4

Actividades

El programa interactivo genera un valor aleatorio de la constante dieléctrica k entre 300 y 400, cada vez que se pulsa el botón titulado Nuevo.

Se introduce un valor de la diferencia de potencial V entre las placas del condensador, en el control titulado d. d. potencial.

Ejemplo:

Leemos en el voltímetro una diferencia de potencial de 80 V y medimos la altura que se eleva el fluido dieléctrico en la regla graduada, h=25.97 mm.

25.97· 10 3 = 1 4π·9· 10 9 k1 2·9.8·1.02· 10 3 ( 80 2· 10 4 ) 2 k=368

Carga constante

Un condensador de placas paralelas está fijado en posición vertical, su parte inferior está en contacto con un dieléctrico líquido. Determinar la altura h del líquido entre las placas, despreciando los efectos capilares, suponiendo que la distancia d entre las placas es mucho menor que las dimensiones lineales de las mismas.

La energía inicial del condensador cargado es

U 0 = 1 2 Q 2 C 0 C 0 = ε 0 ab d

donde a (vertical) y b (horizontal) son las dimensiones de las placas rectangulares cargadas con una carga Q.

Cuando las placas entran en contacto con el dieléctrico líquido se ejerce una fuerza vertical sobre el mismo. La carga Q permanece constante. La disminución de la energía electrostática se compensa con el incremento de la energía potencial gravitatoria ya que no hay energía transferida al sistema desde el exterior.

Calculamos la capacidad C del condensador con una porción h de dieléctrico y su energía Uε.

C= ε kbh d + ε 0 b(ah) d U ε = 1 2 Q 2 C

La energía potencial gravitatoria de la porción de fluido es

U g =ρ( bdh )g h 2 = 1 2 ρbdg h 2

Aplicamos el principio de conservación de la energía

U0= Uε+Ug

h 2 + a (k1) h Q 2 ρg( ε 0 a b 2 ) =0

El campo eléctrico del condensador vacío es

E 0 = σ ε 0 = Q ε 0 ab

donde σ es la densidad de carga de cada una de las placas en C/m2.

h 2 + a (k1) h ε 0 E 0 2 a ρg =0 h= 1 2 a (k1) ( 1+ 1+ 4 ε 0 E 0 2 (k1) 2 ρga )

Habitualmente, la altura h de dieléctrico en el condensador es mucho menor que la altura a de las placas.

h ε 0 (k1) ρg E 0 2

Movimiento del líquido dieléctrico

Consideremos un condensador cilíndrico de longitud L cuyo radio interior es a y exterior b, con L>>b

El condensador en posición vertical se introduce parcialmente en un líquido dieléctrico contenido en un recipiente grande, de tal manera, que la superficie del líquido esté lejos de los extremos del condensador. Siendo L0 la parte sumergida

Se aplica un diferencia de potencial entre las placas del condensador, se observa que el líquido dieléctrico se eleva entre sus placas. En este apartado, calculamos la altura máxima que alcanza el líquido dieléctrico y el tiempo que tarda en alcanzarlo.

En el instante t=0, se concetan las placas del condensador a una batería, en el instante t el líquido se ha elevado x tal como se muestra en la figura

El campo eléctrico en el condensador cilíndrico tiene dirección radial y su módulo es (de acuerdo con la ley de Gauss)

E 1 = λ 1 2π ε 0 r E 2 = λ 2 2πεr

Donde ε0 es la permitividad en el vacío y ε es la permitividad en el dieléctrico. λ1 y λ2 son las densidades lineales de carga libre de la parte superior (vacío) e inferior (dieléctrico)

La diferencia de potencial entre las placas es V

V= a b E 1 dr= λ 1 2π ε 0 ln b a V= a b E 2 dr= λ 2 2πε ln b a

La carga del condensador es

Q= λ 1 ( L( L 0 +x) )+ λ 2 ( L 0 +x )= 2π ε 0 V ln b a ( L( L 0 +x) )+ 2πεV ln b a ( L 0 +x )= 2π ε 0 V ln b a ( L( L 0 +x)+ ε ε 0 ( L 0 +x ) )= 2π ε 0 V ln b a ( L+( k1 )( L 0 +x) ),k= ε ε 0

La capacidad del condensador es

C= Q V = 2π ε 0 ln b a ( L+( k1 )( L 0 +x) )

Alternativamente, obtenemos la capacidad del condensador sumando de las capacidades de un condensador cilíndrico vacío de longitud L-(L0+x) y de otro con dieléctrico de longitud L0+x

C= 2π ε 0 ( L( L 0 +x) ) ln b a + 2πε( L 0 +x ) ln b a = 2π ε 0 ln b a ( L+( k1 )( L 0 +x ) )

Potencial constante

La energía del condensador es

U= 1 2 C V 2

La fuerza sobre el dieléctrico es

F= ( dU dx ) V = π ε 0 ln b a ( k1 ) V 2

Esta fuerza es constante

Ecuación del movimiento del líquido dieléctrico

Se trata de un problema de masa variable

d dt ( m dx dt )=F+ρg L 0 π( b 2 a 2 )mg m=ρπ( b 2 a 2 )( L 0 +x )

La segunda fuerza se justifica del siguiente modo: la presión en la parte inferior de la columna de líquido que hay dentro del condensador es mayor que la presión en la superficie libre, la diferencia es ρgL0, donde ρ es la densidad del líquido dieléctrico. Esto da lugar a una fuerza hacia arriba producto de la diferencia de presión por el área de la sección, ρg L 0 π( b 2 a 2 ) , se le añade a la fuera electrostática F y al peso mg

Obtenemos una ecuación diferencial en términos de la variable x

( L 0 +x ) d 2 x d t 2 + ( dx dt ) 2 = ε 0 ( k1 ) V 2 ρ( b 2 a 2 )ln b a +g L 0 ( L 0 +x )g

Para simplificar la ecuación diferencial, definimos las siguientes magnitudes

a 0 = ε 0 ( k1 ) V 2 ρg( b 2 a 2 )ln b a =g t 0 2 ξ= L 0 +x a 0 ,τ= t t 0 ,α= L 0 a 0

Obtenemos la ecuación diferencial

a 0 ξ a 0 t 0 2 d 2 ξ d τ 2 + a 0 2 t 0 2 ( dξ dτ ) 2 = a 0 g+gα a 0 ξ a 0 g a 0 t 0 2 ξ d 2 ξ d τ 2 + a 0 t 0 2 ( dξ dτ ) 2 =g+gαξg a 0 ξ d 2 ξ d τ 2 + a 0 ( dξ dτ ) 2 =g t 0 2 +g t 0 2 αξg t 0 2 ξ d 2 ξ d τ 2 + ( dξ dτ ) 2 =1+αξ

Para resolver la ecuación diferencial se hace el cambio de variable

p= dξ dτ d 2 ξ d τ 2 = dp dτ = dp dξ dξ dτ =p dp dξ

La ecuación diferencial se convierte en

ξp dp dξ + p 2 =1+αξ ξpdp+( p 2 1α+ξ )dξ=0

El procedimiento para obtener la solución de esta ecuación diferencial se explica en la página titulada Ecuaciones diferenciales (I) en el apartado Ecuaciones diferenciales exactas

Cambiamos el nombre de las variables p por t y ξ por x

ξpdp+( p 2 1α+ξ )dξ=0 ( xt )dt+( t 2 1α+x )dx=0 P( t,x )dt+Q( t,x )dx=0,{ P( t,x )=xt Q( t,x )= t 2 1α+x

Comprobamos que no es diferencial exacta

P x = Q t t2t

Buscamos un factor integrante μ

1 P ( Q t P x )=g(x),μ=exp( g(x)dx ) 1 xt ( 2t2 )= 1 x ,μ=exp( 1 x dx )=exp( lnx )=x

Comprobamos que es diferencial exacta

( μP ) x = ( μQ ) t ( x 2 t ) x = ( x( t 2 1α+x ) ) t 2xt=2xt

Obtenemos la solución de la ecuación diferencial

U= μP(t,x)dt+φ(x)= x·xt·dt+φ(x)= x 2 t 2 2 +φ(x) μQ(t,x)= U x x( t 2 1α+x )=x t 2 + φ x φ x = x 2 ( 1+α )x φ(x)= ( x 2 ( 1+α )x )dx = x 3 3 1+α 2 x 2 +C U(x,t)= x 2 t 2 2 + x 3 3 1+α 2 x 2 +C

Regresando a la notación original, cambiando p por t, x por ξ

1 2 ξ 2 p 2 + ξ 3 3 1+α 2 ξ 2 +C=0

La constante C se determina a partir de las condiciones iniciales, en el instante t=0, x=0, dx/dt=0. O bien, τ=0, ξ=L0/a0=α, p=dξ/dτ=0

α 3 3 1+α 2 α 2 +C=0 C= 1+α 2 α 2 α 3 3

La ecuación que hay que integrar es

1 2 ξ 2 ( dξ dτ ) 2 = 1 3 ( α 3 ξ 3 )+ 1+α 2 ( ξ 2 α 2 ) τ= 1 2 α ξ ξdξ 1+α 2 ( ξ 2 α 2 ) 1 3 ( ξ 3 α 3 )

El líquido dieléctrico alcanza la máxima altura cuando la velocidad es nula, p=dξ/dτ=0

1 3 ( ξ 3 α 3 )= 1+α 2 ( ξ 2 α 2 ) 1 3 ( ξα )( ξ 2 +ξα+ α 2 )= 1+α 2 ( ξα )( ξ+α )

Una solución es ξ=α, que es la altura inicial, las otras dos son

1 3 ( ξ 3 α 3 )= 1+α 2 ( ξ 2 α 2 ) 1 3 ( ξα )( ξ 2 +ξα+ α 2 )= 1+α 2 ( ξα )( ξ+α ) 1 3 ( ξ 2 +ξα+ α 2 )= 1+α 2 ( ξ+α ) 1 3 ξ 2 1 2 ( 1+ α 3 )ξ 1 2 α( 1+ α 3 )=0 ξ m = 1 4 ( α+3+3 α 2 + 10 3 α+1 ) ξ 1 = 1 4 ( α+33 α 2 + 10 3 α+1 )= 1 4 ( 3 α 2 + 10 3 α+1 α3 )

El tiempo Τ que tarda en alcanzar la altura máxima ξm es

Τ= 1 2 t 0 α ξ m ξdξ 1 3 ( ξα )( ξ ξ m )( ξ ξ 1 ) = 3 2 t 0 α ξ m ξdξ ( ξα )( ξ m ξ )( ξ ξ 1 )

Buscamos en la tabla de integrales la solución

u a xdx (ax)(xb)(xc) = 2 ac F( λ,p )+2 a b ac E( λ,p ),a>ub>c p= ab ac ,λ=arcsin( au ab )

En este caso, a es ξm, u es b=α y c es ξ1

λ=arcsin(1)= π 2 E( π 2 ,p )=E( p ) F( π 2 ,p )=K( p ) p= ab ac = ξ m α ξ m +| ξ 1 | Τ= 6 t 0 ( 1 ξ m +| ξ 1 | K( p )+ ξ m α ξ m +| ξ 1 | E( p ) )

E(p) y K(p) son integrales elípticas completas

Nota: la expresión del tiempo Τ difiere de la deducida por los autores del primer artículo citado en las Referencias

Datos

a=1/100; %radio interior
b=2/100; %radio exterior
L=20/100; %longitud del condensador
L0=10/100; %profundidad
V=5000; %batería
k=80; %permitividad relativa del líquido
rho=1000; %densidad del líquido

a0=(k-1)*V^2/(4*pi*9e9*rho*9.8*(b^2-a^2)*log(b/a));
alfa=L0/a0;
t0=sqrt(a0/9.8);
xi_m=(alfa+3+3*sqrt(alfa^2+10*alfa/3+1))/4;
xi_1=abs((3*sqrt(alfa^2+10*alfa/3+1)-alfa-3)/4); %valor absoluto
p2=(xi_m-alfa)/(xi_m+xi_1);
[K,E]=ellipke(p2);
T=sqrt(6)*t0*(K/sqrt(xi_m+xi_1)+xi_m*E*sqrt(xi_m+xi_1)/alfa);
disp([xi_m,a0*xi_m-L0]) %altura máxima
disp([T,t0*T]) %tiempo en alcanzarla

La altura máxima es ξm=13.6195 o bien, xm=0.0167 m=1.67 cm. El tiempo que tarda en alcanzar la altura máxima es Τ=0.6038 o bien, tm=0.0179 s

  13.6195    0.0167
   0.6038    0.0179

Carga constante

La fuerza sobre el dieléctrico es

F= ( dU dx ) Q = Q 2 4π ε 0 ( k1 )ln b a ( L+( k1 )( L 0 +x ) ) 2

Que es una fuerza variable

Ecuación del movimiento del líquido dieléctrico

d dt ( m dx dt )= Q 2 4π ε 0 ( k1 )ln b a ( L+( k1 )( L 0 +x ) ) 2 +ρg L 0 π( b 2 a 2 )mg m=ρπ( b 2 a 2 )( L 0 +x )

Reordenando términos

ρπ( b 2 a 2 ) ( dx dt ) 2 +ρπ( b 2 a 2 )( L 0 +x ) d 2 x d t 2 = Q 2 4π ε 0 ( k1 )ln b a ( L+( k1 )( L 0 +x ) ) 2 +ρg L 0 π( b 2 a 2 )ρπ( b 2 a 2 )( L 0 +x )g ( L 0 +x ) d 2 x d t 2 + ( dx dt ) 2 = Q 2 4π ε 0 ( k1 )ln b a ( L+( k1 )( L 0 +x ) ) 2 ρπ( b 2 a 2 ) +g L 0 ( L 0 +x )g

Para simplificar la ecuación diferencial, definimos las siguientes magnitudes

a ˜ 0 = Q 2 4 π 2 ε 0 ( k1 )ln b a ρ( b 2 a 2 ) L 2 g =g t ˜ 0 2 ξ= L 0 +x a ˜ 0 ,τ= t t ˜ 0 ,α= L 0 a ˜ 0 ,σ=( k1 ) a ˜ 0 L

El resultado es

a ˜ 0 ξ a ˜ 0 t ˜ 0 2 d 2 ξ d t 2 + a ˜ 0 2 t ˜ 0 2 ( dx dt ) 2 = g L 2 a ˜ 0 ( L+( k1 )( L 0 +x ) ) 2 + a ˜ 0 gα a ˜ 0 ξg ξ a ˜ 0 d 2 ξ d τ 2 + a ˜ 0 ( dξ dτ ) 2 = g t ˜ 0 2 ( 1+( k1 ) ξ a ˜ 0 L ) 2 +αg t ˜ 0 2 ξg t ˜ 0 2 ξ d 2 ξ d τ 2 + ( dξ dτ ) 2 = 1 ( 1+σξ ) 2 +αξ

Esta ecuación diferencial no tine solución analítica. Se resuelve aplicando procedimientos numéricos, con las condiciones iniciales, en el instante t=0, x=0, dx/dt=0. O bien, τ=0, ξ=L0/a0=α, p=dξ/dτ=0.

Referencias

Rafael Nardi, Nivaldo A Lemos. Dynamics of a liquid dielectric attracted by a cylindrical capacitor. Eur. J. Phys. 28 (2007) pp. 125-134

I. S. Gradshteyn, I. M. Ryzhik. Table of Integrals, Series, and Products. Seventh Edition. Elsevier (2007). Página 257, 3.132, n° 5