对流场自适应的确定性涡方法

本文系统地研究了随机走步涡团法所面临的主要问题，以统一的途径构造了自适应涡团模型，粘性扩散的确定性算法及确定边界涡量的新方法，数值结果是令人满意的。     

用确定性涡方法模拟二维钝体绕流

以离散涡方法为基础,采用粒子强度交换法模拟涡量扩散方程,对粒子场分布产生扭曲变形时采用粒子强度重分配,数值模拟了突然起动圆柱、风攻角为的方形截面和大海带东桥的流场,计算结果与解析解及其他数值方法计算的结果均吻合的很好,数值结果显示确定性涡方法的高精度的特性使其在模拟钝体扰流问题方面具有很大的优势.     

基于黏性涡域法的二维钝体绕流模拟

针对传统计算流体力学方法需要划分网格且计算耗时长等问题，一种确定性涡方法(黏性涡域法)逐渐发展并日臻成熟。本文利用该方法计算模拟了二维钝体的流场、气动系数及斯特劳哈尔数，并与其他成熟CFD软件计算结果对比来验证模拟精度。黏性涡域法的模拟过程为生成涡量、涡量演变及气动计算。文章首先推导了该方法各过程的基本方程，然后基于Lua语言编制了圆柱和方柱绕流计算程序，最后利用Gnuplot实现流场可视化。利用黏性涡域法模拟不同雷诺数下圆柱绕流、不同无量纲频率和无量纲振幅的振动圆柱气动力时程、不同风向角下的方柱绕流，结果表明：雷诺数为2×10²～1×10⁵时，基于黏性涡域法的圆柱绕流气动计算结果与其他方法的计算结果基本一致；基于该方法的不同风向角下方柱绕流气动计算结果与其他文献结果趋势相同。     

部件性能非确定性对涡轴发动机性能影响量化方法研究

航空发动机各部件性能非确定性的累积叠加会导致发动机总体性能的非确定性。建立了基于蒙特卡洛概率设计方法的非确定性涡轴发动机数学模型,用以量化发动机各部件性能的非确定性对发动机性能的影响。结果表明:各部件性能按正态分布进行非确定性变化时对发动机性能有明显影响,与发动机确定性模型所计算的性能100%达标相比,输出功率的达标概率为49.32%,耗油率的达标概率为50.83%,模型量化结果合理且均方根误差小于3%。发动机在采用限制涡轮前温度的控制规律后,输出功率和耗油率的达标概率进一步降低。     

非定常动静叶干涉作用的二维数值模拟

采用涡方法研究叶轮机内动、静叶间的相互作用,它利用全三维的、通道时间平均的定常解为基础解,并由此给出非定常扰动场的初始解。为计算叶片对扰动场的响应过程,采用拉格朗日方法追踪扰动涡团的对流流动过程,用确定性涡方法来描述流体的粘性扩散过程。用该方法对某型压气机的一级（动、静叶栅距比为１．７４）的动、静叶相互作用,进行了二维有粘问题的数值模拟。数值结果表明,该方法具有较好的收敛性和收敛速度;转子出口的尾迹和各种旋涡等非定常因素在静叶通道中的切割、迁移过程,非定常性对静叶攻角的影响以及脉动压力的振幅衰减特性也在文中作了分析。     

嵌入式多核操作系统确定性研究

为解决嵌入式多核操作系统确定性问题,详细阐述了计算确定性,重点分析了多核结构引入的不确定性问题。通过研究嵌入式多核操作系统的结构,从共享存储的确定性访问、确定性任务调度、确定性通信等方面提出了解决多核结构带来的不确定性问题的方法。     

关于涡量的动力学边界条件

本文讨论了以涡量为变量的流体动力学描述方法中涡量的边界条件问题，并且用有关的数值计算结果对应用各种边界条件的结果进行了验证。     

守恒型可压缩一维湍流方法及其在超声速标量混合层中的应用

一维湍流（ODT）方法是一种能在一维计算域上遵循湍流基本物理规律的湍流建模方法。通过结合确定性和随机性求解方法,能够在一维计算域上准确捕捉到湍流统计规律,且降维建模可显著减小计算量。ODT方法主要被广泛用于不可压湍流和湍流燃烧研究,若要将其拓展用于模拟高速可压缩湍流,需对建模方法进行深度改进。相比于不可压ODT方法,本文基于欧拉参考框架,针对可压缩湍流的特性,将因变量由原始变量改为有利于减小可压缩湍流模拟误差的守恒通量,并加入了组分求解模块。对确定性和随机性求解模块均进行了相应的深度改进,开发出具有标量混合模拟功能的守恒型可压缩ODT方法。在确定性模块中改为求解以守恒通量为变量的一维截断控制方程,在随机性模块中构造一维涡时,将三联映射的作用对象也相应地由原始变量改为守恒通量,并选用了可保证变密度情况下动量守恒的双核变换。通过模拟空间发展超声速平面湍流混合层并将自相似阶段结果与实验结果比对,验证该方法对可压缩剪切湍流场中标量混合的捕捉精度。守恒型可压缩ODT方法模拟得到的速度场和组分场的平均剖面和脉动强度分布与实验结果准确吻合,精度明显优于传统的耦合梯度扩散亚格子模型的大涡模拟方法（LES-GRAD.DIFF.）以及耦合线性涡（LEM）亚格子模型的大涡模拟方法（LES-LEM）,且该方法的降维处理使其在降低计算成本方面具有显著优势。     

某型涡轴发动机性能衰减与部件退化评估

为在某型涡轴发动机持久试验中进行性能健康检查,采用多点匹配方法建立了个体发动机性能衰减计算模型,基于非线性气路分析方法对部件性能退化进行了计算,集合已有经验确定性能偏离部件,完成了某型发动机持久试验中不同阶段的性能衰减评估.结果表明,某型发动机性能衰减及部件性能退化成阶段性分布,前期压气机和燃气涡轮性能变化最快,随后处于相对稳定阶段,动力涡轮性能在目前试验时数范围内基本无变化.该方法的验证为涡轴发动机使用过程中的性能健康与评估提供了有益的经验.      

模拟有限长涡片旋转的一种新数值模型

本文把涡格法中的超收敛性推广到自由涡片的演化问题，涡片的演化由ＢｉｒＫ－ｈｏｆｆ－Ｒｏｔｔ方程所控制，该方程具有强烈的非线性，难以用解析法求解。在数值模拟中通常离散涡方法。本研究表明对自由涡片的旋转问题在精确解和离散涡方法结果之间的误差主要来自点涡近似。使用基于超收格式的新数值模型可以显著地提高数值解的精度。     

关于涡方法的讨论

本文讨论了如何提高二维涡方法的计算精度和减小涡方法计算量的问题。文中指出了改进的途径,同时列举了当前这些方面的研究进展情况。     

尾涡数值模拟的曲涡元技术

连接涡线端点形成涡段并计算其对空间任意点的影响系数是尾迹分析的基本环节,称为涡元技术。本文推出圆弧涡元与其应用程序,并引进抛物拱弧涡元,以代替传统的直涡元。提供的曲涡元程序可方便地用于尾迹分析,算例表明曲涡元在提高精度减少计算量方面的先进性。     

应用 TVD 格式求解 N-S 方程预测透平叶栅二次流损失

应用ＴＶＤ格式和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数紊流模型求解三维Ｎ－Ｓ方程,详细模拟出了一个直列透平叶栅中二次流及其损失在叶栅通道中的产生和发展过程,并给出了叶栅损失的分布状况,其结果深刻揭示了透平叶栅二次流损失的产生发展机制。本文计算得到的叶栅损失的栅距平均值沿叶高的分布曲线与实验值吻合良好,表明格式求解精度较高。     

用于旋翼气动特性计算的一种时间推进自由涡新方法

建立了一个新的用于旋翼尾迹和诱导速度分析的时间推进自由涡计算方法.在该方法中,推导了一个新的二步二阶预测-校正差分算法--"D2PC"算法,并对采用该算法离散表示的涡线主控方程进行修正,以提高尾迹的求解精度;基于Vatistas公式给出涡核模型,并考虑了实际的耗散影响;采用更适合于旋翼气动特性计算的简化升力面桨叶气动模型;同时在自由涡方法中结合了挥舞动力学模型和旋翼配平模型.通过计算的旋翼尾迹几何形状和下洗速度与实验值的对比,验证了本文建立的自由涡方法的有效性.然后利用该方法,计算和分析了在悬停和前飞状态旋翼的尾迹几何形状和下洗诱导速度.     

近地低速飞行时旋翼尾涡系的畸变及其诱导速度

本文用一简化自由涡模型模拟旋翼尾涡系,用一简单涡模拟地面涡,对近地低速飞行时旋翼尾涡系的畸变及其引起的诱导速度变化作了研究。在地面涡和地面的影响下,桨盘附近的尾涡几何位置发生很大变化,导致桨盘处的诱导速度发生很大变化。尾涡畸变在桨盘处引起的诱导速度变化远远大于地面涡和地面在桨盘处直接产生的诱导速度。考虑尾涡畸变后计算出的诱导速度在桨盘处的分布和实验符合很好。     

基于Kirchhoff方法的亚声速平面混合层主涡对并声场分析

为拓展主涡对并气动声场计算分析方法研究,针对亚声速平面混合层主涡对并流动,采用直接数值模拟(DNS)方法求解获取近场声源,利用二维Kirchhoff频域法和时域法外推声辐射远场,并应用Kirchhoff频域法对主涡二次对并声辐射特征进行分析。研究结果表明:Kirchhoff频域法和时域法具有一定差异性,但与主涡对并直接声计算结果对比,显示二者计算精度相当且均具有较高的精确性。通过Kirchhoff频域法应用分析,清晰揭示出主涡二次对并各声模态与主涡卷起、对并过程的对应相关性,利用频域法这种良好的计算分析性能,本文进一步研究了混合层入口扰动相位差对涡并声场的影响,展现出主导声波在声场中的重要作用。     

增升装置微型涡流发生器数值模拟方法研究

 针对数值模拟带微型涡流发生器(VG)的增升装置绕流过程中出现的典型问题,研究相应解决途径.引入面搭接网格技术,避免传统结构网格方法造成的网格数量过大问题;在微型涡流发生器网格上采用特殊边界条件处理以分别对应有/无涡流发生器的状态,排除网格变动对计算结果造成的干扰;研发出了面搭接网格生成技术与特殊边界处理方法结合的求解加装微型涡流发生器的增升装置流动数值模拟方法.研究结果表明:该方法在处理增升装置这类复杂外形、复杂流动问题以及加装涡流发生器的流动控制问题时,具有良好的数值计算精度,能够满足研究需要,为开展涡流发生器精细设计和参数优化的数值模拟研究提供了可靠的技术途径.     

亚声速全机纵横向气动特性计算方法

本文介绍了一种涡格方法,它适用于亚声速飞机的纵横向气动特性计算。该方法中,升力面用布置在中弧面的平均平面上的单涡格层,或用布置在上、下表面的平均平面上的双涡格层模拟,机身则由布置在表面上的四边形涡环或由布置在一组同轴棱柱面上的涡栅来模拟。地面效应用布置镜像涡系的办法处理。文中讨论了机翼尾涡偏斜对平尾、立尾的影响。 该方法计算的结果与有关的理论和实验进行了比较。这些比较表明,此方法计算精度满足工程设计的要求。     

Generalized Kutta–Joukowski theorem for multi-vortex and multi-airfoil flow(a lumped vortex model)

For purpose of easy identification of the role of free vortices on the lift and drag and for purpose of fast or engineering evaluation of forces for each individual body,we will extend in this paper the Kutta–Joukowski(KJ) theorem to the case of inviscid flow with multiple free vortices and multiple airfoils. The major simplification used in this paper is that each airfoil is represented by a lumped vortex,which may hold true when the distances between vortices and bodies are large enough. It is found that the Kutta–Joukowski theorem still holds provided that the local freestream velocity and the circulation of the bound vortex are modified by the induced velocity due to the outside vortices and airfoils. We will demonstrate how to use the present result to identify the role of vortices on the forces according to their position,strength and rotation direction. Moreover,we will apply the present results to a two-cylinder example of Crowdy and the Wagner example to demonstrate how to perform fast force approximation for multi-body and multi-vortex problems. The lumped vortex assumption has the advantage of giving such kinds of approximate results which are very easy to use. The lack of accuracy for such a fast evaluation will be compensated by a rigorous extension,with the lumped vortex assumption removed and with vortex production included,in a forthcoming paper.     

基于高阶WENO格式的旋翼非定常涡流场数值模拟

为提高直升机旋翼黏性涡流场计算流体力学（CFD）的计算精度及降低数值耗散,发展了一套基于五阶加权本质无振荡（WENO）格式的旋翼非定常涡流场数值计算方法。采用运动嵌套网格方法生成围绕前飞状态旋翼的网格系统,以Navier-Stokes方程作为主控方程,湍流模型采用一方程Spalart-Allmaras模型。为了提高旋翼流场中对涡的形成、演化等发展过程的模拟精度,采用高间断分辨率的Roe-WENO格式计算对流通量。对状态变量的插值,选取适当的加权因子,通过多个重构模板的凸组合来构造五阶WENO格式,在交界面附近获得具有五阶精度的状态变量。为了提高计算效率,采用高效的隐式LU-SGS（Lower-Upper Symmetric Gauss-Seidel）推进方法和并行计算策略。最后,运用该方法分别对悬停状态下C-T（Caradonna-Tung）旋翼和UH-60A旋翼以及前飞状态下C-T旋翼和SA349/2旋翼的涡流场及气动特性进行了数值模拟,并将数值结果与传统二阶精度格式的计算值进行了对比。结果表明：在计算时间增长10%～20%的前提下,五阶WENO格式能够捕捉到更精确的涡流场特性,反映了五阶WENO格式在计算旋翼非定常涡流场时具有更高的计算精度与更低的数值耗散特性。