Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Pressure drop in a prototypical 3D magnetohydrodynamic flow across contraction of a fusion blanket manifold

Predicting 3D magnetohydrodynamic (MHD) pressure losses associated with liquid metal flows in complex geometry ducts is required to design breeding blankets for fusion reactors. Of such components, manifolds exhibit major pressure losses. Recent correlations [T. Rhodeset al. Magnetohydrodynamic pressure drop and flow balancing of liquid metal flow in a prototypic fusion blanket manifold. Phys. Fluids. 2018; 30: 057101.] demonstrate promise for predicting pressure drops in electrically insulated inlet manifolds. In the present work, the extend to which these correlations can be applied to outlet manifolds, which feature sudden contractions, is investigated in both viscous-electromagnetic (VE) and inertial-electromagnetic (IE) regimes for Reynolds numbers 50 < Re < 1500, Hartmann numbers 2500 < Ha < 5475, and expansion/contraction ratio 10. Numerical computations have shown that the MHD pressure drops in the eutectic lead-lithium (PbLi) flows in contractions are almost identical to those in the flows featuring expansion with slightly higher magnitudes in the contraction cases. The small discripancy in the MHD pressure drop between contractions and expansions (<8%) suggests that the earlier obtained correlations for the 3D MHD pressure drop in a duct flow with a sudden expansion can also be applied to flows with a sudden contraction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

iharm3D: Vectorized General Relativistic Magnetohydrodynamics

iharm3D is an open-source C code for simulating black hole accretion systems in arbitrary stationary spacetimes using ideal general-relativistic magnetohydrodynamics (GRMHD). It is an implementation of the HARM (“High Accuracy Relativistic Magnetohydrodynamics”) algorithm outlined in Gammie et al. (2003) with updates as outlined in McKinney & Gammie (2004) and Noble et al. (2006). The code is most directly derived from Ryan et al. (2015) but with radiative transfer portions removed. HARM is a conservative finite-volume scheme for solving the equations of ideal GRMHD, a hyperbolic system of partial differential equations, on a logically Cartesian mesh in arbitrary coordinates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Interactions Between Microturbulence and Magnetohydrodynamic Instabilities in Magnetized Plasmas (Final Technical Report)

This document details the work accomplished during the period 04/15/2012-01/31/2020 for the collaborative project “Interactions Between Microturbulence and Magnetohydrodynamic Instabilities in Magnetized Plasmas,” and serves as the final report on the University of California, San Diego part of this collaborative project for Department of Energy Office of Science - Fusion Energy Sciences Grant No. DE- SC0007783. This project focused on the development of analytic and computational tools, seeking to better understand the mutual interactions between small-scale temperature-gradient driven microturbulence and macroscopic magnetohydrodynamic (MHD) instabilities, particularly tearing modes, in magnetically confined plasmas relevant for the generation of fusion energy. This project advanced understanding through numerical simulations of different sets of model equations, corresponding to different microturbulence instabilities. Questions addressed included how the magnitudes, spatial structures, and scaling's of turbulent forcing's and drags impact magnetic island dynamics, as well as quantifying the “region of influence” of an island on the microturbulence dynamics. A key finding of the project is that the microturbulence fluctuations and associated cross-field particle and energy fluxes exhibit spatial asymmetries for sufficiently large islands. The results have improved our fundamental understanding of plasma dynamics, motivated new experimental investigations, and will help inform future fusion predictive modeling efforts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Towards performance portability in the Spark astrophysical magnetohydrodynamics solver in the Flash-X simulation framework

Simulations of core-collapse supernovae, and other astrophysical phenomena, are quintessential extreme-scale computing challenges. For core-collapse supernova simulations to be carried out by the ExaStar project under the Exascale Computing Project umbrella, a robust, efficient, and state-of-the-art magnetohydrodynamics solver is a critical requirement. In Flash-X, the primary software instrument for ExaStar, a new magnetohydrodynamics solver has been designed and implemented from the ground up to achieve accuracy and efficiency for simulations of complex astrophysical flows. This new solver, dubbed Spark, uses high-order spatial reconstruction, Runge-Kutta time integration, and an efficient cell-centered approach to satisfying the divergence-free condition for the magnetic fields. Spark was written to be optimized for data locality in cache hierarchy of CPUs. Since data locality optimizations for cache hierarchy are not directly compatible with those of accelerators, we have taken the approach of using program synthesis to avoid massive amounts of code replication that would be necessary if we were to maintain two different versions of the solver. Our program synthesis relies on a simple key-dictionary approach, implemented in python, that enables us to assemble the version of the solver suitable for the target hardware from code fragments identified by specific keys. In this work, we describe the data locality optimizations of the solver for CPUs and accelerators and the program synthesis tools that enable this portability. We also detail the parallel performance of Spark for both CPUs and accelerators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A general metriplectic framework with application to dissipative extended magnetohydrodynamics

General equations for conservative yet dissipative (entropy producing) extended magnetohydrodynamics are derived from two-fluid theory. Keeping all terms generates unusual cross-effects, such as thermophoresis and a current viscosity that mixes with the usual velocity viscosity. While the Poisson bracket of the ideal version of this model has already been discovered, we determine its metriplectic counterpart that describes the dissipation. This is done using a new and general thermodynamic point of view to derive dissipative brackets, a means of derivation that is natural for understanding and creating dissipative dynamics without appealing to underlying kinetic theory orderings. Finally, the formalism is used to study dissipation in the Lagrangian variable picture where, in the context of extended magnetohydrodynamics, non-local dissipative brackets naturally emerge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Time-dependent relaxed magnetohydrodynamics: Inclusion of cross helicity constraint using phase-space action

In this work, a phase-space version of the ideal magnetohydrodynamic (MHD) Lagrangian is derived from first principles and shown to give a relabeling transformation when a cross-helicity constraint is added in Hamilton's Action Principle. A new formulation of time-dependent “relaxed” magnetohydrodynamics is derived using microscopic conservation of mass and macroscopic constraints on total magnetic helicity, cross helicity, and entropy under variations of density, pressure, fluid velocity, and magnetic vector potential. This gives Euler–Lagrange equations consistent with previous work on both ideal and relaxed MHD equilibria with flow, but generalizes the relaxation concept from statics to dynamics. The application of the new dynamical formalism is illustrated for short-wavelength linear waves, and the interface connection conditions for Multiregion Relaxed MHD (MRxMHD) are derived. The issue of whether E+u×B=0 should be a constraint is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Axial magnetic flux amplification in Hall-magnetohydrodynamic simulations of externally magnetized z-pinches

It is shown that externally magnetized gas-puff z-pinches under compression can exhibit robust axial flux amplification under suitable axial boundary conditions. This effect relies upon the Hall term in the generalized Ohm's law to generate azimuthal currents in the presence of a driving axial electric field. Under dynamic compression, the total current tends to flow in a mostly force-free boundary layer, separating the predominately azimuthal field outside the layer and the predominately axial field inside the layer. The effect only appears to occur if the axial boundaries allow for outflow or absorption of electromagnetic energy. The effect is mitigated by imposing either periodic or conducting axial boundary conditions. A semi-analytic equilibrium analysis agrees with steady-state solutions of the time-dependent electron-magnetohydrodynamic equations and provides an estimate of the scaling of the boundary layer as well as suggesting a scenario for the formation of the boundary layer. When operative, the effect can significantly impede plasma compression due to the increase in axial flux that diffuses into the pre-compressed plasma or through the presence of conductors that inhibit movement of the generated axial flux. Several facilities have noted unusual implosion behavior in z-pinch experiments with applied axial fields that does not appear to be explainable within the standard magnetohydrodynamic model. In conclusion, it is suggested that these experiments can be explained by the axial flux amplification and concentration effect

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Topanga: A kinetic ion plasma code for large-scale ionospheric simulations on magnetohydrodynamic timescales

Topanga is a kinetic ion code developed for simulating large-scale plasma phenomena in the Earth's ionosphere on magnetohydrodynamic timescales. It is a domain-decomposed parallel code that runs on high-performance computing platforms. Features of Topanga include spherical geometry for simplified boundary conditions and computational efficiency; a hybrid plasma model with inertia-less fluid electrons, kinetic ions, and an electric field specified via an Ohm's law; a Maxwell-FDTD (finite difference time domain) plasma model which retains the displacement current in Maxwell's equations and models electron currents in the ionosphere with a tensor conductivity; sponge-layer boundary conditions for absorption of electromagnetic and plasma waves incident on the domain boundaries; and a novel mixed-implicit algorithm for evolving the EM fields inside the Maxwell-FDTD region that is stable over many orders of magnitude in the electron–ion collision frequency. We verify the numerical methods used in Topanga on a pair of test problems. The first test involves modeling a three-dimensional collisionless shock using the hybrid set of equations. The second test involves modeling a spherical TEM mode in vacuum using the Maxwell-FDTD set of equations. Finally, we demonstrate how using the combined set of hybrid and Maxwell-FDTD equations to model the Starfish Prime high-altitude nuclear test recovers a “missing” EM signal on the ground that is not present when using only the hybrid set of equations. The magnitude of this signal in the simulation containing the Maxwell-FDTD region agrees well with the E3a portion of the magnetohydrodynamic electromagnetic pulse from Starfish Prime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Weak magnetohydrodynamic turbulence theory revisited

Two recent papers, P. H. Yoon and G. Choe, Phys. Plasmas 28, 082306 (2021) and Yoon et al., Phys. Plasmas 29, 112303 (2022), utilized in the derivation of the kinetic equation for the intensity of turbulent fluctuations the assumption that the wave spectra are isotropic, that is, the ensemble-averaged magnetic field tensorial fluctuation intensity is given by the isotropic diagonal form, ⟨δBiδBj⟩k=⟨δB2⟩kδij. However, it is more appropriate to describe the incompressible magnetohydrodynamic turbulence involving shear Alfvénic waves by modeling the turbulence spectrum as being anisotropic. That is, the tensorial fluctuation intensity should be different in diagonal elements across and along the direction of the wave vector, ⟨δBiδBj⟩k=12 ⟨δB⊥2⟩k(δij−kikj/k2)+⟨δB∥2⟩k(kikj/k2). In the present paper, we thus reformulate the weak magnetohydrodynamic turbulence theory under the assumption of anisotropy and work out the form of nonlinear wave kinetic equation.

Physics↗
DOE OSTI2022

Modelling general-relativistic plasmas with collisionless moments and dissipative two-fluid magnetohydrodynamics

ABSTRACT Relativistic plasmas are central to the study of black hole accretion, jet physics, neutron star mergers, and compact object magnetospheres. Despite the need to accurately capture the dynamics of these plasmas and the implications for relativistic transients, their fluid modelling is typically done using a number of (overly) simplifying assumptions, which do not hold in general. This is especially true when the mean free path in the plasma is large compared to the system size, and kinetic effects start to become important. Going beyond common approaches used in the literature, we describe a fully relativistic covariant 14-moment based two-fluid system appropriate for the study of electron–ion or electron–positron plasmas. This generalized Israel-Stewart-like system of equations of motion is obtained directly from the relativistic Boltzmann–Vlasov equation. This new formulation can account for non-ideal effects, such as anisotropic pressures and heat fluxes, not present in previous formulations of two-fluid magnetohydrodynamics. We show that a relativistic two-fluid plasma can be recast as a single fluid coupled to electromagnetic fields with (potentially large) out-of-equilibrium corrections. We keep all electron degrees of freedom, which provide self-consistent evolution equations for electron temperature and momentum. The out-of-equilibrium corrections take the form of a collisional 14-moment closure previously described in the context of viscous single fluids. The equations outlined in this paper are able to capture the full two-fluid character of collisionless plasmas found in black hole accretion and flaring processes around compact objects, as well Braginskii-like two-fluid magnetohydrodynamics applicable to weakly collisional plasmas inside accretion discs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Reconnection-driven energy cascade in magnetohydrodynamic turbulence

Magnetohydrodynamic turbulence regulates the transfer of energy from large to small scales in many astrophysical systems, including the solar atmosphere. We perform three-dimensional magnetohydrodynamic simulations with unprecedentedly large magnetic Reynolds number to reveal how rapid reconnection of magnetic field lines changes the classical paradigm of the turbulent energy cascade. By breaking elongated current sheets into chains of small magnetic flux ropes (or plasmoids), magnetic reconnection leads to a previously undiscovered range of energy cascade, where the rate of energy transfer is controlled by the growth rate of the plasmoids. As a consequence, the turbulent energy spectra steepen and attain a spectral index of -2.2 that is accompanied by changes in the anisotropy of turbulence eddies. The omnipresence of plasmoids and their consequences on, for example, solar coronal heating, can be further explored with current and future spacecraft and telescopes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1976

Dissipative, forced turbulence in two-dimensional magnetohydrodynamics

The equations of motion for turbulent two-dimensional magnetohydrodynamic flows are solved in the presence of finite viscosity and resistivity, for the case in which external forces (mechanical and/or magnetic) act on the fluid. The goal is to verify the existence of a magnetohydrodynamic dynamo effect which is represented mathematically by a substantial back-transfer of mean square vector potential to the longest allowed Fourier wavelengths. External forces consisting of a random part plus a fraction of the value at the previous time step are employed, after the manner of Lilly for the Navier-Stokes case. The regime explored is that for which the mechanical and magnetic Reynolds numbers are in the region of 100 to 1000. The conclusions are that mechanical forcing terms alone cannot lead to dynamo action, but that dynamo action can result from either magnetic forcing terms or from both mechanical and magnetic forcing terms simultaneously.

Fyfe, D.↗
NASA NTRS1981

Anisotropic magnetohydrodynamic turbulence in a strong external magnetic field

A strong external dc magnetic field introduces a basic anisotropy into incompressible magnetohydrodynamic turbulence. The modifications that this is likely to produce in the properties of the turbulence are explored for the high Reynolds number case. The conclusion is reached that the turbulent spectrum splits into two parts: an essentially two dimensional spectrum with both the velocity field and magnetic fluctuations perpendicular to the dc magnetic field, and a generally weaker and more nearly isotropic spectrum of Alfven waves. A minimal characterization of the spectral density tensors is given. Similarities to measurements from the Culham-Harwell Zeta pinch device and the UCLA Macrotor Tokamak are remarked upon, as are certain implications for the Belcher and Davis measurements of magnetohydrodynamic turbulence in the solar wind.

Montgomery, D.↗
NASA NTRS1981

Nonlinear magnetohydrodynamic stability

The computer code developed by Bauer et al. (1978) for the study of the magnetohydrodynamic equilibrium and stability of a plasma in toroidal geometry is extended so that the growth rates of instabilities may be estimated more accurately. The original code, which is based on the variational principle of ideal magnetohydrodynamics, is upgraded by the introduction of a nonlinear formula for the growth rate of an unstable mode which acts as a quantitative measure of instability that is important in estimating numerical errors. The revised code has been applied to the determination of the nonlinear saturation, ballooning modes and beta limits for tokamaks, stellarators and torsatrons.

Bauer, F.↗
NASA NTRS1987

Magnetohydrodynamic (MHD) analyses of various forms of activity and their propagation through helio spheric space

Theoretical and numerical modeling of solar activity and its effects on the solar atmosphere within the context of magnetohydrodynamics were examined. Specifically, the scientific objectives were concerned with the physical mechanisms for the flare energy build-up and subsequent release. In addition, transport of this energy to the corona and solar wind was also investigated. Well-posed, physically self-consistent, numerical simulation models that are based upon magnetohydrodynamics were sought. A systematic investigation of the basic processes that determine the macroscopic dynamic behavior of solar and heliospheric phenomena was conducted. A total of twenty-three articles were accepted and published in major journals. The major achievements are summarized.

Wu, S. T.↗
NASA NTRS1988

Magnetohydrodynamic boundary conditions for global models

Boundary conditions in the ionosphere and the upstream solar wind are important in determining the dynamics of global magnetohydrodynamic models of the magnetosphere. It is generally recognized that the orientation of the magnetic field in the upstream solar wind strongly modulates the rate of energy input into the magnetosphere by magnetic reconnection. However, other aspects of the upstream boundary conditions may determine whether the reconnection occurs in a patchy manner, as in flux transfer events, or in a global manner, as in the Paschmann et al. (1979) events. Ionospheric boundary conditions should also affect the reconnection process. For example, ionospheric line-tying can cause x-line motion in the outer magnetosphere. If it is assumed that auroras occur on field lines mapping to x-lines, then auroral motions are different than the local convective motion of the plasma in which they occur. Global magnetohydrodynamic models which incorporate both magnetospheric reconnection and ionospheric convection could be used to investigate the effect of reconnection and convection upon dayside and nightside auroral motions during the course of a magnetic substorm.

Forbes, T. G.↗
NASA NTRS1989

Time evolution from linear to nonlinear stages in magnetohydrodynamic parametric instabilities

The nonlinear evolution of the magnetohydrodynamic (MHD) parametric instability of wave fluctuations propagating along an unperturbed magnetic field is investigated. Both a magnetohydrodynamic perturbation-theoretical approach and a nonlinear MHD simulation are used. It is shown that high harmonic waves are rapidly excited by wave-wave coupling, and that the wave spectrum evolves from a state containing a small number of degrees of freedom in k space to one which contains a large number of degrees of freedom. It is found that the spectral evolution prior to nonlinear saturation is well described by the prturbation theory. During this stage, the ratio of the growth rate of the nth harmonic wave to the linear growth rate of the fundamental wave is n. The nonlinear saturation stage is characterized by a frequency shift of the fundamental wave that destroys the wave-wave resonance condition which, in turn, causes the wave amplitude to cease its growth.

Hoshino, M.↗
NASA NTRS1992

'Reduced' magnetohydrodynamics and minimum dissipation rates

It is demonstrated that all solutions of the equations of 'reduced' magnetohydrodynamics approach a uniform-current, zero-flow state for long times, given a constant wall electric field, uniform scalar viscosity and resistivity, and uniform mass density. This state is the state of minimum energy dissipation rate for these boundary conditions. No steady-state turbulence is possible. The result contrasts sharply with results for full three-dimensional magnetohydrodynamics before the reduction occurs.

Montgomery, David↗