高超声速边界层感受性问题高精度数值模拟

通过五阶WENO格式和六阶对称紧致格式以及三阶TVD R-K结合的方法,对存在强激波和小扰动相互干扰的高超声速边界层感受性问题进行研究.结果表明:此方法能够模拟边界层内不稳定波的产生和发展以及小扰动和激波相互干扰等现象,因此能广泛应用于含激波的感受性等问题的可压缩湍流直接数值模拟以及具有激波和边界层干扰等复杂流场的计算.     

沟槽型面减阻特性数值模拟研究

基于人工压缩性方法,时间项采用隐式离散、对流项采用基于Roe近似Riemann解的迎风格式,研究了三维不可压粘性流动的数值模拟方法,湍流模型采用Baldwin-Barth一方程模型,利用所发展的不可压粘性流数值模拟方法,对典型几类V型、间隔三角型和U型沟槽型面的不可压粘性流场进行了数值模拟计算研究,探索了上述几类沟槽型面参数对湍流减阻特性的影响规律,从表面摩阻、速度场等方面分析了上述几类沟槽型面的减阻作用机理。研究结果表明：上述典型沟槽的减阻效果可以分别达到7．6％,10．5％和5．0％,与国内外相关实验结果符合一致,所发展的计算方法是成功的、可行的。     

高阶气体动理学格式在湍流数值模拟中的应用

本文回顾了高阶气体动理学格式在湍流数值模拟中的应用.与传统的Riemann求解器相比,气体动理学格式可以提供时空耦合的演化过程,这对发展高精度格式十分重要.因此,基于两步四阶时间离散和高精度WE-NO重构,发展了具有四阶时间精度的气体动理学格式.该格式有更高的数值精度和稳定性,并且具有更好的处理复杂流动问题的能力.目前...     

沟槽面湍流减阻数值评估方法

为解决雷诺平均Navier-Stokes(RANS)方法模拟沟槽面减阻的实际应用问题,受Wilcox对粗糙度模化方法的启发,通过深入分析k-ω剪切应力输运(SST)湍流模型近壁区ω值的作用效果,发现在黏性底层内增加壁面ω值可使对应近壁区的湍动能、湍流黏度、雷诺应力均减小,这种对近壁区流动特性的作用效果与真实沟槽面一致。以经典的对称V型沟槽面(高度=间距)减阻实验数据为基础,通过引入减阻效果影响因子——压力梯度与偏航角,建立复杂流动条件下沟槽面几何尺寸到壁面ω值的模化函数方程,将数值模拟结果与实验结果对比,验证了模化沟槽面可以达到与真实沟槽面相同的减阻效果,并依此给出沟槽面减阻在工程实践中的应用方法,以C919翼身组合体巡航构型为例,完成其沟槽面减阻的总体设计与评估。     

低雷诺数NACA0012平面叶栅流场直接数值模拟

采用具有7阶精度的weighted essentially non-oscillatory(WENO)差分格式,直接求解可压缩二维非定常N-S方程组,研究了NACA0012翼型平面叶栅低雷诺数流动的特征.直接模拟及与文献对比的结果表明:叶栅尾缘涡脱落的形成过程与圆柱绕流涡脱落的形成过程非常相似.平面叶栅尾迹区的2阶统计量与孤立翼型尾迹区的2阶统计量具有相同的分布特征,但前者的强度显著大于后者.周期性的涡脱落不仅在上下翼面形成非定常分离,也使得尾迹区某点的总压发生准周期性的变化.随着栅距的减小,翼型上的平均分离位置向前缘移动;尾迹区某点的总压变化频率及其幅值均显著地增加;而且栅距越小,速度脉动2阶统计量反而越大.      

双侧进气突扩燃烧室流场数值模拟

采用有限体积法和SIMPLE算法对冲压发动机突扩燃烧室流场进行了数值模拟,湍流模型采用标准二方程湍流模型,燃烧模型采用EDM模型。给出了流场中速度、速度流函数、湍流强度以及温度的分布。分析了流场特征,并与无燃烧状态下的流场进行了比较。     

高超声速湍流直接数值模拟技术

概述了近年来国内外高超声速湍流直接数值模拟(DNS)技术方面最新的研究进展,主要集中在高精度、高鲁棒性数值方法方面,同时也介绍了近年来典型的高超声速湍流DNS算例。在数值方法方面,主要介绍了高精度激波捕捉格式以及保持计算稳定的数值技术,重点是WENO格式及高阶保单调格式的最新进展。在高超声速湍流DNS算例方面,介绍了压缩性影响、壁温影响、真实气体效应以及高超声速转捩等方面的DNS研究。此外,还简要介绍了作者开发的可压缩高精度计算流体力学软件OpenCFD。     

平面压气机叶栅通道分离流动高精度直接数值模拟

采用高精度有限差分格式直接求解二维Navier-Stokes方程组,数值模拟V103平面压气机叶栅分离流动,数值结果表明:在瞬时流场中,吸力面后部发生流动分离,在分离区前端存在大尺度分离涡,分离涡下游是由二次涡和脱落涡交替形成的涡串,直至叶片尾缘,形成以脱落涡为主结构的尾迹;在时均流场中,吸力面后部存在短分离泡,分离区压力分布存在明显压力平台。与逆压力梯度下平板边界层分离流动相比,瞬时和时均流场结构相似;叶栅通道内无量纲涡脱落频率是前者的两倍。与文献计算结果对比表明:叶片表面压力分布除分离区外吻合很好;非定常计算所得分离区轴向长度比定常计算大41%。在分离区内三个二阶统计量均达到最大值,表明流场强非定常性集中在分离区。     

区域分裂法在湍流流场数值解中的研究

应用区域分裂算法(简称DDM)对湍流流场进行数值研究。对湍流流动控制方程建立一种通用的DDM-1计算格式S_{I,ω,给出重叠区域内变量处理方法。DDM-1的SI,ω计算格式的适用性由突扩圆管分离流流场数值计算所验证,数值计算结果与实验数据比较表明:采用SI,ω计算格式所预测的再附点位置及平均轴向速度、湍动能强度比在整个计算区域内采用单一模型的方法能够更好地与实验值吻合。}      

湍流燃烧的动态自适应直接数值模拟

为了能以较低的计算成本分辨多尺度的湍流燃烧过程,发展了基于动态自适应网格技术的湍流燃烧直接数值模拟方法。在氢气/空气层流预混火焰传播的模拟中,准确计算出了火焰传播速度,表明该方法具有直接求解火焰以及精确捕捉火焰传播过程的能力。在衰减的各向同性均匀湍流的模拟中,借助FFT分析获得了合理的湍动能谱线随湍动能衰减的变化规律,表明该方法能够直接求解湍流。由于在实际应用中,点火发生的位置与低温机理的速率相关,因此在氢气/空气湍流预混火焰中研究了湍流对低温机理的影响,发现湍流输运以及湍流涡团间强烈的拉伸和剪切对低温机理有明显的加速作用,能显著增强放热。     

脊状表面流场数值计算方法研究

根据现有湍流数值模拟方法的特点和脊状表面减阻理论研究的实际需要,研究了基于雷诺平均N-S方程的脊状表面流场数值计算方法,并基于纵向及横向脊状表面的形状特点,提出了不同脊状表面计算域选取及网格剖分的方法。一系列典型横向及纵向脊状表面流场数值计算结果表明,所采用的这套脊状表面数值计算方法不仅能够精确获得脊状表面的减阻特性,还可以获得脊状结构内部及附近流速、压强、涡量等细微流场结构的分布规律。     

基于湍流边界层时均速度分布的脊状表面减阻规律研究

通过对试验测得的脊状表面湍流边界层时均速度分布的分析,获得了脊状表面的减阻范围和基本规律.脊状表面理论零点及壁面摩擦速度的计算采用基于湍流边界层Spalding壁面律公式的最小二乘拟合方法;模型板摩擦阻力计算采用基于边界层动量积分方程的方法.研究表明:脊状结构使得湍流边界层粘性底层增厚,近壁面法向速度梯度降低,过渡层与对数律区上移;V型脊状结构的无量纲宽度s+介于6至18时具有减阻效果;当s+≈12时,减阻量最大,最大值约7.7%.      

微射流作动器三维外流场数值模拟

采用有限体积法求解Faver平均N—s方程，对三维、粘性、非定常、不可压微作动器外流场进行数值模拟。计算结果揭示了微射流产生、发展与耗散的过程，和实验中观察到的现象比较吻合。旋涡对沿展向的不稳定性是涡环破裂、旋涡对的相干性消失的主要原因。二次涡的发现为以后研究微射流旋涡对的细微结构打下良好的基础。     

A novel high-order scheme for numerical simulation of wake flow over helicopter rotors in hover

Accurate prediction of tip vortices is crucial for predicting the hovering performance of a helicopter rotor. A new high-order scheme(we call it WENO-K) proposed by our research group is employed to minimize numerical dissipation and extended to numerical simulation of unsteady compressible viscous flows dominated by tip vortices over hovering rotors. WENO-K is referred to as an adaptively optimized WENO scheme with Gauss-Kriging reconstruction, and its advantage is to reduce dissipation in smoo...     

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

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

马赫数为6的高超声速钝锥湍流边界层空间演化的直接数值模拟

首先求解了来流马赫数为6的零攻角高超声速钝锥边界层的层流基本流场,在选定的计算域入口引入一组有限幅值的T-S波扰动,用高精度差分格式对流动进行了直接数值模拟.引入的扰动触发了转捩,从而得到了空间模式下的湍流边界层.研究了湍流平均场与脉动场的统计特性,给出了相干结构的流动显示图,并对强雷诺比拟的结论进行了检验.