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 2

1 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)
 

The Global Science Journal