An implicit compact scheme solver for two-dimensional multicomponent flows

Abstract A 2D implicit compact scheme solver has been implemented for the vorticity–velocity formulation in the case of nonreacting, multicomponent, axisymmetric, low Mach number flows. To stabilize the discrete boundary value problem, two sets of boundary closures are introduced to couple the velocity and vorticity fields. A Newton solver is used for solving steady-state and time-dependent equations. In this solver, the Jacobian matrix is formulated and stored in component form. To solve the system of linearized equations within each iteration of Newton’s method, preconditioned Bi-CGSTAB is used in combination with a matrix–vector product computed in component form. The almost dense Jacobian matrix is approximated by a partial Jacobian. For the preconditioner equation, the partial Jacobian is approximately factored using several methods. In a detailed study of several preconditioning techniques, a promising method based on ILUT preconditioning in combination with reordering and double scaling using the MC64 algorithm by Duff and Koster is selected. To validate the implicit compact scheme solver, several nonreacting model problems have been considered. At least third order accuracy in space is recovered on nonuniform grids. A comparison of the results of the implicit compact scheme solver with the results of a traditional implicit low order solver shows an order of magnitude reduction of computer memory and time using the compact scheme solver in the case of time-dependent mixing problems.

[1]  Edmond Chow,et al.  An object-oriented framework for block preconditioning , 1998, TOMS.

[2]  M. Smooke,et al.  Error estimate for the modified Newton method with applications to the solution of nonlinear, two-point boundary-value problems , 1983 .

[3]  Iain S. Duff,et al.  On Algorithms For Permuting Large Entries to the Diagonal of a Sparse Matrix , 2000, SIAM J. Matrix Anal. Appl..

[4]  W.Roger Briley,et al.  An overview and generalization of implicit Navier–Stokes algorithms and approximate factorization , 2001 .

[5]  N. Darabiha,et al.  Vector Computers and Complex Chemistry Combustion , 1988 .

[6]  E. Eckert,et al.  Analysis of heat and mass transfer , 1971 .

[7]  F. DUCROSb,et al.  COMPACT FINITE DIFFERENCE SCHEMES ON NON-UNIFORM MESHES . APPLICATION TO DIRECT NUMERICAL SIMULATIONS OF COMPRESSIBLE FLOWS , 1998 .

[8]  W. Schneider.  Flow induced by jets and plumes , 1981, Journal of Fluid Mechanics.

[9]  S. Tsynkov.  Numerical solution of problems on unbounded domains. a review , 1998 .

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

[11]  Iterative modified approximate factorization , 2001 .

[12]  Youcef Saad,et al.  A Basic Tool Kit for Sparse Matrix Computations , 1990 .

[13]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[14]  J. Spurk.  Boundary Layer Theory , 2019, Fluid Mechanics.

[15]  C. Davies,et al.  A novel velocity-vorticity formulation of the Navier-Stokes equations with applications to boundary layer disturbance evolution , 2001 .

[16]  Alexandre Ern.  Vorticity-velocity modeling of chemically reacting flows , 1994 .

[17]  Robert J. Kee,et al.  A FORTRAN COMPUTER CODE PACKAGE FOR THE EVALUATION OF GAS-PHASE, MULTICOMPONENT TRANSPORT PROPERTIES , 1986 .

[18]  Mitchell D. Smooke,et al.  An implicit compact scheme solver with application to chemically reacting flows , 2005 .

[19]  Alexandre Ern,et al.  Vorticity-Velocity Formulation for Three-Dimensional Steady Compressible Flows , 1993 .

[20]  A. J. Baker,et al.  A 3D incompressible Navier–Stokes velocity–vorticity weak form finite element algorithm , 2002 .

[21]  Amable Liñán Martínez,et al.  The virtual origin as a first-order correction for the far-field description of laminar jets , 2002 .

[22]  M. Benzi.  Preconditioning techniques for large linear systems: a survey , 2002 .

[23]  Yousef Saad,et al.  A Flexible Inner-Outer Preconditioned GMRES Algorithm , 1993, SIAM J. Sci. Comput..

[24]  Ayodeji O. Demuren,et al.  HIGHER-ORDER COMPACT SCHEMES FOR NUMERICAL SIMULATION OF INCOMPRESSIBLE FLOWS, PART I: THEORETICAL DEVELOPMENT , 2001 .

[25]  T. Poinsot.  Boundary conditions for direct simulations of compressible viscous flows , 1992 .

[26]  Thomas B. Gatski,et al.  Review of incompressible fluid flow computations using the vorticity-velocity formulation , 1991 .

[27]  新岡 嵩,et al.  Combustion Institute , 2005, Nature.

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

[29]  W. Schneider.  Decay of momentum flux in submerged jets , 1985, Journal of Fluid Mechanics.

[30]  Alexandre Ern,et al.  Detailed Chemistry Modeling of Laminar Diffusion Flames On Parallel Computers , 1995, Int. J. High Perform. Comput. Appl..

[31]  David E. Keyes,et al.  Towards Polyalgorithmic Linear System Solvers for Nonlinear Elliptic Problems , 1994, SIAM J. Sci. Comput..

[32]  H. L. Langhaar,et al.  Steady Flow in the Transition Length of a Straight Tube , 1942 .

[33]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[34]  W. R. Briley,et al.  Solution of the multidimensional compressible Navier-Stokes equations by a generalized implicit method , 1977 .

[35]  Michele Benzi,et al.  Preconditioning Highly Indefinite and Nonsymmetric Matrices , 2000, SIAM J. Sci. Comput..

[36]  Marc S. Ingber,et al.  A Galerkin implementation of the generalized Helmholtz decomposition of vorticity formulations , 2001 .

[37]  Marcelo H. Kobayashi,et al.  A fourth-order-accurate finite volume compact method for the incompressible Navier-Stokes solutions , 2001 .

[38]  Ayodeji O. Demuren,et al.  Higher-Order Compact Schemes for Numerical Simulation of Incompressible Flows , 1998 .

[39]  H. Fasel,et al.  A Compact-Difference Scheme for the Navier—Stokes Equations in Vorticity—Velocity Formulation , 2000 .