Solitary,explosive and periodic solutions for electron acoustic solitary waves with non-thermal hot ions

A theoretical investigation has been made for electron acoustic waves propagating in a system of unmagnetized collisionless plasma consists of a cold electron fluid and ions with two different temperatures in which the hot ions obey the non-thermal distribution. The reductive perturbation method has been employed to derive the Korteweg–de Vries equation for small but finite amplitude electrostatic waves. It is found that the presence of the energetic population of non-thermal hot ions δ, initial normalized equilibrium density of low temperature ions μ and the ratio of temperatures of low temperature ions to high temperature ions β do not only significantly modify the basic properties of solitary structure, but also change the polarity of the solitary profiles. At the critical hot ions density, the KdV equation is not appropriate for describing the system. Hence, a new set of stretched coordinates is considered to derive the modified KdV equation. In the vicinity of the critical hot ions density, neither KdV nor modified KdV equation is appropriate for describing the electron acoustic waves. Therefore, a further modified KdV equation is derived. An algebraic method with computerized symbolic computation, which greatly exceeds the applicability of the existing tanh, extended tanh methods in obtaining a series of exact solutions of the various KdV-type equations, is used here. Numerical studies have been reveals different solutions e.g., bell-shaped solitary pulses, singular solitary “blowup” solutions, Jacobi elliptic doubly periodic wave, Weierstrass elliptic doubly periodic type solutions, in addition to explosive pulses. The results of the present investigation may be applicable to some plasma environments, such as Earth’s magnetotail region.     

Dust acoustic dressed solitons in a four component dusty plasma with nonthermal electron

Nonlinear dust acoustic dressed soliton are studied in a four component dusty plasma. Nonthermal distributions for electrons are considered. The Korteweg–de Vries (KdV) equation is derived by using reductive perturbation technique. A higher order inhomogeneous differential equation is obtained for the higher order correction. The expression for dressed soliton is obtained by the renormalization method. The expressions for higher order correction are determined using a series solution technique.     

Effects of adiabaticity and non-adiabaticity of dust charge fluctuation on the propagation of the dust ion acoustic waves in inhomogeneous dusty plasma

A theoretical investigation has been made for adiabatic positive and negative dust charge fluctuations on the propagation of dust-ion acoustic waves (DIAWs) in a weakly inhomogeneous, collisionless, unmagnetized dusty plasmas consisting of cold positive ions, stationary positively and negatively charged dust particles and isothermal electrons. The reductive perturbation method is employed to reduce the basic set of fluid equations to the variable coefficients Korteweg–de Vries (KdV) equation. Either compressive or rarefactive solitons are shown to exist depending on the critical value of the ion density, which in turn, depends on the inhomogeneous distribution of the ion. The dissipative effects of non-adiabatic dust charge variation has been studied which cause generation of dust ion acoustic shock waves governed by KdV-Burger (KdVB) equation. The results of the present investigation may be applicable to some dusty plasma environments, such as dusty plasma existing in polar mesosphere region.     

Propagation of non-linear ion-sound pulses in dusty plasma waveguides

There are a lot of objects in space associated with dusty plasma inclusions. Such inclusions may bear a prolonged shape and behave as waveguides for ion-sound waves. In the case of space plasmas, the dust particles can possess both negative charge, due to electron attachment, and positive one, due to photoionization. In this paper the propagation of linear and non-linear ion-sound wave pulses in the dusty plasma waveguides, possessing positive charge, is studied. It has been demonstrated that non-linear dynamics of baseband pulse propagation in plasma waveguide possesses essentially non-solitonic behavior. Namely, propagation of a long ion-sound pulse leads to an excitation of a shock-like wave but not a stable localized nonlinear pulse. Also, when a Korteveg–de Vries (KdV) soliton is incident onto the dusty plasma waveguide, some part of the soliton energy is captured by the waveguide and transformed into a multi-pulse structure. Additionally, an interaction of dusty plasma inclusions with KdV soliton can lead to the occurrence of transverse instabilities of the soliton and its eventual destruction.     

Propagation of cylindrical and spherical dust–ion acoustic solitary waves in a relativistic dusty plasma

The properties of cylindrical and spherical dust–ion acoustic solitary waves (DIASW) in an unmagnetized dusty plasma comprising of relativistic ions, Boltzmann electrons, and stationary dusty particles are investigated. Under a suitable coordinate transformation, the cylindrical KdV equation can be solved analytically. The change of the DIASW structure due to the effect of geometry, relativistic streaming factor, ion density and electron temperature is studied by numerical calculation of the cylindrical/spherical Kdv equation. It is noted that with ion pressure the effect of relativistic streaming factor to solitary waves structure is different. Without ion pressure, as the relativistic streaming factor decreases, the amplitude of the solitary wave decreases. However, when the ion pressure is taken into account, the amplitude decreases as the relativistic streaming factor increases and is highly sensitive to relativistic streaming factor. Our results may have relevance in the understanding of astrophysical plasmas.     

Numerical simulation of equatorial plasma bubbles over Cachimbo: COPEX campaign

The problem of day-to-day variability in onset of equatorial spread F (ESF) is addressed using data from the 2002 COPEX observational campaign in Brazil and numerical modeling. The observational results show that for values of virtual height of the F layer base less than 355 km at around 18:35 LT, and for the prereversal peak enhancement of the vertical plasma drift (Vp) less than 30 m/s, the spread-F (ESF) was absent on four nights over Cachimbo (9.5°S, 54.8°W, dip latitude = −2.1°). In this work we analyze the geophysical conditions for the generation of the irregularities by comparing the nights with and without the ESF. In the comparison a numerical code is used to simulate plasma irregularity development in an extended altitude range from the bottom of the equatorial F   layer. The code uses the flux corrected transport method with Boris–Book’s flux limiter for the spatial integration and a predictor–corrector method for the direct time integration of the continuity equation for O⁺${\mathrm{O}}^{+}$ and the SOR (Successive-Over-Relaxation) method for electric potential equation. The code is tested with different evening eastward electric fields (or vertical drifts Vp < 30 m/s and Vp > 30 m/s) in order to study the influence of the prereversal enhancement in the zonal electric field on plasma bubble formation and development. The code also takes into account the zonal wind, the vertical electric field and the collision frequency of ions with neutrals and the amplitude of initial perturbation. The simulation shows a good agreement with the observational results of the ESF. The results of the code suggest that the instability can grow at the F layer bottomside by the Rayleigh–Taylor mechanism only when the Vp > 30 m/s. In the analyzed cases we have considered the competition of other geophysical parameters in the generation of plasma structures.     

On Lambert’s problem and the elliptic time of flight equation: A simple semi-analytical inversion method

In the two-body model, time of flight between two positions can be expressed as a single-variable function and a variety of formulations exist. Lambert’s problem can be solved by inverting such a function. In this article, a method which inverts Lagrange’s flight time equation and supports the problematic 180°$180°$ transfer is proposed. This method relies on a Householder algorithm of variable order. However, unlike other iterative methods, it is semi-analytical in the sense that flight time functions are derived analytically to second order vs. first order finite differences. The author investigated the profile of Lagrange’s elliptic flight time equation and its derivatives with a special focus on their significance to the behaviour of the proposed method and the stated goal of guaranteed convergence. Possible numerical deficiencies were identified and dealt with. As a test, 28 scenarios of variable difficulty were designed to cover a wide variety of geometries. The context of this research being the orbit determination of artificial satellites and debris, the scenarios are representative of typical such objects in Low-Earth, Geostationary and Geostationary Transfer Orbits. An analysis of the computational impact of the quality of the initial guess vs. that of the order of the method was also done, providing clues for further research and optimisations (e.g. asteroids, long period comets, multi-revolution cases). The results indicate fast to very fast convergence in all test cases, they validate the numerical safeguards and also give a quantitative assessment of the importance of the initial guess.     

In situ measurements of ionospheric plasma turbulence over five frequency decades: Heritage flight of the Plasma Local Anomalous Noise Experiment (PLANE)

Observations of ionospheric plasma density and frequency-dependent broadband plasma turbulence made during the heritage flight of the Plasma Local Anomalous Noise Experiment (PLANE) are presented. Rather than record high frequency time series data, the experiment was designed to record Power Spectral Distributions (PSDs) in five decadal frequency bins with upper limits ranging from 1.0 Hz to 10 kHz. Additionally, PLANE was designed distinguish turbulence in the ambient plasma from that local to the spacecraft. The instrument consists of two retarding potential analyzers (RPAs) connected together via a feedback loop to force one analyzer into the I–V trace retardation region at all times. Fluctuations in this measurement are believed to be ambient only as the RPA’s voltage would be too high for locally turbulent plasma to surmount the potential barrier, which is nominally at ram energy. The instrument requires pointing along the spacecraft’s ram velocity vector to make this measurement, thus requiring stabilization in pitch and yaw. During PLANE’s heritage flight, though the satellite’s attitude control system failed early in the mission, plasma data were collected during opportune times in which the instrument rotated into and out of the ram. Observations of plasma density and PSDs of high frequency plasma turbulence were recorded on several occasions. Additionally, a plasma source onboard the satellite was used to generate artificial plasma turbulence, and the PLANE data observed periodic structure presumably associated with the rotation of the spacecraft during these source firings. A brief comparison with other high frequency in situ plasma instruments is presented.     

Interaction of magnetoacoustic solitons in electron-positron plasmas

The interaction between two, four and six magnetoacoustic solitons in electron-positron plasmas are investigated. The extended Poincaré–Lighthill–Kuo (PLK) perturbation method is employed to derived two KdV equations for magnetoacoustic solitons moving towards each other and studied the head-on collision between them and their phase shifts. The Hirota bilinear method is used to have multi-soliton solutions of already derived two KdV equations for right and left moving solitons. The four and six magnetoacoustic solitons solutions of the two KdV equations are obtained to discuss their interaction and phase shifts. It is found that only compressive magnetoacoustic solitons structures are formed in electron-positron plasma. The present study may be useful to understand the collective phenomena related to head-on and overtaking magnetoacoustic solitons interaction in electron-positron plasmas that may occur in a pulsar magnetosphere.     

A lunar cargo mission design strategy using variable low thrust

A design technique for a near optimal, Earth–Moon transfer trajectory using continuous variable low thrust is proposed. For the Earth–Moon transfer trajectory, analytical and numerical methods are combined to formulate the trajectory optimization problem. The basic concept of the proposed technique is to utilize analytically optimized solutions when the spacecraft is flying near a central body where the transfer trajectories are nearly circular shaped, and to use a numerical optimization method to match the spacecraft’s states to establish a final near optimal trajectory. The plasma thruster is considered as the main propulsion system which is currently being developed for crewed/cargo missions for interplanetary flight. The gravitational effects of the 3rd body and geopotential effects are included during the trajectory optimization process. With the proposed design technique, Earth–Moon transfer trajectory is successfully designed with the plasma thruster having a thrust direction sequence of “fixed-varied-fixed” and a thrust acceleration sequence of “constant-variable-constant”. As this strategy has the characteristics of a lesser computational load, little sensitivity to initial conditions, and obtaining solutions quickly, this method can be utilized in the initial scoping studies for mission design and analysis. Additionally, derived near optimal trajectory solution can be used as for initial trajectory solution for further detailed optimization problem. The demonstrated results will give various insights into future lunar cargo trajectories using plasma thrusters with continuous variable low thrust, establishing approximate costs as well as trajectory characteristics.     

Universal scaling laws for fully-developed magnetic field turbulence near and far upstream of the Earth’s bow shock

We analyze the multifractal scaling of the modulus of the interplanetary magnetic field near and far upstream of the Earth’s bow shock, measured by Cluster and ACE, respectively, from 1 to 3 February 2002. The maximum order of the structure function is carefully estimated for each time series using two different techniques, to ensure the validity of our high-order statistics. The first technique consists of plotting the integrand of the pth order structure function, and the second technique is a quantitative method which relies on the power-law scaling of the extreme events. We compare the scaling exponents computed from the structure functions of magnetic field differences with the predictions obtained by the She–Lévêque model of intermittency in anisotropic magnetohydrodynamic turbulence. Our results show a good agreement between the model and the observations near and far upstream of the Earth’s bow shock, rendering support for the modelling of universal scaling laws based on the Kolmogorov phenomenology in the presence of sheet-like dissipative structures.     

MHD modeling of the outer heliosphere: Achievements and challenges

An overview is presented of magnetic-field-related effects in the solar wind (SW) interaction with the local interstellar medium (LISM) and the different theoretical approaches used in their investigation. We discuss the possibility that the interstellar magnetic field (ISMF) introduces north–south and east–west asymmetries of the heliosphere, which might explain observational data obtained by the Voyager 1 and Voyager 2 spacecraft. The SW–LISM interaction parameters that are responsible for the deflection of the interstellar neutral hydrogen flow from the direction of propagation of neutral helium in the inner heliosheath are outlined. The possibility of a strong ISMF, which increases the heliospheric asymmetry and the H–He flow deflection, is discussed. The effect of the combination of a slow-fast solar wind during solar minimum over the Sun’s 11-year activity cycle is illustrated. The consequences of a tilt between the Sun’s magnetic and rotational axes are analyzed. Band-like areas of an increased magnetic field distribution in the outer heliosheath are sought in order to discover regions of possible 2–3 kHz radio emission.     

The International DORIS Service (IDS): Toward maturity

DORIS is one of the four space-geodetic techniques participating in the Global Geodetic Observing System (GGOS), particularly to maintain and disseminate the Terrestrial Reference Frame as determined by International Earth rotation and Reference frame Service (IERS). A few years ago, under the umbrella of the International Association of Geodesy, a DORIS International Service (IDS) was created in order to foster international cooperation and to provide new scientific products. This paper addresses the organizational aspects of the IDS and presents some recent DORIS scientific results. It is for the first time that, in preparation of the ITRF2008, seven Analysis Centers (AC’s) contributed to derive long-term time series of DORIS stations positions. These solutions were then combined into a homogeneous time series IDS-2 for which a precision of less than 10 mm was obtained. Orbit comparisons between the various AC’s showed an excellent agreement in the radial component, both for the SPOT satellites (e.g. 0.5–2.1 cm RMS for SPOT-2) and Envisat (0.9–2.1 cm RMS), using different software packages, models, corrections and analysis strategies. There is now a wide international participation within IDS that should lead to future improvements in DORIS analysis strategies and DORIS-derived geodetic products.     

The realization and convergence analysis of combined PPP based on raw observation

In order to speed up Precise Point Positioning (PPP)’s convergence, a combined PPP method with GPS and GLONASS which is based on using raw observations is proposed, and the positioning results and convergence time have been compared with that of single system. The ionospheric delays and receiver’s Differential Code Bias (DCB) corrections are estimated as unknown parameters in this method. The numerical results show that the combined PPP has not caused significant impacts on the final solutions, but it greatly improved Position Dilution of Precision (PDOP) and convergence speed and enhanced the reliability of the solution. Meanwhile, the convergence speed is greatly influenced by the receiver’s DCB, positioning results in horizontal which are better than 10 cm can be realized within 10 min. In addition, the ionosphere and DCB products can be provided with high precision.     

Model of energetic populations at Ganymede,implications for an orbiter

A model is developed to study the energetic particle populations in Ganymede’s magnetosphere. The main objective is to estimate to what extent the moon could protect an orbiter from radiations. Using Liouville’s theorem, the phase space density of particles coming from Jupiter’s magnetosphere is calculated at any point of Ganymede’s environment. Up to energies of ∼50–100 keV for ions and ∼10–20 MeV for electrons, Ganymede’s magnetic field appears to be able to form distinctive populations as loss-cones over the polar caps and radiation belts. At larger energies, these features are blurred by Larmor radius effects; the moon absorption simply creates a quasi-isotropic layer of ∼500 km thickness where the flux is reduced by ∼40–50%. The predictions are compared to Galileo measurements. In particular, we demonstrate the importance of the moon sweeping in reducing the flux over the polar caps. Interestingly, this can be accounted for by assuming that the particles bouncing between Jupiter and Ganymede are ideally scattered in pitch angle and permanently re-fill the loss-cone, which increases the precipitation on Ganymede’s polar cap. In overall, it is estimated that the radiation dose received by an orbiter of Ganymede will be reduced by more than 50–60% compared to the expected dose at Jupiter/Ganymede distance. This should have a positive impact on the design of a future orbiter of Ganymede.     

GOCE orbit analysis: Long-wavelength gravity field determination using the acceleration approach

The restricted sensitivity of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) gradiometer instrument requires satellite gravity gradiometry to be supplemented by orbit analysis in order to resolve long-wavelength features of the geopotential. For the hitherto published releases of the GOCE time-wise (TIM) and GOCE space-wise gravity field series—two of the official ESA products—the energy conservation method has been adopted to exploit GPS-based satellite-to-satellite tracking information. On the other hand, gravity field recovery from data collected by the CHAllenging Mini-satellite Payload (CHAMP) satellite showed the energy conservation principle to be a sub-optimal choice. For this reason, we propose to estimate the low-frequency part of the gravity field by the point-wise solution of Newton’s equation of motion, also known as the acceleration approach. This approach balances the gravitational vector with satellite accelerations, and hence is characterized by (second-order) numerical differentiation of the kinematic orbit. In order to apply the method to GOCE, we present tailored processing strategies with regard to low-pass filtering, variance–covariance information handling, and robust parameter estimation. By comparison of our GIWF solutions (initials GI for “Geodätisches Institut” and IWF for “Institut für WeltraumForschung”) and the GOCE-TIM estimates with a state-of-the-art gravity field solution derived from GRACE (Gravity Recovery And Climate Experiment), we conclude that the acceleration approach is better suited for GOCE-only gravity field determination as opposed to the energy conservation method.     

Small amplitude ion acoustic solitons in a weakly magnetized plasma with anisotropic ion pressure and kappa distributed electrons

The Zakharov–Kuznetzov (ZK) equation is derived for nonlinear electrostatic waves in a weakly magnetized plasma in the presence of anisotropic ion pressure and superthermal electrons. The anisotropic ion pressure is defined using Chew–Goldberger–Low (CGL) while a generalized Lorentzian (kappa) distribution is assumed for the non-thermal electrons. The standard reductive perturbation method (RPM) is employed to derive the two dimensional ZK equation for the dynamics of obliquely propagating low frequency ion acoustic wave. The influence of spectral index (kappa) of non-thermal electron on the soliton is discussed in the presence of anisotropic ion pressure in plasmas. It is found that ion pressure anisotropy and superthermality of electrons affect both the width and amplitude of the solitary waves. On the other hand the magnetic field is found to alter the dispersive property of the plasma only, and hence the width of the solitons is affected while the amplitude of the solitary waves is independent of external magnetic field. The numerical results are also presented for illustrations.     

The fast debris evolution model

The ‘particles-in-a-box’ (PIB) model introduced by Talent [Talent, D.L. Analytic model for orbital debris environmental management. J. Spacecraft Rocket, 29 (4), 508–513, 1992.] removed the need for computer-intensive Monte Carlo simulation to predict the gross characteristics of an evolving debris environment. The PIB model was described using a differential equation that allows the stability of the low Earth orbit (LEO) environment to be tested by a straightforward analysis of the equation’s coefficients. As part of an ongoing research effort to investigate more efficient approaches to evolutionary modelling and to develop a suite of educational tools, a new PIB model has been developed. The model, entitled Fast Debris Evolution (FADE), employs a first-order differential equation to describe the rate at which new objects ?10 cm are added and removed from the environment. Whilst Talent [Talent, D.L. Analytic model for orbital debris environmental management. J. Spacecraft Rocket, 29 (4), 508–513, 1992.] based the collision theory for the PIB approach on collisions between gas particles and adopted specific values for the parameters of the model from a number of references, the form and coefficients of the FADE model equations can be inferred from the outputs of future projections produced by high-fidelity models, such as the DAMAGE model.     

Studies of the geocenter motion using 16-years DORIS data

An accuracy of geocenter motion estimation is strongly dependent on the geodetic network size and stations distribution over the Earth’s surface. From this point of view DORIS system has an advantage, as its ground network of beacons consists of more than 50 sites, equally distributed over the Earth’s surface. Aiming to study variations of the geocenter movements, the results of DORIS data analysis for the time span 1993.0–2009.0 (inawd06.snx series), performed at the Analysis Centre of the Institute of astronomy of the Russian Academy of Sciences, have been used. DORIS data processing was made with GIPSY/OASIS II software, developed by Jet Propulsion Laboratory and modified for DORIS data processing by Institute Géographique National. Standard deviations of stations coordinates are estimated at the level 0.5–4.0 cm (internal consistency), depending on the number of satellites used in the solution. RMS of estimated components of the DORIS satellites orbits, compared with the solutions of other IDS analysis centres, do not exceed 1–2 cm. Weekly solutions for coordinates have been transformed from free network solutions (inawd06.snx series) to a well defined terrestrial reference frame ITRF2005 with the use of seven parameters of Helmert transformation, which were examined with a view to study variations of the geocenter movements (ina05wd01.geoc time series). In order to estimate linear trend, amplitudes, periods and phases of geocenter variation a method of linear regression was applied. The evaluated amplitudes of annual variations are of the order of 5–7 mm for X and Y components and 27–29 mm for Z component. Semi-annual amplitudes are also noticeable in all components (1–34 mm for X, Y and Z components). Secular trends in the DORIS geocenter coordinates are: −1.2, −0.1 and −0.3 mm/year for X, Y and Z directions respectively.