• A Parallel Domain Decomposition Method for 3D Unsteady Incompressible Flows at High Reynolds Number

    分类: 物理学 >> 普通物理:统计和量子力学,量子信息等 提交时间: 2018-03-15

    摘要: Direct numerical simulation of three-dimensional incompressible ows at high Reynolds number using the unsteady Navier-Stokes equations is challenging. In order to obtain accurate simulations, very #12;ne meshes are necessary, and such simulations are increasingly important for modern engineering practices, such as understanding the behavior around high speed trains, which is the target application of this research. To avoid the time step size constraint imposed by the CFL number and the #12;ne spacial mesh size, we investigate some fully implicit methods, and focus on how to solve the large nonlinear system of equations at each time step on large scale parallel computers. In most of the existing implicit Navier-Stokes solvers, segregated velocity and pressure treatment is employed. In this paper, we focus on the Newton-Krylov-Schwarz method for solving the mono-lithic nonlinear system arising from the fully coupled #12;nite element is cretizationof the Navier-Stokes equations on unstructured meshes. In the subdomain, LU or point-block ILU is used as the local solver. We test the algorithm for some three-dimensional complex unsteady ows, including ows passing a high speed train,on a supercomputer with thousands of processors. Numerical experiments show that the algorithm has superlinear scalability with over three thousand processors for problems with tens of millions of unknowns.