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.
Design of commercial fusion energy technology must be done in a manner safe to the general public. This presentation is intended to cover the major non-radiological categories which must be considered for adequate safety of fusion devices.
Design of commercial fusion energy technology must be done in a manner safe to the general public. This presentation is intended to cover the major non-radiological categories which must be considered for adequate safety of fusion devices.
Wendelstein 7-X (W7-X) (Greifswald, Germany) is an advanced stellarator, which uses the modular coil concept to realize a magnetic configuration optimized for fusion-relevant plasma properties. The magnet system of the machine allows a variation of the rotational transform (iota) at the boundary. In the latest W7-X operational phase a dedicated configuration scan has been performed varying the rotational transform between magnetic configurations with iota = 5/4 and iota = 5/5 at the boundary. This paper presents an overview of the experiments and of the main results with respect to confinement and stability. The main observation is an increase of the plasma energy in several intermediate configurations of the scan when the 5/5-islands are close to the plasma boundary but still inside the last-closed-flux-surface. In addition, these configurations showed marked MHD-activity with a crashing behavior related to the 5/5-islands. The corresponding mode amplitude was correlated with the size of the internal 5/5 islands.
The Magnetized Liner Inertial Fusion (MagLIF) experimental platform [M. R. Gomez et al., Phys. Rev. Lett. 113, 155003 (2014)] represents the most successful demonstration of magneto-inertial fusion (MIF) techniques to date in pursuit of ignition and significant fusion yields. The pressing question remains regarding how to scale MIF concepts like MagLIF to more powerful pulsed-power drivers while avoiding significant changes in physical regimes that could adversely impact performance. In this work, we propose a conservative approach for scaling general MIF implosions, including MagLIF. Underpinning our scaling approach is a theoretical framework describing the evolution of the trajectory and thickness of a thin-walled, cylindrical, current-driven shell imploding on preheated, adiabatic fuel. By imposing that scaled implosions remain self-similar, we obtain a set of scaling rules expressing key target design parameters and performance metrics as functions of the maximum driver current I max . We identify several scaling paths offering unique, complementary benefits and trade-offs in terms of physics risks and driver requirements. Remarkably, when scaling present-day experiments to higher coupled energies, these paths are predicted to preserve or reduce the majority of known performance-degrading effects, including hydrodynamic instabilities, impurity mix, fuel energy losses, and laser-plasma interactions, with notable exceptions clearly delineated. In the absence of α heating, our scaling paths exhibit neutron yield per-unit-length scaling as $\tilde{Y}$∝[I$3\atop{max}$,I${4.14}\atop{max}$] and ignition parameter scaling as χ∝[I max ,I${2.14}\atop{max}$]. By considering the specific physics risks unique to each scaling path, we provide a roadmap for future investigations to evaluate different scaling options through detailed numerical studies and scaling-focused experiments on present-day facilities. Overall, these results highlight the potential of MIF as a key component of the national ignition effort.
A divertor for system supporting plasma for fusion employs permanent magnets providing far weaker magnetic field strength than conventionally used in cryogenic magnetic systems through an ability to place the permanent magnets in close proximity to the plasma containment volume.
Resonances of high energy particles in magnetic confinement devices due to electromagnetic instabilities can strongly modify the distribution, leading to a reduction in fusion power and even discharge termination and particle loss to the device walls through avalanche. The existence of a mode particle resonance depends on properties of the equilibrium, particle trajectories, and perturbation mode harmonic content. Resonance location is a function of particle energy and equilibrium field line helicity. Different methods for finding resonance location and energy dependence are developed. The properties of resonances are discussed using examples from magnetic fusion devices. Here, we show that if mode resonances exist at low particle energy, they very likely also exist at high energy, thus modifying high energy beam particles and fusion products. It is possible for a resonance to appear due to mode induced orbit helicity modification when it is forbidden at low mode amplitude.
Abstract The classical problem of neutral beam ions slowing down in a magnetic mirror geometry is revisited to provide predictive capability for the new Wisconsin HTS Axisymmetric Mirror under development at the University of Wisconsin. A Fokker–Planck model named fast beam ion solver (FBIS) is developed to include the spatial non-uniformity of a physical mirror geometry. The mathematical framework allows for efficient orbit averaging of the pitch-angle scattering operator, and permits a determination of the axial profile of the ambipolar potential confining the electrons. The numerical results from FBIS are consistent with earlier work, but further show how mirror-ratio and a near square magnetic well optimizes the fusion gain. The numerical results are also applied to inform the conceptual design of WHAM++, a low capital-cost breakeven-class magnetic mirror device.
Not Available
Magnetic diffusion plays an important role in inertial confinement fusion with strong magnetic fields. In this paper, we improve a previous analysis of the generation and diffusion of the magnetic field [Morita et al., Phys. Plasmas 25, 094505 (2018)]. For the generation process, we calculate the temporal evolution of the coil current using a self-consistent circuit model. The results show that the peak of the calculated magnetic field is delayed by 1.2 ns compared with that of the incident laser pulse. For the diffusion process, we evaluate the electrical conductivity of warm dense gold over a wide temperature range (300 K–100 eV) by combining the Kubo–Greenwood formula based on a quantum molecular dynamics simulation with the modified Spitzer model. Our simulation shows that the maximum magnetic field (530 T) that penetrates the cone is delayed by 2.5 ns compared with the laser peak. This result is consistent with experiments [Sakata et al., Nat. Commun. 9, 3937 (2018)] that showed that applying a strong magnetic field improved the heating efficiency of fusion fuel.
In the extreme conditions of inertial confinement fusion experiments, heat flow plays a vital role, but local diffusive models frequently break down and overestimate the heat flow. The situation becomes more complicated again in the significant magnetic fields generated during laser–plasma interactions or in magnetized fusion schemes. Accurate non-local and magnetized heat flow computations can be carried out using Vlasov–Fokker–Planck (VFP) simulations, but these are computationally expensive. There is, therefore, significant interest in using faster multi-group models to accurately calculate the non-local heat flow in magnetized plasmas. We benchmark two such multi-group models for calculating the heat flow, M1 and hybrid-AWBS-BGK, against diffusive models and full VFP simulations, before applying the models to realistic example test cases, both magnetized and unmagnetized. We find that the multi-group models generally perform very well for moderate non-localities up to kλmfp∼0.01, but the computational cost increases dramatically. hybrid-AWBS-BGK performs more effectively than M1 at high non-localities, up to kλmfp∼1, due to its adaptive solver and robust P1 closure, but tends to fail in very strong magnetic fields. Both codes are much faster than VFP simulations but are still slow in steep temperature gradients.
High thrust, high specific impulse propulsion systems, such as PuFF (Pulsed Fission-Fusion)offer a new paradigm of speed and safety when exploring the solar system. To realize thesepossibilites, PuFF needs a mature magnetic nozzle design. PuFF’s magnetic nozzle is uniquein that it must provide power in addition to thrust. In the authors’ opinion, the power systemdesign is fairly mature, but nozzle thrust coil design is less so. Previous work has suggested thatthe PuFF magnetic nozzle will behave similarly to a ‘magnetic reflection chamber,’ relying onprinciples of magnetic flux compression to re-direct nozzle flow. However, recent work, usingthe 3D modeling software SPFMax, suggests that plasma motion in the nozzle will be drivenby the Lorentz force. This recent work is limited in that the plasma did not completely leavethe nozzle before the simulation ended. In this work, we aim to extend the SPFMax modelfrom previous work, with a bolstered electromagnetic field solver, so the model can run untilthe plasma is fully evacuated. We hope to confirm previous results regarding the ideal PuFFnozzle design, but we could be surprised. Through improving the PuFF nozzle model, weaim to reduce technical risk with building the nozzle, realize paradigm-shifting performancenumbers, and allow rapid, reliable, and safe exploration of the solar system.
Engineering features are known to cause jets of ablator materials to enter the fuel hot-spot in inertial confinement fusion implosions. The Biermann battery mechanism wraps them in a self-generated magnetic field. Here we show that higher-Z jets have an additional thermoelectric magnetic source term that is not present for hydrogen jets, verified here through a kinetic simulation. It has similar magnitude to the Biermann term. We then include this in an extended magneto-hydrodynamic approach to post process an xRAGE radiation-hydrodynamic implosion simulation. The simulation includes an accurate model for the capsule fill tube, producing a dense carbon jet that becomes wrapped in a 4000 T magnetic field. A simple spherical carbon mix model shows that this insulates the electron heat conduction enough to cause contraction of the jet to an optically thick equilibrium. The denser magnetized jet hydrodynamics could change its core penetration and therefore the final mix mass, which is known to be well correlated with fusion yield degradation. Fully exploring this will require self-consistent magneto-hydrodynamic simulations. Experimental signatures of this self-magnetization may emerge in the high energy neutron spectrum.
A novel wire-array Z-pinch platform has been developed to study the effect of a pulsed magnetic mirror field on a collisional high energy density plasma. The mirror is driven in series with the Z-pinch target and, therefore, scales with driver current. Axial confinement is observed based on visible imaging and axial flow velocity measurements. The presence of axial compressing J xB force is determined indirectly based on Thomson scattering and interferometry measurements and corroborated by three-dimensional extended-MHD simulations. Compared to non-magnetized wire array Z-pinch, a modified pulsed mirror configuration is observed to increase precursor plasma thermal energy density by about 30%. If optimized, such a configuration could potentially improve magnetized liner inertial fusion performance by reducing axial plasma end loss.
A forward analytic model is required to rapidly simulate the neutron time-of-flight (nToF) signals that result from magnetized liner inertial fusion (MagLIF) experiments at Sandia’s Z Pulsed Power Facility. Various experimental parameters, such as the burn-weighted fuel-ion temperature and liner areal density, determine the shape of the nToF signal and are important for characterizing any given MagLIF experiment. Extracting these parameters from measured nToF signals requires an appropriate analytic model that includes the primary deuterium-deuterium neutron peak, once-scattered neutrons in the beryllium liner of the MagLIF target, and direct beamline attenuation. Here, mathematical expressions for this model were derived from the general-geometry time- and energy-dependent neutron transport equation with anisotropic scattering. Assumptions consistent with the time-of-flight technique were used to simplify this linear Boltzmann transport equation into a more tractable form. Models of the uncollided and once-collided neutron scalar fluxes were developed for one of the five nToF detector locations at the Z-Machine. Numerical results from these models were produced for a representative MagLIF problem and found to be in good agreement with similar neutron transport simulations. Twenty experimental MagLIF data sets were analyzed using the forward models, which were determined to only be significantly sensitive to the ion temperature. The results of this work were also found to agree with values obtained separately using a zero scatter analytic model and a high-fidelity Monte Carlo simulation. Finally, inherent difficulties in this and similar techniques are identified, and a new approach forward is suggested.
The formation of the plasma current spike at the end of the thermal quench (TQ) phase is studied systematically, which is found to strongly correlate with the onset of the 1/1 kink-tearing reconnection only in the simulation results presented here. Here, the magnetohydrodynamic activity on the q = 1 surface plays a critical role in the spike formation and the disruption process. In particular, when the safety factor in the magnetic axis q 0 exceeds 1, the plasma major disruption transits into successive minor disruptions, and the start of the TQ phase is delayed.
Experiments on the DIII-D tokamak have identi fied a novel regime in which applied resonant magnetic perturbations (RMPs) increase the particle con finement and overall performance. This work details a robust range of counter-current rotation over which RMPs cause this density pump-in effect for high con finement (H-mode) plasmas. The pump-in is shown to be caused by a reduction of the turbulent transport and to be correlated with a change in the sign of the induced neoclassical transport. Furthermore, this novel reversal of the RMP induced transport has the potential to signifi cantly improve reactor relevant, three-dimensional magnetic con finement scenarios.
The perturbed ion temperature and toroidal flow were measured in rotating Neoclassical Tearing Modes (NTM) in a tokamak for the first time. These toroidally and radially resolved profiles were obtained by impurity ion spectroscopy in a 2,1 NTM in DIII-D. In agreement with drift-kinetic simulations, the electron temperature profile is flat, while the ion temperature gradient is restored across the magnetic island O-point in the presence of fast ions; the perturbed flow has minima in the O-points and maxima at the X-points. Furthermore, these measurements provide the first confirmation of the theoretically expected ion temperature and flow response to a magnetic island needed to predict the NTM on set threshold scaling for ITER and other future tokamaks.
For the purpose of better understanding type-I edge localized mode (ELM) control in ITER with resonant magnetic perturbation (RMP) fields, the plasma response to RMP is computed by a resistive full magneto-hydrodynamic model in toroidal geometry. Five scenarios designed for ITER are considered, ranging from the pre-nuclear to nuclear phases. The plasma response to RMP is quantified by the plasma surface displacement near the X-point of the divertor plasma and at the outboard mid-plane. The optimal coil configurations between two high- Q deuterium-tritium (DT) scenarios (at the same plasma current of 15 MA and the same magnetic field of 5.3 T but different fusion gains, Q = 5 and 10) are predicted to be similar. For the other ITER scenarios with similar edge safety factor q 95 ~ 3 to that of the baseline scenario, the optimal coil phasing is also similar. The optimization results are different for a half-current full-field (7.5 MA/5.3 T) scenario, largely due to the difference in q 95 . The RMP coil currents are also optimized to tailor the core vs edge toroidal torques exerted by the 3D RMP fields on the plasma column. Torque optimization, with various objective functions proposed in the study, is useful for minimizing the side effects of RMP on the plasma core flow in ITER, while still maintaining the ELM control capability. Full utilization of three rows of ELM control coils in ITER is found to be essential to ensure both flexibility and robustness of ELM control, in terms of both linear and quasilinear plasma responses.