Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “magnetohydrodynamical simulations”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Azimuthal magnetic field distribution in gas-puff Z -pinch implosions with and without external magnetic stabilization

An experimental study of the magnetic field distribution in gas-puff Z pinches with and without a preembedded axial magnetic field (B z0 ) is presented. Spatially resolved, time-gated spectroscopic measurements were made at the Weizmann Institute of Science on a 300 kA, 1.6 μs rise time pulsed-power driver. The radial distribution of the azimuthal magnetic field, B θ , during the implosion, with and without a preembedded axial magnetic field of B z0 = 0.26 T, was measured using Zeeman polarization spectroscopy. The spectroscopic measurements of B θ were consistent with the corresponding values of B θ inferred from current measurements made with a B-dot probe. One-dimensional magnetohydrodynamic simulations, performed with the code TRAC-II, showed agreement with the experimentally measured implosion trajectory, and qualitatively reproduced the experimentally measured radial B θ profiles during the implosion when B z0 = 0.26 T was applied. Simulation results of the radial profile of B θ without a preembedded axial magnetic field did not qualitatively match experimental results due to magneto-Rayleigh-Taylor (MRT) instabilities. Furthermore, our analysis emphasizes the importance of MRT instability mitigation when studying the magnetic field and current distributions in Z pinches. Discrepancies of the simulation results with experiment are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2010

Broken Ergodicity in Ideal, Homogeneous, Incompressible Turbulence

We discuss the statistical mechanics of numerical models of ideal homogeneous, incompressible turbulence and their relevance for dissipative fluids and magnetofluids. These numerical models are based on Fourier series and the relevant statistical theory predicts that Fourier coefficients of fluid velocity and magnetic fields (if present) are zero-mean random variables. However, numerical simulations clearly show that certain coefficients have a non-zero mean value that can be very large compared to the associated standard deviation. We explain this phenomena in terms of broken ergodicity', which is defined to occur when dynamical behavior does not match ensemble predictions on very long time-scales. We review the theoretical basis of broken ergodicity, apply it to 2-D and 3-D fluid and magnetohydrodynamic simulations of homogeneous turbulence, and show new results from simulations using GPU (graphical processing unit) computers.

Morin, Lee↗
DOE OSTI2023

Three-dimensional feedback processes in current-driven metal

Using three-dimensional (3D) magnetohydrodynamic simulations, we study how a pit on a metal surface evolves when driven by intense electrical current density j. Redistribution of j around the pit initiates a feedback loop: j both reacts to and alters the electrical conductivity σ, through Joule heating and hydrodynamic expansion, so that j and σ are constantly in flux. Thus, the pit transforms into larger striation and filament structures predicted by the electrothermal instability theory. Both structures are important in applications of current-driven metal: Here, the striation constitutes a density perturbation that can seed the magneto-Rayleigh-Taylor instability, while the filament provides a more rapid path to plasma formation, through 3D j redistribution. Simulations predict distinctive self-emission patterns, thus allowing for experimental observation and comparison.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Observation of a magneto-Rayleigh-Taylor instability in magnetically collimated plasma jets

We present the direct experimental observation of the formation of a diamagnetic cavity and magneto-Rayleigh-Taylor (MRT) instability in a β ≈ 1 high energy density plasma. Proton radiography is used to measure the two dimensional path-integrated magnetic field in a laser-produced plasma propagating parallel to a preimposed magnetic field. Flutelike structures, associated with the MRT instability, are observed to grow at the surface of the cavity, with a measured wavelength of 1.2 mm and growth time of 4 ns. These measurements are in good agreement with predictions of three dimensional magnetohydrodynamic simulations using the GORGON code. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2010

Observations and Simulations of Electron Dynamics Near an Active Neutral Line

Recent observations in the Earth's magnetotail have shown rapid increases in the fluxes of energetic electrons with energies up to 100's of keV associated with dipolarization fronts that propagate into the inner magnetosphere. On August 15, 2001 the four Cluster spacecraft located slightly dawnward of midnight (yGSM approx. -5.4RE) at xGSM approx. -18RE observed a series of earthward propagating dipolarization fronts [Hwang et al., 2010]. At least 6 dipolarization fronts were observed in a 20m interval. Unlike previously reported cases the fluxes of electrons up to 95keV decreased during the passage of the first three fronts over the spacecraft. The energetic electron fluxes increased during the passage of the last three fronts. We have performed a global magnetohydrodynamic simulation of this event using solar wind observations from the ACE satellite to drive the simulation. In the simulation a very complex reconnection system in the near-Earth tail at XGSM approx. -20RE launched a series of earthward propagating dipolarization fronts that are similar to those observed on Cluster. The simulation results indicate that the Cluster spacecraft were just earthward of the reconnection site. In this paper we will present a study of the dynamics of electrons associated with these events by using the large-scale kinetic simulation approach in which we launch a large number of electrons into the electric and magnetic fields from this simulation.

Goldstein, M. L.↗
DOE OSTI2022

Magnetized ablative Rayleigh-Taylor instability in three dimensions

Three-dimensional extended-magnetohydrodynamics simulations of the magnetized ablative Rayleigh-Taylor instability are presented. Previous two-dimensional (2D) simulations claiming perturbation suppression by magnetic tension are shown to be misleading, as they do not include the most unstable dimension. For perturbation modes along the applied field direction, the magnetic field simultaneously reduces ablative stabilization and adds magnetic tension stabilization; the stabilizing term is found to dominate for applied fields > 5 T, with both effects increasing in importance at short wavelengths. Additionally, for modes perpendicular to the applied field, magnetic tension does not directly stabilize the perturbation but can result in moderately slower growth due to the perturbation appearing to be 2D (albeit in a different orientation to 2D inertial confinement fusion simulations). In cases where thermal ablative stabilization is dominant the applied field increases the peak bubble-spike height. Resistive diffusion is shown to be important for short wavelengths and long timescales, reducing the effectiveness of tension stabilization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A robust fourth-order finite-difference discretization for the strongly anisotropic transport equation in magnetized plasmas

We propose a second-order temporally implicit, fourth-order-accurate spatial discretization scheme for the strongly anisotropic heat transport equation characteristic of hot, fusion-grade plasmas. Following Du Toit et al. (2018), the scheme transforms mixed-derivative diffusion fluxes (which are responsible for the lack of a discrete maximum principle) into nonlinear advective fluxes, amenable to nonlinear-solver-friendly monotonicity-preserving limiters. The scheme enables accurate multi-dimensional heat transport simulations with up to seven orders of magnitude of heat-transport-coefficient anisotropies with low cross-field numerical error pollution and excellent algorithmic performance, with the number of linear iterations scaling very weakly with grid resolution and grid anisotropy, and scaling with the square-root of the implicit timestep. We propose a multigrid preconditioning strategy based on a lower-order approximation that renders the scheme efficient and scalable under grid refinement. Several numerical tests are presented that display the expected spatial convergence rates and strong algorithmic performance, including fully nonlinear magnetohydrodynamics simulations of kink instabilities in a Bennett pinch in 2D helical geometry and of ITER in 3D toroidal geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2015

Deep Chandra Observation and Numerical Studies of the Nearest Cluster Cold Front in the Sky

We present the results of a very deep (500 ks) Chandra observation, along with tailored numerical simulations, of the nearest, best resolved cluster cold front in the sky, which lies 90 kpc (19 arcmin) to the north-west of M87. The northern part of the front appears the sharpest, with a width smaller than 2.5 kpc (1.5 Coulomb mean free paths; at 99 per cent confidence). Everywhere along the front, the temperature discontinuity is narrower than 4-8 kpc and the metallicity gradient is narrower than 6 kpc, indicating that diffusion, conduction and mixing are suppressed across the interface. Such transport processes can be naturally suppressed by magnetic fields aligned with the cold front. Interestingly, comparison to magnetohydrodynamic simulations indicates that in order to maintain the observed sharp density and temperature discontinuities, conduction must also be suppressed along the magnetic field lines. However, the northwestern part of the cold front is observed to have a non-zero width. While other explanations are possible, the broadening is consistent with the presence of Kelvin-Helmholtz instabilities (KHI) on length-scales of a few kpc. Based on comparison with simulations, the presence of KHI would imply that the effective viscosity of the intracluster medium is suppressed by more than an order of magnitude with respect to the isotropic Spitzer-like temperature dependent viscosity. Underneath the cold front, we observe quasi-linear features that are approximately 10 per cent brighter than the surrounding gas and are separated by approximately 15 kpc from each other in projection. Comparison to tailored numerical simulations suggests that the observed phenomena may be due to the amplification of magnetic fields by gas sloshing in wide layers below the cold front, where the magnetic pressure reaches approximately 5-10 per cent of the thermal pressure, reducing the gas density between the bright features.

hydrodynamics↗
DOE OSTI2020

Taylor State Merging at SSX: Experiment and Simulation

We describe experiments and simulations of dynamical merging with two Taylor state plasmas in a Swarthmore Spheromak Experiment (SSX) device. Taylor states are formed by magnetized plasma guns at opposite ends of the device. We performed experiments with Taylor states of both senses of magnetic helicity (right-handed twist or left-handed twist). We present results of both counter-helicity merging (one side left-handed, the other right-handed) and co-helicity merging (both sides left-handed). Experiments show significant ion heating, consistent with magnetic reconnection. We suggest that the merged, warm state could be a suitable target for future magneto-inertial fusion experiments. Magnetohydrodynamic simulations of these experiments reveal the structure of the final relaxed, merged state.

Brown, Michael↗
DOE OSTI2025

Integrated simulations of premagnetized and self-magnetizing dynamic screw pinch-driven MagLIF

Magnetically driven implosions such as in magnetized liner inertial fusion (MagLIF) on the Z accelerator suffer from magneto-Rayleigh–Taylor instabilities (MRTI) that dynamically redistribute liner mass during implosion, limiting fusion fuel compression and confinement, which ultimately degrades performance. Driving the implosion with an initially helical drive field that dynamically shifts the direction of the magnetic field surrounding the liner (i.e., a dynamic screw pinch, DSP) is a method proposed to mitigate MRTI in-flight and improve target performance. In DSPs, the axial drive magnetic field component implodes the liner and diffuses through the shocked, melted liner material into the fuel throughout the implosion. Liners can be designed to enable enough axial magnetic flux to diffuse through the liner material to effectively magnetize the fuel region without the need of an initial axial magnetic field (i.e., from external field coils). We present results from three-dimensional radiation-magnetohydrodynamic simulations of MagLIF implosions employing drive magnetic fields composed of axial and azimuthal components (a helical drive field). These simulated DSP-driven MagLIF targets demonstrate improved fuel conditions and thermonuclear yield compared to a traditional MagLIF target implosion. Synthetic x-ray radiography of the imploding liner material and x-ray emission images of the fuel region at the time of peak neutron yield rate indicate superior implosion morphology for DSP-MagLIF implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Spot size measurement of a deuterium–tritium dense plasma focus using neutron radiography

Neutron radiography is a technique uniquely suited to applications in nuclear diagnostics, non-destructive testing, and subcritical experiments. The spatial resolution of neutron radiographs is degraded by optical blur in the imaging system and the neutron source size, where the ideal source is point-like to optimize the point-spread function. A potential neutron source for radiography is the dense plasma focus (DPF), a coaxial Z-pinch that produces thermonuclear and beam-target neutrons. To assess if the source size is suitable for radiography, a neutron imaging system was used to measure the source size of the 4 MA Sodium DPF at the Nevada National Security Site operating with deuterium–tritium gas-fill. The source size was measured using the edge-spread function of tungsten objects, each having a rolled (convex) edge. The spot size was found to be 7–12 mm full-width at half-max (FWHM) assuming a Gaussian source, though comparison is presented for Lorentzian and Bennett distributions. The average FWHM was found to be 8.6 ± 1.2 mm vertically and 10.8 ± 1.2 mm horizontally with respect to the image plane, averaging over varied edges and alignments. The results were sensitive to source alignment and edge metrology, which introduced notable uncertainties. These results are consistent with separate experimental measurements as well as magnetohydrodynamics simulations of this DPF, which suggest that neutron production can originate from pinches ∼5–7 mm off-axis. These results suggest that the DPF should be used for radiography at low magnification (M < 1) where spot size does not dominate spatial blur.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

X-ray imaging and radiation transport effects on cylindrical implosions

Magnetization of inertial confinement implosions is a promising means of improving their performance, owing to the potential reduction of energy losses within the target and mitigation of hydrodynamic instabilities. In particular, cylindrical implosions are useful for studying the influence of a magnetic field, thanks to their axial symmetry. In this work, we present experimental results from cylindrical implosions on the OMEGA-60 laser using a 40-beam, 14.5 kJ, 1.5 ns drive and an initial seed magnetic field of B 0 = 30 T along the axes of the targets, compared with reference results without an imposed B-field. Implosions were characterized using time-resolved x-ray imaging from two orthogonal lines of sight. We found that the data agree well with magnetohydrodynamic simulations, once radiation transport within the imploding plasma is considered. We show that for a correct interpretation of the data in these types of experiments, explicit radiation transport must be taken into account.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Transport coefficient sensitivities in a semi-analytic model for magnetized liner inertial fusion

Performance of magnetized liner inertial fusion (MagLIF) experiments is highly dependent on transport processes including magnetized heat flows and magnetic flux losses. Magnetohydrodynamic simulations used to model these experiments require a choice of model for the transport coefficients, which are the constants of proportionality relating driving terms, such as temperature gradients and currents, to the associated heat and magnetic field transport. The coefficients have been the subject of repeated recalculation using various methods throughout the years. Using a semi-analytic MagLIF model, we compare models for the transport coefficients. The choice of model modifies magnetic-flux losses caused by the Nernst thermoelectric effect and thermal conduction losses. We present simulated results from parameter scans conducted in order to compare the effects of the different models on parameters of interest in MagLIF. In some regions of parameter space, discrepancies of up to 38% are found in integrated quantities like the fusion yield. These results may serve as a guide for experimental validation of the various models, particularly as laser preheat energies and initial axial field strengths are increased on MagLIF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Two-dimensional Thomson scattering measurements of misaligned electron density and temperature gradients and associated Biermann battery produced fields

We present optical Thomson scattering measurements of electron density and temperature in high Mach number laser-driven blast waves in homogeneous gases. Taylor–Sedov blast waves are launched in nitrogen (N2) or helium (He) at pressures between 0.4 mTorr and 10 Torr by ablating a solid plastic target with a high energy laser pulse (10 J, 1012 W/cm2). Experiments are performed at high repetition rate (1 Hz), which allows one-dimensional and two-dimensional Thomson scattering measurements over an area of several cm2 by automatically translating the scattering volume between shots. Electron temperature and density in the blast wave fronts were seen to increase with increasing background gas pressure. Measured electron density and temperature gradients were used to calculate ∂B/∂t∝∇Te×∇ne. The experimentally measured ∂B/∂t showed agreement with the magnetic field probe (B-dot) measurements, revealing that magnetic fields are generated in the observed blast waves via the Biermann battery effect. The results are compared to numerical three-dimensional collisional magnetohydrodynamic simulations performed with FLASH, and are discussed in the context of spontaneous magnetic field generation via the Biermann battery effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Active Control of Alfvén Eigenmodes by Externally Applied 3D Magnetic Perturbations

The suppression and excitation of Alfvén eigenmodes have been experimentally obtained, for the first time, by means of externally applied 3D perturbative fields with different spatial spectra in a tokamak plasma. The applied perturbation causes an internal fast-ion redistribution that modifies the phase-space gradients responsible for driving the modes, determining, ultimately their existence. Hybrid kinetic-magnetohydrodynamic simulations reveal an edge resonant transport layer activated by the 3D perturbative field as the responsible mechanism for the fast-ion redistribution. Finally, the results presented here may help to control fast-ion driven Alfvénic instabilities in future burning plasmas with a significant fusion born alpha particle population.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Nonlinear magnetohydrodynamic modeling of current-drive-induced sawtooth-like crashes in the W7-X stellarator

Sawtooth-like core electron temperature crashes have been observed in W7-X experiments with electron cyclotron current drive. Here we present nonlinear single-fluid magnetohydrodynamic simulations of this phenomenon using the newly developed stellarator modeling capability of the M3D-C 1 code. The near-axis current drive gives rise to two ι = 1 resonances in the equilibrium rotational transform profile so that two consecutive (1, 1) internal kink modes are seen in the simulations. A small-amplitude crash at the inner resonance occurs first, which may correspond to the sawtooth precursors observed in the experiments. A bigger crash at the outer resonance then flattens the core temperature profile, which shows semi-quantitative agreements with experimental measurements on certain metrics such as the crash amplitude and the inversion radius of the temperature change. These results illustrate a likely mechanism of the current-drive-induced sawtooth-like crashes in W7-X and, to some extent, validate the stellarator modeling capability of M3D-C 1 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Impact of strong magnetization in cylindrical plasma implosions with applied B-field measured via x-ray emission spectroscopy

Magnetization is a key strategy for enhancing inertial fusion performance, though accurate characterization of magnetized dense plasmas is needed for a better comprehension of the underlying physics. Measured spectra from imploding Ar-doped D 2 -filled cylinders at the OMEGA laser show distinctive features with and without an imposed magnetic field. A multizone spectroscopic diagnosis leads to quantitative estimates of the plasma conditions, namely revealing a 50% core temperature rise at half mass density when a 30-T seed field is applied. Concurrently, experimental spectra align well with predictions from extended-magnetohydrodynamics simulations, providing strong evidence that the attained core conditions at peak compression are consistent with the impact of a 10-kT compressed field. These results pave the way for the validation of magnetized transport models in dense plasmas and for future magnetized laser implosion experiments at a larger scale. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Magnetized ICF implosions: Scaling of temperature and yield enhancement

This work investigates the impact of an applied magnetic field on the yield and hot-spot temperature of inertial confinement fusion implosions. A scaling of temperature amplification due to magnetization is shown to be in agreement with unperturbed two-dimensional (2D) extended-magnetohydrodynamic simulations. A perfectly spherical hot-spot with an axial magnetic field is predicted to have a maximum temperature amplification of 37%. However, elongation of the hot-spot along field lines raises this value by decreasing the hot-spot surface area along magnetic field lines. A scaling for yield amplification predicts that a magnetic field has the greatest benefit for low-temperature implosions; this is in agreement with simplified 1D simulations, but not 2D simulations where the hot-spot pressure can be significantly reduced by heat-flow anisotropy. Simulations including a P2 drive asymmetry then show that the magnetized yield is a maximum when the capsule drive corrects the hot-spot shape to be round at neutron bang time. An applied magnetic field is also found to be most beneficial for implosions that are more highly perturbed, exceeding the theoretical yield enhancement for symmetric hot-spots. Increasing the magnetic field strength past the value required to magnetize the electrons is beneficial due to the additional suppression of perturbations by magnetic tension.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗