高超声速滚转阻尼导数数值模拟

采用非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程描述物体简谐振动流场,并在Ｅｔｋｉｎ理论下给出绕定轴转动时滚转阻尼导数的计算公式,定常流场的计算采用空间二阶精度的交替方向隐式分解的ＮＮＤ格式,非定常流场的计算采用时、空二阶精度的Ｒｕｎｇｅ－Ｋｕｔｔａ多步格式,采用代数和方法生成物体静、动网格。最后给出高超声速下钝体外形滚转阻尼导数的计算结果,以及滚转力矩系数随瞬时振幅的变化曲线。     

三角翼匀速上仰跨声速非定常流场数值模拟

用有限差分法求解Ｎ－Ｓ方程，数值模拟了细长三角翼作匀速上仰运动时的距声速绕流流场。给出了非定常流场结构，涡破裂特性和非定常升力，阻力特性，研究了上仰运动对动态失速的影响。     

有侧滑角椭圆弹身大攻角状态绕流的数值模拟

用数值模拟方法研究了椭圆细长体在大侧滑角、大攻角状态下的流动特性。数值模拟的出发方程为三维可压缩流的全Ｎ－Ｓ方程，数值格式为ＴＶＤ格式。通过数值模拟，研究了椭圆截面细长体在Ｍ∞＝２．５、侧滑角β＝４５°状态下，攻角从１０°￣３５°范围内流场中的分离特征。     

文章导读

吸气式内外流一体化飞行器动导数数值模拟（147—155，doi：10．7638／kqdlxxb—2015．0020）刘绪，刘伟，周云龙，柴振霞
　　 在标模验证的基础上开展吸气式内外流一体化飞行器 WR—A 动态特性分析，采用 N—S 方程模拟强迫简谐振动的非定常流场，获得了飞行器的直接阻尼导数、加速度导数和旋转导数，并对 WR—A 外形大振幅振动时的进气道性能参数进行了分析。研究表明，在某些情况下反映流动时滞效应的加速度导数占直接阻尼导数的比例最高可以达到40％以上，而在某些情况与直接阻尼导数符号相反。此外，在动导数研究基础上，采用参数化的动导数概念初步开展了非定常气动参数建模研究。     

细长体大攻角绕流的数值模拟

用数值模拟方法研究了细长体大攻角的绕流特性（求解在周向３６０°范围内进行）。数值模拟的出发方程是三维可压流的全Ｎ－Ｓ方程,格式为二阶迎风ＴＶＤ格式。通过数值模拟研究了附近,攻角１９°、３２．５°和５０°状态下,细长体绕流场中的激波、分离流动及表面压力分布等重要特征。     

叶轮机三维非定常流动数值模拟的研究

利用数值模拟手段对压气机内部非定常流场进行了初步研究,在数值模拟中引进了双重时间步方法。研究和讨论了物理时间步长及网格尺寸对计算结果的影响,给出了叶轮机通道中非定常流动的部分流动图画,并对非定常流场进行了初步的分析。     

大迎角振动三角翼跨声速粘性绕流数值分析

本文通过数值求解三维可压缩雷诺平均ＮＳ方程，研究发析了振动三角翼的三维复杂非定常流场，其中对流场项计算格式的空间离散在平行物物理面方向与物面法线方向的空间导数分别采用二阶ＮＮＤ格式和由ＮＮＤ格式高阶插值而得到的三阶迎风格式，时间推进为隐式ＬＵ分解方法，并使用了修正的代数ＢＬ湍流模型模拟高雷诺数的湍流流动。     

基于机器学习的非定常流场网格自适应

现有针对非定常流场数值模拟的网格自适应方法,通常每隔一段时间步就进行一次网格调整,增加了计算复杂度和精度损失。针对这一问题,本文基于间断伽辽金（discontinuity Galerkin, DG）有限元法提出了结合BPNN(backpropagation neural network)和MMPDE(moving mesh partial differential equation)的非定常流场网格自适应方法。该方法首先采用DG有限元法对Navier-Stokes方程进行非定常计算,得到统计意义上的网格间断量;然后以初始网格和间断量训练BPNN回归模型,用于预测任意位置节点的间断量;接着使用MMPDE变分法移动网格节点,使其符合统计意义的间断量分布;最终通过Laplacian网格平滑法保证网格单元质量。圆柱绕流非定常流场算例的验证结果表明,该方法能够在不改变网格拓扑结构和不增加节点数的情况下完成一次性网格自适应,显著提高了非定常流场数值计算的精度和效率。     

耦合POD重构舰面流场的直升机舰面起降数值模拟

为解决舰面非定常流场数据量过大的问题,采用本征正交分解(POD)方法对舰面流场进行重构,发展了一种耦合POD重构流场的直升机舰面起降数值模拟方法。首先采用计算流体力学(CFD)方法计算舰面非定常流场,获得离散数据样本;然后提取流场的POD模态,并截取能够捕捉到原流场主要特征的少量模态对原流场进行重构;最后建立耦合重构舰面流场的直升机高阶飞行动力学模型。以直升机返航进场为例进行数值模拟,并将计算得到的操纵量和飞行状态与飞行试验结果进行对比。结果表明:使用POD方法重构后的舰面流场数据约为原始样本数据的8.5%,且重构流场与原始流场吻合良好,POD方法能够解决舰面非定常流场数据量过大的问题。与飞行试验数据的对比表明,本文方法捕捉到了舰面非定常流场对直升机的影响,可用于直升机舰面起降研究。     

类升力体外形俯仰阻尼特性数值研究

采用有限差分方法求解薄层近似的非定常Navier-Stokes方程，定常流场采用交替方向隐式分解的NND格式，非定常流场采用四步Runge-Kutta方法，在保证时、空二阶精度前提下引入变系数残值光顺技术提高非定常流场计算效率，复杂带翼外形的空间网格通过求解抛物化的椭圆型方程生成，最后在Etkin理论下给出球锥及类升力体外形的俯仰阻尼导数计算结果。     

二维US-FDTD方法的数值稳定和色散分析

研究了二维US-FDTD方法的数值稳定性和数值色散特性。通过对增长因子的计算，证明了US-FDTD方法的无条件稳定性。利用增长因子的相位，推导出了US．FDTD方法的数值色散关系式。分析了US-FDTD方法的数值色散误差。数值分析表明，与ADI-FDTD方法一样，数值色散误差仍然是决定US-FDTD时间步长选取的关键因素。同时发现，数值色散受时间步长及网格大小的影响。     

高阶精度间断Galerkin有限元方法研究

在二维非结构网格上,对高阶精度间断Galerkin有限元方法求解跨音速欧拉方程进行研究。运用间断有限元理论,采用施密特正交化多项式基函数对流场解进行近似描述,使用HLLC近似黎曼解方法计算网格单元边界处的数值通量,积分项通过高斯积分求解,时间推进采用经典四步Runge-Kutta方法,并引入斜率限制器,抑制数值振荡。对NACA0012翼型跨音速无粘流动进行数值模拟,数值结果表明：间断Galerkin有限元方法具有较高的精度,较小的数值耗散和较强的激波捕捉能力。     

轻型燃气轮机燃烧室NOx 排放估算

为了解决NOx排放预估问题,应用数值模拟方法研究轻型燃气轮机燃烧室NOx排放规律。将数值分析结果与文献中经验公式计算数据进行比较与分析,验证数值方法的准确性;基于进口压力、空气质量流量、主燃区温度3个参数拟合NOx排放预估公式。结果表明:数值模拟得到的NOx排放值与文献中经验公式计算出的符合度较高;结合数值计算结果得到的NOx排放预估公式;该公式能有效反映上述3个参数对NOx排放的影响;采用验证后的高精度数值模拟方法建立的NOx预估模型所得结果具有一定的科学性及可靠性,并可降低试验成本。     

矩阵数值秩的一种确定方法及其在控制系统中的应用

一、矩阵的数值秩 矩阵秩的确定在现代控制理论中占有重要的地位。在代数理论中,用行列式或特征值来判断矩阵的秩,由于舍入误差的存在,用计算机判断矩阵秩时,这些方法是不完全适用的。     

跨流域三维复杂绕流问题的气体运动论并行计算

通过研究求解描述跨流域三维绕流问题的Boltzmann模型方程气体运动论耦合迭代数值格式,分析气体运动论数值计算方法的内在并行性;从变量依赖关系、数据通信与并行可扩展性三方面开展基于离散速度空间区域分解计算的研究,发展求解稀薄流到连续流跨流域三维复杂绕流问题的并行算法。通过对不同Knudsen数、不同马赫数、不同攻角三维球体及返回舱绕流的并行计算,计算结果与实验数据和理论分析吻合较好。研究表明,该并行算法负载平衡和并行可扩展性较好,对不同并行计算机系统具有很好适应性,显示直接求解分子速度分布函数的气体运动论计算方法有良好的并行计算特性。     

Numerical investigations on the VFE-2 65-degree rounded leading edge delta wing using the unstructured DLR TAU-Code

The numerical simulation of the flow around a 65° delta wing configuration with rounded leading edges is presented. For the numerical simulation the DLR TAU-Code is used which is based on an unstructured hybrid mesh approach. Within this paper several numerical results are shown, solving the steady RANS equations by different turbulence models. The simulations are carried out within the RTO/AVT-113 task group focusing on experimental and numerical research on delta wing configurations with rounded leading edges. Within this paper the focus is related to the flow topology depending on the angle of attack as well as on Reynolds number effects. Finally the results are compared and verified by experimental data.     

Hoff型夹层板振动分析的边界积分方程法

 本文利用静力问题的基本解建立了Hoff型夹层板振动分析的边界积分方程。在数值实现过程中,边界上和区域内的未知函数均釆用线性插值以提高计算精度。数值结果表明本文方法具有未知量少和计算精度高等优点。对于低阶频率,本文结果与精确解的误差不大于1％。     

消失球方法及其对旋翼厚度噪声的分析研究

对旋翼厚度噪声及其数值模拟方法进行了讨论,着重解决了以下两个问题：（1）通过理论和数值分析,讨论了使用消失球公式计算旋翼旋转噪声时计算结果的有效性;（2）通过计算机模拟,深入分析了旋翼厚度噪声的特征。     

飞行品质仿真计算与试飞数据对比分析

随着飞行控制系统越来越复杂,试飞过程逐渐向数字化过渡。在开展飞行试验前,通过数值模拟对飞机的主要飞行品质参数进行仿真计算,可以为试飞过程提供参考,提高飞行试验的效率。根据某型飞机的风洞试验数据,建立该型飞机的气动力数学模型,基于Flightgear飞行模拟软件搭建功能完备的实时飞行模拟平台,完成飞行品质的数字化试飞。选取典型试飞数据对所建立的气动力模型与数值仿真结果进行验证,结果表明:数值仿真结果与试飞结果具有较好的一致性。     

跨音速喷管化学反应流数值模拟

用时间相关的半隐格式有限差分数值方法求解了化学非平衡反应跨音速喷管流场，在喷管收敛段，流动接近化学平衡状态，控制方程的刚性问题严重，数值积分困难。通过对时间差分项隐式离散、对空间差分项显式离散，流场边界采用参考平面上的特征线法计算等，成功地解决了由于化学反应有限速率带来的数值解不稳定问题。该格式简单、需要计算机存贮空间少。本文完成了一维和轴对称非平衡化学反应喷管流动计算，并与化学平衡流和冻结流的计算结果做了比较。