Characterization of α-particle transport in reactor relevant burning plasmas

Abstract

Nuclear fusion has the potential to provide humanity with a safe, clean, abundant, efficient and reliable energy source for the generations to come, but up to date finding a viable fusion reactor concept remains an ongoing area of research. One of the main difficulties to attain economically viable magnetically controlled thermonuclear fusion reactors is the confinement of α-particles. These α-particles are responsible of sustaining the extreme temperatures required for nuclear reactions, and their loss poses a serious threat to the reactor operational control and to its plasma-facing components. In toroidally shaped fusion devices with a non-uniform magnetic field, α-particles with small parallel velocity become trapped between areas of the high field bouncing between reection points, whose position is highly susceptible to field corrugations. With the exception of symmetric magnetic fields, like those of ideal tokamaks, these socalled trapped α-particles experience non-zero radial average drifts, which might lead to their collisionless losses. There are two principal collisionless mechanisms connecting trapped particle losses with the inhomogeneities of the confining magnetic field. The first is ripple trapping, in which particles fall into local ripples and experience strong radial drifts usually being convective (ballistic). The second mechanism is ripple induced stochastic processes with milder drifts caused by the radial motion of particle reection points, which result either in the banana tip stochastic diffusion or particle transitions, in which the particles change the orbit type near the reection points. While the mitigation of these losses is widely considered in the literature on fusion reactor designs, far too little attention was paid to the statistical characterization of the processes underlying collisionless transport of trapped α-particles, whose nature is generally considered diffusive. This thesis is intended to provide such statistical description and clarify the nature of collisionless trapped α-particle transport for reactor scale configurations in cases of broken symmetry of the underlying magnetic field. To this end, detailed analyses were performed on large ensembles of α-particle trajectories, calculated with the guiding center orbit following Monte Carlo code MOCA for several magnetic configurations: a purely toroidal model with ITER parameters and four quasi-toroidal stellarators loosely based on NCSX project with different levels of magnetic field symmetry. The simulations suggest that while the perfect toroidal magnetic field symmetry of the ITER configuration grants perfect confinement, an increasing departure from quasi-toroidal symmetry leads to faster and larger α-particle losses, most of which belong to particles born with a small parallel velocity in areas of a weak magnetic field on the outer midplane of the configurations. Based on the resulting numerical trajectories, novel techniques were developed capable to calculate the fraction of trapped α-particles and identify the orbit types. Estimates show that about a third of the particles are trapped for ITER, and a fifth for the stellarators, independently on the level of symmetry. These trapped particles mainly follow banana orbits and, to a lesser extent, potato, transition and ripple trapped orbits. Statistical analysis was done for the basic parameters of banana orbits, and it was found that the most probable banana width becomes wider, and that the most probable bouncing time becomes longer as configuration departs from toroidal symmetry. The results of the trapped particle fractions and the most probable bouncing times are in agreement with those obtained by an independent numerical procedure based on the depth of the confining magnetic field and the assumption that α-particles move along the filed lines. To that end, a new figure of merit measuring the level of toroidal symmetry was introduced. The convection velocity and the diffusion transport coefficients were estimated by two methods: using the most probable banana widths and bouncing times, and fitting the time dependence of the moments of the radial probability density functions of banana centers, which were calculated with a new algorithm based on the positions of the reection points. Their statistical analysis suggests that the collisionless transport of trapped α-particles cannot be properly described as being diffusive when the magnetic configuration departs from symmetry. The assumption that ripple-enhanced radial transport of trapped α-particles is diffusive has been extensively used to model experimental data. But it is limited to describe only Gaussian and Markovian transport processes and thus neglects correlations, memory, and spatial effects, that have recently been proved relevant for fusion plasma, especially in cases of turbulent driven transport. In this thesis, the techniques used in characterizing the non-diffusive dynamics of turbulent transport were adapted to study collisionless α-particle neoclassical transport. To build an effective transport model, α-particle trajectories were analyzed with a whole set of tools imported from fractional transport theory. Using the Eulerian propagator and Lagrangian rescale range [R=S] analysis techniques, the Hurst H, the spatial α and the temporal β exponents appearing in fractional transport theory were estimated to describe non-diffusive transport. The results for the ideal toroidally symmetric ITER ripple-less magnetic configuration analyzed by the Lagrangian [R=S] method show an almost zero Hurst exponent pointing out, as expected, to the absence of radial transport. While all perfectly confined trapped α-particles were analyzed for ITER, for the four stellarators, only the particles contributing the most to the losses were considered, i.e. α-particles that get lost in the region with the steepest slope in the loss fraction. The estimated spatial and temporal transport exponents found indicate that the underlying nature of transport is nondiffusive with non-Gaussian and non-Markovian statistics. As the level of toroidal symmetry decreases, the presence of spatial correlations, particularly strong anti-correlations, becomes more pronounced. For all stellarators, there are signs of self-similarity and significant memory effects. The agreement in the Hurst exponents, estimated by both the Lagrangian and Eulerian techniques, shows that as the level of quasitoroidal symmetry increases transport becomes strongly subdiffusive. Although, the validity of the fractional model itself becomes doubtful in the limiting high and low symmetry cases. The work presented in the thesis can be naturally extended to study the validity of the fractional transport model onto other types of confining magnetic fields and various α-particle-related effects, such as collisions, α-particle birth profiles, etc.

La fusión nuclear tiene el potencial de abastecer a la humanidad con una energía segura, limpia, abundante, eficiente y fiable para las generaciones venideras, pero hasta la fecha, encontrar un concepto viable de reactor de fusión es un área de investigación en curso. Una de las mayores dificultades a la hora de conseguir reactores termonucleares económicamente viables es el confinamiento de las partículas α. Estas partículas son las encargadas de mantener las temperaturas extremas que se requieren para las reacciones nucleares y su pérdida supone una seria amenaza para la operación y control del reactor y el de sus componentes en contacto con el plasma. En dispositivos de fusión con forma toroidal y con un campo magnético no uniforme, las partículas α con una velocidad paralela pequeña, quedan atrapadas entre las zonas del campo alto, rebotando entre puntos, cuya posición es muy sensible a las ondulaciones del campo. Exceptuando los campos magnéticos simétricos, como los de los tokamaks ideales, dichas partículas α atrapadas experimentan derivas radiales promedio distintas de cero, que pueden conducir a pérdidas no colisionales. Hay dos mecanismos principales no colisionales que conectan las pérdidas de partículas atrapadas con las heterogeneidades del campo magnético confinante. El primero es el debido al ripple trapping, en el que las partículas caen en ripples locales y experimentan fuertes desviaciones radiales que generalmente son convectivas (balísticas). El segundo mecanismo consiste en procesos estocásticos inducidos por los ripples donde las leves desviaciones causadas por el movimiento radial de los puntos de rebote de las partículas dan como resultado una difusión estocástica al cambiar la trayectoria las partículas cerca de los puntos de reeflxón. Mientras que la mitigación de estas pérdidas ha sido ampliamente estudiada en la literatura relativa a los diseños de reactores de fusión, no se ha prestado mucha atención a la caracterización estadística de los procesos de transporte de las partículas α atrapadas no colisionales, cuya naturaleza en general se ha considerado difusiva. El objetivo de esta tesis es aportar dicha descripción estadística y aclarar la naturaleza del transporte de las partículas α atrapadas no colisionales en configuraciones de tipo reactor cuando hay una ruptura de simetría α en su campo magnético. Con este propósito, realizaremos análisis detallados en muestras amplias de trayectorias de partículas α calculadas con el código Monte Carlo MOCA de seguimiento de órbitas del centro guía, para varias configuraciones magnéticas: un modelo puramente toroidal con parámetros de ITER y cuatro stellarators cuasi-toroidales, inspirados en el proyecto NCSX, con diferentes niveles de simetría magnética. Las simulaciones sugieren que, si bien la perfecta simetría toroidal del campo magnético de la configuración ITER garantiza un confinamiento perfecto, una desviación creciente de la simetría cuasitoroidal conduce a mayores y más rápidas pérdidas de partículas α, la mayoría de las cuales pertenecen a partículas nacidas con un ángulo de ataque pequeño en regiones de campo magnético bajo en la zona externa del plano ecuatorial de las configuraciones. Sobre la base de las trayectorias numéricas resultantes, se desarrollaron nuevas técnicas capaces de calcular la fracción de partículas α atrapadas e identificar los tipos de órbitas. Las estimaciones muestran que alrededor de un tercio de las partículas están atrapadas para ITER, y un quinto para los stellarator, independientemente del nivel de simetría. Estas partículas atrapadas siguen principalmente órbitas de tipo banana y, en menor medida, órbitas potatoes, en tránsito y ripple trapped. En el análisis estadístico de los parámetros básicos de las órbitas banana se encontró que el ancho más probable de las bananas se hace mayor, y que el tiempo de rebote más probable se hace más largo cuando la configuración magnética carece de simetría toroidal. Los resultados de las fracciones de partículas atrapadas y los tiempos de rebote más probables están de acuerdo con los obtenidos por un procedimiento numérico independiente basado en la profundidad del campo magnético confinante y la suposición de que las partículas α se mueven a lo largo de las líneas del campo. Con este propósito, se introdujo una nueva figura de mérito que mide el nivel de simetría toroidal. Las velocidades convectivas y los coeficientes de difusión se estimaron con dos métodos: utilizando los valores más probables de los anchos de las banana y los tiempos de rebote y ajustando la dependencia temporal de la función de distribución radial de probabilidad de los centros de las órbitas de las bananas, que se calcularon con un nuevo algoritmo basado en las posiciones de los puntos de reflexión. Su análisis estadístico sugiere que el transporte de las partículas α atrapadas no colisionales no puede ser descrito como difusivo cuando la configuración magnética pierde la simetría. La suposición de que el transporte radial de partículas α atrapadas inducido por ripple es difusivo se ha utilizado ampliamente para modelar datos experimentales. No obstante, su aplicación se limita a describir solo los procesos de transporte Gaussianos y Markovianos y por lo tanto no tiene en cuenta las correlaciones, la memoria y los efectos espaciales, que recientemente han demostrado ser relevantes en plasmas de fusión, especialmente en casos de transporte turbulento. En esta tesis, las técnicas utilizadas para caracterizar la dinámica no difusiva del transporte turbulento se adaptaron para estudiar el transporte neoclásico de partículas α no colisionales. Para construir un modelo de transporte efectivo, se han analizado las trayectorias de las partículas α con un conjunto de herramientas importadas de la teoría de transporte fraccionario. Usando las técnicas de propagadores Eulerianos y el análisis [R=S] Lagrangiano, se estimaron el exponente de Hurst H y los exponentes espacial α y temporal β que aparecen en la teoría del transporte fraccionario para describir el transporte no difusivo. Los resultados para la configuración magnética sin ripple de ITER con simetría toroidal ideal analizada por el método [R=S] Lagrangiano muestran un exponente de Hurst casi igual a cero que señala, como se esperaba, la ausencia de transporte radial. Mientras que para ITER se analizaron las trayectorias de todas las partículas α atrapadas perfectamente confinadas, para los cuatro stellarators, solo se consideraron las partículas que más contribuyeron a las pérdidas, es decir, partículas α que se pierden en la región con la pendiente más pronunciada en la fracción de pérdidas. Los exponentes espaciales y temporales estimados que se encontraron indican que la naturaleza subyacente del transporte es no difusiva y con estadísticas no Gaussianas y no Markovianas. A medida que disminuye el nivel de simetría toroidal, la presencia de correlaciones espaciales, particularmente fuertes anti-correlaciones, se vuelve más pronunciada. Para todos los stellarators, hay signos de autosimilaridad y efectos significativos en la memoria. El acuerdo entre los exponentes de Hurst, estimados por las técnicas Lagrangiana y Euleriana, muestra que a medida que el nivel de simetría aumenta, el transporte se vuelve fuertemente subdifusivo. Aunque, la validez del modelo fraccionario en sí mismo se vuelve dudosa en los casos limites de alta y baja simetría. El trabajo presentado en esta tesis puede extenderse naturalmente para estudiar la validez del modelo de transporte fraccionario en otros tipos de campos magnéticos confinantes y estudiar varios efectos relacionados con las partículas α, como colisiones, perfiles de nacimiento de las partículas α, etc.