Fuerza de Lorentz sobre una esfera en rotación
De Laplace
Contenido |
1 Introducción
Suponemos una distribución de carga que posee simetría esférica alrededor de un punto central , de forma que su densidad de carga verifica
![\rho(\mathbf{r})=\rho(r')\,](/wiki/images/math/7/f/6/7f63bab59f49dff1efea9dc3210e5ea2.png)
![\mathbf{r}'=\mathbf{r}-\mathbf{r}_C\qquad r'=|\mathbf{r}'|](/wiki/images/math/7/c/4/7c459e098abe06243bd7581edd2e6a6c.png)
Suponemos que esta distribución de carga está localizada, de forma que tiende a cero rápidamente cuando r' crece
Esta distribución de carga se mueve rígidamente, de forma que la velocidad de cada punto puede escribirse como
![\mathbf{v}(\mathbf{r})=\mathbf{v}_C+\vec{\omega}\times\mathbf{r}'](/wiki/images/math/3/9/1/391b27bf1f46e9e3c504f8b53e2cbcc0.png)
Asimismo, esta distribución se encuentra en el seno de un campo electromagnético externo, de forma que cada elemento de carga se encuentra sometido a una fuerza
![\mathrm{d}\mathbf{F}=\rho(\mathbf{E}+\mathbf{v}\times\mathbf{B})\,\mathrm{d}\tau](/wiki/images/math/7/1/6/716b2f9bf12bc72a3f0ef5c0f6982d93.png)
2 Fuerza sobre la distribución
La fuerza neta sobre la distribución de carga será la resultante de las fuerzas diferenciales
![\mathbf{F}=\int \mathrm{d}\mathbf{F}=\int\rho(\mathbf{E}+\mathbf{v}\times\mathbf{B})\,\mathrm{d}\tau](/wiki/images/math/8/d/6/8d6d004fa5b53ecaa230145c5963ac61.png)
Veamos cada contribución por separado.
2.1 Fuerza eléctrica
La fuerza eléctrica sobre la distribución será
![\mathbf{F}_\mathrm{e}=\int \rho\mathbf{E}\,\mathrm{d}\tau](/wiki/images/math/7/f/4/7f4db366424d3f1397c3fb93ea496159.png)
Podemos calcular una expresión aproximada para esta fuerza aplicando que la distribución de carga está localizada en torno a su centro, de manera que podemos sustituir el campo eléctrico por su desarrollo en serie de Taylor en torno al centro de la distribución
![\mathbf{E}(\mathbf{r})=\mathbf{E}_C = \mathbf{r}'\cdot\mathbf{E}_C+\frac{1}{2}\mathbf{r}'\mathbf{r}':\nabla\nabla\mathbf{E}_C + \cdots](/wiki/images/math/0/6/8/06811385ac04ab91fbd464fb18f3096f.png)
donde ,
,
son constantes iguales al valor del campo y sus derivadas sucesivas en el centro de la distribución.
De esta forma la fuerza eléctrica viene dada por la serie
![\mathbf{F}_\mathrm{e}= \int \rho \mathbf{E}_C\,\mathrm{d}\tau+ \int \rho \mathbf{r}'\cdot\nabla\mathbf{E}_C\,\mathrm{d}\tau + \frac{1}{2}\int \rho \mathbf{r}'\mathbf{r}':\nabla\nabla\mathbf{E}_C\,\mathrm{d}\tau+\cdots](/wiki/images/math/4/2/4/424f212b0c7477fbec9542754d5432fa.png)
y, sacando las constantes fuera de las integrales
![\mathbf{F}_\mathrm{e}=Q\mathbf{E}_C+(\mathbf{p}\cdot\nabla)\mathbf{E}+\frac{1}{2}\mathsf{Q}:\nabla\nabla\mathbf{E}_C+\cdots](/wiki/images/math/e/c/7/ec73a450d4be85c708e6565523879d86.png)
donde
![Q=\int\rho\,\mathrm{d}\tau](/wiki/images/math/3/e/0/3e02473fbb303f0d7e55817ee7046c06.png)
![\mathbf{p}=\int\rho\,\mathbf{r}'\mathrm{d}\tau](/wiki/images/math/c/e/b/ceb1d122695e9736443c801da72f0c99.png)
![\mathsf{Q}=\int \rho\mathbf{r}'\mathbf{r}'\,\mathrm{d}\tau](/wiki/images/math/c/c/f/ccfcb9dde370ba08d09ea2ccbd8329f0.png)
Aquí Q es la carga neta de la distribución. es el momento dipolar, el cual se anula debido a la simetría de la distribución.
El término es un tensor. Por la simetría de la distribución, este tensor es realmente un escalar, ya que
![Q_{ik}=\int \rho x'_ix'_k\,\mathrm{d}\tau=0\qquad(i\neq k)](/wiki/images/math/1/c/9/1c9b8347b96d00b6b45bcaf96c0b010c.png)
y
![Q_{ii}=\int \rho {x'_i}^2\,\mathrm{d}\tau = \frac{1}{3}\int \rho r'^2 \,\mathrm{d}\tau = \gamma QR^2 \qquad (i=k)](/wiki/images/math/4/2/8/428c26d8127a13072598b18b0854de33.png)
aquí γ es un coeficiente que depende de la forma concreta de la distribución. Para una esfera cargada uniformemente en volumen vale 1/5, para una cargada en superficie es 1/3, para un decaimiento exponencial vale 4. Reuniendo los dos resultados
![\mathsf{Q}=\gamma Q R^2 \mathsf{I}\,](/wiki/images/math/4/7/a/47ad895d3d0442a9b440b00315b50924.png)
Llevando esto a la expresión de la fuerza, queda
![\mathbf{F}_e = Q\mathbf{E}_C + \frac{\gamma}{2}QR^2\mathsf{I}:\nabla\nabla\mathbf{E}+\cdots](/wiki/images/math/2/7/5/2751cf413e3878f4e9bc7965509ff7dd.png)
pero la contracción del tensor unidad con el gradiente del gradiente es simplemente el laplaciano
![\mathbf{F}_e = Q\mathbf{E}_C + \frac{\gamma}{2}QR^2\nabla^2\mathbf{E}+\cdots](/wiki/images/math/b/d/5/bd573cecdce0bd43b073e2e1b9ec0257.png)
Ahora, para un campo eléctrico externo estático y debido a cargas en posiciones no coincidentes con la distribución (como es el caso del campo de un núcleo sobre un electrón), el laplaciano del campo eléctrico se anula y la fuerza se reduce a
![\mathbf{F}=Q\mathbf{E}_C](/wiki/images/math/0/3/0/0305e90ba33b7c196b5e93d22a2982e9.png)
Hay que señalar que el término en QR2 no es un término cuadrupolar. El término cuadrupolar se refiere a la componente simétrica de traza nula. El término que sale aquí se debe a que el centro de fuerzas eléctricas no coincide con el centro de la distribución.
2.2 Fuerza magnética
La fuerza magnética sobre la distribución será
![\mathbf{F}_\mathrm{m}=\int \rho\mathbf{v}\times\mathbf{B}\,\mathrm{d}\tau](/wiki/images/math/a/1/1/a117d34235c9caf2724fdaafe8ca6649.png)
Sustituyendo la expresión de la velocidad
![\mathbf{F}_\mathrm{m}=\int\rho(\mathbf{v}_C+\vec{\omega}\times\mathbf{r}')\times\mathbf{B}\mathrm{d}\tau](/wiki/images/math/1/4/8/148e47d204e35faf72894d2a86ee22c7.png)
Para el primer término tenemos
![\mathbf{F}_{m1}=\int \rho\mathbf{v}_C\times\mathbf{B}\mathrm{d}\tau = \mathbf{v}_C\times\int\rho \mathbf{B}\mathrm{d}\tau](/wiki/images/math/9/5/c/95c692094b6e37a7c7a68ebbdbee7e36.png)
Aquí podemos hacer el mismo desarrollo que para el campo eléctrico y nos queda
![\mathbf{F}_{m1}=Q\mathbf{v}_C\times\mathbf{B}_C + \frac{\gamma QR^2}{2}\mathbf{v}_C\times(\nabla^2\mathbf{B})](/wiki/images/math/a/0/0/a00b0e6697c9ebb561d66eb5da375958.png)
De nuevo, si este campo magnético se debe a fuentes situadas fuera de la distribución (como el núcleo sobre un electrón), el laplaciano del campo magnético es nulo y este término se reduce a
![\mathbf{F}_{m1}=Q\mathbf{v}_C\times\mathbf{B}_C](/wiki/images/math/7/c/4/7c43799c97c98c6a91c2dbe1fb425a51.png)
Para el segundo término desarrollamos el doble producto vectorial
![\mathbf{F}_{m2}=\int \rho (\vec{\omega}\cdot\mathbf{B})\mathbf{r}'\mathrm{d}\tau-\int \rho \vec{\omega}(\mathbf{r}'\cdot\mathbf{B})\mathrm{d}\tau=\mathbf{F}_{m21}+\mathbf{F}_{m22}](/wiki/images/math/6/3/0/630f114181a245538fb0eee3ab305f11.png)
Aplicando ahora la serie de Taylor para el campo magnético
![\mathbf{B}=\mathbf{B}_C+\mathbf{r}'\cdot\nabla\mathbf{B}_C+\frac{1}{2}\mathbf{r}'\mathbf{r}':\nabla\nabla\mathbf{B}](/wiki/images/math/1/2/c/12c2a7c1db89c7ecca7afd4eb1acea97.png)
obtenemos en primer lugar
![\mathbf{F}_{m21}=\int \rho (\vec{\omega}\cdot\mathbf{B}_C)\mathbf{r}'\mathrm{d}\tau+\int_\rho(\vec{\omega}\cdot(\mathbf{r}'\cdot\nabla\mathbf{B}))\mathbf{r}'\mathrm{d}\tau](/wiki/images/math/1/2/2/1225117ead69952031fb77f7a3836065.png)
En primer término es nulo, por la simetría. Para el segundo, empleando subíndices
![F_{m21i}=\int \rho x_i\omega_j(x_k\partial_k B_j)\mathrm{d}\tau = \gamma QR^2 \delta_{ik} \omega_j\partial_k B_j = \gamma QR^2 \partial_i(\omega_jB_j)\,](/wiki/images/math/0/5/6/05646b1faab59e54e68394ad36612f08.png)
o, en forma vectorial
![\mathbf{F}_{m21} = \gamma QR^2 \nabla(\vec{\omega}\cdot\mathbf{B})](/wiki/images/math/f/5/a/f5af852f7a99dff958f806c5b219a97f.png)
Para tenemos, igualmente
![\mathbf{F}_{m22}=-\int \rho \vec{\omega}(\mathbf{r}'\cdot\mathbf{B})\mathrm{d}\tau =-\vec{\omega}\int \rho\mathbf{r}'\cdot(\mathbf{r}'\cdot\nabla)\mathbf{B}_C)\mathrm{d}\tau](/wiki/images/math/3/b/a/3ba6c45c59e994d2901ebdfc1e9a749e.png)
Empleando de nuevo subíndices
![F_{m22i} = -\omega_i \int \rho x'_jx'_k\partial_k B_j\mathrm{d}\tau = -\omega_i \gamma QR^2 \delta_{jk}\partial_kB_j=-\omega_i \gamma QR^2 \partial_jB_j\,](/wiki/images/math/5/1/1/511047bc28e82da2972f95667383aa12.png)
En forma vectorial
![\mathbf{F}_{m22}=-\vec{\omega}\gamma QR^2 \nabla\cdot\mathbf{B}](/wiki/images/math/4/b/6/4b6a3bdc744f19c8780720b7e4da0cfd.png)
pero la divergencia del campo magnético es siempre nula, por lo que este término se anula. Nos queda por tanto la fuerza magnética
![\mathbf{F}_m = Q\mathbf{v}_C\times\mathbf{B}_C+ \gamma QR^2\nabla(\vec{\omega}\cdot\mathbf{B})](/wiki/images/math/4/f/4/4f480227f1d2f2adea8ee65084c8cf13.png)
2.3 Fuerza total
Sumando las fuerzas eléctrica y magnética obtenemos la fuerza neta
![\mathbf{F}=Q(\mathbf{E}_C+\mathbf{v}_C\times\mathbf{B}_C)+\gamma QR^2 \nabla(\omega\cdot\mathbf{B})](/wiki/images/math/c/3/5/c35a89d2657776c08c5bef5f8117f1cb.png)
Si las fuentes del campo eléctrico y magnético coinciden con lso puntos de la distribución, habrá que añadir los términos correspondientes a los laplacianos de los campos.
3 Momento sobre la distribución
De manera análoga se halla el momento de las fuerzas respecto al centro de la distribución. Este momento será
![\mathbf{M}_C=\int \mathbf{r}'\times\mathrm{d}\mathbf{F}=\int\rho\mathbf{r}'\times(\mathbf{E}+\mathbf{v}\times\mathbf{B})\,\mathrm{d}\tau](/wiki/images/math/4/a/a/4aa9965c9b8414655a13a0d9bb5ba31e.png)
Como con la fuerza, consideraremos cada momento por separado.