Dynamic evolution of a quasi-spherical general polytropic magnetofluid with self-gravity

In various astrophysical contexts, we analyze self-similar behaviours of magnetohydrodynamic (MHD) evolution of a quasi-spherical polytropic magnetized gas under self-gravity with the specific entropy conserved along streamlines. In particular, this MHD model analysis frees the scaling parameter n in the conventional polytropic self-similar transformation from the constraint of n+γ=2 with γ being the polytropic index and therefore substantially generalizes earlier analysis results on polytropic gas dynamics that has a constant specific entropy everywhere in space at all time. On the basis of the self-similar nonlinear MHD ordinary differential equations, we examine behaviours of the magnetosonic critical curves, the MHD shock conditions, and various asymptotic solutions. We then construct global semi-complete self-similar MHD solutions using a combination of analytical and numerical means and indicate plausible astrophysical applications of these magnetized flow solutions with or without MHD shocks.     

Study of inertial kinetic Alfven waves around cusp region

The effect of electron inertia on kinetic Alfven wave has been studied. The expressions for the dispersion relation, growth/damping rate and growth/damping length of the inertial kinetic Alfven wave (IKAW) are derived using the kinetic approach in cusp region. The Vlasov-kinetic theory has been adopted to evaluate the dispersion relation, growth/damping rate and growth/damping length with respect to the perpendicular wave number k_⊥ρ_i (ρ_i is the ion gyroradius) at different plasma densities. The growth/damping rate and growth/damping length are evaluated for different m_e/βm_i, where β is the ratio of electron pressure to the magnetic field pressure, m_{i, e} are the mass of ion and electron, respectively, as I=m_e/βm_i represent boundary between the kinetic and inertial regimes. It is observed that frequency of inertial kinetic Alfven wave (IKAW) ω is decreasing with k_⊥ρ_i and plasma density. The polar cusp is an ideal laboratory for studies of nonlinear plasma processes important for understanding the basic plasma physics, as well as the magnetospheric and astrophysical applications of these processes.     

Variation of the gas and radiation content in the sub-Keplerian accretion disk around black holes and its impact to the solutions

We investigate the variation of the gas and the radiation pressure in accretion disks during the infall of matter to the black hole and its effect to the flow. While the flow far away from the black hole might be non-relativistic, in the vicinity of the black hole it is expected to be relativistic behaving more like radiation. Therefore, the ratio of gas pressure to total pressure (β) and the underlying polytropic index (γ) should not be constant throughout the flow. We obtain that accretion flows exhibit significant variation of β and then γ, which affects solutions described in the standard literature based on constant β. Certain solutions for a particular set of initial parameters with a constant β do not exist when the variation of β is incorporated appropriately. We model the viscous sub-Keplerian accretion disk with a nonzero component of advection and pressure gradient around black holes by preserving the conservations of mass, momentum, energy, supplemented by the evolution of β. By solving the set of five coupled differential equations, we obtain the thermo-hydrodynamical properties of the flow. We show that during infall, β of the flow could vary up to ∼300%, while γ up to ∼20%. This might have a significant impact to the disk solutions in explaining observed data, e.g. super-luminal jets from disks, luminosity, and then extracting fundamental properties from them. Hence any conclusion based on constant γ and β should be taken with caution and corrected.     

Nonextensive Statistical Mechanics: Some Links with Astronomical Phenomena

A variety of astronomical phenomena appear to not satisfy the ergodic hypothesis in the relevant stationary state, if any. As such, there is no reason for expecting the applicability of Boltzmann–Gibbs (BG) statistical mechanics. Some of these phenomena appear to follow, instead, nonextensive statistical mechanics. In the same manner that the BG formalism is based on the entropy S _BG=?k∑ _i p _i ln p _i, the nonextensive one is based on the form S _q=k(1 ?∑ _i p _i ^q)/(q? 1) (with S ₁=S _BG). The stationary states of the former are characterized by an exponential dependence on the energy, whereas those of the latter are characterized by an (asymptotic) power law. A brief review of this theory is given here, as well as of some of its applications, such as the solar neutrino problem, polytropic self-gravitating systems, galactic peculiar velocities, cosmic rays and some cosmological aspects. In addition to these, an analogy with the Keplerian elliptic orbits versus the Ptolemaic epicycles is developed, where we show that optimizing S _q with a few constraints is equivalent to optimizing S _BG with an infinite number of constraints.     

The electron-ion streaming instabilities driven by drift velocities of the order of electron thermal velocity in a nonmagnetized plasma

The kinetic Alfven waves are investigated using Maxwell-Boltzmann-Vlasov equation to evaluate the kinetic dispersion relation and growth/damping rate with magnetic field gradient, density gradient, temperature gradient and velocity gradient with inhomogeneous plasma. The effect of gradient terms is included in the analysis for both the regions k _⊥ ρ _i <1 and k _⊥ ρ _i >1, where k _⊥ is the perpendicular wave number and ρ _i is the ion gyroradius. This study elucidates a possible scenario to account for the particle acceleration and the wave dissipation in inhomogeneous plasmas. This model is able to explain many features observed in plasma sheet boundary layer as well as to evaluate the dispersion relation, growth rate, growth length and damping rate of kinetic Alfven wave. The applicability of this model is assumed for auroral acceleration region, plasma sheet boundary layer and cusp region.     

The calculated and observed diurnal variation of the ionosphere over Millstone Hill on 23–24 March 1970

A numerical model of current F-region theory is use to calculate the diurnal variation of the mid-latitude ionospheric F-region over Millstone Hill on 23–24 March 1970, during quiet geomagnetic conditions. From the solar EUV flux, the model calculates at each altitude and time step primary photoelectron spectra and ionization rates of various ion species. The photoelectron transport equation is solved for the secondary ionization rates, photoelectron spectra, and various airglow excitation rates. Five ion continuity equations that include the effects of transport by diffusion, magnetospheric-ionospheric plasma transport, electric fields, and neutral winds are solved for the ion composition and electron density. The electron and ion temperatures are also calculated using the heating rates determined from chemical reactions, photoelectron collisions, and magnetospheric-ionospheric energy transport. The calculations are performed for a diurnal cycle considering a stationary field tube co-rotating with the Earth; only the vertical plasma drift caused by electric fields perpendicular to the geomagnetic field line is allowed but not the horizontal drift. The boundary conditions used in the model are determined from the incoherent scatter radar measurements of T_e, T_i and O⁺ flux at 800km over Millstone Hill (Evans, 1971a). The component of the neutral thermospheric winds along the geomagnetic field has an important influence on the overall ionospheric structure. It is determined from a separate dynamic model of the neutral thermosphere, using incoherent scatter radar measurements.The calculated diurnal variation of the ionospheric structure agrees well with the values measured by the incoherent scatter radar when certain restrictions are placed on the solar EUV flux and model neutral atmospheric compositions. Namely, the solar EUV fluxes of Hinteregger (1970) are doubled and an atomic oxygen concentration of at least $\text{10}{}^{11}\text{cm}{}^3$ at 120 km is required for the neutral model atmosphere. Calculations also show that the topside thermal structure of the ionosphere is primarily maintained by a flow of heat from the magnetosphere and the night-time F2-region is maintained in part by neutral winds, diffusion, electric fields, and plasma flow from the magnetosphere. The problem of maintaining the calculated night-time ionosphere at the observed values is also discussed.     

Effect of ion collisions on TM models in cylindrical wave guide filled with hot plasma

The expression for damping coefficients (K _i) is derived and discussed numerically, for a cylindrical wave guide, filled with hot collisional and uniaxially magnetised plasma. It is observed that TM modes suffer a very high damping for high values of plasma frequency (w _pe/w = 10) and low values of ion collision frequency (v _i/v _e = 10^?3), where as for low values of plasma frequency (w _pe/w = 0.1) the damping is low. The damping also increases as the ion temperature increases.     

A model for the steady interaction of the solar wind with the Moon

The steady state interaction of the solar wind with the Moon is modeled as a uniform, magnetized, quasi-neutral, collisionless, hypersonic, and hyper-Alfvénic flow of an electronproton plasma past a perfectly ion absorbing, non-magnetized sphere. For the temperature of the electronsT much less than that of the ionsT _i , steady state equations are derived self-consistently from the Vlasov and Maxwell equations by taking advantage of the fact that the ion gyration ratio is small compared to the radius of the Moon, by employing an ordering which requires different scale lengths along the magnetic fieldB and center of mass velocity, and by expanding in a small parameter ? that measures the smallness of ?B terms compared to a dominant term retained. A partial numerical solution is presented and discussed for the limit in which ? is much less than β=(ion pressure/magnetic pressure). In addition, a simple technique is presented whereby the steady state equations can be approximately extended to cases in whichT?T _i for arbitrary ?/β.     

Ion temperature troughs in the equatorial topside ionosphere

The Retarding Potential Analyzer aboard OGO-6 sometimes recorded marked depressions of ion temperature as the satellite crossed the equatorial region. These “T_i troughs” occur at heights between about 700 km and the satellite apogee at 1100 km. At the centre of a trough, close to the dip equator, T_i is frequently 500–1000 K below its value at the northern and southern edges, which are usually 15°–20° in latitude from the centre of the trough. At a given season and local time, the occurrence, symmetry, depth and position of the troughs often vary markedly with longitude. The troughs have no particular association with equatorial troughs of ion concentration N_i.As suggested by Hanson, Nagy and Moffett, the T_i troughs appear to be caused by transequatorial winds that drive F region plasma along geomagnetic field lines. The plasma is adiabatically cooled as it is driven upwards on the “upwind” side of the dip equator, and heated as it descends on the “downwind” side. The available data on the occurrence of troughs at different longitudes, local times and seasons are reasonably consistent with wind directions deduced from Jacchia's model and the OGO-6 thermospheric model of Hedin et al., and with the north-south asymmetries of the tropical 630 nm airglow observed by OGO-4 and OGO-6. Factors determining the latitudinal extent of the troughs are discussed and some questions for further study are listed.     

Total current to cylindrical collectors in collisionless plasma flow

A theory is presented for charged-particle collection by a cylindrical conducting object, such as a spacecraft or an electrostatic probe, which is moving transversely through a collisionless plasma, such as those in the upper atmosphere and space. The calculation is approximate, using symmetric potential profiles which are exact for the infinite-cylinder stationary case. Theoretical current predictions are presented for ratios of collector potential to electron thermal energy eφ_c/kT_e from 0 to ?25, for ion-to-electron temperature ratios T_i/T_c = 1 and 0.5, ratio of collector radius to electron Debye length r_c/λ_D from 0 to 100, and ratio of flow speed to ion thermal speed $\text{S}{}_{\mathrm{i}}\text{= U/(2kT}{}_{\mathrm{i}}\text{/m}{}_{\mathrm{i}}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{)}$ from 0 to 10. Comparisons with existing exact calculations by other authors show that none of these fulfil all of the requirements for nontrivial comparison. Appropriate parameter ranges for future exact calculations are thereby suggested. These are as follows: (a) r_c/λ_D should be large enough that the collector not be in or near orbit-limited conditions; (b) the ratio $\text{S}{}_{\mathrm{i}}{}^2\text{/¦χ}{}_{\mathrm{c,\; i}}\text{¦}$ of ion directed energy to potential energy change in the sheath, should be close to unity or if

$\text{S}{}_{\mathrm{i}}{}^2\text{/¦χ}{}_{\mathrm{c,i}}\text{¦? 1,}\text{then}\text{S}{}_{\mathrm{i}}\text{? 1}$

.     

Effects of ion collisions on the interaction of transverse waves with half-space plasma

For half-space (Z>0), homogeneous, collisonal and warm plasma, the expressions for fields and penetration depth δ/δ _e (in the unit of ion collisionless cold plasma penetration depth, i.e., when v _i =0, υ_0i =) are derived and discussed numerically. It is concluded that the propagation of transverse waves is only slightly affected by the ion collisions and the applied magnetic field when the plasma frequence is greater than the wave frequency (ω _pe >ω). For the case of ω _pe ≤ω, the damping of the wave is not affected by the changes in the ion collision frequency and the ion temperature. However, in this case, the propagation of the wave is drastically affected by the applied magnetic field and the wave damps quickly as the magnetic field strength or the gyrofrequency (Ω _e ) increases.     

Dynamical instability of laminar axisymmetric flows of ideal fluid with stratification

The instability of nonhomentropic axisymmetric flows of ideal fluid with respect to two-dimensional infinitesimal perturbations with the nonconservation of angular momentum is investigated by numerically integrating the differential equations of hydrodynamics. This problem is important in studying the dynamics of astrophysical flows as shear fluid flows around a gravitating center. A complex influence of a nonzero entropy gradient on the instability of sonic and surface gravity modes has been found. In particular, both an increase and a decrease in entropy against the effective gravity g _eff causes the growth of surface gravity modes that are stable at the same parameters for a homentropic flow. At the same time, the growth rate of the sonic instability branches either monotonically increases with increasing rate of decrease in entropy against g _eff or becomes zero at both negative and positive entropy gradients in the unperturbed flow. Calculations also show that growing internal gravity modes appear in the problem with free boundaries under consideration only if the flow is no longer stable with respect to axisymmetric perturbations. In addition, we show that it is improper to specify the entropy distribution in the main flow by a polytropic law with a polytropic index different from the adiabatic value, since the perturbation field does not satisfy the boundary condition at a free boundary in this case.     

Korteweg-de Vries equation for ion acoustic soliton with negative ions in the presence of nonextensive electrons

Korteweg-de Vries (KdV) equation for electrostatic ion acoustic wave in a three component plasma containing positive and negative ions along with the nonextensive electrons is derived. Fast and slow ion acoustic modes which propagate with different velocities are excited. The effects of variation of quantities like q (nonextensive parameter), Q (mass ratio of positive to negative ion), μ (electron to positive ion number density ratio), θ _i (positive ion to electron temperature ratio) and θ _n (negative ion to electron temperature ratio) have been presented for fast and slow ion acoustic modes. Both compressive and rarefactive solitons are observed. It is found that the solitary excitations strongly depend on the mass and density ratios of the positive and negative ions as well as on nonextensive electron parameter.     

Influence of suprathermality on the obliquely propagating dust-acoustic solitary waves in a magnetized dusty plasma

Dust-acoustic (DA) solitary waves are investigated in a magnetized dusty plasma comprising cold dust fluid and kappa-distributed ions and/or electrons. The influence of suprathermal particles, obliqueness, and ion temperature on the DA solitary waves is investigated. We find that only negative DA solitary waves will be excited in this model. Also it is shown that the amplitude of the DA solitary wave decreases with deviation of electrons or ions from Maxwellian distribution via decrease of κ _e or κ _i . The effect of the temperature of the ion decreases with the amplitude and steepness of the solitary wave front.     

Nonplanar shocks in a warm electronegative plasma with electron nonextensivity effects

By employing the reductive perturbation technique, nonlinear cylindrical and spherical Korteweg–de Vries Burgers (KdVB) equation is derived for ion acoustic shock waves in an unmagnetized electronegative plasma. The latter is composed of warm positive and warm negative ions as well as q-distributed nonextensive electrons. Numerically, the modified KdVB equation is solved to examine the impact of nonthermal electrons on the profiles of nonplanar fast ion acoustic shocks. With the help of experimental parameters, it is found that the variations of different quantities, like q (nonextensive parameter), α (the negative-to-positive ion mass ratio), μ (the electron-to-positive ion density ratio) and θ _i (the positive ion-to-electron temperature ratio), η _i0,n0 (the positive/negative ion viscosities) significantly modify the propagation characteristics of nonplanar shocks in electronegative plasmas. The relevance to a laboratory experiment is highlighted, where positive and negative ions are present.     

Investigation of solitary waves in warm plasma for smaller order relativistic effects with variable pressures and inertia of electrons

In the two component relativistic plasmas subject to pressure variation of adiabatic electrons and isothermal ions, both compressive and rarefactive KdV solitons are established in a quite different physical plasma model. It is desirable to define c _s in a new way to substantiate the validity of the model under relativistic effects. The corresponding mathematical condition is also determined, which is a new report of this kind. It is also interesting to report that the relativistic rarefactive solitons cease to exist below some critical ion initial streaming speed v _i0 for a fixed temperature α and electron streaming speed v _e0. Besides, higher initial flux v _i0 of ions under constant temperature is observed to generate higher speed v _i at the passage of time which causes to increase (in relativistic sense) its mass diminishing thereby the growth of soliton amplitudes.     

Ionospheric warming by neutral winds

The effect of frictional heating by means of neutral winds on the ion and electron temperature in the undisturbed ionosphere is studied theoretically by solving a system of basic ionospheric and atmospheric equations. The study shows that both the electron and ion temperatures are increased in the night-time ionosphere through friction. In the region between 150 and 200 km T_i may exceed T₆ by as much as 130°. The increase of T_i due to friction amounts to about 100–200°, depending on the atmospheric model employed in calculating the neutral wind velocity. It is illustrated that frictional heating may be very important for the determination of the neutral temperature from measured ion temperature values.     

Monitoring of energy spectra of particle bursts in the plasma sheet and magnetosheath

Energetic ion (E ? 290 keV) and electron (E_e ? 220 keV) burst intensities were simultaneously monitored at various regions of the plasma sheet and magnetosheath by the CPME JHU/APL instruments on board the IMP-7 and 8 s/c during an extended period from day 250, 1975 to day 250, 1976 when the two spacecraft were closely trailing each other in crossing the geomagnetotail. The energy spectra of the energetic particle populations of different regions in the magnetotail were also computed and monitored simultaneously at the positions of the two spacecraft. The results indicate that the energetic particle intensities are higher and the energy spectra in general considerably softer inside the plasma sheet than the adjacent magnetosheath. The spectral index γ of a power law fit in the computed energy spectrum inside the plasma sheet occasionally exceeds γ > 10 for the ions and γ > 6 for the electrons. Furthermore simultaneous monitoring of particle intensities in the vicinity of the neutral sheet and the high latitude plasma sheet shows higher intensities in the former region. The observations suggest that the energetic particles escape to the magnetosheath from their source inside the plasma sheet by a rigidity dependent process. A dawn-dusk asymmetry in the particle acceleration and escape processes is implied in the observations and discussed in detail.     

Adiabatic indices in a convecting anisotropic plasma

Adiabatic indices for a non-dissipative anisotropic convecting plasma are analyzed, and general expressions for the effective adiabatic index and the partial adiabatic indices parallel (γ_∥) and perpendicular (γ_⟂) to the magnetic field are obtained. It is shown that, in the general case, the value of the effective adiabatic index is not an universal constant and depends on the plasma temperature anisotropy and on the properties of the plasma motion. The values of γ_⟂ and γ_∥ are shown to be independent of the plasma parameters being completely determined by the characteristics of the plasma flow.