Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “magnetic confinement”

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.

At least 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Numerical study of coronal plasma jet formation

A new scenario for solar flare eruption in the coronal holes is analyzed by using MHD stability concepts for a spheromak configuration. The stability properties of a spheromak partially embedded into a conducting surface are studied using three dimensional MHD simulations. In agreement with the analytical theory, a large degree of line-tying stabilizes the spheromak's tilt instability, while the elongation has a destabilizing effect. High-resolution nonlinear simulations also demonstrate current sheet formation at the upper surface of the spheromak, where the tilted magnetic field of the spheromak reconnects with the background magnetic field. The calculated stability threshold and the observed magnetic reconnection support a model of coronal jet eruptions where a dome-like magnetic structure grows through flux emergence on the solar surface, tilts, reconnects, and erupts. Countering the effect from elongation, line-tying strongly stabilizes a spheromak growing from a flux-emergence process, suggesting that to accelerate the onset of eruptive coronal jets, there must be magnetic reconnection at the bottom of the spheromak to detach the structure from the solar surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Further modeling of q 95 windows for the suppression of edge localized modes by resonant magnetic perturbations in the DIII-D tokamak

An improved resonant plasma response model that more accurately captures the physics of the interaction between a tokamak plasma and a resonant magnetic perturbation (RMP) is developed. The model interpolates between the linear and the nonlinear response regimes and takes into account the fact that the slip-frequency is non-zero in the nonlinear regime. The improved model is incorporated into the extended perturbed equilibrium code (EPEC) toroidal asymptotic matching code. The modified EPEC code is used to investigate RMP-induced edge-localized-mode (ELM) suppression in DIII-D H-mode discharge #145380. Somewhat surprisingly, allowing for a finite slip-frequency (i.e., relaxing the so-called no-slip constraint) is found to only slightly facilitate the locking of driven magnetic island chains to the RMP, and, hence, to only slightly facilitate RMP-induced ELM suppression. This is true despite the fact that the nature of non-locked island solutions is radically different when the no-slip constraint is imposed compared to when it is relaxed (in the first case, the widths of the island chains driven at the rational surfaces pulsate, and in the second case, they remain steady). The previously obtained conclusion that the response of a typical H-mode tokamak plasma to an RMP cannot be accurately modeled by linear theory is confirmed. Here, the previously obtained conclusion that the best agreement between theory and observations is achieved by assuming that the natural frequencies of tearing modes, in the absence of the RMP, are determined by the local equilibrium E × B velocity is also confirmed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Comparison of local and global gyrokinetic calculations of collisionless zonal flow damping in quasi-symmetric stellarators

The linear collisionless damping of zonal flows is calculated for quasi-symmetric stellarator equilibria in flux-tube, flux-surface, and full-volume geometry. Equilibria are studied from the quasi-helical symmetry configuration of the Helically Symmetric eXperiment (HSX), a broken symmetry configuration of HSX, and the quasi-axial symmetry geometry of the National Compact Stellarator eXperiment (NCSX). Zonal flow oscillations and long-time damping affect the zonal flow evolution, and the zonal flow residual goes to zero for small radial wavenumber. The oscillation frequency and damping rate depend on the bounce-averaged radial particle drift in accordance with theory. While each flux tube on a flux surface is unique, several different flux tubes in HSX or NCSX can reproduce the zonal flow damping from a flux-surface calculation given an adequate parallel extent. The flux-surface or flux-tube calculations can accurately reproduce the full-volume long-time residual for moderate kx, but the oscillation and damping time scales are longer in local representations, particularly for small kx approaching the system size.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Dissipation in the magnetic turbulence of reversed field pinch plasmas

Reversed field pinch plasmas are subject to tearing instability that creates a broad spectrum of magnetic fluctuations. The dominant fluctuations have poloidal and toroidal mode numbers (m,n) = (1,6 - 10) and can grow to 2 - 3% of the mean magnetic field. Through non-linear coupling this growth culminates in a strong reconnection event and broadening of the magnetic spectrum extending to the ion gyroradius scale. Multiple developments occur during the reconnection stage: ions and electrons are energized, magnetic fluctuation amplitudes increase, plasma flow is halted, and the toroidal magnetic flux increases in a sawtooth-like fashion as the RFP dynamo becomes stronger. Magnetic fluctuations are measured in the plasma edge at multiple radial locations from r/a = 0.75 to 0.96 to assess and characterize the magnetic turbulence. The measured spectrum perpendicular to the mean field, S(k ⊥ ), can be fit to a model spectrum consisting of power law and exponential component with one free parameter that characterizes dissipation. The measured dissipation is much larger than estimated from classical viscous or resistive dissipation, but it is consistent with a flow damping measurement of anomalous viscosity. The measurements show an evolution of the spectrum during which fluctuation power builds up in the smallest wavenumbers and cascades to the larger wavenumber due to the non-linear coupling between the linear (m,n) = (1, > 6) and the non-linear (m,n) = (0, 1) tearing mo

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Growing neoclassical tearing modes seeded via transient-induced-multimode interactions

Nonlinear extended magnetohydrodynamics simulations demonstrating seeding of neoclassical tearing modes (NTMs) via magnetohydrodynamics-transient-induced multimode interactions are presented. Simulations of NTMs are enabled by two recent NIMROD code developments: the implementation of heuristic neoclassical stresses and the application of transient magnetic perturbations (MPs) at the boundary. NTMs are driven unstable by the inherently pressure driven kinetic bootstrap current, which arises due to collisional viscosity between passing and trapped electrons. These simulations use heuristic closures that model the neoclassical electron and ion stresses. NTM growth requires a seed island, which is generated by a transiently applied MP in simulations. The capability is demonstrated using kinetic-based reconstructions with flow of a DIII-D ITER Baseline Scenario discharge (La Haye et al., in Proceedings IAEA FEC, 2020). The applied MP seeds a 2/1 NTM that grows in two phases: a slow growth phase followed by a faster robust growth phase. Additionally, an evolving sequence of higher order core modes are excited at first. Power transfer analysis shows that nonlinear interactions between the core modes and the 2/1 helps drive the initial slow growth. Once the induced 2/1 magnetic island reaches a critical width, the NTM transitions to faster robust growth, which is well described by the nonlinear modified Rutherford equation. Furthermore, this work highlights the role of nonlinear mode coupling in seeding NTMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Influence of anomalous perpendicular transport on linear tearing mode dynamics in tokamak plasmas

The analysis of a previous paper [A. Cole and R. Fitzpatrick, Phys. Plasmas 13, 032503 (2006)] that maps out all of the two-fluid response regimes of a linear tearing layer interacting with an externally generated resonant magnetic perturbation, in a large aspect-ratio tokamak plasma, is generalized to take into account realistic levels of perpendicular particle transport. A new response-regime map is obtained that differs substantially, in many respects, from the old one. The improved analysis is first employed to find all of the two-fluid growth-rate regimes of a non-interacting low-mode-number tearing mode in a large aspect-ratio tokamak plasma. The analysis is then used to determine the scaling of the error-field penetration threshold with machine parameters in large aspect-ratio tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

AC compensation of 3D magnetic diagnostic signals in DIII-D and National Spherical Torus Experiment-Upgrade (NSTX-U) for real-time application

A time domain algorithm has been developed to remove the vacuum pickup generated by both coil current (DC) and induced vessel current (AC) in real time from three dimensional (3D) magnetic diagnostic signals in the National Spherical Torus Experiment-Upgrade (NSTX-U) and DIII-D tokamaks. The possibility of detecting 3D plasma perturbations in real time is essential in modern and future tokamaks to avoid and control MHD instabilities. The presence of vacuum field pickup, due to toroidally asymmetric (3D) coils or to misalignment between sensors and axisymmetric (2D) coils, pollutes the measured plasma 3D field, making the detection of the magnetic field produced by the plasma challenging. Although the DC coupling between coils and sensors can be easily calculated and removed, the AC part is more difficult. Here, an algorithm based on a layered low-pass filter approach for the AC compensation and its application for DIII-D and NSTX-U data is presented, showing that this method reduces the vacuum pickup to the noise level. Comparison of plasma response measurements with and without vacuum compensation shows that accurate mode locking detection and plasma response identification require precise AC and DC compensations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Self-organized magnetic equilibria in tokamak plasmas with very low edge safety factor

Tokamak plasmas often exhibit self-organizing behavior in which internal modes shape the toroidal current density profile, a common example being the sawtooth instability. However, such behavior has not been studied in detail for edge safety factor below 2 due to disruptive kink instabilities that typically prevent operation in this regime. Now, steady tokamak plasmas with an edge safety factor down to 0.8 have been created in the Madison Symmetric Torus, where disruptions are prevented due to a thick, conductive wall and a feedback power supply that sustains the plasma current. Internal measurements and nonlinear magnetohydrodynamic modeling reveal a family of safety factor profiles with a central value clamped near unity as the edge safety factor decreases, indicating current profile broadening through a relaxation process. Furthermore, as the safety factor decreases, the magnetic fluctuations become irregular, and the electron energy confinement time decreases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Implementation of extreme ultraviolet spectroscopy on a sheared-flow-stabilized Z pinch

A diagnostic for extreme ultraviolet spectroscopy was fielded on the sheared-flow-stabilized (SFS) fusion Z-pinch experiment (FuZE-Q) for the first time. The spectrometer collected time-gated plasma emission spectra in the 5–40 nm wavelength (30–250 eV) range for impurity identification, radiative power studies, and for plasma temperature and density measurements. The unique implementation of the diagnostic included fast (10 ns risetime) pulsed high voltage electronics and a multi-stage differential pumping system that allowed the vacuum-coupled spectrometer to collect three independently timed spectra per FuZE-Q shot while also protecting sensitive internal components. Analysis of line emission identifies oxygen (N-, C-, B-, Be-, Li-, and He-like O), peaking in intensity shortly after maximum current (>500 kA). This work provides a foundation for future high energy spectroscopy experiments on SFS Z-pinch devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Calculation of tearing mode stability in an inverse aspect-ratio expanded tokamak plasma equilibrium

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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Resistive hose modes in tokamak runaway electron beams II

Resistive hose instabilities of runaway electron (RE) beams immersed in resistive background plasmas are examined with analytic and numerical calculations. The RE beam-plasma equilibria considered are characteristic of the situation observed post-thermal quench in a tokamak disruption. An analytic linear dispersion relation is presented for the case of a uniform RE current density profile with a sharp boundary in cylindrical geometry. Initial value linear calculations for a more general profile in toroidal geometry find that reducing aspect ratio increases the resistive hose mode growth rates with fixed safety factor profile. Nonlinear calculations in cylindrical and toroidal geometry find that the resistive hose instability-driven fluctuations relax the gradient of the current density profile. In toroidal geometry, changes to the magnetic topology are observed as a result of the resistive hose activity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simulations of the churning mode: Toroidally symmetric plasma convection and turbulence around the X-points in a snowflake divertor

Using a reduced MHD model, extended to include field-aligned thermal conduction, we present numerical simulations of the churning mode (CM): a toroidally symmetric, non-linear plasma vortex in the vicinity of the null points in a snowflake (SF) divertor [D. D. Ryutov et al., Phys. Scr. 89, 088002 (2014)]. Simulations are carried out across a range of inter-null separations, d xx ⁠, and inter-null orientations, θ⁠, primarily in conditions relevant to the MAST-U tokamak. We find that, when d xx is small, the CM induces additional transport across the X-points when β pm ≳8%, where β pm is the ratio of the plasma pressure in the null region to poloidal magnetic pressure at the midplane. This transport also increases approximately linearly as d xx is reduced. A diffusive model of this transport is shown to predict the total transport across the null points, where diffusion coefficients of up to ∼ 10 2 m 2 s −1 centered on a small region around the X-points are used. However, the CM also results in significant changes to the flux surfaces in the null region, which is not captured by this diffusive model. The changes in magnetic geometry mean the fractional exhaust power delivered to each divertor leg is highly sensitive to β pm ⁠, d xx ⁠, and θ⁠. For small values of θ⁠, the CM can induce a change in topology, redirecting exhaust power from a secondary divertor leg on the high field side to one on the low field side. Similar behavior is found in the fraction of exhaust power going to the inner and outer divertor. Such changes in the flux surfaces may not be captured by Grad–Shafranov solvers and so may be a source of error in the magnetic reconstruction of SF experiments. We consistently find that the fractional exhaust power going to a secondary divertor leg on the high field side is small, consistent with SF experiments.

Computer simulation↗
DOE OSTI2025

Simulations of plasma current induced by toroidal field ramping down on tokamaks

In tokamak experiments, the BT ramping down can drive off-axis and parallel inductive current. This approach leads to a decrease in q95 and an increase in normalized beta, βN. The off-axis inductive current also broadens the current profile. Typically, the 1.5D transport code is used to simulate the time evolution of the plasma current profile, which is based on the flux-surface-averaged Faraday’s law. The ONETWO code is one of those transport codes. However, this code cannot simulate the situation of the evolving toroidal field BT. This article proposes a new time-dependent model to take into account the BT ramping down situation. A modified formula of flux-surface-averaged Faraday’s law was derived to consider the effect of the BT ramping rate on the current evolution, and it was implemented in the ONETWO code. Then, the modified ONETWO code was used to simulate the BT ramping down experiment on DIII-D. The simulation result of the plasma current evolution with BT ramping down shows a broader current profile with smaller ohmic current induced by the poloidal field, compared to that without BT ramping down.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Quantum information approach to the implementation of a neutron cavity

Abstract Using the quantum information model of dynamical diffraction we consider a neutron cavity composed of two perfect crystal silicon blades capable of containing the neutron wavefunction. We show that the internal confinement of the neutrons through Bragg diffraction can be modelled by a quantum random walk. Furthermore, we introduce a toolbox for modelling crystal imperfections such as surface roughness and defects. Good agreement is found between the simulation and the experimental implementation, where leakage beams are present, modelling of which is impractical with the conventional theory of dynamical diffraction. Analysis of the standing neutron waves is presented in regards to the crystal geometry and parameters; and the conditions required for well-defined bounces are derived. The presented results enable new approaches to studying the setups utilizing neutron confinement, such as the experiments to measure neutron magnetic and electric dipole moments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

The Multiphysics on Advanced Platforms Project

In 2015, the Lawrence Livermore National Laboratory started development of next-generation multiphysics simulation capabilities for the National Nuclear Security Administration under the Advanced Technologies Development and Mitigation (ATDM) element of the Advanced Simulation and Computing program in collaboration with the Exascale Computing Project (ECP). A key driver for this effort across the NNSA tri-lab was the emergence of advanced high performance computing (HPC) architectures based on heterogeneous compute capabilities, including GPU based systems, as part of the national drive toward exascale computing platforms at multiple Department of Energy (DOE) facilities. Developing a multiphysics code capable of meeting the various simulation needs of the NNSA as defined by the current generation of integrated codes (or ICs), initially developed as part of the Accelerated Strategic Computing Initiative (ASCI) program beginning in 1996, and able to scale to the current 100 petaflop class pre-exascale systems, as well the forthcoming exaflop class computers, is a daunting challenge. To accomplish this ambitious goal, LLNL has embraced two key themes: use of high-order numerical methods and a modular approach to code development. The LLNL next generation effort is organized under the Multi-Physics on Advanced Platforms Project (MAPP). A foundational component of MAPP is the Axom computer science (CS) toolkit which provides infrastructure for the development of modular, performance portable, multi-physics application codes. MARBL is a next-generation application code built on the Axom base to address the modeling needs of the high energy density physics (HEDP) community for simulating high-explosive, magnetic or laser driven experiments such as inertial confinement fusion (ICF), pulsed-power magneto-hydrodynamics (MHD), equation of state (EOS) and material strength studies as part of the NNSA’s stockpile stewardship program (SSP).

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1976

The galactic distribution (in radius and Z) of interstellar molecular hydrogen

Observations of the galactic longitude and latitude distributions of gamma = 2.6 mm CO emission are presented. Analysis of this spectral line data yields the large scale distribution of molecular clouds in the galactic disk and their z-distribution out of the disk. Strong maxima in the number of molecular clouds occur in the galactic nucleus and at galactic radii 4-8 kpc. The peak at 4-8 kpc correlates well with a region of enhanced 100 Mev gamma ray emissivity. This correlation strongly supports the conclusion that the gamma rays are produced as a result of cosmic ray interactions in molecular H2 clouds rather than HI. One important implication of this is that the interstellar magnetic field lines to which cosmic rays are confined must therefore not be excluded from these dense clouds. The width of the cloud layer perpendicular to galactic plane between half density points is 105 + or - 15 pc near the 5.5 kpc peak. The total mass of molecular gas in the interior of the galaxy exceeds that of atomic hydrogen.

Scoville, N. Z.↗
NASA NTRS1983

A static model of chromospheric heating in solar flares

The response of the solar chromosphere to flare processes, namely nonthermal electrons, thermal conduction, and coronal pressure, is modeled. Finite difference methods employing linearization and iteration are used in obtaining simultaneous solutions to the equations of steady-state energy balance, hydrostatic equilibrium, radiative transfer, and atomic statistical equilibrium. The atmospheric response is assumed to be confined to one dimension by a strong vertical magnetic field. A solution is obtained to the radiative transfer equation for the most important optically thick transitions of hydrogen, magnesium, and calcium. The theoretical atmospheres discussed here are seen as elucidating the role of various physical processes in establishing the structure of flare chromospheres. At low coronal pressures, conduction is found to be more important than nonthermal electrons in establishing the position of the transition region. Only thermal conduction can adequately account for the chromospheric evaporation in compact flares. Of the mechanisms considered, only nonthermal electrons bring about significant heating below the flare transition region.

Ricchiazzi, P. J.↗
NASA NTRS1984

Study of alternate optical and fine guidance sensor designs for the space infrared telescope facility (SIRTF)

A unique optical design was developed that compensates for the coma degraded images caused by field chopping in SIRTF. The conic constants of a Cassegrain telescope were altered to compensate for the coma induced by the secondary mirror tilt. The modulation transfer function is essentially independent of secondary mirror tilt, and diffraction limited image quality is maintained over a several arcminute field during chopping. With an untilted secondary mirror, the coma compensated (CC) design has a smaller field than the unchopped Ritchey-Chretien design; but use of relay optics, such as the inverted Cassegrain design developed for the fine guidance sensor (FGS), can increase the CC telescope's field size. A reactionless secondary mirror chopper mechanism that uses superconducting magnets was studied. The heart producing elements are confined to a reaction plate that is not directly viewed by the IR focal plane. A design was also developed for a low moment of inertia, reticulated HIP beryllium secondary mirror consistent with blank fabrication technology and optical finishing requirements.

Wissinger, A.↗