Engineering topics

Klimas, A. J.

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

At least 19 records

NASA NTRS2005

A Mechanism for the Loading-Unloading Substorm Cycle Missing in MHD Global Magnetospheric Simulation Models

Loading and consequent unloading of magnetic flux is an essential element of the substorm cycle in Earth's magnetotail. We are unaware of an available global MHD magnetospheric simulation model that includes a loading- unloading cycle in its behavior. Given the central role that MHD models presently play in the development of our understanding of magnetospheric dynamics, and given the present plans for the central role that these models will play in ongoing space weather prediction programs, it is clear that this failure must be corrected. A 2-dimensional numerical driven current-sheet model has been developed that incorporates an idealized current- driven instability with a resistive MHD system. Under steady loading, the model exhibits a global loading- unloading cycle. The specific mechanism for producing the loading-unloading cycle will be discussed. It will be shown that scale-free avalanching of electromagnetic energy through the model, from loading to unloading, is carried by repetitive bursts of localized reconnection. Each burst leads, somewhat later, to a field configuration that is capable of exciting a reconnection burst again. This process repeats itself in an intermittent manner while the total field energy in the system falls. At the end of an unloading interval, the total field energy is reduced to well below that necessary to initiate the next unloading event and, thus, a loading-unloading cycle results. It will be shown that, in this model, it is the topology of bursty localized reconnection that is responsible for the appearance of the loading-unloading cycle.

Klimas, A. J.↗

NASA NTRS2005

Intermittent Turbulence and SOC Dynamics in a 2-D Driven Current-Sheet Model

Borovsky et al. have shown that Earth's magnetotail plasma sheet is strongly turbulent. More recently, Borovsky and Funsten have shown that eddy turbulence dominates and have suggested that the eddy turbulence is driven by fast flows that act as jets in the plasma. Through basic considerations of energy and magnetic flux conservation, these fast flows are thought to be localized to small portions of the total plasma sheet and to be generated by magnetic flux reconnection that is similarly localized. Angelopoulos et al., using single spacecraft Geotail data, have shown that the plasma sheet turbulence exhibits signs of intermittence and Weygand et al., using four spacecraft Cluster data, have confirmed and expanded on this conclusion. Uritsky et al., using Polar UVI image data, have shown that the evolution of bright, nightside, UV auroral emission regions is consistent with many of the properties of systems in self-organized criticality (SOC). Klimas et al. have suggested that the auroral dynamics is a reflection of the dynamics of the fast flows in the plasma. sheet. Their hypothesis is that the transport of magnetic fludenergy through the magnetotail is enabled by scale-free avalanches of localized reconnection whose SOC dynamics are reflected in the auroral UV emission dynamics. A corollary of this hypothesis is that the strong, intermittent, eddy turbulence of the plasma sheet is closely related to its critical dynamics. The question then arises: Can in situ evidence for the SOC dynamics be found in the properties of the plasma sheet turbulence? A 2-dimensional numerical driven current-sheet model of the central plasma sheet has been developed that incorporates an idealized current-driven instability with a resistive MHD system. It has been shown that the model can evolve into SOC in a physically relevant parameter regime. Initial results from a study of intermittent turbulence in this model and the relationship of this turbulence to the model's known SOC dynamics will be discussed.

Klimas, A. J.↗

NASA NTRS1999

The Role of Self-Organized Criticality in the Substorm Phenomenon and its Relation to Localized Reconnection in the Plasma Sheet

Recent observations of the magnetotail plasma sheet have shown it to be a dynamic and turbulent region. Research has found strong turbulence in the plasma sheet at approximately 20 Earth's Radius tailward of Earth; the turbulence is observed at all activity levels. The existence of strong turbulence in the plasma sheet in the region associated with substorm onset might be thought difficult to reconcile with the coherence and repeatability of the substorm cycle. We review a variety of evidence that strongly suggests the magnetotail is driven, through magnetic flux transfer, into a state of "self-organized criticality" (SOC). It is an important property of physical systems that evolve into SOC that they self-organize into a unique, global dynamic state. This global state is inevitable, and repeatable. In this state, however, small-spatiotemporal-scale system fluctuations are unpredictable and can be only described statistically. This is the basis, we think, for the global coherence and repeatability of the substorm phenomenon in the turbulent plasma sheet. At, or near, substorm onset the plasma sheet can be described by a global SOC state containing significant small scale turbulence. In several recent studies, "sandpile" models were driven into SOC and then shown to reproduce various measures of substorm activity. We discuss the plasma physical foundation of these sandpile models. The evolution of simple continuum plasma sheet models into SOC-like states of many small reconnection events in the turbulent plasma sheet under the will be demonstrated. We view the substorm phenomenon as an avalanche assumption that the plasma sheet is in a SOC state.

Klimas, A. J.↗

NASA NTRS1997

Comparison of interplanetary type 2 radio burst observations by ISEE-3, Ulysses, and WIND with applications to space weather prediction

Interplanetary (IP) type 2 radio bursts are produced by IP shocks driven by solar ejecta, presumably involving shock acceleration of electrons that leads to radio emission. These radio bursts, which can be detected remotely by a sensitive spacecraft radio receiver, provide a method of tracking the leading edge of solar ejecta moving outward from the sun. Consequently, observations of these bursts sometimes provide advance warning of one or more days prior to the onset of geomagnetic activity induced by the solar ejecta. A robust lower limit on the fraction of intense geomagnetic storms, that are preceded by IP type 2 bursts, is provided. It is shown that 41 percent of the geomagnetic storms occurring during the interval September 1978 to February 1983 were preceded by type 2 events in this catalog, and reasons why the fraction is not larger are addressed. Differences in the observing capabilities of the International Sun-Earth Explorer (ISEE) 3, Ulysses, and WIND, to explain why each of these similar spacecraft radio investigations provides a different perspective of IP type 2 emissions are reviewed.

MacDowall, R. J.↗

NASA NTRS1996

Nonlinear ARMA models for the D(st) index and their physical interpretation

Time series models successfully reproduce or predict geomagnetic activity indices from solar wind parameters. A method is presented that converts a type of nonlinear filter, the nonlinear Autoregressive Moving Average (ARMA) model to the nonlinear damped oscillator physical model. The oscillator parameters, the growth and decay, the oscillation frequencies and the coupling strength to the input are derived from the filter coefficients. Mathematical methods are derived to obtain unique and consistent filter coefficients while keeping the prediction error low. These methods are applied to an oscillator model for the Dst geomagnetic index driven by the solar wind input. A data set is examined in two ways: the model parameters are calculated as averages over short time intervals, and a nonlinear ARMA model is calculated and the model parameters are derived as a function of the phase space.

Vassiliadis, D.↗

NASA NTRS1995

On the uniqueness of linear moving-average filters for the solar wind-auroral geomagnetic activity coupling

The relation between the solar wind input to the magetosphere, VB(sub South), and the auroral geomagnetic index AL is modeled with two linear moving-average filtering methods: linear prediction filters and a driven harmonic oscillator in the form of an electric circuit. Although the response of the three-parameter oscillator is simpler than the filter's, the methods yield similar linear timescales and values of the prediction-observation correlation and the prediction Chi(exp 2). Further the filter responses obtained by the two methods are similar in their long-term features. In these aspects the circuit model is equivalent to linear prediction filtering. This poses the question of uniqueness and proper interpretation of detailed features of the filters such as response peaks. Finally, the variation of timescales and filter responses with the AL activity level is discussed.

Vassiliadis, D.↗

NASA NTRS1995

A description of the solar wind-magnetosphere coupling based on nonlinear filters

A nonlinear filtering method is introduced for the study of the solar wind -- magnetosphere coupling and related to earlier linear techniques. The filters are derived from the magnetospheric state, a representation of the magnetospheric conditions in terms of a few global variables, here the auroral electrojet indices. The filters also couple to the input, a representation of the solar wind variables, here the rectified electric field. Filter-based iterative prediction of the indices has been obtained for up to 20 hours. The prediction is stable with respect to perturbations in the initial magnetospheric state; these decrease exponentially at the rate of 30/min. The performance of the method is examined for a wide range of parameters and is superior to that of other linear and nonlinear techniques. In the magnetospheric state representation the coupling is modeled as a small number of nonlinear equations under a time-dependent input.

Vassiliadis, D.↗

NASA NTRS1994

Substorm recurrence during steady and variable solar wind driving: Evidence for a normal mode in the unloading dynamics of the magnetosphere

Farrugia et al (1993) have recently studied substorm activity driven by the passage of an interplanetary magnetic cloud during which the interplanetary magnetic field turned southward for approximately 18 hours. It was shown that both the epsilon and the VB(sub s) parameters varied slowly on the timescale of a substorm but changed considerably over the interval as a whole. The substorm occurrence rate did not reflect the variation in magnetospheric energy loading rate as measured by these parameters but, rather, remained roughly constant with a 50-min average period. Klimas et al. (1992) showed that the Faraday loop analog model of geomagnetic activity predicts this single unloading rate under various constant loading rates. However, various model parameters were adjusted to yield a 1-hour unloading period in agreement with the Bargatze et al. (1985) linear prediction filters and in approximate agreement with the Farrugia et al. (1993) results. It has since been found necessary to add a slow relaxation mechanism to the Faraday loop model to allow for its approach to a ground state during long periods of inactivity. It is proposed that the relaxation mechanism is provided by slow convection of magnetic flux out of the magnetotail to the dayside magnetosphere. In addition, a rudimentary representation of magnetotail-ionosphere coupling has been added to enable comparison of model output to measured AL. The present study is of the modified Faraday loop model response to solar wind input from the Bargatze et al. data set with comparison of its output to concurrent AL. This study has removed the degree of freedom in parameter choice that had earlier allowed adjustments toward the 1-hour unloading period and has, instead, yielded the 1-hour unloading period under various constant loading rates. It is demonstrated that the second peak of the bimodal Bargatze et al. linear prediction filters at approximately equal 1-hour lag and the approximately constant substorm recurrence rate observed by Farrugia et al. can be interpreted as both being due to the existence of a normal unloading recurrence period in the dynamics of the magnetosphere.

Klimas, A. J.↗

NASA NTRS1994

A splitting algorithm for Vlasov simulation with filamentation filtration

A Fourier-Fourier transformed version of the splitting algorithm for simulating solutions of the Vlasov-Poisson system of equations is introduced. It is shown that with the inclusion of filamentation filtration in this transformed algorithm it is both faster and more stable than the standard splitting algorithm. It is further shown that in a scalar computer environment this new algorithm is approximately equal in speed and far less noisy than its particle-in-cell counterpart. It is conjectured that in a multiprocessor environment the filtered splitting algorithm would be faster while producing more precise results.

Klimas, A. J.↗

NASA NTRS1992

A nonlinear dynamical analogue model of geomagnetic activity

Consideration is given to the solar wind-magnetosphere interaction within the framework of deterministic nonlinear dynamics. An earlier dripping faucet analog model of the low-dimensional solar wind-magnetosphere system is reviewed, and a plasma physical counterpart to that model is constructed. A Faraday loop in the magnetotail is considered, and the relationship of electric potentials on the loop to changes in the magnetic flux threading the loop is developed. This approach leads to a model of geomagnetic activity which is similar to the earlier mechanical model but described in terms of the geometry and plasma contents of the magnetotail. The model is characterized as an elementary time-dependent global convection model. The convection evolves within a magnetotail shape that varies in a prescribed manner in response to the dynamical evolution of the convection. The result is a nonlinear model capable of exhibiting a transition from regular to chaotic loading and unloading. The model's behavior under steady loading and also some elementary forms of time-dependent loading is discussed.

Klimas, A. J.↗

NASA NTRS1991

Examination of time-variable input effects in a nonlinear analogue magnetosphere model

The plasma physical analog model (an extension of the damped, harmonic-oscillator dripping faucet model) is employed to consider explicitly the effect of time-varying the inputs. This work is equivalent to considering the effects of northward and southward turnings of the interplanetary magnetic field for various periods of time. It is found that relatively extended episodes (not less than 2 hours) of turned-on input with shorter (about 1 hour) periods of turned-off input lead to model behavior much like the continuously driven case. Going to short input intervals with longer periods of zero input leads to highly irregular and dramatically fluctuating episodes of magnetotail unloading. These results give an insight into the diversity of apparent magnetospheric responses during relatively isolated substorm conditions. This work shows the absolutely critical interdependence (in a nonlinear dynamical system) of input phasing and internal magnetospheric response cycles.

Baker, D. N.↗

NASA NTRS1991

Linear prediction filters for linear and nonlinear modeled geomagnetic activity

It is shown that the Faraday loop analog model of geomagnetic activity exhibits both the directly driven and loading-unloading magnetospheric responses to solar wind input. It is further shown that the directly driven component is a linear response to loading while the loading-unloading response is nonlinear. Linear prediction filters which relate model input to output are discussed. By either allowing or suppressing the loading-unloading model response filters that relate to nonlinear or linear dynamics, respectively, have been computed. Filters that described the directly driven response are finite ranged; they asymptote to zero with increasing lag on a time scale that is fixed by the dissipation rate of the model. Filters that describe the nonlinear total response are infinite ranged; they asymptote with increasing lag to large amplitude periodic oscillations. Some implications of these infinite ranged filters are discussed.

Klimas, A. J.↗

NASA NTRS1991

Indications of low dimensionality in magnetospheric dynamics

Using three separate but related approaches, the question of whether the dynamic response of the magnetosphere to the solar wind input may be described by a low-order system of equations is examined. First, it is determined that the dimension of the subset (the attractor) in the high-dimensional magnetospheric phase space associated with the westward auroral electrojet (AL) index for some of the data sets compiled by Bargatze et al. (1985) is 4.0 + or - 0.2, seemingly independent of activity level. Second, direct modeling of the magnetosphere, considering the bulk properties of the tail plasma, leads to a system of equations that is similar to those previously reported as a dripping faucet model; here, the focus is specifically on the prediction of a natural frequency in this model. Finally, a peak is identified with the predicted frequency in power spectra of AL computed for intervals with both low and high activity. Peaks at other frequencies also appear in the spectra, and such resonances would be expected for a chaotic nonlinear oscillator. Combining these approaches it is concluded that at least some aspects of magnetospheric dynamics may be meaningfully modeled by low-dimensional sets of equations.

Roberts, D. A.↗

NASA NTRS1990

Three-dimensional analytical model for the spatial variation of the foreshock electron distribution function - Systematics and comparisons with ISEE observations

A model which is consistent with the solar wind and shock surface boundary conditions for the foreshock electron distribution in the absence of wave-particle effects is formulated for an arbitrary location behind the magnetic tangent to the earth's bow shock. Variations of the gyrophase-averaged velocity distribution are compared and contrasted with in situ ISEE observations. It is found that magnetic mirroring of solar wind electrons is the most important process by which nonmonotonic reduced electron distributions in the foreshock are produced. Leakage of particles from the magnetosheath is shown to be relatively unimportant in determining reduced distributions that are nonmonotonic. The two-dimensional distribution function off the magnetic field direction is the crucial contribution in producing reduced distributions which have beams. The time scale for modification of the electron velocity distribution in velocity space can be significantly influenced by steady state spatial gradients in the background imposed by the curved shock geometry.

Fitzenreiter, R. J.↗

NASA NTRS1990

The evolution from weak to strong geomagnetic activity - An interpretation in terms of deterministic chaos

An analogue of the magnetosphere developed on the basis of Shaw's (1984) dripping faucet model was used to model the mechanisms of the magnetospheric response to energy transfer from the solar wind. It is demonstrated that geomagnetic activity results from nonlinearly coupled physical processes and that the strength and the nature of the coupling changes dramatically as the magnetosphere is driven harder and harder by increasing energy input. Based on initial results obtained from the model, is is suggested that a chaotic transition takes place in the analogue system as the loading rate is increased beyond a critical value. This model is able to explain many of the features in the results of linear prediction filtering techniques.

Baker, D. N.↗

NASA NTRS1984

Detection of bump-on-tail reduced electron velocity distributions at the electron foreshock boundary

Reduced velocity distributions are derived from three-dimensional measurements of the velocity distribution of electrons in the 7 to 500 eV range in the electron foreshock. Bump-on-tail reduced distributions are presented for the first time at the foreshock boundary consistent with Filbert and Kellogg's proposed time-of-flight mechanism for generating the electron beams. In a significant number of boundary crossings, bump-on-tail reduced distributions were found in consecutive 3 sec measurements made 9 sec apart. It is concluded that, although the beams are linearly unstable to plasma waves according to the Penrose criterion, they persist on a time scale of 3 to 15 sec. Previously announced in STAR as N84-22518

Fitzenreiter, R. J.↗

NASA NTRS1984

Detection of bump-on-tail reduced electron velocity distributions at the electron foreshock boundary

Fitzenreiter, R. J.↗

NASA NTRS1983

A mechanism for plasma waves at the harmonics of the plasma frequency in the electron foreshock boundary

A bump-on-tail unstable reduced velocity distribution, constructed from data obtained at the upstream boundary of the electron foreshock by the GSFC electron spectrometer experiment on the ISEE-1 satellite, is used as the initial plasma state for a numerical integration of the 1D-Vlasov-Maxwell system of equations. The integration is carried through the growth of the instability, beyond its saturation, and well into the stabilized plasma regime. A power spectrum computed for the electric field of the stabilized plasma is dominated by a narrow peak at the Bohm-Gross frequency of the unstable field mode but also contains significant power at the harmonics of the Bohm-Gross frequency. The harmonic power is in sharp peaks which are split into closely spaced doublets. The fundamental peak at the Bohm-Gross frequency is split into a closely spaced triplet. The mechanism for excitation of the second harmonic is shown to be second order wave-wave coupling. Previously announced in STAR as N83-17315

Klimas, A. J.↗