基于动态非结构重叠网格法的直升机前飞非定常流场数值模拟研究

 发展了非结构重叠网格技术，提出适合于直升机旋翼运动的动态非结构重叠网格算法，该算法使得在采用非结构重叠网格模拟直升机流场时实现程序自动处理网格，不需任何人工干预。流动控制方程采用非定常可压缩Euler方程，空间离散采用Jameson有限体积法，时间推进采用双时间步长法，为加速收敛，在双时间步长内循环迭代中采用预处理技术。为验证本发展的预处理解算器以及动态非结构重叠网格的算法的可信性和有效性，开展一系列的数值实验，模拟了不同速度范围的NACA0012翼型绕流,Caradonna旋翼悬停和前飞流场以及Robin直升机机身-旋翼组合的前飞非定常流场，结果表明采用本文算法得到的结果是可信的，网格处理是自动高效的，可以应用于直升机全机流场模拟，为直升机设计提供可信指导数据。     

悬停旋翼无粘流场数值模拟中的多重网格方法

发展了一种基于不同空间离散格式的多重网格算法,并应用于悬停旋翼无粘绕流的Euler方程数值模拟。由于悬停旋翼流场中存在不可压区域,同时旋翼尾涡系统的发展需要较长的时间,使得旋翼流场的收敛速度远低于固定翼流场,因此研究旋翼流场的多重网格算法具有重要意义。空间离散采用了Roe s FDS格式和Jameson中心有限体积格式,时间推进应用了五步Runge-Kutta方法。采用多重网格V循环方式,对一跨声速悬停旋翼无粘流场进行了数值计算。计算结果表明:多重网格算法可以显著加速悬停旋翼无粘流场的数值计算收敛速度;无论在激波分辨率还是在计算精度上,Roe s FDS格式都优于JST格式。     

基于涡限制法的旋翼气动噪声计算

在直升机旋翼流动的数值模拟中,流动求解精度对旋翼气动噪声的预测有着重要影响.采用中心格式有限体积法求解可压缩Euler方程数值模拟旋翼绕流,并采用基于Ffowcs Williams- Hawkings( FW-H)方程的声类比方法计算旋翼气动噪声.为了抑制中心格式的数值耗散,保持机翼翼尖以及旋翼桨尖的尾涡结构,引入了涡...     

用双时间法数值模拟悬停旋翼流场(英文)

本文用改进的双时间法和重叠网格通过求解可压缩的欧拉和Ｎ－Ｓ方程计算了悬停旋翼流场。在双时间法中,我们采用牛顿松弛以加速每个真实时间步中的子迭代。高效的双时间法加上重叠网格使在微机上计算复杂的悬停旋翼流场成为可能,计算结果和实验符合得很好,计算得到的下洗流和载荷证明了我们的方法是可行的。     

基于运动嵌套网格的前飞旋翼绕流N-S方程数值计算

 通过求解 Navier-Stokes方程数值模拟了直升机旋翼前飞非定常流场。为了模拟包括旋转、周期性变距和周期性挥舞的非定常运动,采用了一种能够快速完成重叠网格间流场信息交换的运动嵌套网格方法。空间上采用中心平均的有限体积法进行离散,时间方向应用含隐式子迭代的双时间法推进求解。为了验证程序的正确性,数值计算了一有升力悬停流场,旋翼桨叶表面压力分布的计算值与实验值吻合很好。     

基于隐式算法的悬停旋翼粘性绕流高效CFD分析方法

为显著提高旋翼粘性绕流CFD模拟效率,建立了一套高效的基于隐式LU-SGS算法和OpenMP并行策略的旋翼悬停流场求解方法。首先,求解Poisson方程生成桨叶剖面翼型的贴体正交网格,并通过剖面网格插值、翻折方法生成桨叶C-O型贴体网格;在此基础上,采用基于"扰动衍射"挖洞方法与"Inverse map"相结合的洞边界划定与贡献单元搜索方法,解决了嵌套网格技术中的相关瓶颈问题。然后,以非惯性系下耦合S-A湍流模型的RANS方程为流场主控方程,对流通量采用三阶Roe-MUSCL格式进行离散,时间推进采用高效的隐式LU-SGS方法,同时采用基于数据共享的OpenMP并行策略加速流场求解。最后,运用所建立的方法分别对不同旋翼翼型和悬停状态"Caradonna-Tung"以及UH-60A旋翼的流场及气动特性进行了计算,并给出了根据涡核位置加密网格来提高桨尖涡捕捉精度的方法,同时将计算结果与试验值进行了对比,验证了该方法在旋翼CFD流场模拟中的高效性和高精度特征。     

先进直升机旋翼悬停状态气动性能计算

 为更好地模拟先进直升机旋翼流场,同时准确计算其悬停状态的气动性能,建立了一套基于Navier-Stokes/Euler方程的混合CFD方法。该方法的求解域由两部分组成：一是围绕旋翼桨叶周围的黏性区域,采用可压缩Navier-Stokes方程来模拟旋翼附近的黏性流动和近场尾涡的捕捉;二是离桨叶较远、黏性可以忽略的远场区域,用Euler方程来描述其流动。在该方法中,将三阶逆风格式（MUSCL）与通量差分分裂方法相结合,无需添加人工黏性,因而可有效地减少旋翼尾迹数值耗散。?阌诹鞒》智蠼庖约爸芷谛员呓缣跫氖凳?采用了嵌套网格方法,并给出旋翼网格与背景网格交界面的信息传递方式。应用所建立的计算方法,首先对二维翼型、三维M6机翼的流场进行了数值模拟,以验证计算方法;然后,着重计算了有实验结果可供对比的具有先进气动外形的HELISHAPE 7A模型旋翼和UH-60A直升机旋翼,通过计算得到旋翼表面压力分布、桨叶展向拉力系数分布、桨叶表面细节流动以及气动性能等,进一步验证了该方法的有效性。     

用运动嵌套网格方法数值模拟旋翼前飞非定常流场

通过求解Ｅｕｌｅｒ方程数值模拟了直升机旋翼前飞非定常流场。为了模拟包括旋转、周期性变距和周期性挥舞的非定常运动,本文采用了一种能够快速完成重叠网格间流场信息交换的运动嵌套网格方法。空间采用中心平均的有限体积法进行离散,时间方向采用双时间法隐式推进求解。为了验证方法的正确性,数值模拟了悬停流场,数值结果和实验值吻合很好。尽管前飞状态的计算没有实验验证,但计算结果与其它文献计算结果吻合很好。     

用双时间法数据模拟悬停旋翼流场

酿和改进的双时间法和重叠网格通过求解可压缩的欧拉和Ｎ－Ｓ方程计算了悬停旋翼流场。     

基于重叠网格技术和多重网格算法的悬停旋翼粘性绕流数值模拟

发展了一种基于重叠网格的多重网格算法,并应用于悬停旋翼粘性绕流的Navier-Stokes方程数值模拟.多重网格算法在重叠网格上实施主要面临如下两个问题:粗网格间几何关系的建立,即洞单元、边界单元的确定;洞边界存在的情况下如何进行粗细网格间流场信息的传递.本文针对以上问题给出了相应的解决方案.采用两层网格的V循环,应用有限体积法对一跨声速悬停旋翼粘性流场进行了数值计算.计算结果表明该多重网格方法可以显著加快悬停旋翼粘性流场数值解收敛速度,加速比约为3.     

基于拟压缩方法的叶片扩压器低速流场数值研究

摘要 发展了基于拟压缩方法的求解叶轮机械低马赫数流场的方法,研制了了时间推进方法对叶轮机械全流速范围流场的求解的计算程序,且对拟压缩因子的取值进行了讨论。通过对离心压缩机叶片扩压器的数值计算表明,本文发展的基于拟压缩方法能能够快速正确的预测叶轮机械内的低速流场;引入拟压缩方法是对时间推进方法的一个很好补充,使时间推进方法的应用范围得到了拓展。     

拟压缩性方法数值模拟绕弹体的不可压缩层流脱体涡流动

本文应用在连续生方程中增加压力对时间导数项的拟压缩性方法，数值求解三维定常不可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程。数值求解过程采用隐式的Ｂｅａｍ－Ｗａｒｍｉｎｇ近似因式分解格式。计算绕尖头旋成柱体的大迎角脱体涡流动，得到与实验相符的结果。     

基于时间谱方法的振荡翼型和机翼非定常黏性绕流数值模拟

 传统的双时间方法在非定常计算中长时间的过渡迭代推进求解导致其计算效率相对较低,针对周期性非定常流动问题的流动特征,发展了一种基于离散傅里叶变换的高效时间谱方法,用于求解振荡翼型和机翼的非定常黏性绕流。在时空耦合的雷诺平均Navier-Stokes(RANS)方程的求解中,对流项的离散应用了Roe的通量差分格式,物理时间项的离散方法为时间谱方法,伪时间推进采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式。考虑到湍流的时空耦合效应,时空耦合的Spalart-Allmaras一方程湍流模型的物理时间项同样采用时间谱方法进行离散。为了进一步提高计算效率,当地时间步长和多重网格技术等加速收敛的措施均被采用。算例对俯仰振荡NACA0012翼型和Lann机翼的周期性非定常流场进行了数值计算。结果表明:对于周期性非定常流场的数值模拟,相比于传统的双时间方法,用时间谱方法近似物理时间项,不仅能够提高流场的计算精度,而且更能够大幅度提高计算效率。     

三维方柱不可压缩绕流的大涡模拟计算

用非定常不可压缩流动求解方法对三维方柱绕流问题作了大涡模拟计算研究.求解程序采用虚拟压缩方法和双时间步的时间推进方案,对流项的差分采用高阶精度的数值差分格式.通过与实验和文献结果比较分析,验证了基于虚拟压缩方法的湍流大涡模拟计算方案对于复杂的三维钝体绕流计算是可行的.     

一种新型颤振激励系统的特性研究

在非结构网格的基础上，以欧拉方程为流场的主控方程，采用格心格式空间离散的有限体积法和四步Runge—Kutta显式时间推进的方法．对一种新型颤振激励系统的流场进行了数值研究。研究了该系统的工作原理和激励力特性。研究结果表明，该系统能很好地激励飞机结构。     

圆柱单孔高频压力探针测量高亚音速二维周期流场的方法

本文详细介绍了作者发展的圆柱单孔高频压力探针测量高亚音速二维周期性流场的方法。该方法简化了压缩性影响的修正,测量时置探针于九个方位角。它已被成功地应用于单、多级高速压气机转子出口流场的测量      

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

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

基于拟压缩方法的扩压器低速流场数值研究

发展了基于拟压缩方法的求解叶轮机械低马赫数流场的方法,研制了时间推进方法对叶轮机械全流速范围流场的求解计算程序,且对拟压缩因子的取值进行了讨论.通过对离心压缩机叶片扩压器的数值计算表明,发展的基于拟压缩方法能够快速正确的预测叶轮机械内的低速流场;引入拟压缩方法是对时间推进方法的一个很好补充,使时间推进方法的应用范围得到了拓展.      

后缘喷流对三角翼绕流影响的N-S方程数值分析

本文用拟压缩性方法求解不可压流雷诺平均拟压缩Ｎ－Ｓ方程组,对带有后缘喷流的三角翼粘性绕流进行了数值模拟,求解中采用了Ｂｅａｍ－Ｗａｒｍｉｎｇ隐式近似因子分解格式以及ＭＭＬ代数湍流模型。计算结果说明,后缘喷流使涡核压强降低,使涡核速度增大,从而对三角翼前缘分离涡有稳定作用,并能增大上翼面的负压值和下翼面的正压值,从而可以增加部分升力。计算结果还说明,喷口面积或喷流下偏会使上述作用增强。     

跨音速颤振计算方法研究

采用快捷的动态弹性变形技术生成三维块结构动态贴体网格。结构运动方程采用二级精度的龙格—库塔时间推进。气动力求解则用非定常欧拉方程的双时间有限体积推进,外时间为物理时间,与结构运动方程同步;对每个真实物理时间步采用五步Runge Kutta时间推进进行虚拟定常迭代。针对可压流颤振计算中存在的质量不相似问题,采用变质量、变刚度的方法计算出质量匹配点处的颤振速压值,并根据颤振速压随质量或刚度倍数的变化趋势,得出可用的飞机跨音速颤振速压。对一飞机简化外形的跨音速颤振特性计算得到了合理的计算结果。