Finite amplitude envelope solitons in pair-ion plasmas

The nonlinear coupling between finite amplitude ion thermal waves (ITWs) and quasistationary density perturbations in a pair-ion plasma is considered. A generalized nonlinear Schrödinger equation is derived for the ITW electric field envelope, accounting for large amplitude quasistationary plasma slow motion describing the ITW ponderomotive force. The present theory accounts for the trapping of ITWs in a large amplitude ion density hole. The small amplitude limit is considered and exact analytical solutions are obtained. Finite amplitude solutions are obtained numerically and their characteristics are discussed.

Oblique amplitude modulation of dust-acoustic plasma waves

Theoretical and numerical studies are presented of the nonlinear amplitude modulation of dust-acoustic (DA) waves propagating in an unmagnetized three component, weakly-coupled, fully ionized plasma consisting of electrons, positive ions and charged dust particles, considering perturbations oblique to the carrier wave propagation direction. The stability analysis, based on a nonlinear Schrodinger-type equation (NLSE), shows that the wave may become unstable; the stability criteria depend on the angle theta between the modulation and propagation directions. Explicit expressions for the instability rate and threshold have been obtained in terms of the dispersion laws of the system. The possibility and conditions for the existence of different types of localized excitations have also been discussed.

Modeling relativistic soliton interactions in overdense plasmas: A perturbed nonlinear Schrodinger equation framework

We investigate the dynamics of localized solutions of the relativistic cold-fluid plasma model in the small but finite amplitude limit, for slightly overcritical plasma density. Adopting a multiple scale analysis, we derive a perturbed nonlinear Schrodinger equation that describes the evolution of the envelope of circularly polarized electromagnetic field. Retaining terms up to fifth order in the small perturbation parameter, we derive a self-consistent framework for the description of the plasma response in the presence of localized electromagnetic field. The formalism is applied to standing electromagnetic soliton interactions and the results are validated by simulations of the full cold-fluid model. To lowest order, a cubic nonlinear Schrodinger equation with a focusing nonlinearity is recovered. Classical quasiparticle theory is used to obtain analytical estimates for the collision time and minimum distance of approach between solitons. For larger soliton amplitudes the inclusion of the fifth-order terms is essential for a qualitatively correct description of soliton interactions. The defocusing quintic nonlinearity leads to inelastic soliton collisions, while bound states of solitons do not persist under perturbations in the initial phase or amplitude

Amplitude modulation of quantum-ion-acoustic wavepackets in electron-positron-ion plasmas: Modulational instability, envelope modes, extreme waves

A semirelativistic fluid model is employed to describe the nonlinear amplitude modulation of low-frequency (ionic scale) electrostatic waves in an unmagnetized electron-positron-ion plasma. Electrons and positrons are assumed to be degenerated and inertialess, whereas ions are warm and classical. A multiscale perturbation method is used to derive a nonlinear Schrödinger equation for the envelope amplitude, based on which the occurrence of modulational instability is investigated in detail. Various types of localized ion acoustic excitations are shown to exist, in the form of either bright type envelope solitons (envelope pulses) or dark-type envelope solitons (voids, holes). The plasma configurational parameters (namely, the relativistic degeneracy parameter, the positron concentration, and the ionic temperature) are shown to affect the conditions for modulational instability significantly, in fact modifying the associated threshold as well as the instability growth rate. In particular, the relativistic degeneracy parameter leads to an enhancement of the modulational instability mechanism. Furthermore, the effect of different relevant plasma parameters on the characteristics (amplitude, width) of these envelope solitary structures is also presented in detail. Finally, the occurrence of extreme amplitude excitation (rogue waves) is also discussed briefly. Our results aim at elucidating the formation and dynamics of nonlinear electrostatic excitations in superdense astrophysical regimes.

A Schamel equation for ion acoustic waves in superthermal plasmas

An investigation of the propagation of ion acoustic waves in nonthermal plasmas in the presence of trapped electrons has been undertaken. This has been motivated by space and laboratory plasma observations of plasmas containing energetic particles, resulting in long-tailed distributions, in combination with trapped particles, whereby some of the plasma particles are confined to a finite region of phase space. An unmagnetized collisionless electron-ion plasma is considered, featuring a non-Maxwellian-trapped electron distribution, which is modelled by a kappa distribution function combined with a Schamel distribution. The effect of particle trapping has been considered, resulting in an expression for the electron density. Reductive perturbation theory has been used to construct a KdV-like Schamel equation, and examine its behaviour. The relevant configurational parameters in our study include the superthermality index κ and the characteristic trapping parameter β. A pulse-shaped family of solutions is proposed, also depending on the weak soliton speed increment u₀. The main modification due to an increase in particle trapping is an increase in the amplitude of solitary waves, yet leaving their spatial width practically unaffected. With enhanced superthermality, there is a decrease in both amplitude and width of solitary waves, for any given values of the trapping parameter and of the incremental soliton speed. Only positive polarity excitations were observed in our parametric investigation. 

Electron acoustic soliton energy of the Kadomtsev-Petviashvili equation in the Earth’s magnetotail region at critical ion density

The nonlinear properties of small amplitude electron-acoustic solitary waves (EAWs) in a homogeneous system of unmagnetized collisionless plasma consisted of a cold electron fluid and isothermal ions with two different temperatures obeying Boltzmann type distributions have been investigated. A reductive perturbation method was employed to obtain the Kadomstev-Petviashvili (KP) equation. At the critical ion density, the KP equation is not appropriate for describing the system. Hence, a new set of stretched coordinates
is considered to derive the modified KP equation. Moreover, the solitary solution, soliton energy and the associated electric field at the critical ion density were computed. The present investigation can be of relevance to the electrostatic solitary structures observed in various space plasma environments, such as Earth’s magnetotail region.

Difference equation approach to two-thermocouple sensor characterisation in constant velocity flow environments

Thermocouples are one of the most popular devices for temperature measurement due to their robustness, ease of manufacture and installation, and low cost. However, when used in certain harsh environments, for example, in combustion systems and engine exhausts, large wire diameters are required, and consequently the measurement bandwidth is reduced. This article discusses a software compensation technique to address the loss of high frequency fluctuations based on measurements from two thermocouples. In particular, a difference equation sDEd approach is proposed and compared with existing methods both in simulation and on experimental test rig data with constant flow velocity. It is found that the DE algorithm, combined with the use of generalized total least squares for parameter identification, provides better performance in terms of time constant estimation without any a priori assumption on the time constant ratios of the thermocouples.

Particle Image Velocimetry Investigation of a finite amplitude pressure wave

Particle image velocimetry is used to study the motion of gas within a duct subject to the passage of a finite amplitude pressure wave. The wave is representative of the pressure waves found in the exhaust systems of internal combustion engines. Gas particles are accelerated from stationary to 150 m/s and then back to stationary in 8 ms. It is demonstrated that gas particles at the head of the wave travel at the same velocity across the duct cross section at a given point in time. Towards the tail of the wave viscous effects are plainly evident causing the flow profile to tend towards parabolic. However, the instantaneous mean particle velocity across the section is shown to match well with the velocity calculated from a corresponding measured pressure history using 1D gas dynamic theory. The measured pressure history at a point in the duct was acquired using a high speed pressure transducer of the type typically used for engine research in intake and exhaust systems. It is demonstrated that these are unable to follow the rapid changes in pressure accurately and that they are prone to resonate under certain circumstances.

Exploiting a priori time constant ratio information in difference equation two-thermocouple sensor characterisation

The characterization of thermocouple sensors for temperature measurement in varying-flow environments is a challenging problem. Recently, the authors introduced novel difference-equation-based algorithms that allow in situ characterization of temperature measurement probes consisting of two-thermocouple sensors with differing time constants. In particular, a linear least squares (LS) lambda formulation of the characterization problem, which yields unbiased estimates when identified using generalized total LS, was introduced. These algorithms assume that time constants do not change during operation and are, therefore, appropriate for temperature measurement in homogenous constant-velocity liquid or gas flows. This paper develops an alternative ß-formulation of the characterization problem that has the major advantage of allowing exploitation of a priori knowledge of the ratio of the sensor time constants, thereby facilitating the implementation of computationally efficient algorithms that are less sensitive to measurement noise. A number of variants of the ß-formulation are developed, and appropriate unbiased estimators are identified. Monte Carlo simulation results are used to support the analysis.

Plane wave and Coulomb asymptotics

A simple plane wave solution of the Schrodinger-Helmholtz equation is a quantum eigenfunction obeying both energy and linear momentum correspondence principles. Inclusion of the outgoing wave with scattering amplitude f asymptotic development of the plane wave, we show that there is a problem with angular momentum when we consider forward scattering at the point of closest approach and at large impact parameter given semiclassically by (l + 1/2)/k where l is the azimuthal quantum number and may be large (J. Leech et al., Phys. Rev. Lett. 88. 257901 (2002)). The problem is resolved via non- uniform, non-standard analysis involving the Heaviside step function, unifying classical, semiclassical and quantum mechanics, and the treatment is extended to the case of pure Coulomb scattering.

Modulational instability and localized excitations of dust-ion acoustic waves

Theoretical and numerical studies are carried out of the nonlinear amplitude modulation of dust-ion acoustic waves propagating in an unmagnetized weakly coupled plasma comprised of electrons, positive ions, and charged dust grains, considering perturbations oblique to the carrier wave propagation direction. The stability analysis, based on a nonlinear Schrodinger-type equation, exhibits a wide instability region, which depends on both the angle theta between the modulation and propagation directions and the dust number density n(d). Explicit expressions for the instability increment and threshold are obtained. The possibility and conditions for the existence of different types of localized excitations are also discussed. (C) 2003 American Institute of Physics.

Solving the advection-dispersion equation using the discrete puff particle method

A flexible, mass-conservative numerical technique for solving the advection-dispersion equation for miscible contaminant transport is presented. The method combines features of puff transport models from air pollution studies with features from the random walk particle method used in water resources studies, providing a deterministic time-marching algorithm which is independent of the grid Peclet number and scales from one to higher dimensions simply. The concentration field is discretised into a number of particles, each of which is treated as a point release which advects and disperses over the time interval. The dispersed puff is itself discretised into a spatial distribution of particles whose masses can be pre-calculated. Concentration within the simulation domain is then calculated from the mass distribution as an average over some small volume. Comparison with analytical solutions for a one-dimensional fixed-duration concentration pulse and for two-dimensional transport in an axisymmetric flow field indicate that the algorithm performs well. For a given level of accuracy the new method has lower computation times than the random walk particle method.

A Van der Pol-Mathieu equation for the dynamics of dust grain charge in dusty plasmas

The chaotic profile of dust grain dynamics associated with dust-acoustic oscillations in a dusty plasma is considered. The collective behaviour of the dust plasma component is described via a multi-fluid model, comprising Boltzmann distributed electrons and ions, as well as an equation of continuity possessing a source term for the dust grains, the dust momentum and Poisson's equations. A Van der Pol–Mathieu-type nonlinear ordinary differential equation for the dust grain density dynamics is derived. The dynamical system is cast into an autonomous form by employing an averaging method. Critical stability boundaries for a particular trivial solution of the governing equation with varying parameters are specified. The equation is analysed to determine the resonance region, and finally numerically solved by using a fourth-order Runge–Kutta method. The presence of chaotic limit cycles is pointed out.