×

An implicit compact scheme solver for two-dimensional multicomponent flows. (English) Zbl 1177.76263

Summary: 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.

MSC:

76M20 Finite difference methods applied to problems in fluid mechanics
76R99 Diffusion and convection
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Noskov, M.; Smooke, M. D., An implicit compact scheme solver with application to chemically reacting flows, J Comput Phys, 203, 700-730 (2005) · Zbl 1071.80004
[2] Lele, S. K., Compact finite difference schemes with spectral-like resolution, J Comput Phys, 103, 16-42 (1992) · Zbl 0759.65006
[3] Gamet, L.; Ducros, F.; Nicoud, F.; Poinsot, T., Compact finite difference schemes on non-uniform meshes. Application to direct numerical simulations of compressible flows, Int J Numer Meth Fluids, 29, 159-191 (1999) · Zbl 0939.76060
[4] Bennett, B. A.V.; Smooke, M. D., Unsteady axisymmetric laminar diffusion flames: an application of local rectangular refinement, (Proc 2001 technical meeting of the eastern states section of the Combustion Institute, Hilton Head, SC, December 2-5 (2001), The Combustion Institute), 332-335
[5] Bennett, B. A.V.; Smooke, M. D., Local rectangular refinement solution-adaptive gridding with application to unsteady combustion, (Ninth international conference on numerical combustion, Sorrento, Italy (2002), SIAM and The Combustion Institute), 213-214
[6] Bennett, B. A.V.; Smooke, M. D., Local rectangular refinement solution-adaptive gridding with application to oscillating laminar diffusion flames, (Proc 3rd joint meeting of the US sections of the Combustion Institute, Chicago, IL, March 16-19 (2003), The Combustion Institute), B37
[7] Demuren, A. O.; Wilson, R. V.; Carpenter, M., Higher-order compact schemes for numerical simulation of incompressible flows. Part I: Theoretical development, Numer Heat Transfer Part B, 39, 207-230 (2001)
[8] Pereira, J. M.C.; Kobayashi, M. H.; Pereira, J. C.F., A fourth-order-accurate finite volume compact method for the incompressible Navier-Stokes solutions, J Comput Phys, 167, 217-243 (2001) · Zbl 1013.76054
[9] Saad, Y., A flexible inner-outer preconditioned GMRES algorithm, SIAM J Sci Comput, 14, 2, 461-469 (1993) · Zbl 0780.65022
[10] Knoll, D. A.; Keyes, D. E., Jacobian-free Newton-Krylov methods: a survey of approaches and applications, J Comput Phys, 193, 357-397 (2004) · Zbl 1036.65045
[11] Saad, Y., Iterative Methods for Sparse Linear Systems (2003), SIAM: SIAM Philadelphia · Zbl 1002.65042
[12] Noskov M. An implicit compact scheme solver for modeling steady-state and time-dependent axisymmetric flows with application to laminar diffusion flame simulations. PhD thesis, Mechanical Engineering Department, Yale University, New Haven, 2003.; Noskov M. An implicit compact scheme solver for modeling steady-state and time-dependent axisymmetric flows with application to laminar diffusion flame simulations. PhD thesis, Mechanical Engineering Department, Yale University, New Haven, 2003.
[13] Ern, A.; Smooke, M. D., Vorticity-velocity formulation for three-dimensional steady compressible flows, J Comput Phys, 105, 58-71 (1993) · Zbl 0783.76059
[14] Ern A. Vorticity-velocity modeling of chemically reacting flows. PhD thesis, Mechanical Engineering Department, Yale University, New Haven, 1994.; Ern A. Vorticity-velocity modeling of chemically reacting flows. PhD thesis, Mechanical Engineering Department, Yale University, New Haven, 1994.
[15] Ern, A.; Douglas, C. C.; Smooke, M. D., Detailed chemistry modeling of laminar diffusion flames on parallel computers, Int J Supercomput Appl High Perf Comput, 9, 167-186 (1995)
[16] Gatski, T. B., Review of incompressible fluid flow computations using the vorticity-velocity formulation, Appl Numer Math, 7, 227-239 (1991) · Zbl 0714.76033
[17] Meitz, H. L.; Fasel, H. F., A compact-difference scheme for the Navier-Stokes equations in vorticity-velocity formulation, J Comput Phys, 157, 371-403 (2000) · Zbl 0960.76059
[18] Giovangigli, V.; Darabiha, N., Vector computers and complex chemistry combustion, (Brauner, C. M.; Schmidt-Laine, C., Mathematical modeling in combustion and related topics, NATO Adv Sci Inst Ser E (1988), Martinus Nijhoff Pub.: Martinus Nijhoff Pub. Dordrecht), 491-503
[19] Smooke, M. D., Error estimate for the modified Newton method with applications to the solution of non-linear, 2-point boundary-value problems, J Optimiz Theory Appl, 39, 489-511 (1983) · Zbl 0487.65045
[20] Ern, A.; Giovangigli, V.; Keyes, D. E.; Smooke, M. D., Towards polyalgorithmic linear system solvers for nonlinear elliptic problems, SIAM J Sci Comput, 15, 681-703 (1994) · Zbl 0807.65027
[21] Van der Vorst, H. A., Bi-CGSTAB: a fast and smoothly converging variant of Bi-CG for the solution of non-symmetric linear systems, SIAM J Sci Stat Comput, 13, 631-644 (1992) · Zbl 0761.65023
[22] Wolfram Research, Inc. Mathematica. Version 5.0, Champaign, IL, 2003.; Wolfram Research, Inc. Mathematica. Version 5.0, Champaign, IL, 2003.
[23] Kee RJ, Dixon-Lewis G, Warnatz J, Coltrin ME, Miller JA. A Fortran computer code package for the evaluation of gas-phase multicomponent transport properties. Report SAND86-8246, Sandia National Laboratories, 1986.; Kee RJ, Dixon-Lewis G, Warnatz J, Coltrin ME, Miller JA. A Fortran computer code package for the evaluation of gas-phase multicomponent transport properties. Report SAND86-8246, Sandia National Laboratories, 1986.
[24] Wong, K. L.; Baker, A. J., A 3D incompressible Navier-Stokes velocity-vorticity weak form finite element algorithm, Int J Numer Meth Fluids, 38, 99-123 (2002) · Zbl 1009.76059
[25] Ingber, M. S.; Kempka, S. N., A Galerkin implementation of the generalized Helmholtz decomposition for vorticity formulations, J Comput Phys, 169, 215-237 (2001) · Zbl 1037.76035
[26] Davies, C.; Carpenter, P. W., A novel velocity-vorticity formulation of the Navier-Stokes equations with applications to boundary layer disturbance evolution, J Comput Phys, 172, 119-165 (2001) · Zbl 1065.76573
[27] MacCormack, R. W., Iterative modified approximate factorization, Comput Fluids, 30, 917-925 (2001) · Zbl 0984.76056
[28] Hindmarsh AC, Sloan LJ, Dubois PF. DEC/SOL: solution of dense systems of linear algebraic equations. UCID-30137. Technical report, Lawrence Livermore Laboratory, 1978.; Hindmarsh AC, Sloan LJ, Dubois PF. DEC/SOL: solution of dense systems of linear algebraic equations. UCID-30137. Technical report, Lawrence Livermore Laboratory, 1978.
[29] Chow, E.; Heroux, M. A., An object-oriented framework for block preconditioning, ACM Trans Math Software, 24, 159 (1998) · Zbl 0930.65052
[30] Benzi, M.; Haws, J. C.; Tuma, M., Preconditioning highly indefinite and non-symmetric matrices, SIAM J Sci Comput, 22, 1333-1353 (2000) · Zbl 0985.65036
[31] Duff, I. S.; Koster, J., On algorithms for permuting large entries to the diagonal of a sparse matrix, SIAM J Matrix Anal Appl, 22, 973 (2001) · Zbl 0979.05087
[32] Saad Y. SPARSKIT: a basic tool-kit for sparse matrix computations (version 2), 1998. Available from: http://www-users.cs.umn.edu/ saad/software/SPARSKIT/sparskit.html; Saad Y. SPARSKIT: a basic tool-kit for sparse matrix computations (version 2), 1998. Available from: http://www-users.cs.umn.edu/ saad/software/SPARSKIT/sparskit.html
[33] Benzi, M., Preconditioning techniques for large linear systems: a survey, J Comput Phys, 182, 418-477 (2002) · Zbl 1015.65018
[34] Schlichting, H., Boundary layer theory (1979), McGraw-Hill: McGraw-Hill New York
[35] Eckert, E. R.G.; Drake, R. M., Analysis of heat and mass transfer (1972), McGraw-Hill: McGraw-Hill New York · Zbl 0247.76079
[36] Langhaar, H. L., Steady flow in the transition length of a straight tube, Trans ASME, 64, A55-A58 (1942)
[37] Markatou P. A computational study of methane/air and hydrogen/air combustion above a heated catalytic and non-catalytic surface. PhD thesis, Mechanical Engineering Department, Yale University, New Haven, 1991.; Markatou P. A computational study of methane/air and hydrogen/air combustion above a heated catalytic and non-catalytic surface. PhD thesis, Mechanical Engineering Department, Yale University, New Haven, 1991.
[38] Briley, W. R.; McDonald, H., Solution of the multi-dimensional compressible Navier-Stokes equations by a generalized implicit method, J Comput Phys, 24, 372-397 (1977) · Zbl 0363.76018
[39] Briley, W. R.; McDonald, H., An overview and generalization of implicit Navier-Stokes algorithms and approximate factorization, Comput Fluids, 30, 807-828 (2001) · Zbl 1056.76060
[40] Schneider, W., Flow induced by jets and plumes, J Fluid Mech, 108, 55-65 (1981) · Zbl 0475.76043
[41] Schneider, W., Decay of momentum in submerged jets, J Fluid Mech, 154, 91-110 (1985)
[42] Schneider, W.; Zauner, E.; Bohm, H., The recirculatory flow induced by a laminar axisymmetric jet issuing from a wall, J Fluids Eng, 109, 237-241 (1987)
[43] Poinsot, T. J.; Lele, S. K., Boundary conditions for direct simulations of compressible viscous flows, J Comput Phys, 101, 104-129 (1992) · Zbl 0766.76084
[44] Tsynkov, S. V., Numerical solution of problems on unbounded domains. A review, Appl Numer Math, 27, 465-532 (1998) · Zbl 0939.76077
[45] Revuelta, A.; Sanchez, A. L.; Linan, A., The virtual origin as a first-order correction for the far-field description of laminar jets, Phys Fluids, 14, 1821-1824 (2002) · Zbl 1185.76312
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.