SEARCH · Search NASA

Results for “Magnetohydrodynamics”

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 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Simulation Models for Exploring Magnetic Reconnection

Simulations have played a critical role in the advancement of our knowledge of magnetic reconnection. However, due to the inherently multiscale nature of reconnection, it is impossible to simulate all physics at all scales. For this reason, a wide range of simulation methods have been crafted to study particular aspects and consequences of magnetic reconnection. This article reviews many of these methods, laying out critical assumptions, numerical techniques, and giving examples of scientific results. Plasma models described include magnetohydrodynamics (MHD), Hall MHD, Hybrid, kinetic particle-in-cell (PIC), kinetic Vlasov, Fluid models with embedded PIC, Fluid models with direct feedback from energetic populations, and the Rice Convection Model (RCM).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

A Large Eddy Simulation Study of Flow Turbulence, Alumina Transport, and Bath Temperature Evolution in Conventional Aluminum-Smelting Cell Using OpenFOAM

In this study, a Large Eddy Simulation (LES) of the aluminum-smelting process is performed using OpenFOAM. To understand the coupled behavior of heat transfer, mass transfer, and flow of the smelting process, a multi-physics computational fluid dynamics (CFD) model based on the Eulerian–Eulerian multi-fluid approach is adopted. The model accounts for CO 2 bubble and magnetohydrodynamics (MHD)-driven flow, along with alumina dissolution, transport, and bath temperature evolution. The simulation predictions show small-scale turbulent vortical structures in the anode–cathode space caused by combined effect of MHD and CO 2 bubble-bath interactions and relatively large-scale asymmetric vortices in the inter-anode space caused by the CO 2 bubble-bath interactions. The vortex formation at the edges of the anodes evidently aids in transporting alumina from the central channel to the bottom of the anodes and prevents accumulation of gas bubbles in the periphery of the anode bottom. Symmetric bath cold spots are observed in the vicinity of the feeder. Cold spots are also observed in the anode–cathode distance space below the anode bottom due to the transport of undissolved solid to this region by the flow. The findings from the work are useful in developing and designing alumina-feeding strategy leading to reduced anode effects and smooth operation of the cell. The work also highlights the important flow structures in conventional aluminum-smelting cell.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Three-step LES-C models for flows at high Reynolds numbers

Here, we investigate the need for the second correction step in the recently proposed LES-C models for fluid flows at high Reynolds numbers. These models use a predictor-corrector idea to enhance the efficiency of the existing Large Eddy Simulation models. Different three-step(one defect step, two corrections) LES-C models, based on the Leray-α, ADM and NS-ω LES models, are tested in three different situations. The new Leray-α-C2 model (C2 stands for two correction steps) is applied to the Navier–Stokes equations; the ADC2 is applied to the MagnetoHydroDynamic flow; and the NS-ω-C2 is used in the fluid-fluid interaction problem. We evaluate the effectiveness of the second correction step in all these settings, using qualitative and quantitative numerical tests.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Neutrino transport in general relativistic neutron star merger simulations

Abstract Numerical simulations of neutron star–neutron star and neutron star–black hole binaries play an important role in our ability to model gravitational-wave and electromagnetic signals powered by these systems. These simulations have to take into account a wide range of physical processes including general relativity, magnetohydrodynamics, and neutrino radiation transport. The latter is particularly important in order to understand the properties of the matter ejected by many mergers, the optical/infrared signals powered by nuclear reactions in the ejecta, and the contribution of that ejecta to astrophysical nucleosynthesis. However, accurate evolutions of the neutrino transport equations that include all relevant physical processes remain beyond our current reach. In this review, I will discuss the current state of neutrino modeling in general relativistic simulations of neutron star mergers and of their post-merger remnants. I will focus on the three main types of algorithms used in simulations so far: leakage, moments, and Monte-Carlo scheme. I will review the advantages and limitations of each scheme, as well as the various neutrino–matter interactions that should be included in simulations. We will see that the quality of the treatment of neutrinos in merger simulations has greatly increased over the last decade, but also that many potentially important interactions remain difficult to take into account in simulations (pair annihilation, oscillations, inelastic scattering).

Foucart, Francois (ORCID:0000000346174738)↗

DOE OSTI2021

In-situ magnetic alignment model for additive manufacturing of anisotropic bonded magnets

In this work, we report a mathematical framework which predicts the degree of alignment (DoA) in an in-situ aligned additively manufactured 3D printed bonded magnets. A multiphysics model is developed which couples the harmonious interactions of magnetic particles in a viscous flowing polymer under the presence of an external magnetic field. The hydrodynamic fluid-particle interaction is paired with the magnetophoretic force to predict the particle trajectory and distribution during extrusion through a 3D printer nozzle. Succeeding the force balance, a magnetohydrodynamic torque equilibrium analysis is performed to predict the net-orientation of the magnetic particles as a function of the applied field strength, viscous forces, and particle-to-particle interactions (P2P). Experimental validation of the DoA predictions is performed using 65 vol% Nd-Fe-B+Sm-Fe-N in Nylon12 (DoA exp = 0.620 and DoA theory = 0.686) and 15 vol% Sm-Co in PLA (DoA exp = 0.830 and DoA theory = 0.863). A parametric analysis is performed to analyze the effect of operating and design parameters like alignment field strength, magnetic loading fraction, extrusion load, and particle size. The model predicts a competing behavior between particle-fluid and particle-particle interactions under the presence of an applied field. Lastly, the model provides a framework to efficiently predict the DoA in tandem with a functionalized-magnetic 3D printer and allows the user to adjust the operating parameters according to the desired DoA.

36 MATERIALS SCIENCE↗

DOE OSTI2025

A review of high order strong stability preserving two-derivative explicit, implicit, and IMEX methods

High order strong stability preserving time discretizations ensure the nonlinear non-inner-product strong stability properties of spatial discretizations suited for the stable simulation of hyperbolic PDEs in a wide variety of application areas including fluid dynamics, magnetohydrodynamics, semiconductor devices, electromagnetics, and astrophysics. Over the past decade multiderivative time-stepping have been increasingly used for the time-evolution hyperbolic PDEs, so that the strong stability properties of these methods have become important. In this work we review sufficient conditions for a two-derivative multistage method to preserve the strong stability properties of spatial discretizations in a forward Euler and different conditions on the second derivative. In particular we present the strong stability preserving theory for explicit and implicit two-derivative Runge–Kutta schemes, including a special condition on the second derivative under which these implicit methods may be unconditionally strong stability preserving. This special condition is natural for the stiff component of wide range of plasma physics problems, and can be useful in the context of strong stability preserving implicit-explicit multi-derivative Runge–Kutta schemes, where the time-step restriction is then independent of the stiff term. Lastly, we present the strong stability preserving theory for implicit-explicit multi-derivative general linear methods, and some novel second and third order methods where the time-step restriction is independent of the stiff term.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

On a fully-implicit VMS-stabilized FE formulation for low Mach number compressible resistive MHD with application to MCF

This study presents the development and evaluation of a fully-implicit variational multiscale (VMS) stabilized unstructured finite element (FE) formulation for compressible magnetohydrodynamics (MHD) model, at low Mach number regime. The model describes the dynamics of a compressible conducting fluid in the low Mach number limit in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics phenomena. The specific applications that motivate this study are macroscopic simulations of the longer time-scale stability and disruptions of magnetic confinement fusion (MCF) devices, specifically the ITER tokamak. The discussion considers the development of the VMS FE representation, the structure of the stabilizing terms that deal with significant convective flows, the stabilization of the nearly incompressible response of the fluid flow, and the stabilization of the constraint that enforces the solenoidal involution on the magnetic field. The nonlinear discretized system is solved with scalable preconditioned Newton–Krylov iterative methods, which employs a multiphysics block preconditioning method based on approximate block factorizations and Schur complements. The study presents an evaluation of the VMS method on a 2D cartesian tearing mode instability, and illustrates the scalability of the solvers on MCF relevant problems. A set of results are also presented for longer time-scale stability and disruptions for the ITER tokamak. These include a vertical displacement event (VDE), and a (1,1) internal kink mode. Here, the formulation is demonstrated to be scalable and also reasonably robust with respect to the Lundquist number scaling.

42 ENGINEERING↗

DOE OSTI2020

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

TokaMaker: An open-source time-dependent Grad-Shafranov tool for the design and modeling of axisymmetric fusion devices

Here, in this paper, we present a new static and time-dependent MagnetoHydroDynamic (MHD) equilibrium code, TokaMaker, for axisymmetric configurations of magnetized plasmas, based on the well-known Grad-Shafranov equation. This code utilizes finite element methods on an unstructured triangular grid to enable capturing accurate machine geometry and simple mesh generation from engineering-like descriptions of present and future devices. The new code is designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. A detailed description of the numerical methods of the code, including a novel formulation of the boundary conditions for free-boundary equilibria, and validation of the implementation of those methods using both analytic test cases and cross-code validation is shown. Results show expected convergence across tested polynomial degree for analytic and cross-code test cases.

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↗

DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗

DOE OSTI2026

Fast solvers for tokamak fluid models with PETSc

Multigrid (MG) is widely recognized as a highly effective solver for the model problem, the Laplacian, but textbook MG fails on most problems of interest. MG methods have been applied to complex, real-world applications with careful consideration of the physical model and discretization. In this work we develop the first step in applying MG methods to science and engineering relevant magnetohydrodynamics (MHD) tokamak models in the M3D-C1 (https://m3dc1.pppl.gov) fusion energy science code. The semi-implicit time integrator in M3D-C1 is composed of many linear solves. The implicit advance of the momentum equation is the most challenging and is the focus of this work. The current production solver in M3D-C1 is a block Jacobi (BJ) preconditioner within a Krylov solver, where blocks group degrees of freedom on planes of constant toroidal coordinate. BJ convergence degrades as the number of planes increases due to the spectral properties of the matrix preconditioned with BJ. The partially magnetic field-aligned, regular toroidal grid structure in M3D-C1 is amenable to semi-coarsening geometric MG in the toroidal direction. This paper develops such a solver and demonstrates competitive performance on a runaway electron model of a SPARC (https://cfs.energy/technology/sparc) disruption, and superior robustness on a stellarator model on which the BJ solver fails to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

CAKE: Consistent Automatic Kinetic Equilibrium reconstruction

Magnetohydrodynamic equilibrium reconstructions with kinetic constraints are an essential input of many tokamak stability and transport analysis work- flows. The reconstructions themselves also contain useful information about pressure and current density distributions within the plasma which can be used to drive other calculations. While standardized magnetic-only, and partially kinetic equilibrium reconstruction is widely available and used, re- constructions for well diagnosed machines such as DIII-D, that incorporates a large and heterogeneous set of internal diagnostic measurements yet robustly reaching high numerical accuracy have been a time intensive, manual process. The Consistent Automatic Kinetic Equilibrium reconstruction tool (CAKE) has been developed and implemented at DIII-D for producing low error, kinetically constrained reconstructions without human intervention, and to identify strategies for doing so that would be broadly applicable. Here, the techniques employed to allow CAKE to handle data of limited quality and availability, as well as the full automation of the workflow may assist future efforts to form kinetic constraint profiles in real time plasma control applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Steady state thermo-mechanics and material property definition framework for analyzing DCLL blanket in the fusion nuclear science facility

In this work, a thermo-mechanics model that relies on creating the material property definition framework (MPDF) and multiphysics coupling of the heat transfer and the solid mechanics modules is developed to determine the structural integrity of the recently designed dual cooled lead lithium (DCLL) inboard blanket (IB) for the Fusion Nuclear Science Facility under steady state loads. The MPDF is called to supply fusion relevant neutron irradiation and temperature induced changes in material properties during multiphysics finite element runs, and PbLi temperature profiles are used to approximate Magnetohydrodynamics effect and the nuclear volumetric heating on the PbLi. Neutron irradiation and temperature induced reduction of the yield and ultimate strengths of F82H steel at the first wall (FW) are quantified for one year. A blanket in an assembly with gaps between blanket sectors and another blanket in an assembly with no gaps between blanket sectors, both exposed to radiation damage that lasted for one year are analyzed. Analysis using the elastic ITER structural design criteria for in-vessel components (ITER SDC-IC) design rules and a linear isotropic-hardening-type elastoplastic material model are used where most appropriate. The IB blanket with gaps between blanket sectors will withstand the steady state combined thermal and coolant loads for one year operational period but will fail if no gaps are allowed between blanket sectors. It is recommended that a gap of about 7.62 mm should be provided between IB blanket sectors during assembly which would close up during service, stop neutron streaming, reduce stresses and reduce bending of the FW into the scrape-off layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Model predictive current profile control in tokamaks by exploiting spatially moving electron cyclotron current drives

The plasma control systems in next-generation tokamaks like ITER will balance competing control objectives to achieve the desired level of performance in advanced scenarios while preventing magnetohydrodynamic instabilities and disruptions. During normal tokamak operation, the points of incidence of the electromagnetic waves generated by the electron cyclotron heating and current drives (EC H&CDs) are usually fixed in space. However, the points of incidence can be modified in real-time by changing the angles of the mirrors that reflect the EC H&CD waves. Altering the points of incidence, in turn, varies the ability of the plasma control system to regulate a plasma property. For instance, changing the EC H&CD wave incidence location may place the power demands necessary to achieve a particular plasma target within saturation limits. Therefore, using the EC H&CD deposition location, which is related to the EC H&CD mirror angle, as a supplementary controllable variable may facilitate access to a given target scenario. However, active scenario-control algorithms have not been designed so far to fully exploit this capability in real time. In this work, a model predictive controller that can handle actuation locations as control inputs is developed. In particular, the controller is designed to regulate both the auxiliary powers and the EC H&CD deposition locations in a pre-defined optimal sense to achieve the control objective of attaining and sustaining a target current profile. Further, the proposed controller is tested for a DIII-D tokamak scenario in nonlinear simulations using the Control Oriented Transport SIMulator (COTSIM).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Nonlinear control of safety factor gradient in tokamaks using spatially variable electron cyclotron current drives

Active control of plasma properties may be necessary to achieve stable operation of next-generation tokamaks over large time scales. Such control algorithms can regulate the plasma properties to avoid the onset of magnetohydrodynamic (MHD) instabilities. For instance, the global and local properties of the safety factor profile are linked to the onset of neoclassical tearing modes (NTMs). This work proposes a model-based control approach for deterring/delaying NTMs through active modulation of one of the safety factor properties - the gradient of the safety factor profile at a particular rational safety factor surface. In particular, a novel control-oriented model for the local safety factor gradient is developed. The nonlinear control model is governed by a nonautonomous ordinary differential equation that accounts for a given rational safety factor surface’s spatial variation over time. Further, to improve the controllability of the spatially evolving parameter, the control model treats ECH&CD positions, along with noninductive powers, as controllable variables. A nonlinear control algorithm based on feedback linearization with optimization is synthesized to achieve the objective of regulating the safety factor gradient around a given target. The proposed algorithm allocates optimal ECH&CD positions, in addition to auxiliary powers, at each time instant as the rational safety factor surface drifts to locations with a low control authority. Stability guarantees of the proposed control law are also discussed in this work. The proposed algorithm is tested for a DIII-D tokamak scenario in nonlinear simulations carried out using the Control Oriented Transport SIMulator (COTSIM). Both fixed and moving ECH&CD cases are studied, and their outcomes are compared. Simulation results demonstrate that enthusiastic regulation of the safety factor gradient can be achieved during the ramp-up and flat-top phases of tokamak operation in both fixed and moving ECH&CD cases. However, real-time updates of ECH&CD positions can prevent the saturation of auxiliary powers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Nonlinear control of the minimum safety factor in tokamaks by optimal allocation of spatially moving electron cyclotron current drive

The minimum value of the safety factor profile is related to the magnetohydrodynamic (MHD) stability of the plasma confined in a tokamak. Therefore, active control of the minimum safety factor may mitigate MHD instabilities that can degrade or even terminate plasma confinement. Typically, in most tokamak scenarios, the minimum safety factor evolves spatially with time, i.e., the location at which the safety factor achieves the minimum value changes with time. In addition to the inherent nonlinearities in the minimum safety factor evolution, its spatial variation makes the control design challenging. In particular, complexity in control design may arise from the need for time-dependent nonlinear models that account for spatial variation of the minimum safety factor. Furthermore, the minimum safety factor may drift to locations where the actuator authority is low. The problem of minimum safety factor control with target location tracking and moving electron cyclotron current drive (ECCD) is addressed in this work. A nonlinear time-dependent model that incorporates the spatial variation of the minimum safety factor is presented. A nonlinear controller based on optimal feedback linearization is developed to track a target minimum safety factor. The proposed controller treats the ECCD position as a controllable variable. In other words, the controller prescribes the ECCD position (in addition to the non-inductive powers) in real time based on an optimal criterion that is defined a priori. This work also presents the steps necessary to integrate the minimum safety factor controller with a total energy controller to achieve multiple control objectives simultaneously. In conclusion, the proposed integrated control algorithm is tested using nonlinear simulations in the Control Oriented Transport SIMulator (COTSIM) for a DIII-D tokamak scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

ALEGRA: Finite element modeling for shock hydrodynamics and multiphysics

ALEGRA is a multiphysics finite-element shock hydrodynamics code, under development at Sandia National Laboratories since 1990. Fully coupled multiphysics capabilities include transient magnetics, magnetohydrodynamics, electromechanics, and radiation transport. Importantly, ALEGRA is used to study hypervelocity impact, pulsed power devices, and radiation effects. The breadth of physics represented in ALEGRA is outlined here, along with simulated results for a selected hypervelocity impact experiment.

42 ENGINEERING↗