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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Parallel electric fields in a simulation of magnetotail reconnection and plasmoid evolution

Properties of the electric field component parallel to the magnetic field are investigate in a 3D MHD simulation of plasmoid formation and evolution in the magnetotail, in the presence of a net dawn-dusk magnetic field component. The spatial localization of E-parallel, and the concept of a diffusion zone and the role of E-parallel in accelerating electrons are discussed. A localization of the region of enhanced E-parallel in all space directions is found, with a strong concentration in the z direction. This region is identified as the diffusion zone, which plays a crucial role in reconnection theory through the local break-down of magnetic flux conservation.

Hesse, M.↗
NASA NTRS1987

Magnetic field draping in the Comet Halley coma - Comparison of Vega observations with computer simulations

During the Vega-1 encounter with Comet Halley, the magnetometer observed draping and compression of the interplanetary magnetic field. These are reproduced well by a three-dimensional MHD simulation of the cometary interaction. Rotations in the magnetic field similar to those at closest approach are also observed 2.75 hours earlier. It is suggested that both rotations correspond to the same IMF interval and that the spacecraft had overtaken the plasma and encountered 'older' magnetic field as it penetrated the coma. Analysis of the MHD model indicates that it should take about 3 to 5 hours for a solar wind parcel to pass from the unperturbed solar wind to Vega-1 at closest approach. A simulated magnetic field profile composed of nested sections for different IMF orientations closely resembles the observations. This result supports the hypothesis of layered magnetic orientations in the coma.

Schwingenschuh, K.↗
NASA NTRS1994

Generation of slow-mode waves in front of the dayside magnetopause

Slow-mode waves have been observed to appear frequently in front of the dayside magnetopause. It is found based on two-dimensional global magnetohydrodynamic (MHD) simulations that slow-mode waves are generated through the interaction between the bow shock and various MHD waves (fast-mode, Alfven-mode, or slow-mode waves) in the upstream solar wind. The generated slow-mode waves stay in front of the magnetopause for a long time (over 15 minutes) before the wave energy is convected away tailward. Since various waves are often present in the solar wind, this mechanism may lead to the frequent appearance of slow-mode waves in front of the magnetopause.

Yan, M.↗
NASA NTRS1987

Numerical tests of a modified full implicit continuous Eulerian (FICE) scheme with projected normal characteristic boundary conditions for MHD flows

A numerical method has been developed based on a modified full implicit continuous Eulerian (FICE) scheme and projected normal characteristic boundary conditions for simulating MHD flows which undergo a long process of evolution. An astrophysical flow is chosen for illustration of this procedure, and numerical tests are made to verify the computational stability and physically realistic solution. Three computational tests have been accomplished; they are tests of solving methods, characteristic boundary condition, and time steps. The tests show that the program from the modified FICE scheme with proper boundary conditions and time steps can be made numerically stable for a time long enough to obtain physically plausible solutions.

Wu, S. T.↗
NASA NTRS1989

Parallel electric fields in a simulation of magnetotail reconnection and plasmoid evolution

Properties of the electric field component parallel to the magnetic field (E sub parallel) in a three-dimensional MHD simulation of plasmoid formation and evolution in the magnetotail in the presence of a net dawn-dusk magnetic field component were observed. Particularly emphasized was the spatial location of E(sub parallel), the concept of a diffusion zone and the role of E(sub parallel) in accelerating electrons. A localization of the region of enhanced E(sub parallel) in all space directions with a strong concentration in the z direction was found. This region was identified as the diffusion zone, which plays a crucial role in reconnection theory through the local break-down of magnetic flux conservation. The presence of B(sub y) implies a north-south asymmetry of the injection of accelerated particles into the near-earth region, if the net B(sub y) field is strong enough to force particles to follow field lines through the diffusion region. It is estimated that for a typical net B(sub y) field this should affect the injection of electrons into the near-earth dawn region, so that precipitation into the Northern (Southern) Hemisphere should dominate for duskward (dawnward) net B(sub y). In addition, a spatial clottiness of the expected injection of adiabatic particles which could be related to the appearance bright spots in auroras was observed.

Hesse, Michael↗
NASA NTRS1993

Three-dimensional magnetohydrodynamics of the emerging magnetic flux in the solar atmosphere

The nonlinear evolution of an emerging magnetic flux tube or sheet in the solar atmosphere is studied through 3D MHD simulations. In the initial state, a horizontal magnetic flux sheet or tube is assumed to be embedded at the bottom of MHD two isothermal gas layers, which approximate the solar photosphere/chromosphere and the corona. The magnetic flux sheet or tube is unstable against the undular mode of the magnetic buoyancy instability. The magnetic loop rises due to the linear and then later nonlinear instabilities caused by the buoyancy enhanced by precipitating the gas along magnetic field lines. We find by 3D simulation that during the ascendance of loops the bundle of flux tubes or even the flux sheet develops into dense gas filaments pinched between magnetic loops. The interchange modes help produce a fine fiber flux structure perpendicular to the magnetic field direction in the linear stage, while the undular modes determine the overall buoyant loop structure. The expansion of such a bundle of magnetic loops follows the self-similar behavior observed in 2D cases studied earlier. Our study finds the threshold flux for arch filament system (AFS) formation to be about 0.3 x 10 exp 20 Mx.

Matsumoto, R.↗
NASA NTRS2003

High Performance Parallel Methods for Space Weather Simulations

This is the final report of our NASA AISRP grant entitled 'High Performance Parallel Methods for Space Weather Simulations'. The main thrust of the proposal was to achieve significant progress towards new high-performance methods which would greatly accelerate global MHD simulations and eventually make it possible to develop first-principles based space weather simulations which run much faster than real time. We are pleased to report that with the help of this award we made major progress in this direction and developed the first parallel implicit global MHD code with adaptive mesh refinement. The main limitation of all earlier global space physics MHD codes was the explicit time stepping algorithm. Explicit time steps are limited by the Courant-Friedrichs-Lewy (CFL) condition, which essentially ensures that no information travels more than a cell size during a time step. This condition represents a non-linear penalty for highly resolved calculations, since finer grid resolution (and consequently smaller computational cells) not only results in more computational cells, but also in smaller time steps.

Hunter, Paul↗
NASA NTRS1989

Computer simulation of Alfven resonance in a cylindrical, axially bounded flux tube

The resonant absorption of Alfven waves in an axially bounded cylindrical flux tube is investigated in a dissipative MHD simulation. It is found that in an axially bounded flux tube, in contrast to an infinite periodic model, the resonant frequency is nearly independent of the poloidal component of the magnetic field. This is a consequence of the 'ballooning' structure of the resonant Alfven waves. The scaling with resistivity and viscosity of the width of the resonance layer, the dissipation rate, and the time for steady state absorption to occur, are all in agreement with theory.

Strauss, H. R.↗
NASA NTRS2010

3D Relativistic Magnetohydrodynamic Simulations of Current-Driven Instability: Instability of a Static Column - 1

We have investigated the development of current-driven (CD) kink instability through three-dimensional relativistic MHD simulations. A static force-free equilibrium helical magnetic configuration is considered in order to study the influence of the initial configuration on the linear and nonlinear evolution of the instability. We found that the initial configuration is strongly distorted but not disrupted by the kink instability. The instability develops as predicted by linear theory. In the non-linear regime the kink amplitude continues to increase up to the terminal simulation time, albeit at different rates, for all but one simulation. The growth rate and nonlinear evolution of the CD kink instability depends moderately on the density profile and strongly on the magnetic pitch profile. The growth rate of the kink mode is reduced in the linear regime by an increase in the magnetic pitch with radius and the non-linear regime is reached at a later time than for constant helical pitch. On the other hand, the growth rate of the kink mode is increased in the linear regime by a decrease in the magnetic pitch with radius and reaches the non-linear regime sooner than the case with constant magnetic pitch. Kink amplitude growth in the non-linear regime for decreasing magnetic pitch leads to a slender helically twisted column wrapped by magnetic field. On the other hand, kink amplitude growth in the non-linear regime nearly ceases for increasing magnetic pitch.

Mizuno, Yosuke↗
NASA NTRS2010

Jovian Plasma Torus Interaction with Europa: 3D Hybrid Kinetic Simulation. First results

The hybrid kinetic model supports comprehensive simulation of the interaction between different spatial and energetic elements of the Europa-moon-magnetosphere system with respect to variable upstream magnetic field and flux or density distributions of plasma and energetic ions, electrons, and neutral atoms. This capability is critical for improving the interpretation of the existing Europa flyby measurements from the Galileo orbiter mission, and for planning flyby and orbital measurements, (including the surface and atmospheric compositions) for future missions. The simulations are based on recent models of the atmosphere of Europa (Cassidy etal.,2007;Shematovichetal.,2005). In contrast to previous approaches with MHD simulations, the hybrid model allows us to fully take into account the finite gyro radius effect and electron pressure, and to correctly estimate the ion velocity distribution and the fluxes along the magnetic field (assuming an initial Maxwellian velocity distribution for upstream background ions).Non-thermal distributions of upstream plasma will be addressed in future work. Photoionization,electron-impact ionization, charge exchange and collisions between the ions and neutrals are also included in our model. We consider two models for background plasma:(a) with O(++) ions; (b) with O(++) and S(++) ions. The majority of O2 atmosphere is thermal with an extended cold population (Cassidyetal.,2007). A few first simulations already include an induced magnetic dipole; however, several important effects of induced magnetic fields arising from oceanic shell conductivity will be addressed in later work.

Lipatov, A. S.↗
NASA NTRS2016

Do We Know the Actual Magnetopause Position for Typical Solar Wind Conditions?

We compare predicted magnetopause positions at the subsolar point and four reference points in the terminator plane obtained from several empirical and numerical MHD (magnetohydrodynamics) models. Empirical models using various sets of magnetopause crossings and making different assumptions about the magnetopause shape predict significantly different magnetopause positions (with a scatter greater than 1 Earth radius (R (sub E)) even at the subsolar point. Axisymmetric magnetopause models cannot reproduce the cusp indentations or the changes related to the dipole tilt effect, and most of them predict the magnetopause closer to the Earth than non axisymmetric models for typical solar wind conditions and zero tilt angle. Predictions of two global non axisymmetric models do not match each other, and the models need additional verification. MHD models often predict the magnetopause closer to the Earth than the non axisymmetric empirical models, but the predictions of MHD simulations may need corrections for the ring current effect and decreases of the solar wind pressure that occur in the foreshock. Comparing MHD models in which the ring current magnetic field is taken into account with the empirical Lin et al. model, we find that the differences in the reference point positions predicted by these models are relatively small for B (sub z) equals 0 (note: B (sub z) is when the Earth's magnetic field points north versus Sun's magnetic field pointing south). Therefore, we assume that these predictions indicate the actual magnetopause position, but future investigations are still needed.

Samsonov, A. A.↗
NASA NTRS1990

Global simulation of the solar wind-magnetosphere interaction - The importance of its numerical validity

This paper examines the global interaction of the solar wind with the geomagnetic field, using a newly developed three-dimensional high-precision MHD simulation code which employs the fourth-order Runge-Kutta-Gill scheme in time and the direct finite-difference method in space, and which allows much less numerical error than the conventional code. The simulation reveals in detail the formation processes of the magnetosphere and the bow shock, as well as the plasma sheet formation process. It is shown that, contrary to the conventional understanding, no Kelvin-Helmholtz instability occurs along the magnetopause in an ideal MHD.

Watanabe, K.↗
NASA NTRS1995

Plasmoid formation and evolution in a numerical simulation of a substorm

Plasmoids are thought to occur as a consequence of the formation of a near-Earth neutral line during the evolution of a geomagnetic substorm. Using a 3D, global MHD simulation of the interaction of the Earth's magnetosphere with the solar wind, we initiate a substorm by a southward turning of the Interplanetary Magnetic Field (IMF) after a long period of steady northward field. A large plasmoid is formed and ejected. We show field line maps of its shape and relate its formation time to the progress of the substorm as indicated by the cross polar potential. Because of the large region of closed field in the magnetotail at the time of the substorm, this plasmoid is longer in axial dimension than is typically observed. We compare the simulation results with the type of satellite observations which have been used to argue for the existence of plasmoids or of traveling compression regions (TCRs) in the lobes or magnetosheath. The simulation predicts that plasmoid passage would result in a strong signal in the cross tail electric field.

Slinker, S. P.↗
NASA NTRS2012

Multi-Scale Modeling of Magnetospheric Dynamics

Magnetic reconnection is a key element in many phenomena in space plasma, e.g. Coronal mass Ejections, Magnetosphere substorms. One of the major challenges in modeling the dynamics of large-scale systems involving magnetic reconnection is to quantifY the interaction between global evolution of the magnetosphere and microphysical kinetic processes in diffusion regions near reconnection sites. Recent advances in small-scale kinetic modeling of magnetic reconnection significantly improved our understanding of physical mechanisms controlling the dissipation in the vicinity of the reconnection site in collisionless plasma. However the progress in studies of small-scale geometries was not very helpful for large scale simulations. Global magnetosphere simulations usually include non-ideal processes in terms of numerical dissipation and/or ad hoc anomalous resistivity. Comparative studies of magnetic reconnection in small scale geometries demonstrated that MHD simulations that included non-ideal processes in terms of a resistive term 11 J did not produce fast reconnection rates observed in kinetic simulations. In collisionless magnetospheric plasma, the primary mechanism controlling the dissipation in the vicinity of the reconnection site is nongyrotropic pressure effects with spatial scales comparable with the particle Larmor radius. We utilize the global MHD code BATSRUS and replace ad hoc parameters such as "critical current density" and "anomalous resistivity" with a physically motivated model of dissipation. The primary mechanism controlling the dissipation in the vicinity of the reconnection site in incorporated into MHD description in terms of non-gyrotropic corrections to the induction equation. We will demonstrate that kinetic nongyrotropic effects can significantly alter the global magnetosphere evolution. Our approach allowed for the first time to model loading/unloading cycle in response to steady southward IMF driving. The role of solar wind parameters and ionospheric conductance on the dynamics of the loading/unloading cycle will be discussed.

Kuznetsova, M. M.↗
NASA NTRS2020

A Double Disturbed Lunar Plasma Wake

Under nominal solar wind conditions, a tenuous wake forms downstream of the lunar nightside. However, the lunar plasma environment undergoes a transformation as the Moon passes through the Earth's magnetotail, with hot subsonic plasma causing the wake structure to disappear. We investigate the lunar wake response during a passing coronal mass ejection (CME) on March 8, 2012 while crossing the Earth's magnetotail using both a magnetohydrodynamic (MHD) model of the terrestrial magnetosphere and a three-dimensional hybrid plasma model of the lunar wake. The CME arrives at 1 AU around 10:30 UT and its impact is first detected inside the geomagnetic tail after 11:10 UT by the Acceleration, Reconnection, Turbulence and Electrodynamics of the Moon's Interaction with the Sun (THEMIS-ARTEMIS) satellites in lunar orbit. A global magnetospheric MHD simulation using Wind data for upstream conditions with the OpenGGCM model reveals the magnetosheath compression to the lunar position from 11:20–12:00 UT, accompanied by multiple flux rope or plasmoid-like features developing and propagating tailward. MHD results support plasma changes observed by the THEMIS-ARTEMIS satellites. Lunar-scale simulations using the Amitis hybrid code show a short and misaligned plasma wake during the Moon's brief entry into the magnetosheath at 11:20 UT, with plasma expansion into the void being aided by the higher plasma temperatures. Sharply accelerated flow speed and a compressed magnetic field lead to an enhanced electric field in the lunar wake capable of generating sudden changes to the nightside near-surface electric potential.

A. P. Rasca↗
NASA NTRS1998

Studies of Tenuous Planetary Atmospheres

The final report includes an overall project overview as well as scientific background summaries of dust and sodium in comets, and tenuous atmospheres of Jupiter's natural satellites. Progress and continuing work related to dust coma and tenuous atmospheric studies are presented. Also included are published articles written during the course of the report period. These are entitled: (1) On Europa's Magnetospheric Interaction: An MHD Simulation; (2) Dust-Gas Interrelations in Comets: Observations and Theory; and (3) Io's Plasma Environment During the Galileo Flyby: Global Three Dimensional MHD Modeling with Adaptive Mesh Refinement.

Combi, Michael R.↗
NASA NTRS2017

Solar Activity Across the Scales: From Small-Scale Quiet-Sun Dynamics to Magnetic Activity Cycles

Observations as well as numerical and theoretical models show that solar dynamics is characterized by complicated interactions and energy exchanges among different temporal and spatial scales. It reveals magnetic self-organization processes from the smallest scale magnetized vortex tubes to the global activity variation known as the solar cycle. To understand these multiscale processes and their relationships, we use a two-fold approach: 1) realistic 3D radiative MHD simulations of local dynamics together with high resolution observations by IRIS, Hinode, and SDO; and 2) modeling of solar activity cycles by using simplified MHD dynamo models and mathematical data assimilation techniques. We present recent results of this approach, including the interpretation of observational results from NASA heliophysics missions and predictive capabilities. In particular, we discuss the links between small-scale dynamo processes in the convection zone and atmospheric dynamics, as well as an early prediction of Solar Cycle 25.

Solar↗
NASA NTRS2017

Solar Activity Across the Scales: from Small-Scale Quiet-Sun Dynamics to Magnetic Activity Cycles

Observations as well as numerical and theoretical models show that solar dynamics is characterized by complicated interactions and energy exchanges among different temporal and spatial scales. It reveals magnetic self-­‐organization processes from the smallest scale magnetized vortex tubes to the global activity variation known as the solar cycle. To understand these multiscale processes and their relationships, we use a two-­fold approach: 1) рealistic 3D radiative MHD simulations of local dynamics together with high-­‐resolution observations by IRIS, Hinode, and SDO; and 2) modeling of solar activity cycles by using simplified MHD dynamo models and mathematical data assimilation techniques. We present recent results of this approach, including the interpretation of observational results from NASA heliophysics missions and predictive capabilities. In particular, we discuss the links between small-­‐scale dynamo processes in the convection zone and atmospheric dynamics, as well as an early prediction of Solar Cycle 25.

Magneti↗