Summary of the FESAC Transformative Enabling Capabilities Panel Report
Explore the source record for details and available documents.
Nuclear science and engineering, Physics - Plasma ↗
Engineering topics
Chacon, L.
Publications and source records attributed to Chacon, L..
Wavelet methods for studying the onset of strong plasma turbulence
Recent simulations have demonstrated that coherent current sheets dominate the kinetic-scale energy dissipation in strong turbulence of magnetized plasma. Wavelet basis functions are a natural tool for analyzing turbulent flows containing localized coherent structures of different spatial scales. Here, wavelets are used to study the onset and subsequent transition to fully developed turbulence from a laminar state. Originally applied to neutral fluid turbulence, an iterative wavelet technique decomposes the field into coherent and incoherent contributions. In contrast to Fourier power spectra, finite time Lyapunov exponents, and simple measures of intermittency such as non-Gaussian statistics of field increments, the wavelet technique is found to provide a quantitative measure for the onset of turbulence and to track the transition to fully developed turbulence. The wavelet method makes no assumptions about the structure of the coherent current sheets or the underlying plasma model. Temporal evolution of the coherent and incoherent wavelet fluctuations is found to be highly correlated (a Pearson correlation coefficient of >0.9) with the magnetic field energy and plasma thermal energy, respectively. The onset of turbulence is identified with the rapid growth of a background of incoherent fluctuations spreading across a range of scales and a corresponding drop in the coherent components. This is suggestive of the interpretation of the coherent and incoherent wavelet fluctuations as measures of coherent structures (e.g., current sheets) and dissipation, respectively. The ratio of the incoherent to coherent fluctuations Ric is found to be fairly uniform in the turbulent state across different plasma models and provides an empirical threshold of ~0.1 for turbulence onset. The utility of this technique is illustrated through examples. First, it is applied to the Kelvin–Helmholtz instability from different simulation models including fully kinetic, hybrid (kinetic ion/fluid electron), and Hall MHD simulations. Second, the wavelet diagnostic is applied to the development of turbulence downstream of the bowshock in a global magnetosphere simulation. Finally, the wavelet technique is also shown to be useful as a de-noising method for particle simulations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Kinetic studies of ICF implosions
Here, kinetic effects on inertial confinement fusion have been investigated. In particular, inter-ion-species diffusion and suprathermal ion distribution have been analyzed. The former drives separation of the fuel constituents in the hot reacting core and governs mix at the shell/fuel interface. The latter underlie measurements obtained with nuclear diagnostics, including the fusion yield and inferred ion burn temperatures. Basic mechanisms behind and practical consequences from these effects are discussed.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A curvilinear, fully implicit, conservative electromagnetic PIC algorithm in multiple dimensions
Here, we extend a recently proposed fully implicit PIC algorithm for the Vlasov–Darwin model in multiple dimensions (Chen and Chacón (2015) [1]) to curvilinear geometry. As in the Cartesian case, the approach is based on a potential formulation (Φ, A), and overcomes many difficulties of traditional semi-implicit Darwin PIC algorithms. Conservation theorems for local charge and global energy are derived in curvilinear representation, and then enforced discretely by a careful choice of the discretization of field and particle equations. Additionally, the algorithm conserves canonical-momentum in any ignorable direction, and preserves the Coulomb gauge ∇ • A = 0 exactly. An asymptotically well-posed fluid preconditioner allows efficient use of large cell sizes, which are determined by accuracy considerations, not stability, and can be orders of magnitude larger than required in a standard explicit electromagnetic PIC simulation. We demonstrate the accuracy and efficiency properties of the algorithm with numerical experiments in mapped meshes in 1D-3V and 2D-3V.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Scalable implicit incompressible resistive MHD with stabilized FE and fully-coupled Newton–Krylov-AMG
Here, we discuss that the computational solution of the governing balance equations for mass, momentum, heat transfer and magnetic induction for resistive magnetohydrodynamics (MHD) systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena, as well as the significant range of time- and length-scales that the interactions of these physical mechanisms produce. This paper explores the development of a scalable, fully-implicit stabilized unstructured finite element (FE) capability for 3D incompressible resistive MHD. The discussion considers the development of a stabilized FE formulation in context of the variational multiscale (VMS) method, and describes the scalable implicit time integration and direct-to-steady-state solution capability. The nonlinear solver strategy employs Newton–Krylov methods, which are preconditioned using fully-coupled algebraic multilevel preconditioners. These preconditioners are shown to enable a robust, scalable and efficient solution approach for the large-scale sparse linear systems generated by the Newton linearization. Verification results demonstrate the expected order-of-accuracy for the stabilized FE discretization. The approach is tested on a variety of prototype problems, that include MHD duct flows, an unstable hydromagnetic Kelvin–Helmholtz shear layer, and a 3D island coalescence problem used to model magnetic reconnection. Initial results that explore the scaling of the solution methods are also presented on up to 128K processors for problems with up to 1.8B unknowns on a CrayXK7.
97 MATHEMATICS AND COMPUTING↗