高纬极区离子速度分布函数多项式解及非相干散射谱计算

为更准确描述高纬极区电离层离子分布函数, 分别采用弛豫碰撞模型和麦克斯韦分子碰撞模型描述玻耳兹曼方程的碰撞项, 通过求解两种模型下的输运方程, 分别得到两种模型基于麦克斯韦分布下离子分布函数的13矩近似和基于双麦克斯韦分布下离子分布函数的16矩近似. 进一步根据Sheffield理论, 利用两种模型下离子分布函数的13矩和16矩近似, 计算了非相干散射谱, 并对结果进行对比分析. 结果表明, 相对于弛豫碰撞模型, 麦克斯韦分子碰撞模型能更好地描述电离层E层中离子与中性成分的相互作用. 相对于离子分布函数的13矩近似, 16矩近似更适合描述由于电场增加导致的离子温度各向异性特征.      

高纬极区电离层中离子声波谱线增强的解释

大量的实验研究表明, 在顶部极光区电离层, 利用EISCAT非相干散射雷达和Millstone Hill雷达可以观测到不对称、增强的离子声波谱线. 考虑到低能H+离子束沉降到背景电离层, 以及电离层顶部O+离子的外流事件, 采用含有场向热流项的双麦克斯韦分布描述H+离子束的分布函数, 进而基于离子－离子双流不稳定性理论来解释增强的离子声波谱线. 场向热流的引入可以减小离子声波谱线的不对称现象, 这样得到的结果更符合实际.      

电离层Es离子分布数值建模

从离子连续性方程和动量方程出发,比较全面地考虑了太阳光电离、星光电离、地冕电离、星际背景电离和流星离子流等电离源,综合分析Es层的主要动力学过程和光化学过程,以风剪切理论为基础,研究中性成分碰撞、电场作用、金属离子作用和分子离子作用等机制对离子分布的影响,建立了一维时变电离层Es物理模型.对离子产生率、垂直方向的离子速度在高度和时间上的变化进行仿真和计算,得到了在电场及风剪切作用下的离子密度剖面24h内的变化规律.根据建立的模型,以昆明相干散射雷达观测的径向风结果作为模型输入参考,仿真并得到了电离层Es电子密度的时空分布,据此反演出电离层Es临界频率f₀E_s.与昆明站电离层测高仪同时间观测的情况进行对比,结果比较一致,验证了建立的Es层物理模型正确性.      

基于各向异性分布的非相干散射雷达谱仿真

基于离子在加热条件下呈现各向异性分布的假设,建立离子声波激发条件下非平衡态的三维离子速度分布函数模型,给出任意视线方向上非相干散射谱的计算方法,分析幅度修正因子、离子声速粒子组分、雷达波束方向与磁场方向夹角及温度各向异性对功率谱的影响.研究结果可以解释离子线谱超过一个量级的幅度增强现象和上行离子线强于下行离子线的不对称现象,对加热强扰动电离层条件下的等离子体参数反演具有重要意义.      

热流效应对场向电流确定的影响

给出了任意速度分布函数条件下复频率的介电函数和色散关系的求解方法. 讨论了16矩近似条件下, 场向热流对电子速度分布函数, 介电函数实部与虚部, 以及离子与电子谐振频率和阻尼率的影响, 并与平衡态时麦克斯韦分布等离子体的计算结果进行了比较. 忽略场向热流效应, 正向电子漂移速度条件下, 上行等离子线探测, 会过高估计场向电流; 负向电子漂移速度条件下, 上行等离子线探测会过低估计场向电流. 对上下行等离子线同时探测情况, 忽略热流效应同样会过高估计场向电流.      

近壁区Reynolds应力输运方程的理论模型

在湍流边界层的近壁区,采用对称共振三波的理论模型描述相干结构,根据理论模型对Reynolds应力输运方程逐项进行计算和分析,结果与直接数值模拟符合很好.这不仅在理论上有益于对湍流物理机制的了解,而且展现了一种可能性,即根据相干结构的理论知识来改造湍流模型,使之具有更清晰的物理内涵以提高近壁区的预估精度.     

银河宇宙线的起源和初级质子谱

本文采用无界银河系的爆发点源各向同性弥散的模型研究了初级质子的输运方程.用加速器实验得到的P-P碰撞截面的最新资料和一个合理的扩散系数, 得到了该方程的解, 进而讨论了该解的一些特点以及诸多类似点源对地球附近质子谱的可能贡献.该模型所预期的质子谱能较好地说明实验结果.      

匹配滤波方法在非相干散射雷达测量空间碎片中的应用

利用2010年3月25日欧洲非相干散射雷达的常规电离层探测数据, 分析了实测数据中接收功率突变的现象. 高度发生在500 km以上的功率突变一般不是由于地球物理影响产生的, 而要归因于卫星或空间碎片. 采用匹配滤波方法实现了对空间碎片的探测和提取, 利用相干积累的方法提高了信噪比, 共获得363个碎片, 得到了碎片随高度和时间的分布. 结果表明, 碎片高度分布主要集中在800 km左右, 与欧洲非相干散射组织现有的碎片数据具有一致性.      

时变海面复反射系数建模

基于二维时变海面模型和粗糙面电磁散射的低阶小斜率近似方法,研究了海面复反射系数的时变和统计特性.在电磁散射幅度基础上推导了粗糙面反射系数及其相干和非相干分量,相干分量理论计算公式与经验模型一致.仿真结果表明相干分量幅度与经验模型基本吻合,但相位中存在明显的起伏特性,这种起伏在经验公式中没有体现出来.引入Middleton相位差统计模型为反射系数时变特性建模,仿真数据统计与理论分布吻合较好,且概率密度函数的等高图可以很好地表征时变复反射系数间的相关特性.这种吻合性表明复反射系数在短时间尺度上满足高斯分布,但其幅度和频率被海浪长波调制.      

TVD格式与高超音速气动加热数值模拟

为了在化学非平衡流动中获得准确的流场解以及表面热流分布,将总变差减小TVD(Total Variation Diminishing)格式中的熵修正函数,由各向同性分布改为各向异性分布,同时让熵修正函数中的参数与流场中的压力梯度分布相关.将改进后的熵修正函数运用到高超音速化学非平衡绕流流场的数值模拟中,获得了较使用原有熵修正函数更为准确的流场参数和表面热流分布.采用改进的熵修正函数,可以提高壁面附近的粘性分辨率,降低热流计算结果对壁面附近法向网格尺度的敏感性.      

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

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

离子推力器羽流热效应仿真分析

离子推力器工作时向外喷出的羽流与航天器表面碰撞,会引起敏感材料热变形等热效应,严重时会导致航天任务失败。针对兰州空间技术物理研究所研制的LIPS-200型离子推力器羽流热效应进行了仿真分析。仿真中,使用粒子网格（PIC）方法处理等离子体运动,使用直接模拟蒙特卡罗（DSMC）方法处理粒子间碰撞,使用Maxwell模型处理粒子与壁面的能量交换,对电推进羽流热效应测量中的部分测点进行了数值模拟。结果表明,仿真结果与实验数据符合较好,离子推力器出口轴线上滞止热流仿真值与实验测量值误差小于17.0%。此外,热流计对流场的影响主要集中在热流计附近0.1 m范围内,对整体流场影响较小。      

太阳风质子的波能串级-回旋共振加热机制的初级模型

空间飞船Helios 1和2的观测表明, 由0.3AU, 至1AU, 质子的磁矩T_p⊥/B随日心距离增加而增加, 而纵向绝热不变量T_p∥B2/N_P2则随日心距离增加而下降。这说明质子在垂直于行星际磁场方向上受到加热, 而在平行磁场的方向上受到冷却。以往没有理论能满意地解释上述现象, 本文在Alfvén脉动串级理论的基础上利用回旋波的准线性理论分析研究了串级能量转化为太阳风质子热能的机制, 解释了上述观测事实。      

Improved interpretation of solar wind ion measurements via high-resolution magnetic field data

The Wind   spacecraft’s Faraday cups (FC) continue to produce high-quality, in situ observations of thermal protons (i.e., ionized hydrogen) and α$\alpha$-particles (i.e., fully ionized helium) in the solar wind. By fitting a Wind/FC ion spectrum with a model velocity distribution function (VDF) for each particle species, values for density, bulk velocity, and temperature can be inferred. Incorporating measurements of the background magnetic field from the Wind Magnetic Field Investigation (MFI) allows perpendicular and parallel temperature components to be separated. Prior implementations of this analysis averaged the higher-cadence Wind/MFI measurements to match that of the Wind/FC ion spectra. However, this article summarizes recent and extensive revisions to the analysis software that, among other things, eliminate such averaging and thereby account for variations in the direction of the magnetic field over the time taken to measure the ions. A statistical comparison reveals that the old version consistently underestimates the temperature anisotropy of ion VDF’s: averaging over fluctuations in the magnetic field essentially blurs the perpendicular and parallel temperature components, which makes the plasma seem artificially more isotropic. The new version not only provides a more accurate dataset of ion parameters (which is well suited to the study of microkinetic phenomena), it also demonstrates a novel technique for jointly processing particle and field data. Such methods are crucial to heliophysics as wave-particle interactions are increasingly seen as playing an important role in the dynamics of the solar wind and similar space plasmas.     

COSMIC PLASMA DYNAMO

A new dynamo model based on the polarization of plasma is presented in this paper. From the Maxwell equations in a moving medium, a magnetization vector can be caused with Rongon current. The steady solar magnetic field is solved from the equations. On the assumption that the meridianal flow is ignored, the distribution of magnetic field is put out. In the model, there is no additional parameter considered. The intensity of magnetic field inside the sun ranges from 1-6 T. The surface magnetic field around the pole is in the order of 1 ×10-3 T, at low latitude the calculated surface magnetic field has the order of 1×10-2 T. The maximum magnetic field is around 30° in latitude.     

Current collection through the transport of electrons across magnetic field lines

Collection of electrons by a long conducting cylinder in a flowing plasma is studied by means of numerical simulations. The plasma flow simulates the relative motion between a spacecraft and plasma. The sheath structures and the levels of electron current collections for the cases with and without an ambient magnetic field ( ) are studied. It is found that for the flow perpendicular to the magnetic field, the current is considerably enhanced depending on the relative drift velocity. In the case of a non-zero magnetic field perpendicular to the cylinder axis, the potential structure is a two-dimensional double layer with dimensions L L_|, where L and L_| are the dimensions perpendicular and parallel to , respectively. L is found to be the current limiting radius given by the Parker-Murphy model. For the flow along , the electron current is found to be smaller than that for the flow perpendicular to . This is explained in terms of the potential structures.     

Electromagnetic ion cyclotron resonance heating in the VASIMR

Plasma physics has found an increasing range of practical industrial applications, including the development of electric spacecraft propulsion systems. One of these systems, the Variable Specific Impulse Magnetoplasma Rocket (VASIMR) engine, both applies several important physical processes occurring in the magnetosphere. These processes include the mechanisms involved in the ion acceleration and heating that occur in the Birkeland currents of an auroral arc system. Auroral current region processes that are simulated in VASIMR include lower hybrid heating, parallel electric field acceleration and ion cyclotron acceleration. This paper will focus on using a physics demonstration model VASIMR to study ion cyclotron resonance heating (ICRH). The major purpose is to provide a VASIMR status report to the COSPAR community. The VASIMR uses a helicon antenna with up to 20 kW of power to generate plasma. This plasma is energized by an RF booster stage that uses left hand polarized slow mode waves launched from the high field side of the ion cyclotron resonance. The present setup for the booster uses 2–4 MHz waves with up to 20 kW of power. This process is similar to the ion cyclotron heating in tokamaks, but in the VASIMR the ions only pass through the resonance region once. The rapid absorption of ion cyclotron waves has been predicted in recent theoretical studies. These theoretical predictions have been supported with several independent measurements in this paper. The ICRH produced a substantial increase in ion velocity. Pitch angle distribution studies show that this increase takes place in the resonance region where the ion cyclotron frequency is equal to the frequency on the injected RF waves. Downstream of the resonance region the perpendicular velocity boost should be converted to axial flow velocity through the conservation of the first adiabatic invariant as the magnetic field decreases in the exhaust region of the VASIMR. In deuterium plasma, 80% efficient absorption of 20 kW of ICRH input power has been achieved. No evidence for power limiting instabilities in the exhaust beam has been observed.