Física de Plasmas y Fusión Termonuclear por Confinamiento Magnético: Termodinámica y Estabilidad en el Reactor Tokamak ITER

1. Cinemática Nuclear y Sección Eficaz de la Reacción Deuterio-Tritio

La fusión nuclear se produce cuando dos núcleos ligeros superan la repulsión electrostática colombiana hasta alcanzar distancias del orden del radio nuclear característico ($r_0 \approx 10^{-15} \text{ m}$ o $1 \text{ fm}$), donde la interacción nuclear fuerte, de carácter atractivo y corto alcance, se vuelve dominante.

Entre todas las posibles reacciones exotérmicas entre isótopos de hidrógeno y helio, la que presenta la mayor sección eficaz de dispersión inelástica $(\sigma)$ a temperaturas cinéticas alcanzables industrialmente ($T \sim 10 - 20 \text{ keV}$, equivalentes a $100 - 200 \times 10^6 \text{ K}$) es la fusión de Deuterio ($^2_1\text{H}$ o $\text{D}$) y Tritio ($^3_1\text{H}$ o $\text{T}$):

$^2_1\text{D} + \,^3_1\text{T} \longrightarrow \,^4_2\text{He}\,(3.517\text{ MeV}) + \,^1_0\text{n}\,(14.069\text{ MeV}) + 17.586\text{ MeV}$

1.1 Conservación del Momento Lineal y Reparto de Energía Cinética

El balance energético neto $Q$ liberado por la desintegración de masa se deduce de la relación relativista masa-energía de Einstein:

$Q = \left( m_{\text{D}} + m_{\text{T}} - m_{^4\text{He}} - m_{\text{n}} \right) c^2 = 17.586\text{ MeV}$

En el sistema del centro de masas, asumiendo una velocidad inicial de colisión prácticamente nula en comparación con las velocidades de eyección de los productos relativistas ($v_{\text{rel}} \ll c$), la conservación del momento lineal impone:

$\mathbf{p}_{\alpha} + \mathbf{p}_{\text{n}} = 0 \implies |\mathbf{p}_{\alpha}| = |\mathbf{p}_{\text{n}}| = p$

Dado que la energía cinética no relativista viene dada por $E_k = \frac{p^2}{2m}$:

$E_{\alpha} = Q \left( \frac{m_{\text{n}}}{m_{\alpha} + m_{\text{n}}} \right) \approx 17.586 \times \frac{1.00866}{4.0015 + 1.00866} \approx 3.52\text{ MeV} \quad (20\%)$

$E_{\text{n}} = Q \left( \frac{m_{\alpha}}{m_{\alpha} + m_{\text{n}}} \right) \approx 17.586 \times \frac{4.0015}{4.0015 + 1.00866} \approx 14.07\text{ MeV} \quad (80\%)$

Esta asimetría en el reparto de energía cinética es el fundamento de la arquitectura de ITER:

1. La partícula alfa ($^4\text{He}^{2+}$) posee carga eléctrica positiva $q = +2e$, lo que le permite permanecer atrapada en las líneas helicoidales del campo magnético del Tokamak, cediendo su energía térmica al plasma circundante mediante colisiones colombianas de frenado (calentamiento alfa intrínseco).
2. El neutrón ($^1_0\text{n}$), al ser eléctricamente neutro, escapa libremente de la botella magnética, colisionando con el manto fértil (breeding blanket) de la pared de la vasija de vacío, transfiriendo su energía cinética en forma de calor térmico a un refrigerante de agua a alta presión y transmutando átomos de litio en nuevo tritio combustible.

---

2. Efecto Túnel Cuántico y Pico de Gamow

La barrera de potencial culombiana entre dos núcleos de cargas $Z_1 e$ y $Z_2 e$ a una distancia $r$ viene dada por:

$V_C(r) = \frac{1}{4\pi\varepsilon_0} \frac{Z_1 Z_2 e^2}{r}$

Para el deuterio y el tritio ($Z_1 = Z_2 = 1$) en el radio de contacto nuclear $r_0 \approx 3 \times 10^{-15} \text{ m}$, la energía electrostática requerida para una aproximación clásica es:

$V_C(r_0) \approx \frac{(1.602 \times 10^{-19}\text{ C})^2}{4\pi (8.854 \times 10^{-12}\text{ F/m})(3 \times 10^{-15}\text{ m})} \approx 7.68 \times 10^{-14}\text{ J} \approx 480\text{ keV}$

Si el fenómeno dependiera exclusivamente de la física clásica de colisión, se requeriría una temperatura cinética de plasma de:

$T = \frac{2}{3 k_B} V_C \approx 3.7 \times 10^9\text{ K}$

A temperaturas cinéticas alcanzables en ITER ($T \sim 15\text{ keV} \sim 1.7 \times 10^8\text{ K}$), la penetración de la barrera de Coulomb se produce casi en su totalidad gracias al efecto túnel cuántico mecanocuántico, descrito por la aproximación WKB (Wentzel-Kramers-Brillouin).

2.1 Derivación de la Probabilidad de Penetración y Pico de Gamow

La sección eficaz microscópica nuclear $\sigma(E)$ para una energía de colisión en el centro de masas $E$ se parametriza como:

$\sigma(E) = \frac{S(E)}{E} \exp\left( -2\pi \eta \right)$

donde $S(E)$ es el factor astrofísico que varía suavemente con la energía, y $\eta$ es el parámetro adimensional de Sommerfeld:

$\eta = \frac{Z_1 Z_2 e^2}{4\pi\varepsilon_0 \hbar v} = \sqrt{\frac{E_G}{E}}, \quad E_G = 2 m_{\text{red}} c^2 (\pi \alpha Z_1 Z_2)^2$

donde $E_G$ es la energía de Gamow y $m_{\text{red}} = \frac{m_{\text{D}} m_{\text{T}}}{m_{\text{D}} + m_{\text{T}}}$ es la masa reducida del par colisionante.

Para calcular la reactividad volumétrica $\langle \sigma v \rangle$ en un plasma termalizado con distribución de velocidades de Maxwell-Boltzmann a temperatura $T$:

$\langle \sigma v \rangle = \left( \frac{8}{\pi m_{\text{red}} (k_B T)^3} \right)^{1/2} \int_0^{\infty} S(E) \exp\left( -\frac{E}{k_B T} - \sqrt{\frac{E_G}{E}} \right) dE$

El integrando representa el producto de dos funciones opuestas:

• La cola decreciente de alta energía de Maxwell-Boltzmann: $\exp(-E/k_B T)$.
• La función de transmisión túnel creciente: $\exp(-\sqrt{E_G/E})$.

El máximo del integrando define la ventana o pico de Gamow a la energía óptima de fusión:

$E_0 = \left( \frac{1}{2} \sqrt{E_G} k_B T \right)^{2/3}$

Para $T = 15 \text{ keV}$, $E_0 \approx 66 \text{ keV}$. Esto demuestra cuantitativamente que las reacciones de fusión en ITER no ocurren a la energía térmica media de las partículas, sino en una banda de energía significativamente mayor en la cola maxwelliana.

---

3. El Criterio de Lawson Ampliado y Balance de Potencia Q

En un reactor de fusión de confinamiento magnético estacionario, el estado del plasma se evalúa mediante la conservación global de potencia por unidad de volumen:

$\frac{dW}{dt} = P_{\alpha} + P_{\text{ext}} - P_{\text{rad}} - P_{\text{trans}} = 0$

donde:

• $W = 3 n_e k_B T$ es la densidad de energía térmica almacenada (asumiendo plasma quasineutro $n_e \approx n_i = n$).
• $P_{\alpha} = \frac{1}{4} n^2 \langle \sigma v \rangle E_{\alpha}$ es la potencia de autocalentamiento debida a las partículas alfa confinadas.
• $P_{\text{ext}}$ es la potencia de calentamiento auxiliar inyectada exteriormente (ondas RF y haces neutros).
• $P_{\text{rad}}$ son las pérdidas por radiación de frenado (bremsstrahlung): $P_{\text{brem}} = C_B Z_{\text{eff}} n^2 T^{1/2}$.
• $P_{\text{trans}} = \frac{W}{\tau_E} = \frac{3 n k_B T}{\tau_E}$ son las pérdidas de calor por transporte convectivo y turbulento a través del campo magnético, gobernadas por el tiempo de confinamiento de energía $\tau_E$.

3.1 Definición del Factor de Ganancia de Fusión Q

El factor de amplificación de energía del reactor se define como:

$Q = \frac{P_{\text{fusión}}}{P_{\text{ext}}} = \frac{5 P_{\alpha}}{P_{\text{ext}}}$

• Condición de Breakeven Científico ($Q = 1$): La potencia total generada iguala a la potencia auxiliar inyectada.
• Objetivo de Diseño de ITER ($Q = 10$): Con $P_{\text{ext}} = 50\text{ MW}$, generar $P_{\text{fusión}} = 500\text{ MW}$, lo que implica que el plasma es calentado en dos tercios por sus propias partículas alfa ($P_\alpha = 100\text{ MW} = 2 P_{\text{ext}}$).
• Condición de Ignición ($Q \to \infty$): El plasma se autosostiene indefinidamente sin aporte de calefacción externa ($P_{\text{ext}} = 0$).

3.2 Deducción del Producto Triple de Lawson

Para la condición de ignición ($P_\alpha \ge P_{\text{trans}} + P_{\text{brem}}$), despreciando en primera aproximación las pérdidas de radiación pura en la ventana óptima de $15\text{ keV}$:

$\frac{1}{4} n^2 \langle \sigma v \rangle E_{\alpha} \ge \frac{3 n k_B T}{\tau_E}$

Reordenando para el Producto Triple de Fusión:

$n_e \cdot T \cdot \tau_E \ge \frac{12 k_B T^2}{\langle \sigma v \rangle E_{\alpha}}$

Minimizando el miembro derecho en función de la temperatura para la mezcla D-T, se obtiene el valor crítico absoluto:

$n_e \cdot T \cdot \tau_E \ge 3.0 \times 10^{21} \text{ keV} \cdot \text{s} \cdot \text{m}^{-3} \quad \left( \approx 5.0 \times 10^{28} \text{ m}^{-3}\text{s}\cdot\text{K} \right)$

---

4. Equilibrio Magnetohidrodinámico y la Ecuación de Grad-Shafranov

Un plasma a cientos de millones de grados es un fluido conductor cuasineutro gobernado por las ecuaciones acopladas de Navier-Stokes y las ecuaciones de Maxwell, formulación conocida como Magnetohidrodinámica (MHD) ideal:

$\rho \left( \frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot \nabla)\mathbf{v} \right) = -\nabla p + \mathbf{J} \times \mathbf{B}$

$\nabla \times \mathbf{B} = \mu_0 \mathbf{J}, \quad \nabla \cdot \mathbf{B} = 0$

$\mathbf{E} + \mathbf{v} \times \mathbf{B} = 0 \quad (\text{conductividad infinita, ley de Ohm ideal})$

En estado estacionario y sin flujos macroscópicos de velocidad promedio ($\mathbf{v} = 0$), el equilibrio de fuerzas se reduce estrictamente a:

$\nabla p = \mathbf{J} \times \mathbf{B}$

Esta ecuación indica dos propiedades topológicas fundamentales:

1. $\mathbf{B} \cdot \nabla p = 0$: Las líneas de campo magnético yacen enteramente sobre superficies de presión constante (superficies magnéticas).
2. $\mathbf{J} \cdot \nabla p = 0$: Las líneas de densidad de corriente también fluyen contenidas en las superficies de presión constante.

4.1 Derivación de la Ecuación Bidimensional de Grad-Shafranov

En la geometría axisimétrica toroidal de un Tokamak, usando coordenadas cilíndricas $(R, \phi, Z)$ donde $\phi$ representa el ángulo toroidal de simetría ($\partial/\partial\phi = 0$), el campo magnético divergencia cero $\mathbf{B}$ se puede descomponer unívocamente en componentes poloidal y toroidal:

$\mathbf{B} = \nabla \psi \times \nabla \phi + F(\psi) \nabla \phi = \frac{1}{R} \left( -\frac{\partial \psi}{\partial Z} \hat{\mathbf{R}} + \frac{\partial \psi}{\partial R} \hat{\mathbf{Z}} \right) + \frac{F(\psi)}{R} \hat{\boldsymbol{\phi}}$

donde $\psi(R,Z)$ es la función de flujo poloidal magnético y $F(\psi) = R B_\phi$ es la función de corriente toroidal.

Sustituyendo esta representación en la ley de Ampère $\mathbf{J} = \frac{1}{\mu_0} \nabla \times \mathbf{B}$ e igualando con el gradiente de presión escalar $\nabla p = p'(\psi) \nabla \psi$, se obtiene la celebre Ecuación elíptica no lineal de Grad-Shafranov:

$\Delta^* \psi = -\mu_0 R^2 \frac{dp}{d\psi} - F(\psi) \frac{dF}{d\psi}$

donde el operador elíptico de Grad-Shafranov $\Delta^*$ está definido como:

$\Delta^* \psi \equiv R \frac{\partial}{\partial R} \left( \frac{1}{R} \frac{\partial \psi}{\partial R} \right) + \frac{\partial^2 \psi}{\partial Z^2} = \frac{\partial^2 \psi}{\partial R^2} - \frac{1}{R} \frac{\partial \psi}{\partial R} + \frac{\partial^2 \psi}{\partial Z^2}$

Esta ecuación determina la forma exacta de la vasija y la configuración en D elongada de la sección transversal del plasma en ITER, optimizando el factor de seguridad de línea $q_{95} \approx 3$ para evitar reconexiones magnéticas catastróficas.

---

5. El Límite de Troyon e Inestabilidades de Borde (ELMs)

La eficiencia económica de un Tokamak depende directamente de cuánto plasma caliente puede confinar con un campo magnético dado. Esto se mide mediante el parámetro adimensional Beta del plasma ($\beta$):

$\beta = \frac{\langle p \rangle}{B^2 / (2\mu_0)} = \frac{\text{Presión Cinética del Plasma}}{\text{Presión Magnética del Campo}}$

Si la presión del plasma es demasiado elevada en comparación con la tensión elástica de las líneas de campo magnético, el plasma desarrolla inestabilidades magnetohidrodinámicas destructivas (modos de globo o ballooning modes, y modos de torsión o kink modes).

5.1 El Criterio Empírico de Estabilidad de Troyon

Francis Troyon estableció analíticamente que el valor máximo estable de $\beta$ está rigurosamente acotado por la corriente de plasma y el campo magnético aplicado:

$\beta_{\text{max}} (\%) = \beta_N \left( \frac{I_p [\text{MA}]}{a [\text{m}] \cdot B_T [\text{T}]} \right)$

donde:

• $I_p$ es la corriente toroidal del plasma ($15 \text{ MA}$ en ITER).
• $a$ es el radio menor del toroide ($2.0 \text{ m}$ en ITER).
• $B_T$ es el campo toroidal en el eje magnético ($5.3 \text{ T}$ en ITER).
• $\beta_N$ es el factor de Troyon normalizado.

Para un plasma en régimen estándar de modo de alto confinamiento (Modo H), el límite operativo se sitúa en $\beta_N \approx 1.8 - 2.5$. Superar este límite desencadena una disrupción de plasma: la pérdida completa del confinamiento en menos de $15 \text{ milisegundos}$, liberando energías electromagnéticas de más de $300 \text{ MJ}$ contra las paredes del reactor.

5.2 Modos Localizados en el Borde (ELMs) y su Mitigación

En el borde del plasma se forma una barrera de transporte térmico con gradientes de presión ultra-pronunciados (el pedestal del Modo H). Este pedestal sufre periódicamente inestabilidades explosivas denominadas ELMs (Edge Localized Modes) de Tipo I, que eyectan pulsos de energía de hasta $20 \text{ MW/m}^2$ directamente hacia el desviador (divertor).

Para mitigar los ELMs, ITER implementa dos tecnologías activas:

1. Bobinas RMP (Resonant Magnetic Perturbation): 27 bobinas superconductoras no axisimétricas ubicadas dentro de la vasija de vacío que aplican minúsculas perturbaciones ergódicas en el borde magnético para desestabilizar suavemente el pedestal antes de que acumule energía destructiva.
2. Inyección de Perdigones Criogénicos (Pellet Pacing): Disparadores neumáticos que inyectan pequeños granos esféricos de hidrógeno congelado a $300 \text{ m/s}$ con frecuencias de hasta $30 \text{ Hz}$, forzando micro-ELMs controlados e inofensivos.

---

6. Superconductividad de Nb3Sn y Bobinas Criogénicas a 4.5 K

Para generar el campo toroidal estático de $11.8 \text{ Tesla}$ en el conductor ($5.3 \text{ Tesla}$ en el centro geométrico del plasma, $R_0 = 6.2 \text{ m}$), los electroimanes convencionales de cobre consumirían gigavatios de potencia eléctrica sólo por efecto Joule, haciendo inviable el balance neto de energía.

ITER utiliza 18 bobinas toroidales superconductoras gigantescas de aleación intermetálica de Niobio-Estaño ($\text{Nb}_3\text{Sn}$) y Niobio-Titanio ($\text{NbTi}$):

6.1 Diagrama Crítico del Superconductor de Niobio-Estaño

El estado superconductor sólo existe por debajo de la superficie crítica tridimensional definida por la temperatura $T$, el campo magnético $B$ y la densidad de corriente $J$:

$J_c(T, B) \approx J_{c0} \left( 1 - \frac{T}{T_c(0)} \right)^2 \left( \frac{B}{B_{c2}(T)} \right)^{-1/2} \left( 1 - \frac{B}{B_{c2}(T)} \right)^2$

• Temperatura Crítica: $T_c(0) \approx 18.3 \text{ K}$.
• Campo Crítico Superior: $B_{c2}(0) \approx 29.5 \text{ Tesla}$.
• Régimen Operativo en ITER: El sistema es enfriado con helio líquido supercrítico a $4.5 \text{ Kelvin}$ ($-268.65^\circ\text{C}$) a una presión de $5 \text{ bar}$.

La energía magnética almacenada en el campo toroidal de ITER alcanza los 41 Gigajulios (GJ), equivalente a la energía cinética de un superpetrolero a plena carga o la detonación de 10 toneladas de TNT. Las estructuras de soporte de acero austenítico inoxidable 316LN están diseñadas para resistir fuerzas de centrado hacia el eje de hasta 400 Meganewtons por bobina.

---

7. El Manto Fértil (Breeding Blanket) y el Ciclo Regenerativo del Tritio

El tritio es un isótopo inestable del hidrógeno con un período de semidesintegración radioactiva de $\tau_{1/2} = 12.32 \text{ años}$ mediante desintegración beta pura:

$^3_1\text{T} \longrightarrow \,^3_2\text{He} + e^- + \bar{\nu}_e + 18.6\text{ keV}$

Debido a su corta vida media, las reservas mundiales de tritio civil son extremadamente escasas (menos de $30 \text{ kg}$, producidas principalmente en reactores nucleares de fisión moderados por agua pesada CANDU). Un reactor comercial de $1\text{ GW}$ eléctrico consumirá aproximadamente $55.6\text{ kg}$ de tritio por año. Por tanto, la tecnología de fusión magnética sólo es viable si el reactor es capaz de autosuministrarse mediante transmutación nuclear in situ.

7.1 Reacciones de Regeneración Neutrónica en Litio

La vasija de ITER está revestida interiormente por un manto fértil (breeding blanket) enriquecido con litio metálico y cerámicas de óxido de litio ($^6\text{Li}$ y $^7\text{Li}$), aprovechando el flujo neutrónico de $14.1\text{ MeV}$:

$^6_3\text{Li} + \,^1_0\text{n}\,(\text{térmico}) \longrightarrow \,^4_2\text{He} + \,^3_1\text{T} + 4.784\text{ MeV} \quad (\sigma \approx 940 \text{ barns})$

$^7_3\text{Li} + \,^1_0\text{n}\,(\text{rápido}) \longrightarrow \,^4_2\text{He} + \,^3_1\text{T} + \,^1_0\text{n}' - 2.467\text{ MeV} \quad (\text{endotérmica})$

7.2 Razón de Reproducción de Tritio (TBR)

Para que el ciclo de combustible sea sostenible, la Razón de Reproducción de Tritio (TBR, Tritium Breeding Ratio) debe superar estrictamente la unidad:

$\text{TBR} = \frac{\text{Átomos de Tritio Producidos}}{\text{Átomos de Tritio Consumidos}} \ge 1.05 - 1.15$

Para compensar las fugas neutrónicas por aberturas de diagnóstico y bombeo de vacío, el manto incorpora multiplicadores neutrónicos basados en berilio ($^9\text{Be}$) o aleaciones eutécticas de plomo-litio ($ext{Pb-16Li}$), mediante reacciones de espalación $(n, 2n)$:

$^9_4\text{Be} + \,^1_0\text{n} \longrightarrow 2\,^4_2\text{He} + 2\,^1_0\text{n} - 1.666\text{ MeV}$

---

8. Referencias Bibliográficas y Papers Fundamentales

1. Lawson, J. D. (1957). Some Criteria for a Power Producing Thermonuclear Reactor. Proceedings of the Physical Society. Section B, 70(1), 6-10. DOI: 10.1088/0370-1301/70/1/303.
2. Grad, H., & Rubin, H. (1958). Hydromagnetic Equilibria and Force-Free Fields. Proceedings of the 2nd UN Conference on Peaceful Uses of Atomic Energy, Geneva, Vol. 31, 190-197.
3. Shafranov, V. D. (1966). Plasma Equilibrium in a Magnetic Field. Reviews of Plasma Physics, Vol. 2, Consultants Bureau, New York, 103-151.
4. Troyon, F., Gruber, R., Saurenmann, H., Semenzato, S., & Succi, S. (1984). MHD-Limits to Plasma Beta, and Favoured Profiles in Tokamaks. Plasma Physics and Controlled Fusion, 26(1A), 209-215. DOI: 10.1088/0741-3335/26/1A/319.
5. ITER Physics Expert Groups. (1999). Chapter 1: Introduction and Summary of ITER Physics Basis. Nuclear Fusion, 39(12), 2137-2174. DOI: 10.1088/0029-5515/39/12/301.
6. Wesson, J. (2011). Tokamaks (4th Edition). Oxford University Press, Oxford Master Series in Atomic, Optical and Laser Physics. ISBN: 978-0199592234.
7. Greenwald, M. (2002). Density Limits in Toroidal Plasmas. Plasma Physics and Controlled Fusion, 44(8), R27-R80. DOI: 10.1088/0741-3335/44/8/201.
8. Fasoli, A., et al. (2007). Progress in the ITER Physics Basis: Fast Particles in Tokamak Plasmas. Nuclear Fusion, 47(6), S264-S284.