脉冲等离子体推力器羽流预测计算研究

为了能够将脉冲等离子体推力器成功地运用于空间,需对其羽流进行研究。将一维 MHD双温放电模型的计算结果作为入口条件,运用DSMC（Direct Simulation Monte\|Carlo ）／PIC（Particle in Cell）流体混合算法一体化模拟实验室PPT羽流。验证计算显示该模 型具有一体化模拟脉冲等离子体推力器羽流的能力。对不同初始放电能量下的羽流场进行模 拟,给出了离子、中性粒子、电子温度、轴线上质量流率和出口平面返流质量流率的变化情 况。计算结果显示高放电能量下返流量更大,同时中性粒子在返流中所占比例也越大。
     

Competing effects of surface catalysis and ablation in hypersonic reentry aerothermodynamic environment

Under hypersonic flow conditions, the complicated gas-graphene interactions including surface catalysis and surface ablation would occur concurrently and intervene together with the thermodynamic response induced by spacecraft reentry. In this work, the competing effects of surface heterogeneous catalytic recombination and ablation characteristics at elevated temperatures are investigated using the Reactive Molecular Dynamics(RMD) simulation method. A GasSurface Interaction(GSI) model is establi...     

旋翼桨尖涡生成及演化机理的高精度数值研究

为了细致捕捉直升机旋翼桨尖涡的生成和演化过程,建立了一个基于高精度网格和高阶通量计算格式的旋翼桨尖涡计算流体力学(CFD)求解方法。在该方法中,流场求解选取旋转坐标系下的Navier-Stokes方程为控制方程;空间离散采用迎风Roe格式,并采用低耗散的5阶WENO(Weighted Essentially Non-Osciltatory)格式进行对流通量的计算;时间推进则采用双时间法,在伪时间步上使用隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)推进格式;应用嵌套网格方法实现桨叶网格和背景网格的数据交换。应用所建立的方法对悬停状态的旋翼桨尖涡流场进行了高精度模拟,在桨叶网格上精细地捕捉到了桨尖涡的具体生成过程,在背景网格上捕捉到了更多圈数的桨尖涡尾迹,并对桨尖涡的演化机理进行了深入研究。结果表明:建立的高精度数值方法能够有效地对旋翼桨尖涡的生成和演化过程进行细致模拟;悬停状态下旋翼桨尖气流在上下表面压力梯度的作用下经历了边界层增厚、逐渐卷起形成涡核以及最终脱离桨叶形成桨尖涡的过程。     

一种基于求解椭圆型方程的结构动网格生成方法

基于椭圆型网格生成法,实现了一种简单高效的贴体结构动网格生成方法,可用于具有移动边界问题的非定常流动数值模拟。该方法提出,在网格变形过程中,Poisson方程需要的控制网格间距和正交性的源项可以通过提取已知的静态网格源项直接得到,并在整个动网格生成过程中保持不变。因此,在椭圆型网格生成中需要通过外迭代确定源项的过程可以得到省略,而且该方法不需要人工指定参数。这使得方法具有高效和易于嵌入到已有程序中的特点。数值模拟结果证明,采用这种方法获得的网格能够较好地保持静态网格原有的正交性和光滑性,在相同迭代步数约束下,网格求解效率低于传统弹簧模拟法,但鲁棒性优于弹簧模拟法。     

超微型高速空气轴承的动态特性

对微型转子（半径为1985mm）动态特性进行了数值模拟和实验研究。在数值模拟中,利用ANSYS对空气轴承 微转子系统进行了双向瞬态流固耦合研究;模拟了启动阶段转子在空气轴承调节下的运动过程。结果表明:转子的动态过程为一运动幅度逐渐减小的振荡过程,最后将稳定于某一平衡位置;在轴承供压越大时,转子运动幅度越大,且达到稳定所需时间更长。计算结果发现了轴承间隙宽度与轴承流量的正相关关系。在对应实验中,实验结果与模拟结果很好的相符。由数值模拟和实验研究,说明了在较大止推轴承供压下由于转子轴向运动幅度大,有更大发生碰撞的风险;而在70927Pa的供压下进行启动则更为合适。      

考虑配平特性的超声速客机低声爆气动布局优化研究

降低声爆水平是下一代超声速运输机研制需要解决的关键问题之一。低声爆优化通常使飞行器布局向着机翼后掠角增大、机翼沿机身方向分布范围增大的趋势发展,给飞行器的配平和低速特性带来不利影响。以某超声速客机基本构型为研究对象,建立基于类别/形状函数的翼身组合体参数化建模方法;基于超声速线化理论分析外形几何参数对声爆水平的影响。在此基础上,分别针对机身轮廓、机翼平面形状以及扭转角分布对该构型进行低声爆优化和俯仰力矩特性优化,并采用CFD 方法对优化结果进行校核。结果表明：与基准构型相比,在不显著增加俯仰力矩的基础上,优化构型的阻力降低了19 cts,近场过压显著降低,地面声爆响度降低5.1 PLdB。     

Numerical Simulation of Rotor Flow Field Based on Overset Grids and Several Spatial and Temporal Discretization Schemes

A numerical method based on solutions of Euler/Navier-Stokes (N-S) equations is developed for calculating the flow field over a rotor in hover. Jameson central scheme, van Leer flux-vector splitting scheme, advection upwind splitting method (AUSM) scheme, upwind AUSM/van Leer scheme, AUSM+ scheme and AUSMDV scheme are implemented for spatial discretization, and van Albada limiter is also applied. For temporal discretization, both explicit Runge-Kutta method and implicit lower-upper symmetric Gauss-Seidel (LU-SGS) method are attempted. Simultaneously, overset grid technique is adopted. In detail, hole-map method is utilized to identify intergrid boundary points (IGBPs). Furthermore, aimed at identification issue of donor elements, inverse-map method is implemented. Eventually, blade surface pressure distributions derived from numerical simulation are validated compared with experimental data, showing that all the schemes mentioned above have the capability to predict the rotor flow field accurately. At the same time, vorticity contours are illustrated for analysis, and other characteristics are also analyzed.     

卫星用两种单相流体回路方案对比研究

单相流体回路系统是航天器主动热控技术中的一种重要方式。文章概述了流体回路的工作原理及理论基础，总结了单相流体回路的设计思路，结合某卫星型号，设计了两种流体回路方案，在SINDA／FLUINT中建立数值仿真模型进行分析计算，并比较了两种方案的应用性能。结果显示，将单相流体回路系统应用于卫星，可取得较好的控温效果。     

可压缩流动激波装配在格心型有限体积法中的应用

发展了一种基于格心型有限体积方法(FVM)的激波装配算法。通过定义网格节点属性可以灵活调用激波装配和激波捕捉计算方法。在使用激波装配方法时,激波节点运动速度和下游运动速度通过Rankine-Hugoniot(R-H)关系式获得,同时采用非结构动网格技术描述激波的运动以及调整其他网格节点的位置。流过激波面元的通量为上游单元的基本通量,物理概念更加清晰,通量计算也更为准确。在计算过程中,网格节点属性可以发生变化,以此实现对带有拓扑变化流场的描述。数值试验表明:本文提出的计算方法不但具有较高的计算精度,同时能有效地避免由于捕捉激波而出现的数值问题。     

Numerical prediction and wind tunnel experiment for a pitching unmanned combat air vehicle

The low-speed flowfield for a generic unmanned combat air vehicle (UCAV) is investigated both experimentally and numerically. A wind tunnel experiment was conducted with the Boeing 1301 UCAV at a variety of angles of attack up to 70 degrees, both statically and with various frequencies of pitch oscillation (0.5, 1.0, and 2.0 Hz). In addition, pitching was performed about three longitudinal locations on the configuration (the nose, 35% MAC, and the tail). Solutions to the unsteady, laminar, compressible Navier–Stokes equations were obtained on an unstructured mesh to match results from the static and dynamic experiments. The computational results are compared with experimental results for both static and pitching cases. Details about the flowfield, including vortex formation and interaction, are shown and discussed, including the non-linear aerodynamic characteristics of the vehicle.