SEARCH · Search NASA

Results for “MHD modeling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Fusion gain from cylindrical liner-driven implosions of field reversed configurations

MagLIF experiments [M.R. Gomez et al., Phys. Plasmas 22, 056306 (2015)] on Z have demonstrated the basic principles of Magneto-Inertial Fusion (MIF) for wall confined plasmas. Other MIF schemes have been proposed based on the liner implosion of closed field magnetically confined plasmas such as Field Reversed Configurations (FRCs) [T. P. Intrator et al., Phys. Plasmas 15, 042505 (2008)]. We present a semi-analytical model of liner driven FRC implosions that predicts the fusion gain of such systems. The model predicts a fusion gain near unity for an FRC imploded by a liner driven by the Z Machine. We show that FRCs could be formed and imploded at the Z facility using the AutoMag liner concept [S. A. Slutz et al., Phys. Plasmas 24, 012704 (2017)]. Furthermore, an initial bias magnetic field can be supplied by the external magnets used in MagLIF experiments. The reverse field is then supplied by an AutoMag liner, which has helical conducting paths imbedded in an insulating substance. Experiments [Shipley et al., Phys. Plasmas 26, 052705 (2019)] have demonstrated that AutoMag can generate magnetic fields greater than 30 Tesla inside of the liner. We have performed 2D Radiation MHD simulations of the formation and implosion of an FRC, which are in good agreement with the analytical model. The FRC formation process could be studied on small pulsed power machines delivering about 1 MA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Energy Cost Analysis and Operational Range Prediction Based on Medium- and Heavy-Duty Electric Vehicle Real-World Deployments across the United States

While the market for medium- and heavy-duty battery-electric vehicles (MHD EVs) is still nascent, a growing number of these vehicles are being deployed across the U.S. This study used over 2.3 million miles of operational data from multiple types of MHD EVs across various regions and operating conditions to address knowledge gaps in total cost of ownership and operational range. First, real-world energy cost savings were determined: MHD fleets should experience energy cost savings each year from 2021 to 2035, regardless of vehicle platform, with the greatest savings seen in transit buses (up to USD 4459 annually) and HD trucks (up to USD 3284 annually). Second, to help fleets across various geographies throughout the U.S. assess the suitability of EVs for their year-round operating needs, operational range was modeled using the XGBoost algorithm (R2: 70%) given 22 input features relevant to vehicle efficiency. Finally, this paper recommends (1) that MHD fleets apply energy-saving practices to minimize the impacts of cold temperatures and high congestion levels on vehicle efficiency and range, and (2) that local hauling fleets select trucks with a nominal range nearly double the expected maximum daily range to account for range losses under local, urban driving conditions.

Qiu, Yin (ORCID:0009000900948794)↗

DOE OSTI2025

3D modeling of n = 1 RMP driven heat fluxes on the SPARC tokamak PFCs using HEAT

3D heat flux calculations at the lower outer divertor plate of SPARC using the HEAT code show that 3D fields generated from error field correction coils can lead to enhanced peak heat fluxes up to 15 times larger compared to the axisymmetric case. Previously employed to simulate axisymmetric heat flux on 3D plasma facing components, the HEAT code can now predict 3D heat flux generated by non-axisymmetric plasmas. This is achieved via a new HEAT module which leverages the 3D field line tracing capabilities of MAFOT starting from an M3D-C1 (MHD resistive code) perturbed equilibria. The resulting heat flux is assigned using the magnetic footprint and the heat flux layer model, an extension of the 2D heat flux model also known as the Eich, to 3D non-axisymmetric plasmas. For SPARC, the new capabilities of HEAT are used to calculate the 3D heat loads resulting from n = 1 perturbation fields (with n indicating the toroidal periodicity) applied through a toroidal array of six picture frame coils with different amplitude. The comparison with the unperturbed case shows significant changes in shape and intensity of the heat flux profile. The results show that the application of n = 1 3D field leads to a localized enhancement of the heat flux peak, influenced by the wetted area impacted by the magnetic footprint, and the appearance of a secondary heat flux peak, whose intensity depends on amplitude of the applied 3D field and toroidal location.

3D heat flux↗

DOE OSTI2024

CRADA Final Report: CRADA Number NFE-22-09330 with General Fusion

General Fusion is developing a magnetized target fusion (MTF) approach that involves compressing an initial magnetically confined plasma inside a cavity formed in liquid metal. This approach builds from concepts initially developed under the Linus program at the U.S. Naval Research Laboratory and combines it with advances from compact toroid experiment (CTX) and sustained spheromak plasma experiment (SSPX) in compact toroid plasmas and coaxial Marshall gun systems. Modeling the tokamak during compression is central to designing a successful MTF device. The plasma is formed by coaxial helicity injection in the General Fusion device. Immediately after formation, the plasma has a diverted tokamak configuration with a single null. As the wall moves inwards, the plasma is repelled from the conducting surface and driven inwards by currents induced by its magnetic field in the liquid metal wall. As the liquid metal closes (or bridges) the opening of the coaxial plasma injector, the magnetic field topology alters to remove the null. Due to this, the plasma moves from a diverted to a wall-limited configuration. The liquid metal liner continues to close in and change shape, reducing in radius by a factor of ten at the peak of plasma compression. A model of the MTF plasma must be able to handle this continually varying geometry, and to be predictive, it must faithfully include the real imperfections arising in the process. In this project, we pursued a Monte Carlo approach to closures for MHD by computing kinetic electron trajectories in an MHD plasma background from simulations of GF devices. This requires enhancing the capabilities of the KORC-T code for running large ensembles of kinetic trajectories by porting it to GPU architectures and enabling workflows for large ensembles on OLCF machines. With these capabilities, it is possible to produce a large library of kinetic calculations of electron orbits evolving in plasma configurations spanning the magnetic configurations and plasma density profiles, including non-axisymmetry, arising in the General Fusion’s existing PI3 spherical tokamak device. Using ensembles will capture particles passing a single point in space in a given magnetic configuration, and the entire dataset will cover a range of global magnetic field geometries. By sampling around many starting points, this dataset will capture the spatial dependence of the plasma parameters. From this large dataset, it is possible to produce a reduced model for the kinetic effects not captured in MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Optical design and synthetic analysis of the electron cyclotron emission imaging diagnostic of HL-2M tokamak

The electron cyclotron emission imaging (ECEI) diagnostic is a powerful tool to study the MHD and turbulent transport in magnetically confined fusion plasmas. In this work, the optical system including the local oscillator (LO) coupling and radio frequency (RF) receiving optics, has been designed and analyzed for the HL-2M ECEI diagnostic. The LO optics can illuminate the antenna array and drive the mixer diode to work efficiently, with more than 36% of the beam intensity for the channels which probe the plasma edge relative to those of channels which probe the plasma core. The RF optics aims at guiding the plasma emission signal to the antenna array. To meet different physical requirements, three types of field of view have been achieved for the plasma imaging, with zoom factors of around 1, 1.5 and 2, respectively. The focal surfaces are almost flat, with a maximum off-mapping (defined as the radial distance between the beam waists of the lower-/uppermost antennas and electron cyclotron emission layer) less than 2.5 cm, which can match the electron cyclotron emission layer quite well and suppress the image distortion in the plasma edge region. Furthermore, the impact of possible vibrations and installation error on the focal plane has been evaluated. In addition, the predefined MHD and turbulence perturbations are well reproduced by combining the optical simulation results and synthetic ECEI modeling, which further verifies the good performance of the RF optics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Hints of primordial magnetic fields at recombination and implications for the Hubble tension

Primordial Magnetic Fields (PMFs), long studied as relics of the early Universe, accelerate recombination and have been proposed as a way to relieve the Hubble tension. However, previous studies relied on simplified toy models. Here, in this study, we use recent evaluations of recombination with PMFs, incorporating full magnetohydrodynamic (MHD) simulations and detailed Lyman-alpha radiative transfer, to test PMF-enhanced recombination (bΛCDM) against observational data from the cosmic microwave background (CMB), baryon acoustic oscillations (BAO), and Type Ia supernovae (SN). Focusing on non-helical PMFs with a Batchelor spectrum, we find a preference for presentday total field strengths of approximately 5-10 pico-Gauss. Depending on the dataset combination, this preference ranges from mild (∼ 1.8σ with Planck + DESI) to moderate (∼ 3σ with Planck + DESI + SH0ES-calibrated SN) significance. The bΛCDM has Planck + DESI χ 2 values equal to or better than ΛCDM while predicting a higher Hubble constant. Future high-resolution CMB temperature and polarization measurements will be crucial for confirming or further constraining PMFs at recombination. Field strengths of 5-10 pico-Gauss align closely with those required for cluster magnetic fields to originate entirely from primordial sources, without the need for additional dynamo amplification.

Jedamzik, Karsten [Centre National de la Recherche↗

DOE OSTI2022

2D Magnetohydrodynamic Simulations of the Electrothermal Instability in Metallic Liners

The Virginia Tech (VT) Plasma Dynamics Laboratory Computational (PDCL) and Lawrence Livermore National Laboratory (LLNL) are performing two dimensional (2D) simulations of the electrothermal instability (ETI) using the LLNL multi-physics code Ares. Ares is a multi-physics arbitrary–Lagrangian-Eulerian (ALE) code developed by LLNL and is of particular use in studying magnetohydrodynamic (MHD) instabilities like the ETI due to its resistive MHD, magnetic diffusion, and radiative-hydrodynamics packages. Among its capabilities, it has the ability to model material strength, perform adaptive mesh refinement (AMR), and incorporate a wide variety of equations of state models and conductivity models. The 2D Ares simulation model created by VT-LLNL for studying the development and growth of the electrothermal instability has been configured with initial conditions based on the Mykonos Electrothermal Instability II (METI-II) experiments described by this grant and conducted by team members at the University of Nevada (UNR), the University of New Mexico (UNM), and Sandia National Laboratories. Previously, preliminary 2D Ares simulations of the ETI had been run to approximately 80ns. The rods in these simulations were initiated with sinusoidal perturbations at a similar order of magnitude to those measured on the aluminum rods used for the Mykonos experiment. This model has been improved by increasing the spatial resolution of the simulations and running the simulations further in time. In addition to the simulation run-times extending, the preliminary sinusoidal perturbation has been replaced with a perturbation derived from amplitude measurements by the experimental team, thereby correlating the simulation inputs better to the experimental runs. These new runs are capable of reaching 120ns of simulated time for the uncoated cases and to 200ns the 41 μm coated cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Examining transport and integrated modeling predictive capabilities for negative-triangularity scenarios

This paper investigates the predictive capabilities of TGYRO and TGLF models in assessing the performance of negative triangularity (NT) plasmas compared to positive triangularity (PT) plasmas in fusion devices. TGYRO predicts kinetic profiles, while TGLF analyzes turbulent transport. The study reveals that TGYRO reasonably predicts NT profiles similar to PT, although it overpredicts the high-power scenarios where there is increased experimental MHD activity. TGLF analysis finds reduced linear growth rates in NT and altered flux spectra relative to PT. Additionally, the TGLF SAT0 saturation model is observed to predict high-k transport and a reduction of particle transport with the electron temperature gradient. These findings are further corroborated by core-pedestal modeling using the Stability Transport Equilibrium Pedestal workflow, showing stronger confinement improvements in NT, particularly at higher power densities for the SAT0 saturation model. Furthermore, the study underscores the importance of accurately capturing turbulence saturation mechanisms for NT in order to project its performance accurately in fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Photoionization of seeded combustion products as a method of enhancing the efficiency of magnetohydrodynamic power generators

Here, in this study, we performed an experimental and computational investigation into the feasibility of utilizing photoionization to enhance the electrical conductivity of seeded oxy-fuel combustion products and improve the performance of magnetohydrodynamic (MHD) power generators. We applied a variety of optical and microwave diagnostics to study the ionization and recombination processes of potassium excited by an excimer laser in a high-velocity oxy-fuel free jet. Computational fluid dynamic (CFD) simulations were performed to model the thermophysical properties and species densities of the free jet. The CFD results were validated with position-dependent potassium concentration measurements. Electron recombination exponential lifetimes were measured through time-resolved microwave transmission. The experimental electron lifetimes were compared with lifetimes calculated from CFD-predicted species densities and literature recombination rates. It was determined that K + or O 2 are the most likely recombination partners for photoionized electrons. Time-resolved fluorescence measurements provided evidence of an ionization pathway involving a two-photon ionization of KOH . Finally, a zero-dimensional chemical kinetic model was developed to assess the fundamental viability of inducing a non-equilibrium electron population to provide a net energy return in combustion-driven MHD power generators. We determined that a high energy return is feasible for targeting electrode boundary layers with ultraviolet photoionization. We also found that photoionization could potentially lower the required temperature of the bulk gas flow.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2023

Development of MHD for NekRS

A spectral-element-based formulation of incompressible MHD is presented in the context of the open-source fluid-thermal code, Nek5000/RS. The formulation supports magnetic fields in a solid domain that surrounds the fluid domain. Several steady-state and time-transient model problems are presented as part of the code verification process. Nek5000/RS is designed for large-scale turbulence simulations, which will be the next step with this new MHD capability.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

UNR ETI Computational Investigation (Year 2)

The LANL ASC Arbitrary Lagrange/Eulerian (ALE) code, FLAG, is being used to model experiments conducted under this grant in order to understand the physics limiting high-current conduction and thermal plasma production. This includes the effects of nonlinear magnetic diffusion (NLMD), the magnetic Rayleigh-Taylor (MRT) instability, and the electro-thermal instability (ETI). FLAG is a multi-dimensional, unstructured-grid, multi-material, multi-temperature, shockphysics code with a single-fluid, resistive magnetohydrodynamics (MHD) package. Equations-ofstate (EOS), conductivity, material-strength, and melt can be treated with a variety of analytic or tabular models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Modelling of sawtooth-induced fast ion transport in positive and negative triangularity in TCV

Internal kinks are a common magneto hydro-dynamic (MHD) instability observed in tokamak operation when the q profile in the plasma core is close to unity. This MHD instability impacts both the transport of the bulk plasma (current, particle and energy transport) and minority species, such as fast ions. In TCV (R 0 /a = 0.88 m/0.25 m) the fast ion population is generated in the plasma by neutral beam tangential injection of energies up to 28 keV. TCV features 16 active shaping coils permitting a great flexibility in plasma shape, including negative triangularity (δ) configurations that show surprisingly high confinement. This study focuses on the transport of fast ions induced by sawteeth, by comparing two triangularity cases and simulation results with experimental data. Comparison of two equilibria with opposite δ shows that the fast ion drifts are larger for δ < 0. Furthermore, the sawtooth-induced transport in this case is larger than δ > 0 in similar conditions. Additionally, comparison with experimental data confirms the dominance of the modification of thermal kinetic profiles following the sawtooth crash in explaining drops in the neutron rates and FIDA signals. Additional fast ion diffusion, however, improves the interpretation of the experimental data. For δ < 0, the amplitude of the perturbation better representing the experimental data is larger. Finaly, an exploratory study for 50 keV particles (soon available in TCV) shows that the situation does not worsen for such particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Energy Conversion and Electron Acceleration and Transport in 3D Simulations of Solar Flares

Recent observations and simulations indicate that solar flares undergo extremely complex 3D evolution, making 3D particle transport models essential for understanding electron acceleration and interpreting flare emissions. In this study, we investigate this problem by solving Parker’s transport equation with 3D MHD simulations of solar flares. By examining energy conversion in the 3D system, we evaluate the roles of different acceleration mechanisms, including reconnection current sheet (CS), termination shock (TS), and supra-arcade downflows (SADs). We find that large-amplitude turbulent fluctuations are generated and sustained in the 3D system. The model results demonstrate that a significant number of electrons are accelerated to hundreds of keV and even a few MeV, forming power-law energy spectra. These energetic particles are widely distributed, with concentrations at the TS and in the flare looptop region, consistent with results derived from recent hard X-ray (HXR) and microwave (MW) observations. By selectively turning particle acceleration on or off in specific regions, we find that the CS and SADs effectively accelerate electrons to several hundred keV, while the TS enables further acceleration to MeV. However, no single mechanism can independently account for the significant number of energetic electrons observed. Instead, the mechanisms work synergistically to produce a large population of accelerated electrons. Our model provides spatially and temporally resolved electron distributions in the whole flare region and at the flare footpoints, enabling synthetic HXR and MW emission modeling for comparison with observations. These results offer important insights into electron acceleration and transport in 3D solar flare regions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Non-propagating structures and propagating waves in solar wind turbulence revealed by simulations and observations

Structures and waves are common features of solar wind turbulence at various scales. The interplay between structures and waves is important for processes such as the turbulent energy cascade, plasma heating, and particle scattering. Our understanding of turbulence has been advanced by not only new space missions and numerical simulations, but also techniques that have been developed to interpret the rapidly growing turbulence data. We review basic models of turbulence with a specific focus on the analysis methods for understanding magnetic structures and waves. MHD and kinetic waves in single-spacecraft time series measurements can be identified through mode decomposition or their characteristic polarization signatures. The structures in this paper are considered as zero-frequency, non-propagating or convected modes embedded in the solar wind. The synergy between observations and simulations is most evident in the application of spatial-temporal analysis to multi-spacecraft observation and turbulence simulations. The spatial-temporal analysis has greatly improved our understanding of structures and waves in turbulence. We conclude by discussing prospects for future research.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Beam modulation and bump-on-tail effects on Alfvén eigenmode stability in DIII-D

Beam modulation effects on Alfvén eigenmode stability have been investigated in a recent DIII-D experiment and show that variations in neutral beam modulation period can have an impact on the beam driven Alfvén eigenmode spectrum and resultant fast ion transport despite similar time-averaged input power. The experiment was carried out during the current ramp phase of L-mode discharges heated with sub-Alfvénic 50–80 kV deuterium neutral beams that drive a variety of Alfvén eigenmodes unstable. The modulation period of two interleaved beams with different tangency radii was varied from shot to shot in order to modify the relative time dependent mix of the beam pitch angle distribution as well as the persistence of a bump-on-tail feature near the injection energy (a feature confirmed by imaging neutral particle analyzer measurements). As the beam modulation period is varied from 7 ms to 30 ms on/off (typical full energy slowing down time of τ slow ≈ 50 ms at mid-radius), toroidicity-induced Alfvén eigenmodes (TAEs) located in the outer periphery of the plasma become intermittent and coincident with the more tangential beam. Core mode activity changes from reversed shear Alfvén eigenmodes (RSAEs) to a mix of RSAE and beta-induced Alfvén eigenmodes. Discharges with 30 ms on/off period do not have a persistent bump-on-tail feature, have the lowest average mode amplitude and least fast ion transport. Furthermore, detailed analysis of an individual TAE using TRANSP kick modeling (Monte Carlo evolution of the distribution function with probabilistic 'kicks' by the AEs) and the resistive MHD code with kinetic fast ions, MEGA, find no strong role of energy gradient drive due to bump-on-tail features. Instead, the observed TAE modulation with interleaved beams is likely a pitch angle dependent result combined with slowing down of the tangential beam between pulses. For the conditions investigated, bump-on-tail contributions to TAE drive were found to be 5% or less of the total drive at any given time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Conference Report on the 7th International Symposium on Liquid metals Applications for fusion (ISLA-7)

Supported by the world magnetic fusion research community, a series of International Symposia on Liquid metals Applications for fusion (ISLA) have been held biannually since 2010. The 7th edition (ISLA-7) was held for the period from 12 December through 16 December 2022, at Chubu University located in Kasugai, Aichi, Japan. For the first time in the history of this series of symposia, ISLA-7 was held in a hybrid fashion, due to the COVID-19 situation. The total number of the participants was 60, 34 out of whom attended the symposium in person, and the rest participated online. As to the presentation statistics, 29 papers were presented in person, whereas 21 presentations were delivered online but real-time by the presenters in China, Spain, the UK, and the USA. Both of the presentations delivered in person and online were recorded, and the video has been shared by all participants. These participants represent 11 countries: China, Czech, Italy, Japan, Latvia, Netherlands, Russia, Thailand, the UK, and the USA. All these numbers are among the largest in this series of symposia. Covered by these presentations are; in session-2, program overviews and liquid metal research review; in session-3, liquid metal flows, and MHD issues; in session-4, liquid metal facilities; in sessions-5 and 6, liquid metal experiments and modeling; in session-7, divertor physics and heat flux mitigation; in session-8, plasma and liquid metals interactions; in session-9 liquid metal plasma-facing components, erosion, and wettability. In addition, there were an opening session whereby several opening addresses were delivered and also a closing session where all technical session summaries were presented by the respective session chairs.

divertor↗

DOE OSTI2024

Technical Impacts of Light-Duty and Heavy-Duty Transportation Electrification on a Coordinated Transmission and Distribution System

In this study, we propose a strategy to model the required spatiotemporal charging demand from light-duty (LD) and medium- and heavy-duty (MHD) electric vehicles (EVs) using actual transportation data by mapping the demand for the required EV charging to a realistic and coordinated distribution and transmission electric grid at the predicted times of the day to study their impact on the power system in a variety of load, weather, and EV penetration scenarios. This work is the first study that includes the actual weather data and transportation data with realistic and coordinated distribution and transmission grid data in a large industry-scale level study. The main goal of this study is to identify possible issues and required upgrades in the electric grid, caused by an increase in EV integration. The transmission case study is a large grid with 6717 buses over a Texas footprint, and the distribution grid is over Houston, a city in Texas, covering over three million customers. The resulting overloads and voltage violations experienced in the system are discussed, and required planning upgrades to avoid these issues are suggested.

AC optimal power flow (AC-OPF)↗

DOE OSTI2022

Slow Shock Formation Upstream of Reconnecting Current Sheets

The formation, development, and impact of slow shocks in the upstream regions of reconnecting current layers are explored. Slow shocks have been documented in the upstream regions of magnetohydrodynamic (MHD) simulations of magnetic reconnection as well as in similar simulations with the kglobal kinetic macroscale simulation model. They are therefore a candidate mechanism for preheating the plasma that is injected into the current layers that facilitate magnetic energy release in solar flares. Of particular interest is their potential role in producing the hot thermal component of electrons in flares. During multi-island reconnection, the formation and merging of flux ropes in the reconnecting current layer drives plasma flows and pressure disturbances in the upstream region. These pressure disturbances steepen into slow shocks that propagate along the reconnecting component of the magnetic field and satisfy the expected Rankine–Hugoniot jump conditions. Plasma heating arises from both compression across the shock and the parallel electric field that develops to maintain charge neutrality in a kinetic system. Shocks are weaker at lower plasma β, where shock steepening is slow. While these upstream slow shocks are intrinsic to the dynamics of multi-island reconnection, their contribution to electron heating remains relatively minor compared with that from Fermi reflection and the parallel electric fields that bound the reconnection outflow.

79 ASTRONOMY AND ASTROPHYSICS↗