Simulations of Relativistic Jet Formation in Radio Sources
Explore the source record for details and available documents.
astrophysics astronomy black holes magnetohydrodyn↗
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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Three-Dimensional Modeling of the Solar Wind: From the Coronal Base to the Outer Heliosphere
We have developed a global fully three-dimensional magnetohydrodynamic solar wind model for the region that extends from the coronal base to 100 AU. The simulation domain consists of tree spherical shell subdomains with computational boundaries between them placed at 20 solar radii and 0.3 AU. The location of the first boundary ensures that the flow at the boundary is both supersonic and super-Alfvenic. A steady-state solution in the innermost (coronal) region is obtained by the time-relaxation method. The solution uses a tilted dipole model or solar magnetograms as the boundary condition at the coronal base and includes a flux of Alfven waves in the WKB approximation which provide additional acceleration for the coronal outflow in the open field regions. The intermediate region solution is constructed by the integration of steady-state equations along radius using a marching scheme. The outer region solution (0.3-100 AU) is obtained again by the time relaxation and takes into account turbulence transport and heating as well as heating, flow deceleration, and other effects due to the interstellar pickup protons treated as a separate fluid. We use the model to simulate the global steady-state structure of the solar wind from the coronal base to the heliospheric boundary and compare the results with Ulysses and Voyager observations.
Usmanov, A. V↗
Monte Carlo simulation of particle acceleration at astrophysical shocks
A Monte Carlo code was developed for the simulation of particle acceleration at astrophysical shocks. The code is implemented in Turbo Pascal on a PC. It is modularized and structured in such a way that modification and maintenance are relatively painless. Monte Carlo simulations of particle acceleration at shocks follow the trajectories of individual particles as they scatter repeatedly across the shock front, gaining energy with each crossing. The particles are assumed to scatter from magnetohydrodynamic (MHD) turbulence on both sides of the shock. A scattering law is used which is related to the assumed form of the turbulence, and the particle and shock parameters. High energy cosmic ray spectra derived from Monte Carlo simulations have observed power law behavior just as the spectra derived from analytic calculations based on a diffusion equation. This high energy behavior is not sensitive to the scattering law used. In contrast with Monte Carlo calculations diffusive calculations rely on the initial injection of supra-thermal particles into the shock environment. Monte Carlo simulations are the only known way to describe the extraction of particles directly from the thermal pool. This was the triumph of the Monte Carlo approach. The question of acceleration efficiency is an important one in the shock acceleration game. The efficiency of shock waves efficient to account for the observed flux of high energy galactic cosmic rays was examined. The efficiency of the acceleration process depends on the thermal particle pick-up and hence the low energy scattering in detail. One of the goals is the self-consistent derivation of the accelerated particle spectra and the MHD turbulence spectra. Presumably the upstream turbulence, which scatters the particles so they can be accelerated, is excited by the streaming accelerated particles and the needed downstream turbulence is convected from the upstream region. The present code is to be modified to include a better description of particle scattering (pitch-angle instead of hard-sphere) and as iterative procedure for treating the self-excitation of the MHD turbulence.
Campbell, Roy K.↗
Relaxation processes in a low-order three-dimensional magnetohydrodynamics model
The time asymptotic behavior of a Galerkin model of 3D magnetohydrodynamics (MHD) has been interpreted using the selective decay and dynamic alignment relaxation theories. A large number of simulations has been performed that scan a parameter space defined by the rugged ideal invariants, including energy, cross helicity, and magnetic helicity. It is concluded that time asymptotic state can be interpreted as a relaxation to minimum energy. A simple decay model, based on absolute equilibrium theory, is found to predict a mapping of initial onto time asymptotic states, and to accurately describe the long time behavior of the runs when magnetic helicity is present. Attention is also given to two processes, operating on time scales shorter than selective decay and dynamic alignment, in which the ratio of kinetic to magnetic energy relaxes to values 0(1). The faster of the two processes takes states initially dominant in magnetic energy to a state of near-equipartition between kinetic and magnetic energy through power law growth of kinetic energy. The other process takes states initially dominant in kinetic energy to the near-equipartitioned state through exponential growth of magnetic energy.
Stribling, Troy↗
Simulations of Relativistic Jet Formation in Micro-quasars
Explore the source record for details and available documents.
Astrophysics Astronomy black holes magnetohydrodyn↗
Magnetohydrodynamic Turbulence and the Geodynamo
The ARES Directorate at JSC has researched the physical processes that create planetary magnetic fields through dynamo action since 2007. The "dynamo problem" has existed since 1600, when William Gilbert, physician to Queen Elizabeth I, recognized that the Earth was a giant magnet. In 1919, Joseph Larmor proposed that solar (and by implication, planetary) magnetism was due to magnetohydrodynamics (MHD), but full acceptance did not occur until Glatzmaier and Roberts solved the MHD equations numerically and simulated a geomagnetic reversal in 1995. JSC research produced a unique theoretical model in 2012 that provided a novel explanation of these physical observations and computational results as an essential manifestation of broken ergodicity in MHD turbulence. Research is ongoing, and future work is aimed at understanding quantitative details of magnetic dipole alignment in the Earth as well as in Mercury, Jupiter and its moon Ganymede, Saturn, Uranus, Neptune, and the Sun and other stars.
Shebalin, John V.↗
Forced magnetohydrodynamic turbulence in a uniform external magnetic field
Two-dimensional dissipative MHD turbulence is randomly driven at small spatial scales and is studied by numerical simulation in the presence of a strong uniform external magnetic field. A behavior is observed which is apparently distinct from the inverse cascade which prevails in the absence of an external magnetic field. The magnetic spectrum becomes dominated by the three longest wavelength Alfven waves in the system allowed by the boundary conditions: those which, in a box size of edge 2 pi, have wave numbers (kx, ky) = (1, 1), and (1, -1), where the external magnetic field is in the x direction. At any given instant, one of these three modes dominates the vector potential spectrum, but they do not constitute a resonantly coupled triad. Rather, they are apparently coupled by the smaller-scale turbulence.
Hossain, M.↗
Coronal Heating, Weak MHD Turbulence, and Scaling Laws
Long-time high-resolution simulations of the dynamics of a coronal loop in Cartesian geometry are carried out, within the framework of reduced magnetohydrodynamics (RMHD), to understand coronal heating driven by the motion of field lines anchored in the photosphere. We unambiguously identify MHD anisotropic turbulence as the physical mechanism responsible for the transport of energy from the large scales, where energy is injected by photospheric motions, to the small scales, where it is dissipated. As the loop parameters vary, different regimes of turbulence develop: strong turbulence is found for weak axial magnetic fields and long loops, leading to Kolmogorov-like spectra in the perpendicular direction, while weaker and weaker regimes (steeper spectral slopes of total energy) are found for strong axial magnetic fields and short loops. As a consequence we predict that the scaling of the heating rate with axial magnetic field intensity B, which depends on the spectral index of total energy for given loop parameters, must vary from B3/2 for weak fields to B2 for strong fields at a given aspect ratio. The predicted heating rate is within the lower range of observed active region and quiet-Sun coronal energy losses.
corona↗
Self-Consistent Simulation of the Magnetotail
Kinetic effects may be critically important at various stages of magnetotail evolution prior to and after the break up of substorms. The quasi-adiabaticity of ion dynamics results in the appearance of pronounced off diagonal terms in the pressure tensor (an effect that cannot be reproduced using the standard MHD (magnetohydrodynamic) approach). Only the appearance of these terms could support the formation of thin current sheets that recent experimental data have repeatedly indicated. We present in this paper the results of a self-consistent large-scale kinetic simulation of the Earth's magnetotail. We launch ions from the plasma mantle and self-consistently adjust the background magnetic field by computing the perturbations magnetic field from the particle current. We find that a thin current sheet with a half-thickness of approximately 0.3 R(sub E) quickly develops in the simulation. Also, the magnetotail evolves into a periodic state where the location of the X-line in the tail oscillates between x approximately 40 R(sub E) and x approximately 70 R(sub E) on a 5-minute time scale. Evidence of this periodic motion is also seen in the auroral precipitation profile of ions.
Peroomian, Vahe↗
Plasma Sail Concept Fundamentals
The mini-magnetospheric plasma propulsion (M2P2) device, originally proposed by Winglee et al., predicts that a 15-km standoff distance (or 20-km cross-sectional dimension) of the magnetic bubble will provide for sufficient momentum transfer from the solar wind to accelerate a spacecraft to unprecedented speeds of 50 C80 km/s after an acceleration period of 3 mo. Such velocities will enable travel out of the solar system in period of 7 yr almost an order of magnitude improvement over present chemical-based propulsion systems. However, for the parameters of the simulation of Winglee et al., a fluid model for the interaction of M2P2 with the solar wind is not valid. It is assumed in the magnetohydrodynamic (MHD) fluid model, normally applied to planetary magnetospheres, that the characteristic scale size is much greater than the Larmor radius and ion skin depth of the solar wind. In the case of M2P2, the size of the magnetic bubble is actually less than or comparable to the scale of these characteristic parameters. Therefore, a kinetic approach, which addresses the small-scale physical mechanisms, must be used. A two-component approach to determining a preliminary estimate of the momentum transfer to the plasma sail has been adopted. The first component is a self-consistent MHD simulation of the small-scale expansion phase of the magnetic bubble. The fluid treatment is valid to roughly 5 km from the source and the steady-state MHD solution at the 5 km boundary was then used as initial conditions for the hybrid simulation. The hybrid simulations showed that the forces delivered to the innermost regions of the plasma sail are considerably ( 10 times) smaller than the MHD counterpart, are dominated by the magnetic field pressure gradient, and are directed primarily in the transverse direction.
Khazanov, G. V.↗
Magnetohydrodynamic turbulence in the solar wind
The fluctuations in magnetic field and plasma velocity in solar wind, which possess many features of fully developed magnetohydrodynamic (MHD) turbulence, are discussed. Direct spacecraft observations from 0.3 to over 20 AU, remote sensing radio scintillation observations, numerical simulations, and various models provide complementary methods that show that the fluctuations in the wind parameters undergo significant dynamical evolution independent of whatever turbulence might exist in the solar photosphere and corona. The Cluster mission, with high time resolution particle and field measurements and its variable separation strategies, should be able to provide data for answering many questions on MHD turbulence.
Goldstein, Melvyn L.↗
Magnetohydrodynamic Turbulence and the Geodynamo
Recent research results concerning forced, dissipative, rotating magnetohydrodynamic (MHD) turbulence will be discussed. In particular, we present new results from long-time Fourier method (periodic box) simulations in which forcing contains varying amounts of magnetic and kinetic helicity. Numerical results indicate that if MHD turbulence is forced so as to produce a state of relatively constant energy, then the largest-scale components are dominant and quasistationary, and in fact, have an effective dipole moment vector that aligns closely with the rotation axis. The relationship of this work to established results in ideal MHD turbulence, as well as to models of MHD turbulence in a spherical shell will also be presented. These results appear to be very pertinent to understanding the Geodynamo and the origin of its dominant dipole component. Our conclusion is that MHD turbulence, per se, may well contain the origin of the Earth's dipole magnetic field.
Shebalin, John V.↗
Anomalous transport by magnetohydrodynamic Kelvin-Helmholtz instabilities in the solar wind-magnetosphere interaction
The high latitude, or downstream flank, and dayside low latitude magnetospheric boundaries are modeled in an MHD simulation of Kelvin-Helmholtz instablities in a compressible plasma for parallel and transverse configurations. Detailed attention is given to the nonlinear consequences of the instabilities for several different Alfven and sound Mach numbers in both configurations. Emphasis is given to the anomalous transport of momentum and energy by the Kelvin-Helmholtz instabilities across the magnetospheric boundary, which is crucial in evaluating the instabilities' contribution to magnetospheric convection. It is concluded that the anomalous tangential stress at the magnetospheric boundary caused by the instability may be of the order of 1 percent of the magnetosheath flow momentum adjacent to the boundary, and gives a contribution to the convection potential drop over the polar cap of 10-30 kV.
Miura, A.↗
Spatial and Temporal Signatures of Flux Transfer Events in Global Simulations of Magnetopause Dynamics
Flux transfer events (FTEs) were originally identified by transient bipolar variations of the magnetic field component normal to the nominal magnetopause centered on enhancements in the total magnetic field strength. Recent Cluster and THEMIS multi-point measurements provided a wide range of signatures that are interpreted as evidence for FTE passage (e.g., crater FTE's, traveling magnetic erosion regions). We use the global magnetohydrodynamic (MHD) code BATS-R-US developed at the University of Michigan to model the global three-dimensional structure and temporal evolution of FTEs during multi-spacecraft magnetopause crossing events. Comparison of observed and simulated signatures and sensitivity analysis of the results to the probe location will be presented. We will demonstrate a variety of observable signatures in magnetic field profile that depend on space probe location with respect to the FTE passage. The global structure of FTEs will be illustrated using advanced visualization tools developed at the Community Coordinated Modeling Center
Kuznetsova, Maria M.↗
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.↗
Magnetic Fluxtube Tunneling
We present numerical simulations of the collision and subsequent interaction of two initially orthogonal, twisted, force free field magnetic fluxtubes. The simulations were carried out using a new three dimensional explicit parallelized Fourier collocation algorithm for solving the viscoresistive equations of compressible magnetohydrodynamics. It is found that, under a wide range of conditions, the fluxtubes can 'tunnel' through each other. Two key conditions must be satisfied for tunneling to occur: the magnetic field must be highly twisted with a field line pitch much greater than 1, and the magnetic Lundquist number must be somewhat large, greater than or equal to 2880. This tunneling behavior has not been seen previously in studies of either vortex tube or magnetic fluxtube interactions. An examination of magnetic field lines shows that tunneling is due to a double reconnection mechanism. Initially orthogonal field lines reconnect at two specific locations, exchange interacting sections and 'pass' through each other. The implications of these results for solar and space plasmas are discussed.
Dahlburg, Russell B.↗
Modeling Extreme 'Carrington-Type' Space Weather Events Using Three-dimensional Global MHD Simulations
There is a growing concern over possible severe societal consequences related to adverse space weather impacts on man-made technological infrastructure. In the last two decades, significant progress has been made toward the first-principles modeling of space weather events, and three-dimensional (3-D) global magnetohydrodynamics (MHD) models have been at the forefront of this transition, thereby playing a critical role in advancing our understanding of space weather. However, the modeling of extreme space weather events is still a major challenge even for the modern global MHD models. In this study, we introduce a specially adapted University of Michigan 3-D global MHD model for simulating extreme space weather events with a Dst footprint comparable to the Carrington superstorm of September 1859 based on the estimate by Tsurutani et. al., (2003). Results are presented for a simulation run with "very extreme" constructed/idealized solar wind boundary conditions driving the magnetosphere. In particular, we describe the reaction of the magnetosphere-ionosphere system and the associated induced geoelectric field on the ground to such extreme driving conditions. The model setup is further tested using input data for an observed space weather event of Halloween storm October 2003 to verify the MHD model consistence and to draw additional guidance for future work. This extreme space weather MHD model setup is designed specifically for practical application to the modeling of extreme geomagnetically induced electric fields, which can drive large currents in ground-based conductor systems such as power transmission grids. Therefore, our ultimate goal is to explore the level of geoelectric fields that can be induced from an assumed storm of the reported magnitude, i.e., Dst approx. = −1600 nT.
space weather↗
Electron Drift Resonance in the MHD-Coupled Comprehensive Inner Magnetosphere-Ionosphere Model
Relativistic electrons in the outer radiation belt are highly dynamic and respond to interplanetary solar wind structures interacting with the Earth's magnetic field. A known mechanism dictating electron dynamics is the drift-resonant interaction with ultralow frequency (ULF) waves. The present work simulates the ring current and radiation belt electron populations in the bounce-averaged, kinetic Comprehensive Inner Magnetosphere-Ionosphere model coupled with the Block Adaptive Tree Solar Wind Roe-type Upwind Scheme global magnetospheric magnetohydrodynamic (MHD) code using an idealized ULF wave solar wind density driver. ULF waves generated with 10 min periods (at 1.67 mHz frequencies) in the MHD model are characterized and the corresponding energization of electrons and radial transport of electron phase space density is presented. The drift-resonant electron energy is determined in the simulation and is consistent with the electron resonance conditions in dipolar magnetic fields. The present results will be an important component of understanding inner magnetospheric dynamics and how these inner magnetospheric populations interact with ULF waves resulting from interplanetary solar wind structures.
Komar, C. M.↗