Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical magnetohydrodynamics”

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 199 records · Page 11

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.↗

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↗

Modeling a Carrington-Scale Stellar Superflare and Coronal Mass Ejection from K(exp 1) Cet

Observations from the Kepler mission have revealed frequent superflares on young and active solar-like stars. Superflares result from the large-scale restructuring of stellar magnetic fields, and are associated with the eruption of coronal material (a coronal mass ejection, or CME) and energy release that can be orders of magnitude greater than those observed in the largest solar flares. These catastrophic events, if frequent, can significantly impact the potential habitability of terrestrial exoplanets through atmospheric erosion or intense radiation exposure at the surface. We present results from numerical modeling designed to understand how an eruptive superflare from a young solar-type star, κ(exp 1)Cet, could occur and would impact its astrospheric environment. Our data-inspired, three- dimensional magnetohydrodynamic modeling shows that global-scale shear concentrated near the radial-field polarity inversion line can energize the closed-field stellar corona sufficiently to power a global, eruptive superflare that releases approximately the same energy as the extreme 1859 Carrington event from the Sun. We examine proxy measures of synthetic emission during the flare and estimate the observational signatures of our CME-driven shock, both of which could have extreme space-weather impacts on the habitability of any Earth-like exoplanets. We also speculate that the observed 1986 Robinson-Bopp superflare from κ(exp 1)Cet was perhaps as extreme for that star as the Carrington flare was for the Sun.

Lynch, Benjamin J.↗

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↗

Ion-neutral momentum coupling near discrete high-latitude ionospheric features

A two-dimensional numerical model is developed to study the momentum coupling between the ionosphere and neutral atmosphere in the vicinity of discrete high-latitude features, such as convection channels and plasma density troughs. Based on generalized magnetohydrodynamic equations the model takes account of global pressure gradients, viscous dissipation, ion drag, the Coriolis force, and electrodynamic drifts. Among the findings of an initial steady state investigation are the following: (1) in convection channels, significant shears and rotations of the thermospheric flow can occur below 200 km if a minimum in the electron density profile is present between the E and F regions; (2) in convection channels, the thermospheric wind decreases with height in the F region owing to the effects of horizontal viscosity; and (3) at low altitudes, the boundaries of convection channels may produce Ekman spirals.

St-Maurice, J.-P.↗

FLIP MHD - A particle-in-cell method for magnetohydrodynamics

The fluid-implicit-particle, or 'FLIP' method presently extended to 2D and 3D MHD flow incorporates a Lagrangian field representation and yields a grid magnetic Reynolds number of up to 16 while preserving contact continuities that retain the Galilean invariance of the MHD flow equations. Analytical arguments and numerical examples demonstrate the conservation of mass, momentum, magnetic flux, and energy; 2D calculation results for the illustrative cases of contact discontinuity convection, Rayleigh-Taylor unstable flow.

Brackbill, J. U.↗

Physics of magnetospheric boundary layers

This final report was concerned with the ideas that: (1) magnetospheric boundary layers link disparate regions of the magnetosphere-solar wind system together; and (2) global behavior of the magnetosphere can be understood only by understanding its internal linking mechanisms and those with the solar wind. The research project involved simultaneous research on the global-, meso-, and micro-scale physics of the magnetosphere and its boundary layers, which included the bow shock, the magnetosheath, the plasma sheet boundary layer, and the ionosphere. Analytic, numerical, and simulation projects were performed on these subjects, as well as comparisons of theoretical results with observational data. Other related activity included in the research included: (1) prediction of geomagnetic activity; (2) global MHD (magnetohydrodynamic) simulations; (3) Alfven resonance heating; and (4) Critical Ionization Velocity (CIV) effect. In the appendixes are list of personnel involved, list of papers published; and reprints or photocopies of papers produced for this report.

Cairns, Iver H.↗

MHD Simulations of the Plasma Flow in the Magnetic Nozzle

The magnetohydrodynamic (MHD) flow of plasma through a magnetic nozzle is simulated by solving the governing equations for the plasma flow in the presence of an static magnetic field representing the applied nozzle. This work will numerically investigate the flow and behavior of the plasma as the inlet plasma conditions and magnetic nozzle field strength are varied. The MHD simulations are useful for addressing issues such as plasma detachment and to can be used to gain insight into the physical processes present in plasma flows found in thrusters that use magnetic nozzles. In the model, the MHD equations for a plasma, with separate temperatures calculated for the electrons and ions, are integrated over a finite cell volume with flux through each face computed for each of the conserved variables (mass, momentum, magnetic flux, energy) [1]. Stokes theorem is used to convert the area integrals over the faces of each cell into line integrals around the boundaries of each face. The state of the plasma is described using models of the ionization level, ratio of specific heats, thermal conductivity, and plasma resistivity. Anisotropies in current conduction due to Hall effect are included, and the system is closed using a real-gas equation of state to describe the relationship between the plasma density, temperature, and pressure.A separate magnetostatic solver is used to calculate the applied magnetic field, which is assumed constant for these calculations. The total magnetic field is obtained through superposition of the solution for the applied magnetic field and the self-consistently computed induced magnetic fields that arise as the flowing plasma reacts to the presence of the applied field. A solution for the applied magnetic field is represented in Fig. 1 (from Ref. [2]), exhibiting the classic converging-diverging field pattern. Previous research was able to demonstrate effects such as back-emf at a super-Alfvenic flow, which significantly alters the shape of the magnetic field in both the near- and far-field regions. However, in that work the downstream domain was constrained to a channel of constant cross-sectional area. In the present work we seek to address this issue by modeling the downstream region with a domain that permits free expansion of the plasma, permitting a better evaluation of the downstream effects the applied field has on the plasma. The inlet boundary conditions and applied magnetic field values will also be varied to determine the effect the initial plasma energy content and applied magnetic field energy density have on the near- and far-field plasma properties on the MHD code. This will determine the effect of inlet boundary conditions on the results downstream and address issues related to the restrictive numerical domain previously used.

Smith, T. E. R.↗

Global magnetohydrodynamic simulation of the two-dimensional magnetosphere

The time-dependent magnetohydrodynamic interaction of the solar wind with a two-dimensional dipole magnetic field has been simulated using a novel Lagrangian particle type of MHD code that can treat local low density or vacuum regions without numerical instability. This enables one to simulate the time-dependent magnetic tail. When the solar wind field is southward, a magnetic field line topology consistent with Dungey's model emerges in steady state. The tail, however, is short, and the x-points are only slightly shifted from their vacuum locations, because of strong numerical resistivity. Different configurations resulting from different relative orientations of the solar wind magnetic field and dipole axis are also presented. While the magnetic field is relatively steady, the density and flow in the magnetosheath are turbulent, as are the bow shock and magnetopause; the Kelvin-Helmholtz instability may account for these phenomena. We also model a 'substorm' as the passage of a rotational discontinuity in the solar wind over the dipole. Both 90 and 180 deg shifts to a southward solar wind field cause a violent readjustment of the magnetic tail which eventually settles down to the Dungey configuration.

Leboeuf, J. N.↗

An MHD model of the earth's magnetosphere

It is pointed out that the earth's magnetosphere arises from the interaction of the solar wind with the earth's geomagnetic field. A global magnetohydrodynamics (MHD) model of the earth's magnetosphere has drawn much attention in recent years. In this model, MHD equations are used to describe the solar wind interaction with the magnetosphere. In the present paper, some numerical aspects of the model are considered. Attention is given to the ideal MHD equations, an equation of state for the plasma, the model as an initial- and boundary-value problem, the shock capturing technique, computational requirements and techniques for global MHD modeling, a three-dimensional mesh system employed in the global MHD model, and some computational results.

Wu, C. C.↗

The solar wind flow along the subsolar line in the magnetic barrier and reconnection at the magnetopause

The magnetic barrier is a significant fraction of the magnetosheath region where magnetic forces react back on the plasma flow. We use a magnetohydrodynamic approach to calculate the solar wind flow around the magnetosphere. The use of so-called frozen-in coordinates, where the flow lines and the magnetic field are coordinate axes, is found to be most appropriate. Numerical solutions are presented for the problem of the flow around a blunt body. The results near the magnetopause are used as input parameters for a model of the reconnection of magnetic field lines based on shock-type structures.

Erkaev, N. V.↗

Parallel inhomogeneity and the Alfven resonance. 1: Open field lines

In light of a recent demonstration of the general nonexistence of a singularity at the Alfven resonance in cold, ideal, linearized magnetohydrodynamics, we examine the effect of a small density gradient parallel to uniform, open ambient magnetic field lines. To lowest order, energy deposition is quantitatively unaffected but occurs continuously over a thickened layer. This effect is illustrated in a numerical analysis of a plasma sheet boundary layer model with perfectly absorbing boundary conditions. Consequences of the results are discussed, both for the open field line approximation and for the ensuing closed field line analysis.

Hansen, P. J.↗

Broken Symmetry and Coherent Structure in MHD Turbulence

Absolute equilibrium ensemble theory for ideal homogeneous magnetohydrodynamic (MHD) turbulence is fairly well developed. Theory and Simulation indicate that ideal MHD turbulence non-ergodic and contains coherent structure. The question of applicability real (i.e., dissipative) MHD turbulence is examined. Results from several very long time numerical simulations on a 64(exp 3) grid are presented. It is seen that coherent structure begins to form before decay dominates over nonlinearity. The connection with inverse spectral cascades and selective decay will also be discussed.

Shebalin, John V.↗

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.↗

ON THE ROLE OF INVOLUTIONS IN THE DISCONTINUOUS GALERKIN DISCRETIZATION OF MAXWELL AND MAGNETOHYDRODYNAMIC SYSTEMS

The role of involutions in energy stability of the discontinuous Galerkin (DG) discretization of Maxwell and magnetohydrodynamic (MHD) systems is examined. Important differences are identified in the symmetrization of the Maxwell and MHD systems that impact the construction of energy stable discretizations using the DG method. Specifically, general sufficient conditions to be imposed on the DG numerical flux and approximation space are given so that energy stability is retained These sufficient conditions reveal the favorable energy consequence of imposing continuity in the normal component of the magnetic induction field at interelement boundaries for MHD discretizations. Counterintuitively, this condition is not required for stability of Maxwell discretizations using the discontinuous Galerkin method.

Barth, Timothy↗

Investigating Potential Causes for the Prediction of Spurious Magnetopause Crossings at Geosynchronous Orbit in MHD Simulations

During intense geomagnetic storms, the magnetopause can move in as far as geosynchronous orbit, leaving the satellites in that orbit out in the magnetosheath. Spacecraft operators turn to numerical models to predict the response of the magnetopause to solar wind conditions, but the predictions of the models are not always accurate. This study investigates four storms with a magnetopause crossing by at least one GOES satellite, using four magnetohydrodynamic models at NASA's Community Coordinated Modeling Center to simulate the events, and analyzes the results to investigate the reasons for errors in the predictions. Two main reasons can explain most of the erroneous predictions. First, the solar wind input to the simulations often contains features measured near the L1 point that did not eventually arrive at Earth; incorrect predictions during such periods are due to the solar wind input rather than to the models themselves. Second, while the models do well when the primary driver of magnetopause motion is a variation in the solar wind density, they tend to overpredict or underpredict the integrated Birkeland currents and their effects during times of strong negative interplanetary magnetic field (IMF) Bz, leading to poorer prediction capability. Coupling the MHD codes to a ring current model, when such a coupling is available, generally will improve the predictions but will not always entirely correct them. More work is needed to fully characterize the response of each code under strong southward IMF conditions as it relates to prediction of magnetopause location.

magnetopause↗

Modifications of magnetohydrodynamics as applied to the solar wind

The effect of including the Braginskii viscous stress tensor in magnetohydrodynamics is remarked upon. It is shown that semiquantitative agreement with a recently observed anisotropy in the turbulent solar wind spectrum can be achieved in this way. The modifications of the dynamical equations are simple enough to permit their inclusion in numerical codes. The effects of large 'ion parallel viscosities' also may be significant for plasmas in quite different regimes than the solar wind.

Montgomery, David↗