Engineering topics

Krauss-Varban, D.

Publications and source records attributed to Krauss-Varban, D..

At least 19 records

NASA NTRS1995

Properties of ultra low frequency upstream waves at Venus and Saturn: A comparison

The upstream regions of all planets, except Pluto, have been investigated, using in situ spacecraft measurements and a variety of analysis techniques. The detailed studies at Earth indicate that these waves are generated locally in the magnetically connected solar wind by the interaction with ions backstreaming from the shock. However, since the properties of the solar wind vary with heliocentric distance and since properties of planetary shocks depend on plasma beta, interplanetary magnetic field (IMF) spiral angle and Mach number, the amount of heating, acceleration efficiencies, etc. significantly change with heliocentric distance. In turn the waves seen at each planet propagate not in the same but different (physical) propagation modes. In this paper we compare the ULF wave observations at an outer and an inner planet. We use the results of the ratio, quantites easily derivable with sufficient accuracy at each planet. We use the full electromagnetic dispersion relation for comparison with theoretical predictions.

Orlowski, D. S.↗

NASA NTRS1995

Kinetic structure of intermediate shocks: Implications for the magnetopause

A general study of the structure and stability of intermediate shocks (IS) in an isotropic plasma is presented using a hybrid as well as a resistive Hall MHD code. Special emphasis is put on the question of whether the rotational layers observed at the magnetopause can be intermediate shocks. The shocks are formed dynamically by the interaction between a flowing plasma and a stationary piston. Coplanar ISs (both strong and weak) are found to be stable in a collisionless plasma. The existence of slow shocks in a high beta plasma is also established for the first time. Noncoplanar ISs are found to be time-dependent, evolving toward a rotational discontinuity (RD) after some characteristic time tau which can be quite long (1000 Omega(exp -1), where Omega is the ion gyrofrequency). The value tau is larger the closer the rotation angle is to 180 deg. Rotations larger than 180 deg are found to be unstable, decaying into a state of minimum shear (i.e., rotation angle less than 180 deg). There are various length scales associated with an IS in the kinetic regime. The shortest scale is found to be the length scale over which rotation of the transverse component of the magnetic field takes place. This scale can have a half width as small as one ion inertial length (c/omega(sub p)) for electron sense rotations and 3c/omega(sub p) for ion sense rotations, for an upstream ion beta of unity. Both of these scales are consistent with the observed thickness at the magnetopause and identical to the corresponding RD scales. A detailed study of the mode conversion of the Alfven ion cyclotron waves (A/IC) waves across both slow and intermediate shocks and the resulting downstream wave spectrum are presented. The possibility that the large number of relfected ions observed at the magnetopause may be due to the presence of strong ISs is considered. The identification of strong ISs and their distinction from RDs should be possible in observations due to significant differences that exist between jump conditions and overall structure of the two discontinuities. The jumps in the plasma parameters across a weak IS are typically small. This together with the fact that the weak ISs and RDs have very similar thickness and other overall properties makes the distinction between weak ISs and RDs in the observations largely inconsequential. However, at large noncoplanarity angles the weak IS approaches the RD limit in a relatively short time (approximately less than 100 Omega (exp -1)). Thus, magnetopause rotations with large noncoplanarity angles are most likely either RDs or strong ISs. Finally, direct comparisons between fluid (resistive Hall MHD) and kinetic simulations show that fluid theory is not applicable to study of ISs in a collisionless plasma.

Karimabadi, H.↗

NASA NTRS1995

Kinetic structure of rotational discontinuities: Implications for the magnetopause

Magnetic field rotations in the high ion beta magnetosheath that are part of the magnetopause structure are expected to have only a small normal component. We have studied the properties of rotational discontinuities (RDs) under these conditions, viewed as the limit of weak intermediate shocks (ISs), by performing hybrid simulations with a reflecting wall boundary condition (piston method). With this dynamic formation, the sense and size of rotation are not arbitrarily predetermined, but rather evolve from the given upstream (magnetosheath) and downstream (magnetospheric) boundary conditions, similar to what takes place at the magnetopause. This work focuses on several aspects: the observed minimum shear of RDs, their width, their internal signature, and their relation to ISs in isotropic plasmas. Our simulation results are in agreement with the minimum shear observations, that is, the RDs choose the sense of rotation that corresponds to the minimum angle between the upstream and downstream field vector. The RDs are stable, with a unique scale size. Typical gradient scale half widths are one to four ion inertial lengths with a total width up to ten times of that, in agreement with magnetopause observations. We develop a generalized fluid theory of RDs and discuss the characteristic internal signatures of the rotational layer, comparing the kinetic simulation results to predictions from the generalized fluid theory. The results show that ion inertia, anisotropic pressure, finite Larmor radius effects, nonzero ion heat flux, and reflected ions all contribute to the signatures of RDs on kinetic scales. The RDs may have upstream or downstream wave trains, which become weak for high ion beta and small normal components of the magnetic field. We explain the presence and direction of wave trains in terms of the kinetic properties of the Alfven/ion-cyclotron mode. Away from the RD limit there is a smooth transition to weak intermediate shocks, which have small jumps close to expected Rankine-Hugoniot values. Apart from that, there are few kinetic plasma signatures that distinguish RDs from their neighboring ISs. However, noncoplanar ISs evolve in time into thin RDs. Using the properties of RDs and ISs, we make specific suggestions how these discontinuities can be distinguished observationally in the case of an isotropic plasma.

Krauss-Varban, D.↗

NASA NTRS1995

Electron distributions in quasi-perpendicular shocks: Simulations

One of the outstanding problems of collisionless shock physics is the microphysical description of electron thermalization and associated wave processes. Observations at the Earth's bow shock have shown various features in the electron velocity distribution such as a loss-cone, temperature anisotropies, and beams before a thermalized, flat-topped distribution is reached downstream. It is generally believed that in addition to the macroscopic fields, which determine the overall mapping of the distribution function, electromagnetic whistler waves and lower hybrid waves are responsible for the pitch-angle scattering and thermalization of the electrons. However, a genuine physical explanation requires analysis of the evolution of the electron distribution in the self-consistent wave fields. Particle simulations are ideally suited to do this, and include processes in which both electrons and ions are involved, such as cross-field streaming instabilities. We show results from l-D and 2-D simulations using an implicit code, which allows time steps large compared to the inverse electron plasma frequency. It is shown that in order to capture the correct propagation characteristics of the waves, a large mass ratio between the protons and electrons is generally necessary. Detailed diagnostics of the waves as well as the evolution of the electron velocity distribution will be presented and will be discussed in the frame work of satellite observations.

Pantellini, F.↗

NASA NTRS1995

Linear and nonlinear properties of the ULF waves driven by ring-beam distribution functions

The problem of the exitation of obliquely propagating magnetosonic waves which can steepen up (also known as shocklets) is considered. Shocklets have been observed upstream of the Earth's bow shock and at comets Giacobini-Zinner and Grigg-Skjellerup. Linear theory as well as two-dimensional (2-D) hybrid (fluid electrons, particle ions) simulations are used to determine the properties of waves generated by ring-beam velocity distributions in great detail. The effects of both proton and oxygen ring-beams are considered. The study of instabilities excited by a proton ring-beam is relevant to the region upstream of the Earth's bow shock, whereas the oxygen ring-beam corresponds to cometary ions picked up by the solar wind. Linear theory has shown that for a ring-beam, four instabilities are found, one on the nonresonant mode, one on the Alfven mode, and two along the magnetosonic/whistler branch. The relative growth rate of these instabilities is a sensitive function of parameters. Although one of the magnetosonic instabilities has maximum growth along the magnetic field, the other has maximum growth in oblique directions. We have studied the competition of these instabilities in the nonlinear regime using 2-D simulations. As in the linear limit, the nonlinear results are a function of beam density and distribution function. By performing the simulations as both initial value and driven systems, we have found that the outcome of the simulations can vary, suggesting that the latter type simulations is needed to address the observations. A general conclusion of the simulation results is that field-aligned beams do not result in the formation of shocklets, whereas ring-beam distributions can.

Killen, K.↗

NASA NTRS1995

Wave and ion evolution downstream of quasi-perpendicular bow shocks

Distribution functions of ions heated in quasi-perpendicular bow shocks have a large perpendicular temperature anisotropy that provides free energy for the growth of Alfven ion cyclotron (AIC) waves and mirror waves. Both types of waves have been observed in the Earth's magnetosheath downstream of quasi-perpendicular shocks. We use a two-dimensional hybrid simulations to give a self-consistent description of the evolution of the wave spectra downstream of quasi-perpendicular shocks. Both mirror and AIC waves are identified in the simulated magnetosheath. They are generated at or near the shock front and convected away from it by the sheath plasma. Near the shock, the waves have a broad spectrum, but downstream of the shock, shorter-wavelength modes are heavily damped and only longer-wavelength modes persist. The characteristics of these surviving modes can be predicted with reasonable accuracy by linear kinetic theory appropriate for downstream conditions. We also follow the evolution of the ion distribution function. The shocked ions that provide the free energy for wave growth have a two-component distribution function. The halo is initially gyrophase-bunched and extremely anisotropic. Within a relatively short distance downstream of the shock (of the order of 10 ion inertial lengths), wave-particle interactions remove these features from the halo and reduce the anisotropy of the distribution to near-threshold levels for the mirror and AIC instabilities. A similar evolution has been observed for ions at the Earth's bow shock.

Mckean, M. E.↗

NASA NTRS1995

Two-dimensional structure of slow shocks

The two-dimensional (2-D) structure of switch-off slow magnetosonic shocks is investigated using an electromagnetic hybrid (fluid electrons, kinetic ions) code. It is shown that the basic physical processes occurring at 1-D slow shocks are also operative in 2-D. Specifically, the interaction between the upstream ions and those streaming away from the shock results in the excitation of Alfven/ions-cyclotron (AIC) waves. Depending on the plasma parameters, these waves may either stay in the upstream or convect back into the shock resulting in a non-steady shock behavior which prevents the formation of a trailing wave train. Despite this similarity, some slow shocks which are steady in 1-D are found to be non-steady in 2-D. Fourier analysis of the waves downstream of non-steady shocks identifies them as AIC, demonstrating that the waves remain on the same branch as they convect from upstream into the downstream region.

Omidi, N.↗

NASA NTRS1995

Propagation and damping of broadband upstream whistlers

Previous studies indicated that damping rates of upstream whistlers strongly depend on the details of the electron distribution function. Moreover, detailed analysis of Doppler-shift and whistler dispersion relation indicated that upstream whistlers propagate obliquely in a broad band. In this paper we present results of a kinetic calculation of damping lengths of wide-band whistlers using the sum of 7-drifting bi-Maxwellian electron distributions as a best fit to the International Sun Earth Explorer (ISEE) 1 electron data. For 2 cases, when upstream whistlers are observed, convective damping lengths derived from ISEE magnetic field and ephemeris data are compared with theoretical results. We find that the calculated convective damping lengths are consistent with the data and that upstream whistlers remain marginally stable. We also show that the slope of plasma frame spectra of upstream whistlers, obtained by direct fitting of the observed spectra is between 5 and 7 with a sharp lower frequency cutoff corresponding to a wavelength of about one ion inertial length. When the solar wind velocity is directed largely along the wave normal of the upstream whistlers the polariztion of the right hand waves becomes reversed and low frequencies are switched to high resulting in a peaked spectrum with a strong high frequency cutoff. The overall spectral, wave and particle characteristics, proximity to the shock as well as propagation and damping properties indicate that these waves cannot be generated locally. Instead the observed upstream whistlers arise in the shock ramp most likely by a variety of cross-field drift and/or anisotropy driven instabilities.

Orlowski, D. S.↗

NASA NTRS1995

Experimental studies of the properties of 'simulated' upstream turbulence using a statistical multipoint method

In this report we present a different approach to the multipoint measurement of magnetic fields and plasma. This is called the multi-spacecraft ensemble technique (MET), essentially free of process restrictions, such as linearity and stationarity. We comprehensively discuss the other conditions and limitations intrinsic to this statistical method. We also show the results of the application of the ensemble method to the synthetic data obtained from a hybrid simulation in the region upstream of a quasi-parallel shock. The important implications of the above approach for the CLUSTER mission are discussed.

Orlowski, D. S.↗

NASA NTRS1995

Waves associated with quasi-parallel shocks: Generation, mode conversion and implications

Ions that are energized at quasi-parallel collisionless shocks and move back upstream generate low-frequency waves, largely on the fast/magnetosonic branch. At sufficient Mach number, the waves are convected back into the shock, lead to shock re-formation, and are mode-converted into downstream (magnetosheath) Alfvenic turbulence. Other waves are generated more locally at the interface of the incoming solar wind and the partially thermalized plasma. This paper reviews how recent simulation studies of collisionless shocks in conjunction with linear kinetic theory and proper wave diagnostics have aided in our understanding of the upstream and magnetosheath waves.

Krauss-Varban, D.↗

NASA NTRS1995

Wave and particle evolution downstream of quasi-perpendicular shocks

Distributions of ions heated in quasi-perpendicular bow shocks have large perpendicular temperature anisotropies that provide free energy for the growth of Alfven ion cyclotron (AIC) and mirror waves. These modes are often obsreved in the Earth's magnetosheath. Using two-dimensional hybrid simulations, we show that these waves are produced near the shock front and convected downstream rather than being produced locally downstream. The wave activity reduces the proton anisotropy to magnetosheath levels within a few tens of gyroradii of the shock but takes significantly longer to reduce the anisotropy of He(++) ions. The waves are primarily driven by proton anisotropy and the dynamics of the helium ions is controlled by the proton waves. Downstream of high Mach number shocks, mirror waves compete effectively with AIC waves. Downstream of low Mach number shocks, AIC waves dominate.

Mckean, M. E.↗

NASA NTRS1994

Low-freqency instabilities and the resulting velocity distributions of pickup ions at comet Halley

The interaction between the solar wind and newborn cometary ions is studied using a new analytical theory as well as one- and two-dimensional hybrid simulations. Using the observed parameters upstream of the comet Halley, a detailed study of wave excitation and the resulting particle distributions is presentd. Linear theory as well as simulations show that a variety of modes such as the fast magnetosonic mode, high frequency whistlers and obliquely propagating Alfven ion cyclotron waves can be excited. However, parallel propagating waves are found to be dominant in the wave spectrum and to control the scattering of the pickup ions. Several features of the observed distributions of pickup protons are explained. In particular, it is shown that the observed asymmetric pitch distribution for the pickup protons is due to the small saturation amplitude of the waves for the given parameters. Water group associated waves can lead to energy diffusion and further pitch angle scattering of protons. This effect is most likely to be important in the vicinity of the bow shock of comet Halley where the density of water group ions becomes comparable to that of protons. It is shown that the observed increase in the radius of the proton velocity shell just outside the bow shock can be due to water group waves. The nearby isotropic proton pitch angle distribution observed by Neugebauer et al. (1989) just outside the bow shock may, however, be related to the presence of a rotational discontinuity which has been identified in the magnetic field data. Just outside the bow shock, simulations show that parallel propagating water group waves can steepen with attached whistler wave packets. The steepening process at parallel propagation is a transient effect, in an important contrast to the case of steepening at oblique angles. The smaller beam densitites at comet Halley appears to be the main reason not only why waves at comet Halley have smaller amplitudes but also why oblique, steepening magnetosonic waves have not been detected at comet Halley, whereas they have been seen at comet Giacobini-Zinner. Halley, whereas they have been seen at comet Giacobini-Zinner.

Karimabadi, H.↗

NASA NTRS1994

Mode properties of low-frequency waves: Kinetic theory versus Hall-MHD

In fluid theory, the ordering of low-frequency modes in a homogeneous plasma is based on the phase velocity, since modes do not intersect each other in dispersion diagrams as a function of wavenumber or other parameters. In linear kinetic theory, modes cross each other. Thus a consistent and useful classification should be based on the physical properties of the modes instead. This paper attempts such a classification by documeting the dispersion and general mode properties of the low-frequency waves (omega much less than (OMEGA(sub ci) OMEGA(sub ce) (exp 1/2)), where OMEGA(sub ci), OMEGA(sub ce) are the cyclotron frequencies of the ions and electrons, respectively) in kinetic theory, and by comparing them to the results of two-fluid theory. Kinetic theory gives a seperate Alfven/ion-cyclotron (A/IC) wave with phase speed Omega/k approximately = v(sub A) cos theta for omega much less than OMEGA(sub ci), where v(sub A) is the Alfven velocity and theta the angle of propagation between wave vector k and background magnetic field B(sub o). For a given wavenumber, the magnetosonic mode is a double-valued solution with a singular point in theta, beta parameter space, where beta is the ratio of thermal pressure to magnetic pressure. It is shown that a branch cut starting at the singular point theta approximately 30 deg, beta approximately 3 and leading to larger beta gives a practical and consitent seperation of this double-valued magnetosonic solution. Selection of this branch cut results in a moderately damped fast/magnetos onic and a heavily damped slow/sound wave. A comprehensive review of the polarization, compressibility and other mode properties is given and shown to be consistent with the selected branch cut. At small wavenumbers, the kinetic mode properties typically start to deviate significantly from their fluid counterparts at beta approximately 0.5. At larger beta, there is no longer a consistent correspondence between the fluid and kinetic modes. Kinetic theory also dictates the use of different mode properties to distinguish between them in observational data. For example, the phase between the density and magnetic field perturbation may become useless at high beta, whereas the direction of the magnetic field perturbations with respect to k and B(sub o) remains a useful characteristic. Two quantities based on this characteristic are suggested and are shown to be useful also to distinguish between the mirror mode and A/IC waves in a plasma with temperature anisotropy.

Krauss-Varban, D.↗

NASA NTRS1994

Electron acceleration at nearly perpendicular collisionless shocks. 3: Downstream distributions

Spacecraft observations at the Earth's bow shock and at interplanetary shocks have established that the largest fluxes of accelerated suprathermal electrons occur in so-called shock spike events immediately downstream of the shock ramp. Previous theoretical efforts have mainly focused on explaining upstream energetic electron beams. Here we investigate the general motion and acceleration of energetic electrons in a curved, nearly perpendicular shock by numerically integrating the orbits of solar wind halo electrons in shock fields generated by a hybrid simulation (core electron fluid and kinetic ions). Close to the angle Theta(sub Bn) = 90 degs between the upstream magnetic field and shock normal, the calculations result in a (perpendicular) temperature increase proportional to the magnetic field ratio and give the highest phase space densities in the overshoot. For a steep distribution, the temperature change can correspond to an enhancement of the distribution by several orders of magnitude. These results are in agreement with predictions from adiabatic mapping. With smaller angles Theta(sub Bn), the overshoot and downstream densities fall off quickly, because the adiabatic energy gain is less and fewer electrons transmit. The shock curvature also leads to an accumulation of electrons close to 90 degs. Without pitch angle scattering, energization is only significant within a few (approximately 5 to 10 degs) degrees of the point of tangency. However, shock spike events appear to be observed more easily and farther away from 90 degs. Given that over a region of several degrees around 90 degs the theory gives enhancements of up to approximately 4 orders of magnitude, such electrons could in principle account for the typically observed enhancements of 1 to 2 orders of magnitude, if they were distributed over Theta(sub Bn). To test the idea that scattering could efficiently redistribute the energetic electrons, we have conducted test particle simulations in which artificial pitch angle scattering is included. We find that such a process can indeed be very efficient and can explain observations of shock spike events far away from Theta(sub Bn) is approximately 90 degs. The scattering naturally leads to much higher phase space densities at smaller Theta(sub Bn) than what a local one-dimensional mapping would predict and thus can account for an observed discrepancy with adiabatic theory stated in the literature.

Krauss-Varban, D.↗

NASA NTRS1994

Growing 'Alfvenic' modes in the upstream region of Saturn

Recent studies of low-frequency electromagnetic waves upstream of the Saturn bow shock have shown that these waves, in contrast to those at Earth, are observed not in one, but in at least two, distinct frequency bands. The results of wave mode identification based on the Hall-magnetohydrodynamic (MHD) model of plasma and observed wave polarization suggested that these waves propagate in the high beta intermediate mode. However, the underlying instability was not unambiguously determined. In the present paper we use the full electromagnetic dispersion relation derived from linear Vlasov theory in order to examine which of the plasma modes, with observed properties, are unstable in an isotropic Maxwellian plasma in the presence of backstreaming proton beams consistent with Voyager 2 observations at Saturn. As a result we find that the unstable 'Alfvenic' beam mode, as well as resonant and non-resonant fast magnetosonic modes have properties consistent with the data. Moreover, we find that in contrast to the Earth's upstream waves, at Saturn no 'kinetic' normal mode can account for the observed magnetic polarization.

Orlowski, D. S.↗

NASA NTRS1994

A test of the Hall-MHD model: Application to low-frequency upstream waves at Venus

Early studies suggested that in the range of parameter space where the wave angular frequency is less than the proton gyrofrequency and the plasma beta, the ratio of the thermal to magnetic pressure, is less than 1 magnetohydrodynamics provides an adequate description of the propagating modes in a plasma. However, recently, Lacombe et al. (1992) have reported significant differences between basic wave characteristics of the specific propagation modes derived from linear Vlasov and Hall-magnetohydrodynamic (MHD) theories even when the waves are only weakly damped. In this paper we compare the magnetic polarization and normalization magnetic compression ratio of ultra low frequency (ULF) upstream waves at Venus with magnetic polarization and normalized magnetic compression ratio derived from both theories. We find that while the 'kinetic' approach gives magnetic polarization and normalized magnetic compression ratio consistent with the data in the analyzed range of beta (0.5 less than beta less than 5) for the fast magnetosonic mode, the same wave characteristics derived from the Hall-MHD model strongly depend on beta and are consistent with the data only at low beta for the fast mode and at high beta for the intermediate mode.

Orlowski, D. S.↗

NASA NTRS1993

Propagation characteristics of waves upstream and downstream of quasi-parallel shocks

The propagation characteristics of waves upstream and downstream of quasi-parallel shocks are investigated by using 2D hybrid simulations. At low Alfven Mach numbers, M(A) below about 2, the shock is initially associated with upstream phase-standing whistlers. At later times, backstreaming ions excite longer-wavelength whistlers via the right-hand resonant ion/ion instability. These waves propagate along the magnetic field at a group velocity no smaller than the upstream flow speed, so that the waves remain in the upstream region. At higher MA (above about 3), these waves are convected back into the shock, causing its reformation and downstream perturbations. Shock transmitted waves mode-convert into Alfven/ion-cyclotron waves which have a wave vector along the shock normal (pointing upstream) and convect downstream. The 2D simulation results confirm our earlier suggestion that the upstream waves should be field aligned, and that their convection into the downstream is associated with linear mode conversion into the Alfven/ion-cyclotron branch.

Krauss-Varban, D.↗

NASA NTRS1993

Structure and length scales of rotational discontinuities

Results are presented of a series of 1D hybrid (kinetic ions and fluid electrons) simulations, using a nondynamic method to form rotational discontinuities (RDs) at an angle of 60 deg between the normal direction and the upstream magnetic field. Ion kinetic effects are found to be important with regard to the stability of the current layer, its thickness and scaling, and its dependence on initial conditions. A detailed examination of the downstream wave train and gradient scale size showed that RDs with alpha1 = +/-180 deg are stable within the framework of the simulation. With proper initialization, electron-sense rotations of -270 deg were found that do not disintegrate immediately but rather appear to consist of a superposition of a solitary wave with -360-deg rotation and a 90-deg ion-sense rotation.

Krauss-Varban, D.↗