Sajben跨声速扩压器分离流动中湍流模式数值研究

为研究湍流模式对激波/湍流边界层干扰内流流动的影响,提高数值计算准确度,使用SA,SST k-ω,非线性EASM k-ω,Gao-Yong四个湍流模式对Sajben扩压器内激波/湍流边界层干扰流动进行了数值计算。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。文中展示了四个湍流模式计算得到的壁面压力、速度剖面、摩阻系数等分布。计算值与实验值符合很好,四个湍流模式总体上能够较好地模拟扩压器内激波/湍流边界层干扰复杂分离流动。Gao-Yong湍流模式对分离区内的压力、速度型的模拟更加准确,而非线性EASM k-ω模式对分离再附点位置计算最理想。     

基于非交错有限体积格式和 k-ω 湍流模型模拟叶栅湍流分离流动

发展了一种求解任意坐标系下二维Ｎ－Ｓ方程的非交错有限体积法和一种高阶ＴＶＤ对流格式。采用一种ｋ－ω湍流模型和几种不同的ｋ－ω模型对叶栅湍流分离流动进行了详细的分析。数值结果与得到的实验结果比较表明,此方法不仅精度高,而且具有较好的数值稳定性。     

对 k-ω 湍流模型的改进

为了更适合于较复杂几何外形的数值计算，本文提出仅引入１个参数来改进ｋ－ω湍流模型并采用了一种双网格系统的有限体积格式来求解Ｎ－Ｓ方程与改进的的ｋ－ω模型方程。通过平板附面层的计算结果表明了用本文提出的对ｋ－ω湍流模型的改进对控制转捩是有效的。     

基于非交错有限体积式和k—w湍流模型模拟叶栅湍流分离流动

发展了一种求解任意坐标系下二维Ｎ－Ｓ方程的非交错有限体积法和一种高阶ＴＶＤ对流格式。采用一种ｋ－ω湍流模型和几种不同的ｋ－ω模型对叶栅湍流分离流动进行了详细的分析。数值结果与得到的实验结果比较表明,此方法不仅精度高,而且具有较好的数值稳定性。     

基于非结构/混合网格的脱体涡模拟算法

为了提高二阶精度有限体积算法的湍流数值模拟能力,在原始Roe格式基础上建立了与脱体涡模拟(DES)方法相匹配的二阶混合耗散自适应格式,根据流场信息自动调节格式耗散,并分别基于Spalart-Allmaras一方程湍流模型和k-ω剪切应力输运(SST)两方程湍流模型,发展了基于非结构/混合网格的DES方法。采用该方法计算了雷诺数为3 900的圆柱绕流和NACA 0021翼型60°大迎角分离流两个典型算例,通过与试验数据以及其他数值结果的对比验证了该方法的可行性。同时开展了不同数值格式、湍流模型的对比分析,研究结果表明:采用混合格式的DES算法能够解析更小尺度的湍流涡结构、计算数据更接近试验值;本文的DES类算法受其基准湍流模型影响较小。     

两种湍流模型在跨声速绕流计算的应用研究

湍流模型在对复杂流场的数值计算中起着非常重要的作用.本文采用SA一方程湍流模型和SST k-ω二方程湍流模型,通过解耦求解雷诺平均N-S方程和湍流模型方程实现对亚跨超声速湍流流场的数值模拟.对NACA0012翼型和ONERA-M6机翼跨声速绕流流场进行了计算,对压力分布和激波位置与实验结果进行了细致的比较,并分析了不同离散格式、不同网格疏密及壁面函数对计算结果的影响,在计算过程中这两个模型体现出了较好的简捷性和健壮性.     

高精度IPDG湍流模拟及通量格式数值特性北大核心CSCD

为更高精度数值求解复杂湍流问题,在内罚间断伽辽金(internal penalty discontinuous Galerkin,IPDG)方法框架内发展了基于SST k-ω模型的湍流模拟方法,通过对亚/跨/超声速工况下流场的数值计算与湍流特征捕捉,验证了方法的适用性,进而系统分析了AUSM、Lax-F、HLL和Roe 4种通量格式在IPDG湍流模拟中的数值特性。结果表明:AUSM格式在超声速工况下脱体激波面“褶皱”现象明显;Lax-F格式和HLL格式数值耗散大,在激波解析方面精度较低,且Lax-F格式在激波脚后诱导产生流动分离;Roe格式具备宽速域适用性,计算精度较高且能精确解析激波结构。     

跨声速下迎风格式的比较研究

研究了不同迎风格式的耗散机理,通过求解N-S方程,比较了其耗散性能.采用SST,EASM k-ω湍流模型,分析了湍流数值模拟中不同格式的计算精度,重点比较了迎风格式抑制壁面振荡的能力.选取跨声速收敛-扩张喷管内流场算例,从不同角度考察了格式的计算性能,数值计算结果和试验结果吻合,获得了不同格式耗散性能的比较结果.     

超声速横喷干扰湍流场数值模拟

本文探讨了超声速平板横向喷流干扰湍流流场的数值模拟方法。分别采用Roe、Harten-Yee及AUSM 格式对三维雷诺平均Navier-Stokes方程进行离散,湍流模拟应用Baldwin-Lomax代数模型,讨论了B-L模型中Fy最大值的不同选取方式。分析了数值格式、网格划分对计算结果的影响。将来流马赫数为5、喷流马赫数为3的实验结果与计算结果进行比较,选择了与实验符合较好的数值模拟方法,为进一步利用数值计算方法研究喷流干扰流场特性奠定基础。     

基于火焰面模型的超声速湍流燃烧数值模拟研究

在非结构混合网格有限体积框架下,将RANS方程和层流火焰面模型相结合,开展超声速湍流燃烧流动数值模拟研究,时间和空间离散分别采用LU-SGS格式和HLLC格式。用SST k-ω模型模拟湍流,FlameMaster 3.9生成火焰面数据库,并采用β-PDF和δ-PDF分布来考察湍流和化学反应的相互作用。采用该数值方法,对DLR氢燃料超燃燃烧室和北航乙烯燃料气动斜坡超燃燃烧室开展了数值模拟研究,并将数值计算结果与实验数据进行对比,二者符合较好。研究结果表明:在主流中心支板顺流喷射超燃燃烧室内,燃烧主要发生在支板后的燃料/空气混合层内,并且火焰核心落在恰当混合物分数Zst附近区域。      

复杂几何形状喷管内外三维流场的数值模拟

研究了球面收敛矩形扩张的二元喷管（ＳＣＦＮ）流场的数值模拟,为了有效地捕获激波与滑流面及较好地模拟分离,采用了Ｒｏｅ格式、ＭＵＳＣＬ方法及Ｖｅｎｄａｔａｋｒｉｓｈｎａｎ限制器、Ｓｐａｌａｒｔ－Ａｌｌｍａｒａｓ方程紊流模型。研究中给出了两个喷管模型四个流动状态的数值模拟结果。结果表明,除了有严重分离的流动情况以外,该算法能得到相当满意的计算结果。     

旋转离心叶轮内湍流的非结构化网格数值计算

给出了同位非结构化网格数值求解旋转离心叶轮内湍流流场的SIMPLE算法。网格划分采用改进的阵面推进法。控制方程采用通用的对流扩散方程并利用有限控制容积法离散。湍流模型采用标准的k-ε模型。形成的代数方程组用带有预条件器的共轭梯度平方法（CGS）求解。最后利用所给出的算法对一实验用的旋转离心叶轮进行计算。计算结果和实验数据基本相符,说明该算法合理,编程正确,可作为流场预测的理论工具。     

Gao-Yong理性湍流方程

对近十年来Gao-Yong理性湍流方程的研究进行综述,给出矢量形式的可压湍流封闭方程组的推导过程。侧偏统计平均保留了湍流脉动量的一阶统计信息,在引入加权漂移速度对称性及正交各向异性后,导出了漂移流的连续方程、动量方程及机械能方程,最后依据湍流物理的唯象论,使用动量传输链概念模化封闭了整个方程组。方程组不含任何经验系数,不使用壁面函数,保留了NS方程的均化的非线性特性。其级数形式的能量方程与非线性现象多尺度层次现象相对应,具备了描述湍流统计平均流动及拟序结构流动的双重功能。已进行了大量算例验证,验证结果证明了Gao-Yong湍流方程对广泛范围的复杂湍流问题的适定性。     

任意非正交曲线坐标系中湍流的数值计算

本文旨在运用张量分析理论,结合成熟的SIMPLE算法,发展一种适用于复杂几何形状内湍流流场数值求解的方法。假设流动为二维单连通域内的等温、不可压、定常、无化学反应、单组分湍流流动,采用双方程湍流模型。由雷诺平均的湍流控制方程的基本形式[1,2],结合张量分析理论。      

湍流的耗散及弥散相互作用理论

 湍流具有涡团散裂、耗散和弥散的特性。根据Prigogine倡导的耗散结构理论,推导了表征耗散与弥散相互作用的新的湍流控制方程组。其特点是:用稳定性分析得到湍流动能产生项,再根据广义熵增原理推出并列存在的,分别适用于强弱涡量的二个湍流动量方程。运用该理论已成功地计算了一些典型的湍流问题:湍流边界层中的马蹄涡拟序结构、钝体涡尾区的湍流能量逆转、湍流涡团散裂弛豫及各向异性分布。文中给出了部分算例。     

单列叶栅紊流流场有限分析数值研究

本文将有限分析方法用于曲线座标系上紊流 N-S方程数值计算,研究了高雷诺数时叶栅粘性紊流流场。有限分析方法在网格单元上对非线性偏微分程进行线化处理,在解析边界条件下求出线化偏微分方程的解析解,以此解构造离散方程。有限分析方法能够根据对流的方向和大小自动改变格式系数,并具有数值扩散小、精度高和稳定性好等优点。本文以 k-ε紊流模型模化紊流,以壁面函数方法处理近壁区流动参数。      

N-S方程有限分析数值研究

 本文将有限分析方法用于曲线座标系上紊流N-S方程的数值计算。分别计算了单列和串列叶栅内部流场,计算中采用了k-ε紊流模型和壁面函数。计算结果与试验结果相比较,吻合程度令人满意。有限分析方法在网格单元上对N-S方程进行线化处理,以解析边界条件作为约束,得出解析解,在解析解基础上构造离散代数方程。有限分析方法的最大特点是可以适应对流速度大小和方向,自动调整格式系数,因而具有数值扩散小和稳定性高等优点。     

有侧壁干扰的机翼半模型风洞实验的N-S方程数值模拟

应用三维可压、雷诺平均Navier-Stokes方程数值模拟了机翼半模实验风洞侧壁干扰和三维机翼半模与安装侧壁结合部流场。计算采用中心有限体积多步Runge-Kutta时间步长格式,将Baldwin-Lomax两层代数紊流模型扩展应用于三维拐角区流动。用本文方法计算了3个算例 ,并与国外的有关实验数据进行比较,计算和实验吻合良好。     

离心叶轮内稀释固液两相流动的数值模拟

 对固液两相流动进行了数值模拟。液相和固相分别采用 k- ε湍流模型和代数颗粒流( AP)湍流模型并考虑了液固两相间的相互作用。用此模型结合 SIMPLE- C算法对离心叶轮内的稀释固液两相流动进行了计算并研究了颗粒对流场的影响。计算结果同试验结果基本一致     

多级轴流压气机的径向掺混

 随着压气机向着更高的级负荷和更低的展弦比方向发展,径向掺混问题在气动分析和气动设计中越来越重要,当前对决定径向掺混的主要因素众说不一。本文用严格的理论,得到了用统一形式表示含二次流和湍流在内的各个径向掺混因素的多级轴流压气机通流基本方程组,并讨论了二维通流流场湍流掺混与实际流场三维剪切结构的关系。基本方程组的数值结果与试验结果以及其它掺混理论结果比较表呀,该方程组能够满意地模拟多级轴流压气机中的径向掺混特性。