中国科学院机构知识库网格
Chinese Academy of Sciences Institutional Repositories Grid
复杂可压缩流动高精度数值模拟方法研究

文献类型:学位论文

作者何康
答辩日期2024-11
文献子类博士
授予单位中国科学院大学
授予地点北京
导师李新亮 ; 刘洪伟
关键词混合动理学 WGVC方法 WENO方法 动态非定常流动 DP-LUR 双时间步方法 隐式大涡模拟
学位专业流体力学
其他题名Study of High Order Numerical Simulation Methods for Complex Compressible Flows
英文摘要

复杂可压缩流动中往往同时具有多尺度结构和间断结构,对数值格式和计算程序的设计都提出了苛刻的要求。对数值格式而言,一方面,多尺度结构要求数值格式具有良好的频谱特性(低耗散和低色散),以准确解析解的振幅和相位。另一方面,间断结构则要求数值格式具有足够的耗散,以捕捉间断结构而不产生数值振荡。对计算程序而言,实际流动的复杂性和工况的多样性要求计算程序不断拓展功能,以满足工程和学科发展的需求。本文针对复杂可压缩流动的高精度数值方法和计算程序开展了深入研究,主要研究内容和结论如下:

(1)提出了一种用于可压缩流动模拟的混合动理学通量WGVC(Weighted Group Velocity Control)-WENO(Weighted Essentially Non-Oscillatory)方法。该方法继承了文献\cite{liu2015}中提出的混合动理学方法,考虑了气体分子的自由传输和碰撞效应,并结合了文献\cite{he2014}中的WGVC-WENO重构技术。具有以下特点:在光滑区域,碰撞相关的通量起主要贡献,且具有较小耗散性,并采用WGVC重构方法;在非光滑区域附近,无碰撞的KFVS(Kinetic Flux Vector Splitting)通量起主要作用,具有较大耗散性,并与WENO重构技术相匹配。数值实验表明,该方法比文献\cite{liu2015}和其他FVS(Flux Vector Splitting)方法\cite{steger1981flux,harten1983upstream,kao2004lax}分辨率更高,耗散性更小。利用多个一维和二维算例,验证了该方法不仅具有高精度和低耗散特性,还展现出较强的鲁棒性和激波捕捉能力。

(2)为满足多尺度结构的分辨率要求,提出了一种优化的WGVC格式,该方法在保持精度的同时优化了数值格式的谱特性。通过设计波包的光滑度量因子和非线性加权方案,在低波数范围内实现精度控制,在中波数范围内实现群速度控制,最终显著提升了整体的谱特性。此外,通过将优化的WGVC格式嵌入到激波捕捉格式中,如WENO/TENO(Targeted Essentially Non-Oscillatory)格式,可以在增强格式激波捕捉能力的同时,保留了WGVC格式在中低波数的谱特性。理论和数值实验证明,优化的WGVC格式具有保持精度、色散和耗散误差小等优点,非常适合于湍流/激波边界层干扰等复杂流动问题的数值模拟。

(3)为实现动态非定常流动问题的高精度求解,基于OpenCFD-SC,开发了一款动态非定常高精度计算程序,程序利用刚性动网格和超限插值动网格实现了网格运动,利用局部时间步长方法和DP-LUR(Data-Parallel Lower-Upper Relaxation)方法实现高效的隐式并行计算。随后利用NACA0012翼型的俯仰振荡跨声速流动作为测试算例,对程序的正确性和收敛性进行了验证。计算结果与实验和他人数值模拟的结果基本一致,所开发的程序能够有效捕捉复杂的动态流动特征,为动态非定常流动问题的高精度数值模拟提供了研究基础。

(4)利用所开发的程序,开展了不同固定攻角和动态俯仰条件下的RAE2822超临界翼型的隐式大涡模拟研究。研究深入探讨了激波与湍流边界层的相互作用,揭示了随攻角变化激波强度和位置的运动规律,尤其在$5^\circ$攻角下,激波强度显著增强,且位置后移。动态俯仰运动中,激波呈现非标准的周期性运动,升力和力矩系数出现明显滞后。通过对比不同俯仰方向的流动特性和湍流特性,发现上仰过程加强了激波下游的流向速度脉动和正应力,削弱了其他方向的速度脉动及雷诺应力,而下俯过程则相反。通过对比不同攻角下的流动特性,发现动态流场特性受平均攻角影响很大,在合适的平均攻角和运动频率下,将有助于维持流动的稳定。

语种中文
源URL[http://dspace.imech.ac.cn/handle/311007/97814]  
专题力学研究所_高温气体动力学国家重点实验室
推荐引用方式
GB/T 7714
何康. 复杂可压缩流动高精度数值模拟方法研究[D]. 北京. 中国科学院大学. 2024.

入库方式: OAI收割

来源:力学研究所

浏览0
下载0
收藏0
其他版本

除非特别说明,本系统中所有内容都受版权保护,并保留所有权利。