旋翼表面非定常压力脉动计算的三维自由尾迹非定常面元法

利用最新发展的高分辨率三维自由尾迹非定常面元方法计算多叶片直升机旋翼的三维非定常气动力，将风洞实验数据与旋翼在悬停和出现严重桨涡干扰的低速前飞下滑状态下的计算结果进行了比较，在悬停时，还与欧拉方法的结果进行了比较，本方法与欧拉方法的结果以及风洞实验结果都吻合很好。     

基于椭圆积分的单支点半柔壁喷管大挠度成型半解析求解与验证

为提升跨超声速风洞单支点半柔壁(SJSF)喷管分析精度和效率,发展了基于椭圆积分的柔性壁板大挠度变形半解析求解方法.以0.3 m低温风洞单支点半柔壁喷管为研究对象,建立型面组件力学模型,基于欧拉梁理论推导柔性壁板型面曲线的微分方程,给出了椭圆积分的表达形式.利用Maple软件编制了计算程序,并以此针对马赫数为1.15和...     

跨声速三维非线性洞壁干扰的数值计算

分别以固壁条件和洞壁附近的压力分布模拟各类实壁和透气壁试验段的洞壁边界条件，利用Ｅｕｌｅｒ方程和Ｎ－Ｓ方程数值求解模型在风洞中的绕流场，得出洞壁干扰对跨声速模型绕流和气动力的影响，初步的研究结果表明，该方法能较有效地模拟模型在跨声速风洞中的绕流场，经洞壁干扰修正后的ＧＢＭ－０４Ａ模型在０．６ｍ风洞中的试验结果与无干扰参考结果吻合较好     

战斗机弹射救生系统数值模拟研究

对弹射救生系统这种几何外形复杂、不规则非流线型钝体,以Jameson有限体积法为基础发展了一套基于高质量混合网格的Euler方程解算器,在Euler方程解算器基础上,采用S-A湍流模型发展了一套同样基于混合网格的N-S方程解算器,初步建立了适合弹射救生系统这类几何外形复杂物体的数值模拟技术平台,并采用多种加速收敛措施和合适的耗散项模型使得所发展的解算器在软件实现上鲁棒性很强。通过选取合适的计算起始迎角,对弹射救生系统成功地进行了大迎角和大侧滑角绕流计算,获得了与风洞试验比较吻合的计算结果,并且可以清楚地观察流场中马赫数分布、弹射救生系统各部位压力分布等流场细节,对弹射救生系统空气动力特性分析具有重要的指导意义。该计算平台可作为风洞试验的合理补充、加强和完善。     

结冰风洞中液滴相变效应数值模拟

为明晰结冰风洞中液滴相变效应,发展了基于Euler法的气液两相传质传热耦合流动计算方法,探索了3 m×2 m结冰风洞主试验段构型相变效应对液滴传热过程的影响,开展了参数影响研究,评估了试验段内液滴过冷状态。结果表明：构型内液滴经历了先蒸发后凝结两个阶段,蒸发效应促进了准一维传热阶段液滴温度的下降趋势,使液滴温度趋于湿球温度,而凝结效应则抑制了三维收缩阶段液滴温度的下降趋势,进而增大了试验段液气温差,影响液滴过冷状态;增大初始相对湿度和试验段气流速度,会导致蒸发效应减弱而凝结效应增强,进而增强了相对湿度和气流速度对液滴过冷状态的影响程度,与此相反,初始液滴尺寸的增加,则会导致蒸发效应和凝结效应均减弱,进而减弱了液滴尺寸对液滴过冷状态的影响程度;在典型工况下,小尺寸液滴（液滴直径范围为40~100 μm）在高风速时（试验段气流速度大于164 m/s）将偏离过冷状态（液气温差超过3℃）。     

双欧法与四元数法的应用比较

对解决欧拉方程奇异性的双欧拉和四元数法进行了对比研究，四元数法从理论上讲比较完美，但实际应用中存在较大的累积计算误差，从而影响计算精度；双欧法利用正、反欧拉方程间精华区倒挂关系进行分区交替运算，把精华区扩展到全域，不仅根除了奇异性，而且计算误差小。因此，对于解决欧拉方程奇异性来讲，双欧法要优于四元数法。     

基于欧拉方程的二维振荡机翼非定常气动设计反命题方法

建立了基于欧拉方程组的二维振荡机翼非定常气动设计反命题方法的数学模型.通过一系列变换,将物理时空中的求解域转换成映象坐标系中的规范区域;并导出了映象坐标下的欧拉方程,结合非定常反命题的边界条件,便可用有限差分方法求解.本方法引入映象坐标系解决反命题边界形状的不确定性,并能利用欧拉方程现有的各种差分格式于反命题求解.     

Experience in the use of computational aerodynamics to predict store release characteristics

In the early days, store separation tests were conducted in a hit or miss fashion—the stores would be dropped from the aircraft at gradually increasing speeds until the store came close to or sometimes actually hit the aircraft. In some cases, this led to loss of the aircraft, and made some test pilots reluctant to participate in store separation flight test programs.During the 1960's, the Captive Trajectory System (CTS) method for store separation wind tunnel testing was developed. The CTS provided a considerable improvement over the hit or miss method, and became widely used in aircraft/store integration programs prior to flight-testing. However, since fairly small-scale models had to be used in the wind tunnel tests, in many cases the wind tunnel predictions did not match the flight test results. No mechanism was then in place to resolve the wind tunnel/flight test discrepancies.During this same time frame Computational Fluid Dynamics (CFD) had finally matured to the point of providing a trajectory solution for a store in an aircraft flowfield. However, since the computational tools were necessarily (due to computer resource limitations) limited to linear techniques, and since most store separation problems occur at transonic speeds, these tools had limited application.Recent advances in computer resources have greatly improved the capability of CFD to predict store release characteristics. Instead of using linear or approximate schemes, time dependent Euler and Navier Stokes trajectories could be computed in a reasonable time frame.Three international CFD Challenges, held during the last decade of the 20th century, have shown that CFD can not only match wind tunnel test data, but also predict flight test trajectories for complex stores at transonic speeds. It appears that CFD has matured to the point that it can be usefully integrated into aircraft/store compatibility programs.     

螺旋桨流场Euler/NS数值解对比分析

采用多块点对点网格生成技术,生成螺旋桨流场空间计算网格。针对不同前进比,分别求解Euler和NS方程对螺旋桨流场进行数值模拟,并与实验结果进行了对比分析。研究发现无粘和有粘的计算结果均接近于实验结果,因此基于Euler方程的流场数值模拟更适用于工程设计。     

基于CFD的带控制舵导弹的动导数计算

介绍了基于计算流体动力学(CFD)的阻尼导数预测方法以及其在复杂带控制舵导弹外形的俯仰/滚转阻尼导数预测中的应用.方法采用求解三维非定常ALE(arbitrary Lagrangian-Eulerian)形式RANS(Reynolds-averaged Navier-Stokes)方程加湍流模型的方法得到模型强迫运动的气动力时间历程,再通过绕模型非定常运动过程得到的气动力迟滞曲线计算得到俯仰及滚转阻尼导数.为了验证采用数值方法的正确性,针对基本带翼导弹外形,采用强迫模型按正弦形式振动,计算得到了超声速范围内马赫数从1.5到2.5的俯仰及滚转稳定性参数,通过与风洞试验及弹道数据结果的比较发现,两者吻合较好.      

带控制舵弹体气动特性分析

本文用数值模拟方法研究了带控制舵弹体的气动特性。该问题外形复杂，在控制舵和弹身之间会出现明显的缝隙。为了解决此类外形的气动力计算问题，我们发展了结构网格和非结构网格的混合网格技术。这样一方面充分发挥了贴体结构网格质量好、效率高的优点，另一方面充分发挥了非结构网格几何灵活性好、易于模拟复杂外形的优点。数值计算与风洞实验相当吻合，说明本方法具有良好的精度。最后对两种外形进行大量计算并进行气动特性分析，给出了静稳定特性和配平特性的研究结果。     

HB-2 high-velocity correlation model at high angles of attack in supersonic wind tunnel tests

Responding to a need for experimental data on a standard wind tunnel model at high angles of attack in the supersonic speed range, and in the absence of suitable reference data, a series of tests of two HB-2 standard models of different sizes was performed in the T-38 trisonic wind tunnel of Vojnotehnickˇi Institut(VTI), in the Mach number range 1.5–4.0, at angles of attack up to+30°. Tests were performed at relatively high Reynolds numbers of 2.2 millions to 4.5 millions(based on model forebody diameter). Results were compared with available low angle of attack data from other facilities, and, as a good agreement was found, it was assumed that, by implication, the obtained high angle of attack results were valid as well. Therefore, the results can be used as a reference database for the HB-2 model at high angles of attack in the supersonic speed range, which was not available before. The results are presented in comparison with available reference data, but also contain data for some Mach numbers not given in other publications.     

马赫数2～4连续可调风洞数值模拟及静态标定试验

型面旋转连续可调风洞可以获得出口马赫数连续变化的出口流场,其结构简单、易于调节、响应迅速,成为目前国内外的研究热点,具有很好的发展前景。本文针对某优化设计的出口马赫数2~4的变马赫数风洞喷管,进行了全流场三维数值仿真校核,并进行了风洞流场静态标定试验,获得了该连续可调风洞在不同马赫数下的流场品质及流场均匀区大小。结果显示:该连续可调风洞在马赫数2~4下的出口流场均匀性良好,流场品质满足固定几何风洞流场的国军标要求;试验得到的不同马赫数工况下均匀区大小均大于理想菱形区的2/3,均匀区马赫数标准差在马赫数2~3工况下小于0.01,在马赫数3.5工况下为0.013 6,说明该连续可调风洞可以实现马赫数的连续变化,且具有良好的流场品质,可用于后续的风洞试验。     

二维折角管道内超声速流动的多重网格计算

本文应用Euler模型和隐式TVD格式计算了二维折角管道内的超声速流动问题,为了有效地模拟激波和膨胀波系的相互干扰及壁面的反射,本文试将多重网格法和TVD格式结合起来进行数值求解。考虑到目前绝大部分文章采用V-循环多重网格法,我们对多重网格法中的几种多层次循环形式(V-循环、W-循环)进行了比较性的数值计算,旨在从中发现合适的多层次循环求解的形式,以期充分发挥多重网格法在Euler方程求解过程中的效益。同时,本文也是对多重网格法和TVD格式结合求解的一次初步尝试。     

小迎角轴对称体跨声速绕流的一种Euler解法

应用小迎角展开的方法将小迎角下旋成体绕流Ｅｕｌｅｒ方程进行摄动展开，对零介和一阶方程分别采用隐式ＮＮＤ格式进行了计算，与实验进行比较得到了较为满意的结果。该方法采用二维的计算得到了三维的结果，大大提高了计算的效率。     

利用转移矩阵进行姿态控制

传统欧拉角控制方法在处理大姿态角机动时往往出现奇异和交联问题。针对这一问题对欧拉角的跳变现象进行了讨论,并从飞行器实际姿态与需求姿态之间的空间关系入手,给出了一种消除欧拉角跳变现象的有效方法－转移矩阵法,仿真计算表明,该方法在大型航天器的姿态控制中可以有效地消除欧拉角奇异、交联和跳跃现象,在欧拉角跳变的情况下仍然能够实现对飞行器的控制。     

解决欧拉方程奇异性的方法探讨

提出了一种解决欧拉方程奇异性的工程算法,从飞行器运动方程入手,推导出具体计算公式,并进行了数值计算验证。该方法无论是理论分析还是工程计算都证明是可行的,没有发现原理缺陷和方法误差,不会产生奇异性,在飞行器仿真、模拟和姿态控制中具有普遍的应用和参考价值。     

高升力下二维构型阻力准确测量研究

翼型风洞试验阻力测量常使用尾迹流场测量积分求取阻力的方法,但各积分公式均建立在一定的假设基础上,有一定适用范围。在多段翼型流场N-S方程数值模拟和风洞试验的基础上,研究高升力情况下低速风洞阻力精确测量技术。通过N-S方程数值模拟求解多段翼型绕流场,分析尾迹流场的特点和常规风洞试验阻力计算公式推导时所作假设,提出新的更为准确的型阻计算公式;利用多段翼型绕流的数值模拟结果,积分表面压力和摩擦力求得翼型的气动特性,并利用计算得到的尾迹流场信息按照常规和新提出的风洞试验型阻计算公式计算阻力,将三者进行比较,检验提出的新型阻计算公式的准确性;通过风洞试验检验数值模拟得到的流场特点和新型阻计算公式。研究表明：在高升力条件下,传统型阻计算公式有很大的局限性,必须进行改进;提出的考虑尾迹区流动特点的新型阻计算公式能够得到更准确的阻力值。     

大型煤气柜风荷载的风洞试验及数值模拟

对某大型煤气柜进行了风洞测压试验及风压数值模拟。分析了试验模型表面风压分布及其脉动特性,并同数值计算结果、规范条文中类似断面结构的风压分布作对比。结果表明：风洞试验中由于结构表面分布的工字钢及表面粗糙度的处理,雷诺数效应对表面风压分布影响并不明显,但对表面绕流场分离区的风压值有一定影响。结构的均方根升力、阻力系数在频域表现为宽带谱;采用基于雷诺平均的RNGk-ε湍流模型能较准确地模拟表面平均风压分布,其计算结果同样可为结构抗风设计提供参考;在不同量级雷诺数下数值模拟得到的平均风压分布能反映出雷诺数效应的影响。