流向磁场作用下二维磁流体槽道湍流直接模拟

对低磁雷诺数近似下流向磁场作用的二维磁流体槽道湍流进行了直接数值模拟(DNS,Direct Numerical Simulation),给出了 Re=10 000时不同磁相互作用数下二维磁流体槽道湍流的近壁速度分布、湍流脉动速度均方根、雷诺应力等统计量的变化,并与中性流体二维槽道湍流进行了比较.结果表明:流向磁场作用会导致对数区上移;雷诺应力最大值随磁相互作用数呈线性变化;随着磁相互作用数的增大,下壁面平均涡量的时间演化由拟周期性向周期性转变,且脉动周期逐渐增大并当流动层流化后下壁面平均涡量成为常值.      

一种新的烧蚀热响应算法在质量估算中的应用

针对预测飞行器再入过程中气动热流引发的烧蚀热响应导致热防护罩表层材料质量损耗的问题,研究了热防护罩的几何模型和烧蚀质量估算的方法,通过建立三自由度再入轨迹动力学方程,应用修正的牛顿流体理论计算气动系数,以及Detra-Kemp-Riddell和Tauber-Sutton理论计算驻点热流密度和热辐射,利用一维非线性热传导方程模拟了碳化材料的烧蚀过程,提出了基于Newton-Raphson和TDMA的烧蚀热响应算法估计飞行器热防护罩质量损耗的方法。通过分析,实现了再入全过程热防护材料烧蚀深度连续动态变化的预测,能够有效替代热平衡积分法,估算的烧蚀质量为优化热防护罩的几何模型和再入轨迹提供了参考依据。     

火星再入飞行器风洞试验与真实飞行之间相关性的探讨

由于风洞试验条件限制,难以完全模拟火星再入飞行器真实飞行环境,因此需要建立火星再入飞行器风洞条件与真实飞行之间的关联关系。基于国外文献公开数据,采用数值方法和对比分析方法探讨了类"探路者号"外形的火星再入飞行器的风洞试验与真实飞行之间的外推方法。结果表明,在高焓空气风洞和常规空气风洞试验条件下,可以将模型驻点附近的无量纲压力和压力系数作为相关性参数,将风洞条件与飞行条件相关联起来,但是不能直接利用风洞试验的热流、无量纲热流和Stanton数作为关联参数;在高焓CO₂风洞试验条件下,可以利用模型驻点附近的无量纲压力、压力系数和Stanton数作为外推参数,但是不能直接将风洞试验的热流、无量纲热流作为相关性参数,将风洞条件下的风洞数据通过外推获取飞行条件下飞行器的性能参数。      

球头锥-柱再入飞行器的动力环境预示

考察了超声速和高超声速流场中一类球头锥-柱再入飞行器表面稳定压力和脉动压力的分布特性,利用CFX(Computation Fluid Explorer)软件建立球头锥-柱的外流场域网格,对其进行定常、非定常流场分析,给出一套预测表面稳定压力和脉动压力环境的工程新方法.将该方法与理论和经验估算方法进行了对比分析,验证了该方法的可行性和适用性.同时利用该方法研究了马赫数、来流攻角等因素对再入飞行器表面脉动压力环境的影响.利用ANSYS大型软件建立了再入飞行器导弹头部的动力学分析模型,根据确定的再入飞行器表面脉动压力环境和随机振动分析理论,采用有限元谱分析方法计算得到该飞行器弹头部件的随机振动响应加速度的均方根值和功率谱密度曲线,计算结果满足工程应用的要求.      

高温燃气风洞加热特性仿真分析

提出了利用亚音速高温燃气流进行近空间高超飞行器热环境地面模拟的试验方案,在试验装置试验段,通过高温高速的燃气流引射低速的冷气流,达到仅使飞行器头锥驻点附近区域产生局部高温而其余头锥蒙皮表面低温的目的.对某型高超音速飞行器的头锥利用高温燃气进行加热并利用CFD(Computational Fluid Dynamics)方法对13种模拟工况进行了数值仿真分析.将数值模拟计算结果与飞行器在高超音速飞行状态下对应机体部位气动热的理论计算值进行了对比,证实了亚声速高温燃气热环境模拟方法的可行性,为高温燃气地面模拟设备技术方案论证提供了依据.      

等离子体磁流体发电研究进展

简要分析了等离子体磁流体发电系统的工作原理和发电过程,阐述了目前磁流体发电研究中的重点和关键问题,从数值模拟和实验研究两方面回顾了国内外的研究情况和研究进展。分析认为等离子体磁流体动力学将会受到越来越多的重视,且将会推动航空航天技术的进步。     

高超声速再入钝头体表面热流计算

基于跟踪流线的轴对称比拟法,采用纯工程算法、Euler数值计算与边界层内工程算法相结合的方法,对高超声速再入钝头体的表面热流进行了计算,并将计算结果与风洞实验数据进行了对比,两者吻合较好,验证了两种工程算法在计算高超声速飞行器热环境方面的正确性.将两种工程算法与数值求解N-S方程进行对比,表明工程算法在迎风面的热流计算方面有较高的精度,节约了计算时间,很好地满足高超声速飞行器概念研究和设计的需要.      

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

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

乘波构型飞行器磁流体进气道一体化概念设计

采用准一维的磁流体方程组对设计状态下内磁流体发生器进行参数设计,利用斜激波理论设计磁流体进气道一体化的外形.以该外形为基点,对安装五对电极的内磁流体发生器采用五方程模型进行三维数值模拟,并通过调节其磁作用数以满足准一维假设下设计的进气道外形.研究表明,磁流体进气道最小出口流动均匀参数可以作为判别内发生器最优磁作用数的准则,在该磁作用数下内磁流体发生器可以有效地提高进气道的性能.      

超声速磁流体管道流动的数值模拟

应用可同时适用于描述磁流体MHD (MagnetoHydroDynamic)强化超燃冲压发动机磁流体发电器与磁流体加速器的,包含磁流体N-S(Navier-Stokes)方程与磁感应方程的可压缩磁流体方程组对超声速磁流体管道流动进行了数值模拟研究,以验证本项工作所提出的对2组方程进行交替迭代的分裂式算法.给出了算例的计算结果,并与无磁场及恒磁场下相同管道的流动计算结果进行了对比分析.     

可重复使用运载器递归积分滑模姿态控制方法

针对存在模型不确定和外部干扰的可重复使用运载器再入段姿态控制问题,提出了一种基于自适应滑模干扰观测器的递归积分滑模控制方法。首先,基于可重复使用运载器再入段姿态运动模型,建立了面向控制的模型;其次,设计了自适应滑模干扰观测器,以精准估计和补偿由模型不确定和外部干扰构成的复合干扰;然后,基于递归思想设计了一种新型递归积分滑模控制器,利用Lyapunov稳定性理论证明了闭环系统的有限时间稳定性;最后,数值仿真结果验证了该方法具有较强的鲁棒性和较快的收敛速度。     

基于ESO的可重复使用飞行器积分滑模控制器设计

考虑参数不确定和外界干扰对姿态控制的影响,开展基于扩张状态观测器(ESO,extended states observer)的可重复使用飞行器姿态控制研究.首先,对飞行器运动模型进行描述,在合理假设条件下,得到面向控制器设计的纵向和横侧向姿态模型;其次,分别针对纵向和横侧向姿态模型设计积分滑模控制器(ISMC,integral sliding mode control);然后,引入ESO对耦合作用、参数不确定及外界干扰进行估计,并在控制器中进行补偿;最后,六自由度仿真结果表明,本研究给出的控制策略能够实现对给定制导指令的稳定跟踪,确保安全再入飞行.     

可重复使用飞行器再入姿态的区间二型自适应模糊滑模控制设计

针对具有强非线性、多变量耦合特性的可重复使用飞行器（RLV），同时考虑模型参数不确定性和外界干扰对飞行器再入姿态跟踪的影响，提出了一种基于区间二型自适应模糊滑模的姿态控制方法。首先，建立飞行器再入动态模型，并基于反步思想将控制模型转化为姿态角和角速率相关子系统。其次，将模型参数不确定性和外界干扰视作子系统非线性项的一部分。再次，采用区间二型模糊系统逼近子系统非线性项，并结合自适应技术和滑模控制方法分别设计虚拟控制量和实际控制量。此外，引入一阶低通滤波器用以处理子系统虚拟控制律。通过Lyapunov方法的分析证明了闭环控制系统的稳定性，且飞行器姿态跟踪误差可收敛于原点附近的小邻域。最后，利用飞行器的数值仿真验证了所设计控制方法能有效跟踪飞行器参考指令，且对外界干扰有较强的鲁棒性。      

基于遗传算法的航天器再入动态终迹圈研究

为实现航天器安全准确地沿最优飞行轨迹再入,对航天器再入飞行过程中任意时刻飞行状态,在满足所有约束条件下,计算航天器所能达到的地面最大着陆区域,以判断该时刻飞行器能否到达预定着陆点。文章采用改进的遗传算法对航天器再入飞行轨迹进行优化,结合再入动态终迹圈计算方法,实现再入动态终迹圈的仿真计算。通过对仿真结果的对比分析,得到再入动态终迹圈及相应再入飞行轨迹的特性,对航天器再入飞行制导及航迹规划具有一定的借鉴意义。     

Application of ADER Scheme in MHD Simulation

The Arbitrary accuracy Derivatives Riemann problem method(ADER) scheme is a new high order numerical scheme based on the concept of finite volume integration,and it is very easy to be extended up to any order of space and time accuracy by using a Taylor time expansion at the cell interface position.So far the approach has been applied successfully to flow mechanics problems.Our objective here is to carry out the extension of multidimensional ADER schemes to multidimensional MHD systems of conservation laws by calculating several MHD problems in one and two dimensions: (ⅰ) Brio-Wu shock tube problem,(ⅱ) Dai-Woodward shock tube problem,(ⅲ) Orszag-Tang MHD vortex problem.The numerical results prove that the ADER scheme possesses the ability to solve MHD problem,remains high order accuracy both in space and time,keeps precise in capturing the shock.Meanwhile,the compared tests show that the ADER scheme can restrain the oscillation and obtain the high order non-oscillatory result.     

空间等离子体压力各向异性对磁场重联的影响

基于二维时变可压缩磁流体动力学模拟,数值研究了等离子体压力各向异性对磁场重联的影响,发现一个小的压力各向异性(P_⊥=1.02P_//)即可大大加速磁场重联的发展,这可能是由于磁镜不稳定性与撕裂模不稳定性共同起作用.在P_⊥P_//的情况下,撕裂模受到抑制,电流片中不能形成大型磁岛.      

A CLASS OF TWO-STEP TVD MACCORMACK TYPE NUMERICAL SCHEME FOR MHD EQUATIONS

In this paper, a new numerical scheme of Total Variation Diminishing (TVD) Mac-Cormack type for MagnetoHydroDynamic (MHD) equations is proposed by taklng into account of the characteristics such as convergence, stability, resolution. This new scheme is established by solving the MHD equations with a TVD modified MacCormack scheme for the purpose of developing a scheme of quick convergence as well as of TVD property. To show the validation, simplicity and practicability of the scheme for modelling MHD problems, a self-similar Cauchy problem with the discontinuous initial data consisting of constant states, and the collision of two fast MHD shocks, and two-dimensional Orszag and Tang‘s MHD vortex problem are discussed with the numerical results conforming to the existing results obtained by the Roe type TVD, the high-order Godunov scheme, and Weighted Essentially Non-Oscillatory (WENO) scheme. The numerical tests show that this two-step TVD MacCormack numerical scheme for MHD system is of robust operation in the presence of very strong waves, thin shock fronts, thin contact and slip surface discontinuities.