SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
The tearing mode stability of an inverse aspect-ratio expanded tokamak plasma equilibrium of general shape is investigated using asymptotic matching techniques. Particular emphasis is placed on the conservation of toroidal electromagnetic angular momentum. The TJ code, which is a specific implementation of the results of the investigation, is described.
This paper presents a sparsified Fourier neural operator for coupled time-dependent partial differential equations (ST-FNO) as an efficient machine learning surrogate for fluid and particle-based fusion codes such as NIMROD (Non-Ideal Magnetohydrodynamics with Rotation - Open Discussion) and GTC (Gyrokinetic Toroidal Code). ST-FNO leverages the structures in the governing equations and utilizes neural operators to represent Green's function-like numerical operators in the corresponding numerical solvers. Once trained, ST-FNO can rapidly and accurately predict dynamics in fusion devices compared with first-principle numerical algorithms. In general, ST-FNO represents an efficient and accurate machine learning surrogate for numerical simulators for multi-variable nonlinear time-dependent partial differential equations, with the proposed architectures and loss functions. The efficacy of ST-FNO has been demonstrated using quiescent H-mode simulation data from NIMROD and kink-mode simulation data from GTC. The ST-FNO H-mode results show orders of magnitude reduction in memory and central processing unit usage in comparison with the numerical solvers in NIMROD when computing fields over a selected poloidal plane. The ST-FNO kink-mode results achieve a factor of 2 reduction in the number of parameters compared to baseline FNO models without accuracy loss.
A new frequency-stepped Doppler backscattering (DBS) system has been integrated into a real-time steerable electron cyclotron heating launcher system to simultaneously probe local background turbulence (f < 10 MHz) and high-frequency (20–550 MHz) density fluctuations in the DIII-D tokamak. The launcher allows for 2D steering (horizontally and vertically) over wide angular ranges to optimize probe location and wavenumber response. The vertical steering can be optimized during a discharge in real time. The new DBS system employs a programmable frequency synthesizer with adjustable dwell time as a source to launch either O or X-mode polarized millimeter waves. This system can step in real-time over the entire E-band frequency range (60–90 GHz). This combination of capabilities allows for the diagnosis of the complex internal spatial structure of high power (>200 kW) helicon waves (476 MHz) injected from an external antenna during helicon current drive experiments in DIII-D. Broadband density fluctuations around the helicon frequency are observed during real-time scans of measurement location and wavenumber during these experiments. Analysis indicates that these broadband high-frequency fluctuations are a result of backscattering of the DBS millimeter-wave probe beam from plasma turbulence modulated by the helicon wave. It is observed that background turbulence is effectively locally “tagged” with the helicon wave electric field, forming images of the turbulent spectrum in the overall density fluctuation spectrum that appear as high-frequency sidebands of the turbulence. These observations of background turbulence and high-frequency fluctuations open up the possibility of monitoring local helicon wave amplitude by comparing the high-frequency signal amplitude to the simultaneously measured background turbulence. In combination with the real-time measurement location and wavenumber scanning capabilities (offered by real-time frequency-stepping and steering), this allows rapid determination of the spatial distribution of the helicon wave power during steady-state plasma operation.
It is possible to produce a ponderomotive effect in a plasma system without time-varying fields, if the plasma flows over spatial oscillations in the field. This can be achieved by superimposing a spatially oscillatory perturbation on a guide field, then setting up an electric field perpendicular to the guide field to drive flow over the perturbation. However, subtle distinctions in the structure of the resulting electric field can entirely change the behavior of the resulting ponderomotive force. Previous work has shown that, in slab models, these distinctions can be explained in terms of the polarization of the effective wave that appears in the co-moving frame. Here, we consider what happens to this picture in a cylindrical system, where the transformation to the co-moving (rotating) frame is not inertial. It turns out that the non-inertial nature of this frame transformation can lead to counterintuitive behavior, partly due to the appearance of parallel (magnetic-field-aligned) electric fields in the rotating frame even in cases where none existed in the laboratory frame. Apart from the academic interest of this study, the practical impact lies in being better able to anticipate the antenna configuration on the plasma periphery of a cylindrical plasma that will lead to optimal ponderomotive barrier formation in the interior plasma.
Considerable progress has been made in deuterium-tritium-layered implosion experiments on the OMEGA Laser System, bringing the prospects for thermonuclear ignition in direct-drive configurations with megajoule-class lasers closer to reality. Doing so has required navigating the balance between improved 1D performance and multidimensional stability. Using statistical modeling based on over 350 cryogenic implosions to identify various degradation mechanisms, and combined with multidimensional simulations and experimental techniques such as target offsets to combat residual flows, core conditions have repeatably been achieved that extrapolate to the burning-plasma state when scaled to 2.15 MJ of symmetric laser illumination. Using high implosion velocities (> 450 km/s) and moderately high adiabats (~5), these experiments produced record-high scaled Lawson parameters in direct drive equal to 89 ± 2% of that required for ignition with expected yields of up to 1.5 ± 0.2 MJ. To improve these results still further, focused physics studies are performed to improve physics understanding and identify routes to even greater performance. Recent studies include investigations into the impact of mounting features, laser imprint, reduced fuel temperatures, and greater on-target intensities through subscale experiments. This manuscript gives a summary of the cryogenic direct-drive program on the OMEGA laser, including routes taken to achieve the current best performance, the status of recent focused physics investigations, and future designs—such as target solutions to laser imprint and reducing vapor density to increase convergence—that are expected lead to the demonstration of hydro-equivalent ignition on OMEGA.
Understanding the physics of the onset and evolution of tearing modes (TMs) in tokamak plasmas is crucial for high-performance steady-state operations. The onset condition of (m,n) = (2,1) tearing stability, $Δ′$ > $Δ_c^′$> 0, has been studied in DIII-D steady-state hybrid discharges with accurate equilibrium reconstruction and well-measured plasma profiles. Here, m is the poloidal mode number and n is the toroidal mode number. The onset of n = 1 TMs is carefully identified using the Mirnov array to determine in-phase signals for the onset time and using electron cyclotron emission measurements to trace back from the saturated island width in the experiments for both initial island width and its growth rate. The tearing stability index $Δ′$ at TM onset is calculated in multiple ways, using the modified Rutherford equation and using the PEST-III or resistive DCON magnetohydrodynamic codes, and the results are in reasonable agreement with each other. The calculated $Δ′$ is compared with an analytical formula of the tearing stability threshold $Δ_c^′$ to determine the condition for the mode onset, $Δ′$ > $Δ_c^′$. Finally, the local gradient effect of the plasma current and the pressure profiles on the tearing stability is investigated in terms of the difference between $Δ′$ and $Δ_c^′$ through a series of modeled equilibria with the parametric variations. A database of experimental mode instability onsets is compatible with the simulated regions of instability obtained from the parametric variations.
A subgrid electron–temperature–gradient (ETG) model is demonstrated here which averages local electron-scale turbulence over intermediate scales in space and time to include in global ion–temperature–gradient (ITG) simulations. This approach results in ion-scale equations, which incorporate the electron heat transport from ETG turbulence and the effects of electron-scale turbulence on the ion scale. Flux-tube ETG Cyclone Base Case simulations are performed at different radial locations and a kinetic form of the flux is added to global ion-scale (IS) simulations as a source term. Analytic radial profiles of ETG heat flux are constructed and compared to flux-tube results at multiple radial locations. Different ratios of ITG to ETG heat flux levels are considered and the results of capturing ETG heat transport in global ion-scale simulations are discussed. The potential coupling of the ETG streamer potential and intermediate-scale zonal flows to the ion scale is further addressed.
Experimental analysis at DIII-D shows that small edge localized modes (ELMs) deposit a larger fraction of their energy to the first wall, compared to type-I ELMs in similar magnetic configuration and input power. The energy (λQ) and particle (λΓ) flux decay lengths in the scrape-off layer (SOL) are up to 3 and 5 times larger, respectively, for small ELMs than for larger type-I ELMs. Transport dynamics of ELM filaments in the SOL are found to be related to divertor conditions, where high divertor collisionality, typical for partially detached plasmas, is associated with increased cross field ELM radial fluxes. Results show that a sufficiently large outer wall gap and/or limiters might be needed in future scenarios to protect the first wall, if operating with small ELMs and a cold divertor. This might also have implications for RF heating in future devices, where the coupling efficiency is dependent to some degree to the outer-wall gap.
This work presents the hardware design and first results from a newly commissioned correlation reflectometer radiometer diagnostic that measures the cross-phase angle between electron density and temperature fluctuations in ASDEX Upgrade plasmas. This diagnostic employs cross correlations between signals measured by a tunable, continuous wave, perpendicular incidence, fluctuation reflectometer, and a 24-channel radiometer sharing the same line of sight. Novel measurements in the pedestal of a helium H-mode plasma with small edge localized modes show changes in the cross-phase angle between the electron density and temperature fluctuations from ∼90° to 120°, suggesting changes in the properties of the turbulence driving transport in the plasma edge.
Here, this article describes errata in 'Initial characterization of electron temperature and density profiles in PEGASUS spherical tokamak discharges driven solely by local helicity injection,' Phys. Plasmas 28, 102504 (2021) that correct reported electron temperature values and their implications.
Magnetically driven implosions involving a liner collapsing onto a target are inherently vulnerable to the magneto-Rayleigh–Taylor instability (MRTI). Among the various approaches proposed to achieve fusion conditions within the target, the staged Z-pinch (SZP) concept employs a high-Z liner, the advantages of which remain an active area of investigation. Consequently, ongoing design optimization efforts are essential, while critical physical processes such as magnetic field (B-field) diffusion and radiation transport increase the complexity of required simulations. In this study, we utilize the new capabilities of the FLASH code to simulate in 2D, for the first time using this code, staged Z-pinch configurations (designated SZP1 and SZP1*), focusing on their stability. First, a comparison of simulation results with theoretical predictions of MRTI growth provide new insights into the often-overlooked influences of high mesh resolution and initial perturbation seeding on instability dynamics. These findings then lead to a discussion on potential improvements for future SZP experiments. These include the use of an axial B-field stabilization and the optimization of radiation transport processes during the implosion. The results presented herein establish the framework for simulating multi-dimensional Z pinches using the FLASH code and pave the way for the development of innovative experimental configurations leveraging its advanced simulation capabilities.
Magnetized liner inertial fusion (MagLIF) experiments have demonstrated fusion-relevant ion temperatures up to 3.1 keV and thermonuclear production of up to 1.1 × 1013 deuterium–deuterium neutrons. This performance was enabled through platform development that provided increases in applied magnetic field, coupled preheat energy, and drive current. Advanced coil designs with internal reinforcement enabled an increase from 10 to 20 T. An improved laser pulse shape, beam smoothing, and thinner laser entrance foils increased preheat energy coupling from less than 1 to 2.3 kJ. A redesign of the final transmission line and load region increased peak load current from 16 to 20 MA. The wider range of input parameters was leveraged to study target performance trends with preheat energy, applied magnetic field, and peak load current. Ion temperature and neutron yield generally followed trends in two-dimensional clean Lasnex calculations. Stagnation performance improved with peak load current when other input parameters were also increased such that convergence was maintained. This dataset suggests that reducing convergence to less than 30 would improve predictability of target performance. Lasnex was used to identify a simulation-optimized scaling path, which suggests 10+ kJ of fusion yield is possible on the Z facility with achievable input parameters. This path also indicates >10 MJ could be generated through volume burn on a future facility with a path to high yield (>200 MJ) using cryogenic dense fuel layers. The newly developed MagLIF platform enables exploration of both this simulation optimized scaling path and a recently developed similarity-scaling path.
Dense plasma focii (DPFs) are appealing as energy efficient sources of short pulses of ions, neutrons, and x rays. The output of these sources is expected to scale with input current (I4), but has been shown to drop at the MA level [S. K. H. Auluck, “On the failure of neutron yield scaling in the dense plasma focus,” Phys. Plasmas 30, 080701 (2023)]. New results on the MegaJOuLe Neutron Imaging Radiography DPF showed neutron yield production in agreement with the input current scaling beyond the previously observed drop. This work provides insight into the pinch formation on a DPF and reports on the two different mechanisms leading to neutron generation inside a DPF using a combination of kinetic simulations and experimental data. A combination of particle-in-cell (PIC) and 1D shock theory results are used to describe the pinch formation and disassembly and the corresponding thermonuclear and beam-target mechanisms. The temporal evolution of the pinch column predicted by the PIC simulations shows qualitative agreement with the experimental data from plasma photon emission as well as temporal neutron pulse shapes. In MJ-class DPFs, both thermonuclear and beam-target mechanisms can occur over the course of the implosion and contribute to the total neutron production. Hence the neutron source size of a DPF will change throughout the implosion. Experimental neutron radiographs show the increase in source size as the pinch breaks apart, in agreement with simulation's prediction.
We apply the global gyrokinetic code GENE-3D to perform kinetic-ballooning-mode (KBM) turbulence simulations in the Wendelstein 7-X (W7-X) stellarator using experimentally relevant parameters. In contrast to the flux-tube approach, we observe saturation of the turbulent fluctuations at finite amplitude. Analyzing the turbulent characteristics, we identify the presence of weakly driven KBM turbulence hybridizing with ion-temperature-gradient and trapped-electron-mode turbulence in the scenario under consideration, showing that it is possible that KBM activity can be driven in W7-X, opening up the path for studying this type of turbulence, which is highly relevant for future fusion power plants.
The radiated power distribution is a crucial aspect of heat transport, heat load mitigation, and plasma exhaust performance. An imaging video bolometer camera has been installed to measure the plasma radiated power in the divertor region of the Wendelstein 7-X (W7-X) stellarator. This diagnostic offers a wide-angle (40° × 68°) sampling of the plasma volume in both the poloidal and toroidal directions. The field-of-view is covered with a large number (> 500) of bolometer channels, providing imaging capability. The diagnostic design is introduced here together with its data analysis procedure. A set of laboratory experiments is performed to assess the thermal properties of the gold absorber foil and their spatial uniformity. Following installation, a heat source originating from the inertially cooled front of the diagnostic is identified and filtered out. The discharge data indicate a satisfactory signal-to-noise ratio as well as spatiotemporal resolution. These represent the first toroidally resolved images of the line-integrated radiated power in the W7-X island divertor. The diagnostic was then upgraded with a thinner platinum absorber and adjusted mirrors. Early data from the most recent experimental campaign employing the upgraded design show a considerable improvement in the diagnostic performance with more bolometer channels (> 1400), extended coverage, and increased spatial resolution.
A kinetic ion-induced electron emission (IIEE) model for general applications is developed to obtain the emitted electron energy spectrum for a distribution of ion impacts on a metallic surface. We assume an ionization cascade mechanism and use empirical models for the ion and electron stopping powers. The emission spectrum and the secondary electron yield (SEY) are validated for a variety of materials. The IIEE model is used to study the effect of IIEE on the plasma-material interactions of Z-pinch electrodes. Un-magnetized Boltzmann-Poisson simulations are performed for a Z-pinch plasma doubly bounded by two biased copper electrodes with and without IIEE at bias potentials from 0 to 9 kV. At the anode, the SEY decreases from 0 to 1 kV, but then increases at higher bias potentials. At the cathode, the SEY is much larger due to higher energy ion bombardment and grows with bias potential. As the bias potential increases, the emitted cathode electrons are accelerated to higher energies into the domain, collisionally heating the plasma. Above 1 kV, the heating is strong enough to increase the plasma potential. Despite SEY greater than 1, only a classical sheath forms as opposed to a space-charge limited or inverse sheath due to the emitted electron flux not reaching the space charge current saturation limits. Furthermore, the current in the emissionless cases saturates to a value lower than experiment. With IIEE, the current does not saturate and continues to increase with the 4 kV case, matching most closely with the experiment.
The aim of this study is to analyze the stability of helical Alfvén eigenmodes (HAEs) in TJ-II discharges and the stabilizing effect of the energetic particles generated by the neutral beam injector (NBI) on pressure gradient-driven modes (PGDMs). HAE and PGDM stability is studied using the linear version of the gyro-fluid code FAR3d and the continuous structure by the STELLGAP code. First, Alfvén eigenmode (AE) and PGDM activity observed in the experiments is reproduced by the simulations, identifying unstable m/n = 4/7 − 2/3 and 7/12 − 5/8 HAEs triggered around ρ = 0.66 showing a frequency of 209 and 204 kHz, respectively, as well as 5/3 PGDM. Next, a parametric study is performed with respect to the thermal ion density and iota profile in the middle-outer plasma region to verify the robustness of the simulation results with respect to the uncertainty of experimental profiles. The analysis confirms that experimental uncertainty does not cause large deviations in the simulation results, showing the destabilization of the same HAEs for all the configurations tested. The simulations also indicate the decay of the 5/3 PGDM growth rate as the energetic particle (EP) population in the plasma increases, consistent with the experiment. Stability analysis of the n = 3, 7, 11, n = 5, 9, 13, n = 6, 10, 14, and n = 8, 12 helical families is performed with respect to the NBI operational regime for different EP energies, β as well as deposition profiles. The most unstable configuration is the radially localized on-axis NBI operation (stiff EP density profile gradients nearby the magnetic axis). Using the simulation model that reproduces the observed Alfvén activity, we extend the study to analyze NBI performance within a theoretical framework. It shows that increasing NBI voltage (which raises EP energy) leads to a degradation in NBI performance for a given power (related to EP β and their density). To achieve better NBI operation, higher voltage must be balanced with lower injection power, ensuring stable AEs while keeping the same EP β.
This paper presents the validation of the GENE-KNOSOS-Tango framework for recovering both the steady-state plasma profiles in the considered radial domain and selected turbulence trends in a stellarator. This framework couples the gyrokinetic turbulence code GENE, the neoclassical transport code KNOSOS, and the transport solver Tango in a multi-timescale simulation feedback loop. Ion-scale kinetic-electron and electron-scale adiabatic-ion flux-tube simulations were performed to evolve the density and temperature profiles for four OP1.2b W7-X scenarios. The simulated density and temperature profiles showed good agreement with the experimental data using a reasonable set of boundary conditions. Equally important was the reproduction of observed trends for several turbulence properties, such as density fluctuations and turbulent heat diffusivities. Key effects were also touched upon, such as electron-scale turbulence and the neoclassical radial electric field shear. The validation of the GENE-KNOSOS-Tango framework enables credible predictions of physical phenomena in stellarators and reactor performance based on a given set of edge parameters.