E×B guiding-centre motion in three electrostatic waves: statistics of separatrix crossings

1998 ◽  
Vol 59 (4) ◽  
pp. 761-772
Author(s):  
B. WEYSSOW
Keyword(s):  
Phase Space ◽  
Low Frequency ◽  
Integration Scheme ◽  
Symplectic Integrator ◽  
Third Order ◽  
Special Configuration ◽  
Electrostatic Waves ◽  
Hamiltonian Theory ◽  
Numerical Integration Methods ◽  
Adiabatic Motion

E×B guiding-centre (GC) motion in a special configuration of three low-frequency electrostatic waves can be considered as a paradigmatic Hamiltonian system for studying adiabatic motion and separatrix crossings. A peculiarity of this system is that a single initial condition gives rise to two stroboscopic phase-space trajectories. According to the classical Hamiltonian theory, the proportion of points on the stroboscopic trajectories is a function of the time evolution of the surfaces enclosed by the separatrices in the phase space. This behaviour is qualitatively observed in test-particle numerical experiments. The ability of numerical integration methods like the ‘classical’ fourth-order Runge–Kutta integration scheme or a third-order symplectic integrator to reproduce the statistics is analysed.

scholarly journals Double layers in the downward current region of the aurora

2003 ◽  
Vol 10 (1/2) ◽  
pp. 45-52 ◽  
Author(s):  
R. E. Ergun ◽  
L. Andersson ◽  
C. W. Carlson ◽  
D. L. Newman ◽  
M. V. Goldman
Keyword(s):  
Magnetic Field ◽  
Electron Beam ◽  
Phase Space ◽  
Electric Field ◽  
Electric Fields ◽  
Double Layers ◽  
Parallel Electric Field ◽  
Electrostatic Waves ◽  
Current Region ◽  
Decisive Evidence

Abstract. Direct observations of magnetic-field-aligned (parallel) electric fields in the downward current region of the aurora provide decisive evidence of naturally occurring double layers. We report measurements of parallel electric fields, electron fluxes and ion fluxes related to double layers that are responsible for particle acceleration. The observations suggest that parallel electric fields organize into a structure of three distinct, narrowly-confined regions along the magnetic field (B). In the "ramp" region, the measured parallel electric field forms a nearly-monotonic potential ramp that is localized to ~ 10 Debye lengths along B. The ramp is moving parallel to B at the ion acoustic speed (vs) and in the same direction as the accelerated electrons. On the high-potential side of the ramp, in the "beam" region, an unstable electron beam is seen for roughly another 10 Debye lengths along B. The electron beam is rapidly stabilized by intense electrostatic waves and nonlinear structures interpreted as electron phase-space holes. The "wave" region is physically separated from the ramp by the beam region. Numerical simulations reproduce a similar ramp structure, beam region, electrostatic turbulence region and plasma characteristics as seen in the observations. These results suggest that large double layers can account for the parallel electric field in the downward current region and that intense electrostatic turbulence rapidly stabilizes the accelerated electron distributions. These results also demonstrate that parallel electric fields are directly associated with the generation of large-amplitude electron phase-space holes and plasma waves.

A pseudo third order symplectic integrator

2005 ◽  
Vol 29 (1) ◽  
pp. 92-103
Author(s):  
Liu Fu-yao ◽  
Wu Xin ◽  
Lu Ben-kui
Keyword(s):  
Symplectic Integrator ◽  
Third Order

Multirate Integration for Multibody Dynamic Analysis With Decomposed Subsystems

1999 ◽  
Author(s):  
Sung-Soo Kim ◽  
Jeffrey S. Freeman
Keyword(s):  
Dynamic Analysis ◽  
Multibody System ◽  
Equations Of Motion ◽  
Low Frequency ◽  
Inertia Force ◽  
Integration Scheme ◽  
Integration Algorithm ◽  
Multibody Dynamic ◽  
Higher Order Derivative ◽  
Derivative Information

Abstract This paper details a constant stepsize, multirate integration scheme which has been proposed for multibody dynamic analysis. An Adams-Bashforth Moulton integration algorithm has been implemented, using the Nordsieck form to store internal integrator information, for multirate integration. A multibody system has been decomposed into several subsystems, treating inertia coupling effects of subsystem equations of motion as the inertia forces. To each subsystem, different rate Nordsieck form of Adams integrator has been applied to solve subsystem equations of motion. Higher order derivative information from the integrator provides approximation of inertia force computation in the decomposed subsystem equations of motion. To show the effectiveness of the scheme, simulations of a vehicle multibody system that consists of high frequency suspension motion and low frequency chassis motion have been carried out with different tire excitation forces. Efficiency of the proposed scheme has been also investigated.

scholarly journals Electron phase-space hole transverse instability at high magnetic field

2019 ◽  
Vol 85 (5) ◽  
Author(s):  
I. H. Hutchinson
Keyword(s):  
Phase Space ◽  
Low Frequency ◽  
Force Balance ◽  
Complex Frequency ◽  
Electron Hole ◽  
Particle In Cell ◽  
Planar Electron ◽  
Wave Phenomena ◽  
Slow Growing ◽  
Infinite Particle

Analytic treatment is presented of the electrostatic instability of an initially planar electron hole in a plasma of effectively infinite particle magnetization. It is shown that there is an unstable mode consisting of a rigid shift of the hole in the trapping direction. Its low frequency is determined by the real part of the force balance between the Maxwell stress arising from the transverse wavenumber $k$ and the kinematic jetting from the hole’s acceleration. The very low growth rate arises from a delicate balance in the imaginary part of the force between the passing-particle jetting, which is destabilizing, and the resonant response of the trapped particles, which is stabilizing. Nearly universal scalings of the complex frequency and $k$ with hole depth are derived. Particle in cell simulations show that the slow-growing instabilities previously investigated as coupled hole–wave phenomena occur at the predicted frequency, but with growth rates 2 to 4 times greater than the analytic prediction. This higher rate may be caused by a reduced resonant stabilization because of numerical phase-space diffusion in the simulations.

Optimal Third-Order Symplectic Integration Modeling of Seismic Acoustic Wave Propagation

2020 ◽  
Vol 110 (2) ◽  
pp. 754-762 ◽  
Author(s):  
Chuan Li ◽  
Jianxin Liu ◽  
Bo Chen ◽  
Ya Sun
Keyword(s):  
Wave Propagation ◽  
Wave Equation ◽  
Acoustic Wave ◽  
Seismic Wave ◽  
Truncation Error ◽  
Seismic Wave Propagation ◽  
Symplectic Integrator ◽  
Third Order ◽  
The Third ◽  
2D And 3D

ABSTRACT Seismic wavefield modeling based on the wave equation is widely used in understanding and predicting the dynamic and kinematic characteristics of seismic wave propagation through media. This article presents an optimal numerical solution for the seismic acoustic wave equation in a Hamiltonian system based on the third-order symplectic integrator method. The least absolute truncation error analysis method is used to determine the optimal coefficients. The analysis of the third-order symplectic integrator shows that the proposed scheme exhibits high stability and minimal truncation error. To illustrate the accuracy of the algorithm, we compare the numerical solutions generated by the proposed method with the theoretical analysis solution for 2D and 3D seismic wave propagation tests. The results show that the proposed method reduced the phase error to the eighth-order magnitude accuracy relative to the exact solution. These simulations also demonstrated that the proposed third-order symplectic method can minimize numerical dispersion and preserve the waveforms during the simulation. In addition, comparing different central frequencies of the source and grid spaces (90, 60, and 20 m) for simulation of seismic wave propagation in 2D and 3D models using symplectic and nearly analytic discretization methods, we deduce that the suitable grid spaces are roughly equivalent to between one-fourth and one-fifth of the wavelength, which can provide a good compromise between accuracy and computational cost.

Sign in / Sign up

Export Citation Format

Share Document