Simulation study of motion of charged particles trapped in Earth’s magnetosphere

This article aims to understand the motion of the charged particles trapped in the Earth’s inner magnetosphere. The emphasis is on identifying the numerical scheme, which is appropriate to characterize the trajectories of the charged particles of different energies that enter the Earth’s magnetosphere and get trap along the magnetic field lines. These particles perform three different periodic motions, namely: gyration, bounce, and azimuthal drift. However, often, the gyration of the particle is ignored, and only the guiding center of the particle is traced to reduce the computational time. It is because the simulation of all three motions (gyro, bounce, and drift) together needed a robust numerical scheme, which has less numerical dissipation. We have developed a three-dimensional test particle simulation model in which the relativistic equation of motion is solved numerically using the fourth and sixth-order Runge-Kutta methods. The stability of the simulation model is verified by checking the conservation of total kinetic energy and adiabatic invariants linked with each type of motion. We found that the sixth-order Runge-Kutta method is suitable to trace the complete trajectories of both proton and electron of a wide energy range, 5 keV to 250 MeV for L = 2 – 6. We have estimated the bounce and drift periods from the simulations, and they are found to be in good agreement with the theory. The study implies that a simulation model with sixth-order Runge-Kutta method can be applied to the time-vary, non-analytical form of magnetic configuration in future studies to understand the dynamics of charged particles trapped in Earth’s magnetosphere.     

Ion- and electron-acoustic solitons and double layers in multi-component space plasmas

A general model for the ion- and electron-acoustic solitons and double layers in a multi-component unmagnetized plasma consisting of background electrons, counter-streaming electron beams and ions is discussed. The model is based on the multi-fluid equations and the Poisson equation, and uses the Sagdeev pseudo-potential techniques. For identical counter-streaming electron beams and depending upon the plasma parameters, three types of solutions, namely, ion-acoustic, slow and fast electron-acoustic soliton/double layer, are possible. Generally, the ion acoustic solitons have positive potentials, slow-electron acoustic solitons have negative potentials and fast electron-acoustic solitons and double layers can have either positive or negative potentials depending on the core electron density. As beam speed is increased, first ion-acoustic and then slow electron-acoustic solitons disappear. At large beam speed, only fast electron-acoustic solitons/double layers survive. The results may be relevant to the observations of the electrostatic solitary waves (ESWs) observed in the Earth’s magnetosphere.     

Electron acoustic solitary waves in the Earth’s magnetotail region

Satellite observations have revealed solitary potential structures in the Earth’s magnetotail region. These structures have both positive (compressive) and negative (rarefactive) electrostatic potentials. In this paper we study the electron-acoustic solitary waves (EASWs) in an unmagnetized plasma consisting of cold plasma electrons and isothermal ions with two different temperatures. Using the reductive perturbation method, the nonlinear evolution of such structures is studied. The numerical computations are performed to study the role of two temperature ions in the generation of EASWs. In this case, the model supports the existence of both positive and negative electrostatic potentials with bipolar pulses. The electric field associated with these positive and negative solitary structures are numerically computed. The present study could be useful to construe the compressive and rarefactive electric field bipolar pulses associated with the BEN type emissions in the magnetospheric regions where the electron beams are not present.     

Detumbling and reorientation maneuvers and stabilization of NASASCOLE system

The questions of rotational maneuver and vibration stabilization of the NASA Spacecraft Control Laboratory Experiment (SCOLE) system is considered. The mathematical model of the SCOLE system includes the rigid body dynamics as well as the elastic dynamics representing transverse and torsional deformations of the elastic beam connecting the orbiter and end body (reflector). For the rotational maneuver, a new control law (orbiter control law) is derived using an orbiter input torque vector. Detumbling and reorientation maneuvers of the SCOLE system are accomplished using this control law; however, this excites the elastic modes of the beam. The orbiter control law asymptotically linearizes the flexible dynamics. Using the linearized model, a linear feedback control law is designed for vibration suppression. An observer is designed for estimating the state variables using sensor outputs which are also used for the synthesis of the control law. Simulation results are presented to show that in the closed-loop system detumbling and reorientation maneuvers can be accomplished and the effect of control and observation spillover is insignificant