Search NASA⌕ Search

Engineering topics

Wang, Xueyi

Publications and source records attributed to Wang, Xueyi.

Magnetic Reconnection Inside Solar Wind Rotational Discontinuity During Its Interaction With the Quasi-Perpendicular Bow Shock and Magnetosheath

Here, using a three-dimensional global hybrid simulation, we investigate the formation and evolution of ion-scale magnetic reconnection inside an interplanetary rotational discontinuity (RD) owing to its interaction with the quasi-perpendicular (Q-⟂) bow shock and the magnetosphere. The interplanetary magnetic field (IMF) is initially predominantly northward, while it changes to purely southward across the RD. A significantly thinned RD current layer with a width ~ ion skin depth and normal magnetic field B n ≈0 is formed by a shock compression process as the RD interacts with the Q-⟂ shock. Magnetic reconnection thereupon takes place inside the thin RD current layer, where Hall magnetic and electric fields, reconnection electric field, and high-speed ion outflow jets are identified. Simultaneously, flux ropes form with an extension of a few ion inertial lengths. As the RD is transmitted into the magnetosheath, multiple reconnection sites lead to the formation of longer flux ropes. Moreover, magnetosheath reconnection is also found at new reconnection sites inside the RD. No reconnection is found in the Q-⟂ shock alone outside the RD. The flux ropes propagate poleward and tailward in the magnetosheath. In addition, magnetopause reconnection takes place under the southward IMF on the sunward side of the RD after the RD passes through the magnetopause.

79 ASTRONOMY AND ASTROPHYSICS↗

Global Hybrid Simulations of Interaction Between Interplanetary Rotational Discontinuity and Bow Shock/Magnetosphere: Can Ion-Scale Magnetic Reconnection be Driven by Rotational Discontinuity Downstream of Quasi-Parallel Shock?

Ion-scale magnetic reconnection has been observed downstream of the terrestrial quasi-parallel (Q-∥) shock. Whether it is driven by interplanetary discontinuities or turbulent Q-∥ shock, however, is unclear. Using three-dimensional global hybrid simulation, we investigate the generation of magnetic reconnection downstream of the Q-∥ shock, while an interplanetary rotational discontinuity (RD) is launched to the bow shock. Cases with various solar wind Alfvén Mach numbers, M A = 3.0 to 8, and propagation directions of the RD are presented. The propagation direction n is assumed to be in the GSE xz plane and pointing earthward, with n = (-sin(θ 12 /2),0,-cos(θ 12 /2)), where θ 12 is the angle between the upstream (B1) and downstream (B2) magnetic fields across the transmitted RD. It is found that magnetic reconnection occurs inside the RD downstream of the Q-∥ shock, forming flux ropes extending along the dawn-dusk direction about tens of ion inertial lengths. Large-amplitude low-frequency waves originated from the Q-∥ shock lead to the bending and squeezing of the field lines around the RD, which play an important role in triggering reconnection inside the RD. As the RD impacts the dayside magnetopause, magnetopause reconnection takes place between the field lines behind the RD and geomagnetic field lines. Nevertheless, no reconnection is found downstream of the Q-∥ shock itself or outside the RD in the magnetosheath. The existent and structure of reconnection in the magnetosheath are found to strongly depend on the parameters M A and n. Our simulation shows that ion-scale magnetic reconnection is driven by an external driver in the form of the compression of an RD around the bow shock and in the magnetosheath, rather than caused by the turbulent Q-∥ shock alone.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetopause Reconnection as Influenced by the Dipole Tilt Under Southward IMF Conditions: Hybrid Simulation and MMS Observation

Abstract Using a three‐dimensional (3‐D) global‐scale hybrid code, the Magnetospheric Multiscale (MMS) reconnection event around 02:13 UT on 18 November 2015, highlighted in the Geospace Environment Modeling (GEM) Dayside Kinetic Challenge, is simulated, in which the interplanetary magnetic field (IMF) points southward and the geomagnetic field has a −27° dipole tilt angle. Strong southward plasma jets are found near the magnetopause as a result of the dayside reconnection. Our results indicate that the subsolar magnetopause reconnection X line shifts from the subsolar point toward the Northern Hemisphere due to the effect of the tilted geomagnetic dipole angle, consistent with the MMS observation. Subsequently, the reconnection X lines or sites and reconnection flux ropes above the equator propagate northward along the magnetopause. The formation and global distribution of the X lines and the structure of the magnetopause reconnection are investigated in detail with the simulation. Mirror mode waves are also found in the middle of the magnetosheath downstream of the quasi‐perpendicular shock where the plasma properties are consistent with the mirror instability condition. As a special outcome of the GEM challenge event, the spatial and temporal variations in reconnection, the electromagnetic power spectra, and the associated D‐shaped ion velocity distributions in the simulated reconnection event are compared with the MMS observation.

Guo, Zhifang↗

Generation of kinetic Alfvén waves in dayside magnetopause reconnection: A 3-D global-scale hybrid simulation

In this paper, we perform a three-dimensional (3-D) global-scale hybrid simulation to investigate the generation of kinetic Alfvén waves (KAWs) in reconnection at the dayside magnetopause. In the simulation, the magnetopause reconnection takes place due to a self-consistent interaction between the solar wind and Earth's magnetosphere, in which the solar wind is assumed to carry a southward interplanetary magnetic field. As multiple X-line reconnection occurs, shear Alfvénic perturbations are generated from reconnection, which propagate toward polar regions along the magnetic field lines. Our detailed analysis of the electromagnetic characteristics reveals the existence of kinetic Alfvén waves on the ion kinetic scales in these perturbations, in which the parallel electric field is generated. The 3-D structures and the Poynting fluxes in the shear Alfvén waves/KAWs in the global magnetopause are presented. It is shown that these KAWs carry electromagnetic energy into the cusp, which may lead to aurora brightening.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetosheath Reconnection Before Magnetopause Reconnection Driven by Interplanetary Tangential Discontinuity: A Three‐Dimensional Global Hybrid Simulation With Oblique Interplanetary Magnetic Field

Abstract Terrestrial dayside dynamics associated with a southward turning, oblique interplanetary magnetic field (IMF) carried by an interplanetary tangential discontinuity (TD) is investigated by performing a three‐dimensional global‐scale hybrid simulation systematically for cases in which the incoming solar wind TD possesses various magnetic field rotation angles ΔΦ = 90° to 180° and half widths w = 2 d i 0 to w = 30 d i 0 , where d i 0 is the ion inertial length in the solar wind. Overall, the TD is compressed while being transmitted into the magnetosheath, with different compression processes downstream of the quasi‐parallel (Q‐||) and quasi‐perpendicular (Q‐ ⊥ ) shocks. It is found that magnetosheath reconnection may take place downstream of both the Q‐|| and Q‐ ⊥ shocks due to interaction of the directional TD with the bow shock and magnetopause, but the existence of magnetosheath reconnection depends on w and ΔΦ . Magnetosheath flux ropes are formed through three‐dimensional patchy reconnection in the thinned current sheet, with a longer rope lengths under a larger ΔΦ . There exists a dawn‐dusk asymmetry in the spatial extent of the flux ropes, which becomes more significant as ΔΦ decreases. When ΔΦ decreases to 90° , no reconnection flux ropes are found. Magnetopause reconnection is initiated when the magnetic fluxes with a southward turning IMF (on the sunward side) reach the magnetopause, and the magnetopause flux ropes can be mixed with the magnetosheath ones. Our simulation demonstrates that the effects of a southward turning of the IMF may not be a simple field direction change that leads to reconnection only at the magnetopause.

Guo, Zhifang↗

Particle-in-cell simulations of magnetically driven reconnection using laser-powered capacitor coils

In this study, we propose an experimental scheme to fulfill magnetically driven reconnections. Here, two laser beams are focused on a capacitor-coil target and then strong currents are wired in two parallel circular coils. Magnetic reconnection occurs between the two magnetic bubbles created by the currents in the two parallel circular coils. A two-dimensional particle-in-cell simulation model in the cylindrical coordinate is used to investigate such a process, and the simulations are performed in the ( r , z ) plane. The results show that with the increase of the currents in the two coils, the associated magnetic bubbles expand and a current sheet is formed between the two bubbles. Magnetic reconnection occurs when the current sheet is sufficiently thin. A quadrupole structure of the magnetic field in the θ direction ( B θ ) is generated in the diffusion region and a strong electron current along the r direction ( J e r ) is also formed due to the existence of the high-speed electron flow away from the X line in the center of the outflow region. Finally, because the X line is a circle along the θ direction, the convergence of the plasma flow around r = 0 will lead to the asymmetry of J e r and B θ between the two outflow regions of magnetic reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ion acceleration and heating by kinetic Alfvén waves associated with magnetic reconnection

In a previous study on the generation and signatures of kinetic Alfv en waves (KAWs) associated with magnetic reconnection in a current sheet revealed that KAWs are a common feature during reconnection [Liang et al. J. Geophys. Res.: Space Phys. 121, 6526 (2016)]. In this paper, ion acceleration and heating by the KAWs generated during magnetic reconnection are investigated with a three-dimensional (3-D) hybrid model. It is found that in the outflow region, a fraction of inflow ions are accelerated by the KAWs generated in the leading bulge region of reconnection, and their parallel velocities gradually increase up to slightly super-Alfv enic. As a result of waveparticle interactions, an accelerated ion beam forms in the direction of the anti-parallel magnetic field, in addition to the core ion population, leading to the development of non-Maxwellian velocity distributions, which include a trapped population with parallel velocities consistent with the wave speed. We then heat ions in both parallel and perpendicular directions. In the parallel direction, the heating results from nonlinear Landau resonance of trapped ions. In the perpendicular direction, however, evidence of stochastic heating by the KAWs is found during the acceleration stage, with an increase of magnetic moment μ. The coherence in the T $\perp$ ion temperature and the perpendicular electric and magnetic fields of KAWs also provides evidence for perpendicular heating by KAWs. The parallel and perpendicular heating of the accelerated beam occur simultaneously, leading to the development of temperature anisotropy with the perpendicular temperature T $\perp$ >T $\parallel$ temperature. The heating rate agrees with the damping rate of the KAWs, and the heating is dominated by the accelerated ion beam. In the later stage, with the increase of the fraction of the accelerated ions, interaction between the accelerated beam and the core population also contributes to the ion heating, ultimately leading to overlap of the beams and an overall anisotropy with T $\perp$ >T $\parallel$ .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Foreshock wave interaction with the magnetopause: Signatures of mode conversion

Our previous hybrid simulation under a radial interplanetary magnetic field (IMF) and a supercritical solar wind Mach number has shown that foreshock compressional waves originated from the quasi–parallel (Q–∥) shock are mode converted to kinetic Alfvén waves (KAWs) at the Alfvén resonance surface of the subsolar magnetopause. In this paper, three–dimensional global dayside mode conversion is investigated for cases under various solar wind conditions using the global hybrid model. The global patterns and propagations of KAWs are distinguished and presented. Under a near–critical Mach number ( M A =3), KAW structures due to mode conversion exhibit a feature of broader excitation regions in the magnetopause boundary layer (MPBL) compared to supercritical Mach number ( M A =5) shocks. For cases with an oblique IMF with supercritical Mach numbers ( M A =5), the amplitude of magnetosheath compressional waves is larger at the quasi–parallel shock (Q–∥) than at the quasi–perpendicular (Q–⊥) shock. Downstream of the Q–∥ shock, there is a general trend that the perturbations of density ( N ) and magnetic field ( B ) change from predominantly in–phase in the magnetosheath to antiphase near the MPBL. While downstream of the Q–⊥ shock, an antiphase relation between N and B is dominant throughout the magnetosheath and magnetopause except near the shock transition. The compressional drivers are found to reach an extended region of the magnetopause due to the combined effects of wave propagation in the plasma frame and flow convection, leading to a broad region of mode conversion at the magnetopause. Subsequently, the resulting KAWs can be carried to the regions downstream of the Q–⊥ shock owing to the flow convection at the magnetopause. Furthermore the KAWs propagate poleward along the geomagnetic field lines and meanwhile are carried tailward by the ambient flows, and they are more intense in the downstream of Q–∥ shocks than downstream of Q–⊥ shocks.

79 ASTRONOMY AND ASTROPHYSICS↗

Interactions of electrons with two lower hybrid waves

The effects of perturbed orbits on the interactions of electrons with two lower hybrid waves, one of which is resonant with electrons at a low phase velocity (v p1 = 3.8V the , where v p1 is the wave phase velocity and V the is the electron thermal speed) while the other is off-resonant at a high phase velocity (v p2 = 5.5V the ), have been studied by using the particle simulation code based on the gyro-kinetic electron and fully-kinetic ion (GeFi) model [Lin et al., Plasma Phys. Controlled Fusion 47, 657 (2005)]. When the amplitude of the off-resonant wave is sufficiently small so that the resonances of these two waves do not overlap, the variation of the resonant wave amplitude is similar to that predicted by O'Neil's theory [O'Neil, Plasma Fluid 8, 12 (1965)]. With the increasing amplitude, the two resonances overlap and large scale chaos emerges. As a result, the damping of the resonant wave can be enhanced, which is due to that the trapped electron orbits are significantly perturbed by the off-resonant wave. The diffusion process gives rise to the enhanced damping. In conclusion, when the overlap is sufficiently large, the damping of the off-resonant wave and the oscillatory behavior of the wave amplitude are observed. In addition, the resonant plateau in the distribution function can be broadened due to the change in the chaotic region boundaries as the electron perturbed orbits are taken into account.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

3D electrostatic gyrokinetic electron and fully kinetic ion simulation of lower-hybrid drift instability of Harris current sheet

The eigenmode stability properties of three-dimensional lower-hybrid-drift-instabilities (LHDI) in a Harris current sheet with a small but finite guide magnetic field have been systematically studied by employing the gyrokinetic electron and fully kinetic ion (GeFi) particle-in-cell (PIC) simulation model with a realistic ion-to-electron mass ratio m i /m e . In contrast to the fully kinetic PIC simulation scheme, the fast electron cyclotron motion and plasma oscillations are systematically removed in the GeFi model, and hence one can employ the realistic m i /m e . The GeFi simulations are benchmarked against and show excellent agreement with both the fully kinetic PIC simulation and the analytical eigenmode theory. Our studies indicate that, for small wavenumbers, ky, along the current direction, the most unstable eigenmodes are peaked at the location where $\vec{k}$• $\vec{B}$ =0, consistent with previous analytical and simulation studies. Here, $\vec{B}$ is the equilibrium magnetic field and $\vec{k}$ is the wavevector perpendicular to the nonuniformity direction. As ky increases, however, the most unstable eigenmodes are found to be peaked at $\vec{k}$ •$\vec{B}$ ≠0. Additionally, the simulation results indicate that varying m i /m e , the current sheet width, and the guide magnetic field can affect the stability of LHDI. Simulations with the varying mass ratio confirm the lower hybrid frequency and wave number scalings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Particle simulations of mode conversion between slow mode and fast mode in lower hybrid range of frequencies

The propagation and mode conversion of lower hybrid waves in an inhomogeneous plasma are investigated by using the nonlinear δf algorithm in a two-dimensional particle-in-cell simulation code based on the gyrokinetic electron and fully kinetic ion (GeFi) scheme [Lin et al., Plasma Phys. Controlled Fusion 47, 657 (2005)]. The characteristics of the simulated waves, such as wavelength, frequency, phase, and group velocities, agree well with the linear theoretical analysis. It is shown that a significant reflection component emerges in the conversion process between the slow mode and the fast mode when the scale length of the density variation is comparable to the local wavelength. The dependences of the reflection coefficient on the scale length of the density variation are compared with the results based on the linear full wave model for cold plasmas. It is indicated that the mode conversion for the waves with a frequency of 2.45 GHz (ω ∼ 3ω{sub LH}, where ω{sub LH} represents the lower hybrid resonance) and within Tokamak relevant amplitudes can be well described in the linear scheme. As the frequency decreases, the modification due to the nonlinear term becomes important. For the low-frequency waves (ω ∼ 1.3ω{sub LH}), the generations of the high harmonic modes and sidebands through nonlinear mode-mode coupling provide new power channels and thus could reduce the reflection significantly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Generation of kinetic Alfven waves in the high-latitude near-Earth magnetotail: A global hybrid simulation

In this paper, effects of a fast flow in the tail plasma sheet on the generation of kinetic Alfven waves (KAWs) in the high-latitude of the near-Earth magnetotail are investigated by performing a two-dimensional (2-D) global-scale hybrid simulation, where the plasma flow is initialized by the E×B drift near the equatorial plane due to the existence of the dawn-dusk convection electric field. It is found that firstly, the plasma sheet becomes thinned and the dipolarization of magnetic field appears around (x,z)=(−10.5R{sub E},0.3R{sub E}), where R{sub E} is the radius of the Earth. Then, shear Alfven waves are excited in the plasma sheet, and the strong earthward flow is braked by the dipole-like magnetic field. These waves propagate along the magnetic field lines toward the polar regions later. Subsequently, KAWs with k{sub ⊥}≫k{sub ∥} are generated in the high-latitude magnetotail due to the existence of the non-uniformity of the magnetic field and density in the polar regions. The ratio of the electric field to the magnetic field in these waves is found to obey the relation (δE{sub z})/(δB{sub y} )∼ω/k{sub ∥} of KAWs. Our simulation provides a mechanism for the generation of the observed low-frequency shear Alfven waves in the plasma sheet and kinetic Alfven waves in the high-latitude near-Earth magnetotail, whose source is suggested to be the flow braking in the low-latitude plasma sheet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗