基于有限差分-谱方法的分数阶Burgers流体的非稳态驻点流动

研究了分数阶Burgers流体通过拉伸平板的非稳态驻点流动问题。将分数阶导数引入Burgers流体模型可以更好地模拟流动过程,但也增加了模型的复杂性和求解难度。首次运用有限差分-谱方法求解分数阶Burgers流体模型,离散格式构造简单有效。采用谱方法对控制方程中的空间项进行离散,利用有限差分方法分别结合L-1和L-2算法离散控制方程中的时间项,给出了两种离散格式,并且通过构造数值算例证明了离散格式的收敛性。结果表明,在靠近平板处,速度随着分数阶导数的增加而减小,而无穷远处的流体速度呈现出相反的趋势,体现了分数阶导数的记忆特性。此外,雷诺数越小,流体的粘度越大,导致流体速度越大。由于松弛时间参数的松弛特性,靠近平板处松弛时间参数对速度分布有抑制作用,远离平板处松弛时间促进流体流动。

基于非结构/混合网格的高阶精度格式研究进展

尽管以二阶精度格式为基础的计算流体力学(CFD)方法和软件已经在航空航天飞行器设计中发挥了重要的作用,但是由于二阶精度格式的耗散和色散较大,对于湍流、分离等多尺度流动现象的模拟,现有成熟的CFD软件仍难以给出满意的结果,为此CFD工作者发展了众多的高阶精度计算格式.如果以适应的计算网格来分类,一般可以分为基于结构网格的有限差分格式、基于非结构/混合网格的有限体积法和有限元方法,以及各种类型的混合方法.由于非结构/混合网格具有良好的几何适应性,基于非结构/混合网格的高阶精度格式近年来备受关注.本文综述了近年来基于非结构/混合网格的高阶精度格式研究进展,重点介绍了空间离散方法,主要包括k-Exact和ENO/WENO等有限体积方法,间断伽辽金(DG)有限元方法,有限谱体积(SV)和有限谱差分(SD)方法,以及近来发展的各种DG/FV混合算法和将各种方法统一在一个框架内的CPR(correction procedure via reconstruction)方法等.随后简要介绍了高阶精度格式应用于复杂外形流动数值模拟的一些需要关注的问题,包括曲边界的处理方法、间断侦测和限制器、各种加速收敛技术等.在综述过程中,介绍了各种方法的优势与不足,其间介绍了作者发展的基于"静动态混合重构"的DG/FV混合算法.最后展望了基于非结构/混合网格的高阶精度格式的未来发展趋势及应用前景.

风雷软件LES开发设计与验证

基于国家数值风洞风雷软件开源框架,设计开发了LES(Large Eddy Simulation)湍流模型,主要包括Fourier谱/有限差分方法 LES求解器和有限体积/有限差分方法 LES求解器。简要介绍了采用的不可压缩流动求解的投影法、Fourier谱/有限差分混合方法、亚格子模型等理论方法,给出设计的软件框架和计算流程,尤其说明采用的松/紧2种耦合模式。通过数值模拟不可压缩槽道湍流、亚临界雷诺数圆柱绕流、NACA0012临界攻角的低频振荡算例,验证求解器的计算精度和复杂湍流模拟能力。基于风雷开源框架设计的开源LES模型,具备高精度数值格式、亚格子模型、湍流统计等通用模块,可为国内学者提供一个LES湍流模拟研究的开放平台。

三维尾迹型流动中的大尺度旋涡位错形成

应用紧致有限差分———Fourier谱方法求解三维不可压Navier Stokes方程组 ,数值研究尾迹型流动的三维演化和旋涡位错 .计算中在来流剖面中引入局部展向非均匀性 .数值结果表明流动的不稳定性发展导致三维旋涡涡街流场 .局部非均匀性引起涡街间 (沿展向 )在频率、相位及强度上的差异 .在非均匀区中产生并形成大尺度链状旋涡位错结构 .用数值模拟详细描述了旋涡位错的形成与特征 .