Search NASA⌕ Search

SEARCH · Search NASA

Results for “MHD simulation”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Velocity shear generation of solar wind turbulence

A two-dimensional incompressible MHD spectral code is used to show that shear-driven turbulence is a possible means for producing many observed properties of the evolution of the magnetic and velocity fluctuations in the solar wind and, in particular, the evolution of the cross helicity ('Alfvenicity') at small scales. It is shown that large-scale shear can nonlinearly produce a cascade to smaller scale fluctuations even when the linear Kelvin-Helmholtz mode is stable, and that a roughly power law inertial range is established by this process. The evolution found is similar to that seen in some other simulations of MHD turbulence.

Roberts, D. A.↗
NASA NTRS2009

2D and 3D Numerical Simulations of Flux Cancellation

Cancellation of magnetic flux in the solar photosphere and chromosphere has been linked observationally and theoretically to a broad range of solar activity, from filament channel formation to CME initiation. Because this phenomenon is typically measured at only a single layer in the atmosphere, in the radial (line of sight) component of the magnetic field, the actual processes behind this observational signature are ambiguous. It is clear that reconnection is involved in some way, but the location of the reconnection sites and associated connectivity changes remain uncertain in most cases. We are using numerical modeling to demystify flux cancellation, beginning with the simplest possible configuration: a subphotospheric Lundquist flux tube surrounded by a potential field, immersed in a gravitationally stratified atmosphere, spanning many orders of magnitude in plasma beta. In this system, cancellation is driven slowly by a 2-cell circulation pattern imposed in the convection zone, such that the tops of the cells are located around the beta=1 level (i.e., the photosphere) and the flows converge and form a downdraft at the polarity inversion line; note however that no flow is imposed along the neutral line. We will present the results of 2D and 3D MHD-AMR simulations of flux cancellation, in which the flux at the photosphere begins in either an unsheared or sheared state. In all cases, a low-lying flux rope is formed by reconnection at the polarity inversion line within a few thousand seconds. The flux rope remains stable and does not rise, however, in contrast to models which do not include the presence of significant mass loading.

Karpen, Judith T.↗
NASA NTRS2009

20 and 3D Numerical Simulations of Flux Cancellation

Cancellation of magnetic flux in the solar photosphere and chromosphere has been linked observationally and theoretically to a broad range of solar activity, from filament channel formation to CME initiation. Because this phenomenon is typically measured at only a single layer in the atmosphere, in the radial (line of sight) component of the magnetic field, the actual processes behind this observational signature are ambiguous. It is clear that reconnection is involved in some way, but the location of the reconnection sites and associated connectivity changes remain uncertain in most cases. We are using numerical modeling to demystify flux cancellation, beginning with the simplest possible configuration: a subphotospheric Lundquist flux tube surrounded by a potential field, immersed in a gravitationally stratified atmosphere, spanning many orders of magnitude in plasma beta. In this system, cancellation is driven slowly by a 2-cell circulation pattern imposed in the convection zone, such that the tops of the cells are located around the beta= 1 level (Le., the photosphere) and the flows converge and form a downdraft at the polarity inversion line; note however that no flow is imposed along the neutral line. We will present the results of 2D and 3D MHD-AMR simulations of flux cancellation, in which the flux at the photosphere begins in either an unsheared or sheared state. In all cases, a lOW-lying flux rope is formed by reconnection at the polarity inversion line within a few thousand seconds. The flux rope remains stable and does not rise, however, in contrast to models which do not include the presence of significant mass loading.

Karpen, Judith T.↗
NASA NTRS2014

Simulations of Emerging Magnetic Flux. II. The Formation of Unstable Coronal Flux Ropes and the Initiation of Coronal Mass Ejections

We present results from three-dimensional magnetohydrodynamic simulations of the emergence of a twisted convection zone flux tube into a pre-existing coronal dipole field. As in previous simulations, following the partial emergence of the sub-surface flux into the corona, a combination of vortical motions and internal magnetic reconnection forms a coronal flux rope. Then, in the simulations presented here, external reconnection between the emerging field and the pre-existing dipole coronal field allows further expansion of the coronal flux rope into the corona. After sufficient expansion, internal reconnection occurs beneath the coronal flux rope axis, and the flux rope erupts up to the top boundary of the simulation domain (approximately 36 Mm above the surface).We find that the presence of a pre-existing field, orientated in a direction to facilitate reconnection with the emerging field, is vital to the fast rise of the coronal flux rope. The simulations shown in this paper are able to self-consistently create many of the surface and coronal signatures used by coronal mass ejection (CME) models. These signatures include surface shearing and rotational motions, quadrupolar geometry above the surface, central sheared arcades reconnecting with oppositely orientated overlying dipole fields, the formation of coronal flux ropes underlying potential coronal field, and internal reconnection which resembles the classical flare reconnection scenario. This suggests that proposed mechanisms for the initiation of a CME, such as "magnetic breakout," are operating during the emergence of new active regions.

corona↗
NASA NTRS1991

Oblique MHD cosmic-ray modified shocks: Two-fluid numerical simulations

We present the first results of time dependent, two-fluid, cosmic-ray (CR) modified, MHD shock simulations. The calculations were carried out with a new numerical code for 1-D ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they appear to have converged to dynamical steady states. We also present results from an initial time dependent simulation which extend the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescale than parallel shocks.

Frank, Adam↗
NASA NTRS1985

On the formation of coronal cavities

The formation of a coronal cavity and its relation to a quiescent prominence is studied theoretically. The stability of condensation modes of a plasma in the coronal streamer model (Steinolfson et al., 1982) is considered using a two-dimensional time-dependent ideal MHD numerical simulation. It is found that a plasma with beta = 0.5 is unstable but one with beta = 4 is stable because the density enhancement of the plasma trapped by the closed fields increases with the strength of the magnetic field. The means by which condensation modes can produce a coronal cavity and/or initiate the formation of a prominence (depending on the field configuration) are discussed. It is argued that prominence and cavity material is all supplied from the chromospheric level in the form of spicules.

An, C.-H.↗
NASA NTRS1992

Constraints on cosmic ray propagation in the galaxy

The goal was to derive a more detailed picture of magnetohydrodynamic turbulence in the interstellar medium and its effects on cosmic ray propagation. To do so, radio astronomical observations (scattering and Faraday rotation) were combined with knowledge of solar system spacecraft observations of MHD turbulence, simulations of wave propagation, and modeling of the galactic distribution to improve the knowledge. A more sophisticated model was developed for the galactic distribution of electron density turbulence. Faraday rotation measure data was analyzed to constrain magnetic field fluctuations in the ISM. VLBI observations were acquired of compact sources behind the supernova remnant CTA1. Simple calculations were made about the energies of the turbulence assuming a direct link between electron density and magnetic field variations. A simulation is outlined of cosmic ray propagation through the galaxy using the above results.

Cordes, James M.↗
NASA NTRS1993

Magnetospheric response to solar wind dynamic pressure variations: Interaction of interplanetary tangential discontinuities with the bow shock

Some magnetic impulse events observed in the polar region are related to vortices associated with plasma convection in the ionosphere. Recent analyses of satellite and ground data suggest that the interaction of solar wind dynamic pressure pulses and the magnetosphere may lead to the formation of velocity vortices in the magnetopause boundary layer region. This can in turn lead to the presence of vortices in the polar ionosphere. However, before reaching the Earth's magnetopause, these interplanetary pressure pulses must interact with and pass through the bow shock. A variation of the solar wind dynamic pressure may be associated with shocks, magnetic holes, or tangential discontinuities (TDs) in the interplanetary medium. We study the interaction of interplanetary TDs with the Earth's bow shock (BS) using both theoretical analysis and MHD computer simulations. It is found that as a result of the collision between a TD and the BS, the jump in the solar wind dynamic pressure associated with the TD is significantly modified, the bow shock moves, and a new fast shock or fast rarefaction wave, which propagates in the downstream direction, is excited. Our theoretical analysis shows that the change in the plasma density across the interplanetary TD plays the most important role in the collision process. In the case with an enhanced dynamic pressure behind the interplanetary TD, the bow shock is intensified in strength and moves in the earthward direction. The dynamic pressure jump associated with the transmitted TD is generally reduced from the value before the interaction. A fast compressional shock is excited ahead of the transmitted TD and propagates toward the Earth's magnetosphere. For the case in which the dynamic pressure is reduced behind the interplanetary TD, the pressure jump across the transmitted TD is substantially weakened, the bow shock moves in the sunward direction, and a rarefaction wave which propagates downstream is excited. We also simulate and discuss the interaction of a pair of tangential discontinuities, which may correspond to a magnetic hole, with the BS.

Wu, Bor-Han↗
NASA NTRS1995

Origin of coronal mass ejection and magnetic cloud: Thermal or magnetic driven?

A fundamental problem in Solar-Terrestrial Physics is the origin of the solar transient plasma output, which includes the coronal mass ejection and its interplanetary manifestation, e.g. the magnetic cloud. The traditional blast wave model resulted from solar thermal pressure impulse has faced with challenge during recent years. In the MHD numerical simulation study of CME, the authors find that the basic feature of the asymmetrical event on 18 August 1980 can be reproduced neither by a thermal pressure nor by a speed increment. Also, the thermal pressure model fails in simulating the interplanetary structure with low thermal pressure and strong magnetic field strength, representative of a typical magnetic cloud. Instead, the numerical simulation results are in favor of the magnetic field expansion as the likely mechanism for both the asymmetrical CME event and magnetic cloud.

Zhang, Gong-Liang↗
NASA NTRS1998

Chinks in Solar Dynamo Theory: Turbulent Diffusion, Dynamo Waves and Magnetic Helicity

In this first year of our investigation we explored the role of compressibility and stratification in the dissipation of magnetic fields. The predictions of Mean Field Electrodynamics have been questioned because of the strong feedback of small scale magnetic structure on the velocity fields. In 2-D, this nonlinear feedback results in a lengthening of the turbulent decay time. In 3-D alpha-quenching is predicted. Previous studies assumed a homogeneous fluid. This first year we present recent results from 2-D compressible MHD decay simulations in a highly stratified atmosphere that more closely resembles to solar convection zone. We have applied for NCCS T3E time to assist in the performance of our 3-D calculations.

DeLuca, E. E.↗
NASA NTRS2008

Simulations of Solar Wind Turbulence

Recently we have restructured our approach to simulating magnetohydrodynamic (MHD) turbulence in the solar wind. Previously, we had defined a 'virtual' heliosphere that contained, for example, a tilted rotating current sheet, microstreams, quasi-two-dimensional fluctuations as well as Alfven waves. In this new version of the code, we use the global, time-stationary, WKB Alfven wave-driven solar wind model developed by Usmanov and described in Usmanov and Goldstein [2003] to define the initial state of the system. Consequently, current sheets, and fast and slow streams are computed self-consistently from an inner, photospheric, boundary. To this steady-state configuration, we add fluctuations close to, but above, the surface where the flow become super-Alfvenic. The time-dependent MHD equations are then solved using a semi-discrete third-order Central Weighted Essentially Non-Oscillatory (CWENO) numerical scheme. The computational domain now includes the entire sphere; the geometrical singularity at the poles is removed using the multiple grid approach described in Usmanov [1996]. Wave packets are introduced at the inner boundary such as to satisfy Faraday's Law [Yeh and Dryer, 1985] and their nonlinear evolution are followed in time.

Goldstein, Melvyn L.↗
NASA NTRS2017

Building Reliable Forecasts of Solar Activity

Solar ionizing radiation critically depends on the level of the Sun’s magnetic activity. For robust physics-based forecasts, we employ the procedure of data assimilation, which combines theoretical modeling and observational data such that uncertainties in both the model and the observations are taken into account. Currently we are working in two major directions: 1) development of a new long-term forecast procedure on time-scales of the 11-year solar cycle, using a 2-dimensional mean-field dynamo model and synoptic magnetograms; 2) development of 3-dimensional radiative MHD (Magnetohydrodynamic) simulations to investigate the origin and precursors of local manifestations of magnetic activity, such as the formation of magnetic structures and eruptive dynamics.

Kitiashvili, Irina↗
NASA NTRS2017

Realistic 3D Radiative Modeling of Turbulent Structure of Moderate-Mass Stars and Sun

Understanding the turbulent dynamics of the Sun and stars is a critical element for interpreting observed processes and phenomena on different scales and for predicting extreme events such as flares and super flares. High-­‐resolution observations of the Sun and high-­‐fidelity radiative MHD numerical simulations have substantially advanced our understanding of solar and stellar local dynamics and magnetism from the upper convection zone to the atmosphere and corona. However, global modeling of the Sun with such a high degree of realism is currently not affordable due to the extremely high computational cost of resolving the scales in the convection zone. The physics of the deep solar dynamics can be effectively addressed through modeling more massive solar-­‐type stars where the convection zone is shallower and the convective overturning time is much shorter than those on the Sun. We present recent 3D realistic simulation results of moderate-­‐mass stars and discuss links between solar and stellar dynamics, such as the multiscale structure of granulation, convective overshoot, and others. In particular, these simulations have provided better understanding of the dynamics of the tachocline (the overshoot layer at the bottom of the convection zone) and have explained long-­‐ standing results from helioseismology.

Realistic↗
NASA NTRS2022

STITCH: A Subgrid-Scale Model for Energy Buildup in the Solar Corona

The solar corona routinely exhibits explosive activity, in particular coronal mass ejections and their accompanying eruptive ares, that have global-scale consequences. These events and their smaller counterparts, coronal jets, originate in narrow, sinuous lament channels. The key processes that form and evolve the channels operate on still smaller spatial scales and much longer time scales, culminating in a vast separation of characteristic lengths and times that govern these explosive phenomena. In this article, we describe implementation and tests of an efficient subgrid-scale model for generating eruptive structures in magnetohydrodynamics (MHD) coronal simulations. STITCH {STatistical InjecTion of Condensed Helicity { is a physics-based, reduced representation of helicity condensation: a process wherein small-scale vortical surface convection forms ubiquitous current sheets, and pervasive reconnection across the sheets mediates an inverse cascade of magnetic helicity and free energy, thereby forming the lament channels. We have developed a formalism, STITCH, that abstracts these complex processes into a single term in Ohm's law and the induction equation that directly injects tangential magnetic flux into the low corona. We show that our approach is in very good agreement with a full helicity-condensation calculation that treats all of the dynamics explicitly, while enabling substantial reductions in temporal duration and spatial resolution. In addition, we illustrate the flexibility of STITCH at forming localized lament channels and at energizing complex surface flux distributions that have sinuous boundaries. STITCH is simple to implement and computationally ecient, making it a powerful technique for physics-based modeling of solar eruptive events.

J T Dahlin↗
NASA NTRSDate not supplied

Variety in the Variability of Accreting Supermassive Binary Black Holes

Accreting supermassive binary black holes are key multi-messenger sources for LISA, yet are challenging to simulate realistically as solving the radiation magnetohydrodynamics (MHD) equations over the full dynamic spatio-temporal range of the problem is computationally infeasible at present. We report on our collaboration's progress simulating these systems using general relativistic MHD simulations and dynamic GR. In order to cover a larger temporal range in one set of simulations, we constrain our view to the circumbinary disk region and measure the effects binary mass ratio, accretion disk size, and black hole spin have on the structure and variability of the accretion flow. We particularly emphasize how these parameters influence the overdensity feature, which orbits the binary near the edge of the cavity, since it is responsible for most of the electromagnetic emission's variability and variability is a key signature of a system being a binary. Extending to smaller length scales, we will report on simulations following accretion all the down to the event horizons so that we can begin to investigate how black hole spin affects mini-disk dynamics, accretion rate, and jet power.

Scott C Noble↗