Search NASA⌕ Search

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.

At least 361 records · Page 20

Filament formation due to photospheric shear

Recently, Wu et al. (1983, 1986) developed a three-dimensional, time-dependent axisymmetric ideal magnetohydrodynamic (MHD) model which has been used to study lower solar atmospheric response due to photospheric motion. They have shown that various observed solar features can be formed from these numerical simulation experiments. Some work currently in progress which has resulted from a specific numerical simulation experiment is reported. These results may elucidate the possible physical mechanism of formation of filament.

Wu, S. T.↗

Vlasov Simulation of Electrostatic Solitary Structures in Multi-Component Plasmas

Electrostatic solitary structures have been observed in the Earth's magnetosheath by the Cluster spacecraft. Recent theoretical work has suggested that these solitary structures are modeled by electron acoustic solitary waves existing in a four-component plasma system consisting of core electrons, two counter-streaming electron beams, and one species of background ions. In this paper, the excitation of electron acoustic waves and the formation of solitary structures are studied by means of a one-dimensional electrostatic Vlasov simulation. The present result first shows that either electron acoustic solitary waves with negative potential or electron phase-space holes with positive potential are excited in four-component plasma systems. However, these electrostatic solitary structures have longer duration times and higher wave amplitudes than the solitary structures observed in the magnetosheath. The result indicates that a high-speed and small free energy source may be needed as a fifth component. An additional simulation of a five-component plasma consisting of a stable four-component plasma and a weak electron beam shows the generation of small and fast electron phase-space holes by the bump-on-tail instability. The physical properties of the small and fast electron phase-space holes are very similar to those obtained by the previous theoretical analysis. The amplitude and duration time of solitary structures in the simulation are also in agreement with the Cluster observation.

magnetohydrodynamics (MHD)↗

Hydrodynamic Stability and Magnetic Reconnection in Disks and Stars

The purpose of this grant is to study parametric instability. The simplest example of parametric instability is a harmonic oscillator with a periodic modulation of the spring constant. If the modulation frequency is close to twice the natural frequency of the oscillator, the amplitude of oscillation tends to grow exponentially. The growth rate is proportional to the strength of the modulation, but it also depends upon the closeness to resonance of the two frequencies, and upon natural damping rate or "Q" of the oscillator. Parametric instabilities are very common in physics. A familiar example is a jogger's ponytail--normally a very strongly damped pendulum, it can be destabilized by the variation in effective gravity during the jogger's stride. Observation confirms that the period of the pendulum is half that of the jogger's vertical motion. In astrophysics, parametric instability may occur by external tidal forcing, or by interaction among eigenmodes. In the latter case, an energetic eigenmode may destabilize modes of half its frequency, provided some weak nonlinearity exists to couple them. Under a previous Astrophysical Theory grant (NAGW-2419), the PI discovered a parametric instability of tidally forced disks such as the accretion disks in cataclysmic variables and X ray binaries [2]. The destabilized modes are tightly-wound, incompressible, three-dimensional waves analogous to g-modes and r-modes in stars. Later work has confirmed our analysis [4]. It was hoped that these modes might provide a source of turbulence and angular momentum transport in accretion disks. However, a follow-up investigation of this instability by local numerical simulations, although confirming the analytically estimated growth rates, found negligible angular momentum flux [3]. Other work, partly supported by the ATP, now strongly indicates that the transport mechanism in such disks is magnetohydrodynamic turbulence [6]. Nevertheless, the parametric mechanism may truncate the outer edges of disks in close binaries [2], and it may be important in disks of very low ionization such as protostellar disks, or even cataclysmic-variable disks in quiescence where the MHD mechanism may be ineffective [5]. All analyses up to 1996 were done in a local approximation where the orbital frequency, shear rate, and tidal field were treated as constants. The locally computed growth rate turns out to depend strongly on radius, and it was unclear how to average these local rates to obtain the correct global rate. This is a critical issue for accretion disks in close binaries, because the local growth rate is comparable to the orbital frequency towards the outer edge of the disk but decreases rapidly inwards. Paper #1 examined this issue in a simplified global model where the destabilizing terms vary with position. We found that the global growth rate is essentially equal to the maximum local rate, provided that the latter is smoothed over a radial range equal to the distance that the destabilized wave propagates at its group speed in one growth time. Thus, in an accretion disk, waves would grow rapidly in the outer parts but would propagate both inwards and outwards at a maximum group speed of order the disk thickness divided by the orbital period.

Goodman, Jeremy↗

A simulation study of impulsive penetration of solar wind irregularities into the magnetosphere at the dayside magnetopause

The impulsive penetration processes that occur when a plasma irregularity in the magnetosheath, modeled as a field-aligned filament, impinges on the dayside magnetopause, are investigated via a 2D magnetohydrodynamic code. It is found that, when the magnetic fields of the magnetosheath and magnetosphere are parallel or antiparallel, even plasma filaments with small initial velocities can eventually penetrate through the magnetopause and enter the magnetosphere. If the magnetic fields are not aligned, penetration will occur only if the initial kinetic energy density of the plasma filament is more than 50 times larger than the magnetic energy density due to the transverse component of the magnetic field.

Ma, Z. W.↗

Comparing Magnetopause Predictions From Two MHD Models During A Geomagnetic Storm and A Quiet Period

Magnetopause location is an important prediction of numerical simulations of the magnetosphere, yet the models can err, either under-predicting or over-predicting the motion of the boundary. This study compares results from two of the most widely used magnetohydrodynamic (MHD) models, the Lyon–Fedder–Mobarry (LFM) model and the Space Weather Modeling Framework (SWMF), to data from the GOES 13 and 15 satellites during the geomagnetic storm on 22 June 2015, and to THEMIS A, D, and E during a quiet period on 31 January 2013. The models not only reproduce the magnetopause crossings of the spacecraft during the storm, but they also predict spurious magnetopause motion after the crossings seen in the GOES data. We investigate the possible causes of the over-predictions during the storm and find the following. First, using different ionospheric conductance models does not significantly alter predictions of the magnetopause location. Second, coupling the Rice Convection Model (RCM) to the MHD codes improves the SWMF magnetopause predictions more than it does for the LFM predictions. Third, the SWMF produces a stronger ring current than LFM, both with and without the RCM and regardless of the LFM spatial resolution. During the non-storm event, LFM predicts the THEMIS magnetopause crossings due to the southward interplanetary magnetic field better than the SWMF. Additionally, increasing the LFM spatial grid resolution improves the THEMIS predictions, while increasing the SWMF grid resolutions does not.

magnetohydrodynamics↗

Resonant absorption of phase-mixed Alfven surface waves in ideal and resistive magnetohydrodynamics - Initial-value problem

An initial-value approach to the study of Alfven surface waves, in which the linearized MHD equations are solved numerically in time and space, is considered. A disturbance of a particular wavelength is excited in a nonuniform plasma and its temporal evolution is simulated with numerical solutions of the ideal and resistive MHD equations in slab geometry. This procedure generates a continuous spectrum and therefore simulates phase-mixing and subsequent resonant absorption. Energy in the phase-mixed surface waves resistively dissipates with the absorption time and width scaling as resistivity to the -1/3 and 1/6 powers, respectively.

Steinolfson, R. S.↗

Constraints on cosmic ray propagation in the galaxy

The goal was to derive a more detailed picture of magnetohydrodynamic turbulence in the interstellar medium and its effects on cosmic ray propagation. To do so, radio astronomical observations (scattering and Faraday rotation) were combined with knowledge of solar system spacecraft observations of MHD turbulence, simulations of wave propagation, and modeling of the galactic distribution to improve the knowledge. A more sophisticated model was developed for the galactic distribution of electron density turbulence. Faraday rotation measure data was analyzed to constrain magnetic field fluctuations in the ISM. VLBI observations were acquired of compact sources behind the supernova remnant CTA1. Simple calculations were made about the energies of the turbulence assuming a direct link between electron density and magnetic field variations. A simulation is outlined of cosmic ray propagation through the galaxy using the above results.

Cordes, James M.↗

Self-similar evolution of the nonlinear magnetic buoyancy instability

A new type of self-similar solution of ideal magnetohydrodynamics (MHD) in the nonlinear stage of the undular model (k parallel to B) of the magnetic buoyancy instability (the ballooning instability in fusion plasma physics or the Parker instability in astrophysics) is found through MHD simulation and theory. The linear theory developed agrees well with the simulation in the early (linear) stage. The nonlinear stages of the instability in the simulation show the self-similar evolution. One of the solutions obtained from the nonlinear analysis has the characteristics of nonlinear instability in Lagrangian coordinates; the fluid velocity and the Alfven speed on each magnetic loop increase exponentially with time, because the loop is evacuated by the field-aligned motion of matter resulting from gravitational acceleration. In the later stage of the nonlinear evolution, the solution property changes from exponential to power-law time dependence. The latter corresponds to a force-free expansion solution. The later saturation of the velocity increment is also discussed.

Shibata, K.↗

Shock Structure and Magnetic Fields Generation Associated with Relativistic Jets Unmagnetized Pair Plasma

Using 3D and 2D particle-in-cell simulations we investigate a shock structure, magnetic field generation, and particle acceleration associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized pair plasma. The simulations use long computational grids which allow to study the formation and dynamics of the system in a spatial and temporal way. We find for the first time a relativistic shock system comparable to a predicted magnetohydrodynamic shock structure consisting of leading and trailing shocks separated by a contact discontinuity. Strong electromagnetic fields resulting from the Weibel two-stream instability are generated in the trailing shock where jet matter is thermalized and decelerated. We analyze the formation and nonlinear development through saturation and dissipation of those fields and associated particle acceleration. In the AGN context the trailing shock corresponds to the jet shock at the head of a relativistic astrophysical jet. In the GRB context this trailing shock can be identified with the bow shock driven by relativistic ejecta. The strong electromagnetic field region in the trailing shock provides the emission site for the hot spot at the leading edge of AGN jets and for afterglow emission from GRBs.

Niemiec, J.↗

The Fundamental Physical Processes Producing and Controlling Stellar Coronal/ Transition-Region/Chromospheric Activity and Structure

Our LTSA grant supports a long-term collaborative investigation of stellar activity. The project involves current NASA spacecraft and supporting ground-based telescopes, will make use of future missions, and utilizes the extensive archives of IUE, ROSAT, HST, and EUVE. Our interests include observational work (with a nonnegligible groundbased component); specialized processing techniques for imaging and spectral data; and semiempirical modeling, ranging from optically-thin emission measure studies to simulations of optically-thick resonance lines. Collaborations with our cool-star colleagues here in Boulder (at JILA and the High Altitude Observatory) provide access to even broader expertise, particularly on the solar corona, convection, and magnetohydrodynamic phenomena (including "dynamo" theories). The broad-brush of our investigation include the following: (1) where do coronae occur in the Hertzsprung-Russell diagram? (2) the winds of coronal stars: hot, cool, or both? (3) age, activity, rotation relations; (4) atmospheric inhomogeneities; and (5) heating mechanisms, subcoronal flows and flares. Our observation task has been to map the global properties of chromospheres and coronae in the H-R diagram and conduct detailed studies of key objects.

Ayres, Thomas R.↗

Presentation Summary: State of the AGN: Progress Toward Understanding Black Hole Accretion Processes

Accretion of plasma onto black holes power some of the most powerfulsystems in the cosmos. Supermassive black holes at the centers ofgalaxies represent the high-mass limit of these objects and so accountfor the most luminous accretors. As a result, their influence spansvast spatial and temporal scales of cosmic phenomena: intraclusterheating, intergalactic media, galactic feedback and star formation,kiloparsec-scale jets/outflows, variability over time scales of minutesto centuries, and luminous multi-wavelength electromagnetic emissionextending all the way down to its event horizon. Their intrigue isheightened by the fact that they lie at the intersection of variousphysical laws---e.g., general relativistic gravity,magnetohydrodynamics, radiation, high-energy particle physics,thermodynamics, and photo-ionization---which all must be reconciled toarrive at a fundamental understanding and probe for new physics, liketests of general relativity. These physics ingredients must beincorporated into simulations performed somehow on dynamical scalesranging from that of the event horizon to parsec-scales. Fortunately,new computational and theoretical techniques---such as GPU computing,radiation transport, and and novel gridding techniques---are enablingprogress to larger scales, more degrees of freedom, and even to binarysystems. Some of the topics we will survey include recent progress onsimulating the relationship between the disk-jet interaction, howtilted black holes behave, radiation-dominated flow, and how binaryAGN affect the standard picture of black hole accretion. Along theway, we will highlight how new technologies have enabled thesescientific rewards. Future directions and open questions will beprovided to inspire discussion and interaction during the session.

Black Holes↗

Study of radiatively sustained cesium plasmas for solar energy conversion

The results of a study aimed at developing a high temperature solar electric converter are reported. The converter concept is based on the use of an alkali plasma to serve as both an efficient high temperature collector of solar radiation as well as the working fluid for a high temperature working cycle. The working cycle is a simple magnetohydrodynamic (MHD) Rankine cycle employing a solid electrode Faraday MHD channel. Research milestones include the construction of a theoretical model for coupling sunlight in a cesium plasma and the experimental demonstration of cesium plasma heating with a solar simulator in excellent agreement with the theory. Analysis of a solar MHD working cycle in which excimer laser power rather than electric power is extracted is also presented. The analysis predicts a positive gain coefficient on the cesium-xenon excimer laser transition.

Palmer, A. J.↗

On the Origin of Near-Radial Magnetic Fields in the Heliosphere: Numerical Simulations

Deviations from the direction of the "Parker spiral" can be found in in situ measurements of the interplanetary magnetic field on essentially all scales. One intriguing subset is the intervals of near-radial magnetic field, lasting for many hours. Some such intervals are obviously associated with coronal mass ejections, while others appear to be embedded within the ambient solar wind. Most occur on declining speed profiles, such that, when mapped back to the Sun, an entire radial field interval appears to have been launched at approximately the same time. It has been proposed that these events are the result of abrupt, semipermanent speed decreases on these field lines close to the Sun, and that such speed changes might be due to interchange reconnection. In this study, we use a three-dimensional, time-dependent magnetohydrodynamic model to assess to what extent this can account for near-radial magnetic fields observed relatively far out in the heliosphere. We find that sudden speed drops on the trailing portions of high-speed flows can produce strongly underwound (that is near radial) field lines in the heliosphere, although significantly larger speed gradients are required than are typically observed. Moreover, the simulations also reproduce the decreases in density, temperature, and magnetic field strength that are also commonly observed within these events. The question of what produces the abrupt speed drops remains to be answered.

Mikic, Zoran↗

Large-scale computations in fluid mechanics; Proceedings of the Fifteenth Summer Seminar on Applied Mathematics, University of California, La Jolla, CA, June 27-July 8, 1983. Parts 1 & 2

Papers are presented on such topics as the use of semi-Lagrangian advective schemes in meteorological modeling; computation with high-resolution upwind schemes for hyperbolic equations; dynamics of flame propagation in a turbulent field; a modified finite element method for solving the incompressible Navier-Stokes equations; computational fusion magnetohydrodynamics; and a nonoscillatory shock capturing scheme using flux-limited dissipation. Consideration is also given to the use of spectral techniques in numerical weather prediction; numerical methods for the incorporation of mountains in atmospheric models; techniques for the numerical simulation of large-scale eddies in geophysical fluid dynamics; high-resolution TVD schemes using flux limiters; upwind-difference methods for aerodynamic problems governed by the Euler equations; and an MHD model of the earth's magnetosphere.

Engquist, B. E.↗

State of the AGN: Progress Toward Understanding Black Hole Accretion Processes

Accretion of plasma onto black holes power some of the most powerful systems in the cosmos. Supermassive black holes at the centers of galaxies represent the high-mass limit of these objects and so account for the most luminous accretors. As a result, their influence spans vast spatial and temporal scales of cosmic phenomena: intracluster heating, intergalactic media, galactic feedback and star formation, kiloparsec-scale jets/outflows, variability over time scales of minutes to centuries, and luminous multi-wavelength electromagnetic emission extending all the way down to its event horizon. Their intrigue is heightened by the fact that they lie at the intersection of various physical laws---e.g., general relativistic gravity, magnetohydrodynamics, radiation, high-energy particle physics, thermodynamics, and photo-ionization---which all must be reconciled to arrive at a fundamental understanding and probe for new physics, like tests of general relativity. These physics ingredients must be incorporated into simulations performed somehow on dynamical scales ranging from that of the event horizon to parsec-scales. Fortunately, new computational and theoretical techniques---such as GPU computing, radiation transport, and and novel gridding techniques---are enabling progress to larger scales, more degrees of freedom, and even to binary systems. Some of the topics we will survey include recent progress on simulating the relationship between the disk-jet interaction, how tilted black holes behave, radiation-dominated flow, and how binary AGN affect the standard picture of black hole accretion. Along the way, we will highlight how new technologies have enabled these scientific rewards. Future directions and open questions will be provided to inspire discussion and interaction during the session.

LISA↗