Fast unsteady flow computations with a Jacobian-free Newton–Krylov algorithm

Despite the advances in computer power and numerical algorithms over the last decades, solutions to unsteady flow problems remain computing time intensive. Especially for large Reynolds number flows, nonlinear multigrid, which is commonly used to solve the nonlinear systems of equations, converges s...

Full description

Saved in:
Bibliographic Details
Published in:Journal of computational physics 2010-12, Vol.229 (24), p.9201-9215
Main Authors: Lucas, Peter, van Zuijlen, Alexander H., Bijl, Hester
Format: Article
Language:English
Subjects:
Algorithms
Computational fluid dynamics
Computational techniques
Exact sciences and technology
Fluid flow
Higher order implicit time integration
Jacobian-free Newton–Krylov
Mathematical methods in physics
Mathematical models
Nonlinearity
Physics
Reynolds number
Stiffness due to various flow and grid parameters
Turbulence
Turbulent flow
Unsteady flow
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Despite the advances in computer power and numerical algorithms over the last decades, solutions to unsteady flow problems remain computing time intensive. Especially for large Reynolds number flows, nonlinear multigrid, which is commonly used to solve the nonlinear systems of equations, converges slowly. The stiffness induced by the large aspect ratio cells and turbulence is not tackled well by this solution method. In previous work we showed that a Jacobian-free Newton–Krylov ( jfnk) algorithm, preconditioned with an approximate factorization of the Jacobian that approximately matches the target residual operator, enables a speed up of a factor of 10 compared to standard nonlinear multigrid for two-dimensional, large Reynolds number, unsteady flow computations. The goal of this paper is to demonstrate that the jfnk algorithm is also suited to tackle the stiffness induced by the maximum aspect ratio, the grid density, the physical time step and the Reynolds number. Compared to standard nonlinear multigrid, speed ups up to a factor of 25 are achieved.
ISSN:0021-9991
1090-2716
DOI:10.1016/j.jcp.2010.08.033