激波聚焦诱导气液两相爆震燃烧的数值模拟

对以激波聚焦和增加障碍物方式诱导煤油-空气气液两相爆震燃烧的过程进行了数值模拟.采用欧拉-拉格朗日方法建立了脉冲爆震发动机(PDE)中气液两相流的喷射、雾化、掺混过程.研究发现环形爆震波在爆震管凹腔内经过反射、汇聚后能够引燃可燃混合物.而在障碍物处,激波的反射和再反射聚焦能够形成高温高压点(2700K,25MPa),产生局部爆炸,有助于形成稳定的脉冲爆震燃烧(波面速度为1900m/s,温度为2 950K),有效地缩短由缓燃向爆震转变(DDT)距离至0.45m.     

几种常用格式和湍流模型在高速流动中的适用性研究

针对高超声速飞行器中存在的强激波、激波/边界层干扰、分离、湍流等复杂流动现象,对比分析了当前计算流体力学中的主要空间离散格式及湍流模型,发现不同格式对强激波的分辨率基本相同,Roe和LDE(low diffusion E-CUSP(convective upwind split pressure))格式对摩擦因数和传热系数的模拟优于其他格式;S-A(Spalart-Allmaras)一方程湍流模型计算的摩擦因数比k-ωSST(shear stress transport)两方程湍流模型高10%左右,而后者预测的分离区约为前者的2倍,且分离点靠前.     

温度对液体中超声速度的影响

本文导出了超声波速度同液体压缩系数及密度关系，研究了压缩系数及密度与温度的关系，进而研究了温度对声速及声时的影响，用实验测量了不同的液体成份下的声时同温度的关系，从而提出了进行液体成份分析所必需的参量，对用超声波分析液体成份提供了理论依据，对提高产品质量及实现工业在线检测有一定的实用价值。     

相对论性离子束（环束）流引起的电磁不稳定性

在相对论性框架下，研究了冷离子流产生的右旋共振模，非共振模电磁不稳定性的特征结果表明，非共振模不稳定性增长率受离子流初始参数影响较大，而右旋共振模的增长率则受离子流初始参数影响较小。此结果可应用于太阳高能离子流通过低日冕区和地球磁层辐射带俘获的高能质子情形。     

Solar cycle and latitude dependence of high-beta suprathermal plasma conditions in interplanetary space between 1.3 and 5.4 AU

The analysis of energetic particles and magnetic field measurements from the Ulysses spacecraft has shown that in a series of events, the energy density contained in the suprathermal tail particle distribution is comparable to or larger than that of the magnetic field, creating conditions of high-beta plasma. In this work we analyze periods of high-beta suprathermal plasma occurrences (β_ep > 1) in interplanetary space, using the ratio (β_ep) of the energetic particle (20 keV to ∼5 MeV) and magnetic field energy densities from measurements covering the entire Ulysses mission lifetime (1990–2009) in order to reveal new or to reconfirm some recently defined interesting characteristics. The main key-results of the work are summarized as follows: (i) we verify that high-beta events are detected within well identified regions corresponding mainly to the vicinity of shock surfaces and magnetic structures, and associated with energetic particle intensity enhancements due to (a) reacceleration at shock-fronts and (b) unusually large magnetic field depressions. (ii) We define three considerable features for the high-beta events, concentrated on the next points: (a) there is an appreciable solar-activity influence on the high-beta events, during the maximum and middle solar-cycle phase, (b) the annual peak magnitude and the number of occurrences of high events are well correlated with the sunspot number, (c) the high-beta suprathermal plasma events present a spatial distribution in heliographic latitudes (HL) up to ∼±80°, and a specific important concentration on the low (−25° ? HL < −6°, 6° < HL ? 25°) and median (−45° ? HL < −25°, 25° < HL ? 45°) latitudes. We also reconfirm by a statistical analysis the results of Marhavilas and Sarris (2011), that the high-beta suprathermal plasma (β_ep > 1) events are characterized by a very large parameter β_ep (up to 1732.5), a great total duration (406 days) and a large percentage of the Ulysses-mission lifetime (which is equal to 6.34% of the total duration with usable measurements, and 11.3% of the duration in presence of suprathermal particles events).     

Plasma and electromagnetic wave simulations of meteors

Every day billions of meteoroids impact and disintegrate in the Earth’s atmosphere. Current estimates for this global meteor flux vary from 2000 to 200,000 tons per year, and estimates for the average velocity range between 10 km/s and 70 km/s. The basic properties of this global meteor flux, such as the average mass, velocity, and chemical composition remain poorly constrained. We believe much of the mystery surrounding the basic parameters of the interplanetary meteor flux exists for the following reason, the unknown sampling characteristics of different radar meteor observation techniques, which are used to derive or constrain most models. We believe this arises due to poorly understood radio scattering characteristics of the meteor plasma, especially in light of recent work showing that plasma turbulence and instability greatly influences meteor trail properties at every stage of evolution. We present our results on meteor plasmas simulations of head echoes using particle in cell (PIC) ions, which show that electric fields strongly influence early stage meteor plasma evolution, by accelerating ions away from the meteoroid body. We also present the results of finite difference time domain electromagnetic simulations (FDTD), which can calculate the radar cross section of the simulated meteor plasmas. These simulations have shown that the radar cross section depends in a complex manner on a number of parameters. These include the angle between radar and meteor entry, a large dependence on radar frequency, which shows that for a given meteor plasma size and density, the reflectivity as a function of probing radar frequency varies, but typically peaks below 100 MHz.     

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.     

Height and time variability of planetary wave activity

The height–season and year-to-year regularities of parameters of first and second spatial harmonics determine the structure of the stratosphere and mesosphere circulation and its variability. In the period 1992–2002 at heights 0–55 km, the amplitudes and phases of the first and second spatial harmonics in the field of temperature, geopotential height, zonal and meridional wind were calculated by the method of harmonic decomposition. Dispersion (standard or mean square deviation) of their day-to-day and year-to-year variations was calculated by their wavelength constants. Height and season patterns of variability have been estimated. The difference in height–longitude variability for wave numbers m = 1 and 2 has been discovered. At the same time, the intensity of wave disturbances for m = 1 is less than for m = 2 excluding the polar areas, where a significant variability appears at the heights 0–55 km. There is also a tendency for the intensity of year-to-year variations to decrease in comparison with day-to-day variations. In cold and warm periods the amplitude of perturbation waves with m = 2 both for day-to-day and year-to-year variations is greater than for waves with m = 1. Transient height areas in the interval of 20–30 km are more distinct for day-to-day variations of polar area.     

The role of radial perpendicular diffusion and latitude dependent acceleration along the solar wind termination shock

A numerical model, based on Parker’s transport equation, describing the modulation of anomalous cosmic rays and containing diffusive shock acceleration is applied. The role of radial perpendicular diffusion at the solar wind termination shock, and as the dominant diffusion coefficient in the outer heliosphere, is studied, in particular the role it plays in the effectiveness of the acceleration of anomalous protons and helium when its latitude dependence is changed. It is found that the latitudinal enhancement of radial perpendicular diffusion towards the heliospheric poles and along the termination shock has a prominent effect on the acceleration of these particles. It results in a ‘break’ in the energy spectrum for anomalous protons at ∼6.0 MeV, causing the spectral index to change from E^−1.38 to E^−2.23, but for anomalous helium at ∼3.0 MeV, changing the spectral index from E^−1.38 to E^−2.30. When approaching the simulated TS, the changes in the modulated spectra as they unfold to a ‘steady’ power law shape at energies below 50 MeV are much less prominent as a function of radial distances when radial perpendicular diffusion is increased with heliolatitude.     

Kinetic Excitation Mechanisms for ION-Cyclotron Kinetic Alfvén Waves in Sun-Earth ConnectionI

We study kinetic excitation mechanisms for high-frequency dispersive Alfvén waves in the solar corona, solar wind, and Earth's magnetosphere. The ion-cyclotron and Cherenkov kinetic effects are important for these waves which we call the ion-cyclotron kinetic Alfvén waves (ICKAWs). Ion beams, anisotropic particles distributions and currents provide free energy for the excitation of ICKAWs in space plasmas. As particular examples we consider ICKAW instabilities in the coronal magnetic reconnection events, in the fast solar wind, and in the Earth's magnetopause. Energy conversion and transport initiated by ICKAW instabilities is significant for the whole dynamics of Sun-Earth connection chain, and observations of ICKAW activity could provide a diagnostic/predictive tool in the space environment research. This revised version was published online in August 2006 with corrections to the Cover Date.