Search NASA⌕ Search

Engineering topics

Lyon, J. G.

Publications and source records attributed to Lyon, J. G..

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The inviscid incompressible limit of Kelvin–Helmholtz instability for plasmas

The Kelvin–Helmholtz Instability (KHI) is an interface instability that develops between two fluids or plasmas flowing with a common shear layer. KHI occurs in astrophysical jets, solar atmosphere, solar flows, cometary tails, planetary magnetospheres. Two applications of interest, encompassing both space and fusion applications, drive this study: KHI formation at the outer flanks of the Earth’s magnetosphere and KHI growth from non-uniform laser heating in magnetized direct-drive implosion experiments. Here, we study 2D KHI with or without a magnetic field parallel to the flow. We use both the GAMERA code, which solves the compressible Euler equations, and the STRATOSPEC code, which solves the Navier-Stokes equations under the Boussinesq approximation, coupled with the magnetic field dynamics. GAMERA is a global three-dimensional MHD code with high-order reconstruction in arbitrary nonorthogonal curvilinear coordinates, which is developed for a large range of astrophysical applications. STRATOSPEC is a three-dimensional pseudo-spectral code with an accuracy of infinite order (no numerical diffusion). Magnetized KHI is a canonical case for benchmarking hydrocode simulations with extended MHD options. An objective is to assess whether or not, and under which conditions, the incompressibility hypothesis allows to describe a dynamic compressible system. For comparing both codes, we reach the inviscid incompressible regime, by decreasing the Mach number in GAMERA, and viscosity and diffusion in STRATOSPEC. Here, we specifically investigate both single-mode and multi-mode initial perturbations, either with or without magnetic field parallel to the flow. The method relies on comparisons of the density fields, 1D profiles of physical quantities averaged along the flow direction, and scale-by-scale spectral densities. We also address the triggering, formation and damping of filamentary structures under varying Mach number or Atwood number, with or without a parallel magnetic field. Comparisons show very satisfactory results between the two codes. The vortices dynamics is well reproduced, along with the breaking or damping of small-scale structures. We end with the extraction of growth rates of magnetized KHI from the compressible regime to the incompressible limit in the linear regime assessing the effects of compressibility under increasing magnetic field. The observed differences between the two codes are explained either from diffusion or non-Boussinesq effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

Pseudobreakups During January 10, 1997

The January 10, 1997 interplanetary pressure pulse (observed at 0053 UT at Wind) caused a dayside aurora, as seen in Polar Ultraviolet Imager (UVI) data, that propagates tailward and to lower L.

interplanetary pressure pulse polar ultraviolet im↗
NASA NTRS1996

Magnetic Field Orientation Effects on the Standoff Distance of Earth's Bow Shock

Three-dimensional, global MHD simulations of solar wind flow onto a prescribed magnetopause obstacle are used to show that a bow shock's nose location a(sub s), and the relative subsolar magnetosheath thickness Delta(sub ms)/a(sub mp) are strong functions of the IMF cone angle theta (between v(sub sw) and B(sub sw)) and the Alfven Mach number M(sub A). For a given M(sub A) the shock is more distant for higher theta (restricted to the interval 0-90deg by symmetries), while a(sub s)/a(sub mp) and Delta(sub ms/a(sub mp) increase with decreasing M(sub A) for theta greater than or approximately 20deg but decrease with decreasing M(sub A) for theta approximately Odeg. Large differences in Delta(sub ms/a(sub mp) are predicted between theta = Odeg and 90deg at low M(sub A), with smaller differences remaining even at M(sub A) approximately 10. The theta = Odeg results confirm and extend the previous work of Spreiter and Rizzi [1974]. The simulations show that successful models for the subsolar shock location cannot subsume the dependences on M(sub A) and theta into a sole dependence on M(ms). Instead, they confirm a recent prediction [Cairns and Grabbe, 1994] that a(sub s)/a(sub mp) and Delta(sub ms)/a(sub mp) should depend strongly on theta and M(sub A) for M(sub A) less than or approximately 10 (as well as on other MHD variables). Detailed comparisons between theory and data remain to be done. However, preliminary comparisons show good agreement, with distant shock locations found for low M(sub A) and large theta greater than or approximately 45deg and closer locations found for theta less than or approximately 20deg even at M/A approximately 8.

Cairns, Iver H.↗
NASA NTRS1995

MHD Simulations of Earth's Bow Shock at low Mach Numbers: Standoff Distances

Global, three-dimensional, ideal MHD simulations of Earth's bow shock are reported for low Alfven Mach numbers M(sub A) and quasi-perpendicular magnetic field orientations. The simulations use a hard, infinitely conducting magnetopauause obstacle, with axisymmetric three-dimensional location given by scaled standard model, to directly address previous gasdynamic (GD) and field-aligned MHD (FA-MHD) work. Tests of the simulated shocks' density jumps X for 1.4 approx. less than MA approx. less than 10 and the high M(sub A) shock location, and reproduction of the GD relation between magnetosheath thickness and X for quasi-gasdynamic MHD runs with M(sub A) much greater than M(sub s), confirm that the MHD code is working correctly. The MHD simulations show the standoff distance a(sub s), increasing monotonically with decreasing M(sub A). Significantly larger a(sub s), are found at low M(sub A) than predicted by GD and phenomenological MHD models and FA-MHD simulations, as required qualitatively by observations. The GD and FA-MHD predictions err qualitatively, predicting either constant or decreasing a(sub s), with decreasing M(sub A). This qualitative difference between quasi- perpendicular MHD and FA-MHD simulations is direct evidence for a(sub s), depending on the magnetic field orientation Theta. The enhancement factor over the phenomenological MHD predictions at MA approx. 2.4 agrees quantitatively with one observatiorial estimate. A linear relationship is found between the magnetosheath thickness and X, modified both quantitatively and intrinsically by MHD effects from the GD result. The MHD and GD results agree in the high M(sub A) limit. An MHD theory is developed for a(sub s), restricted to sufficiently perpendicular Theta and high sonic Mach numbers M(sub s). It explains the simulation results with excellent accuracy. Observational and further simulation testing of this MHD theory, and of its predicted M(sub A), Theta, and M(sub s) effects, is desirable.

Cairns, Iver H.↗
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 NTRS1995

Topological structure of the magnetotail as a function of interplanetary magnetic field direction

Magnetic reconnection between the interplanetary magnetic field (IMF) and the geomagnetic field is thought to play a major role in the transfer of solar wind momentum and energy to the magnetosphere. As the angle between the IMF and the geomagnetic field is changed at the bow of the magnetosphere, the topological record of the location of the reconnection region should be recorded in the magnetosheath and on the magnetopause along the flanks of the tail, because the super fast flow freezes strong magnetic gradients formed in the bow reconnection regions into the plasma downstream. In this report, we present results from a three-dimensional, magnetohydrodynamic (MHD), global numerical simulation code for the location of the separatrix between unconnected IMF magnetosheath field lines and reconnected field lines which penetrate the magnetopause and connect to the polar ionosphere. The angle between the IMF direction and the line where the separatrix crosses the magnetopause is shown to be a sensitive function of the IMF clock angle. We also explain how this behavior can be used to derive an approximate relation for the dependence of the cross-polar voltage on the IMF clock angle. We conclude with a note of caution concerning the importance of physical boundary conditions in magnetoplasma simulations.

Fedder, J. A.↗
NASA NTRS1995

The Earth's magnetosphere is 165 R(sub E) long: Self-consistent currents, convection, magnetospheric structure, and processes for northward interplanetary magnetic field

The subject of this paper is a self-consistent, magnetohydrodynamic numerical realization for the Earth's magnetosphere which is in a quasi-steady dynamic equilibrium for a due northward interplanetary magnetic field (IMF). Although a few hours of steady northward IMF are required for this asymptotic state to be set up, it should still be of considerable theoretical interest because it constitutes a 'ground state' for the solar wind-magnetosphere interaction. Moreover, particular features of this ground state magnetosphere should be observable even under less extreme solar wind conditions. Certain characteristics of this magnetosphere, namely, NBZ Birkeland currents, four-cell ionospheric convection, a relatively weak cross-polar potential, and a prominent flow boundary layer, are widely expected. Other characteristics, such as no open tail lobes, no Earth-connected magnetic flux beyond 155 R(sub E) downstream, magnetic merging in a closed topology at the cusps, and a 'tadpole' shaped magnetospheric boundary, might not be expected. In this paper, we will present the evidence for this unusual but interesting magnetospheric equilibrium. We will also discuss our present understanding of this singular state.

Fedder, J. A.↗
NASA NTRS1994

MHD simulations of the magnetosheath

Ideal magnetohydrodynamic (MHD) simulations of the dayside magnetosheath are reported. The calculations show a depletion layer near the magnetopause and a new stagnation flow structure. The stagnation flow appears to be a generalization of the stagnation line for ideal MHD. For Mach numbers typical of the solar wind at the Earth, there is no evidence of the slow mode waves recently observed, but for low Mach numbers prominent slow mode waves are seen in the magnetosheath.

Lyon, J. G.↗
NASA NTRS1992

Preliminary study of the CRRES magnetospheric barium releases

Preliminary theoretical and computational analyses of the Combined Release and Radiation Effects Satellite (CRRES) magnetospheric barium releases are presented. The focus of the studies is on the evolution of the diamagnetic cavity which is formed by the barium ions as they expand outward, and on the structuring of the density and magnetic field during the expansion phase of the releases. Two sets of simulation studies are discussed. The first set is based upon a 2D ideal MHD code and provides estimates of the time and length scales associated with the formation and collapse of the diamagnetic cavity. The second set uses a nonideal MHD code; specifically, the Hall term is included. This additional term is critical to the dynamics of sub-Alfvenic plasma expansions, such as the CRRES barium releases, because it leads to instability of the expanding plasma. Detailed simulations of the G4 and G10 releases were performed. In both cases the expanding plasma rapidly structured: the G4 release structured at time t less than about 3 s and developed scale sizes of about 1-2 km, while the G10 release structured at time t less than about 22 s and developed scale sizes of about 10-15 km. It is also found that the diamagnetic cavity size is reduced from those obtained from the ideal MHD results because of the structure. On the other hand, the structuring allows the formation of plasma blobs which appear to free stream across the magnetic field; thus, the barium plasma can propagate to larger distances traverse to the magnetic field than the case where no structuring occurs. Finally, a new normal mode of the system was discovered which may be excited at the leading edge of the expanding barium plasma.

Huba, J. D.↗
NASA NTRS1991

Reconnection voltage as a function of IMF clock angle

Magnetic reconnection between the IMF and the geomagnetic field is thought to play a major role in the transfer of solar wind momentum and energy to the magnetosphere. Both analytic modeling and analysis of geophysical data have shown that this coupling process should be a sensitive function of the clock angle of the IMF. Results are presented from a three-dimensional, MHD, global numerical simulation code for the reconnection voltage between the closed geomagnetic field and the IMF as a function of the IMF clock angle. These results are consistent with a sin(theta/2) functional behavior.

Fedder, J. A.↗
NASA NTRS1987

The solar wind-magnetosphere-ionosphere current-voltage relationship

The global current-voltage relationship for the solar wind-magnetosphere-ionosphere system (SW-M-I) is investigated, restricting the study to strong southward IMF. The dynamo presently identified is on open field lines, and it operates at close to short circuit and at a fraction of the available power output. Control of the dynamo by ionospheric conductivity is discussed. Implications of the simulation results, including the relationship between open and closed field dynamos, the effect of solar wind conditions, the control of reconnection on the bow, the size of the open field line region, and the effects of increased auroral conductivity, are discussed.

Fedder, J. A.↗
NASA NTRS1986

The effect of different resistivity models on magnetotail dynamics

Two-dimensional, time-dependent MHD simulations of the interaction between the solar wind and the earth's magnetosphere have been performed to study magnetotail dynamics with varying forms of anomalous resistivity. In general, the resulting models conform to the neutral line model proposed for substorms with near-earth x points and high-speed tailward flows occurring in the magnetotail. However, in the case where Joule heating is included in the physical description of the system, the near-earth x point never moves far down tail, and high-speed tailward flows never stop. Only in the case where there is no Joule heating does the x point move down tail. Simultaneously, the high-speed tailward flows cease. These results indicate that the mechanism of energy dissipation can have an important effect on reconnection processes and the global magnetospheric dynamics.

Lyon, J. G.↗
NASA NTRS1986

Numerical simulations of comets - Predictions for Comet Giacobini-Zinner

Simulations of Comet Giacobini-Zinner's interaction with solar wind are described and results are presented. The simulations are carried out via the numerical solution of the ideal MHD equations as an initial value problem in a uniform solar wind. The calculations are performed on a Cartesian mesh centered at the comet. Results reveal that the first significant modifications of the solar wind along the ISEE/ICE trajectory will occur 100,000 km from the solar wind comet axis.

Fedder, J. A.↗
NASA NTRS1984

Fast collisionless tearing in an anisotropic neutral sheet

The collisionless tearing mode in a neutral sheet is studied in the presence of ion temperature (Ti) anisotropy, using Vlasov description for both ions and electrons. It is found that the growth rate of the instability is significantly enhanced if the ratio of Ti perpendicular to the equilibrium magnetic field to Ti parallel to the field is greater than one. For typical magnetotail parameters with modest temperature anisotropy, it is shown that the linear e-folding time is reduced to a small fraction of the time delays believed to precede the onset of reconnection. This enhancement of the growth rate is due to the Lorentz force acting on the ions that cross the neutral plane, traversing beyond the conventional electron-tearing layer.

Chen, J.↗
NASA NTRS1981

Computer simulation of a geomagnetic substorm

A global two-dimensional simulation of a substormlike process occurring in earth's magnetosphere is presented. The results are consistent with an empirical substorm model - the neutral-line model. Specifically, the introduction of a southward interplanetary magnetic field forms an open magnetosphere. Subsequently, a substorm neutral line forms at about 15 earth radii or closer in the magnetotail, and plasma sheet thinning and plasma acceleration occur. Eventually the substorm neutral line moves tailward toward its presubstorm position.

Lyon, J. G.↗