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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS2014

Onset of Reconnection in the near Magnetotail: PIC Simulations

Using 2.5-dimensional particle-in-cell (PIC) simulations of magnetotail dynamics, we investigate the onset of reconnection in two-dimensional tail configurations with finite Bz. Reconnection onset is preceded by a driven phase, during which magnetic flux is added to the tail at the high-latitude boundaries, followed by a relaxation phase, during which the configuration continues to respond to the driving. We found a clear distinction between stable and unstable cases, dependent on deformation amplitude and ion/electron mass ratio. The threshold appears consistent with electron tearing. The evolution prior to onset, as well as the evolution of stable cases, are largely independent of the mass ratio, governed by integral flux tube entropy conservation as imposed in MHD (magnetohydrodynamics). This suggests that ballooning instability in the tail should not be expected prior to the onset of tearing and reconnection. The onset time and other onset properties depend on the mass ratio, consistent with expectations for electron tearing. At onset,we found electron anisotropies T⊥⁄ T∥ (bottom tail divided by parallel tail) equals 1.1-1.3, raising growth rates and wavenumbers. Our simulations have provided a quantitative onset criterion that is easily evaluated in MHD simulations, provided the spatial resolution is sufficient. The evolution prior to onset and after the formation of a neutral line does not depend on the electron physics, which should permit an approximation by MHD simulations with appropriate dissipation terms.

PIC↗

NASA NTRS1984

Numerical Simulations of Convectively Excited Gravity Waves

Magneto-convection and gravity waves are numerically simulated with a nonlinear, three-dimensional, time-dependent model of a stratified, rotating, spherical fluid shell heated from below. A Solar-like reference state is specified while global velocity, magnetic field, and thermodynamic perturbations are computed from the anelastic magnetohydrodynamic equations. Convective overshooting from the upper (superadiabatic) part of the shell excites gravity waves in the lower (subadiabatic) part. Due to differential rotation and Coriolis forces, convective cell patterns propagate eastward with a latitudinally dependent phase velocity. The structure of the excited wave motions in the stable region is more time-dependent than of the convective motions above. The magnetic field tends to be concentrated over giant-cell downdrafts in the convective zone but is affected very little by the wave motion in the stable region.

Glatzmaier, G. A.↗

NASA NTRS1986

Formation of active region and quiescent prominence magnetic field configurations

To investigate the formation of prominences, researchers studied chromospheric mass injection into an overlying coronal dipole magnetic field using a 2-D ideal magnetohydrodynamic (MHD) numerical model. Researchers propose that active region prominences are formed by chromospheric plasmas injected directly into the overlying coronal magnetic field and that quiescent prominences are formed by plasmas evaporated at the interface between spicules and corona. Hence, for the simulation of an active region prominence magnetic field we inject the mass from one side, but use a symmetric mass injection to form a quiescent prominence field configuration. Researchers try to find optimum conditions for the formation of Kippenhahn-Schuluter(K-S)type field configuration for stable support of the injection plasmas. They find that the formation of K-S type field configuration by mass injection requires a delicate balance between injection velocity, density, and overlying magnetic fields. These results may explain why a prominence does not form on every neutral line.

An, C.-H.↗

NASA NTRS1995

Non-local transport in turbulent MHD convection

The nonlocal non-diffusive transport of passive scalars in turbulent magnetohydrodynamic (MHD) convection is investigated using transilient matrices. These matrices describe the probability that a tracer particle beginning at one position in a flow will be advected to another position after some time. A method for the calculation of these matrices from simulation data which involves following the trajectories of passive tracer particles and calculating their transport statistics, is presented. The method is applied to study the transport in several simulations of turbulent, rotating, three dimensional compressible, penetrative MDH convection. Transport coefficients and other diagnostics are used to quantify the transport, which is found to resemble advection more closely than diffusion. Some of the results are found to have direct relevance to other physical problems, such as the light element depletion in sun-type stars. The large kurtosis found for downward moving particles at the base of the convection zone implies several extreme events.

Miesch, Mark↗

NASA NTRS2016

Modeling the Ionosphere-Thermosphere Response to a Geomagnetic Storm Using Physics-based Magnetospheric Energy Input: OpenGGCM-CTIM Results

The magnetosphere is a major source of energy for the Earth's ionosphere and thermosphere (IT) system. Current IT models drive the upper atmosphere using empirically calculated magnetospheric energy input. Thus, they do not sufficiently capture the storm-time dynamics, particularly at high latitudes. To improve the prediction capability of IT models, a physics-based magnetospheric input is necessary. Here, we use the Open Global General Circulation Model (OpenGGCM) coupled with the Coupled Thermosphere Ionosphere Model (CTIM). OpenGGCM calculates a three-dimensional global magnetosphere and a two-dimensional high-latitude ionosphere by solving resistive magnetohydrodynamic (MHD) equations with solar wind input. CTIM calculates a global thermosphere and a high-latitude ionosphere in three dimensions using realistic magnetospheric inputs from the OpenGGCM. We investigate whether the coupled model improves the storm-time IT responses by simulating a geomagnetic storm that is preceded by a strong solar wind pressure front on August 24, 2005. We compare the OpenGGCM-CTIM results with low-earth-orbit satellite observations and with the model results of Coupled Thermosphere-Ionosphere-Plasmasphere electrodynamics (CTIPe). CTIPe is an up-to-date version of CTIM that incorporates more IT dynamics such as a low-latitude ionosphere and a plasmasphere, but uses empirical magnetospheric input. OpenGGCMCTIM reproduces localized neutral density peaks at approx. 400 km altitude in the high-latitude dayside regions in agreement with in situ observations during the pressure shock and the early phase of the storm. Although CTIPe is in some sense a much superior model than CTIM, it misses these localized enhancements. Unlike the CTIPe empirical input models, OpenGGCM-CTIM more faithfully produces localized increases of both auroral precipitation and ionospheric electric fields near the high-latitude dayside region after the pressure shock and after the storm onset, which in turn effectively heats the thermosphere and causes the neutral density increase at 400 km altitude.

Connor, Hyunju K.↗

NASA NTRS1986

Models for application of radiation boundary condition for MHD waves in collapse calculations

The problem of reflection of magnetohydrodynamic (MHD) waves at the boundary of a numerical grid has to be resolved in order to obtain reliable results for the end state of the (isothermal) collapse of a rotating, magnetic protostellar cloud. Since the goal of investigating magnetic braking in collapse simulations is to see if the transport of angular momentum via alfven waves is large enough to solve the angular momentum problem an approximation that artificially suppresses large amplitudes in the MHD waves can be self-defeating. For this reason, four alternate methods of handling reflected waves where no assumptions are made regarding the amplitudes of the waves were investigated. In order to study this problem (of reflection) without interference from other effects these methods were tried on two simpler cases. The four methods are discussed.

Vanajakshi, C. T.↗

NASA NTRS2000

MHD-Based Specification of Magnetotail Plasma and Fields: Possibilities and Limitations

Magnetohydrodynamics (MHD) constitutes the simplest comprehensive and self-consistent formulation of the properties of space plasmas. As such, it has been applied with large success to the dynamics of solar system plasmas. For the nightside region of the Earth, the magnetotail, MHD simulations have led to new understanding of the structure and dynamics of the plasmas, in response to both changes in boundary conditions as well as internal dynamical processes. As a result, substantial knowledge of the structure and dynamics of the nightside region have been accumulated. In basic MHD conservation laws on magnetospheric structure, as well as the role of non-MHD processes in the initiation and evolution of dynamical processes of the magnetotail. This presentation will consist of three parts, the first of which addresses the basic constraints on magnetotail structure as well as their consequences for magnetotail specification and forecasting. We will then review some examples of magnetotail structural changes brought about by solar-wind-like boundary conditions. Last, we consider the role of non-MHD processes in magnetotail specification and forecasting. The emphasis here will be on inclusion of these processes into MHD models, and their impact on the overall structure and dynamics.

Hesse, M.↗

NASA NTRS2023

Mshpy23: A User-Friendly, Parameterized Model of Magnetosheath Conditions

Lunar Environment heliospheric X-ray Imager (LEXI) and Solar wind – Magnetosphere - Ionosphere Link Explorer (SMILE) will observe magnetosheath and its boundary motion in soft X-rays for understanding magnetopause reconnection modes under varioussolar wind conditions after their respective launches in 2024 and 2025. Magnetosheath conditions, namely, plasma density, velocity, and temperature, are key parameters for predicting and analyzing soft X-ray images from the LEXI and SMILE missions. We developed a user-friendly model of magnetosheath that parameterizes number density, velocity, temperature, and magnetic field by utilizing the global Magnetohydrodynamics (MHD) model as well as the pre-existing gas-dynamic and analytic models. Using this parameterized magnetosheath model, scientists can easily reconstruct expected soft X-ray images and utilize them for analysis of observed images of LEXI and SMILE without simulating the complicated global magnetosphere models. First, we created an MHD based magnetosheath model by running a total of 14 OpenGGCM global MHD simulations under 7 solar wind densities (1, 5, 10, 15, 20, 25, and 30cm−3) and 2 interplanetary magnetic field BZ components (± 4nT), and then parameterizing the results in new magnetosheath conditions. We compared the magnetosheath model result with THEMIS statistical data and it showed good agreement with a weighted Pearson correlation coefficient greater than 0.77, especially for plasma density and plasma velocity. Second, we compiled a suite of magnetosheath models incorporating previous magnetosheath models (gas-dynamic, analytic), and did two case studies to test the performance. The MHD based model was comparable to or better than the previous models while providing self consistency among the magnetosheath parameters. Third, we constructed a tool to calculate a soft X-ray image from any given vantage point, which can support the planning and data analysis of the aforementioned LEXI and SMILE missions. A release of the code has been uploaded to a Github repository.

Magnetosheath↗

NASA NTRS2001

Fluid Aspects of Solar Wind Disturbances Driven by Coronal Mass Ejections

Transient disturbances in the solar wind initiated by coronal eruptions have been modeled for many years, beginning with the self-similar analytical models of Parker and Simon and Axford. The first numerical computer code (one-dimensional, gas dynamic) to study disturbance propagation in the solar wind was developed in the late 1960s, and a variety of other codes ranging from simple one-dimensional gas dynamic codes through three-dimensional gas dynamic and magnetohydrodynamic codes have been developed in subsequent years. For the most part, these codes have been applied to the problem of disturbances driven by fast CMEs propagating into a structureless solar wind. Pizzo provided an excellent summary of the level of understanding achieved from such simulation studies through about 1984, and other reviews have subsequently become available. More recently, some attention has been focused on disturbances generated by slow CMEs, on disturbances driven by CMEs having high internal pressures, and disturbance propagation effects associated with a structured ambient solar wind. Our purpose here is to provide a brief tutorial on fluid aspects of solar wind disturbances derived from numerical gas dynamic simulations. For the most part we illustrate disturbance evolution by propagating idealized perturbations, mimicking different types of CMEs, into a structureless solar wind using a simple one-dimensional, adiabatic (except at shocks), gas dynamic code. The simulations begin outside the critical point where the solar wind becomes supersonic and thus do not address questions of how the CMEs themselves are initiated. Limited to one dimension (the radial direction), the simulation code predicts too strong an interaction between newly ejected solar material and the ambient wind because it neglects azimuthal and meridional motions of the plasma that help relieve pressure stresses. Moreover, the code ignores magnetic forces and thus also underestimates the speed with which pressure disturbances propagate in the wind.

Gosling, J. T.↗

NASA NTRS2003

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The PI (Cranmer) and Co-I (A. van Ballegooijen) made significant progress toward the goal of building a "unified model" of the dominant physical processes responsible for the acceleration of the solar wind. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a one-dimensional and time-independent model of plasma heating and acceleration. The accomplishments in the report period are thus divided into these two categories: 1a. Focused Study of Kinetic MHD Turbulence. We have developed a model of magnetohydrodynamic (MHD) turbulence in the extended solar corona that contains the effects of collisionless dissipation and anisotropic particle heating. A turbulent cascade is one possible way of generating small-scale fluctuations (easy to dissipate/heat) from a pre-existing population of low-frequency Alfven waves (difficult to dissipate/heat). We modeled the cascade as a combination of advection and diffusion in wavenumber space. The dominant spectral transfer occurs in the direction perpendicular to the background magnetic field. As expected from earlier models, this leads to a highly anisotropic fluctuation spectrum with a rapidly decaying tail in the parallel wavenumber direction. The wave power that decays to high enough frequencies to become ion cyclotron resonant depends on the relative strengths of advection and diffusion in the cascade. For the most realistic values of these parameters, though, there is insufficient power to heat protons and heavy ions. The dominant oblique waves undergo Landau damping, which implies strong parallel electron heating. We thus investigated the nonlinear evolution of the electron velocity distributions (VDFs) into parallel beams and discrete phase-space holes (similar to those seen in the terrestrial magnetosphere) which are an alternate means of heating protons via stochastic interactions similar to particle-particle collisions. 1b. Focused Study of the Multi-Mode Detailed Balance Formalism. The PI began to explore the feasibility of using the "weak turbulence," or detailed-balance theory of Tsytovich, Melrose, and others to encompass the relevant physics of the solar wind. This study did not go far, however, because if the "strong" MHD turbulence discussed above is a dominant player in the wind's acceleration region, this formalism is inherently not applicable to the corona. We will continue to study the various published approaches to the weak turbulence formalism, especially with an eye on ways to parameterize nonlinear wave reflection rates. 2. Building the Unified Model Code Architecture. We have begun developing the computational model of a time-steady open flux tube in the extended corona. The model will be "unified" in the sense that it will include (simultaneously for the first time) as many of the various proposed physical processes as possible, all on equal footing. To retain this generality, we have formulated the problem in two interconnected parts: a completely kinetic model for the particles, using the Monte Carlo approach, and a finite-difference approach for the self-consistent fluctuation spectra. The two codes are run sequentially and iteratively until complete consistency is achieved. The current version of the Monte Carlo code incorporates gravity, the zero-current electric field, magnetic mirroring, and collisions. The fluctuation code incorporates WKJ3 wave action conservation and the cascade/dissipation processes discussed above. The codes are being run for various test problems with known solutions. Planned additions to the codes include prescriptions for nonlinear wave steepening, kinetic velocity-space diffusion, and multi-mode coupling (including reflection and refraction).

Cranmer, Steven R.↗

NASA NTRS1994

Coronal heating by the resonant absorption of Alfven waves: The effect of viscous stress tensor

The time-dependent linearized magnetohydrodynamics (MHD) equations for a fully compressible, low-beta, viscoresistive plasma are solved numerically using an implicit integration scheme. The full viscosity stress tensor (Braginskii 1965) is included with the five parameters eta(sub i) i = 0 to 4. In agreement with previous studies, the numerical simulations demonstrate that the dissipation on inhomogeneities in the background Alfven speed occurs in a narrow resonant layer. For an active region in the solar corona the values of eta(sub i) are eta(sub o) = 0.65 g/cm/s, eta(sub 1) = 3.7 x 10(exp -12) g/cm/s, eta(sub 2) = 4 eta(sub 1), eta(sub 3) = 1.4 x 10(exp -6) g/cm/s, eta(sub 4) = 2 eta(sub 3), with n = 10(exp 10)/cu cm, T = 2 x 10(exp 6) K, and B = 100 G. When the Lundquist number S = 10(exp 4) and R(sub 1) much greater than S (where R(sub 1) is the dimensionless shear viscous number) the width of the resistive dissipation layer d(sub r) is 0.22a (where a is the density gradient length scale) and d(sub r) approximately S(exp -1/3). When S much greater than R(sub 1) the shear viscous dissipation layer width d(sub r) scales as R(sub 1)(exp -1/3). The shear viscous and the resistive dissipation occurs in an overlapping narrow region, and the total heating rate is independent of the value of the dissipation parameters in agreement with previous studies. Consequently, the maximum values of the perpendicular velocity and perpendicular magnetic field scale as R(sub 1)(exp -1/3). It is evident from the simulations that for solar parameters the heating due to the compressive viscosity (R(sub 0) = 560) is negligible compared to the resistive and the shear viscous (R(sub 1)) dissipation and it occurs in a broad layer of order a in width. In the solar corona with S approximately equals 10(exp 4) and R(sub 1) approximately equals 10(exp 14) (as calculated from the Braginskii expressions), the shear viscous resonant heating is of comparable magnitude to the resistive resonant heating.

Ofman, L.↗

NASA NTRS2023

The Effect of the Ambient Solar Wind Medium on a CME-driven Shock and the Associated Gradual Solar Energetic Particle Event

We present simulation results of a gradual solar energetic particle (SEP) event detected on 2021 October 9 by multiple spacecraft, including BepiColombo (Bepi) and near-Earth spacecraft such as the Advanced Composition Explorer (ACE). A peculiarity of this event is that the presence of a high-speed stream (HSS) affected the low-energy ion component (≲5 MeV) of the gradual SEP event at both Bepi and ACE, despite the HSS having only a modest solar wind speed increase. Using the EUHFORIA (European Heliospheric FORecasting Information Asset) magnetohydrodynamic model, we replicate the solar wind during the event and the coronal mass ejection (CME) that generated it. We then combine these results with the energetic particle transport model PARADISE (PArticle Radiation Asset Directed at Interplanetary Space Exploration). We find that the structure of the CME-driven shock was affected by the nonuniform solar wind, especially near the HSS, resulting in a shock wave front with strong variations in its properties such as its compression ratio and obliquity. By scaling the emission of energetic particles from the shock to the solar wind compression at the shock, an excellent match between the PARADISE simulation and in situ measurements of ≲5 MeV ions is obtained. Our modeling shows that the intricate intensity variations observed at both ACE and Bepi were influenced by the nonuniform emission of energetic particles from the deformed shock wave and demonstrates the influence of even modest background solar wind structures on the development of SEP events.

Solar energetic particles↗

NASA NTRS1992

The theory of sunspots

This review covers the present state of our theoretical understanding of the physics of sunspots, along with the principal observational results that need to be explained. The topics covered range from the detailed structure of an individual sunspot to the broad connection between sunspots and the global solar magnetic field and the solar cycle. Our aim is to give a critical discussion of the theoretical ideas and models without presenting mathematical details. After outlining the historical development of the basic concepts associated with the magnetohydrodynamic theory of sunspots, we discuss recent treatments of their properties and structure, placing special emphasis on developments that have occurred within the last ten years. There have been remarkable improvements in the theoretical modelling of sunspots, led by new ideas and by more elaborate and realistic numerical simulations. At the same time, new observations have raised new theoretical questions or caused old ones to be reconsidered. In particular, measurements of oscillations in and around sunspots have opened up the new field of sunspot seismology, while recent high-resolution observations have forced us to rethink the structure of a sunspot penumbra.

Thomas, John H.↗

NASA NTRS2002

Interplanetary Magnetic Field Control of the Entry of Solar Energetic Particles into the Magnetosphere

We have investigated the entry of energetic ions of solar origin into the magnetosphere as a function of the interplanetary magnetic field orientation. We have modeled this entry by following high energy particles (protons and 3 He ions) ranging from 0.1 to 50 MeV in electric and magnetic fields from a global magnetohydrodynamic (MHD) model of the magnetosphere and its interaction with the solar wind. For the most part these particles entered the magnetosphere on or near open field lines except for some above 10 MeV that could enter directly by crossing field lines due to their large gyroradii. The MHD simulation was driven by a series of idealized solar wind and interplanetary magnetic field (IMF) conditions. It was found that the flux of particles in the magnetosphere and transport into the inner magnetosphere varied widely according to the IMF orientation for a constant upstream particle source, with the most efficient entry occurring under southward IMF conditions. The flux inside the magnetosphere could approach that in the solar wind implying that SEPs can contribute significantly to the magnetospheric energetic particle population during typical SEP events depending on the state of the magnetosphere.

Richard, R. L.↗

NASA NTRSDate not supplied

Toward Hall MHD Power Generation and Drag Augmentation Using a Coaxial-Electrode Configuration During Hypersonic Entry

Planetary EDL requires that a spacecraft must reduce its relative velocity from multiple km/s to zero for safe landing. The entry phase requires slowing down through an intense hypersonic flight regime (Mach 25 to Mach 5). Magnetohydrodynamics (MHD) is a next-generation concept with a wide untapped set of technologies that can expand planetary exploration capabilities and create more robust entry vehicles by utilizing in-situ resources, increasing control authority, and reducing heat flux. There is therefore a need for an experimentally validated simulation tool for MHD technologies in hypersonic entry.

MHD drag↗

NASA NTRS1994

Enhanced cosmic-ray acceleration rates in highly inclined astrophysical shocks

The theory of cosmic-ray acceleration by subluminal magnetohydrodynamics (MHD) astrophysical shocks according to the test particle approximation is extended to highly oblique shocks where flow speeds which appear nonrelativistic in the shock rest frame appear on transformation to the de Hoffmann-Teller or E identically equal to 0 frame to have upstream flow speeds approaching c. Monte Carlo simulation shows that relative to the predictions of diffusion theory, as the upstream E identically equal to 0 frame velocity approaches c, flatter spectra and faster acceleration rates occur. These spectral and acceleration time changes are similar to those found for relativistic, parallel MHD shocks and may affect all nonthermal active galactic nuclei emission from relativistic electrons. Also there is an additional means of increasing the upper limit of the cosmic-ray spectrum expected from active galactic nuclei, although the approximation used may not accurately reproduce the spectral shape.

Lieu, R.↗

NASA NTRS1987

An adaptive grid, unsteady model for two-dimensional magnetohydrodynamic (MHD) flow

An adaptive grid finite difference method for solving multi-dimensional, time-dependent, magnetohydrodynamic (MHD) equations is developed. The method is capable of solving problems that include high gradients due to geometry, propagation of shock waves and unsteady boundary conditions. The grid generation technique is based on variational principles with direct control over grid concentration, smoothness and skewness. An example for a two-dimensional MHD simulation of the propagation of a solar-flare-generated shock wave in solar wind flow in the heliographic equatorial plane is selected for illustration of this method.

Panitchob, Supat↗

NASA NTRS1988

The role of magnetohydrodynamics in heliospheric space plasma physics research

Magnetohydrodynamics (MHD) is a fairly recent extension of the field of fluid mechanics. While much remains to be done, it has successfully been applied to the contemporary field of heliospheric space plasma research to evaluate the 'macroscopic picture' of some vital topics via the use of conducting fluid equations and numerical modeling and simulations. Some representative examples from solar and interplanetary physics are described to demonstrate that the continuum approach to global problems (while keeping in mind the assumptions and limitations therein) can be very successful in providing insight and large scale interpretations of otherwise intractable problems in space physics.

Dryer, Murray↗