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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2024

Exploring Standing and Reflected Slow-mode Waves in Flaring Coronal Loops: A Parametric Study Using 2.5D MHD Modeling

Abstract Recent observations of reflected propagating and standing slow-mode waves in hot flaring coronal loops have spurred our investigation into their underlying excitation and damping mechanisms. To understand these processes, we conduct 2.5D magnetohydrodynamic (MHD) simulations using an arcade active region model that includes a hot and dense loop. Our simulations allow for in-depth parametric investigations complementing and expanding our previous 3D MHD modeling results. We excite these waves using a large-amplitude, flow pulse applied at one footpoint of the loop in two distinct models as motivated by observations from the Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA). The first model (Model 1) incorporates classical compressive viscosity coefficient, while the second model (Model 2) adopts a 10-times enhanced viscosity coefficient. We obtain the following major results: (1) Our 2.5D MHD simulations reinforce previous conclusions derived from 1D and 3D MHD models that significantly enhanced viscosity is crucial for the rapid excitation of standing slow waves with damping times consistent with observations by Wang et al. (2015). (2)We uncover that nonlinearity in Model 1 delays the conversion of a reflected propagating wave into a standing wave. In contrast, Model 2 exhibits a much weak influence of nonlinearity on the excitation time of standing waves, thanks to the suppression of these effects by enhanced viscosity. (3) Our results reveal that the transverse temperature structure holds more influence on wave behavior than the density structure. In Model 1, increased loop temperature contrast significantly enhances wave trapping within the structure, mitigating the impact of temperature-dependent viscous damping. Conversely, in Model 2, the impact of temperature structure on wave behavior weakens in comparison to the effect of viscosity. (4) Model 1 displays evident nonlinear coupling to the fast and kink magnetoacoustic waves and pronounced wave leakage into the corona. Model 2 exhibits significantly weaker effects in this regard. Analyzing three observed wave events by SDO/AIA aligns with Model 2 predictions, providing further support for the substantial viscosity increase. Our 2.5D study unravels the complex interplay of wave-flow phenomena and nonlinear processes in coronal loops, extending our previous 1D modeling results to incorporate more realistic loop geometry. This provides insights into scenarios where 3D effects may be neglected, thereby enhancing our understanding of the intricate dynamics of the solar corona.

Tongjiang Wang↗
NASA NTRS1986

The synchrotron-maser theory of type II solar radio emission processes - The physical model and generation mechanism

A theory is proposed to explain the generation mechanism of type II solar radio bursts. It is suggested that the shock wave formed at the leading edge of a coronal transient can accelerate electrons. Because of the nature of the acceleration process, the energized electrons can possess a 'hollow-beam' type distribution function. When the electron beam propagates along the ambient magnetic field to lower altitudes and attains larger pitch angles, a synchrotron-maser instability can set in. This instability leads to the amplification of unpolarized or weakly polarized radiation. The present discussion incorporates a model which describes the ambient magnetic field and background plasma by means of MHD simulation. The potential emission regions may be located approximately, according to the time-dependent MHD simulation. Since the average local plasma frequency in the source region can be evaluated from the MHD model, the frequent drift associated with the radiation may be estimated. The result seems to be in good agreement with that derived from observations.

Wu, C. S.↗
NASA NTRS2016

Disk Emission from Magnetohydrodynamic Simulations of Spinning Black Holes

We present the results of a new series of global, three-dimensional, relativistic magnetohydrodynamic (MHD) simulations of thin accretion disks around spinning black holes. The disks have aspect ratios of H/R approx. 0.05 and spin parameters of a/M = 0, 0.5, 0.9, and 0.99. Using the ray-tracing code Pandurata, we generate broadband thermal spectra and polarization signatures from the MHD simulations. We find that the simulated spectra can be well fit with a simple, universal emissivity profile that better reproduces the behavior of the emission from the inner disk, compared to traditional analyses carried out using a Novikov-Thorne thin disk model. Finally, we show how spectropolarization observations can be used to convincingly break the spin-inclination degeneracy well known to the continuum-fitting method of measuring black hole spin.

accretion↗
NASA NTRS2022

Electron Anisotropies in Magnetotail Dipolarization Events

Anisotropic electron distributions can be a major source of free energy for the generation of microinstabilities and waves. Here we investigate specifically anisotropies of the suprathermal electrons associated with dipolarization events in the geomagnetic tail. The investigation is based on an MHD simulation of magnetotail reconnection, flow bursts and dipolarization. This simulation is used to trace test particles in the electromagnetic fields of the MHD simulation and investigate their acceleration and contributions to flux increases. The simulated velocity distributions yield anisotropies which can be dominantly parallel (“cigar”-shaped), perpendicular (“pancake”-shaped) to the magnetic field, or a combination thereof (e.g., “rolling-pin”-shaped), depending on location, relative timing, and energy.

electron anisotropies↗
NASA NTRS1997

Numerical simulation of MHD turbulence in three dimensions

The evolution of Alfvenic turbulence in 3D spherical geometry can now be studied. In simulations, a fast stream is sandwiched between two slower streams. The inflow is both supersonic and superAlfvenic. Alfven waves entering the box are convected into the medium and interact nonlinearly with the velocity shear and with any structures (i.e., flux tubes) that might be present. These initial simulations suggest that velocity shear, even in spherical geometry, is able to drive a turbulent cascade which results in approximately Kolmogoroff-like power spectra.

Goldstein, M. L.↗
NASA NTRS2013

Characterization of the Optical and X-ray Properties of the Northwestern Wisps in the Crab Nebula

We have studied the wisps to the northwest of the Crab pulsar as part of a multi-wavelength campaign in the visible and in X-rays. Optical observations were obtained using the Nordic Optical Telescope in La Palma and X-ray observations were made with the Chandra X-ray Observatory. The observing campaign took place from October 2010 until September 2012. About once per year we observe wisps forming and peeling off from (or near) the region commonly associated with the termination shock of the pulsar wind. We find that the exact locations of the northwestern wisps in the optical and in X-rays are similar but not coincident, with X-ray wisps preferentially located closer to the pulsar. This suggests that the optical and X-ray wisps are not produced by the same particle distribution. It is also interesting to note that the optical and radio wisps are also separated from each other (Bietenholz et al. 2004). Our measurements and their implications are interpreted in terms of a Doppler-boosted ring model that has its origin in MHD modeling. While the Doppler boosting factors inferred from the X-ray wisps are consistent with current MHD simulations of PWNe, the optical boosting factors are not, and typically exceed values from MHD simulations by about a factor of 4.

Weisskopf, M. C.↗
NASA NTRS1991

Large-scale gas dynamical processes affecting the origin and evolution of gaseous galactic halos

Observations of galactic halo gas are consistent with an interpretation in terms of the galactic fountain model in which supernova heated gas in the galactic disk escapes into the halo, radiatively cools and forms clouds which fall back to the disk. The results of a new study of several large-scale gas dynamical effects which are expected to occur in such a model for the origin and evolution of galactic halo gas will be summarized, including the following: (1) nonequilibrium absorption line and emission spectrum diagnostics for radiatively cooling halo gas in our own galaxy, as well the implications of such absorption line diagnostics for the origin of quasar absorption lines in galactic halo clouds of high redshift galaxies; (2) numerical MHD simulations and analytical analysis of large-scale explosions ad superbubbles in the galactic disk and halo; (3) numerical MHD simulations of halo cloud formation by thermal instability, with and without magnetic field; and (4) the effect of the galactic fountain on the galactic dynamo.

Shapiro, Paul R.↗
NASA NTRS1994

Towards an MHD theory for the standoff distance of Earth's bow shock

An MHD theory is developed for the stand-off distance a(sub s) of the bow shock and the thickness delta(sub ms) of the magnetosheath, using the empirical Spreiter et al. relation delta(sub ms) = kX and the MHD density ratio X across the shock. The theory includes as special cases the well-known gasdynamic theory and associated phenomenological MHD-like models for delta(sub ms) and a(sub s). In general, however, MHD effects produce major differences from previous models, especially at low Alfven (M(sub A)) and sonic (M(sub S)) Mach numbers. The magnetic field orientation, M(sub A), M(sub S) and the ratio of specific heats gamma are all important variables of the theory. Three principal conclusions are reached. First, the gasdynamic and phenomenological models miss important dependances on field orientation and M(sub S) and generally provide poor approximations to the MHD results. Second, changes in field orientation and M(sub S) are predicted to cause factor of approximately 4 changes in delta(sub ms) at low M(sub A). Third, using Spreiter et al.'s value for k in the MHD theory leads to maximum a(sub s) values at low M(sub A) and nominal M(sub S) that are much smaller than observations and MHD simulations require. Resolving this problem requires either the modified Spreiter-like relation and larger k found in recent MHD simulations and/or breakdown in the Spreiter-like relation at very low M(sub A).

Cairns, Iver H.↗
NASA NTRS1994

Towards an MHD Theory for the Standoff Distance of Earth's Bow Shock

A magnetohydrodynamic (MHD) theory is developed for the standoff distance a(s) of the bow shock and the thickness Delta(ms) of the magnetosheath, using the empirical Spreiter et al. relation Delta(ms) = kX and the MHD density ratio X across the shock. The theory includes as special cases the well-known gasdynamic theory and associated phenomenological MHD-like models for Delta(ms) and As. In general, however, MHD effects produce major differences from previous models, especially at low Alfev (Ma) and Sonic (Ms) Mach numbers. The magnetic field orientation Ma, Ms and the ratio of specific heats gamma are all important variables of the theory. In contrast, the fast mode Mach number need play no direct role. Three principle conclusions are reached. First the gasdynamic and phenomenological models miss important dependences of field orientation and Ms generally provide poor approximations to the MHD results. Second, changes in field orientation and Ms are predicted to cause factor of approximately 4 changes in Delta(ms) at low Ma. These effects should be important when predicting the shock's location or calculating gramma from observations. Third, using Spreiter et al.'s value for k in the MHD theory leads to maxima a(s) values at low Ma and nominal Ms that are much smaller than observations and MHD simulations require. Resolving this problem requires either the modified Spreiter-like relation and larger k found in recent MHD simulations and/or a breakdown in the Spreiter-like relation at very low Ma.

Carins, Iver H.↗
NASA NTRS2006

Plasma Circulation in the Magnetosphere

We investigate the global structure and dynamics of plasma circulation produced by prototypical solar wind disturbances of the interplanetary magnetic field and dynamic pressure. We track the global circulation and energization of solar wind, polar wind, and auroral wind plasmas throughout the magnetosphere, until they precipitate or escape into the downstream solar wind. We use the full equations of motion of the plasma ions within fields produced by a global MHD simulation of the dynamic solar wind interaction. We use the dynamic hot plasma density and Poynting energy flux specified at the inner boundary of the MHD simulation as drivers of conjugate ion outflow fluxes using local empirical relations obtained from the FAST and Polar missions. Birkeland currents computed by the MHD code are used to derive a field-parallel potential drop from a Knight-like relation [as modified by Lyons and Evans, 1980]. This potential drop is applied to each ion as an initial bulk energy, added to a thermal heating driven by the locally incident Poynting flux. The solar wind pressure increase case (B(sub Y) = 5; B(sub z) = 0 nT) produces an immediate substorm owing to compression of pre-existing plasmas. The SB(sub z), interval (embedded in NB(sub z)) produces a substorm after about one hour of development. Both disturbances enhance the auroral wind flux and heavy ion pressure of the magnetosphere substantially, with complex dynamic structuring by auroral acceleration vortexes and dynamic reconnection. Comparisons are made with observations during disturbed periods including the Halloween 2003 super-storm and other periods.

Moore, T. E.↗
NASA NTRS2014

Energetic Electrons in Dipolarization Events: Spatial Properties and Anisotropy

Using the electromagnetic fields of an MHD simulation of magnetotail reconnection, flow bursts, and dipolarization, we further investigate the acceleration of electrons to suprathermal energies. Particular emphasis is on spatial properties and anisotropies as functions of energy and time. The simulation results are compared with Time History of Events and Macroscale Interactions during Substorms observations. The test particle approach successfully reproduces several observed injection features and puts them into a context of spatial maps of the injection region(s): a dominance of perpendicular anisotropies farther down the tail and closer to the equatorial plane, an increasing importance of parallel anisotropy closer to Earth and at higher latitudes, a drop in energy fluxes at energies below approximately 10 keV, coinciding with the plasma density drop, together with increases at higher energy, a triple peak structure of flux increases near 0 deg, 90 deg, and 180 deg, and a tendency of flux increases to extend to higher energy closer to Earth and at lower latitudes. We identified the plasma sheet boundary layers and adjacent lobes as a main source region for both increased and decreased energetic electron fluxes, related to the different effects of adiabatic acceleration at high and low energies. The simulated anisotropies tend to exceed the observed ones, particularly for perpendicular fluxes at high energies. The most plausible reason is that the MHD simulation lacks the effects of anisotropy-driven microinstabilities and waves, which would reduce anisotropies.

anisotropy↗
NASA NTRS2015

The Role of the Hall Effect in Global Structure and Dynamics of Planetary Magnetospheres: Ganymede as a Case Study

We present high-resolution Hall MHD simulations of Ganymede's magnetosphere demonstrating that Hall electric fields in ion-scale magnetic reconnection layers have significant global effects not captured in resistive MHD simulations. Consistent with local kinetic simulations of magnetic reconnection, our global simulations show the development of intense field-aligned currents along the magnetic separatrices. These currents extend all the way down to the moon's surface, where they may contribute to Ganymede's aurora. Within the magnetopause and magnetotail current sheets, Hall J x B forces accelerate ions to the local Alfven speed in the out-of-plane direction, producing a global system of ion drift belts that circulates Jovian magnetospheric plasma throughout Ganymede's magnetosphere. We discuss some observable consequences of these Hall-induced currents and ion drifts: the appearance of a sub-Jovian 'double magnetopause' structure, an Alfvenic ion jet extending across the upstream magnetopause, and an asymmetric pattern of magnetopause Kelvin-Helmholtz waves.

Magnetic Reconnection↗
NASA NTRS2001

Modeling Magnetospheric Sources

We have used global magnetohydrodynamic, simulations of the interaction between the solar wind and magnetosphere together with single particle trajectory calculations to investigate the sources of plasma entering the magnetosphere. In all of our calculations solar wind plasma primarily enters the magnetosphere when the field line on which it is convecting reconnects. When the interplanetary magnetic field has a northward component the reconnection is in the polar cusp region. In the simulations plasma in the low latitude boundary layer (LLBL) can be on either open or closed field lines. Open field lines occur when the high latitude reconnection occurs in only one cusp. In the MHD calculations the ionosphere does not contribute significantly to the LLBL for northward IMF. The particle trajectory calculations show that ions preferentially enter in the cusp region where they can be accelerated by non-adiabatic motion across the high latitude electric field. For southward IMF in the MHD simulations the plasma in the middle and inner magnetosphere comes from the inner (ionospheric) boundary of the simulation. Solar wind plasma on open field lines is confined to high latitudes and exits the tailward boundary of the simulation without reaching the plasma sheet. The LLBL is populated by both ionospheric and solar wind plasma. When the particle trajectories are included solar wind ions can enter the middle magnetosphere. We have used both the MHD simulations and the particle calculations to estimate source rates for the magnetosphere which are consistent with those inferred from observations.

Walker, Raymond J.↗
NASA NTRS1988

The effect of the solenoidal condition on the numerical magnetohydrodynamic simulation of coronal dynamics

Several MHD simulations of coronal dynamics performed more than 10 years ago violated the solenoidality condition for an 'open' magnetic field topology. Using an improved code, consideration is given to the effect on the physical validity of the numerical simulation (for a representative pulse disturbance) for the case when solenoidality is deliberately violated as compared to the case when it is preserved. It is found that the error incurred in this specific case ('open' topology) in the energy density and in the plasma density profiles is rather small and, hence, does not invalidate the earlier conclusions concerning mass and wave motion.

Wu, S. T.↗
NASA NTRS2011

Enhanced Spectral Anisotropies Near the Proton-Cyclotron Scale: Possible Two-Component Structure in Hall-FLR MHD Turbulence Simulations

Recent analysis of the magnetic correlation function of solar wind fluctuations at 1 AU suggests the existence of two-component structure near the proton-cyclotron scale. Here we use two-and-one-half dimensional and three-dimensional compressible MHD models to look for two-component structure adjacent the proton-cyclotron scale. Our MHD system incorporates both Hall and Finite Larmor Radius (FLR) terms. We find that strong spectral anisotropies appear adjacent the proton-cyclotron scales depending on selections of initial condition and plasma beta. These anisotropies are enhancements on top of related anisotropies that appear in standard MHD turbulence in the presence of a mean magnetic field and are suggestive of one turbulence component along the inertial scales and another component adjacent the dissipative scales. We compute the relative strengths of linear and nonlinear accelerations on the velocity and magnetic fields to gauge the relative influence of terms that drive the system with wave-like (linear) versus turbulent (nonlinear) dynamics.

Ghosh, Sanjoy↗
NASA NTRS2004

Fitting Flux Ropes to a Global MHD Solution: A Comparison of Techniques

Flux rope fitting (FRF) techniques are an invaluable tool for extracting information about the properties of a subclass of CMEs in the solar wind. However, it has proven difficult to assess their accuracy since the underlying global structure of the CME cannot be independently determined from the data. In contrast, large-scale MHD simulations of CME evolution can provide both a global view as well as localized time series at specific points in space. In this study we apply 5 different fitting techniques to 2 hypothetical time series derived from MHD simulation results. Independent teams performed the analysis of the events in "blind tests", for which no information, other than the time series, was provided. F rom the results, we infer the following: (1) Accuracy decreases markedly with increasingly glancing encounters; (2) Correct identification of the boundaries of the flux rope can be a significant limiter; and (3) Results from techniques that infer global morphology must be viewed with caution. In spite of these limitations, FRF techniques remain a useful tool for describing in situ observations of flux rope CMEs.

Riley, Pete↗
NASA NTRS2004

Fitting Flux Ropes to a Global MHD Solution: A Comparison of Techniques

Flux rope fitting (FRF) techniques are an invaluable tool for extracting information about the properties of a sub-class of coronal mass ejections (CMEs) in the solar wind. However, it has proven difficult to assess their accuracy since the underlying global structure of the CME cannot be independently determined from the data. In contrast, large-scale MHD simulations of CME evolution can provide both a global view as well as localized time series at specific points in space. In this study we apply five different fitting techniques to two hypothetical time series derived from MHD simulation results. Independent teams performed the analysis of the events in "blind tests", for which no information, other than time series, was provided. From the results, we infer the following: 1) Accuracy decreases markedly with increasingly glancing encounters; 2) Correct identification of the boundaries of the flux rope can be a significant limiter; 3) Results from techniques that infer global morphology must be viewed with caution. In spite of these limitations, FRF techniques remain a useful tool for describing in situ observations of flux rope CMEs.

Riley, Pete↗
NASA NTRS1989

Slow shock characteristics as a function of distance from the X-line in the magnetotail

Both particle and MHD simulations are performed to study the characteristics of slow shocks in the magnetotail. The particle simulations indicate that switch-off shocks exhibit large amplitude rotational wave trains, while magnetotail slow shocks with an intermediate Mach number M(An) less than M(c) of about 0.98 do not display such rotational wave trains. The MHD simulations show that the spontaneous reconnection process in the near-earth plasma sheet leads to the formation of a pair of slow shocks tailward of the reconnection line (X-line). The properties of slow shocks are found to vary as a function of the distance from X-line due to the formation of plasmoid. Slow shocks in most regions of magnetotail are found to be nonswitch-off shocks with M(An) of less than 0.98. The present results are used to discuss the lack of large amplitude rotational wave trains at slow shocks in the deep magnetotail.

Lee, L. C.↗