Surface tension-driven flows in horizontal liquid metal layers

The problem of surface tension-driven flows in horizontal liquid layers has been studied experimentally, and theoretically by direct numerical simulation and small perturbation analysis. We focus our attention on situations in which the depth of the fluid (liquid tin; small Prandtl number, Pr=0.015) is small enough to ensure the predominance of the surface tension forces over those due to the buoyancy. The surface velocity has been experimentally obtained for liquid tin layer with various aspect ratio (length to height) in the range 5<A<83. The thermal gradients are ranged from 5 to 40°K/cm. In the numerical study, the Navier-Stokes and energy equations are solved by an efficient finite difference technique. The parameters governing the flow behaviour in the liquid are varied to determine their effects on thermocapillary convection: the Reynolds number 10<Re<210⁴ and the aspect ratio 2<A<25; with Pr kept constant at Pr=0.015. The linear eigenequation resulting from small spatial disturbances of the Couette flow solution is solved using an Tau-Chebyshev approximation. A notable feature of the theoretical study is the totally different end circulations. In the region near the cold wall a multicell structure is evident. This agrees with the eigensolution which is of complex type, indicating spatial periodicity. In the hot wall region the flow is accelerated to reach the velocity value for the fully-developed Couette flow which is reached under conditions such as Re/A<20. The transition from viscous to boundary layer regime occurs for a critical value (Re/A)_c of nearly about 200, as deduced from the numerical and experimental results.     

两种用于NS-DSMC耦合方法的连续失效参数对比

NS-DSMC（Navier Stokes-Direct Simulation Monte Carlo）耦合方法是计算连续-稀薄跨流域流动的主要方法,应用过程中如何确定连续流域和稀薄流域的界面是此方法的关键问题之一,界面位置通常通过连续失效参数来判定.为合理选择连续失效参数,对目前广泛使用的两种连续失效参数Kn_Q和B参数进行了理论上的分析和比较,表明虽然他们建立的出发点不同,但在数学形式上具有一定相似性.通过圆柱绕流问题的数值试验进一步验证了两种连续失效参数具有一定对应关系.Kn_Q和B两种连续失效参数在NS-DSMC耦合方法中应用效果相差不大,但Kn_Q 的阈值0.05适用范围较广,而B参数的阈值对不同流动问题会有变化.      

基于分数阶Maxwell模型的振荡流传热特性分析

经求解基于分数阶Maxwell模型的粘弹性流体在周期性振荡压力梯度下,圆直管内的运动方程和能量方程,得到了振荡管流换热的速度分布、温度分布以及热扩散系数的解析解形式.通过对无量纲热扩散系数的分析可知,影响粘弹性流体管内振荡流轴向换热的无量纲参数有:Womersley数Wo、Deborah数De、无量纲振幅Δx/R和流体普朗特数Pr.分数阶Maxwell模型振荡流传热也存在粘弹性流体流动中存在的共振现象,且共振峰的数量随De数的减小而增加,发生共振的起始频率随De数的减小而降低.共振峰值出现的位置即频率值与Pr和无量纲振幅Δx/R无关.     

INTERBALL statistical study of ion flow fluctuations in the plasma sheet

We use ion distribution measurements with CORALL instrument on-board the INTERBALL/Tail spacecraft to study plasma flows in the mid-tail (−9> X> −27 R_E) plasma sheet. Three velocity components computed every 2 minutes exhibit two types of velocity variations: Earthward bursty bulk flows (BBFs) and random flow fluctuations. Their properties are in a good agreement with the observations of the ISEE-2 spacecraft (Borovsky et al., 1997). The INTERBALL/Tail spacecraft configuration favors measurements of V_z component, in contrast to previous experiments in which only V_x and V_y were measured reliably. In the outer part of the plasma sheet V_y and V_z fluctuations were close to each other (variances σ(V_y) and σ(V_z) were about 160 and 110 km/s, respectively), but in the inner part at the dusk flank amplitude of V_y fluctuations increased and was 2 times higher than that of V_z component. This asymmetry of fluctuations should be taken into account during modern theoretical analysis and simulations.     

Calcium abundance measurements using the Yohkoh BCS

Soft X-ray observations by SMM and other spacecraft have shown that the abundance of certain elements in solar corona varies from flare to flare. In this study, observations made by the Yohkoh Bragg Crystal Spectrometer (BCS) in helium-like Ca XIX have been analysed, and Ca abundance determined for 177 flares observed during the first four years of the mission (1991–1995). The average abundance of Ca relative to H for all flares is A_Ca = (3.64±0.39) × 10⁻⁶. As with an earlier study of SMM data, the abundance is found to be enhanced compared to the photosphere ((2.24±0.10) × 10⁻⁶), and with only minor variation from flare to flare. However, the absolute value and range of values determined by this study is smaller than in the previous study; these differences are discussed.     

Auroral signatures of transient processes in the outer magnetosphere

We present an analysis of sporadic and recurrent injections of magnetospheric ions in the midnight auroral oval during substorms and of the associated ionospheric ion outflows. The source of plasma sheet precipitating ions is determined using a simple method, based on the measured relation between the ion inverse velocity and time (l = v × t). This method is applied here to two typical passes of the Interball-Auroral (IA) satellite at distances of 3 R_E above the auroral regions. Substorm related ion injections are shown to be mainly due to time of flight effects. In contrast with particle trajectory computations (Sauvaud et al., 1999), the inverse velocity method does not require magnetic and electric field models and can thus be used systematically for the detection of time of flight dispersed ion structures (TDIS). This allowed us to build a large database of TDIS events and to perform a statistical analysis of their spatial distribution. For the cases presented here the source region of the injected ions is found at radial distances from 18 to 30 R_E near the equatorial magnetosphere. At Interball altitudes ( 3 R_E), ion injections detected at the poleward boundary of the nighside auroral oval are associated with shear Alfvén waves superimposed over large-scale quasi-static current structures. We show that the most poleward TDIS are collocated with a large outflow of ionospheric H⁺ and O⁺ displaying pitch-angle distributions peaked in the pitch-angle range 90°–120°. These ions are thus accelerated perpendicularly to the magnetic field not only in the main auroral acceleration region but also up to at least 3 R_E. The expanding auroral bulge thus constitutes a significant source of H⁺ and O⁺ ions for the mid-tail magnetosphere.     

Plasma convection in the low-latitude mantle of the near-earth magnetosphere

A detailed analysis of Geotail observations on June 1, 1997 provides evidence for the appearance of the plasma mantle at low-latitudes, Z_GSM 3 R_E for southward interplanetary magnetic field (IMF), and clarifies the characteristics of the plasma flow in the low-latitude mantle. The tailward flow speed in the mantle is 80 160 km/s, smaller than the flow speed (300 km/s) in the adjacent magnetosheath. The component of the velocity perpendicular to the magnetic field lines is predominantly poleward at speeds up to 60 km/s, and the convection tends to be stronger for |;B_Y|<|;B_Z| than for IMF |B_Y|>|B_Z|;. This polewardconvection can be interpreted as being driven by a J × B force operating on the field lines.     

On the adjust of minor constituents density in middle atmosphere during solar proton events

Profiles of O₃ partial pressure and of other minor atmospheric constituents (NO, NO₂, HCL, HF and H₂O), observed in the middle atmosphere during Solar proton events (20.04.1998; 05.04.2000), were analysed. Conclusions were drawn that under SCR impact a short-term O₃ partial pressure increase and destruction of some freon constituents took place.     

Inconsistency of magnetic field and plasma velocity variations in the distant plasma sheet: Violation of the “frozen-in” criterion?

Measurements of the magnetic field and low energy plasma by the GEOTAIL spacecraft have been used to study the relationship between variations of the plasma velocity and of the magnetic field in the distant (100–200 R_E) and middle (40–80 R_E) tail. The analysis was carried out separately for the tail lobes and the plasma sheet. It is shown that the absolute values of the magnetic field and plasma velocity, as well as their corresponding components (V_X and B_X, V_Y and B_Y, V_Z and B_Z), are linearly connected in the tail lobes. In the plasma sheet, however, the plasma velocity and the magnetic field do not seem to be related to one another. The distant plasma sheet seems to be in a regime of turbulence. The diffusion coefficients estimated from our data set of the velocity parameters in the plasma sheet are in good agreement with the theoretical predictions of Antonova and Ovchinnikov (1996, 1999).     

Intensity profiles of the multiple scattering light near the cometary nucleus

The multiple scattering of solar radiation in the cometary atmosphere is treated with the method of successive scattering. Referring to in situ measurements of comet Halley about the size and spatial distributions of dust, the optical thickness τ₁ of dust has been estimated, i.e. τ₁=0.03 at wavelength λ=0.62μm in a quiet time, but τ₁=0.3 when the outbursts/jets occur. In the derivation of τ₁, optical properties of dust including a mixing ratio of absorbing to silicate grains, are determined based on the polarimetry of P/Halley at λ=0.62μm observed during the phase angles over Nov. 1985 to May 1986 at the Dodaira Station of Tokyo Astronomical Observatory.

It is found that a temporary enhancement of τ₁ leads an increase of the upward reflected intensity when the surface albedo A of the nucleus is less than 0.04, but the reverse is true when A>0.04. On the other hand, the intensity of the downward radiation at the surface of the nucleus always decreases as an increase of τ₁.     

Convective motions in liquid metals for material processings. Numerical simulation of oscillation regimes and effect of rotation

The paper is concerned with the numerical simulation and the analysis of some kinds of flow regimes which can develop in Bridgman and Czochralski systems for material processings. The flows in the liquid phase are investigated considering two-dimensional and axisymmetric models. The time-dependent regimes were studied for a zero-Prandtl-number fluid layer confined inside a two-dimensional cavity of aspect ratio (length-to-height) A=4, involving a stress-free upper surface and submitted to a horizontal temperature gradient. The range of Grashof number was varied up to the conditions at which the flow goes from oscillatory to chaotic type behaviours. The combined influence of the temperature gradients and of the rotations of the crucible and of the seed/crystal was investigated for a Czochralski model. The axisymmetric regimes were studied for a Pr_m=0.015 liquid melt confined inside a cylindrical crucible of aspect ratio (height-to-radius) A_m=2, and coupled to a viscous encapsulant liquid layer (10<Pr_e<1200) of aspect ratio A_e=0.5. A number of steady and (transient) time-dependent flow patterns are identified.     

Numerical simulation of three-dimensional reconnection due to the instability of collisionless current sheets

Based on analytical calculations we have currently argued that spontaneous reconnection through thin collisionless current sheets is an essentially three-dimensional (3 D) process (Büchner, 1996 a, b). Since 3 D kinetic PIC codes have become available, the three dimensional nature of the collisionless current sheet decay are now illustrated by numerical simulations (Büchner and Kuska, 1996; Pritchett and Coroniti, 1996; Zhu and Winglee, 1996). While the latter two claim a coupling to a longer wavelength kink mode as a main factor, destabilizing thin current sheets in 3 D, our simulations have revealed that even shorter scale perturbations in the current direction suffice to destabilize thin sheets very quickly. Since past simulation runs, however, were limited to mass ratios near unity, the influence of the electrons was not treated adequately. We have now investigated the stability of thin collisionless current sheets including 64 times lighter negatively charged particles. We can now show that while the two-dimensional tearing instability slows down for M = M_p/m_e = 64, the three-dimensional current sheet decay is a much faster process — practically as fast as the mass ratio M = 1 3 D sheet decay, even without kinking the sheet. We further conclude that, unlike the two-dimensional tearing instability, the three-dimensional decay of thin current sheets is not controlled by the electrons. For a sheet width comparable with the ion inertial length, we also recovered signatures of the Hall effect as predicted by Vasyliunas (1975) in the mass ratio M = 64 case. The ion inertial length seems to be the critical scale at which the sheet starts to decay.     

Collisionless reconnection in the magnetotail

Recent developments regarding collisionless reconnection in current sheets with a finite normal magnetic field component (B_z) are reviewed. In 2-D x, z configurations the ion tearing mode is stabilized by the electron compressibility. When the y dependence is included, cross-field current instabilities can be excited. Of these, the drift kink mode appears to be particularly important. 3-D electromagnetic particle simulations indicate that this mode can act as the precursor to the growth of tearing modes and subsequent reconnection.     

LMXB luminosity function and connection to globular cluster in early type galaxies

We derive bias-corrected X-ray luminosity functions (XLFs) of LMXBs detected in 14 E and S0 galaxies observed with Chandra. After correcting for incompleteness, the individual XLFs are statistically consistent with a single power-law. A break at or near L_X,Eddington , as previously reported, is not required in any individual case. The combined XLF with a reduced error, however, suggests a possible break at L_X = 5 × 10³⁸ erg s⁻¹, which may be consistent with the Eddington luminosity of neutron stars with the largest possible mass (3 M), or of He-enriched neutron star binaries. We confirm that the total X-ray luminosity of LMXBs is correlated with the the near-IR luminosities, but the scatter exceeds that expected from measurement errors. The scatter in L_X(LMXB)/L_K appears to be correlated with the specific frequency of globular clusters, as reported earlier.

We cross-correlate X-ray binaries with globular clusters determined by ground-based optical and HST observations in 6 giant elliptical galaxies. With the largest sample reported so far (300 GC LMXBs with a 5:2 ratio between red and blue GCs), we compare their X-ray properties, such as X-ray hardness, XLF and L_X/L_B and find no statistically significance difference between different groups of LMXBs. Regardless of their association with GCs, both GC and field LMXBs appear to follow the radial profile of the optical halo light, rather than that of more extended GCs. This suggests that while metallicity is a primary factor in the formation of LMXBs in GCs, there may be a secondary factor (e.g., encounter rate) playing a non-negligible role.     

IMF effect on ionospheric trough occurrence at equinoxes

Previous observations have shown that there is a relationship between the F region trough and both B_z and B_y components of the interplanetary magnetic field (IMF). Since IMF governs the polar cap convection, we investigate here if this relationship can be explained by means of polar cap convection. The study is limited to equinox seasons. The poleward and equatorward edges of the trough are determined from satellite tomographic observations and their locations are plotted in magnetic coordinates together with the convection pattern given by Papitashvili and Rich [Papitashvili, V.O., Rich, F.J. High-latitude ionospheric convection models derived from DMSP ion drift observations and parameterized by the IMF strength and direction. J. Geophys. Res. 107, 2002, doi:10.1029/2001JA000264] using IMF measurements coincident with trough observations. The results indicate a close relationship between the troughs and convection. Most of the troughs are seen within the dusk cell and the pattern of trough observations rotates with the convection pattern, when B_y changes its sign. More dayside troughs are observed when B_z is negative than in the opposite case, i.e. fast convective flow favours the dayside trough occurrence. Nightside troughs are observed more frequently when B_y is negative. In both evening and morning sectors the trough is situated close to the edges of convection cells, which partly contradicts previous results showing that the troughs are associated with the convection reversal. It is concluded that plasma convection has an important role in trough generation, although the effect of a strong electric field and other mechanisms like precipitation certainly have a role of their own.     

The collective gyration of a heavy ion cloud in a magnetized plasma

In both the ionospheric barium injection experiments CRIT I and CRIT II, a long-duration oscillation was seen with a frequency close to the gyro frequency of barium and a time duration of about one second. A model for the phenomenon which was proposed for the CRIT 1 experiment is here compared to the results from CRIT II which made a much more complete set of measurements. The model follows the motion of a low-β ion cloud through a larger ambient plasma. The internal field of the model is close to antiparallel to the injection direction v_i but slightly tilted towards the self-polarization direction E_P = −v_ixB. As the ions move across the magnetic field, the space charge is continuously neutralized by magnetic-field aligned electron currents from the ambient ionosphere, drawn by the divergence in the perpendicular electric field. These currents give a perturbation of the magnetic field related to the electric field perturbation by ΔE/ΔB ≈ V_a. The model predictions agree quite well with the observed vector directions, field strengths, and decay times of the electric and magnetic fields in CRIT II. The possibility to extend the model to the active region, where the ions are produced in this type of self-ionizing injection experiments, is discussed.     

PANS方法在双圆柱绕流数值模拟中的性能分析

PANS方法以传统的雷诺平均Navier-Stokes(RANS)方程为模板,引入控制参数对原方程进行修改,是一类能够有效模拟分离湍流流动的混合方法.采用基于Menter 剪切应力输运(SST)湍流模型的PANS方法,选取双圆柱绕流作为研究对象, 分析双圆柱流动结构,验证方法的准确性与可靠性,考察PANS方法中模型参数选取对于计算结果的影响,评估不同类型PANS方法的性能.研究表明:模化湍动能比例这一模型参数对计算结果有显著影响,且计算域内统一的模型参数取值难以处理复杂流动问题.在此基础上使方法中模化湍动能比例取值随流场物理信息和网格尺度变化,所得结果与SST DES计算结果及实验结果吻合较好,表明该方法适用于复杂湍流数值模拟.      

Large geomagnetic storms of extreme solar event periods in solar cycle 23

During extreme solar events such as big flares or/and energetic coronal mass ejections (CMEs) high energy particles are accelerated by the shocks formed in front of fast interplanetary coronal mass ejections (ICMEs). The ICMEs (and their sheaths) also give rise to large geomagnetic storms which have significant effects on the Earth’s environment and human life. Around 14 solar cosmic ray ground level enhancement (GLE) events in solar cycle 23 we examined the cosmic ray variation, solar wind speed, ions density, interplanetary magnetic field, and geomagnetic disturbance storm time index (D_st). We found that all but one of GLEs are always followed by a geomagnetic storm with D_st  −50 nT within 1–5 days later. Most(10/14) geomagnetic storms have D_st index  −100  nT therefore generally belong to strong geomagnetic storms. This suggests that GLE event prediction of geomagnetic storms is 93% for moderate storms and 71% for large storms when geomagnetic storms preceded by GLEs. All D_st depressions are associated with cosmic ray decreases which occur nearly simultaneously with geomagnetic storms. We also investigated the interplanetary plasma features. Most geomagnetic storm correspond significant periods of southward B_z and in close to 80% of the cases that the B_z was first northward then turning southward after storm sudden commencement (SSC). Plasma flow speed, ion number density and interplanetary plasma temperature near 1 AU also have a peak at interplanetary shock arrival. Solar cause and energetic particle signatures of large geomagnetic storms and a possible prediction scheme are discussed.     

Effect of HF emission of the topside sounder transmitter aboard the COSMOS-1809 satellite on the ionospheric plasma

The experiment on investigation of effect of the HF emission (300 W) by the dipole antenna on the ionospheric plasma was carried out onboard the COSMOS-1809 satellite (1987). The sounder accelerated particles (SAP) at the electron cyclotron harmonics n · ω_cs and in the frequency region of antenna resonance were detected by the charged particle spectrometer.     

From solar nanoflares to stellar giant flares: Scaling laws and non-implications for coronal heating

In this study we explore physical scaling laws applied to solar nanoflares, microflares, and large flares, as well as to stellar giant flares. Solar flare phenomena exhibit a fractal volume scaling, V(L)  L^1.9, with L being the flare loop length scale, which explains the observed correlation between the total emission measure EM_p and flare peak temperature T_p in both solar and stellar flares. However, the detected stellar flares have higher emission measures EM_p than solar flares at the same flare peak temperature T_p, which can be explained by a higher electron density that is caused by shorter heating scale height ratios s_H/L ≈ 0.04–0.1. Using these scaling laws we calculate the total radiated flare energies E_X and thermal flare energies E_T and find that the total counts C are a good proxy for both parameters. Comparing the energies of solar and stellar flares we find that even the smallest observed stellar flares exceed the largest solar flares, and thus their observed frequency distributions are hypothetically affected by an upper cutoff caused by the maximum active region size limit. The powerlaw slopes fitted near the upper cutoff can then not reliably be extrapolated to the microflare regime to evaluate their contribution to coronal heating.