高分辨率格式在非定常无粘可压缩流动中的应用研究

本文发展了一种计算非定常无粘可压缩流动的高分辨率格式。通量计算采用AUSMPW+Rie-mann求解器,其中使用了基于单调性的压力加权函数f,并与Roe Riemann求解器进行比较。为提高精度,使用五阶WENO格式进行左右状态插值,且时间推进采用三阶TVD Runge-Kutta方法。用该方法对移动接触间断问题、Sod激波管问题、二维激波反射问题和双马赫反射问题进行了数值计算,数值结果表明基于五阶WENO插值的AUSMPW+格式有很高的激波及接触间断分辨能力,并比基于五阶WENO插值的Roe格式有更高的计算效率。     

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法.外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格-库塔方法.采用四阶龙格-库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率.使用非结构动态网格技术对积冰后的翼型网格进行快速调整.在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性.     

重力波波包在时变背景场中传播的数值研究

采用二阶算子分裂格式对高斯型重力波波包在时变背景温场和时变背景风场中的非线性传播过程分别进行了数值模拟和分析.研究表明,时变背景场对波包空间结构演变影响不明显,各物理量的偏振关系与线性理论偏振关系较一致;在时变背景温场中传播,水平群速度和垂直群速度的变化趋势与线性理论结果相反,波频率的变化也与线性理论色散方程计算结果有一定的差别,这表明了线性重力波理论的局限性;在时变背景风场中传播,波频率的变化比时变背景温场中更显著.传播路径与在稳恒背景场中比较,更接近线性理论计算出的路径.      

一般状态方程可压缩多流体流动的PPM方法

采用流体体积分数的混合型多流体数值模型和高精度多介质Piecewise Parabolic Method(PPM)算法模拟一般形式Mie-Grüneisen状态方程的复杂多流体可压缩流动.根据多流体接触界面无振荡原则设计高精度计算格式,采用双波近似和两层迭代算法求解一般状态方程的Riemann问题,利用Strang分裂方法将一维算法直接推广到高维情况,不需要动能校正.二维多流体接触间断平移问题的模拟结果表明本算法在接触间断附近压力和速度没有振荡,界面数值耗散很小,被控制在两到三个网格之内.本文还模拟了二维Riemann问题和激波同多流体交界面的相互作用,算例结果表明本文方法可以有效地处理接触间断、激波等物理问题,且具有耗散小精度高的特点.     

跨声速压气机数值模拟格式效应研究

针对跨声速叶轮机械复杂内部粘性流动,开发了全三维高精度RANS数值求解器。分别采用FVS格式与AUSM+格式结合多种限制器对Rotor 67跨声速压气机转子进行数值模拟,分析和比较了各格式的计算效果并与试验结果对比。结果表明:与FVS格式相比,AUSM+格式的数值粘性更小,边界层的模拟精度更高;Hemker限制器的综合表现最优,Van Albada限制器的粘性分辨率略低于Hemker限制器,Minmod限制器对流动分离现象的捕捉能力较差,Van Leer限制器容易引入色散误差。     

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

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

沿流向微结构沟槽流场直接数值模拟

采用高阶格式对覆有V型对称沟槽表面的槽道湍流流动进行了直接数值模拟,数值方法采用有限差分法。为精确求解沟槽壁面的湍流流动,对流项的离散采用7阶WENO（Weighted Essentially Non-Oscillatory）格式;时间推进采用分数步时间推进与低耗散、低色散Runge-Kutta方法（LDDRK方法）结合的格式;黏性项的离散采用6阶中心格式。模拟的雷诺数为5 000（基于槽道高度的1/2）,计算的沟槽宽度范围为13~44,沟槽斜壁与水平面夹角为60°。数值模拟结果表明,与平板相比,沿流向沟槽表面的阻力最高降低了9%。数据分析发现出现减阻效果时,沟槽减少了近壁面处顺流向涡的数目,并且减阻机理与微沟槽阻碍大尺度流向涡与沟槽壁面的直接碰撞,使沟槽表面湍流脉动得到削弱有关。     

特定条件下高阶WENO格式计算结果误差

通过比较一阶迎风格式和五阶WENO格式模拟激波、接触间断、曲线坐标系下的均匀流和激波正规反射等4个简单流场得到的数值结果,发现WENO格式模拟的激波和接触间断在从初始间断变化成数值过渡区的过程中出现的非物理波动比一阶迎风格式的结果更加明显,流场结构也更加复杂;同时,由坐标变换而产生的几何诱导误差和边界近似模型误差也明显...     

隐式多重网格法求解叶轮机械三维跨声速湍流流场

本文提出了一种高效率求解叶轮机械三维跨声速湍流流场的隐工多重网格算法。它以矢通量分裂格式为基础，采用ＦＡＳ型多重网格加速技术。     

基于高阶耗散紧致格式的GMRES方法收敛特性研究

计算效率较低是当前限制高阶精度计算方法应用的重要因素。为了提高高阶精度混合型耗散紧致格式（HDCS）的计算效率,发展了适合多块对接网格的广义最小残值（GMRES）方法,并利用GMRES方法开展了HDCS格式的加速收敛研究。首先研究了GMRES的预处理方法、CFL数和内层迭代步数对HDCS数值模拟收敛特性的影响,计算结果显示：点松弛方法是一种高效的预处理方法;CFL数对计算收敛速度影响较大;GMRES方法存在最优的内层迭代步数。利用GMRES方法完成了NACA 0012翼型绕流、NLR 7301翼型绕流和DLR-F4翼身组合体绕流的数值模拟,并与其他隐式时间推进方法进行了对比,GMRES方法计算更加稳定,并且计算效率相对LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法可以提高5倍以上。研究结果表明,本文发展的GMRES方法在多块对接网格中具有良好的计算稳定性,计算结果的残差可以收敛到更低的量级,并且可以较大幅度地提高高阶精度数值模拟的计算效率。