Sub-iteration leads to accuracy and stability enhancements of semi-implicit schemes for the Navier-Stokes equations

We present an iterative semi-implicit scheme for the incompressible Navier-Stokes equations, which is stable at CFL numbers well above the nominal limit. We have implemented this scheme in conjunction with spectral discretizations, which suffer from serious time step limitations at very high resolution. However, the approach we present is general and can be adopted with finite element and finite difference discretizations as well. Specifically, at each time level, the nonlinear convective term and the pressure boundary condition - both of which are treated explicitly in time - are updated using fixed-point iteration and Aitken relaxation. Eigenvalue analysis shows that this scheme is unconditionally stable for Stokes flows while numerical results suggest that the same is true for steady Navier-Stokes flows as well. This finding is also supported by error analysis that leads to the proper value of the relaxation parameter as a function of the flow parameters. In unsteady flows, second- and third-order temporal accuracy is obtained for the velocity field at CFL number 5-14 using analytical solutions. Systematic accuracy, stability, and cost comparisons are presented against the standard semi-implicit method and a recently proposed fully-implicit scheme that does not require Newton's iterations. In addition to its enhanced accuracy and stability, the proposed method requires the solution of symmetric only linear systems for which very effective preconditioners exist unlike the fully-implicit schemes.

[1]  Steven J. Ruuth,et al.  Implicit-explicit methods for time-dependent partial differential equations , 1995 .

[2]  G. Karniadakis,et al.  Suppressing vortex-induced vibrations via passive means , 2009 .

[3]  S. Orszag,et al.  High-order splitting methods for the incompressible Navier-Stokes equations , 1991 .

[4]  André Robert,et al.  A stable numerical integration scheme for the primitive meteorological equations , 1981 .

[5]  Bruce M. Irons,et al.  A version of the Aitken accelerator for computer iteration , 1969 .

[6]  J. Butcher.  The numerical analysis of ordinary differential equations: Runge-Kutta and general linear methods , 1987 .

[7]  Jie Shen,et al.  Velocity-Correction Projection Methods for Incompressible Flows , 2003, SIAM J. Numer. Anal..

[8]  J. Verwer,et al.  Numerical solution of time-dependent advection-diffusion-reaction equations , 2003 .

[9]  J. Brandts.  [Review of: W. Hundsdorfer, J.G. Verwer (2003) Numerical Solution of Time-Dependent Advection-Diffusion-Reaction Equations] , 2006 .

[10]  S. Dong,et al.  Strong and auxiliary forms of the semi-Lagrangian method for incompressible flows , 2005 .

[11]  M. Zamir,et al.  The Physics of Pulsatile Flow , 2000, Biological Physics Series.

[12]  Fotis Sotiropoulos,et al.  Curvilinear immersed boundary method for simulating fluid structure interaction with complex 3D rigid bodies , 2008, J. Comput. Phys..

[13]  W. Wall,et al.  Fixed-point fluid–structure interaction solvers with dynamic relaxation , 2008 .

[14]  David J. Silvester,et al.  Implicit algorithms and their linearization for the transient incompressible Navier-Stokes equations , 1997 .

[15]  G. Karniadakis,et al.  A direct numerical simulation study of flow past a freely vibrating cable , 1997, Journal of Fluid Mechanics.

[16]  M. A. Casarin,et al.  Low-energy basis preconditioning for elliptic substructured solvers based on unstructured spectral/ hp element discretization , 2001 .

[17]  George Em Karniadakis,et al.  Frequency selection and asymptotic states in laminar wakes , 1989, Journal of Fluid Mechanics.

[18]  R. F. Warming,et al.  An Implicit Factored Scheme for the Compressible Navier-Stokes Equations , 1977 .

[19]  Howard C. Elman,et al.  Fast Nonsymmetric Iterations and Preconditioning for Navier-Stokes Equations , 1996, SIAM J. Sci. Comput..

[20]  Jie Shen,et al.  An unconditionally stable rotational velocity-correction scheme for incompressible flows , 2010, J. Comput. Phys..

[21]  Francis X. Giraldo,et al.  A spectral element semi-Lagrangian (SESL) method for the spherical shallow water equations , 2003 .

[22]  G. Karniadakis,et al.  Spectral/hp Element Methods for CFD , 1999 .

[23]  Spencer J. Sherwin,et al.  Parallel performance of the coarse space linear vertex solver and low energy basis preconditioner for spectral/hp elements , 2009, Parallel Comput..

[24]  Ekkehard Ramm,et al.  Accelerated iterative substructuring schemes for instationary fluid-structure interaction , 2001 .

[25]  M. O. Deville,et al.  Numerical Evaluation of the Accuracy and Stability Properties of High-order Direct Stokes Solvers with or without Temporal Splitting , 2006, J. Sci. Comput..