Higher order time integration schemes for the unsteady Navier-Stokes equations on unstructured meshes

The efficiency gains obtained using higher-order implicit Runge-Kutta (RK) schemes as compared with the second-order accurate backward difference schemes for the unsteady Navier-Stokes equations are investigated. Three different algorithms for solving the nonlinear system of equations arising at each time step are presented. The first algorithm (nonlinear multigrid, NMG) is a pseudo-time-stepping scheme which employs a nonlinear full approximation storage (FAS) agglomeration multigrid method to accelerate convergence. The other two algorithms are based on inexact Newton's methods. The linear system arising at each Newton step is solved using iterative/Krylov techniques and left preconditioning is used to accelerate convergence of the linear solvers. One of the methods (LMG) uses Richardson's iterative scheme for solving the linear system at each Newton step while the other (PGMRES) uses the generalized minimal residual method. Results demonstrating the relative superiority of these Newton's method based schemes are presented. Efficiency gains as high as 10 are obtained by combining the higher-order time integration schemes such as fourth-order Runge-Kutta (RK64) with the more efficient inexact Newton's method based schemes (LMG).

[1]  D. P. Young,et al.  GMRES acceleration of computational fluid dynamics codes , 1985 .

[2]  William J. Rider,et al.  Physics-Based Preconditioning and the Newton-Krylov Method for Non-equilibrium Radiation Diffusion , 2000 .

[3]  Dimitri J. Mavriplis,et al.  On Convergence Acceleration Techniques for Unstructured Meshes , 1998 .

[4]  L. Trefethen,et al.  Numerical linear algebra , 1997 .

[5]  Eli Turkel,et al.  PRECONDITIONING-SQUARED METHODS FOR MULTIDIMENSIONAL AERODYNAMICS , 1997 .

[6]  G. Dahlquist.  A special stability problem for linear multistep methods , 1963 .

[7]  Mark D. Sanetrik,et al.  Time-accurate Navier-Stokes calculations with multigrid acceleration , 1993 .

[8]  Chang-Hsien Tai,et al.  Design of optimally smoothing multistage schemes for the Euler equations , 1992 .

[9]  E. Hairer,et al.  Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems , 1993 .

[10]  Charles H. K. Williamson,et al.  A SERIES IN 1/√Re TO REPRESENT THE STROUHAL–REYNOLDS NUMBER RELATIONSHIP OF THE CYLINDER WAKE , 1998 .

[11]  R. Henderson.  Nonlinear dynamics and pattern formation in turbulent wake transition , 1997, Journal of Fluid Mechanics.

[12]  Dimitri J. Mavriplis,et al.  Directional Agglomeration Multigrid Techniques for High-Reynolds Number Viscous Flows , 1998 .

[13]  Hester Bijl,et al.  Time Integration Schemes for the Unsteady Navier-stokes Equations , 2001 .

[14]  A. Jameson.  Time dependent calculations using multigrid, with applications to unsteady flows past airfoils and wings , 1991 .

[15]  D. Mavriplis.  Multigrid Strategies for Viscous Flow Solvers on Anisotropic Unstructured Meshes , 1997 .

[16]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[17]  Dimitri J. Mavriplis,et al.  MULTIGRID TECHNIQUES FOR UNSTRUCTURED MESHES , 1995 .

[18]  Dimitri J. Mavriplis,et al.  Large Eddy and Detached Eddy Simulations Using an Unstructured Multigrid Solver , 2001 .

[19]  R. Dembo,et al.  INEXACT NEWTON METHODS , 1982 .

[20]  Gene H. Golub,et al.  Matrix computations , 1983 .

[21]  M. Carpenter,et al.  Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations , 2003 .

[22]  D. Mavriplis.  An assessment of linear versus non-linear multigrid methods for unstructured mesh solvers , 2001 .

[23]  G. Raithby,et al.  A multigrid method based on the additive correction strategy , 1986 .

[24]  J. W. Ruge,et al.  4. Algebraic Multigrid , 1987 .

[25]  V. Venkatakrishnan,et al.  IMPLICIT SCHEMES AND PARALLEL COMPUTING IN UNSTRUCTURED GRID CFD , 1995 .

[26]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .