Engineering Papers⌕ Search

Engineering topics

Omidi, N.

Publications and source records attributed to Omidi, N..

At least 37 records · Page 2

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

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

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

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

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

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

Magnetic structure of slow shocks and the associated ion dissipation

The kinetic structure of slow magnetosonic shocks obtained from 1- and 2-D hybrid (fluid electrons, particle ions) simulations is discussed. The magnetic structure of the shock varies considerably with the Mach number. At low Mach numbers, the shock may be laminar with no associated waves, whereas, at high Mach numbers the shock may be unsteady with Alfvenic turbulence in the upstream and the downstream regions. The nature of ion dissipation also varies with the Mach number. At weak laminar shocks, ion dissipation consists of adiabatic heating and reflection of a small number of ions back into the upstream. At stronger shocks, particle scattering by Alfven/ion-cyclotron waves in the upstream contributes to the heating which continues well into the downstream region.

Omidi, N.↗

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

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

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

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

Hybrid simulations of the effects of interstellar pickup hydrogen on the solar wind termination shock

Hybrid (kinetic ions/fluid electrons) plasma simulations are used to study the effects of a population of energetic interstellar pickup hydrogen ions on the solar wind termination shock. The pickup hydrogen is treated as a second ion species in the simulations, and thus the effects of the pick-ups on the shock, as well as the effects of the shock on the pickups, are treated in a fully self-consistent manner. For quasi-perpendicular shocks with 10-20 percent pickup hydrogen the pickup ions manifest themselves in a small foot ahead of the shock ramp caused by pickup ion reflection. For oblique shocks with smaller angles between the field and the shock normal, a large fraction of the pickup ions are reflected and move back upstream where they excite large amplitude magnetosonic waves which steepen into shocklets. These backstreaming pickup ions may provide advance warning of a spacecraft encounter with the termination shock.

Liewer, P. C.↗

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

Electromagnetic ion/ion cyclotron instability - Theory and simulations

Linear theory and 1D and 2D hybrid simulations are employed to study electromagnetic ion/ion cyclotron (EMIIC) instability driven by the relative streaming of two field-aligned ion beams. The characteristics of the instability are studied as a function of beam density, propagation angle, electron-ion temperature ratios, and ion beta. When the propagation angle is near 90 deg the EMIIC instability has the characteristics of an electrostatic instability, while at smaller angles electromagnetic effects play a significant role as does strong beam coupling. The 2D simulations point to a narrowing of the wave spectrum and accompanying coherent effects during the linear growth stage of development. The EMIIC instability is an important effect where ion beta is low such as in the plasma-sheet boundary layer and upstream of slow shocks in the magnetotail.

Winske, D.↗

Kinetic structure of slow shocks - Effects of the electromagnetic ion/ion cyclotron instability

The structure of slow magnetosonic shocks in the low beta regime is analyzed with attention given to ion heating and the effects of waves upstream of the electromagnetic ion/ion cyclotron (EMIIC) instability. Shock formation is assessed by means of three methods - a relaxation method and two based on dynamic flow interactions - to determine the effects of initialization and boundary conditions on the formation. Good solutions are found with the piston method and the similar flow-flow method in which the plasma is injected from two boundaries to form two slow shocks. Plasma parameters and shock normal angle are found to be the key variables dictating the structure of the magnetosonic shocks. Four unique classes of resultant shock structures are described in which classical, steady, or nonsteady behavior is found. The analysis also yields insight into the relationship between EMIIC instability and ion dissipation.

Omidi, N.↗

Hybrid simulations of intermediate shocks - Coplanar and noncoplanar solutions

A hybrid code is used to investigate the kinetic structure and stability of subfast intermediate shocks (IS) formed dynamically by the interaction between a flowing plasma and a stationary piston. Results of the kinetic simulation of noncoplanar ISs are compared with predictions of the MHD theory. The relevance of the results of the study to observations of the magnetopause is discussed.

Karimabadi, H.↗

Rotational discontinuities in anisotropic plasmas

The kinetic structure of rotational discontinuities (RDs) in anisotropic plasmas with T-perpendicular/T-parallel is greater than 1 is investigated by using a 1D electromagnetic hybrid code. To form the RD, a new approach is used where the plasma is injected from one boundary and reflected from the other, resulting in the generation of a traveling fast shock and an RD. The results show that across the RD both the magnetic field strength and direction, as well as the plasma density change. Given that such a change can also be associated with an intermediate shock, the Rankine-Hugoniot relations are used to confirm that the observed structures are indeed RDs. It is found that the thickness of RDs is a few ion inertial lengths and is independent of the rotation angle. Also, the preferred sense of rotation is in the electron sense; however, RDs with a rotation angle larger than 180 deg are found to be unstable, changing their rotation to a stable ion sense.

Omidi, N.↗