Search NASA⌕ Search

SEARCH · Search NASA

Results for “magnetohydrodynamical simulations”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Pulsed Plasma Accelerator Modeling

This report presents the main results of the modeling task of the PPA project. The objective of this task is to make major progress towards developing a new computational tool with new capabilities for simulating cylindrically symmetric 2.5 dimensional (2.5 D) PPA's. This tool may be used for designing, optimizing, and understanding the operation of PPA s and other pulsed power devices. The foundation for this task is the 2-D, cylindrically symmetric, magnetohydrodynamic (MHD) code PCAPPS (Princeton Code for Advanced Plasma Propulsion Simulation). PCAPPS was originally developed by Sankaran (2001, 2005) to model Lithium Lorentz Force Accelerators (LLFA's), which are electrode based devices, and are typically operated in continuous magnetic field to the model, and implementing a first principles, self-consistent algorithm to couple the plasma and power circuit that drives the plasma dynamics.

Goodman, M.↗
NASA NTRS2010

Interchange Reconnection and Coronal Hole Dynamics

We investigate the effect of magnetic reconnection between open and closed fields, often referred to as interchange reconnection, on the dynamics and topology of coronal hole boundaries. The most important and most prevalent three-dimensional topology of the interchange process is that of a small-scale bipolar magnetic field interacting with a large-scale background field. We determine the evolution of such a magnetic topology by numerical solution of the fully three-dimensional MHD equations in spherical coordinates. First, we calculate the evolution of a small-scale bipole that initially is completely inside an open field region and then is driven across a coronal hole boundary by photospheric motions. Next the reverse situation is calculated in which the bipole is initially inside the closed region and driven toward the coronal hole boundary. In both cases, we find that the stress imparted by the photospheric motions results in deformation of the separatrix surface between the closed field of the bipole and the background field, leading to rapid current sheet formation and to efficient reconnection. When the bipole is inside the open field region, the reconnection is of the interchange type in that it exchanges open and closed fields. We examine, in detail, the topology of the field as the bipole moves across the coronal hole boundary and find that the field remains well connected throughout this process. Our results, therefore, provide essential support for the quasi-steady models of the open field, because in these models the open and closed flux are assumed to remain topologically distinct as the photosphere evolves. Our results also support the uniqueness hypothesis for open field regions as postulated by Antiochos et al. On the other hand, the results argue against models in which open flux is assumed to diffusively penetrate deeply inside the closed field region under a helmet streamer. We discuss the implications of this work for coronal observations.

MHD (MAGNETOHYDRODYNAMICS)↗
NASA NTRS2022

MMS Observations of Storm-Time Magnetopause Boundary Layers in the Vicinity of the Southern Cusp

During a storm-time interval around winter solstice, observations by the Magnetospheric Multi-Scale (MMS) Mission show multiple distinct magnetopause boundary layers (BLs) in the vicinity of the southern cusp. The microphysics of the solar wind-magnetosphere interaction during storm times are not well understood, because the observations are relatively lacking. This event enables the opportunity to probe the storm-time magnetopause, and observations support that MMS was near a reconnection site equatorward of the southern cusp, suggesting active reconnection in close proximity to closed magnetic flux regions in the BL. The Grid Agnostic magnetohydrodynamics (MHD) for Extended Research Applications global MHD simulation shows evidence for transient secondary reconnection sites near the southern cusp, demonstrating mechanisms to form closed field line regions of the BL.

Brandon L. Burkholder↗
NASA NTRS2023

Investigating Potential Causes for the Prediction of Spurious Magnetopause Crossings at Geosynchronous Orbit in MHD Simulations

During intense geomagnetic storms, the magnetopause can move in as far as geosynchronous orbit, leaving the satellites in that orbit out in the magnetosheath. Spacecraft operators turn to numerical models to predict the response of the magnetopause to solar wind conditions, but the predictions of the models are not always accurate. This study investigates four storms with a magnetopause crossing by at least one GOES satellite, using four magnetohydrodynamic models at NASA's Community Coordinated Modeling Center to simulate the events, and analyzes the results to investigate the reasons for errors in the predictions. Two main reasons can explain most of the erroneous predictions. First, the solar wind input to the simulations often contains features measured near the L1 point that did not eventually arrive at Earth; incorrect predictions during such periods are due to the solar wind input rather than to the models themselves. Second, while the models do well when the primary driver of magnetopause motion is a variation in the solar wind density, they tend to overpredict or underpredict the integrated Birkeland currents and their effects during times of strong negative interplanetary magnetic field (IMF) Bz, leading to poorer prediction capability. Coupling the MHD codes to a ring current model, when such a coupling is available, generally will improve the predictions but will not always entirely correct them. More work is needed to fully characterize the response of each code under strong southward IMF conditions as it relates to prediction of magnetopause location.

magnetopause↗
NASA NTRS1991

Numerical simulation of self-field MPD thrusters

A fully two-dimensional magnetohydrodynamics code has been developed to predict self-field, steady-state MPD thruster performance. The governing equations are outlined, and methods of solution are presented. Model predictions are compared with experimental data for two thruster geometries over a range of discharge currents and mass flow rates. Model limitations are evaluated, and issues concerning quasi-steady versus steady-state thruster comparisons are discussed.

La Pointe, Michael R.↗
NASA NTRS2000

Global MHD Simulations of Space Plasma Environments: Heliosphere, Comets, Magnetospheres of Plants and Satellites

Magnetohydrodynamics (MHD) provides an approximate description of a great variety of processes in space physics. Accurate numerical solutions of the MHD equations are still a challenge, but in the past decade a number of robust methods have appeared. Once these techniques made the direct solution of MHD equations feasible, a number of global three-dimensional models were designed and applied to many space physics objects. The range of these objects is truly astonishing, including active galactic nuclei, the heliosphere, the solar corona, and the solar wind interaction with planets, satellites, and comets. Outside the realm of space physics, MHD theory has been applied to such diverse problems as laboratory plasmas and electromagnetic casting of liquid metals. In this paper we present a broad spectrum of models of different phenomena in space science developed in the recent years at the University of Michigan. Although the physical systems addressed by these models are different, they all use the MHD equations as a unifying basis.

Kabin, K.↗
NASA NTRS2018

Simulated Encounters of the Parker Solar Probe with a Coronal-Hole Jet

Solar coronal jets are small, transient, collimated ejections most easily observed in coronal holes (CHs). The upcoming Parker Solar Probe (PSP) mission provides the first opportunity to encounter CH jets in situ near the Sun and examine their internal structure and dynamics. Using projected mission orbital parameters, we have simulated PSP encounters with a fully three-dimensional magnetohydrodynamic (MHD) model of a CH jet. We find that three internal jet regions, featuring different wave modes and levels of compressibility, have distinct identifying signatures detectable by PSP. The leading Alfvén wave front and its immediate wake are characterized by transAlfvénic plasma flows with mild density enhancements. This front exhibits characteristics of a fast switch-on MHD shock, whose arrival is signaled by the sudden onset of large-amplitude transverse velocity and magnetic-field oscillations highly correlated in space and time. The trailing portion is characterized by supersonic but subAlfvénicout flows of dense plasma with uncorrelated velocity and magnetic-field oscillations. This compressible region contains most of the jet's mass. The volume between the immediate wake and dense jet, the remote wake,mixes and transitions the characteristics of the two other regions. In addition to probing each region separately, we also simulate a corotational PSP-jet encounter. In this scenario, the simulated spacecraft hovers over the jet producing CH, as may occur during the mission's corotational phases, sampling each jet region in turn. We estimate that PSP will encounter numerous CH jets over the lifetime of the mission.

Roberts, Merrill A.↗
NASA NTRSDate not supplied

Performance Analysis of Magnetohydrodynamic Drag Modulation for Actively Controlled Aerocapture at Neptune

Missions to the Ice Giants are a top priority for flagship missions this coming decade. However, a fully propulsive orbital insertion into these planets requires an immense amount of fuel, taking a significant portion of the spacecraft mass and restricting the scientific payload. To mitigate this, aerocapture has been heavily investigated. Although simulations have shown that aerodynamically controlled aerocapture can successfully insert into an orbit around both Ice Giants, the deep atmospheric pass required necessitates a complex, mass expensive, and sometimes prohibitive thermal protection system. Magnetohydrodynamic drag modulation serves as a potential alternative control method for aerocapture which could not only save propellant mass compared to fully propulsive orbital insertion, but also save thermal protection system mass compared to conventional aerocapture methods. Both aerodynamically controlled and magnetohydrodynamically controlled aerocapture methodologies were simulated in NASA Langley’s high-fidelity six degree-of-freedom flight dynamics code, the Program to Optimize Simulated Trajectories II. Each method was simulated to identical missions to Neptune using a numerical predictor-corrector algorithm to optimize the control towards the target orbit. The results concluded that magnetohydrodynamic drag modulation can successfully capture and performs on par with aerodynamic drag modulation while significantly reducing the heatflux and aeroshell complexity.

Danny N Nguyen↗
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 NTRS1995

Solar wind interaction with small bodies. 1: Whistler wing signatures near to Gaspra and Ida

Data from Galileo's two asteroid flybys reveal magnetic fluctuations that we interpret as perturbations of the solar wind magnetic field caused by an interaction with the nearby asteroid. The scale sizes of the bodies (approximately 14 km for Gaspra and approximately 30 km for Ida) are intermediate between the ion and electron gyroradii, which implies that the asteroid-imposed perturbations propagate in the whistler mode. Special properties of the whistler mode include phase phase velocities that can exceed the solar wind speed and confinement of the disturbance to directions nearly aligned with the magnetic field. These features of the interaction impose a structure on the solar wind disturbance that differs greatly from the forms familiar for either magnetized or unmagnetized bodies of magnetohydrodynamic spatial scales. We examine both data and computer simulations of the interaction with special attention to what can be inferred about the interaction itself from analysis of the data.

Wang, Z.↗
NASA NTRS2007

Broken Symmetry and Coherent Structure in MHD Turbulence

Absolute equilibrium ensemble theory for ideal homogeneous magnetohydrodynamic (MHD) turbulence is fairly well developed. Theory and Simulation indicate that ideal MHD turbulence non-ergodic and contains coherent structure. The question of applicability real (i.e., dissipative) MHD turbulence is examined. Results from several very long time numerical simulations on a 64(exp 3) grid are presented. It is seen that coherent structure begins to form before decay dominates over nonlinearity. The connection with inverse spectral cascades and selective decay will also be discussed.

Shebalin, John V.↗
NASA NTRS2023

Auroral, Ionospheric and Ground Magnetic Signatures of Magnetopause Surface Modes

Surface waves on Earth's magnetopause have a controlling effect upon global magnetospheric dynamics. Since spacecraft provide sparse in situ observation points, remote sensing these modes using ground-based instruments in the polar regions is desirable. However, many open conceptual questions on the expected signatures remain. Therefore, we provide predictions of key qualitative features expected in auroral, ionospheric, and ground magnetic observations through both magnetohydrodynamic theory and a global coupled magnetosphere-ionosphere simulation of a magnetopause surface eigenmode. These show monochromatic oscillatory field-aligned currents (FACs), due to both the surface mode and its non-resonant Alfvén coupling, are present throughout the magnetosphere. The currents peak in amplitude at the equatorward edge of the magnetopause boundary layer, not the open-closed boundary as previously thought. They also exhibit slow poleward phase motion rather than being purely evanescent. We suggest the upward FAC perturbations may result in periodic auroral brightenings. In the ionosphere, convection vortices circulate the poleward moving FAC structures. Finally, surface mode signals are predicted in the ground magnetic field, with ionospheric Hall currents rotating perturbations by approximately (but not exactly) 90° compared to the magnetosphere. Thus typical dayside magnetopause surface modes should be strongest in the East-West ground magnetic field component. Overall, all ground-based signatures of the magnetopause surface mode are predicted to have the same frequency across L-shells, amplitudes that maximize near the magnetopause's equatorward edge, and larger latitudinal scales than for field line resonance. Implications in terms of ionospheric Joule heating and geomagnetically induced currents are discussed.

M. O. Archer↗
NASA NTRS1983

Reconnection rates, small scale structures and simulations

The study of reconnection in the context of one fluid, two dimensional magnetohydrodynamics (MHD), with spatially uniform constant density, viscosity and resistivity is though to retain most of the physics important in reconnection. Much of the existing reconnection literature makes use of this approach. This discussion focuses on attempts to determine the properties of reconnection solutions to MHD as precisely as possible without regard to the intrinsic limitations of the model.

Matthaeus, W. H.↗
NASA NTRS2011

Magnetohydrodynamic Turbulence and Reconnection in the Magnetotail

We have used a global MHD simulation with high spatial resolution to investigate the origin and properties of turbulence in the plasma sheet. In this simulation we imposed a steady southward IMF with a magnitude of 5 nT at the upstream simulation boundary for more than three hours followed by ninety minutes of northward IMF of the same magnitude. The solar wind number density was 20 cm-3, the thermal pressure was 20 pPa, and the velocity was 500 km/s in the x direction. The moderately high dynamic pressure confined the magnetotail to the high-grid resolution region. Even for these nominal solar wind parameters and steady driving the plasma sheet became turbulent. The power spectral densities and probability distribution functions computed from the simulations were comparable to those obtained from spacecraft observations. The largest scale vortices were associated with reconnection outflows and, in the southward IMF case, with the diversion of high speed flows in the near-Earth region. Both time and space domain analyses revealed that there were three scales present, the large scale of the driving processes, the intermediate inertial scale and the dissipative scale.

El-Alaoui, Mostafa↗
NASA NTRS2016

Coronal Magnetic Field Profiles from Shock-CME Standoff Distances

Coronagraphs observe coronal mass ejections (CMEs) and driven shocks in white light images.From these observations the shocks speed and the shocks standoff distance from the CMEs leading edge can be derived. Using these quantities, theoretical relationships between the shocks Alfvenic Mach number MA and standoff distance, and empirical radial profiles for the solar wind velocity and number density, the radial magnetic field profile upstream of the shock can be calculated. These profiles cannot be measured directly. We test the accuracy of this method for estimating the radial magnetic field profile upstream of the shock by simulating a sample CME that occurred on 29 November 2013 using the three-dimensional (3-D) magnetohydrodynamic Block-Adaptive-Tree-Solar wind-Roe-Upwind-Scheme code, retrieving shock-CME standoff distances from the simulation, and comparing the estimated and simulated radial magnetic field profiles. We find good agreement between the two profiles (within +/-30%) between 1.8 and 10R.Our simulations confirm that a linear relationship exists between the standoff distance and the inverse compression ratio at the shock. We also find very good agreement between the empirical and simulated radial profiles of the number density and speed of the solar wind and inner corona.

Schmidt, J. M.↗
NASA NTRS2006

GRMHD Simulations of Jet Formation with RAISHIN

We have developed a new three dimensional general relativistic magnetohydrodynamic (GRMHD) code, RAISHIN, using a conservative, high-resolution shock capturing scheme. Numerical fluxes are calculated using the Harten, Lax, & van Leer (HLL) approximate Riemann solver scheme. The flux-interpolated, constrained transport scheme is used to maintain a divergence-free magnetic field. We describe code performance on some test problems in both special and general relativity. Our new GRMHD code has proven to be accurate to second order and has successfully passed several numerical test problems including highly relativistic and magnetized tests in both special and general relativity. We have performed several simulations of non-rotating and rotating black hole systems with a geometrically thin accretion disk. The simulations show the formation of jets driven by the Lorentz force and the gas pressure. It appears that the rotating black hole creates an additional faster, and more collimated outflow inside a broader, slower outflow that is also generated by the rotating accretion disk around a non-rotating black hole. The kinematic jet structure could thus be a sensitive function of black hole rotation.

Mizuno, Y.↗
NASA NTRS1994

Nonlinear evolution of a large-amplitude circularly polarized Alfven wave: Low beta

The nature of turbulent cascades arising from the parametric instabilities of a monochromatic field-aligned large-amplitude circularly polarized Alfven wave is investigated via direct numerical simulation for the case of low plasma Beta and no wave dispersion. The magnetohydrodynamic code permits nonlinear couplings in the parallel direction to the ambient magnetic field and one perpendicular direction. Compressibility is included in the form of a polytropic equation of state. Anisotropic turbulent cascades, similar to those found in early incompressible two-dimensional simulations, occur after nonlinear saturation of the parallel propagating decay instability. The turbulent spectrum can be divided into three regimes: the lowest wave numbers are dominated by lower sideband remnants of the parametric process, intermediate wave numbers display nearly incompressible dynamics, and the highest wave numbers are dominated by acoustic turbulence.

Ghosh, S.↗