SEARCH · Search NASA

Results for “Numerical 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 253 records · Page 14

All sources NASA NTRS DOE OSTI

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↗

NASA NTRS2002

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗

DOE OSTI2025

Ten-moment fluid modeling of the Weibel instability

We investigate the one-dimensional non-relativistic Weibel instability through the capture of anisotropic pressure tensor dynamics using an implicit 10-moment fluid model that employs the electromagnetic Darwin approximation. The results obtained from the 10-moment model are compared with an implicit particle-in-cell simulation. The linear growth rates obtained from the numerical simulations are in good agreement with the theoretical fluid and kinetic dispersion relations. The fluid dispersion relations are derived using Maxwell’s equations and the Darwin approximation. We also show that the magnetohydrodynamic approximation can be used to model the Weibel instability if one accounts for an anisotropic pressure tensor and unsteady terms in the generalised Ohm’s law. In addition, we develop a preliminary theory for the saturation magnetic field strength of the Weibel instability, showing good agreement with the numerical results.

Kuldinow, D. A. (ORCID:0000000319730196)↗

NASA NTRS1979

Numerical study of an explosion in a non-homogeneous medium with and without magnetic fields

A version of the two step Lax-Wendroff difference method with second order accuracy is used to seek solutions of the unsteady magnetohydrodynamic (MHD) equations for the study of an explosion in a non-homogeneous medium with and without magnetic fields. The explosion is initiated by introducing a finite amount of energy in a small volume of gas which leads to an instantaneous increase of gas temperature. Two types of magnetic field configurations (i.e. open and closed) are considered to illustrate the dependence and differences of magnetic and a non-magnetic flow motion. Numerical results show the strong dependence of induced MHD flow field on the magnetic field configuration and strength. The development and decay of the fast and slow MHD shock waves, as well as the ordinary gasdynamic shock waves, are also represented well by the present computer simulations. In conclusion, we have demonstrated that the numerical scheme we have used is a reliable one for studying 2-dimensional MHD problems in non-homogeneous medium.

Han, S. M.↗

DOE OSTI2024

Numerical minimization of neoclassical poloidal viscosity for supersonic equilibria in tokamak geometry

Extensive experimental evidence has shown that the presence of poloidal flow in tokamaks can dramatically improve transport properties. However, theory indicates that poloidal flows are damped by poloidal viscosity, thus necessitating external drivers, such as neutral beam injection or radio frequency heating. In this work, ideal magnetohydrodynamic equilibria are calculated via the FORTRAN code FLOW [Guazzotto et al., Phys. Plasmas 11, 604 (2004)] and a postprocessor is used to estimate the neoclassical poloidal viscosity. The equilibrium inputs, which correspond to intuitive physical quantities, are then numerically optimized to reduce a viscosity figure of merit. We present supersonic equilibria in tokamak geometry with minimized neoclassical poloidal viscosities for various velocity free function inputs, plasma aspect ratios, and collisionality regimes. Benchmarks are made against an analytic theory as well as a classical expression of poloidal viscosity. Numerical confirmation of the analytic theory is obtained in the high aspect ratio and high collisionality limit. Good agreement is also seen near the plasma core and edge, with discrepancies arising in the intermediate region. Outside of these limits, rotation input function profiles are found that provide ∼order of magnitude improvements over the analytic theory, with additional progress being made toward predictions for tokamak-relevant equilibria.

Physics↗

DOE OSTI2023

Transition from wall modes to multimodality in liquid gallium magnetoconvection

Coupled laboratory-numerical experiments of Rayleigh-Bénard convection in liquid gallium subject to a vertical magnetic field are presented. The experiments are carried out in two cylindrical containers with diameter-to-height aspect ratio Γ = 1.0 and 2.0 at varying thermal forcing (Rayleigh numbers 10 5 ≲ Ra ≲ 10 8 ) and magnetic field strength (Chandrasekhar numbers 0 ≲ Ch ≲ 3 × 10 5 ). Laboratory measurements and numerical simulations confirm that magnetoconvection in our finite cylindrical tanks onsets via nondrifting wall-attached modes, in good agreement with asymptotic predictions for a semi-infinite domain. With increasing supercriticality, the experimental and numerical thermal measurements and the numerical velocity data reveal transitions between wall mode states with different azimuthal mode numbers and between wall-dominated convection to wall and interior multimodality. These transitions are also reflected in the heat transfer data, which combined with previous studies connect onset to supercritical turbulent behaviors in liquid metal magnetoconvection over a large parameter space. The gross heat transfer behaviors between magnetoconvection and rotating convection in liquid metals are compared and discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

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

The overall goal of this project was to develop a better understanding between small-scale microturbulence and large-scale magnetohydrodynamic (MHD) instabilities, particularly magnetic islands and tearing modes, in fusion relevant magnetically confined plasmas. This project advanced understanding through numerical simulations of different model equations. Questions addressed included how the magnitudes, spatial structures, and scalings of turbulent forcings and drags impact island dynamics, as well as what the "region of influence" of a tearing mode on the turbulence dynamics. A key finding of the project is that the microturbulence (and associated transport) exhibits 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↗

NASA NTRS1980

An adaptive computation mesh for the solution of singular perturbation problems

In singular perturbation problems, control of zone size variation can affect the effort required to obtain accurate, numerical solutions of finite difference equations. The mesh is generated by the solution of potential equations. Numerical results for a singular perturbation problem in two dimensions are presented. The mesh was used in calculations of resistive magnetohydrodynamic flow in two dimensions.

Brackbill, J. U.↗

NASA NTRS1983

Long-time states of inverse cascades in the presence of a maximum length scale

It is shown numerically, both for the two-dimensional Navier-Stokes (guiding-center plasma) equations and for two-dimensional magnetohydrodynamics, that the long-time asymptotic state in a forced inverse-cascade situation is one in which the spectrum is completely dominated by its own fundamental. The growth continues until the fundamental is dissipatively limited by its own dissipation rate.

Hossain, M.↗

NASA NTRS1983

Long-time states of inverse cascades in the presence of a maximum length scale

Hossain, M.↗

NASA NTRS1999

Practical Aspects of Stabilized FEM Discretizations of Nonlinear Conservation Law Systems with Convex Extension

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the Galerkin least-squares (GLS) and the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the POE system. Central to the development of the simplified GLS and DG methods is the Degenerative Scaling Theorem which characterizes right symmetrizes of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobean matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler, Navier-Stokes, and magnetohydrodynamic (MHD) equations. Linear and nonlinear energy stability is proven for the simplified GLS and DG methods. Spatial convergence properties of the simplified GLS and DO methods are numerical evaluated via the computation of Ringleb flow on a sequence of successively refined triangulations. Finally, we consider a posteriori error estimates for the GLS and DG demoralization assuming error functionals related to the integrated lift and drag of a body. Sample calculations in 20 are shown to validate the theory and implementation.

Barth, Timothy↗

DOE OSTI2025

A scalable multidimensional fully implicit solver for Hall magnetohydrodynamics

We propose an optimally performant fully implicit algorithm for the Hall magnetohydrodynamics (HMHD) equations based on multigrid-preconditioned Jacobian-free Newton-Krylov methods. HMHD is a challenging system to solve numerically because it supports stiff fast dispersive waves. The preconditioner is formulated using an operator-split approximate block factorization (Schur complement), informed by physics insight. We use a vector-potential formulation (instead of a magnetic field one) to allow a clean segregation of the problematic $\nabla$ x $\nabla$ x operator in the electron Ohm's law subsystem. This segregation allows the formulation of an effective damped block-Jacobi smoother for multigrid. We demonstrate by analysis that our proposed block-Jacobi iteration is convergent and has the smoothing property. The resulting HMHD solver is verified linearly with wave propagation examples, and nonlinearly with the GEM challenge reconnection problem by comparison against another HMHD code. We demonstrate the excellent algorithmic and parallel performance of the algorithm up to 16384 MPI tasks in two dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Novel approach to general curvilinear coordinates for plasma fluid applications

In general geometry, plasma fluid equations include nonlinear geometric sources associated with fictitious forces, which pose significant challenges to computer simulations. We reformulate the plasma fluid hierarchy to rigorously preserve geometry and conservation properties critical to numerical simulations, while concealing the geometric sources. In their discrete form, the reformulated models conserve mass, angular momentum, and energy naturally, by simple analogy with the continuum equations. These conservation properties have minimal requirements in discrete space, namely, the anti-symmetry of the first derivative and the orthogonality of the scalar and cross products. By decoupling magnetic geometry, coordinate systems, and numerical discretization, this enables maximum flexibility while preserving physics fidelity. As a testbed, we apply the novel representation to the resistive magnetohydrodynamic system, which involves a complete set of curvilinear operations. We verify the correctness of the approach using steady state liquid metal flows and the classic Orszag–Tang vortex.

Halpern, Federico D. [General Atomics, San Diego, ↗

NASA NTRS2007

Broken Symmetries and Magnetic Dynamos

Phase space symmetries inherent in the statistical theory of ideal magnetohydrodynamic (MHD) turbulence are known to be broken dynamically to produce large-scale coherent magnetic structure. Here, results of a numerical study of decaying MHD turbulence are presented that show large-scale coherent structure also arises and persists in the presence of dissipation. Dynamically broken symmetries in MHD turbulence may thus play a fundamental role in the dynamo process.

Shebalin, John V.↗

DOE OSTI2023

Collapsar Black Holes Are Likely Born Slowly Spinning

Collapsing stars constitute the main black hole (BH) formation channel, and are occasionally associated with the launch of relativistic jets that power γ-ray bursts (GRBs). Thus, collapsars offer an opportunity to infer the natal (before spin-up/down by accretion) BH spin directly from observations. We show that once the BH saturates with a large-scale magnetic flux, the jet power is dictated by the BH spin and mass accretion rate. Core-collapse simulations by Halevi et al. and GRB observations favor stellar density profiles that yield an accretion rate of $\dot m\approx 10^{-2} M_{\odot}s^{-1}$, weakly dependent on time. This leaves the spin as the main factor that governs the jet power. By comparing the jet power to characteristic GRB luminosities, we find that the majority of BHs associated with jets are likely born slowly spinning with a dimensionless spin of a ≃ 0.2, or a ≃ 0.5 for wobbling jets, with the main uncertainty originating in the unknown γ-ray radiative efficiency. This result could be applied to the entire core-collapse BH population, unless an anticorrelation between the stellar magnetic field and angular momentum is present. In a companion paper, Jacquemin-Ide et al., we show that regardless of the natal spin, the extraction of BH rotational energy leads to spin-down to a ≲ 0.2, consistent with gravitational-wave observations. We verify our results by performing the first 3D general-relativistic magnetohydrodynamic simulations of collapsar jets with characteristic GRB energies, powered by slowly spinning BHs. We find that jets of typical GRB power struggle to escape from the star, providing the first numerical indication that many jets fail to generate a GRB.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1992

A test suite for magnetohydrodynamical simulations

A collection is presented of MHD problems that will provide researchers who have interests in modeling MHD phenomena with a battery of tests (a 'test suite') for calibrating their numerical algorithms. Use of these tests will provide a common reference for comparison of different numerical MHD algorithms. The test suite includes both 1D and 2D problems. Taken together, these problems test the abilities of a numerical method to propagate accurately all of the MHD wave families in moving and stationary media, to capture MHD shocks and contact discontinuities, and to model accurately Lorentz force terms in multiple dimensions. An example solution of each test problem is presented. Diagnostic convergence-testing procedures that provide a quantitative evaluation of MHD algorithms are demonstrated.

Stone, James M.↗

NASA NTRS1992

Coronal shock waves

The formation of magnetohydrodynamic (MHD) shocks near the leading edge of Coronal Mass Ejections (CME's) is considered using results from both shock theory and numerical simulations. By using just the shock theory results and a simplified model corona, the following is shown: a slow shock front should be concave upward (away from the solar surface); a configuration containing intermediate shocks should also be concave upward near the center of the CME and concave downward at some distance from the CME centerline; and a fast shock configuration should be concave downward. These results are verified with numerical simulations of the MHD equations for the propagation of shocks through a quiescent coronal streamer. By examining the characteristic wave speeds in the ambient streamer, it is estimated that the various shock configurations should form for CME speeds within the following speed increments: slow shocks, 200 to 300 km/sec; intermediate shocks, 300 to 900 km/sec; and fast shocks, greater than 900 km/sec.

Steinolfson, Richard S.↗

NASA NTRS2001

Simple Analytic Expressions for the Magnetic Field of a Circular Current Loop

Analytic expressions for the magnetic induction (magnetic flux density, B) of a simple planar circular current loop have been published in Cartesian and cylindrical coordinates [1,2], and are also known implicitly in spherical coordinates [3]. In this paper, we present explicit analytic expressions for B and its spatial derivatives in Cartesian, cylindrical, and spherical coordinates for a filamentary current loop. These results were obtained with extensive use of Mathematica "TM" and are exact throughout all space outside of the conductor. The field expressions reduce to the well-known limiting cases and satisfy V · B = 0 and V x B = 0 outside the conductor. These results are general and applicable to any model using filamentary circular current loops. Solenoids of arbitrary size may be easily modeled by approximating the total magnetic induction as the sum of those for the individual loops. The inclusion of the spatial derivatives expands their utility to magnetohydrodynamics where the derivatives are required. The equations can be coded into any high-level programming language. It is necessary to numerically evaluate complete elliptic integrals of the first and second kind, but this capability is now available with most programming packages.

Simpson, James C.↗