A parallel Newton-Krylov-Schur flow solver for the Navier-Stokes equations using the SBP-SAT approach

This paper presents a three-dimensional Newton-Krylov flow solver for the NavierStokes equations which uses summation-by-parts (SBP) operators on multi-block structured grids. Simultaneous approximation terms (SAT’s) are used to enforce the boundary conditions and the coupling of block interfaces. The discrete equations are solved iteratively with an inexact Newton method. The linear system of each Newton iteration is solved using a Krylov subspace iterative method with an approximate-Schur parallel preconditioner. The algorithm is validated against an established two-dimensional flow solver. Additionally, results are presented for laminar flow around the ONERA M6 wing, as well as low Reynolds number flow around a sphere. Using 384 processors, the solver is capable of obtaining the steady-state solution (reducing the flow residual by 12 orders of magnitude) on a 4.1 million node grid around the ONERA M6 wing in 4.2 minutes. Convergence to 3 significant figures in force coefficients is achieved in 83 seconds. Parallel scaling tests show that the algorithm scales well with the number of processors used. The results show that the SBP-SAT discretization, solved with the parallel Newton-Krylov-Schur algorithm, is an efficient option for three-dimensional Navier-Stokes solutions, with the SAT’s providing several advantages in enforcing boundary conditions and block coupling.

[1]  Jan Nordström,et al.  Weak and strong wall boundary procedures and convergence to steady-state of the Navier-Stokes equations , 2012, J. Comput. Phys..

[2]  D. Gottlieb,et al.  A Stable and Conservative Interface Treatment of Arbitrary Spatial Accuracy , 1999 .

[3]  Jan Nordström,et al.  High-order finite difference methods, multidimensional linear problems, and curvilinear coordinates , 2001 .

[4]  Jing Gong,et al.  A stable and conservative high order multi-block method for the compressible Navier-Stokes equations , 2009, J. Comput. Phys..

[5]  B. Strand.  Summation by parts for finite difference approximations for d/dx , 1994 .

[6]  F. White.  Viscous Fluid Flow , 1974 .

[7]  Jason E. Hicken,et al.  Parallel Newton-Krylov Solver for the Euler Equations Discretized Using Simultaneous-Approximation Terms , 2008 .

[8]  Ch. Hirsch,et al.  Fundamentals Of Computational Fluid Dynamics , 2016 .

[9]  David W. Zingg,et al.  Efficient Newton-Krylov Solver for Aerodynamic Computations , 1998 .

[10]  D. Gottlieb,et al.  Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: methodology and application to high-order compact schemes , 1994 .

[11]  David W. Zingg,et al.  Higher-order spatial discretization for turbulent aerodynamic computations , 2001 .

[12]  Jan Nordström,et al.  Boundary and Interface Conditions for High-Order Finite-Difference Methods Applied to the Euler and Navier-Stokes Equations , 1999 .

[13]  William W. Willmarth,et al.  Some experimental results on sphere and disk drag , 1971 .

[14]  T. Pulliam.  Efficient solution methods for the Navier-Stokes equations , 1986 .

[15]  David W. Zingg,et al.  A High-Order Parallel Newton-Krylov Flow Solver for the Euler Equations , 2009 .

[16]  Paul Orkwis,et al.  Jacobian update strategies for quadratic and near-quadratic convergence of Newton and Newton-like implicit schemes , 1993 .

[17]  P. Roache.  Fundamentals of computational fluid dynamics , 1998 .

[18]  S. Taneda.  Experimental Investigation of the Wake behind a Sphere at Low Reynolds Numbers , 1956 .

[19]  H. Kreiss,et al.  Finite Element and Finite Difference Methods for Hyperbolic Partial Differential Equations , 1974 .

[20]  David E. Keyes,et al.  Aerodynamic applications of Newton- Krylov-Schwarz solvers , 1995 .

[21]  T. Papanastasiou,et al.  Viscous Fluid Flow , 1999 .

[22]  A. Jameson,et al.  Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes , 1981 .

[23]  Magnus Svärd,et al.  Steady-State Computations Using Summation-by-Parts Operators , 2005, J. Sci. Comput..

[24]  Magnus Svärd,et al.  Stable and accurate schemes for the compressible Navier-Stokes equations , 2008, J. Comput. Phys..

[25]  Jan S. Hesthaven,et al.  A Stable Penalty Method for the Compressible Navier-Stokes Equations: III. Multidimensional Domain Decomposition Schemes , 1998, SIAM J. Sci. Comput..

[26]  Masha Sosonkina,et al.  Distributed Schur Complement Techniques for General Sparse Linear Systems , 1999, SIAM J. Sci. Comput..

[27]  Magnus Svärd,et al.  A stable high-order finite difference scheme for the compressible Navier-Stokes equations: No-slip wall boundary conditions , 2008, J. Comput. Phys..

[28]  B. V. Leer,et al.  Experiments with implicit upwind methods for the Euler equations , 1985 .

[29]  Magnus Svärd,et al.  A stable high-order finite difference scheme for the compressible Navier-Stokes equations, far-field boundary conditions , 2007, J. Comput. Phys..

[30]  David W. Zingg,et al.  A Three-Dimensional Multi-Block Newton-Krylov Flow Solver for the Euler Equations , 2005 .

[31]  Thomas H. Pulliam,et al.  Comparison of Several Spatial Discretizations for the Navier-Stokes Equations , 1999 .