用N－S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二级时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型￣［１］。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

用N—S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二维时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

细长体大攻角绕流的数值模拟

用数值模拟方法研究了细长体大攻角的绕流特性（求解在周向３６０°范围内进行）。数值模拟的出发方程是三维可压流的全Ｎ－Ｓ方程,格式为二阶迎风ＴＶＤ格式。通过数值模拟研究了附近,攻角１９°、３２．５°和５０°状态下,细长体绕流场中的激波、分离流动及表面压力分布等重要特征。     

高阶精度格式WCNS在三角翼大攻角模拟中的应用研究

采用5阶精度的加权紧致非线性格式(WCNS-E-5)数值模拟了65°三角翼的大攻角绕流流场,主要目的是考核高阶精度格式WCNS在大攻角旋涡流动方面以及跨声速流场的激波附面层干扰、涡破裂位置的模拟能力,重点研究不同网格规模和湍流模型对尖前缘三角翼涡系之间的相互作用的影响。通过求解任意坐标系下的雷诺平均N-S方程,采用5阶精度的加权紧致非线性格式(WCNS-E-5)和多块对接结构网格技术,两种湍流模型分别是一方程SA和两方程SST湍流模型,在与相应试验结果对比的基础上,详细研究了WCNS-E-5格式在跨声速大攻角旋涡流动中的表现,以及不同网格规模、两种湍流模型对主涡二次涡相互作用、涡破裂位置和表面压力分布的影响。本文的研究结果表明,高阶精度格式WCNS-E-5能成功应用于三角翼的跨声速大攻角流动,网格规模的增加进一步提高流场分辨率,SST湍流模型相对SA湍流模型在三角翼大攻角流动中具有更好的适用性。     

翼型大攻角绕流的数值模拟

以求解二维N-S方程数值模拟NACA0012翼型大攻角状态的可压绕流特性;N-S方程是在贴体坐标系中给出的,以代数方法生成C型网格系统。采用LU-ADI格式计算,为提高格式的稳定性在隐式和显式部分分别添加了2阶和4阶人工粘性项。应用BaldwinLomax湍流二层代数模型模拟了大攻角时前缘分离涡的形成,旋涡对流及其非定常现象。在某些Mach数和攻角下NACA0012翼型的湍流解具有周期性。通过比较,本文数值计算结果同实验及国外相应的数值计算结果基本吻合。     

尖拱弹身大攻角N-S方程数值计算

本文采用有限体积的NND格式,数值模拟了弹身超音速有攻角绕流,分析了周向网格和法向网格对弹身背风区分离流态及气动性能的影响,计算与实验结果比较,两者相当一致.     

基于无网格算法的化学非平衡流数值模拟

为了验证AUSMPW+格式在无网格算法的精度等以及无网格算法能否用于复杂化学非平衡流数值模拟中,首先将AUSMPW+格式推广到无网格算法中,给出了AUSMPW+在无网格算法中的具体形式。其次,将无网格算法应用于带化学反应的多组分气体Euler方程的求解,给出无网格条件下求解化学非平衡流的具体过程。最后,对NACA0012跨声速绕流场、高超声速弹丸绕流场、氢气/氧气诱导燃烧流场、球锥高超声速绕流场、楔体斜爆轰流场进行了数值模拟。数值模拟结果与相关文献计算结果和实验结果吻合较好,表明了AUSMPW+格式在无网格算法中表现较高精度。同时,采用无网格算法能较好地模拟复杂化学非平衡流场,正确分辨复杂的物理现象,为化学非平衡流数值模拟提供一种新的算法,拓宽了无网格算法的应用范围。     

机动飞行器中等攻角高超声速无粘绕流数值模拟

本文采用ＮＮＤ差分格式，利用推进－迭代方法，数值模拟了机动飞行器的中等攻角高超声速无粘绕流流场，发展了一种适用于处理中等攻下背风面流场的技术，大量计算表明，结果准确、可靠。气动力力系数和壁面压力分布与实验数据一致，流场波系结构模拟正确。     

跨声速钝体绕流数值模拟研究

通过数值求解ＮＳ方程，进行了三方面工作（１）格式和算法对返回舱跨声速流场特性的影响研究；（２）网格形态对捕捉声速弱激波的有效性研究；（３）模拟了无支杆和有支杆返回舱跨声大攻角绕流流场，给出了后球体上的跨声速激波和后支杆上的第二道激波。     

有侧滑角椭圆弹身大攻角状态绕流的数值模拟

用数值模拟方法研究了椭圆细长体在大侧滑角、大攻角状态下的流动特性。数值模拟的出发方程为三维可压缩流的全Ｎ－Ｓ方程，数值格式为ＴＶＤ格式。通过数值模拟，研究了椭圆截面细长体在Ｍ∞＝２．５、侧滑角β＝４５°状态下，攻角从１０°￣３５°范围内流场中的分离特征。     

分区对接网格在跨流域气体运动论统一算法的应用研究

为了模拟跨流域复杂物形绕流或多体流场,需要在位置空间建立能描述整个流场特性的网格系统,相应的数值计算方法也需要有所改进。本文采用分区对接网格处理复杂物形绕流或多体干扰流场,得到能准确描述流场特征的贴体网格系统。在此基础上数值求解考虑转动非平衡影响的Boltzmann-Rykov模型方程,在速度空间应用离散速度坐标法、位置空间采用数值差分NND格式。网格对接面上考虑分布函数的传递,包含网格点信息和网格点上速度的信息。通过计算得到复杂物形的流动特性。相关算例表明网格分区之间的流场参数光滑过渡。与文献及DSMC结果的对比证实本文基于分区对接网格的气体运动论统一算法是可靠的。     

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

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

分区对接网格在跨流域气体运动论统一算法的应用研究北大核心CSCD

为了模拟跨流域复杂物形绕流或多体流场,需要在位置空间建立能描述整个流场特性的网格系统,相应的数值计算方法也需要有所改进。本文采用分区对接网格处理复杂物形绕流或多体干扰流场,得到能准确描述流场特征的贴体网格系统。在此基础上数值求解考虑转动非平衡影响的Boltzmann-Rykov模型方程,在速度空间应用离散速度坐标法、位置空间采用数值差分NND格式。网格对接面上考虑分布函数的传递,包含网格点信息和网格点上速度的信息。通过计算得到复杂物形的流动特性。相关算例表明网格分区之间的流场参数光滑过渡。与文献及DSMC结果的对比证实本文基于分区对接网格的气体运动论统一算法是可靠的。     

超音速来流中侧向喷流干扰流场的数值模拟

通过求解三维可压缩N S方程组,对尖锥/圆柱组合体模型在超音速绕流中的侧向喷流干扰流场进行了模拟,数值方法采用有限体积法,差分格式采用NND格式。分析了干扰流场的结构,研究了干扰效应对气动性能的影响,得到来流Ma∞=3 3和4 5时力放大因子和干扰力矩系数随攻角变化的规律,计算结果与实验一致。     

非结构混合网格鲁棒自适应技术

数值格式、湍流模型和计算网格是影响CFD数值模拟精度的3个主要因素。结合流场信息的网格自适应技术具备动态优化计算网格的能力,被NASA列为未来CFD发展的一项关键技术。本文针对非结构混合网格,发展了网格单元分布优化、表面网格几何投影和空间网格协调匹配3项关键技术,建立了高鲁棒性几何保真的网格自适应系统。首先,为了提高自适应方法的鲁棒性和通用性,发展了基于标准面网格的多面体网格单元分布优化方法。其次,发展了仅依赖表面网格信息的局部曲面重构技术,采用参数点映射方法实现了新增表面网格点的几何投影,消除了自适应系统对几何CAD系统的依赖。再次,采用改进的距离函数方法实现了空间网格与投影后表面网格的快速匹配。最后,结合基于流场特征的自适应探测器,采用二阶格式的有限体积方法,开展了30P30N三段翼绕流和三角翼大迎角绕流的网格自适应数值模拟。结果表明,通过网格自适应对网格单元的分布进行优化后,流场求解的收敛性和模拟精度都得到了显著提高。     

用非结构动网格技术计算确定再入飞行器配平攻角

提出了用数值计算直接确定飞行器配平攻角的自由振动法.该方法采用非结构动网格技术,使飞行器在气动力作用下绕通过质心的轴线自由振动,最终衰减到配平攻角.构造了一种MUSCL类型有限体积格式,求解考虑了动网格效应的三维非定常Euler方程,采用弹簧近似技术实现非结构网格的变形运动.作为例证,计算确定了载人飞船返回舱的配平攻角,计算结果与实验结果接近.该方法外形适应能力强,效率较高,具有较好的应用前景.     

头部几何微扰对流动不对称性影响的数值研究

使用低耗散的Roe格式,数值模拟了头部几何微扰对细长体不对称绕流的影响.计算发现,周向网格的加密本身也是一种微扰动,会对流场产生很大的影响.如果再加上头部几何微扰,随其高度的变化,侧向力分布显示出类似于几何扰动周向角变化时产生的双稳态现象.在对应的攻角下,头部扰动的方式决定了头部分离的方向,进而决定了后体的流场发展情况.     

低雷诺数下二维翼型绕流的流场特性分析

采用高精度有限差分格式,对低雷诺数下二维翼型绕流进行了直接数值模拟,计算了雷诺数为1.0×104,NACA0012翼型0°,4°以及10°攻角下的流场。计算结果表明:在0°和4°攻角条件下,翼型绕流尾迹区的统计特性相似,0°攻角下的统计量值具有很好的对称性;在距翼型尾缘0.3弦长以后的尾迹区,旋涡排列成类似涡街的结构,涡量的极值、压力的极小值都位于旋涡中心,沿着流向,涡量极值的绝对值逐渐减小,压力的极小值逐渐增大。10°攻角下,翼型上表面从前缘开始分离,尾迹区统计分析结果所得图象与0°和4°完全不同,数值上较后者结果大;在翼型尾缘处,涡量的卷吸,高压力区域的形成,是旋涡脱落的条件,正向和反向旋涡的交替脱落,形成了类似涡街的结构。      

超声速有翼飞行器的无粘分区数值模拟

本文采用ＭａｃＣｏｒｍａｃｋ两步显式格式，用空间推进法求解Ｅｕｌｅｒ方程来对有翼飞行器的超声速绕流流场进行数值模拟，应用特征边条件处理场面边界及装配激波的计算，在翼缘及表面斜率不连续处，采用局部解析处理，提高了算法的可靠性，复杂外形的计算采用分区处理技术，网格由代数形式生成，在大攻角的计算中，通过对横流速度型的修正改善了背风面压力分布特性，解决了推进中断问题，同时模拟出背风面涡。     

含侧喷流弹丸三维绕流干扰流场数值分析

从全Navier-Stokes方程出发，应用三维隐式有限体积TVD格式，在超声速和零攻角条件下，分别对无侧喷流和含侧喷流的低阻远程弹丸三维粘性绕流流场进行了数值模拟，为避免分区计算而用代数方法生成了弹丸绕流场O型网络，计算结果得到了清晰的流场波系结构，并从流场波系结构、弹体和弹底表面压力分布等角度将两者进行了对比与分析，表明侧喷流不仅提供增程推力，而且改变了弹丸阻力特性，其结果可用于研究弹丸的火箭增程技术。