Search NASA⌕ Search

SEARCH · Search NASA

Results for “magnetohydrodynamics (MHD)”

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
NASA NTRS1981

Parametric analysis of closed cycle magnetohydrodynamic (MHD) power plants

A parametric analysis of closed cycle MHD power plants was performed which studied the technical feasibility, associated capital cost, and cost of electricity for the direct combustion of coal or coal derived fuel. Three reference plants, differing primarily in the method of coal conversion utilized, were defined. Reference Plant 1 used direct coal fired combustion while Reference Plants 2 and 3 employed on site integrated gasifiers. Reference Plant 2 used a pressurized gasifier while Reference Plant 3 used a ""state of the art' atmospheric gasifier. Thirty plant configurations were considered by using parametric variations from the Reference Plants. Parametric variations include the type of coal (Montana Rosebud or Illinois No. 6), clean up systems (hot or cold gas clean up), on or two stage atmospheric or pressurized direct fired coal combustors, and six different gasifier systems. Plant sizes ranged from 100 to 1000 MWe. Overall plant performance was calculated using two methodologies. In one task, the channel performance was assumed and the MHD topping cycle efficiencies were based on the assumed values. A second task involved rigorous calculations of channel performance (enthalpy extraction, isentropic efficiency and generator output) that verified the original (task one) assumptions. Closed cycle MHD capital costs were estimated for the task one plants; task two cost estimates were made for the channel and magnet only.

Owens, W.↗
NASA NTRS1984

A full-implicit-continuous-Eulerian (FICE) scheme for multidimensional transient magnetohydrodynamic (MHD) flows

A full implicit continuous Eulerian (FICE) scheme is developed for solving multidimensional transient MHD flow problems. The physical system under consideration is a general case of a transient MHD flow in which an initial steady state is subject to a finite amplitude disturbance. The governing equations are described, their finite difference formulation is presented and the FICE algorithm is given. The boundary conditions are treated by classifying them into physical and computational ones. The usefulness of the FICE algorithm is demonstrated using a physical example concerning the dynamical response of the static solar atmosphere due to a representative photospheric disturbance.

Wu, S. T.↗
NASA NTRS1986

Magnetohydrodynamics (MHD) modelling of flare energy buildup, the energy release phase, and its propagation into heliospheric space

Solar flare energy buildup at the photospheric level and energy release and transport into heliospheric space are examined using a composite MHD model. A four phase composite MHD model is described. An example demonstrating the applicability of the model is presented; the model was applied to the active region AR 2372. The limitations of this composite MHD model approach to analyzing solar flare energy buildup are discussed.

Wu, S. T.↗
NASA NTRS1988

Magnetohydrodynamic (MHD) modelling of solar active phenomena via numerical methods

Numerical ideal MHD models for the study of solar active phenomena are summarized. Particular attention is given to the following physical phenomena: (1) local heating of a coronal loop in an isothermal and stratified atmosphere, and (2) the coronal dynamic responses due to magnetic field movement. The results suggest that local heating of a magnetic loop will lead to the enhancement of the density of the neighboring loops through MHD wave compression. It is noted that field lines can be pinched off and may form a self-contained magnetized plasma blob that may move outward into interplanetary space.

Wu, S. T.↗
NASA NTRS1991

Shear-induced instability and arch filament eruption - A magnetohydrodynamic (MHD) numerical simulation

A situation wherein a bipolar magnetic field embedded in a stratified solar atmosphere undergoes symmetrical shear motion at the footpoints is investigated via a 2D (nonplanar) MHD simulation. It was found that the vertical plasma flow velocities grow exponentially, leading to a new type of global MHD instability. The growth rate increases almost linearly until it reaches the same order of magnitude as the Alfven speed. Then a nonlinear MHD instability occurs beyond this point. It was found that the central loops are pinched by opposing Lorentz forces, and the outer closed loops stretch upward with the vertically-rising mass flow. The nonlinear dynamical shearing instability is illustrated by a numerical example that is given for three different values of the plasma beta that span several orders of magnitude.

Wu, S. T.↗
NASA NTRS2008

Recovering Photospheric Velocities from Vector Magnetograms by using a Three-Dimensional, Fully Magnetohydrodynamic (MHD) Model

We introduce a numerical simulation method for recovering the photospheric velocity field from the vector magnetograms. The traditional method is local correlation tracking (LCT) which is based on measuring the relative displacements of features in blocks of pixels between successive white-light images or magnetograms. Within this method, there are a variety of implementations. One of recently developed implementations is induction local correlation tracking (ILCT) as described by Welsch et al. (2004). They employ the normal component of magnetic induction equation as a constraint to assure consistent solutions. Our numerical method uses the fully three-dimensional MHD equations to recover the photospheric velocity field with individual vector magnetograms. We compare our method to the ILCT method using NOAA AR8210 as an example. The differences and similarities are discussed in detail.

Wang, A. H.↗
DOE OSTI2022

An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD

The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Additionally, therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (mfem.org). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in Chacón et al. (2002) to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime (≥ 10 7 ) with the successful resolution of plasmoid instabilities and thin current sheets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Identification of Kelvin-Helmholtz generated vortices in magnetised fluids

The Kelvin-Helmholtz Instability (KHI), arising from velocity shear across the magnetopause, plays a significant role in the viscous-like transfer of mass, momentum, and energy from the shocked solar wind into the magnetosphere. While the KHI leads to growth of surface waves and vortices, suitable detection methods for these applicable to magnetohydrodynamics (MHD) are currently lacking. A novel method is derived based on the well-established λ-family of hydrodynamic vortex identification techniques, which define a vortex as a local minimum in an adapted pressure field. The J × B Lorentz force is incorporated into this method by using an effective total pressure in MHD, including both magnetic pressure and a pressure-like part of the magnetic tension derived from a Helmholtz decomposition. The λ MHD method is shown to comprise of four physical effects: vortical momentum, density gradients, fluid compressibility, and the rotational part of the magnetic tension. A local three-dimensional MHD simulation representative of near-flank magnetopause conditions (plasma β’s 0.5 – 5 and convective Mach numbers M f ∼ 0.4) under northward interplanetary magnetic field (IMF) is used to validate λ MHD . Analysis shows it correlates well with hydrodynamic vortex definitions, though the level of correlation decreases with vortex evolution. Overall, vortical momentum dominates λ MHD at all times. During the linear growth phase, density gradients act to oppose vortex formation. By the highly nonlinear stage, the formation of small-scale structures leads to a rising importance of the magnetic tension. Compressibility was found to be insignificant throughout. Finally, a demonstration of this method adapted to tetrahedral spacecraft observations is performed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Cholla-MHD: An Exascale-capable Magnetohydrodynamic Extension to the Cholla Astrophysical Simulation Code

Abstract We present an extension of the massively parallel, GPU native, astrophysical hydrodynamics code Cholla to magnetohydrodynamics (MHD). Cholla solves the ideal MHD equations in their Eulerian form on a static Cartesian mesh utilizing the Van Leer + constrained transport integrator, the HLLD Riemann solver, and reconstruction methods at second and third order. Cholla’s MHD module can perform ≈260 million cell updates per GPU-second on an NVIDIA A100 while using the HLLD Riemann solver and second order reconstruction. The inherently parallel nature of GPUs combined with increased memory in new hardware allows Cholla’s MHD module to perform simulations with resolutions ∼500 3 cells on a single high-end GPU (e.g., an NVIDIA A100 with 80 GB of memory). We employ GPU direct Message Passing Interface to attain excellent weak scaling on the exascale supercomputer Frontier, while using 74,088 GPUs and simulating a total grid size of over 7.2 trillion cells. A suite of test problems highlights the accuracy of Cholla’s MHD module and demonstrates that zero magnetic divergence in solutions is maintained to round off error. We also present new testing and CI tools using GoogleTest, GitHub Actions, and Jenkins that have made development more robust and accurate and ensure reliability in the future.

Astronomy & Astrophysics↗
NASA NTRS1984

Sialon Electrodes and Insulators for MHD Device

Rectangular magnetohydrodynamic (MHD) channel structure for electrical power generation designed using pure sialon ceramic for insulating portion of structure and metal-bearing sialon cermet for conducting portion.

Phillips, W. M.↗
NASA NTRS2010

Broken Ergodicity in MHD Turbulence

Ideal magnetohydrodynamic (MHD) turbulence may be represented by finite Fourier series, where the inherent periodic box serves as a surrogate for a bounded astrophysical plasma. Independent Fourier coefficients form a canonical ensemble described by a Gaussian probability density function containing a Hermitian covariance matrix with positive eigenvalues. The eigenvalues at lowest wave number can be very small, resulting in a large-scale coherent structure: a turbulent dynamo. This is seen in computations and a theoretical explanation in terms of 'broken ergodicity' contains Taylor s theory of force-free states. An important problem for future work is the case of real, i.e., dissipative flows. In real flows, broken ergodicity and coherent structure are still expected to occur in MHD turbulence at the largest scale, as suggested by low resolution simulations. One challenge is to incorporate coherent structure at the largest scale into the theory of turbulent fluctuations at smaller scales.

Shebalin, John V.↗
DOE OSTI2023

Ideal MHD induced temperature flattening in spherical tokamaks

This paper extends the analysis first presented in Jardin et al. [Phys. Rev. Lett. 128, 245001 (2022)] to more thoroughly examine the stability of spherical torus equilibrium to ideal magnetohydrodynamic (MHD) infernal modes and their nonlinear consequences. Here we demonstrate that in a 3D resistive magnetohydrodynamic (MHD) simulation of a NSTX discharge, anomalous transport can occur due to these instabilities. Here we generate a family of equilibrium of differing β and use this to show that these instabilities could explain the experimentally observed flattening of the electron temperature profile at modest β. The modes studied in this paper are found to occur with poloidal mode number m and toroidal mode number n when the ratio m/ n is in the range of 1.2–1.5, when the central safety factor is in this range or slightly lower, and when the central region has very low magnetic shear. Our analysis gives some insight as to why the unstable linear growth rates are oscillatory functions of the toroidal mode number n. We present a simulation of an initially stable configuration that passes through a stability boundary at a critical β as it is heated. We also show that a particular NSTX discharge is unstable to these modes over a timescale of several hundred ms. We conclude that these modes must be taken into account when performing predictive modeling. An appendix shows that similar modes can be found in R / a = 4 tokamaks for certain q-profiles and β values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

Solar-Driven Liquid-Metal MHD Generator

Liquid-metal magnetohydrodynamic (MHD) power generator with solar oven as its heat source has potential to produce electric power in space and on Earth at high efficiency. Generator focuses radiation from Sun to heat driving gas that pushes liquid metal past magnetic coil. Power is extracted directly from electric currents set up in conducting liquid. Using solar energy as fuel can save considerable costs and payload weight, compared to previous systems.

Hohl, F.↗
NASA NTRS1986

Laser-powered MHD generators for space application

Magnetohydrodynamic (MHD) energy conversion systems of the pulsed laser-supported detonation (LSD) wave, plasma MHD, and liquid-metal MHD (LMMHD) types are assessed for their potential as space-based laser-to-electrical power converters. These systems offer several advantages as energy converters relative to the present chemical, nuclear, and solar devices, including high conversion efficiency, simple design, high-temperature operation, high power density, and high reliability. Of these systems, the Brayton cycle liquid-metal MHD system appears to be the most attractive. The LMMHD technology base is well established for terrestrial applications, particularly with regard to the generator, mixer, and other system components. However, further research is required to extend this technology base to space applications and to establish the technology required to couple the laser energy into the system most efficiently. Continued research on each of the three system types is recommended.

Jalufka, N. W.↗
NASA NTRS1997

Global Magnetohydrodynamic Modeling of the Solar Corona

Under this contract, we have continued our investigations of the large scale structure of the solar corona and inner heliosphere using global magnetohydrodynamic (MHD) simulations. These computations have also formed the basis for studies of coronal mass ejections (CMES) using realistic coronal configurations. We have developed a technique for computing realistic magnetohydrodynamic (MHD) computations of the solar corona and inner heliosphere. To perform computations that can be compared with specific observations, it is necessary to incorporate solar observations into the boundary conditions. We have used the Wilcox Solar Observatory synoptic maps (collected during a solar rotation by daily measurements of the line-of-sight magnetic field at central meridian) to specify the radial magnetic field (B,) at the photosphere. For the initial condition, we use a potential magnetic field consistent with the specified distribution of B, at the lower boundary, and a wind solution consistent with the specified plasma density and temperature at the solar surface. Together this initial condition forms a (non-equilibrium) approximation of the state of the solar corona for the time-dependent MHD computation. The MHD equations are then integrated in time to steady state. Here we describe solutions relevant to a recent solar eclipse, as well as Ulysses observations. We have also developed a model configuration of solar minimum, useful for studying CME initiation and propagation.

Linker, Jon A.↗
DOE OSTI2022

Numerical Analysis of Liquid Metal MHD Flow and Heat Transfer for Open-Surface Li Divertor in FNSF

Within the ongoing U.S.-based program on the development of liquid metal plasma-facing components, numerical simulations and analyses are performed to address the feasibility of the open-surface Li divertor. In the previous scoping studies (Smolentsev, 2021), heat-removal capabilities of the divertor were assessed using a simplified flow model for a slug-type velocity profile and constant flow thickness. Here, new analyses take into account forces acting on the flowing Li layer. Three reduced-order mathematical models are applied under the conditions of the U.S. Fusion Nuclear Science Facility (FNSF) to access magnetohydrodynamic (MHD) flow development effects, velocity distribution, and surface waves: 1) fully developed MHD flow; 2) quasi-2-D developing MHD flow; and 3) multiphase MHD flow. The obtained results for MHD flows and the surface heat flux computed with the plasma code scrape-off layer plasma simulation for ITER (SOLPS-ITER) are then used as input data to compute the temperature distribution in the divertor by solving the convection–diffusion energy equation.

Smolentsev, Sergey↗
NASA NTRS1994

MHD simulations of the magnetosheath

Ideal magnetohydrodynamic (MHD) simulations of the dayside magnetosheath are reported. The calculations show a depletion layer near the magnetopause and a new stagnation flow structure. The stagnation flow appears to be a generalization of the stagnation line for ideal MHD. For Mach numbers typical of the solar wind at the Earth, there is no evidence of the slow mode waves recently observed, but for low Mach numbers prominent slow mode waves are seen in the magnetosheath.

Lyon, J. G.↗
NASA NTRS2009

Cavities of Weak Magnetic Field Strength in the Wake of FTEs: Results from Global Magnetospheric MHD Simulations

We use the global magnetohydrodynamic (MHD) code BATS-R-US to model multipoint observations of Flux Transfer Event (FTE) signatures. Simulations with high spatial and temporal resolution predict that cavities of weak magnetic field strength protruding into the magnetosphere trail FTEs. These predictions are consistent with recently reported multi-point Cluster observations of traveling magnetopause erosion regions (TMERs).

Kuznetsova, M. M.↗