Engineering Papers⌕ Search

Engineering topics

Fedder, J. A.

Publications and source records attributed to Fedder, J. A..

At least 19 records

Global Response to Local Ionospheric Mass Ejection

We revisit a reported "Ionospheric Mass Ejection" using prior event observations to guide a global simulation of local ionospheric outflows, global magnetospheric circulation, and plasma sheet pressurization, and comparing our results with the observed global response. Our simulation framework is based on test particle motions in the Lyon-Fedder-Mobarry (LFM) global circulation model electromagnetic fields. The inner magnetosphere is simulated with the Comprehensive Ring Current Model (CRCM) of Fok and Wolf, driven by the transpolar potential developed by the LFM magnetosphere, and includes an embedded plasmaspheric simulation. Global circulation is stimulated using the observed solar wind conditions for the period 24-25 Sept 1998. This period begins with the arrival of a Coronal Mass Ejection, initially with northward, but later with southward interplanetary magnetic field. Test particles are launched from the ionosphere with fluxes specified by local empirical relationships of outflow to electrodynamic and particle precipitation imposed by the MIlD simulation. Particles are tracked until they are lost from the system downstream or into the atmosphere, using the full equations of motion. Results are compared with the observed ring current and a simulation of polar and auroral wind outflows driven globally by solar wind dynamic pressure. We find good quantitative agreement with the observed ring current, and reasonable qualitative agreement with earlier simulation results, suggesting that the solar wind driven global simulation generates realistic energy dissipation in the ionosphere and that the Strangeway relations provide a realistic local outflow description.

Moore, T. E.↗

How Ionospheric Ions Populate the Magnetosphere during a Magnetic Storm

Ionospheric oxygen ions have been observed throughout the magnetosphere, from the plasma sheet to the ring current region. I t has been found that the O+ /H+ density ratio in the magnetosphere increases with geomagnetic activity and varies with storm phases. During the magnetic storm in late September to earIy October 2002, Cluster was orbiting in the plasma sheet and ring current regions. At prestorm time, Cluster observed high H+ density and low O+ density in the plasma sheet and lobes. During the storm main phase, 0+ density has increased by 10 times over the pre-storm level. Strong field-aligned beams of O+ were observed in the lobes. O+ fluxes were significantly reduced in the central plasma sheet during the storm recovery. However, 0+ was still evident on the boundaries of the plasma sheet and in the lobes. In order to interpret the Cluster observations and to understand how O+ ions populate the magnetosphere during a magnetic storm, we model the storm in early October 2002 using our global ion kinetic simulation (GIK). We use the LFN global simulation model to produce electric and magnetic fields in the outer magnetosphere, the Strangeway outflow scaling with Delcourt ion trajectories to include ionospheric outflows, and the Fok inner magnetospheric model for the plasmaspheric and ring current response to all particle populations. We find that the observed composition features are qualitatively reproduced by the simulations, with some quantitative differences that point to future improvements in the models.

Fok, Mei-Ching↗

Plasma Sheet Circulation Pathways

Global simulations of Earth's magnetosphere in the solar wind compute the pathways of plasma circulation through the plasma sheet. We address the pathways that supply and drain the plasma sheet, by coupling single fluid simulations with Global Ion Kinetic simulations of the outer magnetosphere and the Comprehensive Ring Current Model of the inner magnetosphere, including plasmaspheric plasmas. We find that the plasma sheet is supplied with solar wind plasmas via the magnetospheric flanks, and that this supply is most effective for northward IMF. For southward IMF, the innermost plasma sheet and ring current region are directly supplied from the flanks, with an asymmetry of single particle entry favoring the dawn flank. The central plasma sheet (near midnight) is supplied, as expected, from the lobes and polar cusps, but the near-Earth supply consists mainly of slowly moving ionospheric outflows for typical conditions. Work with the recently developed multi-fluid LFM simulation shows transport via plasma "fingers" extending Earthward from the flanks, suggestive of an interchange instability. We investigate this with solar wind ion trajectories, seeking to understand the fingering mechanisms and effects on transport rates.

Moore, Thomas E.↗

Plasma Circulation in the Magnetosphere

We investigate the global structure and dynamics of plasma circulation produced by prototypical solar wind disturbances of the interplanetary magnetic field and dynamic pressure. We track the global circulation and energization of solar wind, polar wind, and auroral wind plasmas throughout the magnetosphere, until they precipitate or escape into the downstream solar wind. We use the full equations of motion of the plasma ions within fields produced by a global MHD simulation of the dynamic solar wind interaction. We use the dynamic hot plasma density and Poynting energy flux specified at the inner boundary of the MHD simulation as drivers of conjugate ion outflow fluxes using local empirical relations obtained from the FAST and Polar missions. Birkeland currents computed by the MHD code are used to derive a field-parallel potential drop from a Knight-like relation [as modified by Lyons and Evans, 1980]. This potential drop is applied to each ion as an initial bulk energy, added to a thermal heating driven by the locally incident Poynting flux. The solar wind pressure increase case (B(sub Y) = 5; B(sub z) = 0 nT) produces an immediate substorm owing to compression of pre-existing plasmas. The SB(sub z), interval (embedded in NB(sub z)) produces a substorm after about one hour of development. Both disturbances enhance the auroral wind flux and heavy ion pressure of the magnetosphere substantially, with complex dynamic structuring by auroral acceleration vortexes and dynamic reconnection. Comparisons are made with observations during disturbed periods including the Halloween 2003 super-storm and other periods.

Moore, T. E.↗

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

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

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

'Skidding' of the CRRES G-9 barium release

A simulation study and experimental data of the CRRES G-9 ionospheric barium release are presented. The simulation study is based on a 2D electrostatic code that incorporates time-dependent coupling to the background plasma. It is shown that the densest portion of the barium ion cloud 'skids' about 15 km within the first three seconds following the release, consistent with the optical data analyses.

Huba, J. D.↗

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

A parametric study of the solar wind interaction with comets

The Naval Research Laboratory's magnetohydrodynamic simulation code is used to simulate the solar wind interction with comet Halley for two different outgassing rates and several different solar wind states. The magnetic field is more strongly draped for fast solar wind conditions than slow. For higher mass loading rates, the tail becomes wider and contains more magnetic flux. The visual appearance of the comet differs for the case in which the interplanetary magnetic field lies in the plane of the sky from the case when it lies along the line of sight. The ion tail appears shorter in the latter case. Thus variation in the IMF direction can cause significant changes in the appearance of comets. The comet also creates a large momentum flux deficit in the solar wind with a narrow enhanced region within it corresponding to the ion tail.

Russell, C. T.↗

The visual appearance of comets under varying solar wind conditions

Three-dimensional MHD simulations have been performed for four different sets of solar wind conditions and cometary outgassing rates appropriate to the Halley encounters. Even though the simulations are single fluid calculations, it is possible to separate the solar wind and cometary ions using the divergenceless nature of the solar wind ions. The cometary ion density is then integrated along the line-of-sight from the observer through the comet to determine how the comet would look to a distant observer under these different conditions. In general, comet tails appear longer when the interplanetary magnetic field lies in the plane of the sky rather than along the line-of-sight. Also, the tail shrinks as the speed of the solar wind increases and/or the mass loading rate decreases.

Russell, C. T.↗

Simulations of plasma structure evolution in the high-latitude ionosphere

The Naval Research Laboratory has recently developed a two-dimensional inertial, electrostatic code which has been successfully applied to the development and evolution of ionospheric structure driven by plasma instabilities. This code models the ionosphere and magnetosphere as a set of horizontal two-dimensional layers which are coupled by the vertical magnetic field lines at high latitudes. It is shown that the development of instability-generated structure can be strongly dependent on this coupling. For example, the influence of magnetospheric coupling on the E x B gradient drift instability is to retard the instability's growth and to isotropize density irregularities. The influence of ionospheric coupling on the Kelvin-Helmholtz instability is to retard its growth and to suppress vortex formation. This paper presents the results of numerical simulations of these instabilities and discusses their application to high-latitude ionospheric structure.

Huba, J. D.↗

Nonlinear evolution of the Kelvin-Helmholtz instability in the high-latitude ionosphere

The first numerical simulations of the nonlinear evolution of the electrostatic Kelvin-Helmholtz (K-H) instability with ionospheric Perderson conductivity coupling are presented. It is found that the K-H instability develops in a distinctly different manner in the nonlinear regime with Pedersen coupling than without it. Pedersen coupling effects, in conjunction with a neutral wind and density gradient, are shown to result in an increased time scale for K-H instability wave growth, to inhibit K-H vortex formation, to lead to nonlinear structures which can be described as 'breaking waves', and to generate, in the nonlinear regime, small-scale turbulence by means of secondary instabilities growing on primary waves. The spatial power spectra of the electrostatic potential and density fluctuations are computed, and differences with and without Pedersen effects are reported.

Keskinen, M. J.↗

Effect of possible passage through Halley's magnetic tail on geomagnetic activity

It is found that there occurred a geomagnetic storm on May 18/19, 1910 that cannot be dismissed as a recurrent storm. The period of storminess is close to the period of time at which the magnetic and plasma tail of Halley should have passed over the earth. The signature is worldwide and of the form and magnitude expected if the solar wind were shielded from the earth by the comet. Comparison with computer simulations applicable closer to the comet suggests that this interpretation is feasible. If this interpretation is correct, the 'lobes' of the Halley tail have deeper lows at 24 x 10 to the 6th km than in the present model at 7 x 10 to the 6th km downstream and are at least twice as wide.

Russell, C. T.↗

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

Characteristics of cometary picked-up ions in a global model of Giacobini-Zinner

Energetic ions observed during the International Cometary Explorer (ICE) spacecraft flyby of comet Giacobini-Zinner provide information about both the constitution of comets and the plasma physical processes associated with their interaction with the solar wind. In this investigation the details of ion 'pickup,' in the limit where small-scale fluctuations in the plasma and magnetic field are neglected, are modeled by following the motion of a large number of initially cold, heavy (mass 18) ions in a global magnetohydrodynamic model of the local plasma and magnetic field. The results indicate how the background or macroscopic velocity and magnetic field structure of the comet can affect the average spatial and spectral characteristics of the observed cometary ions. These effects, which occur by virtue of forces associated with the compression and the curvature of the magnetic field in the presence of the stagnating plasma flow, can explain the double maxima in the time series of the energetic ion flux observed along the ICE trajectory.

Kimmel, C. D.↗

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

Near-tail reconnection as the cause of cometary tail disconnections

In a cometary tail disconnection event the plasma tail appears to separate from the coma and to accelerate away from it. As this occurs a new tail begins to form. It is proposed that these disconnections arise in a manner analogous to geomagnetic substorms, i.e., by the formation of a strongly reconnecting region in the near tail that forms a magnetic island in the coma and ejects the plasma tail by strengthening the magnetic 'slingshot' within the tail. This reconnection process may be triggered by several different processes, such as interplanetary shocks or variations in the Alfven Mach number.

Russell, C. T.↗