一种滤波SST方法在翼型深失速模拟中的应用

 为了提高原始剪切应力输运(SST)湍流模型对于分离流动的求解精度,将大涡模拟(LES)中的滤波因子和SST方程相结合构造出一种滤波SST方法,利用湍流尺度对流场求解区域进行划分,近壁面附近的稳态流动由湍流模型控制,远壁面采用LES方法进行模拟。与传统混合RANS/LES方法相比,该方法的特点是:滤波因子的选取不再依赖于网格尺度,可以有效地降低网格诱导分离现象发生的概率。采用该方法对NACA0021翼型深失速特性进行了仿真研究,对比了非定常雷诺平均Navier-Stokes(URANS)方法和SST-DES方法,从仿真结果可以看出滤波SST方法有效地激活了分离区域的脉动,充分展现了分离的三维特性;同时算例求解结果证明该方法的精度高于URANS方法,与试验结果吻合较好,显示其具有一定的工程应用价值。     

T型尾翼风洞颤振实验保护装置绕流特性分析

采用SST两方程湍流模型,通过求解非定常Navier-Stokes（N-S）方程,对T型尾翼风洞实验流场进行了模拟,分析了保护装置对T型尾翼风洞实验流场的影响,研究了保护装置几何外形和保护装置基座后移距离对流场影响。通过对平尾气动力的分析以及对非定常流场的对比,可以得出：采用NACA系列翼型对基座进行气动整流后,基座两侧局部超声速区显著减小,局部激波减弱甚至消失,流场品质得到改善。且采用NACA0010翼型对基座修形后的结果最理想。随着保护装置基座后移距离的增加,平尾气动力均方根值和波动值先是急剧减小,达到0.85倍平均气动弦长后开始有所增大,在2.45～4.05倍平均气动弦长范围基本不再变化,稳定到单独T型尾翼模型相应系数1倍左右。此结论对T型尾翼风洞颤振实验保护装置设计具有一定的指导意义。     

三维翼面结冰过程及其影响数值研究

为得到三维机翼及翼身尾构型在结冰条件大气中的结冰过程和影响,采用结构化网格对M6机翼和翼身尾构结冰过程进行数值模拟,对飞机结冰过程的数值模拟采用欧拉两相流法求解水滴撞击特性,基于Messinger热力学模型模拟冰形增长过程,对其水滴撞击特性、结冰形状及结冰后对气动特性影响进行研究.分析结果可得:机翼展向,随着弦长变小,结冰量逐渐增多,尤其对大展弦比机翼更为明显.同时机身前端迎风面也是比较容易积冰的区域.     

基于 PIV 速度场测量重构压强场的研究进展

利用PIV测量得到的速度场数据重构空间压强场是一种新颖的压强测量方法。目前国外的一些仿真计算和风洞实验已经证明了该方法的可行性和有效性。本文首先详细介绍了基于PIV速度场测量重构压强场的基本原理--两种压强梯度计算方法（拉格朗日方法和欧拉方法）和两种压强积分方法（平面 Poisson 法和直接空间积分法）,然后从速度场测量、压强梯度计算和压强分布计算3个方面综述了基于PIV速度场测量重构空间压强场的关键技术及相关的研究进展,最后从 PIV 速度测量的改进、参数的优化设置、算法的改进与创新、探索并完善3D压强分布计算、可压缩流动条件下的压强重构技术等5个方面探讨了该方法的发展方向,以期引起国内同行对该技术的高度重视并为其进一步发展提供一定的参考。     

Interaction of single-pulse laser energy with bow shock in hypersonic flow

Pressure sensing and schlieren imaging with high resolution and sensitivity are applied to the study of the interaction of single-pulse laser energy with bow shock at Mach 5. An Nd:YAG laser operated at 1.06 lm, 100 mJ pulse energy is used to break down the hypersonic flow in a shock tunnel. Three-dimensional Navier–Stokes equations are solved with an upwind scheme to simulate the interaction. The pressure at the stagnation point on the blunt body is measured and calculated to examine the pressure variation during the interaction. Schlieren imaging is used in conjunction with the calculated density gradients to examine the process of the interaction. The results show that the experimental pressure at the stagnation point on the blunt body and schlieren imaging fit well with the simulation. The pressure at the stagnation point on the blunt body will increase when the transmission shock approaches the blunt body and decrease with the formation of the rarefied wave. Bow shock is deformed during the interaction. Quasi-stationary waves are formed by high rate laser energy deposition to control the bow shock. The pressure and temperature at the stagnation point on the blunt body and the wave drag are reduced to 50%, 75% and 81% respectively according to the simulation. Schlieren imaging has provided important information for the investigation of the mechanism of the interaction.     

Nb基合金中氧原子扩散及团簇化的第一性原理

利用第一性原理计算,研究了氧原子在Nb中的扩散现象以及施加应变的影响,并对氧原子在Nb-Ti、Nb-Zr合金中的团簇化进行了计算分析。发现在单向应变下,在平行于应变方向上,氧原子的扩散势垒随着应变的增大不断减小,当应变为10%时,扩散激活能从初始的0.92eV降低至0.5eV左右;在垂直于应变方向上,氧原子更容易扩散至被拉长的八面体间隙中,即单向应变下氧原子会在Nb中沿应变方向聚集排列,这会在实际材料中导致应力的集中。计算研究了氧原子在Nb合金中的团簇化问题,发现在纯Nb中,氧原子倾向于均匀分布而不是形成氧原子对;添加合金元素Ti和Zr后,随着氧原子浓度的升高,氧原子对将在合金元素周围形成,Zr周围氧原子对的形成可使体系的能量降低0.29eV。     

基于变步长CFD/RBD方法的旋转弹轨迹仿真

采用变步长方法实现非定常Navier-Stokes方程和刚体六自由度运动方程的耦合求解,并对美国ARL旋转弹试验模型的飞行轨迹进行了仿真。旋转弹自由飞行过程中受到空气阻力作用,转速下降,气动力变化趋缓。因此轨迹仿真时先选取较小的时间步长,随着转速的下降逐渐增大步长,这样可以减少计算时间。刚体六自由度运动方程的数值求解采用变步长Adams预估校正法。将定步长和变步长两种方法应用于旋转弹实验模型轨迹仿真,并将计算得到的旋转弹空间位置、姿态角与试验结果进行了对比,结果显示:方法能较精确地仿真旋转弹飞行轨迹。在相同初始步长条件下,变步长方法较定步长方法计算时间减半;在相同计算时间条件下,变步长方法较定步长方法精度更高。     

控制系统中的矩阵方程的数值解

系统分析与设计中遇到的关键问题之一是求解某些矩阵微分方程和矩阵代数方程,而它们的求解往往是困难的。本文就这些方程的解法总结了近十几年来国内外的一些较好妁数值方法,並给出了数例。     

彗尾剪切磁流体的修正B—KdV方程及其非线性特征

本文运用二维约化摄动方程推导出彗尾剪彻磁流体的修正B-KdV方程,并定性讨论了其非线性性态。主要结论是:(1)推导出彗尾剪切磁流体修正的B-KdV方程,(2)给出彗尾宏观非线性波动与剪切背景的制约关系,(3)定性地解释了彗尾等离子体的一些复杂结构,(4)给出太阳风与彗尾等离子体相互作用的一种输运机制。     

圆楔体轴向共形缝隙阵的分析

本文提出一种严格分析圆楔体轴向共形缝隙阵的方法。在直接求解法的基础上,运用矩量法对空腔探针馈电方式的缝隙阵的互耦进行了分析,并得到了辐射场量表达式。最后举例对圆楔体轴向共形二元缝隙阵的方向图进行了数值计算并由实验证实其结果。