Engineering topics

Karimabadi, H.

Publications and source records attributed to Karimabadi, H..

At least 19 records

NASA NTRS2016

Turbulence Heating ObserveR: - Satellite Mission Proposal

The Universe is permeated by hot, turbulent, magnetized plasmas. Turbulent plasma is a major constituent of active galactic nuclei, supernova remnants, the intergalactic and interstellar medium, the solar corona, the solar wind and the Earths magnetosphere, just to mention a few examples. Energy dissipation of turbulent fluctuations plays a key role in plasma heating and energization, yet we still do not understand the underlying physical mechanisms involved. THOR is a mission designed to answer the questions of how turbulent plasma is heated and particles accelerated, how the dissipated energy is partitioned and how dissipation operates in different regimes of turbulence. THOR is a single-spacecraft mission with an orbit tuned to maximize data return from regions in near-Earth space magnetosheath, shock, foreshock and pristine solar wind featuring different kinds of turbulence. Here we summarize the THOR proposal submitted on 15 January 2015 to the Call for a Medium-size mission opportunity in ESAs Science Programme for a launch in 2025 (M4). THOR has been selected by European Space Agency (ESA) for the study phase.

Vaivads, A.↗

NASA NTRS2014

The Link Between Shocks, Turbulence, and Magnetic Reconnection in Collisionless Plasmas

Global hybrid (electron fluid, kinetic ions) and fully kinetic simulations of the magnetosphere have been used to show surprising interconnection between shocks, turbulence and magnetic reconnection. In particular collisionless shocks with their reflected ions that can get upstream before retransmission can generate previously unforeseen phenomena in the post shocked flows: (i) formation of reconnecting current sheets and magnetic islands with sizes up to tens of ion inertial length. (ii) Generation of large scale low frequency electromagnetic waves that are compressed and amplified as they cross the shock. These 'wavefronts' maintain their integrity for tens of ion cyclotron times but eventually disrupt and dissipate their energy. (iii) Rippling of the shock front, which can in turn lead to formation of fast collimated jets extending to hundreds of ion inertial lengths downstream of the shock. The jets, which have high dynamical pressure, 'stir' the downstream region, creating large scale disturbances such as vortices, sunward flows, and can trigger flux ropes along the magnetopause. This phenomenology closes the loop between shocks, turbulence and magnetic reconnection in ways previously unrealized. These interconnections appear generic for the collisionless plasmas typical of space, and are expected even at planar shocks, although they will also occur at curved shocks as occur at planets or around ejecta.

MAGNETOSPHERE↗

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

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

A new approach to the linear theory of single-species tearing in two-dimensional quasi-neutral sheets

We have developed the linear theory of collisionless ion tearing in a two-dimensional magnetotail equilibrium for a single resonant species. We have solved the normal mode problem for tearing instability by an algorithm that employs particle-in-cell simulation to calculate the orbit integrals in the Maxwell-Vlasov eigenmode equation. The results of our single-species tearing analysis can be applied to ion tearing where electron effects are not included. We have calculated the tearing growth rate as a function of the magnetic field component B(sub n) normal to the current sheet for thick and thin current sheets, and we show that marginal stability occurs when the normal gyrofrequency Omega(sub n) is comparable to the Harris neutral sheet growth rate. A cross-tail B(sub y) component has little effect on the growth rate for B(sub y) approximately = B(sub n). Even in the limit B(sub y) much greater than B(sub n), the mode is strongly stabilized by B(sub n). We report than random pitch angle scattering can overcome the stabilizing effect of B(sub n) and drive the growth rate up toward the Harris neutral sheet (B(sub n) = 0) value when the pitch angle diffusion rate is comparable to Omega(sub n).

Brittnacher, M.↗

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

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 NTRS1995

Physics of intermediate shocks: A review

Intermediate shocks (ISs) lead to a transition from super-Alfvenic to sub-Alfvenic flow and are different from slow and fast shocks in that an IS rotates the component of the magnetic field tangent to the shock plane by 180 deg. Another peculiarity of ISs is that for the same upstream conditions an IS can have two different downstream states. There also exist a second class of ISs which rotate the magnetic field by an angle other than 180 deg. Due to their noncoplanar nature they cannot be time-stationary and are referred to as time-dependent intermediate shocks (TDIS). The existence of ISs has been the subject of much controversy over the years. Early studies questioned the physical reality of ISs. However, the studies of ISs found a new impetus when C.C. Wu showed that ISs do exist and are stable within the resistive MHD framework. In this paper, after a brief historical overview of the subject, we will review the latest developments in the study of ISs. In particular, we will address the questions of stability and structure of ISs and the relationship between ISs and other discontinuities. One of the recent developments has been the finding that ISs can be unsteady, reforming in time. Details of this process will be discussed. Finally, we examine the effect of anisotropy on the resolutions and discuss the relevance of ISs to the observed field rotations at the Earth's magnetopause.

Karimabadi, H.↗

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

On the energy principle and ion tearing in the magnetotail

We re-examine the use of the energy principle as applied to the tearing instability in the magnetotail. We demonstrate that when a magnetic field component normal to the current sheet is present, electron pitch-angle diffusion (PAD) either by micro-turbulence or by chaotic orbits cannot remove the strong stabilization of the tearing mode caused by electron compressibility. We find that our conclusions are in agreement with those of Pellat et al. (1991), who argued on the basis of canonical P(sub y) conservation that the stabilization of the ion tearing mode cannot be removed by the introduction of PAD. Our results are at variance with those of Kuznetsova and Zelenyi (1991), who argued that the application of the energy principle used by Pellat et al. (1991) is incorrect, and that tearing is in fact unstable in the limit of strong PAD. We show that the disagreement between these two studies can be traced to an incorrect orbit evaluation first introduced by Coroniti (1980) and subnsequently used by Kuznetsova and Zelenyi (1991).

Brittnacher, M.↗

NASA NTRS1992

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

NASA NTRS1992

Physics of pitch angle scattering and velocity diffusion. I - Theory

A general theory for the pitch angle scattering and velocity diffusion of particles in the field of a spectrum of waves in a magnetized plasma is presented. The test particle theory is used to analyze the particle motion. The form of diffusion surfaces is examined, and analytical expressions are given for the resonance width and bounce frequency. The resonance widths are found to vary strongly as a function of harmonic number. The resulting diffusion can be quite asymmetric with respect to pitch angle of 90 deg. The conditions for the onset of pitch angle scattering and energy diffusion are explained in detail. Some of the known shortcomings of the standard quasi-linear theory are also addressed, and ways to overcome them are shown. In particular, the often stated quasi-linear gap at 90 deg is found to exist only under very special cases. For instance, oblique wave propagation can easily remove the gap. The conditions for the existence of the gap are described in great detail. A new diffusion equation which takes into account the finite resonance widths is also discussed. The differences between this new theory and the standard resonance broadening theory is explained.

Karimabadi, H.↗

NASA NTRS1992

Anomalous ion mixing within a Kelvin-Helmholtz vortex in a collisionless plasma

Anomalously fast ion mixing is observed in a hybrid code simulation (particle ions and a massless electron fluid) of the Kelvin-Helmholtz instability in a collisionless plasma. While the traditional view predicts that the ion mixing occurs on a time scale longer than the roll-up time T(v) of the Kelvin-Helmholtz vortex by a factor of (velocity shear scale length)/(ion Larmor radius), the simulation results show that the mixing occurs within a time scale comparable to T(v). It is shown that this fast mixing is due to the scattering of ions by unsteady electromagnetic fields.

Terasawa, T.↗

NASA NTRS1991

Two-dimensional simulations of the ion/ion acoustic instability and electrostatic shocks

A newly developed 2D electrostatic code with particle ions and Boltzmann electrons is used to investigate the details of the ion/ion acoustic instability and the structure of electrostatic shocks. The simulation results show that, for the parameters relevant to the plasma sheet boundary layer, the saturation mechanism of the ion/ion acoustic instability is ion trapping. It is also shown that the 2D structure of electrostatic shocks is considerably different from that suggested by previous 1D simulations. The main reason for this difference is the presence of shock reflected ions, which through the ion/ion acoustic instability lead to the generation of large amplitude waves in the upstream region propagating obliquely to the shock normal. These waves play an important role in the shock dissipation process.

Karimabadi, H.↗

NASA NTRS1991

Collisionless reconnection in two-dimensional magnetotail equilibria

A two-dimensional particle simulation model based on the Darwin approximation to Maxwell's equations for studying collisionless reconnection in the magnetotail has been developed. Simulations of the pure ion tearing mode in a thin current sheet with normal B(z) field component demonstrate that in this limit this mode grows more slowly than expected based on previous analytic estimates. The saturation level of the tearing instability greatly surpasses estimates based on a simple trapping argument. The effect of the normal field component on the evolution of the tearing instability is considered. It is found that a normal field of even a few percent on axis strongly inhibits the growth of the instability.

Pritchett, P. L.↗

NASA NTRS1990

Particle orbits in two-dimensional equilibrium models for the magnetotail

Assuming that there exist an equilibrium state for the magnetotail, particle orbits are investigated in two-dimensional kinetic equilibrium models for the magnetotail. Particle orbits in the equilibrium field are compared with those calculated earlier with one-dimensional models, where the main component of the magnetic field (Bx) was approximated as either a hyperbolic tangent or a linear function of z with the normal field (Bz) assumed to be a constant. It was found that the particle orbits calculated with the two types of models are significantly different, mainly due to the neglect of the variation of Bx with x in the one-dimensional fields.

Karimabadi, H.↗

NASA NTRS1990

Particle acceleration by a wave in a strong magnetic field - Regular and stochastic motion

A general theory for the acceleration of a charged particle by a coherent wave of arbitrary polarization, propagation angle, and phase velocity in the presence of a uniform and strong magnetic field is presented. It is shown that the Hamiltonian surfaces are topologically open for waves with parallel phase velocity, equal to or larger than the speed of light. The trapping width is found to be a strong function of the index of refraction (N), and for N = 1 the trapping width increases as a function of the harmonic number. A novel mechanism for coherently accelerating particles to unlimited energies is presented.

Karimabadi, H.↗