Engineering PapersSearch

Engineering topics

Schunk, R. W.

Publications and source records attributed to Schunk, R. W..

At least 37 records · Page 2

Plasmaspheric wind

Observational evidence is presented indicating that beyond L = 1.7-2 plasma corotating in the plasmasphere expands continuously with a small outward directed bulk velocity perpendicular to geomagnetic field lines. A numerical simulation of plasmaspheric flux tube drift motion is presented in support of such an outward plasma expansion. The maximum expansion velocity inside the plasmasphere is determined by the maximum value of the plasma interchange velocity, which is inversely proportional to the value of the integrated Pedersen conductivity.

Lemaire, J.

The flow of plasma in the solar terrestrial environment

The scientific goals of the program are outlined, and some of the papers submitted for publication within the last six months are briefly highlighted. Some of the topics covered include ionosphere-magnetosphere coupling, polar cap arcs, polar wind, convection vortices, ionosphere-plasmasphere coupling, and the validity of macroscopic plasma flow models.

Schunk, R. W.

Ionospheric simulation compared with Dynamics Explorer observations for November 22, 1981

The present study uses an extensive DE-2 data base to both constrain inputs to a time-dependent ionospheric model (TDIM) for a simulation of the ionosphere and then check the simulated densities. The investigation was carried out for both a quiet period and a substorm period. The quiet-day study produced very good agreement between modeled and observed electron densities in the topside ionosphere with two significant exceptions: First, across the polar region the DE-2 LANG densities showed fine structure in addition to the overall regional density morphology. Second, a surprising discrepancy arose in the presunrise and midlatitude trough. The TDIM densities were an order of magnitude lower than those observed by DE-2. The substorm study showed remarkably good agreement with the observed densities.

Sojka, J. J.

Semikinetic and generalized transport models of the polar and solar winds

In order to gain a better understanding of the relative merits of the transport and kinetic approaches to modeling thermal space plasmas, a comparison is presented of a transport (bi-Maxwellian-based 16-moment equations) and a semikinetic description of supersonic flow in the solar wind and also of both supersonic and subsonic flows in the polar wind for 'steady state' conditions. The study shows: remarkable agreement between the two models for supersonic collisionless flows, even for the higher-order moments; the inadequacy of the semikinetic approach for modeling subsonic flows; and the superiority of the 16-moment transport over the semikinetic approach for modeling the solar wind. Further evidence that the bi-Maxwellian-based transport equations are a useful tool for studying 'thermal' space plasmas that develop non-Maxwellian features is provided.

Demars, H. G.

Comparison of the USU ionospheric model with the UCL-Sheffield coupled thermospheric-ionospheric model

Several physical models of the high-latitude ionosphere have been developed that describe the time-dependent evolution of the E- and F-region plasma density. The models require a variety of inputs, including solar EUV fluxes, magnetospheric convection, auroral precipitation, and neutral atmosphere. Of specific relevance to this study is how the neutral atmosphere is incorporated into the ionospheric models. For the USU ionospheric model, the neutral atmosphere is the MSIS 1986 empirical model, while for the UCL-Sheffield coupled thermospheric-ionospheric model the neutral atmosphere is computed simultaneously with the ionosphere. Both models were run for similar solar and magnetospheric conditions (solar maximum, moderate geomagnetic activity, and winter solstice). Solar maximum conditions ensured a strong coupling between the ionosphere and thermosphere, which provided the possibility of a large ionospheric difference between the two physical models. This was further enhanced by choosing winter conditions so that the densities were not dominated by sunlight. The comparison of the two models indicated that both models predict the same morphological features with similar ionospheric densities, generally within about 30 percent.

Sojka, J. J.

Approaches to ionospheric modelling, simulation and prediction

The ionosphere is a complex, multispecies, anisotropic medium that exhibits a significant variation with time, space, season, solar cycle, and geomagnetic activity. In recent years, a wide range of models have been developed in an effort to describe ionospheric behavior. The modeling efforts include: (1) empirical models based on extensive worldwide data sets; (2) simple analytical models for a restricted number of ionospheric parameters; (3) comprehensive, 3D, time-dependent models that require supercomputers; (4) spherical harmonic models based on fits to output obtained from comprehensive numerical models; and (5) ionospheric models driven by real-time magnetospheric inputs. In an effort to achieve simplicity, some of the models have been restricted to certain altitude or latitude domains, while others have been restricted to certain ionospheric parameters, such as the F-region peak density, the auroral conductivity, and the plasma temperatures. The current status of the modeling efforts is reviewed.

Schunk, R. W.

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable, and another one of our goals is to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulation). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, R. W.

Field-aligned current associated with a distorted two-cell convection pattern during northward interplanetary magnetic field

The influence of the ionospheric conductance on the field-aligned current associated with a distorted two-cell convection pattern during northward IMF was investigated using the Heppner-Maynard (1987) convection model and the Utah State University conductivity model described by Rasmussen and Schunk (1987). Results show that the variation of the ionospheric conductivity distribution can significantly affect the features of the field-aligned current for northward IMF, where matching or mismatching between the conductance gradient and the convection electric field plays a key role. It was found that the increase of the field-aligned current in the polar cap observed during summer is mainly due to the increasing contribution from the Pedersen current, and that the increase of the field-aligned current in both the oval region and the evening-midnight sector during the active aurora period is mainly due to the increasing contribution from the Hall current.

Zhu, L.

On the discontinuity in kinetic solutions of the collisionless polar wind

The use of a Maxwellian ion velocity distribution in kinetic solutions of the collisionless polar wind regime results in a velocity distribution and associated moments that display sharp discontinuities at the boundary, which quantitatively affects the solution at high altitudes. This paper investigates the causes of this discontinuity, using both analytical and numerical approaches; it is shown that the choice of appropriate boundary conditions will eliminate the discontinuity.

Khoyloo, A.

Latitudinal dynamics of steady solar wind flows

In many previous studies, it has been assumed that the streamlines of the solar wind are in a purely radial direction beyond a few solar radii. Therefore, the nonradial modulations by the magnetic force and gas pressure is neglected at large heliocentric distances. The work reported in this paper includes a two-dimensional magnetohydrodynamic (MHD) study of the latitudinal structure of solar wind flows. The dynamical effects of the interplanetary magnetic field (IMF) are investigated by MHD solutions and by comparing them with both one-dimensional and two-dimensional hydrodynamic (HD) solutions. The corotational effect of the sun and the influence of transient events in the solar atmosphere are neglected. With such a simplification, the latitudinal structure of the steady solar wind is controlled by the dynamic balance of the magnetic force and the gas pressure. Several different patterns of the evolution of the solar wind latitudinal structure are exhibited in the calculation. The results indicate the existence of a proton number density maximum at the magnetic neutral line, whether or not there is a density maximum or minimum at the inner boundary. The drift motion of magnetic field lines toward the magnetic neutral line enhances the magnetic field strength around the neutral sheet, which may provide a possible explanation of the discrepancy between the measured IMF at 1 AU and that extrapolated from the photospheric magnetic field by the current source-surface modeling.

Yang, W.-H.

Comparison of semikinetic and generalized transport models of the polar wind

A rigorous comparison is made of a semikinetic model and the bi-Maxwellian-based 16-moment transport equations for polar wind conditions. The comparisons show remarkably close agreement in the corresponding predictions for the altitudinal variation of the various plasma parameters. These results indicate that the bi-Maxwellian-based transport equations are a powerful tool for modeling a wide range of thermal plasma flows throughout the solar-terrestrial environment. In the collisionless limit, the results indicate that the 16-moment transport equations are equivalent to the kinetic models commonly used in solar and polar wind studies.

Demars, H. G.

Solutions to bi-Maxwellian transport equations for radial solar wind beyond 28 R(S)

This paper presents solar wind solutions for radial flow between 28 solar radii and 1 AU using the bi-Maxwellian-based 16-moment set of transport equations. In addition to the number density, drift velocity, and parallel and perpendicular temperatures, the 16-moment equations account for the transport of both longitudinal and transverse thermal energies as well as stress. Also, using the 16-moment approximation for the distribution function and assuming plasma parameter values characteristic of the solar wind, contour plots are generated for the proton velocity distribution function. It is shown how the shape of these plots depends on various macroscopic plasma parameters.

Demars, H. G.

Plasma expansion characteristics of ionized clouds in the ionosphere - Macroscopic formulation

A macroscopic plasma expansion model, based on a numerical solution of the type-dependent nonlinear coupled continuity and momentum equations for background O(+) ions and several released ion species, was used to examine plasma expansion characteristics of ionized clouds in the ionosphere. Information is obtained on Ba(+), Li(+), and Ba(+) - Li(+) clouds; different cloud sizes (Gaussian half widths of 0.1, 0.5, and 1 km); cloud/background ion density ratios covering two orders of magnitude; electron/ion temperature ratios of 1, 5, and 10; and several cloud-background relative velocities along B (0, 1, 2, 4, 6, and 8 km/sec). The macroscopic expansion features obtained were found to be in general agreement with those obtained from the small-scale numerical simulations.

Schunk, R. W.

Model and observation comparison of the universal time and IMF by dependence of the ionospheric polar hole

The polar ionospheric F-region often exhibits regions of marked density depletion. These depletions have been observed by a variety of polar orbiting ionospheric satellites over a full range of solar cycle, season, magnetic activity, and universal time (UT). An empirical model of these observations has recently been developed to describe the polar depletion dependence on these parameters. Specifically, the dependence has been defined as a function of F10.7 (solar), summer or winter, Kp (magnetic), and UT. Polar cap depletions have also been predicted /1, 2/ and are, hence, present in physical models of the high latitude ionosphere. Using the Utah State University Time Dependent Ionospheric Model (TDIM) the predicted polar depletion characteristics are compared with those described by the above empirical model. In addition, the TDIM is used to predict the IMF By dependence of the polar hole feature.

Sojka, J. J.

Particle-in-cell simulations of sheath formation around biased interconnectors in a low-earth-orbit plasma

The interaction between satellite solar arrays and the LEO plasma is presently studied with particle-in-cell simulations in which an electrical potential was suddenly applied to the solar cell interconnector. The consequent temporal response was followed for the real O(+)-electron mass ratio in the cases of 100- and 250-V solar cells, various solar cell thicknesses, and solar cells with secondary electron emission. Larger applied potentials and thinner solar cells lead to greater initial polarization surface charges, and therefore longer discharging and shielding times. When secondary electron emission from the cover glass is brought to bear, however, the potential structure is nearly planar, allowing constant interaction between plasma electrons and cover glass; a large fraction of the resulting secondary electrons is collected by the interconnector, constituting an order-of-magnitude increase in collected current.

Thiemann, H.

The longitude dependence of the dayside F region trough - A detailed model-observation comparison

The nighttime main F-region trough extends into the sunlit afternoon sector. This trough feature exhibits both a strong magnetic-activity dependence and a longitude (UT) dependence. Whalen (1987), using IGY ionosonde data, showed that both of these effects are readly extracted from f0F2 observations. This study shows that the longitude effect is the same as that contained in the Utah State University time-dependent ionospheric model. It arises from the offset of the geomagnetic axis from the geographic axis. The magnetic-activity dependency is associated with the westward convection in the afternoon sector. It is also contained in the ionospheric model via the empirical magnetospheric convection model.

Sojka, J. J.

A test of convection models for IMF Bz north

The Utah State University Ionospheric Model was run to obtain diurnally reproducible ionospheric densities and temperatures for summer and winter conditions using both distorted two-cell and three-cell convection patterns. Differences due to the different convection patterns manifest themselves in the depth and location of polar holes in the F-region electron density. While the total depth of the model holes is a characteristic of the diurnally reproducible pattern, the features appear and are recognizable within 0.5 h. Langmuir probe data from 41 DE-2 passes, during which the IMF Bz component was northward, have been qualitatively checked against the model predictions. The cross polar cap electron density profiles of a large majority of the passes more closely conform to the distorted two-cell runs for both polarities of the IMF By component. This test can be generalized to rule out proposed convection patterns based on the presence/absence and position of polar electron density holes.

Maynard, N. C.

Solar wind proton velocity distributions - Comparison of the bi-Maxwellian based 16-moment expansion with observations

The purpose of this paper is to study the possible types of velocity distributions that can be obtained from the bi-Maxwellian based 16-moment expansion of the distribution function, assuming macroscopic parameter values characteristic of the range of solar wind conditions. While previous studies also took heat flow into account, the theoretical expansions for f and the definitions of the physical moments adopted in these studies were different from those used in this paper. The choice of the 16-moment expansion and corresponding moment definitions was motivated by the fact that this is the correct generalization of the widely-used Maxwellian-based 13-moment expansion to the case where the zeroth-order distribution is a bi-Maxwellian. It is found that most of the features characteristic of solar wind proton distributions can be reproduced with the 16-moment distribution, including the appearance of secondary peaks. It is also shown how each of the physically significant velocity moments affects the shape of the distribution function.

Demars, H. G.