磁控进气道二维性能计算

超燃冲压发动机进气道通常都是针对某一特定的马赫数进行设计的,当飞行马赫数不在设计马赫数时发动机性能急剧降低,磁控进气道的应用很有可能解决这一技术难题.对磁控进气道内的流动情况进行了二维数值模拟.飞行马赫数大于设计马赫数时采用磁控进气道可以调整激波位置使激波回到设计点,并减小燃烧室入口处的马赫数.分析了霍尔效应对磁控进气道性能的影响,结果表明分段电极可以有效减小霍尔效应的不利影响.      

基于磁流体控制的高超声速进气道黏性效应

建立引入电磁源项的二维低磁雷诺数磁流体动力学(MHD)方程组,对高超声速二维前体/进气道黏性流场进行了数值模拟.在给出了进气道高于设计马赫数的非设计工况下黏性流场的基本特征基础上,进一步分析了施加MHD控制对进气道黏性效应的影响.结果表明:施加MHD控制可以有效抑制非设计工况下内进气道表面的附面层分离,改善上壁面的热状况,平衡上、下壁面之间的热负担;黏性作用下,进气道流场及性能参数随磁感强度的变化规律与无黏模型计算结果存在较大差别,对磁流体控制的高超声速进气道研究不可忽略黏性的影响.      

非结构网格高超声速化学非平衡MHD流场数值模拟

 在非结构网格上对考虑化学非平衡效应的二维高超声速磁流体绕钝头体流动进行了数值模拟。控制方程由二维理想磁流体动力学（MHD）方程和组元连续方程两部分组成,化学动力学模型为5组元17反应模型。MHD方程空间离散采用AUSM格式,时间推进采用显式5步龙格-库塔格式,并通过弱耦合的方式与化学反应控制方程结合在一起。计算模型为二维钝头体,外加磁场为偶极子场,磁场源位于钝头体内部。在高超声速来流条件下,对有、无磁场干扰,是否考虑化学反应下的4种工况进行了数值计算,得到了满意的结果,并与有限的参考文献进行了对比。结果表明本文发展的方法能准确地模拟考虑化学非平衡效应的高超声速MHD流场。     

基于电子束电离的磁流体力学进气道流动控制数值模拟

采用基于电子束电离的磁流体力学(MHD)控制系统,对高超声速流场附面层,以及非设计状态下的高超声速进气道流场的磁流体控制进行了深入研究.控制方程为低磁雷诺数Navier-Stokes方程,采用等离子体动力学模型与电子束模型模拟空气电离过程.研究结果表明:①电子束电离能有效提高流场的电导率,增强磁场对流场的控制效率;②基于电子束诱导电离的MHD控制系统能有效地控制高超声速流场的附面层,但其控制效率跟电子束能量大小相关;③基于电子束诱导电离的MHD控制系统能有效地改变非设计状态下高超声速飞行器的斜激波结构,使进气道重新满足Shock-on-lip(SOL)条件,但进气道的总压恢复系数以及流量将会降低.      

三维磁流体动力学管道流动加减速控制数值研究

阐述了磁流体动力学(MHD)控制流场作用机理.在小磁雷诺数条件下,运用数值模拟方法,对不同外加电磁场条件下三维MHD管道流动的流场情况进行研究,得出不同磁场、电场等MHD参数作用下MHD加速器的性能.对MHD加速器的应用前景进行了展望.计算结果显示,电场取6000V/m,磁场为0.92T 情况下,流场的加速性能可达13.46%,并且可以通过减小磁场或增大电场进一步提高加速性能;在仅添加0.92T磁场条件下,速度减速可达16.35%.     

Numerical analysis of deflection control of a gas plasma jet based on magnetohydrodynamic staggered electrode configuration

To control the deflection of the gas plasma jet, a new analytical method is proposed based on the Magnetohydrodynamic (MHD) technique. Based on the typical MHD power generation model, the applied voltage is applied to the staggered electrodes, that is, a pair of electrodes on the same side wall are connected to generate an axial current in the channel. Under the action of the magnetic field perpendicular to the direction of the flow, the plasma is subjected to electromagnetic forces perpendicular to these two directions, and the jet is deflected. The computational model including the Navier-Stokes equations coupled with electromagnetic source terms, the electric potential equation and Ohm’s law is solved. The deflection of the gas jet under the action of an electromagnetic field is observed, and the maximum deflection angle is about 14.8°. The influences of the electric field, magnetic field, and conductivity on the jet deflection are studied. Results show that although the influences of these three factors on the deflection are similar, and the effect of increasing the electric field strength is slightly greater, priority should be given to increasing the magnetic field strength from the perspective of reducing energy consumption. The Stuart number is introduced to assess the ability of electromagnetic force to control jet deflection. When the electromagnetic parameters are constant, this solution provides better control of low-density and low-speed fluid flows. The calculation results show that using the staggered electrode method configuration is feasible in terms of controlling the deflection of a plasma jet deflection.     

二维超燃冲压发动机磁控进气道的数值模拟

针对飞行马赫数大于设计马赫数的情况,采用二维磁流体动力学方程对磁控进气道进行了数值模拟研究。数值模拟结果表明磁流体装置的电磁作用可以使非设计马赫数下进气道激波满足SOL(shock on lip)条件,并使出口处的流动变得均匀。分析了磁流体作用位置、宽度和深度等关键控制参数对该类进气道性能的影响,计算结果表明,磁流体作用区域越靠近飞行器前缘,而且在纵向越深,则进气道出口处的流动越均匀,但流率会有所下降;若磁流体作用区域较宽,则需较小的磁场就能使非设计马赫数下进气道的激波结构满足SOL条件。     

基于正交实验设计的螺旋通道变流器水动力性能数值分析

变流器是螺旋通道磁流体推进器推进通道中重要的水力部件之一,变流器的设计影响着推进器的水力性能和推进性能。提出了基于四阶贝塞尔曲线的螺旋通道变流器三维建模方法,明确了变流器的结构参数;采用正交实验设计方法,选用五因素四水平正交设计表基于CFD方法分析了变流器结构参数对螺旋通道水动力学性能的影响及其变化规律。分析结果表明螺旋通道导流器会大幅减小水力损失,整流器会增加水力损失,但能提高整流效果。分析结果可以为螺旋通道磁流体推进器的水力优化及推进性能提升提供有意义的指导作用。     

Investigation of MHD power generation with supersonic non-equilibrium RF discharge

Magnetohydrodynamic (MHD) power generation with supersonic non-equilibrium plasma is demonstrated. Capacitively coupled radio frequency (RF) discharge (6 MHz, maximum continual power output of 200 W) was adopted to ionize the Mach number 3.5 (650 m/s), 0.023 kg/m3 airflow. In a MHD channel of 16 mm × 10 mm × 20 mm, MHD open voltage of 10 V is realized in the magnetic field of 1.25 T, and power of 0.12 mW is extracted steadily and con-tinuously in the magnetic field of 1 T. The reasons for limited power generation are proposed as:low conductivity of RF discharge; large touch resistance between MHD electrode and plasma;strong current eddies due to flow boundary layer. In addition, the cathode voltage fall is too low to have obvious effects on MHD power generation.     

MHD控制激波诱导边界层分离的机理

激波边界层相互作用是高超声速飞行器面临的重大问题,激波入射到平板引起的边界层分离是其中最具代表性的一种。用加权的3阶ENN格式计算了小磁雷诺数近似的MHD方程,研究了MHD控制层流边界层分离的机理。数值结果显示,通过局部电离空气并施加洛伦兹力,能使分离点向下游移动,分离区尺寸减小,从而抑制和缓解由于激波-边界层相互作用而引起的分离。     

Thermal damage analysis of aircraft composite laminate suffered from lightning swept stroke and arc propagation

Lightning strike is a complicated process involving multi-field coupling. In order to investigate the thermal damage mechanism of carbon fiber reinforced composites subject to lightning swept stroke, a complete numerical method is presented. Numerical model of lightning discharge is established based on Magneto Hydro Dynamics (MHD) and calculated by FLUENT secondary development technology. Considering aerodynamic flow effect, channel formation and evolution process during lightning discharge is analyzed for lightning current waveform A. Thermal-electric Coupling model is presented according to Radial Basis Function (RBF) interpolation theory, which is implemented by compiling program to make lightning current and heat energy inject into composite laminate. Consequently, damage mechanism of composite laminate under lightning swept stroke is studied based on the coupled numerical model and element deletion method. Ablation damage morphology of composite laminate is analyzed to understand plasma expansion and reattachment in arc root. The results show that aerodynamic flow makes the lightning channel move fast and composite laminate is deteriorated due to thermal damage.     

磁零点及通道宽度对霍尔推力器性能影响的数值研究

磁场位形和通道尺度会改变霍尔推力器等离子体放电过程,影响推力器的宏观放电特性。为分析磁场和通道宽度对推力器放电性能的影响规律,本文针对霍尔推力器轴对称通道结构和放电物理过程建立2D3V物理模型,采用粒子模拟方法研究了霍尔推力器磁零点磁场位形不同通道宽度的电势、粒子数密度、电子温度、电离速率、比冲及推功比的变化规律,结果表明：在具有磁零点磁场位形下,随着通道宽度增加,通道出口处电势降增加,加速区缩短,离子径向速度减少,壁面腐蚀降低;当磁零点位置在内壁面,推力器通道宽度由14 mm增加到16 mm时,推力器比冲和推功比增大,推力器放电效率提高;当磁零点位置在通道中轴线或外壁面,且通道宽度大于14 mm时,推力器比冲增大,推功比减小,推力器效率下降。     

A Class of TVD Type Combined Numerical Scheme for MHD Equations With a Survey About Numerical Methods in Solar Wind Simulations

It has been believed that three-dimensional, numerical, magnetohydrodynamic (MHD) modelling must play a crucial role in a seamless forecasting system. This system refers to space weather originating on the sun; propagation of disturbances through the solar wind and interplanetary magnetic field (IMF), and thence, transmission into the magnetosphere, ionosphere, and thermosphere. This role comes as no surprise to numerical modelers that participate in the numerical modelling of atmospheric environments as well as the meteorological conditions at Earth. Space scientists have paid great attention to operational numerical space weather prediction models. To this purpose practical progress has been made in the past years. Here first is reviewed the progress of the numerical methods in solar wind modelling. Then, based on our discussion, a new numerical scheme of total variation diminishing (TVD) type for magnetohydrodynamic equations in spherical coordinates is proposed by taking into account convergence, stability and resolution. This new MHD model is established by solving the fluid equations of MHD system with a modified Lax-Friedrichs scheme and the magnetic induction equations with MacCormack II scheme for the purpose of developing a combined scheme of quick convergence as well as of TVD property. To verify the validation of the scheme, the propagation of one-dimensional MHD fast and slow shock problem is discussed with the numerical results conforming to the existing results obtained by the piece-wise parabolic method (PPM). Finally, some conclusions are made. This revised version was published online in August 2006 with corrections to the Cover Date.     

超声速氩气流磁流体发电初步实验研究(英文)

利用激波风洞,采用氦气驱动氩气,在平衡接触面运行方式下得到高温气体,通过在低压段注入电离种子K2CO3粉末,实现高温条件下导电流体的产生,开展了超声速氩气流磁流体功率提取初步实验研究。在喷管入口总压0.32MPa、总温6504K,磁场强度约0.5T、喷管出口气流速度1959m/s的条件下,对分段磁流体功率提取通道电极的感应电压和短路电流进行了测量,实验测量结果与理论计算相吻合,并由电压电流计算得出了平均电导率约20S/m左右,在负载系数为0.5的情况下,磁流体功率提取通道最大的功率密度可达4.7971MW/m3,最大焓提取率为0.34%。最后分析并给出了气体状态参数T1,M1,T2,M2的测试原理与方法。     

NND格式在多维理想磁流体方程组中的应用

采用修正的四步Runge-Kutta方法求解三维一般曲线坐标系下的理想磁流体方程组,为克服数值振荡,加特征型NND格式进行后处理.特征型NND格式推广到求解三维磁流体(MHD)问题需要知道雅可比通量的左右特征矩阵,在具体计算时需要克服矩阵的奇性.本文用三维程序采用推广的特征NND格式计算了一维MHD激波管和二维(MHD)喷管流动,计算结果表明,特征NND格式保持了TVD格式高精度的优点,又具有计算简单的特点,在包括强弱间断等复杂波系的定常和非定常MHD流场数值模拟中是成功的.     

Low-Frequency Waves and Instabilities in Collisionless Plasmas

We present a theoretical overview of low-frequency waves and instabilities in collisionless, multi-component plasmas with gyrotropic ( ) thermal pressure. We show that the complete dispersion relation can be obtained in the framework of a mixed magnetohydrodynamic (MHD)-kinetic formalism, which uses the MHD mass, momentum, and induction equations, together with the kinetically corrected version of the double-adiabatic equations of state. The complete dispersion relation contains not only the three standard modes (fast, slow, and Alfvén) from double-adiabatic MHD, but also the mirror mode from kinetic theory. We examine the stability properties of these four modes, firstly in the case of a uniform medium, and secondly in the case of a stratified and rotating medium. We also discuss the connections with the quasi-interchange modes (interchange and translation) often referred to in the context of magnetospheric physics.     

The future of space plasma simulation

Space plasma simulation is a subject which is in its infancy, but which is already having an important impact on space science. Its growth is being spurred onward by the continuing increase in capacity (speed and memory) of computers and by advances in the sophistication of numerical models. These advances are making it possible to simulate more realistic situations using more complex models. Already significant three-dimensional MHD calculations of the magnetosphere and its interaction with the solar wind have been carried out. In addition multi-dimensional particle simulations are illuminating many of the microscopic physics processes which go on (instabilities, saturation levels and wave nonlinearity, shock structure, etc.). Notwithstanding these advances, the surface has only been scratched; many challenges and opportunities are provided by simulation both for the space physicist and the model builders (also for computer designers). In MHD models more physics need to be included (Hall effect, gyroviscosity, accurate models of boundaries, how do we put microscopic physics effects into macroscopic codes, etc.). For model builders correct treatments of systems containing a large range of important space and time scales, magnetic field strengths, Alfven wave velocity, etc. present real challenges. What are the best ways to diagnose these complex models and obtain meaningful information? What quantities should be looked at? How should they be displayed? A discussion of the promises, the prospects, and the challenges of the above topics will be given with examples taken from recent work.     

基于HLLC格式的理想MHD方程数值模拟

针对理想MHD(magnetohydrodynamics)方程数值求解困难的问题,基于原始的HLLC (Harten-Lax-Van Leer Contact wave)近似黎曼解方法,发展出一种新的适用于理想MHD问题的通量计算格式。控制方程采用有限体积法离散,时间推进采用隐式的LU-SGS（lower-upper symmetric Gauss-Seidel）格式,并且采用双曲型散度清除技术来抑制磁场散度的累积。通过一维激波管问题的数值模拟表明,HLLC-MHD格式能准确的分辨并捕捉复杂的磁流体力学波系结构,保证数值计算的精度,最大的数值计算误差不超过10%;通过二维的Rotor问题的数值模拟研究表明,HLLC-MHD格式能够应用于多维理想MHD问题的数值模拟,并且能够准确捕捉多维情况下磁流场中的阿尔文波;通过对比有、无散度清除的Rotor问题计算结果表明,双曲型散度清除技术可以将磁场散度峰值从50降低到2,有效抑制高磁场散度区域的散度累积,同时也会将误差传播到低散度区域,并且引起边界处散度的累积,影响计算的稳定性。