发动机进气道水滴撞击特性分析

针对某发动机进气道的防冰系统设计,对进气道的水滴撞击特性进行了研究.在对进气道流场进行计算的基础上,采用差分法对水滴运动方程进行了数值求解,得到了水滴运动轨迹,从而确定了水滴的撞击极限、总收集系数和局部水收集系数等水滴撞击特性参数;此外,还研究了飞行高度、飞行速度及水滴半径对水滴撞击特性的影响;研究发现水滴撞击极限、总收集系数和局部水收集系数随飞行速度的减小、飞行高度的增加或者水滴半径的增加而增大.这些为进气道的防冰系统设计奠定了基础.     

超声行波驱动的玻璃表面液滴运动数值模拟

针对液滴铺展和移动的动力学行为在工业生产和微流控芯片等领域有着重要作用,提出了一种基于超声行波理论的弹性体平板驱动模型,利用压电陶瓷的逆压电效应在弹性体玻璃产生超声行波从而驱动液滴运动。借助多物理场软件COMSOL建立液滴模型,首先进行行波分析,验证了驱动液滴的可行性。在0~60 ms中,液滴在超声行波的驱动下进行收缩-铺展的正弦振荡运动。然后通过液滴内部流场结构分析发现,当液滴半径铺展到最大后开始收缩时,液滴与基底接触面处的速度首先发生变化,表明液滴内部速度场的变化对接触线是否发生移动有着重要作用。液滴内部流场存在一个类似于椭圆形的漩涡,说明液滴运动不是单纯由于收缩-铺展而引起的平动,而是滚动着朝前运动。最后分别探讨了液滴移动速度与驱动电压、驱动频率以及动力黏度的关系,结果表明液滴移动速度受动力黏度影响较为显著。     

三相线钉扎的大液滴蒸发实验

利用实践十号返回式科学实验卫星蒸发对流箱,开展了三相线处于钉扎状态且接触半径大于毛细长度的无水乙醇大滴在加热PTFE表面蒸发的地基科学实验.实验发现,液滴体积随时间线性递减,但钉扎大液滴蒸发过程中没有出现恒定接触角(CCA)阶段.与小液滴蒸发的恒定接触半径(CCR)阶段相同,大液滴的平均蒸发速率也与初始体积无关,表明受...     

基于红外热像仪对珠状凝结换热过程实验

通过非接触式测温装置红外热像仪,观测冷凝液滴的形成、生长以及液滴之间汇聚、滚落过程.实验发现在珠状冷凝过程中,其表面温度变化显著.随着冷凝过程的进行,小液滴逐渐汇聚成大液滴,其液滴单位面积换热量明显下降,液滴表面温度逐渐降低.微纳米级小液滴在珠状冷凝换热中具有至关重要的作用.液滴在洁净表面上的珠状凝结和液滴在汇合后自清理表面上的珠状凝结表现出不同的温度分布,即珠状凝结汇聚后的自清理表面一直被前驱膜所覆盖.      

蒸气在铜表面冷凝的可视化实验

通过建立蒸气冷凝过程的可视化实验装置,在地面重力条件下对蒸气在铜表面的冷凝情况进行实验.采用高速摄像仪对蒸气在铜表面上发生的冷凝过程进行观测,并对比分析了H₂O,R141b以及FC-72蒸气的冷凝情况.研究发现水蒸气在冷凝过程中一些大液滴在合并周围小液滴后会发生回弹现象,从而清理掉周边区域的小液滴.随着冷凝的进行,被清理掉小液滴的区域又会重新生成小液滴.这种回弹现象相当于延长了滴状冷凝的时间,从而起到强化冷凝换热效果.R141b及FC-72在光滑铜表面上均未发生滴状冷凝.在加大过冷度后,冷凝器内残留的空气还会在冷凝液膜上冷凝成小水滴,产生滴膜共存现象.水蒸气在铜表面冷凝过程中,滴状冷凝转换成膜状冷凝后热流量降低约20%,因此需要对铜表面进行改性处理,降低其表面能,从而强化冷凝换热.      

Time dependent charging of layer clouds in the global electric circuit

There is much observational data consistent with the hypothesis that the ionosphere-earth current density (J_z) in the global electric circuit, which is modulated by both solar activity and thunderstorm activity, affects atmospheric dynamics and cloud cover. One candidate mechanism involves J_z causing the accumulation of space charge on droplets and aerosol particles, that affects the rate of scavenging of the latter, notably those of Cloud Condensation Nuclei (CCN) and Ice Forming Nuclei (IFN) (  and ). Space charge is the difference, per unit volume, between total positive and total negative electrical charge that is on droplets, aerosol particles (including the CCN and IFN) and air ions. The cumulative effects of the scavenging in stratiform clouds and aerosol layers in an air mass over the lifetime of the aerosol particles of 1–10 days affects the concentration and size distribution of the CCN, so that in subsequent episodes of cloud formation (including deep convective clouds) there can be effects on droplet size distribution, coagulation, precipitation processes, and even storm dynamics.     

The iron Kα-line diagnostics of the rotating black hole metric in Seyfert galaxies

Observations of Seyfert galaxies in X-ray region reveal the wide emissive lines in their spectra, which can arise in inner parts of accretion disks, where the effects of General Relativity (GR) must be taken into account. A spectrum of a solitary emission line of a hot spot in Kerr accretion disk is simulated depending on the radial coordinate r and the angular momentum a = J/M of a black hole, under the assumption of equatorial circular motion of a hot spot. It is shown that the characteristic two-peak line profile with the sharp edges arises at a large distance (about r ≈ (3–10) r_g). The inner regions emit the line, which is observed with one maximum and extremely wide red wing. We present results of simulations for the isothermal and Shakura–Sunayev disks.     

Transport across complex magnetic environment

It is well known that the rich, highly complicated field configurations of the Heliospheric Magnetic Field has an essential effect on the cross-field propagation of charged particles of relatively low rigidity. In this process the field line separation is the most significant factor of the magnetic structure. The aim of this paper is to develop a framework, in which the particle transport can be investigated for various types of field-line separation. It turns out that the presence of any diffusive cross-field motion of the particles, even if it is arbitrarily small, can wash out most of the details belonging to a specific field line separation model. In that case the overall cross-field motion is forced to become diffusive, irrespectively of the type of the field line spreading, whether this spreading is diffusive, or even non-diffusive in nature. For the various separation models the expressions derived are very similar.     

Free librations of the two-layer Moon and the possibilities of their detection

Based on the forthcoming second stage of the Japanese Lunar mission ILOM (2013), when an optical telescope will be set on the surface near one of the Lunar poles, the possibility to detect free Lunar modes (Chandler-like wobble and free-core nutation) is considered. The difference between the Lunar Eulerian and Chandler-like wobble is explained. The terms “arbitrary libration” and “free libration” are discussed. The geometrical and physical interpretations of the free polar motion over the Lunar surface are considered from the viewpoints of Lunar surface-based observations and the Lunar Navigation Almanac. The dependencies of the free libration period on the core’s radius, density, and ellipticity are modelled and discussed.     

Resonant motion, ballistic escape, and their applications in astrodynamics

A special set of solutions governing the motion of a particle, subject to the gravitational attractions of the Earth, the Moon, and, eventually, the Sun, is discussed in this paper. These solutions, called resonant orbits, correspond to a special motion where the particle is in resonance with the Moon. For a restricted set of initial conditions the particle performs a resonance transition in the vicinity of the Moon. In this paper, the nature of the resonance transition is investigated under the perspective of the dynamical system theory and the energy approach. In particular, using a new definition of weak stability boundary, we show that the resonance transition mechanism is strictly related to the concept of weak capture. This is shown through a carefully computed set of Poincaré surfaces, at different energy levels, on which both the weak stability boundary and the resonant orbits are represented. It is numerically demonstrated that resonance transitioning orbits pass through the weak stability boundaries. In the second part of the paper the solar perturbation is taken into account, and the motion of the resonant orbits is studied within a four-body dynamics. We show that, for a wide class of initial conditions, the particle escapes from the Earth–Moon system and targets an heliocentric orbit. This is a free ejection called a ballistic escape. Astrodynamical applications are discussed.     

Dust in magnetised plasmas: Basic theory and some applications

In this paper we review the theory of charged test particle motion in magnetic fields. This theory is then extended to charged $\text{dust}$ particles, for which gravity and charge fluctuations play an important role. It is shown that systematic drifts perpendicular to the magnetic field and stochastic transport effects may then have to be considered - none of which occur in the case of atomic particles (with the exception of charge exchange reactions). Some applications of charged dust particle transport theory to planetary rings are then briefly discussed.     

磁扰期间的三维漂移壳分离

基于磁层粒子动力学理论,首先对比了计算漂移壳分离的引导中心法和磁力线追踪法,计算表明两种方法的计算结果一致.然后分别采用T89c和T96磁层磁场模式,用磁力线追踪法数值计算了不同初始位置(≤9Re)、不同初始投掷角、不同Kp指数和不同太阳风压力下,带电粒子的漂移壳分离.计算结果揭示了漂移壳分离随初始位置、投掷角、Kp指数和太阳风压力的变化.其具体特征如下. (1)随着径向距离的增大,漂移壳分离效应愈加显著,由正午出发的粒子将被稳定捕获,而午夜出发的径向距离≥7Re的部分大投掷角粒子将沿磁层顶逃逸. (2)正午出发的粒子,漂移到午夜时其漂移壳随投掷角减小向外排列;午夜出发的粒子,漂移到正午时其漂移壳随投掷角增大排列; 90°投掷角粒子在磁赤道面的漂移壳沿着磁场等值线排列. (3)漂移壳分离随Kp指数和太阳风压力增大变得显著,且随这两种扰动参数的变化特征和趋势是基本相似的.      

Multi-case Study of Current Sheet Flapping Motions Induced by Non-adiabatic Ions

In this paper, the y-component of magnetic field line curvature in the plasma sheet was analyzed, and two kinds of shear structures of the flapping current sheet were found, i.e. symmetric and antisymmetric. The alternating bending orientations of the guiding field are exactly corresponding to alternating north-south asymmetries of the bouncing ion population in the sheet center. Those alternating asymmetric plasma sources consequently induce the current sheet flapping motion as a driver. In addition, a substantial particle population with downward motion was observed in the center of a bifurcated current sheet. This population is identified as the quasi-adiabatic particles, and provides a net current opposite to the conventional cross-tail current.      

Simulation study of motion of charged particles trapped in Earth’s magnetosphere

This article aims to understand the motion of the charged particles trapped in the Earth’s inner magnetosphere. The emphasis is on identifying the numerical scheme, which is appropriate to characterize the trajectories of the charged particles of different energies that enter the Earth’s magnetosphere and get trap along the magnetic field lines. These particles perform three different periodic motions, namely: gyration, bounce, and azimuthal drift. However, often, the gyration of the particle is ignored, and only the guiding center of the particle is traced to reduce the computational time. It is because the simulation of all three motions (gyro, bounce, and drift) together needed a robust numerical scheme, which has less numerical dissipation. We have developed a three-dimensional test particle simulation model in which the relativistic equation of motion is solved numerically using the fourth and sixth-order Runge-Kutta methods. The stability of the simulation model is verified by checking the conservation of total kinetic energy and adiabatic invariants linked with each type of motion. We found that the sixth-order Runge-Kutta method is suitable to trace the complete trajectories of both proton and electron of a wide energy range, 5 keV to 250 MeV for L = 2 – 6. We have estimated the bounce and drift periods from the simulations, and they are found to be in good agreement with the theory. The study implies that a simulation model with sixth-order Runge-Kutta method can be applied to the time-vary, non-analytical form of magnetic configuration in future studies to understand the dynamics of charged particles trapped in Earth’s magnetosphere.     

Energization of charged particle in a time-dependent chaotic magnetic field with an implication of the production of seed particles in solar energetic particle events

We investigate the acceleration of charged particles in a time-dependent chaotic magnetic field in this work. In earlier works, it has been demonstrated that in an asymmetric wire-loop current systems (WLCSs), the magnetic field is of chaotic in nature. Furthermore, observations also showed that there exist time-varying current loops and current filaments in solar corona. It is therefore natural to conceive that the magnetic field on the solar surface is chaotic and time-dependent. Here, we develop a numerical model to study the acceleration process of charged particles in a time-varying chaotic magnetic field that is generated by an ensemble of 8 WLCSs. We found that the motion of energetic particles in the system is of diffusive in nature and a power law spectrum can quickly develop. The mechanism examined here may serve as an efficient pre-acceleration mechanism that generates the so-called seed particles for diffusive shock acceleration at a coronal mass ejection (CME) driven shock in large solar energetic particle (SEP) events.     

The dynamical environments analysis of surface particles for different shaped asteroids

Asteroids are coming to be a popular topic in the areas of astrophysical studies and deep space exploration recently. However, surface dynamics of asteroids is still a difficult problem. This paper aims at the motion analysis of surface particles for different asteroids. The dynamical analysis method of particles’ movement is given for three parts: global motion trend, local motion trend and static analysis. A dimensionless parameter ζ is defined to distinguish the predominant term to determine the distribution of effective potential. Three kinds of common asteroids: spheroidal asteroid, spindle-shaped asteroid and dumbbell-shaped asteroid are all discussed for those three parts with different parameter ζ. The motion trend of particles on the surface of each kind of asteroid is given. The static analysis of surface particles for different asteroids is also illustrated. Based on them, some common rules for different shaped asteroids are revealed. This paper could not only provide a reference for asteroid exploration missions but also be meaningful for the research of morphologic evolution of asteroids.     

背景压强对电推进羽流场影响的数值模拟

真空舱内背景压强是电推力器地面试验过程中影响工作性能评估和羽流场参数诊断的重要参数。针对LIPS-200型离子推力器羽流场参数的数值仿真中采用的背景压强建立方法进行了仿真分析。仿真中采用混合粒子网格(PIC)方法和直接模拟蒙特卡罗(DSMC)方法处理羽流场中等离子体运动和粒子间碰撞,分别采用虚拟粒子和计算粒子建立压强的方式,对电推进羽流场进行了数值模拟,并与绝对真空环境进行对比分析。结果表明:背景压强的存在导致中性粒子和电荷交换离子数密度较绝对真空环境高1个量级以上。虚拟粒子可大幅提高计算效率,获得的流场中电荷交换离子分布与计算粒子结果相近,但中性粒子分布相差较大,虚拟粒子无法表征壁面及真空泵的影响。      

基于光场成像的燃烧颗粒粒径与速度测量方法

为实现固体推进剂中金属燃料动态燃烧的精细化表征,提出一种基于光场成像的固体推进剂中燃烧颗粒粒径与速度三维动态同步测量方法。构建固体推进剂燃烧光场成像测量系统,获得金属颗粒动态燃烧过程的光场图像,通过重聚焦算法得到燃烧颗粒场的重聚焦图像;进一步开展光场相机标定实验,获取深度与最佳重聚焦参数的关系曲线;利用图像分割算法对重聚焦图像中的颗粒进行识别和定位,获取燃烧颗粒的粒径大小,并结合三维粒子跟踪技术对颗粒的轨迹和速度进行重构;开展推进剂药条燃烧实验验证研究。结果表明：光场成像技术能够获得不同深度的颗粒信息,跟踪颗粒的动态燃烧过程,并实现对燃烧金属颗粒粒径与三维速度的同步测量。     

A fine mist of very small comet dust particles

A comet nucleus considered as an aggregate of interstellar dust would produce a mist of very finely divided (radius ～ 0.01 μm) particles of carbon and metal oxides accompanying the larger dust grains. These small particles which are very abundant in the interstellar dust size spectrum would provide substantial physical effects because of their large surface area. They may show up strongly in particle detectors on the Halley probes. A strong basis for serious consideration of these particles comes from the other evidence that interstellar dust grains are the building blocks of comets; e.g. (1) the explanation of the “missing” carbon in comets; (2) The S₂ molecule detection which suggests that the comet solid ice materials have been previously subjected to ultraviolet radiation (as are interstellar grains) before aggregation into the comet; (3) the predicted dust to gas ratio.     

带电粒子在近地球区运动的计算方法

基于磁层粒子动力学理论,根据单粒子方法和偶极磁场模型,在近地球区(L＜10),详细讨论了带电粒子的运动特征以及适合不同能量的带电粒子的计算方法,定量分析了磁层中的电场和磁场对各种能量的带电粒子漂移运动的影响.结果表明,对于能量低于105 eV的电子和102 eV的质子,宜采用引导中心近似方法;对于能量介于105～108 eV之间的电子和102～108 eV之间的质子,可在部分区域内采用引导中心近似方法,若运动区在10 Re内,则只能采用变步长的轨道法;而对于引导中心方法,粒子能量高于105 eV时可以忽略电场漂移的影响,粒子能量低于103 eV时不必考虑磁场漂移的影响,从而简化了引导中心方程组,提高了数值计算效率.