Large eddy simulations of incompressible turbulent flows using parallel computing techniques
Large eddy simulations of incompressible turbulent flows using parallel computing techniques
复制标题
DOI：
10.1002/fld.1560
复制
发表时间：
2008-04
影响因子：
1.8
通讯作者：
Anup Gokarn;F. Battaglia;R. Fox;James C. Hill;J. Réveillon
中科院分区：
文献类型：
--
作者：
Anup Gokarn;F. Battaglia;R. Fox;James C. Hill;J. Réveillon
本文介绍了在分布式存储机上用大涡模拟（LES）求解不可压缩高雷诺数湍流的详细过程。滤波后的Navier-Stokes方程采用部分交错变量排列进行离散化，并采用有限差分网格进行求解。空间导数采用二阶中心差分格式和六阶紧化格式。时间导数采用三阶低存储龙格-库塔法。数值方案的验证首先通过模拟驱动空腔流动和向后台阶的流动进行。然后对Smagorinsky动态亚网格湍流模型进行了通道内流动的验证。用文献中的相关数据验证了仿真结果。由于LES计算量大，求解器采用消息传递接口并行化。利用一种高效的平行线性方程求解器求解椭圆压力泊松方程。对该并行程序进行了复杂流态流动解的测试，并将初步结果与实验数据进行了比较。在多达256个处理器的并行集群上测试了相同几何形状的程序性能。这项工作中的新方法是对湍流的LES使用部分交错的可变排列，从而避免了可能掩盖子网格效应的任何形式的人工耗散的需要。版权所有©2007 John Wiley & Sons， Ltd
This paper presents a detailed procedure to solve incompressible high Reynolds number turbulent flows using large eddy simulations (LES) on distributed memory machines. The filtered Navier–Stokes equations are discretized using a partial‐staggered variable arrangement and solved using a finite difference grid. A second‐order central difference scheme and sixth‐order compact scheme are employed for the spatial derivatives. A third‐order low storage Runge–Kutta method is used for the temporal derivatives. Validation of the numerical scheme is performed first by simulating a driven cavity flow and flow over a backward‐facing step. The dynamic Smagorinsky subgrid turbulence model is then validated for flow in a channel. Simulations are validated with relevant data available in literature. Since LES is computationally expensive, the solver is parallelized using message passing interface. An efficient parallel linear equation solver is utilized for solving the elliptical pressure Poisson equation. The parallel program is tested for solutions of flow in a complex flow configuration and preliminary results are compared with experimental data. Performance of the program for the same geometry is tested on a parallel cluster up to 256 processors. The novel approach in this work is the use of a partial‐staggered variable arrangement for LES of turbulent flows, obviating the need for any form of artificial dissipation that might mask the subgrid effect on the solution. Copyright © 2007 John Wiley & Sons, Ltd.