基于DUGKS的槽道湍流大涡模拟研究  

LARGE EDDY SIMULATION OF TURBULENT CHANNEL FLOW USING DUGKS-LES

作　　者：郭琳 沈杰 彭程 郭照立[3] 王连平 牛中国[5] Guo Lin;Shen Jie;Peng Cheng;Guo Zhaoli;Wang Lian-Ping;Niu Zhong-Guo(Department of Mechanics and Aerospace Engineering,Southern University of Science and Technology,Shenzhen 518055,Guangdong,China;Key Laboratory of High-efficiency and Clean Mechanical Manufacture,Ministry of Education,School of Mechanical Engineering,Shandong University,Jinan 250100,China;Institute of Interdisciplinary Research for Mathematics and Applied Science,School of Energy and Power Engineering,Huazhong University of Science and Technology,Wuhan 430074,China;Guangdong Provincial Key Laboratory of Turbulence Research and Applications,Center for Complex Flows and Soft Matter Research,Southern University of Science and Technology,Shenzhen 518055,Guangdong,China;Low Speed High Reynolds Number Aeronautical Key Laboratory,AVIC Aerodynamics Research Institute,Harbin 150001,China)

机构地区：[1]南方科技大学力学与航空航天工程系,广东深圳518055 [2]山东大学机械工程学院,山东大学高效洁净机械制造教育部重点实验室,济南250100 [3]华中科技大学能源与动力工程学院,华中科技大学数学与应用学科交叉创新研究院,武汉430074 [4]南方科技大学复杂流动及软物质研究中心,广东省湍流基础研究与应用重点实验室,广东深圳518055 [5]中国航空工业空气动力研究院低速高雷诺数航空重点实验室,哈尔滨150001

出　　处：《力学学报》2024年第12期3394-3408,共15页Chinese Journal of Theoretical and Applied Mechanics

基　　金：国家自然科学基金(U2241269,T2250710183);广东省湍流基础研究与应用重点实验室项目(2023B1212060001)资助。

摘　　要：基于介观数值模拟方法的大涡模拟(large eddy simulation,LES)是研究湍流的重要手段之一.离散统一气体动理学方法(discrete unified gas kinetic scheme,DUGKS)是一种相对较新的介观方法,其作为湍流大涡模拟工具的研究还处于起步阶段.文章通过耦合标准离壁Smagorinsky亚格子应力模型和Musker壁面模型,实现了DUGKS-LES介观大涡模拟方法.推导了在DUGKS中实现大涡模拟所需的有效总应力和修正分布函数非平衡态的关系式;通过壁面函数逆求解壁面剪切应力时,设计了可以保证收敛性的快速牛顿迭代法,并利用逆问题的近似解进一步提高了迭代效率.最后,分别对摩擦雷诺数为180,395和950的槽道湍流进行大涡模拟,并与相关文献中介观大涡模拟方法的结果进行了比较.结果表明,DUGKS-LES方法在准确性上具有优势,以谱方法模拟结果为基准,流向速度、雷诺应力和脉动速度等湍流统计量的均方误差均可以控制在11%以内.在DUGKS框架下,与文献中利用自适应亚格子应力模型相比,文章结合Musker壁面模型的处理方法,可以得到更精确的流向速度和雷诺切应力.The mesoscopic-scheme-based large eddy simulation(LES)has emerged as an alternative approach to study turbulent flows.The discrete unified gas kinetic scheme(DUGKS)is a relatively newly developed mesoscopic method,and the realization of LES in DUGKS is still in its early stages.In this paper,a DUGKS-LES algorithm is proposed by coupling the standard off-wall Smagorinsky subgrid model and the Musker wall model.Firstly,the relation between the effective total stress and the non-equilibrium distribution function for LES within the DUGKS framework is derived.When solving the wall shear stress using the wall function,a rapid Newton iteration method is designed to ensure the convergence,by incorporating the approximate solution of the inverse problem.Finally,large eddy simulations of the turbulent channel flows at friction Reynolds numbers of 180,395 and 950 are performed to validate the proposed DUGKS-LES method.Compared with the results from the spectral method,the mean square errors of turbulence statistics such as streamwise velocity,Reynolds stress and fluctuating velocity are maintained below 11%.In the DUGKS framework,when compared with the wall-adapting subgrid models used in the literature,the proposed DUGKS-LES method incorporating the Musker wall model can achieve more precise results in terms of mean streamwise velocity and Reynolds shear stress.

关 键 词：离散统一气体动理学方法 大涡模拟 壁面模型 槽道湍流 

分 类 号：O357.5[理学—流体力学]