Self-similar solutions to stellar winds with cosmic rays

A two-fluid model is used to study the time evolution of stellar winds including the dynamical effect of cosmic rays. Neglecting the diffusion of cosmic rays, we seek self-similar solutions to spherically symmetric winds with a termination shock. The velocity upstream of the shock is taken to be zero. Physical solutions are those that can connect the shock to the star with the velocity approaching zero at the star. Two parameters govern the behaviour of the solutions, namely, the ratio of the upstream sound speed to the shock speed (in an inertial frame) and the gravitation potential of the central star. In some parameter regimes, no physical solution is possible.     

Investigation of natural and artificial stimulation of the ionospheric Alfvén resonator at high latitude

A brief review is provided of recent progress in understanding the ionospheric Alfvén resonator (IAR) at high latitude. Firstly, naturally occurring resonances of the IAR as detected by pulsation magnetometers in the auroral zone at Sodankylä and in the polar cap at Barentsburg are considered. The characteristics of the IAR in the two regions are broadly similar, although the effects of solar illumination are less clear at the higher latitudes. Secondly we review recent attempts to stimulate the IAR through high-power radio frequency experiments both in the auroral zone at Tromsø with the European Incoherent SCATter (EISCAT) heater, and within the polar cap at Longyearbyen with the Space Plasma Exploration by Active Radar (SPEAR) facility. In the auroral zone at, Tromsø the stimulated IAR has been observed by ground-based magnetometers, and through electron acceleration observed on the FAST spacecraft. At SPEAR in the polar cap, the stimulated IAR has been investigated, with ground magnetometers, with the first results indicative of a positive detection.     

Shear Flow Instabilities in Low-Beta Space Plasmas

We study instabilities driven by a sheared plasma flow in the low-frequency domain. Two unstable branches are found: the ion-sound mode and the kinetic Alfvén mode. Both instabilities are aperiodic. The ion-sound instability does not depend on the plasma β (gas/magnetic pressure ratio) and has a maximum growth rate of about 0.1 of the velocity gradient dV ₀/dx. On the other hand, the kinetic Alfvén instability is stronger for larger β and dominates the ion-sound instability for β > 0.05. Possible applications for space plasmas are shortly discussed.     

Shock Wave Interaction with the Magnetopause

The magnetopause is in continuous motion and shock waves and impulsive acceleration events can occur. As an example, we show that the interaction of an interplanetary shock with the bow shock can generate a shock wave that after passing through the magnetosheath can interact with the magnetopause. In fluid dynamics, when a shock wave encounters a fluid discontinuity, the interface may become unstable and form bubbles and spikes. We consider this Richtmyer-Meshkov instability in magnetohydrodynamics. At the dayside magnetopause, the instability tends to be stabilized by the magnetic field. However, the shock wave interaction can initiate magnetic field reconnection for the southward IMF, which may be important in strong interplanetary shock events. This revised version was published online in August 2006 with corrections to the Cover Date.     

Simulating planetary wave propagation to the upper atmosphere during stratospheric warming events at different mountain wave scenarios

Parameterization schemes of atmospheric normal modes (NMs) and orographic gravity waves (OGWs) have been implemented into the mechanistic Middle and Upper Atmosphere Model (MUAM) simulating atmospheric general circulation. Based on the 12-members ensemble of runs with the MUAM, a composite of the stratospheric warming (SW) has been constructed using the UK Met Office data as the lower boundary conditions. The simulation results show that OGW amplitudes increase at altitudes above 30 km in the Northern Hemisphere after the SW event. At altitudes of about 50 km, OGWs have largest amplitudes over North American and European mountain systems before and during the composite SW, and over Himalayas after the SW. Simulations demonstrate substantial (up to 50–70%) variations of amplitudes of stationary planetary waves (PWs) during and after the SW in the mesosphere-lower thermosphere of the Northern Hemisphere. Westward travelling NMs have amplitude maxima not only in the Northern, but also in the Southern Hemisphere, where these modes have waveguides in the middle and upper atmosphere. Simulated variations of PW and NM amplitudes correspond to changes in the mean zonal wind, EP-fluxes and wave refractive index at different phases of the composite SW events. Inclusion of the parameterization of OGW effects leads to decreases in amplitudes (up to 15%) of almost all SPWs before and after the SW event and their increase (up to 40–60%) after the SW in the stratosphere and mesosphere at middle and high northern latitudes. It is suggested that observed changes in NM amplitudes in the Southern Hemisphere during SW could be caused by divergence of increased southward EP-flux. This EP-flux increases due to OGW drag before SW and extends into the Southern Hemisphere.     

一种有效的消波装置

对无限长渠道中小振幅入射波通过两块间距为ｌ,透水性能不同的多孔透水板的反射和透射进行了研究。应用分离变量法得到了反射波和透射波的波形。结果表明,当ｌ为入射波半波长的倍数时,反射波和透射波振幅最大,消波效果差,最佳消波率为５０％。当ｌ为入射波１／４波长加半波长的倍数时,反射波和透射波振幅最小,消波效果好,最佳消波率大于５０％。此外,还得到了无反射波的条件。     

递推阻尼最小二乘及其在INS/双星组合中应用

由于双星系统的误差模型未知,且稳定性较差,采用Kalman滤波实现惯导(INS)/双星组合导航时,其滤波特性较差;而此时最小二乘法是实际工程应用中常用的方法,但易产生参数爆发现象.因此,提出了一种基于多输入-多输出(MI-MO)系统的递推阻尼最小二乘法,并将其应用到INS/双星组合导航系统中.仿真结果表明,利用递推阻尼最小二乘法能达到和递推最小二乘法同样的精度,且能明显抑制参数爆发现象.      

A self-similar solution of a curved shock wave and its time-dependent force variation for a starting flat plate airfoil in supersonic flow

The problem of aeroelasticity and maneuvering of command surface and gust wing interaction involves a starting flow period which can be seen as the flow of an airfoil attaining suddenly an angle of attack. In the linear or nonlinear case, compressive Mach or shock waves are generated on the windward side and expansive Mach or rarefaction waves are generated on the leeward side. On each side, these waves are composed of an oblique steady state wave, a vertically-moving one-dimensional unsteady wave, and a secondary wave resulting from the interaction between the steady and unsteady ones. An analytical solution in the secondary wave has been obtained by Heaslet and Lomax in the linear case, and this linear solution has been borrowed to give an approximate solution by Bai and Wu for the nonlinear case. The structure of the secondary shock wave and the appearance of various force stages are two issues not yet considered in previous studies and has been studied in the present paper. A self-similar solution is obtained for the secondary shock wave, and the reason to have an initial force plateau as observed numerically is identified. Moreover, six theoretical characteristic time scales for pressure load variation are determined which explain the slope changes of the time-dependent force curve.     

Kink mode m = 1 in magnetic tube with discontinuous magnetic field

In this paper, we study conditions of realization and stability of kink modes with azimuthal wave numbers $m=\pm 1$ in cylindrical plasma flex with twisted magnetic field and homogeneous current along its axis. We assume permanent axial magnetic field both inside and outside the flex, surrounded by currentless plasma. Azimuthal magnetic field decreases inversely proportional to the distance from the boundary beyond the flex. We derived dispersion equations for stable and unstable modes in approximation of “thin” plasma flex. The analysis of equations has been provided for the case of discontinuous axial magnetic field on flex’s boundary. Conditions of propagation of wave modes have been defined. It was shown, that unstable modes can be implemented in certain interval of longitudinal wavenumbers. Results can be applied for the interpretation of solar magnetic tubes behavior, using measurements, provided by spacecrafts.     

斜激波总压损失率极小值理论解与物理意义

以法向马赫数作为激波强度表征量,对斜激波关系式进行重新推导,得到了穿过激波总压损失率极小值的理论解。控制方程表达式为激波角对物面角的线性函数。依据斜激波总压损失率极小值解析公式,首先,绘制了针对超声速流总压损失率应用的楔形角-激波角-马赫数的斜激波效率图。其次,通过生成斜激波三维总压损失率等值线图,呈现了总压损失率在楔形角-激波角-特征马赫数空间上的分布规律。此外,利用斜激波效率图,揭示了等总压损失率条件下马赫数与激波角的对称双解现象。