Multi-GPU parallel computation of unsteady incompressible flows using kinetically reduced local Navier–Stokes equations

Abstract Numerical simulations of 2D doubly periodic shear layers and 3D decaying homogeneous isotropic turbulence are presented using Kinetically Reduced Local Navier–Stokes (KRLNS) equations that is applicable to the unsteady incompressible flows without the need for sub-iterations and is capable of capturing the correct transient behavior. To achieve high accuracy, the KRLNS equations is discretized with higher order central difference approximations and 4-stage Runge–Kutta method. The results are compared with the solutions obtained by Lattice Boltzmann method (LBM) and pseudo-spectral method (PSM), which is the standard approach for this problem. Parallel computations are carried out on multiple GPUs (Tesla K40), maximum 4 GPUs available, based on the domain decomposition method and the speedup of the KRLNS equations is investigated. It is found that all three methods can capture the transient flow fields of unsteady incompressible flow and a large speedup for the KRLNS equations is obtained.

[1]  Nobuyuki Satofuka,et al.  Comparison of Local Computational Approaches for Unsteady Viscous Incompressible Flows , 2014, Modeling, Simulation and Optimization for Science and Technology.

[2]  Nobuyuki Satofuka,et al.  Simulation of doubly periodic shear layers using Kinetically Reduced Local Navier–Stokes equations on a GPU , 2013 .

[3]  Pietro Asinari,et al.  Artificial compressibility method revisited: Asymptotic numerical method for incompressible Navier-Stokes equations , 2010, J. Comput. Phys..

[4]  Effect of collision and velocity model of lattice Boltzmann model on three-dimensional turbulent flow simulation , 2017 .

[5]  Nobuyuki Satofuka,et al.  Higher order numerical simulation of unsteady viscous incompressible flows using kinetically reduced local Navier–Stokes equations on a GPU , 2015 .

[6]  Kazuhiro Nakahashi,et al.  Parallel Processing of the Building-Cube Method on a GPU Platform , 2011 .

[7]  Max Grossman,et al.  Professional CUDA C Programming , 2014 .

[8]  A. Chorin.  A Numerical Method for Solving Incompressible Viscous Flow Problems , 1997 .

[9]  Ulrich Rüde,et al.  Performance modeling and analysis of heterogeneous lattice Boltzmann simulations on CPU-GPU clusters , 2015, Parallel Comput..

[10]  M. Minion,et al.  Performance of Under-resolved Two-Dimensional Incompressible Flow , 1995 .

[11]  Jonathan R Clausen,et al.  Entropically damped form of artificial compressibility for explicit simulation of incompressible flow. , 2012, Physical review. E, Statistical, nonlinear, and soft matter physics.

[12]  D. Wolf-Gladrow.  Lattice-Gas Cellular Automata and Lattice Boltzmann Models: An Introduction , 2000 .

[13]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[14]  Shigeo Kida,et al.  Kolmogorov similarity in freely decaying turbulence , 1987 .

[15]  Santosh Ansumali,et al.  Kinetically reduced local Navier-Stokes equations: an alternative approach to hydrodynamics. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[16]  S. Ansumali,et al.  Kinetically reduced local Navier-Stokes equations for simulation of incompressible viscous flows. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.