Entrar Página Discusión Historial Go to the site toolbox

Desarrollo Multipolar magnético

De Laplace

(Diferencias entre revisiones)
(Estimación del error cometido)
(Eliminado todo el contenido de la página)
 
(Una edición intermedia no se muestra.)
Línea 1: Línea 1:
-
__TOC__
 
-
==Introducción==
 
-
Supongamos que tenemos una distribución de corriente estacionaria que ocupa una pequeña región del espacio y queremos hallar el campo en puntos alejados.
 
-
Cuando se dice “una pequeña región del espacio” se entiende que comparada con la distancia al punto de observación. Matemáticamente:
 
-
 
-
<center><math>\delta=\frac{\mathrm{max}(|\mathbf{r}'|)}{r}\ll 1</math></center>
 
-
 
-
Como con el [[Desarrollo multipolar eléctrico|campo eléctrico]], la idea del desarrollo multipolar es hacer un cálculo aproximado, más sencillo que la integral exacta (la cual puede ser imposible de calcular) mediante el empleo de una serie de Taylor. Partimos de la expresión del [[Potencial vector magnético|potencial vector]] para el caso de una espira
 
-
 
-
<center><math>\mathbf{A}=\frac{\mu_0I}{4\pi}\oint\frac{\mathrm{d}\mathbf{r}'}{|\mathbf{r}-\mathbf{r}'|}</math></center>
 
-
 
-
Aplicando el desarrollo del binomio de Newton
 
-
 
-
<center><math>\frac{1}{|\mathbf{r}-\mathbf{r}'|}=\frac{1}{r}+\frac{\mathbf{r}\cdot\mathbf{r}'}{r^3}+\mathrm{O}(\delta^2)</math></center>
 
-
 
-
resulta (de forma no trivial) la expresión aproximada para el potencial vector
 
-
 
-
<center><math>\mathbf{A}=\frac{\mu_0}{4\pi}\frac{\mathbf{m}\times\mathbf{r}}{r^3}+\cdots</math></center>
 
-
 
-
donde
 
-
 
-
<center><math>\mathbf{m}=\frac{I}{2}\oint \mathbf{r}'\times\mathrm{d}\mathbf{r}'</math></center>
 
-
 
-
es el denominado ''momento dipolar magnético'' de la espira.
 
-
 
-
A una distribución de corriente que produce este potencial vector (y el campo magnético resultante) se le denomina ''[[dipolo magnético]]''. Podemos imaginar un dipolo magnético como una pequeña espira de corriente, si bien el concepto es más general:
 
-
 
-
* Una distribución de corriente de volumen confinada a una pequeña región del espacio también produce el campo de un dipolo magnético, siendo su momento dipolar
 
-
 
-
<center><math>\mathbf{m}= \frac{1}{2}\int \mathbf{r}'\times\mathbf{J}\,\mathrm{d}\tau'</math></center>
 
-
 
-
:En particular, este resultado incluye el caso de una distribución de carga en rotación.
 
-
 
-
* Lo mismo ocurre para una distribución de corriente superficial
 
-
 
-
<center><math>\mathbf{m}= \frac{1}{2}\int \mathbf{r}'\times\mathbf{K}\,\mathrm{d}S'</math></center>
 
-
 
-
* La partículas elementales (electrones, protones,&hellip;) producen un campo magnético dipolar, en el que el momento magnético es proporcional al llamado espín de la partícula. Este momento dipolar no está asociado a ninguna corriente &ldquo;clásica&rdquo;, pero su comportamiento es análogo.
 
-
 
-
===Demostración detallada===
 
-
La obtención de la expresión aproximada del potencial vector no es elemental.
 
-
 
-
Para empezar hay que justificar el desarrollo de <math>1/|\mathbf{r}-\mathbf{r}'|</math>. Para ello observamos que se puede escribir como
 
-
 
-
<center><math>\frac{1}{|\mathbf{r}-\mathbf{r}'|}=\frac{1}{\sqrt{(\mathbf{r}-\mathbf{r}')\cdot(\mathbf{r}-\mathbf{r}')}}=\left(r^2-2\mathbf{r}\cdot\mathbf{r}'+r^{'2}\right)^{-1/2}</math></center>
 
-
 
-
y sacando factor común un <math>r^2</math> de la raíz
 
-
 
-
<center><math>\frac{1}{|\mathbf{r}-\mathbf{r}'|}=(r^2)^{-1/2}\left(1-2\frac{\mathbf{r}\cdot\mathbf{r}'}{r^2}+\frac{r^{'2}}{r^2}\right)^{-1/2} = \frac{1}{r}\left(1+\left(-2\frac{\mathbf{r}\cdot\mathbf{r}'}{r^2}+\frac{r^{'2}}{r^2}\right)\right)^{-1/2}</math></center>
 
-
 
-
Hasta aquí no hay aproximación alguna. Observamos que en el último factor tenemos 1 más algo mucho más pequeño que la unidad (pues <math>r'\ll r</math>). La fórmula general del binomio de Newton nos dice que si <math>x\ll 1</math>
 
-
 
-
<center><math>(1+x)^n = 1 + n x + \mathrm{O}(x^2)\,</math></center>
 
-
 
-
Aplicando esto al resultado anterior
 
-
 
-
<center><math>\frac{1}{|\mathbf{r}-\mathbf{r}'|} = \frac{1}{r}\left(1+\left(-2\frac{\mathbf{r}\cdot\mathbf{r}'}{r^2}+\frac{r^{'2}}{r^2}\right)\right)^{-1/2} = \frac{1}{r}\left(1-\frac{1}{2}\left(-2\frac{\mathbf{r}\cdot\mathbf{r}'}{r^2}+\frac{r^{'2}}{r^2}\right)+\mathrm{O}(\delta^2)\right)</math></center>
 
-
 
-
pero de hecho, el segundo de los dos sumandos del paréntesis también es de orden <math>\delta^2</math>, por lo que podemos despreciarlo y reducir el desarrollo a
 
-
 
-
<center><math>\frac{1}{|\mathbf{r}-\mathbf{r}'|} = \frac{1}{r}\left(1+\frac{\mathbf{r}\cdot\mathbf{r}'}{r^2}+\mathrm{O}(\delta^2)\right)=
 
-
\frac{1}{r}+\frac{\mathbf{r}\cdot\mathbf{r}'}{r^3}+\cdots</math></center>
 
-
 
-
El segundo paso es sustituir esto en la expresión del potencial vector. Nos queda
 
-
 
-
<center><math>\mathbf{A}=\frac{\mu_0I}{4\pi}\left(\oint\frac{\mathrm{d}\mathbf{r}'}{r}+\oint\frac{(\mathbf{r}\cdot\mathrm{r}')\mathrm{d}\mathbf{r}'}{r^3}+\cdots\right)</math></center>
 
-
 
-
Estas integrales son sobre la variable <math>\mathbf{r}'</math>, así que <math>\mathbf{r}</math> es una constante en ellas y puede ser extraído en la medida de lo posible
 
-
 
-
<center><math>\mathbf{A}=\frac{\mu_0I}{4\pi}\left(\frac{1}{r}\oint\mathrm{d}\mathbf{r}'+\frac{1}{r^3}\oint(\mathbf{r}\cdot\mathrm{r}')\mathrm{d}\mathbf{r}'+\cdots\right)</math></center>
 
-
 
-
La primera de estas dos integrales es el desplazamiento neto al recorrer una curva cerrada, por lo que se anula identicamente,
 
-
 
-
<center><math>\oint \mathrm{d}\mathbf{r}'=\left.\mathbf{r}'\right|_{\mathbf{r}_0}^{\mathbf{r}_0}=\mathbf{0}</math></center>
 
-
 
-
Lo que nos dice este resultado, de nuevo, es que el campo magnético de corrientes estacionarias, no posee término monopolar, esto es, que el campo de corrientes no equivale al campo de cargas magnéticas (''[[monopolo magnético|monopolos]]'').
 
-
 
-
El segundo término no se puede integrar de forma inmediata. Para conseguir extraer <math>\mathbf{r}</math> de la integral podemos usar el [[Generalización del teorema de Stokes|teorema vectorial]], generalización del teorema de Stokes
 
-
 
-
<center><math>\oint \mathrm{d}\mathbf{r}' \phi = \int \mathrm{d}\mathbf{S}\times\nabla\phi</math></center>
 
-
 
-
que, en nuestro caso da
 
-
 
-
<center><math>\oint \mathrm{d}\mathbf{r}'(\mathbf{r}\cdot\mathrm{r}') = \int \mathrm{d}\mathbf{S}\times\nabla'(\mathbf{r}\cdot\mathrm{r}')</math></center>
 
-
 
-
pero, como se [[Algunas identidades vectoriales|puede demostrar]]
 
-
 
-
<center><math>\nabla'(\mathbf{r}\cdot\mathbf{r}')=\mathbf{r}</math></center>
 
-
 
-
por lo que la integral se convierte en
 
-
 
-
<center><math>\oint \mathrm{d}\mathbf{r}'(\mathbf{r}\cdot\mathrm{r}') = \int \mathrm{d}\mathbf{S}\times\mathbf{r}=\mathbf{S}\times\mathbf{r}</math></center>
 
-
 
-
y, aplicación la expresión del [[vector superficie]]
 
-
 
-
<center><math>\mathbf{S}=\frac{1}{2}\oint \mathbf{r}'\times \mathrm{d}\mathbf{r}'</math></center>
 
-
 
-
nos queda finalmente
 
-
 
-
<center><math>\mathbf{A}=\frac{\mu_0}{4\pi}\,\frac{\mathrm{m}\times\mathbf{r}}{r^3}+\cdots</math>{{qquad}}{{qquad}}<math>\mathbf{m}=I\mathbf{S}=\frac{I}{2}\oint \mathbf{r}'\times \mathrm{d}\mathbf{r}'</math></center>
 
-
 
-
La extensión al caso de una distribución de corriente volumétrica o superficial requiere aun más cálculo vectorial, por lo que nos limitaremos a señalar que puede obtenerse siguiento las reglas usuales de transformación
 
-
 
-
<center><math>\int_\Gamma(\ldots)I\mathrm{d}\mathbf{r}'\to \int_S(\ldots)\mathbf{K}\mathrm{d}S'\to\int_\tau(\ldots)\mathbf{J}\mathrm{d}\tau'</math></center>
 
-
 
-
==Ejemplos==
 
-
===Espira plana===
 
-
Para una espira plana, el [[Vector_superficie#Caso_de_una_curva_plana|vector superficie]] es
 
-
 
-
<center><math>\mathbf{S}=S\mathbf{n}</math></center>
 
-
 
-
Siendo <math>S</math> el área de la porción de plano limitada por la curva y n la normal a éste (con el sentido asignado por la regla de la mano derecha.
 
-
 
-
Para una espira circular
 
-
 
-
<center><math>\mathbf{m}=I\pi R^2\,\mathbf{n}</math></center>
 
-
 
-
===Espira alabeada===
 
-
Para una espira no contenida en un plano el [[Vector_superficie#Caso_de_una_curva_alabeada_2|vector superficie]] es
 
-
 
-
<center><math>\mathbf{S}=S_x\mathbf{u}_x+S_y\mathbf{u}_y+S_z\mathbf{u}_z</math></center>
 
-
 
-
con <math>S_i</math> las áreas de las proyecciones sobre los planos coordenados.
 
-
 
-
Como ejemplo, consideremos una espira alabeada formada por tres cuadrantes de circunferencia, de radio <math>R</math> que van de <math>R\mathbf{u}_x</math> a <math>R\mathbf{u}_y</math>, de ahí a <math>R\mathbf{u}_z</math>, y de vuelta a <math>R\mathbf{u}_x</math>. Tal como se ve en artículo sobre el [[vector superficie]], para esta superficie tenemos
 
-
 
-
<center><math>\mathbf{m}=I\mathbf{S}=\frac{I\pi R^2}{4}(\mathbf{u}_x+\mathbf{u}_y+\mathbf{u}_z)</math></center>
 
-
 
-
===Esfera cargada rotante===
 
-
{{ac|Momento magnético de una esfera en rotación}}
 
-
Un caso particular de corriente de volumen localizada en el espacio es el de una distribución de carga en rotación. Este modelo puede ayudar a describir el comportamiento de una partícula cargada, como un electrón, que se caracteriza por un momento angular (el espín), equivalente en ciertos aspectos a una rotación.
 
-
 
-
Supongamos una esfera de radio <math>R</math>, con una carga <math>q</math> distribuida uniformemente en el volumeny que gira en torno a un eje (que tomaremos como eje Z) con velocidad angular <math>\mathbf{w}</math>. Para este sistema la densidad de corriente es, en esféricas
 
-
 
-
<center><math>\mathbf{J}=\rho_0\mathbf{v}=\rho_0\mathbf{w}\times\mathbf{r}=\rho_0\omega r\,\mathrm{sen}\,\theta\,\mathbf{u}_\varphi</math>{{qquad}}{{qquad}}<math>\rho_0=\frac{q}{4\pi R^3/3}</math></center>
 
-
 
-
llevando esto a la integral resulta el momento magnético
 
-
 
-
<center><math>\mathbf{m}=\frac{q R^2}{5}\mathbf{w}</math></center>
 
-
 
-
==Campo magnético==
 
-
Al igual que el potencial vector, el campo magnético posee un desarrollo multipolar , que se puede obtener hallando el rotacional de cada sumando del potencial vector
 
-
 
-
<center><math>\mathbf{B}=\nabla\times\mathbf{A}=\nabla\times\left(\frac{\mu_0}{4\pi}\,\frac{\mathbf{m}\times\mathbf{r}}{r^3}\right)+\cdots</math></center>
 
-
 
-
Hallando el [[Dipolo_magnético#Demostración|rotacional del término dipolar]] queda
 
-
 
-
<center><math>\mathbf{B}=\frac{\mu_0}{4\pi}\,\frac{3(\mathbf{m}\cdot\mathbf{r})\mathbf{r}-r^2\mathbf{m}}{r^5}+\cdots</math></center>
 
-
 
-
==Estimación del error cometido==
 
-
El error cometido al aproximar el campo magnético exacto por su desarrollo multipolar puede estimarse a partir del primer término despreciado. El error relativo suele ser de orden <math>\delta</math> o más pequeño, esto es, que si la distancia del punto de observación a la distribución es 5 veces el tamaño de esta distribución, el error será del orden de 1/5 = 20%.
 
-
 
-
Por ejemplo, para una espira circular el [[Campo magnético de una espira circular|campo en los puntos de su eje]] es
 
-
 
-
<center><math>\mathbf{B}(z)=\frac{\mu_0I R^2\mathbf{u}_z}{2(R^2+z^2)^{3/2}}</math></center>
 
-
 
-
A una altura <math>z=5R</math> este campo vale
 
-
 
-
<center><math>\mathbf{B}(5R)=\frac{\mu_0I R^2\mathbf{u}_z}{2(R^2+25R^2)^{3/2}}=\frac{1}{52\sqrt{26}}\,\frac{\mu_0I}{R}\simeq 0.003771 \frac{\mu_0I\mathbf{u}_z}{R}</math></center>
 
-
 
-
Empleando el desarrollo multipolar, la espira puede aproximarse por un dipolo de momento
 
-
<center>
 
-
<math>\mathbf{m}=I\pi R^2\mathbf{u}_z</math></center>
 
-
 
-
y el campo producido por él en el eje
 
-
 
-
<center><math>\mathbf{B}\simeq\frac{\mu_0}{4\pi}\,\frac{3(\mathbf{m}\cdot\mathbf{r})\mathbf{r}-r^2\mathbf{m}}{r^5}=\frac{\mu_0IR^2\mathbf{u}_z}{2z^3}</math></center>
 
-
 
-
que a una altura <math>z = 5R</math> vale
 
-
 
-
<center><math>\mathbf{B}= \frac{1}{250}\,\frac{\mu_0I\mathbf{u}_z}{R}=0.004\,\frac{\mu_0I\mathbf{u}_z}{R}</math></center>
 
-
 
-
El error relativo cometido es
 
-
 
-
<center><math>\varepsilon = \frac{|\mathbf{B}_\mathrm{apr}-\mathbf{B}_\mathrm{ex}|}{|\mathbf{B}_\mathrm{ex}|}= \frac{52\sqrt{26}-250}{250}=6\%</math></center>
 
-
 
-
vemos que el error es de solo el 6%, en vez del 20% previsto. La razón es que por la simetría del sistema el siguiente término se anula y hay que ir a un orden superior, quedando un error del orden de <math>\delta^2</math>.
 
-
 
-
Además de proporcionarnos un resultado con un error pequeño, el desarrollo multipolar nos permite hallar el valor del campo para todos los puntos del espacio, y no solo para los puntos del eje.
 
-
 
-
==Momentos de orden superior==
 
-
Como en el desarrollo multipolar eléctrico, existen situaciones en las que el momento dipolar se anula y es preciso ir a órdenes superiores del desarrollo. Así tenemos el cuadrupolo magnético, el octupolo, etc.
 
-
 
-
Un ejemplo de sistema cuadrupolar lo forman dos espiras muy próximas por las cuales circulan corrientes en sentidos opuestos. En este caso, los campos dipolares se anulan mutuamente, por lo que es necesario recurrir a un orden superior del desarrollo.
 
-
[[Categoría:Campo magnético de corrientes estacionarias]]
 

última version al 18:23 4 abr 2009

Herramientas:

Herramientas personales
TOOLBOX
LANGUAGES
licencia de Creative Commons
Esta página fue modificada por última vez el 18:23, 4 abr 2009. - Esta página ha sido visitada 1.225 veces. - Aviso legal - Acerca de Laplace