Commun. Comput. Phys., 13 (2013), pp. 1408-1431. A Jacobian-Free Newton Krylov Implicit-Explicit Time Integration Method for Incompressible Flow Problems Samet Y. Kadioglu 1*, Dana A. Knoll 21 Department of Mathematical Engineering, Yildiz Technical University, 34210 Davutpasa-Esenler, Istanbul, Turkey; and Fuels Modeling and Simulation Department, Idaho National Laboratory, P.O. Box 1625,MS 3840, Idaho Falls, ID 83415, USA. 2 Theoretical Division, Los Alamos National Laboratory, P.O. Box 1663, MS B-216, Los Alamos, NM, 87545, USA. Received 2 March 2012; Accepted (in revised version) 18 July 2012 Available online 8 October 2012 doi:10.4208/cicp.020312.180712a Abstract 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 Notice: Undefined variable: pac in /var/www/html/issue/abstract/readabs.php on line 164 Key words: Incompressible flow, Navier-Stokes equations, IMEX method, JFNK method, preconditioner. *Corresponding author. Email: samet.kadioglu@inl.gov (S. Y. Kadioglu), nol@lanl.gov (D. A. Knoll)