非结构聚合多重网格法流场数值模拟研究

 在非结构网格上应用聚合多重网格技术来加速 Euler方程的收敛过程。设计了一种高效率的网格聚合方法。采用四重网格,在每一层网格上采用有限体积法进行计算。通过对多段翼型绕流的数值求解验证了多重网格加速收敛的高效性。     

高阶精度方法下的湍流生成项对低速流动数值模拟的影响研究

基于雷诺平均的Navier-Stokes方程和结构网格技术,采用五阶空间离散精度的WCNS格式,开展了SST两方程模型不同湍流生成项对低速流动数值模拟的计算分析。主要目的是为高阶精度格式在复杂外形上的应用提供技术支撑。计算模型包含了低速NLR7301两段翼型和Trap Wing高升力构型,研究内容主要包括不同湍流生成项对收敛历程、边界层湍流粘性系数分布、边界层速度分布、压力系数分布、气动特性的影响。在与试验数据对比的基础上,计算结果表明:对于低速二维流动,不同湍流生成项对收敛历程有比较明显的影响,对附面层湍流粘性系数分布和速度型影响不明显,不同湍流生成项主要影响主翼前缘的吸力峰值,进而影响升力系数和压差阻力系统;对于低速三维流动,不同湍流生成项对低速流动的收敛特性影响不明显,对翼梢涡的模拟精度有比较明显的影响,进而影响翼梢站位的压力分布和总体气动特性。     

径向基函数在三维Euler方程数值计算中的应用

基于三维Euler方程,对非结构四面体网格给出了一种基于紧支径向基函数重构的ENO型有限体积格式,方法的主要思想是先对每一个四面体单元构造插值径向基函数,而在计算交界面的流通量采用高斯积分公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法。最后用该格式对一些典型算例进行了数值模拟,结果表明该方法计算速度快,对间断有很好的分辨能力。     

非结构动网格生成的弹簧-插值联合方法

 提出了基于弹簧法的两点改进,以解决边界发生大变形时的非结构网格变形问题。为了使边界运动引起的网格变形能更好地由边界传递到内部网格中,提出了一种新的基于Delaunay网格插值的弹簧倔强系数逐层改进方法。为了提高弹簧法的计算效率,引入背景网格和直接插值方法,提出了弹簧-插值法。弹簧-插值法首先生成计算域的背景网格(粗网格),然后由弹簧法求解边界运动引起的背景网格变形,最后利用变形后的背景网格直接插值得到计算网格的变形。算例结果表明：改进后的方法一方面有效地提高了动网格的变形能力和变形后的网格质量;另一方面通过降低弹簧法的求解规模,显著地提高了动网格的变形效率。     

高阶精度方法下的湍流生成项对跨声速流动数值模拟的影响研究

利用五阶空间离散精度的WCNS格式和多块结构网格技术,通过求解雷诺平均NS方程,开展了SST两方程模型不同湍流生成项组合方式对跨声速流动数值模拟影响的计算分析。研究的主要目的是为高阶精度格式在复杂外形上的工程应用提供技术支撑。计算模型采用了RAE2822超临界翼型和DLR-F6翼身组合体构型。研究内容主要包括不同湍流生成项对残差收敛历程、边界层湍流粘性系数分布、边界层速度分布、压力系数分布以及模型整体气动力特性的影响。不同湍流生成项组合方式的流场计算结果还与风洞试验数据进行了对比。研究结果表明:对于小迎角不存在明显分离的跨声速流动,不同湍流生成项对流场的高精度计算结果的影响很小,可以不用考虑。     

基于SST湍流模型的模拟SRM内流场数值仿真

为准确模拟固体火箭发动机燃烧室内流场,采用基于格心的迎风型有限体积法求解定常雷诺平均Navier-Stokes方程,在空间离散方法上,采用AUSM-PW矢通量分裂格式,时间推进采用三阶三步TVD型Runge-Kutta显式方法,将Ment-er F R提出的SST( shear-stress-transport)湍流模型及其改进形式用于燃烧室湍流流场的数值模拟,并将计算结果与Wilcox的和Spalart-Allmaras湍流模型进行了对比。结果表明, Menter F R的SST湍流模型计算的燃烧室内的径向速度分布与实验值吻合得最好,最大误差约为5.1%;计算的燃烧室内湍流强度分布与实验的规律一致,而其余湍流模型计算的结果与实验值有很大差异。     

基于Transition SST模型的高雷诺数圆柱绕流数值研究

为了研究高雷诺数下圆柱绕流边界层的转捩现象和圆柱尾迹近壁区的流动特征,首先通过在典型雷诺数下采用Transition SST四方程转捩模型模拟圆柱绕流得到的结果与实验结果及采用SST k-ω两方程湍流模型模拟结果的对比分析,验证了Transition SST模型在模拟高雷诺数下圆柱绕流的优越性,并较为准确地预测出了圆柱绕流边界层的转捩现象及尾迹近壁区的流动特征。然后分别对亚临界区、临界区、超临界区和过临界区的圆柱绕流问题进行了数值模拟,分析了不同雷诺数下圆柱绕流的流场结构及圆柱表面压力系数、摩擦力系数的变化规律,研究了圆柱绕流近壁区的流动特征、边界层转捩的流动机理、转捩位置及其随雷诺数的变化规律。结果表明,亚临界区二维圆柱绕流边界层发生层流分离,无分离泡和转捩现象;临界区和超临界区二维圆柱绕流边界层先产生了分离泡现象,之后流动发生了转捩并在转捩后发生湍流分离;过临界区二维圆柱绕流边界层流动在转捩之后发生湍流分离,无分离泡现象;在临界区、超临界区和过临界区,二维圆柱绕流边界层转捩位置随雷诺数增大向前驻点移动。     

基于Euler方程解摄动的超声速旋转导数计算方法

基于笛卡尔网格和采用Euler方程基本流场解,通过Riemann不变量摄动,对旋转扰动引起的物面法向速度构建摄动关系,发展了一种超声速旋转导数快速计算方法。文中给出了具体的推导过程,可以看出由此得到的超声速旋转导数,仅与定常欧拉解基本流相关。因此,一旦CFD计算得到定常欧拉解,就可一次性地快速计算出感兴趣的超声速旋转导数。为了验证该计算方法,本文选用国际动导数标准模型Basic Finner算例进行了计算,计算结果与试验数据、文献数据具有很好的一致性,展示出本文方法计算耗时少、精度高、工程实用性强等优点,且该算法基于笛卡尔网格,特别适合复杂气动外形飞行器超声速旋转导数的快速计算。     

浸入式贴体网格边界方法

提出一种基于重叠网格的浸入式边界方法,用以模拟复杂的三维黏性跨声速流动.该方法采用一套固定的笛卡儿正交网格,用以主流的求解;采用一套可以移动的贴体等距面网格,用以拟合或者离散物面的作用力,通过空间插值,实现两套网格重叠部分的信息传递.分析了浸入式贴体网格边界方法的优势,介绍了贡献单元的寻找策略和物理通量的插值方法.流场...     

N-S方程数值模拟方法研究及边界层中相干结构的非线性演化

建立了一套三维不可压流动的直接数值模拟方法。该算法在x及y向构造了基于非等间距网格的紧致有限差分格式和非线性项的迎风紧致型格式,在z向采用Fourier谱方法。并研究、提出了三维、非定常流体运动下游边界的无反射出流条件。另外,本文还构建了多涡结构的初始扰动理论模型,利用直接数值模拟的方法研究了边界层中多涡结构的非线性演化特点。结果表明多涡结构非线性演化中的许多特点,如雷诺应力分布、高剪切层的形成、喷射和扫掠等与壁湍流中相干结构的发展规律和现象十分相似。同时,通过对边界层中多涡结构的非线性演化问题的直接数值模拟,验证了该算法具有计算精度高,稳定性好,收敛速度快,出流边界影响小等优点。