Efecto pellizco en el plasma: cómo el tokamak atrae partículas hacia el centro por sí solo
Contexto
Este trabajo se sitúa en la periferia de las investigaciones sobre turbulencia y transporte en el plasma de borde de los tokamaks. La mayoría de los estudios vecinos examinan el transporte turbulento de calor y partículas hacia afuera, fijando la densidad en la frontera interna o utilizando modelos simplificados. Aquí, por primera vez en una simulación fluida global con un perfil de densidad autoconsistente, se muestra que el pellizco surge de manera natural, sin suposiciones artificiales. Esto conecta el trabajo con cuestiones de larga data en la física del confinamiento, como la formación del pedestal de densidad en el modo H.
Métodos
Se utilizó el código GDB, un modelo electromagnético tridimensional global basado en las ecuaciones de Braginskii reducidas por deriva. Se simuló un plasma con parámetros cercanos a los de Alcator C-Mod: radio mayor de 68,5 cm, radio menor de 22 cm, campo de 5,3 T. La región de cálculo abarcaba toda la capa anular toroidal desde 5,3 cm dentro de la última superficie cerrada hasta 3,3 cm en la capa de desprendimiento (scrape-off layer). El perfil de densidad inicial era plano, la fuente de partículas estaba localizada cerca de la separatriz y las fuentes de calor en la frontera interna. La simulación duró 8,4 ms de tiempo de transporte, lo que permitió que el perfil de densidad alcanzara un estado cuasiestacionario. Para analizar los mecanismos, se realizaron cálculos adicionales con la conductividad térmica electrónica y los flujos de calor transversales desactivados.
Resultados
La simulación mostró que, a partir del perfil inicial plano, en unos pocos milisegundos se forma un perfil con un pico de densidad en el centro de la región de cálculo. El flujo de partículas hacia el interior existe durante toda la evolución, pero su naturaleza cambia. En la etapa temprana, cuando el gradiente de densidad es positivo (la densidad aumenta hacia afuera), domina el flujo turbulento asociado con la termodifusión electrónica. El análisis lineal de un modelo simplificado mostró que, cuando la relación de gradientes de densidad y temperatura es negativa (η<0), las ondas de deriva crean un desfase entre las fluctuaciones de densidad y potencial que produce un flujo negativo (hacia adentro). En la etapa tardía, cuando el perfil de densidad se vuelve plano, la contribución turbulenta se debilita, pero surge un flujo E×B de equilibrio debido a la asimetría poloidal de la densidad y el potencial. Esta asimetría es generada por el flujo de calor de Pfirsch-Schlüter, pero la teoría neoclásica clásica predice un flujo total nulo debido a la respuesta de Boltzmann de los electrones. En el modelo GDB se incluyen efectos no ideales (resistividad, inercia de electrones, inducción electromagnética) que crean un pequeño desfase poloidal δs ≈ −0.08π entre la densidad y el potencial. Este desfase conduce a un flujo finito hacia adentro que explica cuantitativamente la acumulación de densidad. Los cálculos de control con los flujos de calor transversales desactivados mostraron que el pico de densidad se disipa, confirmando el papel de la asimetría poloidal.
Implicaciones
Este trabajo se diferencia de los estudios vecinos en que no impone el perfil de densidad desde el exterior, sino que permite que evolucione de manera autoconsistente. Es la primera demostración, en el marco de un modelo fluido global del borde, de que el pellizco de partículas surge naturalmente y conduce a un pico de densidad. Se muestra que la imagen neoclásica clásica es incompleta: tener en cuenta las no idealidades en la ley de Ohm abre un canal de transporte adicional. Esto cambia la comprensión del balance de partículas en el plasma de borde y proporciona una base para predicciones más precisas de los perfiles de densidad.
Desarrollo futuro
El tema podría desarrollarse hacia la inclusión de una geometría más realista (por ejemplo, elongación, divertor) y efectos cinéticos para plasmas calientes. Es importante verificar si el pellizco E×B no ideal descubierto se mantiene en modelos girocinéticos y con menor colisionalidad. También es prometedor investigar el papel de la ionización asimétrica poloidal de neutros, que en tokamaks reales podría intensificar o debilitar el pellizco.
Impacto
Los resultados influirán en el modelado del confinamiento del plasma en tokamaks y stellarators, especialmente en la predicción de perfiles de densidad y la transición a modos de confinamiento mejorado. También son útiles para interpretar experimentos de inyección de gas y de pellets.
Próximos pasos
Los próximos pasos incluyen la comparación con perfiles de densidad experimentales en Alcator C-Mod y otras máquinas, así como la extensión del modelo para incluir correcciones cinéticas y geometría real.
Problemas abiertos clave
Este trabajo está directamente relacionado con el problema no resuelto de la formación del pedestal de densidad en el modo H y con la cuestión de qué mecanismos determinan el pico de densidad en diferentes regímenes de confinamiento. También arroja luz sobre el antiguo enigma de por qué en el modo I no hay barrera para partículas, aunque sí la hay para el calor.
En la simulación, las partículas 'nadan contra la corriente': se mueven desde el borde hacia el centro, aunque el gradiente de densidad inicialmente apunta hacia afuera. Este transporte contraintuitivo es posible gracias a la correlación entre las fluctuaciones de densidad y el campo eléctrico, no a la difusión ordinaria.
La idea de que el plasma se comprima espontáneamente hacia el centro recuerda al concepto de 'pellizco' en la ciencia ficción, como en la novela 'Mundo Anillo' de Larry Niven, donde los campos magnéticos contienen el plasma en un reactor.
Cifras clave
- duración de la simulación: 8.4 ms
- desfase δs: ≈ −0.08π
- pico de densidad: 1.6×1020 m^-3
- velocidad de crecimiento de la densidad: en 4 ms el perfil se vuelve plano, en 8 ms alcanza el pico
- relación de gradientes η: de −2 a +2