Magnetohidrodinámica (MHD)

Confinamiento de plasmas para reactores nucleares por fusión.

I. I NTRODUCCI ÓN elementos pesados resultantes tienen una masa ligeramente

menor que los elementos de fusión y esta diferencia de masa
La búsqueda de una fuente de energı́a a gran escala y limpia, resulta en liberación de grandes cantidades de energı́a.
ha llevado a la creación de reactores nucleares; instalación Un conjunto de partı́culas con energı́a promedio de 5KeV
capaz de iniciar, controlar y mantener las reacciones nucleares tiene una temperatura efectiva de casi 40 millones Kelvin, a
en cadena La liberación de energı́a en los isotopos radiactivos una temperatura de unos cuantos miles de grados, el gas se
que ocurren de forma natural es demasiado lenta para ser una disocia; es decir, los electrones se desprenden de sus átomos,
fuente de energı́a de uso práctico. No obstante, a través del y dichos electrones junto con los núcleos separados se mueven
proceso de fisión nuclear es posible una liberación de energı́a aleatoriamente. A este conjunto de materia ionizada se le
mucho más rápida. conoce como plasma.
Pero la fisión no es la única reacción nuclear liberadora de
energı́a. La fusión de núclidos ligeros, como el hidrógeno, Las reacciones de fusión incontroladas en forma de ex-
en elementos pesados también es un proceso productor plosiones termonucleares que usan deuterio y tritio, donde el
de energı́a. El problema central en la búsqueda de energı́a plasma no está confinado, ya se han logrado en la Tierra. Por
de fusión controlada ha sido diseñar una configuración lo tanto, no hay duda de que la fusión termonuclear funciona y
de campo magnético que permita contener el plasma en es un reacción altamente exoenergética. [1]
un equilibrio mas o menos estable. Este trabajo se centra
en el estudio de problemas magnetohidrodinámicos en
dispositivos de confinamiento magnético tipo tokamak, los P LASMA COMO FLUIDO MHD.
cuales tienen forma tipo toroidal. La investigación de la fusión
Para entender algunas de las caracterı́sticas del plasma, resulta
nuclear controlada a sido estudiada por más de 50 años, aún util su descripción como un fluido, en la que se consideran al-
encontrándose en experimentación. gunos aspectos macroscospicos sin tener en cuenta de manera
detallada la distribución de velocidades microscospicas. Ası́
En los experimentos de confinamiento magnético orientados a un plasma se puede describir como dos fluidos magnetohidro-
fusión nuclear se busca mantener un plasma en las condiciones dinámicos de iones y electrones con densidades de masa ρmi
adecuadas para que se produzcan suficientes reacciones de ρme , densidad de carga ρ, densidad de corriente j, velocidades
fusión de forma controlada. En tales condiciones el plasma se de flujo Vi Ve , y presiones pi , pe ; las cuales pueden ser expre-
encuentra totalmente ionizado y se puede modelar como un sadas en función de promedios en el espacio de velocidades,
fluido conductor. El marco teórico para describir el comporta- usando la función de distribución en el espacio de velocidades
miento macroscópico de un plasma de esta naturaleza es la fi (r, v,t), el número de densidad de iones se expresa como [2]:
magnetohidrodinámica (MHD). La magnetohidrodinámica
estudia la dinámica de fluidos conductores de electricidad en
Z
ni (r,t) = fi (r, v,t)dv (1)
presencia de campos eléctricos y magnéticos; especı́ficamente
efectos que aparecen por la interacción entre el movimiento
El transporte en plasmas es uno de los temas más importantes
del fluido y un campo magnético cualquiera que este presente.
en la investigación de fusión. La ecuación de transporte para
. Algunos ejemplos incluyen plasmas, metales lı́quidos y
particulas es:
el agua salada. Estudiando con la MHD plasmas calientes,
tales como los plasmas en astrofı́sica y termonucleares de ∂
origen natural (estrellas) o artificial (tokamaks). En el caso n(r,t) + 5 · [n(r,t)V (r,t)] = 0 (2)
∂t
ideal, el plasma se trata como un fluido homogéneo, como
un conductor perfecto, por lo que posee una conductividad Un plasma consiste en muchas partı́culas cargadas en movi-
eléctrica infinita y con una viscosidad nula. miento y tiene muchos grados de libertad magnetohidrodinámi-
cos, ası́ como grados de libertad en el espacio de velocidad. El
calentar un plasma aumenta la energı́a cinética de las partı́culas
R EACCIONES DE FUSI ÓN Y PLASMA .
cargadas, pero al mismo tiempo puede inducir fluctuaciones
En las reacciones de fusión dos núcleos ligeros chocan en los campos eléctrico y magnético, que a su vez aumentan la
aproximándose lo suficiente para que entre en juego la fuerza difusión anómala; siendo muy importante determinar si algún
nuclear, fundiéndose para formar un núcleo más pesado. Los modo perturbado en particular es estable o inestable . Si la

1
Magnetohidrodinámica MHD

desviación del estado de equilibrio es pequeña, se puede uti- al equilibrio termodinámico, la inestabilidad se desarrolla.
lizar una aproximación linealizada, ası́ estas inestabilidades Una configuración de confinamiento será estable cuando
pueden ser descritas por ecuaciones magnetohidrodinámicas cualquier pequeña deformación de la configuración del plasma
lineales. [3] y del campo magnético produzcan fuerzas restauradoras que
eliminen la deformación, conduciendo al sistema al estado
C ONFINAMIENTO MAGN ÉTICO . inicial no perturbado.

En la aproximación magnetohidrodinámica la condición de Tomando la fig.1 con la región del plasma en x < 0 y
equilibrio magnetohidrostático es: la región del vacı́o en x > 0, la aceleración g se aplica en
la dirección x. Los iones y los electrones se desplazan en
~ 1 ~ ~
~
5P = (J ∧ B) (3) direcciones opuestas entre sı́, debido a la aceleración, con
c velocidades de deriva, y la separación de carga debida a la
deriva de iones y electrones opuestos produce un campo
donde J~ es la densidad de corriente en el plasma y ~B la in-
eléctrico. El flujo magnético originalmente dentro pero cerca
ducción magnética total, que incluye al campo externo y el
del lı́mite del plasma se desplaza de manera que está fuera del
generado por el propio plasma. La ecuación 3 es la base del
lı́mite, mientras que el flujo externo se mueve para llenar la
confinamiento magnético de plasmas, expresando que la fuerza
depresión que queda en el lı́mite; este tipo de inestabilidad
ejercida por el plasma para expandirse puede ser expresada
se conoce como inestabilidad de intercambio. Fenómenos
por la fuerza de Lorentz y el plasma puede ser confinado por
similares ocurren en la hidrodinámica cuando un fluido
el campo magnético.
denso es soportado contra la gravedad por un fluido de
La complejidad se manifiesta en que en la ecuación 3, J~ y ~B
menor densidad. Por lo tanto, la inestabilidad del intercambio
no son independientes (~B es la suma del campo externo más
~ Reescribiendo 3, considerando un plasma también se llama inestabilidad de Rayleigh-Taylor, y es la
el generado por J).
principal inestabilidad que presenta el confinamiento inercial.
encerrado en una cavidad cilı́ndrica, donde el eje del cilindro y
campo magnético están orientados en el eje z:

B2
5⊥ (P + )=0 (4)
8π
Con 5⊥ el gradiente en dirección transversa al campo. Tal
condición de equilibrio resulta equivalente a la condición de
equilibrio hidrostática, con la excepción que la presión consta
de dos términos (la presión mecánica P y la magnética B2 /8π).

C ONFINAMIENTO I NERCIAL
Esta alternativa consiste en alojar la mezcla gaseosa en una
pequeña bola o cuenta esférica, la cual es bombardeada desde Figura 1: Derivas de iones y electrones, y campo eléctrico
todas direcciones por haces de rayos láser de alta potencia, resultante por el intercambio de inestabilidad
depositando una gran cantidad de energı́a. La inercia mantie-
ne confinado el plasma tiempo suficiente para que ocurra la Los problemas de estabilidad del plasma pueden estudiarse con
reacción de fusión aproximaciones infinitesimales de las perturbaciones desde el
estado de equilibrio. De la ley de Ohm y las ecuaciones de
I NESTABILIDADES Maxwell se puede obtener:
 
Existe una variedad de inestabilidades en el modelo MHD 4πσ 1 1
5 ×(5 × B) = − B + 5 ×(V × B) (5)
ideal que pueden afectar el confinamiento de un plasma, c ct c
inestabilidades que no desaparecen si se expande el modelo
para considerar efectos microscópicos pero el modelo MHD o
c2
ideal sirve como punto de partida para el estudio de todo ∂t B = 5 × (V × B) + 52 B (6)
dispositivo de confinamiento de plasmas. 4πσ
El exceso de energı́a libre que tiene un plasma confinado Por lo que el campo magnético en un punto en el plasma se
respecto a uno que no lo está desencadena inestabilidades; el puede cambiar por la convección de fluidos.
equilibrio termodinámico corresponde al mı́nimo absoluto de
energı́a libre. Las inestabilidades magnetohidrodinámicas se
S ISTEMAS T OROIDALES
deben a que un plasma confinado puede convertir parte de su
energı́a interna en movimiento organizado a través del campo En un sistema toroidal las lı́neas de campo magnético se cie-
magnético; si esto disminuye la energı́a del sistema y lo acerca rran sobre si mismas guiando a las partı́culas cargadas sobre

2
Magnetohidrodinámica MHD

trayectorias dentro del toro, pero un campo toroidal generado T OKAMAK .

por solo una componente Bφ es inadecuado ya que el plasma
El tokamak es una máquina experimental diseñada para aprove-
es desplazado hacia la pared exterior del toro, destruyendo el
char la energı́a de la fusión, donde el plasma tiene un confina-
confinamiento. La solución consiste en modificar el campo
miento toroidal. El campo toroidal utilizado para estabilizarse
toroidal a modo de obtener lı́neas de fuerza helicoidales, de tal
contra las inestabilidades de MHD es lo suficientemente fuerte
forma que el desplazamiento de cargas en dirección vertical
como para satisfacer la condición de Kruskal-Shafranov.
sea virtualmente anulado; esto se consigue añadiendo además
del campo toroidal, una componente poloidal.

Figura 3: Esquema del sistema de campo magnético en un

tokamak.

Las lineas de campo magnético en geometrı́as tan cerradas se

ven forzadas a seguir superficies toroidales, y las partı́culas
de plasma inicialmente se agrupan en espiral a lo largo de las
lineas de campo; no obstante los campos toroidales simples
no pueden confinar el plasma en equilibrio, por lo que se debe
Figura 2: Geometrı́a toroidal emplear un esquema estabilizador. Actualmente este tipo de
reactor es la base del proyecto internacional ITER.

Figura 4: Tokamak [4]

3
Magnetohidrodinámica MHD

E CUACI ÓN DE G RAD -S HAFRANOV R EFERENCIAS

El equilibrio ideal del tokamak se describe mediante soluciones [1] M. Kikuchi, Fusion Physics. International Atomic Energy
asimétricas al conjunto de las siguientes ecuaciones: Agency, 2012.
5 ×~B = µ µ0~j (7)
[2] D. R. Nicholson, Introduction to Plasma Theory. Wiley
5 ·~B = 0 (8) series in Plasma Physics, 1983.
5 p = ~j × ~B (9) [3] K. Miyamoto, Plasma Physics and Controlled Nuclear
Resulta conveniente dividir los campos vectoriales de densidad Fusion. Springer, 3 ed., 2005.
de corriente y campo magnético en partes poloidales y toroida-
les; ~j = ~j pol − ~jtor , ~B = ~B pol + ~Btor , expresando las funciones [4]
de flujo en termino de componentes poloidales, en coordenadas
cilindricas R, φ , Z. (ver Fig. 2)

~j pol = 1 ~eφ × 5F (10)

R

~B pol = − 1 ~eφ × 5ψ (11)

R
Ahora, el producto punto de 9 con 5ψ da la parte no trivial de
la fuerza de balance:
dP
5ψ · 5P = (5ψ)2
dψ
= 5ψ · (~j pol × ~Btor ) + 5(~jtor × ~B pol ) (12)
 
µ0 1
= (5ψ)2 − 2 F 0 + jφ
R R
la cual constituye una ecuación para la densidad de corriente
toroidal ~jφ . Colocando esto dentro de la componene toroidal
de la ley de Ampere, se puede llegar a una ecuación diferencial
parcial semilineal, para el flujo de la función ψ [1].

∂ 2ψ 1 ∂ ψ ∂ 2ψ 2 dP µ02 dF 2
− + = −µ 0 R jφ = −µ0 R −
∂ R2 R ∂ R ∂ Z 2 ψ 2 fψ
(13)
en equilibrio, la presión del plasma y la función de campo
toroidal dependen de ψ. El equilibrio MHD del plasma queda
determinado en términos de la ecuación de Grad- Shafranov
una vez queden especificadas P y F, una vez especificadas se
deben imponer las condiciones de contorno en ψ para incluir
efecto de las corrientes que circulan por las bobinas externas.

II. C ONCLUSIONES
La búsqueda de reactores nucleares por fusión nuclear hoy
en dı́a se encuentra en investigación, y como se muestra el
problema principal para la creación de estos es la forma de
confinar el plasma; si bien existen diferentes formas de con-
finarlo como lo es la inercial y la magnética, en ambas se
utiliza en lo macroscópico el marco teórico de la magnetohi-
drodinámica, la cual nos proporciona para el caso la ecuación
de Grad-Shafranov, la cual nos dará la forma espacial del plas-
ma y la configuración magnética en equilibrio. Estos trabajos
de investigación siguen aún en proceso existiendo diversas
asociaciones para llevar a cabo esta misión y la mayorı́a se
centra en el confinamiento magnético, como es el caso del
ITER que trabaja en la creación de un tokamak.