基于最小二乘无网格方法的大涡模拟研究  

Large eddy simulating based on least square meshless method

作　　者：王园丁[1] 谭俊杰[1] 蔡晓伟[1] 任登凤[1] 

机构地区：[1]南京理工大学能源与动力工程学院,江苏南京210094

出　　处：《华中科技大学学报（自然科学版）》2013年第4期103-108,共6页Journal of Huazhong University of Science and Technology(Natural Science Edition)

基　　金：国家自然科学基金资助项目(11072114)

摘　　要：对基于最小二乘无网格方法的大涡模拟进行研究,采用Smagorinsky亚格子模型封闭滤波后的N-S控制方程.针对无网格方法的特点提出适合无网格方法的空间滤波尺度求解形式,用AUSM+格式求解数值通量,最小二乘法拟合空间导数,三阶龙格-库塔法进行时间迭代.在此基础之上对Re=22 000状态下的二维方柱绕流问题进行数值模拟与分析.讨论了初始条件对计算结果的影响,统计了斯特劳哈尔、平均升阻力系数、方柱后方中心线上的平均速度及其均方根值以及不同位置处的平均速度及其脉动值沿Y方向的分布,并将统计结果与实验值及其他数值模拟结果进行对比,结果符合较好.Large eddy simulation （LES） based on the least square meshless method was researched. The subgrid scale model of Smagorinsky was employed to close the Favre-filtered Navier-Stokes equa- tions. The space filter length of the meshless method was proposed. The AUSM＋scheme was used to calculate the flux, the spatial derivative was fitted by the least square method, and the third order Runge-Kuta scheme was used for the time advance. Then, two dimensional flow past a square cylinder at Re= 22 000 was performed with uniform inflow condition. The influence of the initial condition was discussed. Some global quantities, such as Strouhal number, the mean values of drag and lift coeffi- cient were calculated. The mean and root-mean-square value （rms） values of velocity along the center line of the square cylinder and the distribution of the mean and rms values of velocity along the Y di- rection were also investigated. The results agree well with the experiments and other numerical simu- lation values.

关 键 词：大涡模拟 无网格 最小二乘 方柱绕流 数值模拟 湍流 

分 类 号：O35[理学—流体力学]