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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Advances and Challenges in Observations and Modeling of the Global-Sun Dynamics and Dynamo

Computational heliophysics has shed light on the fundamental physical processes inside the Sun, such as the differential rotation, meridional circulation, and dynamo-generation of magnetic fields. However, despite the substantial advances, the current results of 3D MHD simulations are still far from reproducing helioseismic inferences and surface observations. The reason is the multi-scale nature of the solar dynamics, covering a vast range of scales, which can-not be solved with the current computational resources. In such a situation, significant progress has been achieved by the mean-field approach, based on the separation of small-scale turbulence and large-scale dynamics. The mean-field simulations can reproduce solar observations, qualitatively and quantitatively, and uncover new phenomena. However, they do not reveal the complex physics of large-scale convection, solar magnetic cycles, and the magnetic self-organization that causes sunspots and solar eruptions. Thus, developing a synergy of these approaches seems to be a necessary but very challenging task.

SMD↗
NASA NTRSDate not supplied

Flare Ribbon Signatures of Reconnection Plasmoids

Solar flares are explosive space weather events that rapidly convert stored magnetic energy into bulk motion, plasma heating, and particle acceleration. Understanding the structure and dynamics of the magnetic reconnection that powers flares is critical for predicting the energy release. In particular, the amount of energy transferred to energetic particles is thought to be highly dependent on whether the reconnection is primarily turbulent (e.g., plasmoid dominated) or instead laminar. We present high-resolution MHD simulations of three-dimensional reconnection in an eruptive flare and compare the results to recent data. Although flare reconnection is challenging to observe directly in the corona, highly detailed constraints on its dynamics can be obtained from observations of flare ribbons that track the chromospheric footpoints of newly reconnected field lines. The analogues of flare ribbons in our simulations are identified by tracking discontinuous changes in field-line magnetic connectivity due to the reconnection. In our highest-resolution calculations, we find that these ribbon analogues are highly structured and exhibit many ‘whorl’ patterns that are linked to turbulent plasmoids in the reconnecting current sheet. Flare ribbon fine structure therefore reveals crucial information about the fundamental turbulent vs. laminar nature of the reconnection. We discuss the implications of these results for understanding reconnection-driven energy release throughout the heliosphere.

J. T. Dahlin↗
NASA NTRS2022

Impact of 3D Structure on Magnetic Reconnection

Results from 2.5D and 3D studies of the onset and development of the tearing instability are presented, using high fidelity resistive MHD simulations. A limited parameter study of the strength of the reconnecting field (or shear angle) was performed. An initially simple 1D equilibrium was used, consisting of a modified force-free current sheet, with periodic boundary conditions in all directions. In all cases, the linear and nonlinear evolution led to a primary current sheet between two large flux ropes. The global reconnection rate during this later stage was analyzed in all simulations. It was found that in 2.5D the primary current sheet fragmented owing to plasmoids, and as expected, the global reconnection rate, calculated using multiple methods, increases with the strength of the reconnecting field owing to a stronger Alfvén speed. In 3D, the presence of interacting oblique modes of the tearing instability complicates the simple 2.5D picture, entangling the magnetic field of the inflow and introducing a negative effect on the reconnection rate. The two competing effects of stronger Alfvén speed and entangling, which both increase with the strength of the reconnecting field, resulted in a decrease in the reconnection rate with increasing reconnecting field. For all simulations, the 3D rates were less than in 2.5D but suggest that as one goes to weak reconnecting field (or strong guide field) the system becomes more 2.5D-like and the 2.5D and 3D rates converge. These results have relevance to situations like nanoflare heating and flare current sheets in the corona.

Lars K S Daldorff↗
NASA NTRS2022

Beyond Basic Drag in Interplanetary CME Modeling: Effects of Solar Wind Pileup and High-Speed Streams

Coronal mass ejections (CMEs) cause severe space weather effects throughout our solar system. As a fast CME propagates through interplanetary space, it accumulates solar wind materials at its front. This pileup of materials, or CME-driven sheath, can be important in determining the geoeffectiveness of a CME. We take an existing arrival time model that includes expansion and deformation of the CME flux rope (ANTEATR; Kay & Gopalswamy, 2018, https://doi.org/10.1029/2018JA025780; Kay & Nieves-Chinchilla, 2021a, https://doi.org/10.1029/2020JA028911) and add a pileup procedure (PUP) as a physics-based approach to modeling the CME-driven sheath. ANTEATR-PUP solves the Rankine-Hugoniot equations for an oblique shock to determine the shock speed and sheath density, magnetic field, and temperature. The extra sheath mass affects the background drag calculation. Additionally, ANTEATR can now use any 1D profile for the background solar wind as opposed to the simple empirical models it previously relied upon. We present initial results from ANTEATR-PUP and compare with previous ANTEATR findings. Using results from an MHD simulation, we explore the effects of interactions with a static high-speed stream (HSS) on the CME's and sheath's interplanetary evolution. The drag forces essentially disappear while a CME remains within the HSS, but reappear stronger once the CME exits. The HSS-CME interaction produces the largest changes in the CME and sheath properties at 1 au when it occurs either close to the Sun near the inner simulation boundary at 0.1 au or right before the CME reaches 1 au. We estimate that these changes could significantly affect the geoeffectiveness.

C Kay↗
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 NTRSDate not supplied

Transients and Dayside Magnetospheric Interactions

Observations at and deep within the Earth’s dayside magnetosphere provide considerable evidence indicating that the solar wind-magnetosphere interaction is often unsteady. Bursty reconnection resulting in the formation of flux transfer events and transient auroral brightenings is common. Fluctuations in solar wind parameters may trigger some bursts of reconnection, while other bursts may result from intrinsic magnetopause instabilities. Large amplitude solitons often travel along the dayside magnetopause. Deeper within the dayside magnetosphere, quasi-periodic transient compressions energize radiation belt electrons and ring current ions. Most of these perturbations occur behind the quasi-parallel bow shock, indicating that they are driven by density structures generated by kinetic processes within the foreshock. This presentation compares and contrasts the amplitudes and occurrence patterns of density and dynamic pressure perturbations intrinsic to the solar wind with those generated within the foreshock, specifies how they move across the magnetosphere, examines recent global MHD simulation results for their interaction with the bow shock, magnetosheath, and magnetopause. It then outlines new paths toward quantifying the significance of bursty reconnection to the overall interaction, via new imaging and analysis techniques.

D. G. Sibeck↗
NASA NTRS2024

Direct Imaging of Magnetohydrodynamic Wave Mode Conversion Near A 3D Null Point on the Sun

Mutual conversion of various kinds of magnetohydrodynamic (MHD) waves can have profound impacts on wave propagation, energy transfer, and heating of the solar chromosphere and corona. Mode conversion occurs when an MHD wave travels through a region where the Alfvén and sound speeds are equal (e.g., a 3D magnetic null point). Here we report the direct extreme ultraviolet (EUV) imaging of mode conversion from a fast-mode to a slow-mode MHD wave near a 3D null point using Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA) observations. An incident fast EUV wavefront associated with an adjacent eruptive flare propagates laterally through a neighboring pseudostreamer. Shortly after the passage of the fast EUV wave through the null point, a slow-mode wave appears near the null that propagates upward along the open structures and simultaneously downward along the separatrix encompassing the fan loops of the pseudostreamer base. These observations suggest the existence of mode conversion near 3D nulls in the solar corona, as predicted by theory and MHD simulations. Moreover, we observe decaying transverse oscillations in both the open and closed structures of the pseudostreamer, along with quasiperiodic type III radio bursts indicative of repetitive episodes of electron acceleration.

Pankaj Kumar↗
NASA NTRS2024

The Role of Magnetic Reconnection in the Structure and Dynamics of Fast Coronal Mass Ejections

Both observations and models of flare-associated coronal mass ejections (CMEs) suggest that magnetic reconnection in an ejection's wake substantially increases the net, outward Lorentz force acceler­ating the CME. A stronger outward force can cause a feedback loop, driving further magnetic reconnection in a "reconnective instability." The flux accretion model captures this by relating reconnected flux, ∆ rec , and magnetic field strength, B CME , to increased outward Lorentz force, ∆F r . To better understand recon­nection's role in CME dynamics, we analyze two snapshots from a 2.5D, MHD simulation of a breakout eruption. Outward Lorentz forces increase substantially as reconnection proceeds, caused primarily by "flank currents," which flow just inside the boundary of the rising ejection's wake and parallel to its axis. This model's reconnection jet also alters the ejection's internal structure, an effect that could be sought in observations. Analyzing reconnection-induced Lorentz forces in 3D simulations could provide additional insights into CME dynamics.

Magnetic Reconnection↗
NASA NTRS1995

Time-dependent simulation of oblique MHD cosmic-ray shocks using the two-fluid model

Using a new, second-order accurate numerical method we present dynamical simulations of oblique MHD cosmic-ray (CR)-modified plane shock evolution. Most of the calculations are done with a two-fluid model for diffusive shock acceleration, but we provide also comparisons between a typical shock computed that way against calculations carried out using the more complete, momentum-dependent, diffusion-advection equation. We also illustrate a test showing that these simulations evolve to dynamical equilibria consistent with previously published steady state analytic calculations for such shocks. In order to improve understanding of the dynamical role of magnetic fields in shocks modified by CR pressure we have explored for time asymptotic states the parameter space of upstream fast mode Mach number, M(sub f), and plasma beta. We compile the results into maps of dynamical steady state CR acceleration efficiency, epsilon(sub c). We have run simulations using constant, and nonisotropic, obliquity (and hence spatially) dependent forms of the diffusion coefficient kappa. Comparison of the results shows that while the final steady states achieved are the same in each case, the history of CR-MHD shocks can be strongly modified by variations in kappa and, therefore, in the acceleration timescale. Also, the coupling of CR and MHD in low beta, oblique shocks substantially influences the transient density spike that forms in strongly CR-modified shocks. We find that inside the density spike a MHD slow mode wave can be generated that eventually steepens into a shock. A strong layer develops within the density spike, driven by MHD stresses. We conjecture that currents in the shear layer could, in nonplanar flows, results in enhanced particle accretion through drift acceleration.

Frank, Adam↗
NASA NTRS2003

Numerical Simulation of Turbulent MHD Flows Using an Iterative PNS Algorithm

A new parabolized Navier-Stokes (PNS) algorithm has been developed to efficiently compute magnetohydrodynamic (MHD) flows in the low magnetic Reynolds number regime. In this regime, the electrical conductivity is low and the induced magnetic field is negligible compared to the applied magnetic field. The MHD effects are modeled by introducing source terms into the PNS equation which can then be solved in a very efficient manner. To account for upstream (elliptic) effects, the flowfields are computed using multiple streamwise sweeps with an iterated PNS algorithm. Turbulence has been included by modifying the Baldwin-Lomax turbulence model to account for MHD effects. The new algorithm has been used to compute both laminar and turbulent, supersonic, MHD flows over flat plates and supersonic viscous flows in a rectangular MHD accelerator. The present results are in excellent agreement with previous complete Navier-Stokes calculations.

Kato, Hiromasa↗
NASA NTRS1976

MHD solution of interplanetary disturbances generated by simulated velocity perturbations

An MHD time-dependent numerical simulation, restricted to the solar equatorial plane, is used to demonstrate the interplanetary disturbances caused by several simplified coronal holes. Each 'hole' is assumed to have a configuration such that the higher solar wind velocity produced within their longitudinal extent is Gaussian over a 7-day period at the inner boundary (0.3 AU) of the numerical simulation. A second, twin coronal hole is assumed to rotate on the solar disk behind its predecessor. It is shown that the first coronal hole-produced interplanetary shock ensemble is overtaken by the second ensemble because of the higher velocity, lower density environment into which the latter propagates. A number of features predicted by MHD similarity theory are confirmed by the numerical simulation. These features include (1) strong azimuthal magnetic and plasma density compression, accompanied with average temperature depression, at the contact surface between forward and reverse shock ensembles, and (2) increasing spatial separation distance between forward and reverse shocks.

Dryer, M.↗
NASA NTRS2016

Extended Magnetohydrodynamics with Embedded Particle-in-Cell Simulation of Ganymede's Magnetosphere

We have recently developed a new modeling capability to embed the implicit particle-in-cell (PIC) model iPIC3D into the Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme magnetohydrodynamic (MHD) model. The MHD with embedded PIC domains (MHO-EPIC) algorithm Is a two-way coupled kinetic-fluid model. As one of the very first applications of the MHD-EPIC algorithm, we simulate the Interaction between Jupiter's magnetospherlc plasma and Ganymede's magnetosphere. We compare the MHO-EPIC simulations with pure Hall MHD simulations and compare both model results with Galileo observations to assess the Importance of kinetic effects In controlling the configuration and dynamics of Ganymede's magnetosphere. We find that the Hall MHD and MHO-EPIC solutions are qualitatively similar, but there are significant quantitative differences. In particular. the density and pressure inside the magnetosphere show different distributions. For our baseline grid resolution the PIC solution is more dynamic than the Hall MHD simulation and it compares significantly better with the Galileo magnetic measurements than the Hall MHD solution. The power spectra of the observed and simulated magnetic field fluctuations agree extremely well for the MHD-EPIC model. The MHO-EPIC simulation also produced a few flux transfer events (FTEs) that have magnetic signatures very similar to an observed event. The simulation shows that the FTEs often exhibit complex 3-0 structures with their orientations changing substantially between the equatorial plane and the Galileo trajectory, which explains the magnetic signatures observed during the magnetopause crossings. The computational cost of the MHO-EPIC simulation was only about 4 times more than that of the Hall MHD simulation.

Toth, Gabor↗
NASA NTRS2012

The Effect of Magnetic Topology on the Escape of Flare Particles

Magnetic reconnection in the solar atmosphere is believed to be the driver of most solar explosive phenomena. Therefore, the topology of the coronal magnetic field is central to understanding the solar drivers of space weather. Of particular importance to space weather are the impulsive Solar Energetic particles that are associated with some CME/eruptive flare events. Observationally, the magnetic configuration of active regions where solar eruptions originate appears to agree with the standard eruptive flare model. According to this model, particles accelerated at the flare reconnection site should remain trapped in the corona and the ejected plasmoid. However, flare-accelerated particles frequently reach the Earth long before the CME does. We present a model that may account for the injection of energetic particles onto open magnetic flux tubes connecting to the Earth. Our model is based on the well-known 2.5D breakout topology, which has a coronal null point (null line) and a four-flux system. A key new addition, however, is that we include an isothermal solar wind with open-flux regions. Depending on the location of the open flux with respect to the null point, we find that the flare reconnection can consist of two distinct phases. At first, the flare reconnection involves only closed field, but if the eruption occurs close to the open field, we find a second phase involving interchange reconnection between open and closed. We argue that this second reconnection episode is responsible for the injection of flare-accelerated particles into the interplanetary medium. We will report on our recent work toward understanding how flare particles escape to the heliosphere. This work uses high-resolution 2.5D MHD numerical simulations performed with the Adaptively Refined MHD Solver (ARMS).

Antiochos, S. K.↗
NASA NTRS2012

A Model for SEP Escape

Magnetic reconnection in the solar atmosphere is believed to be the driver of most solar explosive phenomena. Therefore, the topology of the coronal magnetic field is central to understanding the solar drivers of space weather. Of particular importance to space weather are the impulsive Solar Energetic particles that are associated with some CME/eruptive flare events. Observationally, the magnetic configuration of active regions where solar eruptions originate appears to agree with the standard eruptive flare model. According to this model, particles accelerated at the flare reconnection site should remain trapped in the corona and the ejected plasmoid. However, flare-accelerated particles frequently reach the Earth long before the CME does. We present a model that may account for the injection of energetic particles onto open magnetic flux tubes connecting to the Earth. Our model is based on the well-known 2.5D breakout topology, which has a coronal null point (null line) and a four-flux system. A key new addition, however, is that we include an isothermal solar wind with open-flux regions. Depending on the location of the open flux with respect to the null point, we find that the flare reconnection can consist of two distinct phases. At first, the flare reconnection involves only closed field, but if the eruption occurs close to the open field, we find a second phase involving interchange reconnection between open and closed. We argue that this second reconnection episode is responsible for the injection of flare-accelerated particles into the interplanetary medium. We will report on our recent work toward understanding how flare particles escape to the heliosphere. This work uses high-resolution 2.5D MHD numerical simulations performed with the Adaptively Refined MHD Solver (ARMS).

Antiochos, S. K.↗
NASA NTRS1993

Structure of the dayside reconnection layer in resistive MHD and hybrid models

Numerical simulations were performed to investigate the structure of the reconnection layer at the dayside magnetopause. Two typical cases are examined in detail; both are asymmetric in magnetic field and plasma density. In case 1, the guide fields in the magnetosheath and in the magnetosphere are set at zero and thus the tangential magnetic fields on the two sides of the initial current sheet are exactly antiparallel. In case 2, the angle between the tangential magnetic fields on the two sides of the initial current sheet is 145 deg. The results obtained from a resistive MHD model and from a hybrid model are found to be different. In the MHD simulation of case 1, a 2-4 intermediate shock is found to bound the reconnection layer on the magnetosheath side, while an Alfven wave pulse bounds the reconnection layer on the magnetospheric side. In case 2, it is found that a time-dependent intermediate shock (TDIS) bounds the reconnection layer on the magnetosheath side, with a slow expansion wave propagating behind. With the MHD simulations, in the general case in which the tangential magnetic fields on the two sides of the initial current sheet are not exactly antiparallel, a rotational discontinuity across which the tangential magnetic field rotates, a large angle is found to bound the reconnection layer on the magnetosheath side.

Lin, Y.↗
NASA NTRS1995

Analytic MHD Theory for Earth's Bow Shock at Low Mach Numbers

A previous MHD theory for the density jump at the Earth's bow shock, which assumed the Alfven M(A) and sonic M(s) Mach numbers are both much greater than 1, is reanalyzed and generalized. It is shown that the MHD jump equation can be analytically solved much more directly using perturbation theory, with the ordering determined by M(A) and M(s), and that the first-order perturbation solution is identical to the solution found in the earlier theory. The second-order perturbation solution is calculated, whereas the earlier approach cannot be used to obtain it. The second-order terms generally are important over most of the range of M(A) and M(s) in the solar wind when the angle theta between the normal to the bow shock and magnetic field is not close to 0 deg or 180 deg (the solutions are symmetric about 90 deg). This new perturbation solution is generally accurate under most solar wind conditions at 1 AU, with the exception of low Mach numbers when theta is close to 90 deg. In this exceptional case the new solution does not improve on the first-order solutions obtained earlier, and the predicted density ratio can vary by 10-20% from the exact numerical MHD solutions. For theta approx. = 90 deg another perturbation solution is derived that predicts the density ratio much more accurately. This second solution is typically accurate for quasi-perpendicular conditions. Taken together, these two analytical solutions are generally accurate for the Earth's bow shock, except in the rare circumstance that M(A) is less than or = 2. MHD and gasdynamic simulations have produced empirical models in which the shock's standoff distance a(s) is linearly related to the density jump ratio X at the subsolar point. Using an empirical relationship between a(s) and X obtained from MHD simulations, a(s) values predicted using the MHD solutions for X are compared with the predictions of phenomenological models commonly used for modeling observational data, and with the predictions of a modified phenomenological model proposed recently. The similarities and differences between these results are illustrated using plots of X and a(s) predicted for the Earth's bow shock. The plots show that the new analytic solutions agree very well with the exact numerical MHD solutions and that these MHD solutions should replace the corresponding phenomenological relations in comparisons with data. Furthermore, significant differences exist between the standoff distances predicted at low M(A) using the MHD models versus those predicted by the new modified phenomenological model. These differences should be amenable to observational testing.

Grabbe, Crockett L.↗