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

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

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

为了研究磁流体动力学(Magnetohydrodynamics:MHD)加速边界层对激波-湍流边界层相互作用的影响,用高阶有限差分法求解了小磁雷诺数近似的MHD湍流方程。其中,无粘通量采用WENN格式离散、粘性通量采用Roe平均中心差分离散,时间采用半隐式推进,并采取追赶法求解。计算给出了湍流、电场、磁场和电导率等参数对边界层分离的影响,数值结果显示:在同样的逆压梯度下,湍流边界层分离能更快地趋于稳态流场,且分离区比层流小;通过施加洛仑兹力加速,边界层速度型面变得更加饱满、位移厚度减小、分离点和再附点向激波与固壁的交点靠近,分离区尺寸减小甚至最终被消除。     

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

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

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

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

三维等离子体MHD气动热环境数值模拟

电磁流动控制技术是一个多学科交叉融合的重要研究方向,在高超声速飞行器气动特性优化、气动热环境减缓、边界层转捩和等离子体分布等流动控制方面显示出广阔的应用前景。考虑高超声速飞行器绕流流场中发生的离解、复合、电离和置换等化学反应,气体分子振动能激发以及化学非平衡效应,耦合电磁场作用并基于低磁雷诺数假设,通过数值模拟求解三维非平衡Navier-Stokes流场控制方程和Maxwell电磁场控制方程,建立磁场与三维等离子体流场耦合数值模拟方法及程序,采用典型算例进行考核。在此基础上,开展不同条件下磁场对再入三维等离子体流场以及气动热环境影响分析。研究表明:建立的高超声速飞行器的等离子体流场与磁场耦合计算方法及程序,其数值模拟结果与文献符合,外加磁场使飞行器头部弓形激波外推,磁场强度越强,激波面外推距离越大;不同磁场强度环境下,流场中温度峰值大小略有变化,变化幅度较小;磁场对绝大部分区域的热流有减缓作用,作用的大小与飞行高度、马赫数以及磁场的配置紧密相关;当前的计算条件下,飞行的高度越高,磁场的作用越明显。     

磁流体加速技术对脉冲爆轰发动机增推影响的二维数值研究

为了研究磁流体加速技术对脉冲爆轰发动机增推过程的影响,利用含磁流体运动的气液两相爆轰模型,同时采用CE/SE方法对含磁流体加速装置的脉冲爆轰发动机点火、爆轰及排气过程进行数值研究。模拟带磁流体加速器的脉冲爆轰发动机流场,揭示磁流体加速器的增推过程,分析磁流体加速器在不同磁感应强度及管长条件下对脉冲爆轰发动机推进性能的影响。计算结果表明:磁流体加速装置能够实现对脉冲爆轰发动机的推力补偿作用。磁感应强度及加速段管长增加,脉冲爆轰发动机的平均推力得到显著提高。加速段管长0.4m,磁感应强度2T时,平均推力可提高30.457%,此时推力提升率达到最大。     

一种新的用于MHD方程逆风格式的Jacobian矩阵分裂方法

针对理想MHD方程,提出了一种新的基于MacCormack算法的雅可比矩阵分裂方法,克服了原有方法稳定性差的问题,并成功地应用于理想MHD方程的求解.控制方程在非结构混合网格上进行空间离散,其中对流项采用本文发展的逆风向量分裂格式,并引入了双曲型磁场散度清除技术,时间推进为显式5步龙格-库塔方法.对MHD激波管流动和带均匀磁场干扰的二维高超声速钝头体绕流流场进行了数值模拟,得到了与参考文献相吻合的数值结果,表明本文发展的数值分裂方法可以有效地捕捉MHD流场的流动特征,并且具有比MacCormack方法更高的稳定性和计算精度.     

MHD加速器模式磁控进气道的优化设计

为提高超燃冲压发动机进气道在非设计状态下的性能,对磁控进气道进行了研究。采用二维磁流体动力学(MHD)模型对加速器模式的磁控进气道进行了数值模拟和参数优化。分析了电磁作用使空气流率增加的原因,选取了一组优化的设计参数进行数值模拟,确定了磁流体关键参数与进气道主要性能参数的匹配原则。分析表明磁场的大小和方向以及电磁作用的位置对进气道性能有重要影响;唇口附近及上方的电磁作用对增加空气流率起到了关键的作用,磁流体加速器可以显著增加进气道的空气捕获率和压缩比,但由于不可逆效应总压恢复系数会减小。研究结果表明,当飞行马赫数小于设计马赫数时加速器模式的磁控进气道可以提高进气道的性能。     

磁场作用下的TIG电弧行为分析

在磁流体动力学理论的基础上,针对外加纵向磁场作用下的钨极氩弧焊(TIG)电弧在旋转圆柱坐标系中建立了不同电流密度的数学模型.通过Fluent用户自定义函数进行编程开发,对其流场和电磁场进行耦合求解,得到了外加纵向磁场作用下的TIG电弧热场和流场的分布规律.结果 表明,纵向磁场作用下的TIG电孤温度场呈现空心钟罩形态,并...     

双楔面诱导的斜爆轰波阵面的磁流体控制

基于高超声速飞行器爆轰推进发展需要,数值研究了由双楔面诱导的斜爆轰波类型及其磁流体(MHD)控制的可行性,探讨了后楔倾角变化对双楔斜爆轰流场与MHD控制的影响.结果表明:对于前、后楔倾角分别为15°与20°的双楔面诱导的斜爆轰,Lorentz力可以将不同来流马赫数条件下的稳定斜爆轰波阵面恢复到设计位置,并可以使不稳定斜爆轰波趋于稳定,但是,无法将不稳定斜爆轰波恢复到设计位置.而当后楔倾角稍微增加时,若爆轰流场稳定,同样可对斜爆轰波阵面进行有效控制.但当后楔倾角大于25°时,稳定斜爆轰波流场失稳,此时MHD无法控制流场稳定性.但若增加气体反应速率,虽然斜爆轰波不稳定,MHD控制却可使其趋于稳定,但无法将爆轰波阵面恢复到设计位置.      

三维磁流体强化超燃冲压发动机数值模拟

建立了三维磁流体强化超燃冲压发动机内部黏性流场的求解模型.针对马赫数为6设计了联合应用磁控进气道和磁流体能量旁路的磁流体强化超燃冲压发动机模型.针对该模型进行了数值模拟研究,分析其中的三维流场结构、电参数分布规律以及能量转换特性.结果表明:当飞行马赫数为8时,磁控进气道的应用能够使头部压缩激波回到唇口,使分离区消失,内进气道中的流动恢复到设计状态.磁流体能量旁路可有效降低燃烧室入口处的马赫数,从而改善发动机性能.其中发生器中的流动参数和电参数的分布比较理想,效果显著;而加速器要取得显著的加速效果则需要大量的能量输入.在加速器中,电极附近焦耳耗散严重,导致局部高温以及流动的复杂性,性能不够理想.      

超声速磁流体加速实验及一维模型分析

为了获得负载系数、电导率等参数变化对超声速磁流体加速效果的影响规律,利用激波风洞,采用氩气与碳酸钾作为工质,电容提供电能的方式,在磁感应强度为0.5T的条件下,进行了不同电容充电电压下的超声速磁流体加速实验研究,并对一维定常理想分段法拉第型磁流体加速模型进行了分析.通过实验获得了不同电容充电电压下#10电极间的电压、电流、负载系数、电导率及#20电极开路电压等数据,在300, 400V电容充电电压下,气流速度分别增加11.4%和24.0%,在500V电容充电电压下气流速度减小11.1%.实验及模型分析得出不同的负载系数会使超声速磁流体处于加速或减速的不同状态,而电导率会影响注入总能量的大小,使磁流体流动的速度梯度大小发生改变.      

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...     

磁流体湍流中的结构和能量传输

磁流体湍流中存在多种相干结构,包括电流和涡旋结构。文章主要对近几年有关磁流体湍流相干结构及其相关的能量传输的一些工作进行了介绍。不同于中性流体湍流中的管状结构,磁流体中的相干结构多呈现为片状,强电流片附近也存在强涡片结构。磁流体湍流中的能量级串使能量跨尺度从大尺度传递至小尺度;与中性流体湍流能量传输主要集中在相近尺度上不同,磁流体湍流除了动能传输,还有磁能传输,以及动能与磁能之间的跨尺度互相转换。研究表明,磁流体湍流中的能量传输及耗散具有强间歇性,集中在较小区域范围,且与相干结构存在关联。在存在背景磁场的情况下,磁流体湍流中的相干结构在背景磁场方向上被拉长,而同样方向上的能量传输受到抑制。在高马赫数的情况下,激波的产生会使能量传输明显增强,同时带来很强的间歇性,结构函数标度律远远偏离线性标度关系,并且随着阶数的增加,标度指数会达到饱和值。     

Emerging Parameter Space Map of Magnetic Reconnection in Collisional and Kinetic Regimes

In large-scale systems of interest to solar physics, there is growing evidence that magnetic reconnection involves the formation of extended current sheets which are unstable to plasmoids (secondary magnetic islands). Recent results suggest that plasmoids may play a critical role in the evolution of reconnection, and have raised fundamental questions regarding the applicability of resistive MHD to various regimes. In collisional plasmas, where the thickness of all resistive layers remain larger than the ion gyroradius, simulations results indicate that plasmoids permit reconnection to proceed much faster than the slow Sweet-Parker scaling. However, it appears these rates are still a factor of ～10× slower than observed in kinetic regimes, where the diffusion region current sheet falls below the ion gyroradius and additional physics beyond MHD becomes crucially important. Over a broad range of interesting parameters, the formation of plasmoids may naturally induce a transition into these kinetic regimes. New insights into this scenario have emerged in recent years based on a combination of linear theory, fluid simulations and fully kinetic simulations which retain a Fokker-Planck collision operator to allow a rigorous treatment of Coulomb collisions as the reconnection electric field exceeds the runaway limit. Here, we present some new results from this approach for guide field reconnection. Based upon these results, a parameter space map is constructed that summarizes the present understanding of how reconnection proceeds in various regimes.     

Electrodynamic coupling of the magnetosphere and ionosphere

The auroral zone ionosphere is coupled to the outer magnetosphere by means of field-aligned currents. Parallel electric fields associated with these currents are now widely accepted to be responsible for the acceleration of auroral particles. This paper will review the theoretical concepts and models describing this coupling. The dynamics of auroral zone particles will be described, beginning with the adiabatic motions of particles in the converging geomagnetic field in the presence of parallel potential drops and then considering the modifications to these adiabatic trajectories due to wave-particle interactions. The formation of parallel electric fields can be viewed both from microscopic and macroscopic viewpoints. The presence of a current carrying plasma can give rise to plasma instabilities which in a weakly turbulent situation can affect the particle motions, giving rise to an effective resistivity in the plasma. Recent satellite observations, however, indicate that the parallel electric field is organized into discrete potential jumps, known as double layers. From a macroscopic viewpoint, the response of the particles to a parallel potential drop leads to an approximately linear relationship between the current density and the potential drop.The currents flowing in the auroral circuit must close in the ionosphere. To a first approximation, the ionospheric conductivity can be considered to be constant, and in this case combining the ionospheric Ohm's Law with the linear current-voltage relation for parallel currents leads to an outer scale length, above which electric fields can map down to the ionosphere and below which parallel electric fields become important. The effects of particle precipitation make the picture more complex, leading to enhanced ionization in upward current regions and to the possibility of feedback interactions with the magnetosphere.Determining adiabatic particle orbits in steady-state electric and magnetic fields can be used to determine the self-consistent particle and field distributions on auroral field lines. However, it is difficult to pursue this approach when the fields are varying with time. Magnetohydrodynamic (MHD) models deal with these time-dependent situations by treating the particles as a fluid. This class of model, however, cannot treat kinetic effects in detail. Such effects can in some cases be modeled by effective transport coefficients inserted into the MHD equations. Intrinsically time-dependent processes such as the development of magnetic micropulsations and the response of the magnetosphere to ionospheric fluctuations can be readily treated in this framework.The response of the lower altitude auroral zone depends in part on how the system is driven. Currents are generated in the outer parts of the magnetosphere as a result of the plasma convection. The dynamics of this region is in turn affected by the coupling to the ionosphere. Since dissipation rates are very low in the outer magnetosphere, the convection may become turbulent, implying that nonlinear effects such as spectral transfer of energy to different scales become important. MHD turbulence theory, modified by the ionospheric coupling, can describe the dynamics of the boundary-layer region. Turbulent MHD fluids can give rise to the generation of field-aligned currents through the so-called -effect, which is utilized in the theory of the generation of the Earth's magnetic field. It is suggested that similar processes acting in the boundary-layer plasma may be ultimately responsible for the generation of auroral currents.     

MHD Flow Control of Oblique Shock Waves Around Ramps in Low-temperature Supersonic Flows

This article is devoted to experimental study on the control of the oblique shock wave around the ramp in a low-temperature supersonic flow by means of the magnetohydrodynamic(MHD) flow control technique. The purpose of the experiments is to take advantage of MHD interaction to weaken the oblique shock wave strength by changing the boundary flow characteristics around the ramp. Plasma columns are generated by pulsed direct current(DC) discharge, the magnetic fields are generated by Nd-Fe-B rare-earth permanent magnets and the oblique shock waves in supersonic flow are generated by the ramp. The Lorentz body force effect of MHD interaction on the plasma-induced airflow velocity is verified through particle image velocimetry(PIV) measurements. The experimental results from the supersonic wind tunnel indicate that the MHD flow control can drastically change the flow characteristics of the airflow around the ramp and decrease the ratio of the Pitot pressure after shock wave to that before it by up to 19. 66%, which leads to the decline in oblique shock wave strength. The oblique shock waves in front of the ramp move upstream by the action of the Lorentz body force. The discharge characteristics are analyzed and the MHD interaction time and consumed energy are determined with the help of the pulsed DC discharge images. The interaction parameter corresponding to the boundary layer velocity can reach 1. 3 from the momentum conservation equation. The velocity of the plasma column in the magnetic field is much faster than that in the absence of magnetic field force. The plasma can strike the neutral gas molecules to transfer momentum and accelerate the flow around the ramp.     

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.