基于声扰动方程的气动噪声传播积分计算方法

传统声比拟思想假设声传播介质均匀,所获得的远场噪声积分方法无法考虑近场流动非均匀性引起的声折射效应,声场预测结果往往存在较大的数值误差。为解决上述问题,本文采用声扰动方程描述声波运动,并利用声学互易原理,在频域内构建关联近场声源与远场声压的伴随格林函数,从而将非均匀介质折射效应纳入伴随格林函数。基于伴随格林函数和声扰动方程源项的数值结果,利用积分方程计算气动噪声的空间传播。结合上述方法对二维圆柱和NACA0012翼型绕流噪声开展了数值预测研究,远场声压预测结果与CFD直接计算值吻合。数值研究表明,本文提出的基于声扰动方程的气动噪声传播积分计算方法可以考虑非均匀介质对声传播的影响,能够提高气动噪声的预测精度。     

基于混合RANS/LES方法与FW-H方程的气动声学计算研究

在气动噪声的数值计算过程中,非定常流动的求解精度对声学计算结果有着重要的影响.以机体噪声计算的标准算例双圆柱绕流气动噪声问题为研究对象,采用基于非线性k-ε湍流模型的限制数值尺度(LNS)方法对双圆柱绕流进行了数值模拟,将计算得到的气动特性和流动特征与相应的试验结果进行了对比分析.为了求解远场观测点处的气动噪声,在精确求解双圆柱绕流流动的基础上结合基于FW-H(Ffowcs Williams-Hawkings)方程的声类比方法进行数值计算,并通过圆柱体的展向相关性对计算结果进行了修正,将得到的最终结果与相应的声学试验结果进行了对比,两者吻合良好,表明该数值方法是准确、可靠的.     

围壳形状对潜艇流致噪声的影响计算

针对潜艇在航行中产生的流噪声现象,建立了适合于潜艇绕流过程中流噪声计算的三维数值方法和计算程序。其基本思路为:采用多块技术与SIMPLE方法求解RANS方程计算定常流场,在稳态流场的基础上,采用迭代时间推进法和DES湍流模型求解非定常流场,再根据瞬态流场结果采用Lighthill声比拟方法获得声观测点随时间变化的声压。对尺度放大的潜艇标模SUBOFF以及三种改进的围壳外形进行了流致噪声比较计算,获得了不同围壳外形的声压频谱和总声压特性以及典型模型的表面脉动压力分布规律。结果表明:把潜艇围壳后缘形状从矩形变为梯形,可控制围壳的尾流摆动,能够有效降低潜艇流噪声,最大的降噪效果可达5dB左右。     

非紧致气动噪声半空间传播基本解的边界积分方法

为了获得低马赫数流动诱发的非紧致气动噪声在半空间内传播的基本解,结合复等效源方法和边界元方法建立了半空间精确格林函数的边界积分方程,当半空间边界为刚度型阻抗边界时可避免奇异性积分。同时基于等效源方法提出一种半空间二维非紧致圆柱声散射模型,推导了静止介质中声散射基本解的理论表达式。对静止介质中的二维圆柱声散射,数值解在研究的频率与观察点处与理论解一致。采用数值方法计算二维圆柱绕流诱发的半空间声场基本解,结果显示半空间边界强化了声源的声辐射能力,小于马赫数0.2的介质运动对声传播的影响可以忽略。     

非结构网格高超声速化学非平衡MHD流场数值模拟

 在非结构网格上对考虑化学非平衡效应的二维高超声速磁流体绕钝头体流动进行了数值模拟。控制方程由二维理想磁流体动力学（MHD）方程和组元连续方程两部分组成,化学动力学模型为5组元17反应模型。MHD方程空间离散采用AUSM格式,时间推进采用显式5步龙格-库塔格式,并通过弱耦合的方式与化学反应控制方程结合在一起。计算模型为二维钝头体,外加磁场为偶极子场,磁场源位于钝头体内部。在高超声速来流条件下,对有、无磁场干扰,是否考虑化学反应下的4种工况进行了数值计算,得到了满意的结果,并与有限的参考文献进行了对比。结果表明本文发展的方法能准确地模拟考虑化学非平衡效应的高超声速MHD流场。     

基于比拟理论计算圆柱/翼型的气动噪声

采用计算流体力学(CFD)与"声比拟"相结合的方法,通过两个步骤模拟流动的声学远场。第一步,在包含所有声源的近场区域内,通过求解URANS方程获得控制面处的非定常流场参数;第二步,采用基于可穿透数据面的FW-H方法模拟声学远场。该方法与传统的半经验方法相比,计算精确,易于工程实现,并且可以计算非线性噪声。根据气动声学计算积分公式,构造了高阶的时间离散格式,采用高斯积分公式计算圆柱/翼型的气动噪声。在时域和频域上计算了观测点处的声压及声压级随时间步数的变化情况,其结果与国外实验结果对比取得了较好的一致性,此外研究了翼型尾流部分界面对计算精度的影响。     

NND格式在航天飞机头部段N—S方程求解中的应用

本文采用实质上为二阶精度的无波动、无自由参数的耗散差分格式(NND格式)求解了非定常的二维和三维完全N-S方程,对类似于航天飞机头部外形的超声速粘性气体绕流流场进行了模拟。在计算中,使用了代数网格生成技术,对脱体激波采用了装配法。作为验证算例的二维圆柱绕流计算结果同文献[4]的数据进行了比较。二维凹陷外形绕流的计算结果同无粘流结果进行了对比,由于粘性流动出现了分离,二者的壁面压力是有差异的。最后,对三维外形的绕流也进行了模拟,计算结果提供了流场的细节。     

基于多重网格方法的跨声速颤振数值模拟研究

使用多重网格的方法对二维翼型和三维机翼的跨声速颤振进行了数值求解.流场控制方程为N-S方程,湍流模型采用SA模型网格.计算网格采用结构网格,空间离散采用中心格式,使用双时间法进行时间推进.流场与结构之间的数据通过径向基函数(RBF)插值方法进行交换.通过耦合求解流场控制方程和结构运动方程得到Isogai(case A)机翼模型和Agard445.6机翼的跨声速颤振边界.通过与使用单重网格时得到的结果相比较得出,多重网格方法能够节省计算时间,提高颤振分析计算效率.     

自动编码器在流场降阶中的应用

自动编码器作为一种压缩算法,在数据降维和去噪等方面有着广泛实践,有条件作为一种降阶方法在流场识别与数据处理方面得到应用。文章中以圆柱绕流为例,首先对圆柱后速度场建立了编码模型,用来对原始数据进行降维和特征提取,之后将编码后的数据与流场特征量相关联,建立了由流场编码回归圆柱表面压力系数的神经网络,探索了降维后数据的应用。结果表明,自动编码得到的结果能够承载原始速度场的主要信息,解码后速度场与原速度场测试均方根误差小于0.02,压力回归测试均方根误差可小于0.1。说明自动编码器能够作为一种流场的特征提取和降阶方法,在未来得到更广泛的应用。     

一种用于飞机结冰的海绵冰模型及应用

基于两相流理论,采用欧拉法对飞机结冰进行数值模拟,在空气流场速度分布的基础上,求解水滴流场,得到水滴收集率和撞击极限.结果表明:欧拉法比拉格朗日法在计算水滴收集效率和撞击极限方面有优势,尤其是在较大的水滴平均等效直径情况下.对翼型的冰形和圆柱表面的温度预测数据和实验数据进行对比,圆柱冰形厚度误差在15.0%之内,计算结果与实验结果吻合较好,同时表明结冰模型是有效可行的.      

基于非结构网格CFD技术的旋翼气动噪声计算方法研究

将基于非结构网格技术的旋翼流场CFD计算方法与基于FW-H和Kirchhoff方程的声学方法相结合,建立了一套既适合于直升机旋翼厚度、载荷和桨-涡干扰噪声,又适合于跨声速高速脉冲噪声的综合计算模型。为提高旋翼流场及桨叶表面气动载荷计算的精度,主控方程的求解采用了三维可压非定常的N-S方程,网格划分则使用非结构运动嵌套网格方法。在噪声计算中,通过FW-H方法计算旋翼的厚度噪声和载荷噪声,并选取能够包含流场非线性区域的旋转面作为Kirchhoff积分面,由Kirchhoff方法计算包含四极子项的高速脉冲噪声。应用该模型,以AH-1旋翼为算例,计算了不同飞行状态下的旋翼气动噪声,并与可得到的试验结果进行比较,验证了方法的有效性。然后,着重对两种声学方法对计算结果的影响进行了对比研究,并分析了旋翼厚度噪声、载荷噪声和四极子噪声的特性。     

多重网格加速的LUSGS算法用于喷管流场数值模拟

采用有限体积法,二阶精度插值的Roe格式对三维流场N-S方程进行空间离散;采用LUSGS隐式时间推进法和多重网格技术加快计算收敛速度;采用独特的显式算法克服了两方程湍流模型刚性对流场计算收敛速度的影响,并给出了该方法的理论依据。通过轴对称喷管内流场、二元喷管内外流场的数值模拟验证了采用以上方法在计算效率上的提高。最后将RNG模型计算结果与SA模型计算结果及实验数据做对比分析。     

五种湍流涡粘模型在二维方柱绕流数值模拟中的对比研究

为研究不同雷诺时均Navier-Stokes(RANS)模型在求解钝体绕流场的差异,采用五种两方程湍流涡粘模型计算了二维方柱的绕流场(雷诺数Re=22 000),得到了不同湍流模型在不同网格离散方案下的方柱气动力与流场特性。结合以往文献数据,通过对比不同湍流模型计算结果的差异,揭示了不同湍流模型的特点。研究结果表明:网格离散方案对方柱绕流数值模拟结果有重要影响;RNGk-ε湍流模型的计算效率最高,SSTk-ω湍流模型的计算效率最低;Standard k-ε湍流模型的计算结果准确程度整体弱于其余湍流模型;RNGk-ε湍流模型、Realizable k-ε湍流模型与Standard k-ω湍流模型的计算结果大致相当,较接近大涡模拟结果;SSTk-ω湍流模型的模拟结果优于其余湍流模型,其尾流区速度场与试验结果吻合较好;两方程湍流模型计算的二维方柱绕流结果与试验结果以及大涡模拟结果相比,存在不可忽略的差异。     

带下反桨尖旋翼气动噪声数值分析

采用高精度、高效的CFD(计算流体力学)方法求解旋翼流场,在获得准确的流场信息的基础上,噪声计算采用基于声类比法的FW-H方程进行求解。利用算例对气动噪声预估方法进行了验证,针对带不同下反角度桨尖旋翼的悬停状态和前飞状态的气动噪声进行了对比计算分析,着重开展了下反角为0°、20°和45°的三副旋翼的近场噪声及地面噪声计算分析,结果表明,下反桨尖具有降噪效果,能够在一定程度上抑制旋翼气动噪声。     

浸入边界法模拟小圆柱对主圆柱涡脱落的抑制

使用浸入边界法研究了小圆柱对主圆柱涡脱落的抑制.方法中使用非贴体笛儿尔网格,易于处理包含复杂边界的流动问题.采用离散附加力直接加入边界条件方法对虚单元进行重构,使边界条件在浸入边界上精确满足.使用隐式分步法解二维非定常不可压Navier-Stokes方程,通过速度和压力解耦提高计算效率.数值模拟单圆柱绕流及不同位置小圆柱和主圆柱的流动干扰,通过分析流场涡结构和升、阻力系数,得到小圆柱对主圆柱涡脱落的延迟和抑制作用.计算结果与已有实验结论和数值结果对比,计算误差不超过5%,说明浸入边界法可以简单有效地处理圆柱涡脱落抑制这类流动干扰问题.     

烧蚀环境下的圆柱绕流计算模型

基于多学科理论建立了高温高压气流环境下圆柱烧蚀、剥蚀的数理模型。利用离散涡方法计算流场与圆柱表面压力分布;采用三方程烧蚀模型计算热化学控制机制下的烧蚀速率,并将其与扩散控制机制下的结果相比较以确定烧蚀量;引入颗粒轨道模型求解剥蚀颗粒的运动。基于该模型对高温高压燃烧室内圆柱形烧蚀试件的绕流场、烧蚀量以及剥蚀颗粒的运动轨迹进行了编程计算分析。研究表明:低雷诺数条件下烧蚀环境对绕流流场的影响不大,而较高雷诺数条件下烧蚀绕流流场与无烧蚀绕流流场的涡街分布存在显著差别;在低雷诺数条件下流场对烧蚀速率的影响甚微,较高雷诺数条件下烧蚀速率略大于无流场烧蚀情况;涡量的存在决定着剥蚀颗粒的运动轨迹与分布,进而影响尾流热场与烧蚀速率。该模型可为热防护材料的烧蚀实验提供参考,并为烧蚀、剥蚀过程的非线性分析提供数据支持。     

跨声速条件下轴对称收扩喷管内外流场的数值研究

采用有限体积法,二阶精度插值的Roe格式对三维流场N-S方程进行空间离散;采用LUSGS隐式时间推进法和多重网格技术加快计算收敛速度;采用不同湍流模型、壁面函数对轴对称喷管内外流场进行三维数值模拟.计算结果与实验数据对比表明,采用涡粘修正的RNG模型,内流场使用增强壁面函数,外流场使用标准壁面函数可以得到准确的计算结果.研究表明:在亚声速条件下,喷管外壁面压力分布随落压比的变化不大;超声速条件下,喷管外壁面激波随着落压比的增大前移.      

三维副翼铰链力矩计算

 数值模拟了带副翼偏转构型的二维、三维亚、跨、超声速来流情况的绕流流场,采用有限体积方法和AUSM+迎风格式数值求解N-S方程,四步龙格-库塔时间推进,湍流模型为Baldwin-Lomax和Spalart-Allmaras湍流模型。二维模拟压强分布与实验数据以及CFD软件模拟结果吻合较好。采用交错对接网格,数值模拟了不同副翼偏角和缝隙、来流马赫数、迎角和雷诺数的三维构型流场。通过舵面铰链力矩与实验数据的对比分析表明了该数值方法模拟此类复杂流场的可靠性和舵面铰链力矩计算的有效性。发现并研究了亚声速来流时铰链力矩随迎角的反向变化趋势,初步分析了副翼的缝隙和雷诺数效应。     

面向稀薄流非线性本构预测的机器学习方法

稀薄非平衡流域内连续介质假设已经失效,主要围绕Boltzmann方程及模型方程对稀薄非平衡流开展理论与计算研究,统一气体动理论格式（UGKS）是其中一种代表性方法。在稀薄非平衡流数值模拟中,Navier-Stokes （N-S）方程连续介质假设已经失效,不能有效描述流场非平衡特征。UGKS方法虽然计算精度高,但速度空间离散导致计算效率低下,多维高速条件下数值计算难以开展。基于数据驱动的思想,在N-S方程与UGKS方法的研究基础上发展出了一种稀薄非平衡流非线性本构关系求解方法（DNCR）。该方法以N-S与UGKS求解器获得的流场数值模拟计算结果作为训练数据集,基于流场特征参数采用极端随机树算法生成机器学习模型,对预测流场中线性黏性应力项与热流项进行非线性修正,并耦合非线性本构关系求解宏观守恒方程得到目标状态稀薄非平衡流动数值解。针对DNCR方法中所采用的机器学习方法-极端随机树模型,通过二维顶盖驱动方腔流算例对高维非线性建模涉及的特征参数选取、参数调优开展了相关验证工作,选取若干典型状态对极端随机树模型的泛化性能开展研究,并评估了相关模型与方法的计算精度与计算效率。     

多级轴流压气机特性的一种计算方法

本描述了多级轴流压气机特性计算的一种方法，流场计算采用流线曲率法解完全径向平衡方程。本在特性计算中引进了环壁附面损失模型，考虑了叶片去亏损力，二次流和径向间隙漏流等影响。采用了新的流场计算模型，流场计算稳定，收敛较快。并对计算结果与试验数据进行了对比，认为这种计算方法可以用于轴流压气机特性计算。