A Jacobian-Free Newton Krylov Implicit-Explicit Time Integration Method for Incompressible Flow Problems

We have introduced a fully second order IMplicit/EXplicit (IMEX) time integration technique for solving the compressible Euler equations plus nonlinear heat conduction problems (also known as the radiation hydrodynamics problems) in Kadioglu et al., J. Comp. Physics [22, 24]. In this paper, we study the implications when this method is applied to the incompressible Navier-Stokes (N-S) equations. The IMEX method is applied to the incompressible flow equations in the following manner. The hyperbolic terms of the flow equations are solved explicitly exploiting the well understood explicit schemes. On the other hand, an implicit strategy is employed for the non-hyperbolic terms. The explicit part is embedded in the implicit step in such a way that it is solved as part of the non-linear function evaluation within the framework of the Jacobian-Free Newton Krylov (JFNK) method [8, 29, 31]. This is done to obtain a self-consistent implementation of the IMEX method that eliminates the potential order reduction in time accuracy due to the specific operator separation. We employ a simple yet quite effective fractional step projection methodology (similar to those in [11,19,21,30]) as our preconditioner inside the JFNK solver. We present results from several test calculations. For each test, we show second order time convergence. Finally, we present a study for the algorithm performance of the JFNK solver with the new projection method based preconditioner. AMS subject classifications: 01-08, 35Q35, 76Dxx

[1]  Robert L. Street,et al.  A comparison of staggered and non-staggered grid Navier--Stokes solutions for the 8:1 cavity natural convection flow , 2005 .

[2]  S. Armfield.  Finite difference solutions of the Navier-Stokes equations on staggered and non-staggered grids , 1991 .

[3]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[4]  Yong Jung Kim.  A MATHEMATICAL INTRODUCTION TO FLUID MECHANICS , 2008 .

[5]  Steven J. Ruuth,et al.  Implicit-Explicit Methods for Time-Dependent PDE''s , 1993 .

[6]  Dana A. Knoll,et al.  A fully second order implicit/explicit time integration technique for hydrodynamics plus nonlinear heat conduction problems , 2010, J. Comput. Phys..

[7]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .

[8]  Jie Shen,et al.  An overview of projection methods for incompressible flows , 2006 .

[9]  Steven J. Ruuth,et al.  Implicit-explicit Runge-Kutta methods for time-dependent partial differential equations , 1997 .

[10]  Chi-Wang Shu,et al.  Strong Stability-Preserving High-Order Time Discretization Methods , 2001, SIAM Rev..

[11]  H. B. Keller,et al.  Driven cavity flows by efficient numerical techniques , 1983 .

[12]  A. Chorin.  Numerical Solution of the Navier-Stokes Equations* , 1989 .

[13]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[14]  J. W. Thomas.  Numerical Partial Differential Equations: Finite Difference Methods , 1995 .

[15]  Dana A. Knoll,et al.  Multiphysics Analysis of Spherical Fast Burst Reactors , 2009 .

[16]  Jim E. Morel,et al.  The Coupling of Radiation and Hydrodynamics , 1999 .

[17]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

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

[19]  M. Sussman,et al.  A Coupled Level Set and Volume-of-Fluid Method for Computing 3D and Axisymmetric Incompressible Two-Phase Flows , 2000 .

[20]  Stanley Osher,et al.  A second order primitive preconditioner for solving all speed multi-phase flows , 2005 .

[21]  Boyce E. Griffith,et al.  An accurate and efficient method for the incompressible Navier-Stokes equations using the projection method as a preconditioner , 2009, J. Comput. Phys..

[22]  J. W. Thomas.  Numerical Partial Differential Equations , 1999 .

[23]  B. P. Leonard,et al.  The ULTIMATE conservative difference scheme applied to unsteady one-dimensional advection , 1991 .

[24]  Chi-Wang Shu,et al.  High Order Weighted Essentially Nonoscillatory Schemes for Convection Dominated Problems , 2009, SIAM Rev..

[25]  Rick M. Rauenzahn,et al.  A second order self-consistent IMEX method for radiation hydrodynamics , 2010, J. Comput. Phys..

[26]  William J. Rider,et al.  On consistent time-integration methods for radiation hydrodynamics in the equilibrium diffusion limit: low-energy-density regime , 2001 .

[27]  M. Sussman,et al.  A Second Order JFNK-based IMEX Method for Single and Multi-phase Flows , 2011 .

[28]  Yousef Saad,et al.  Hybrid Krylov Methods for Nonlinear Systems of Equations , 1990, SIAM J. Sci. Comput..

[29]  Steven J. Ruuth.  Implicit-explicit methods for reaction-diffusion problems in pattern formation , 1995 .

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

[31]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[32]  C. Bruneau,et al.  The 2D lid-driven cavity problem revisited , 2006 .

[33]  M. Sussman.  A second order coupled level set and volume-of-fluid method for computing growth and collapse of vapor bubbles , 2003 .

[34]  J. Strikwerda.  Finite Difference Schemes and Partial Differential Equations , 1989 .

[35]  G. Batchelor,et al.  An Introduction to Fluid Dynamics , 1968 .

[36]  D. Drikakis.  Bifurcation phenomena in incompressible sudden expansion flows , 1997 .

[37]  C. Fletcher.  Computational techniques for fluid dynamics , 1992 .

[38]  H. Kreiss,et al.  Time-Dependent Problems and Difference Methods , 1996 .

[39]  Chi-Wang Shu,et al.  Total variation diminishing Runge-Kutta schemes , 1998, Math. Comput..

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

[41]  W. Rider,et al.  High-Resolution Methods for Incompressible and Low-Speed Flows , 2004 .

[42]  P. Colella,et al.  A second-order projection method for the incompressible navier-stokes equations , 1989 .

[43]  Chia-Jung Hsu.  Numerical Heat Transfer and Fluid Flow , 1981 .

[44]  C. Kelley.  Solving Nonlinear Equations with Newton's Method , 1987 .

[45]  P. Moin,et al.  Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations , 1984 .

[46]  P. Wesseling.  Principles of Computational Fluid Dynamics , 2000 .