A Semi-implicit Numerical Scheme for Reacting Flow: II. Stiff, Operator-Split Formulation

Journal of Computational Physics

Volumn 154, Issue 2, 1999, Pages 428-467

  Knio, Omar M ^a   Najm, Habib N ^b   Wyckoff, Peter S ^b  

^a Johns Hopkins University   (United States)

^b SANDIA NATIONAL LABORATORIES   (United States)

Author keywords

Chemistry; Flow; Implicit; Operator splitting; Projection; Reacting; Stiff

Indexed keywords

COMPUTATIONAL FLUID DYNAMICS; EQUATIONS OF MOTION; GAS DYNAMICS; GASES; HEAT CONVECTION; INTEGRATION; LINEAR EQUATIONS; NONLINEAR EQUATIONS;

FLOW; IMPLICIT; OPERATOR SPLIT; OPERATOR-SPLITTING; PROJECTION; REACTING; REACTING FLOWS; SECOND ORDERS; STIFF; VELOCITY FIELD;

DIFFUSION IN LIQUIDS;

EID: 0000392923     PISSN: 00219991     EISSN: None     Source Type: Journal    
DOI: 10.1006/jcph.1999.6322     Document Type: Article

References (90)

1. R. Aiken (Ed.), Stiff Computation (Oxford Univ. Press, London, 1985).
2. E. Hairer and G. Wanner, Solving Ordinary Differential Equations, Vol. II, Stiff and Differential-Algebraic Problems (Springer-Verlag, Berlin, 1996).
3. M. Frenklach, H. Wang, M. Goldenberg, G. P. Smith, D. M. Golden, C. T. Bowman, R. K. Hanson, W. C. Gardiner, and V. Lissianski, GRImech - An Optimized Detailed Chemical Reaction Mechanism for Methane Combustion, Top. Rep. GRI-95/0058, GRI, 1995.
4. H. N. Najm, P. S. Wyckoff, and O. M. Knio, A semi-implicit numerical scheme for reacting flow. I Stiff chemistry, J. Comput. Phys. 143(2), 381 (1998).
5. M. D. Smooke, I. K. Puri, and K. Seshadri, A comparison between numerical calculations and experimental measurements of the structure of a counterflow diffusion flame burning diluted methane in diluted air, in Proceedings, Twenty-First Symposium (International) on Combustion (Combustion Institute, Pittsburgh, PA, 1986), pp. 1783-1792.
6. C. F. Curtiss and J. O. Hirschfelder, Integration of stiff equations, Proc. Nat. Acad. Sci. U.S.A. 38, 235 (1952).
7. C. W. Gear, Numerical Initial Value Problems in Ordinary Differential Equations (Prentice-Hall, Englewood Cliffs, NJ, 1971).
8. A. C. Hindmarsh, GEAR: Ordinary Differential Equation System Solver, Technical Report UCID-30001, Rev. 3, Lawrence Livermore National Laboratory, Livermore, CA, 1974.
9. A. C. Hindmarsh, GEARB: Solution of Ordinary Differential Equations Having Banded Jacobian, Technical Report UCID-30059, Rev. 2, Lawrence Livermore National Laboratory, Livermore, CA, 1977.
10. A. C. Hindmarsh, LSODE and LSODI: Two new initial value ordinary equation solvers, ACM-SIGNUM 15(4), 10 (1980).
11. A. C. Hindmarsh, ODEPACK: A systematized collection of ODE solvers, in Scientific Computing, edited by R. Stepleman (North-Holland, Amsterdam, 1983), pp. 55-64.
12. P. N. Brown, G. D. Byrne, and A. C. Hindmarsh, VODE: A variable coefficient ODE solver, SIAM J. Sci. Statist. Comput. 10, 1038 (1989).
13. P. J. Vanderhouven and B. P. Sommeijer, Analysis of parallel diagonally implicit iteration of Runge-Kutta methods, Appl. Numer. Math. 11, 169 (1993).
14. N. H. Cong, A parallel DIRK method for stiff initial-value problems, J. Comput. Appl. Math. 54, 121 (1994).
15. X. Zhong, Additive semi-implicit Runge-Kutta methods for computing high-speed nonequilibrium reactive flows, J. Comput. Phys. 128, 19 (1996).
16. J. J. B. Baeza, F. P. Plá, and G. R. Ramos, Stiffness-adaptive Taylor method for the integration of non-stiff and stiff kinetic models, J. Comput. Chem. 13, 810 (1992).
17. P. Sun, D. P. Chock, and S. L. Winkler, An implicit-explicit hybrid solver for a system of stiff kinetic equations, J. Comput. Phys. 115, 515 (1994).
18. D. P. Chock, S. L. Winkler, and P. Sun, Comparison of stiff chemistry solvers for air quality modeling, Environ. Sci. Technol. 28, 1882 (1994).
19. D. Dabdub and J. H. Seinfeld, Extrapolation techniques used in the solution of stiff ODEs associated with chemical kinetics of air quality models, Atmos. Environ. 29, 403 (1995).
20. S. Elliot, R. P. Turco, and M. Z. Jacobson, Tests on combined projection/forward differencing integration for stiff photochemical family systems at long time step, Computers Chem. 17, 91 (1993).
21. E. Hesstvedt, O. Hov, and I. S. A. Isaksen, Quasi-steady-state approximations in air pollution modeling: Comparison of two numerical schemes for oxidant prediction, Int. J. Chem. Kinetics 10, 971 (1978).
22. T. R. Young and J. P. Boris, A numerical technique for solving stiff ordinary differential equations associated with the chemical kinetics of reactive flow problems, J. Phys. Chem. 81, 2424 (1977).
23. J. G. Verwer and M. van Loon, An evaluation of explicit pseudo-steady-state approximation schemes for stiff ODE systems from chemical kinetics, J. Comput. Phys. 113, 347 (1994).
24. W. Gong and H.-R. Cho, A numerical scheme for the integration of the gas-phase chemical rate equations in three-dimensional atmospheric models, Atmos. Environ. 27A, 2591 (1993).
25. R. D. Saylor and G. D. Ford, On the comparison of numerical methods for the integration of kinetic equations in atmospheric chemistry and transport models, Atmos. Environ. 29(19), 2585 (1995).
26. H. H. Robertson, The solution of a set of reaction rate equations, in Numerical Analysis: an Introduction, edited by J. Walsh (Academic Press, New York, 1966), pp. 178-182.
27. K. Radhakrishnan, New integration techniques for chemical kinetic rate equations. I. Efficiency comparison, Combust. Sci. Technol. 46, 59 (1986).
28. Y. D'Angelo and B. Larrouturou, Comparison and analysis of some numerical schemes for stiff complex chemistry problems, RAIRO Modél. Math. Anal. Numér. 29, 259 (1995).
29. J. G. Verwer, Gauss-Seidel iteration for stiff ODES from chemical kinetics, SIAM J. Sci. Comput. 15, 1234 (1994).
30. J. G. Verwer, J. G. Blom, and E. J. Spee, A comparison of stiff ODE solvers for atmospheric chemistry problems, Atmos. Environ. 30, 49 (1996).
31. C. J. Aro, CHEMSODE: A stiff ODE solver for the equations of chemical kinetics, Comput. Phys. Commun. 97, 304 (1996).
32. R. Courant, K. O. Friedrichs, and H. Lewy, Über die Partiellen Differenzengleichungen der Mathematischen Physik, Math. Ann. 100, 32 (1928); Translation, On the partial difference equations of mathematical physics, IBM J. Res. Develop. 11, 215 (1967).
33. R. Courant, K. O. Friedrichs, and H. Lewy, Über die Partiellen Differenzengleichungen der Mathematischen Physik, Math. Ann. 100, 32 (1928); Translation, On the partial difference equations of mathematical physics, IBM J. Res. Develop. 11, 215 (1967).
34. D. A. Anderson, J. C. Tannehill, and R. H. Pletcher, Computational Fluid Mechanics and Heat Transfer (Hemisphere, New York, 1984).
35. C. A. J. Fletcher, Computational Techniques for Fluid Dynamics (Springer-Verlag, Berlin, 1988), Vol. I.
36. J. M. Varah, Stability restrictions on second order, three level finite difference schemes for parabolic equations, SIAM J. Numer. Anal. 17(2), 300 (1980).
37. U. M. Ascher, S. J. Ruuth, and B. T. R. Wetton, Implicit-explicit methods for time-dependent partial differential equations, SIAM J. Numer. Anal. 32(3), 797 (1995).
38. W. Hundsdorfer and J. G. Verwer, A note on splitting errors for advection-reaction equations, Appl. Numer. Math. 18, 191 (1995).
39. J. G. Verwer, J. G. Blom, and W. Hundsdorfer. An implicit-explicit approach for atmospheric transport-chemistry problems, Appl. Numer. Math. 20, 191 (1996).
40. P. Sun, A pseudo-non-time-splitting method in air quality modeling, J. Comput. Phys. 127, 152 (1996).
41. J. G. Verwer, E. J. Spee, J. G. Blom, and W. H. Hundsdorfer, A second order Rosenbrock method applied to photochemical dispersion problems, SIAM J. Sci. Comput. 20, 1456 (1999).
42. J. Frank, W. Hundsdorfer, and J. G. Verwer, On the stability of implicit-explicit linear multistep methods, Appl. Numer. Math. 25, 193 (1997).
43. O. Knoth and R. Wolke, Implicit-explicit Runge-Kutta methods for computing atmospheric reactive flows, Appl. Numer. Math. 28, 327 (1998).
44. J. G. Verwer, and D. Simpson, Explicit methods for stiff ODEs from atmospheric chemistry, Appl. Numer. Math. 18, 413 (1995).
45. O. Knoth and R. Wolke, An explicit-implicit numerical approach for atmospheric chemistry-transport modeling, Atmos. Environ. 32, 1785 (1998).
46. H. N. Najm, O. M. Knio, P. H. Paul, and P. S. Wyckoff, A study of flame observables in premixed methane-air flames, Combust. Sci. Technol. 140(1-6), 369 (1998).
47. W. H. Hundsdorfer, Numerical Solution of Advection-Diffusion-Reaction Equations, Lecture Notes, Report NM-N9603, CWI, Amsterdam, 1996; http://info4u.cwi.nl.
48. E. J. Spee, Coupling advection and chemical kinetics in a global atmospheric test model, in Air Pollution III, edited by H. Power et al. (Computational Mechanics, Southampton-Boston, 1995), Vol. 1, pp. 319-326.
49. J. G. Verwer, J. G. Blom, M. van Loon, and E. J. Spee, A comparison of stiff ODE solvers for atmospheric chemistry problems, Atmos. Environ. 30, 49 (1995).
50. E. J. Spee, P. M. de Zeeuw, J. G. Verwer, J. G. Blom, and W. H. Hundsdorfer, Vectorization and Parallelization of a Numerical Scheme for 3D Global Atmospheric Transport-Chemistry Problems, Report NM-R9620, CWI, Amsterdam, 1996; http://info4u.cwi.nl.
51. E. J. Spee, J. G. Verwer, P. M. de Zeeuw, J. G. Blom, and W. Hundsdorfer, A numerical study for global atmospheric transport-chemistry problems, Math. Comput. Simulation 48, 177 (1998).
52. L. A. Khan and P. L.-F. Liu, An operator splitting algorithm for coupled one-dimensional advection-diffusion-reaction equations, Comput. Methods Appl. Mech. Eng. 127, 181 (1995).
53. Q. Sheng, Solving linear partial differential equations by exponential splitting, IMA J. Numer. Anal. 9, 199 (1989).
54. J. P. Wright, Numerical instability due to varying time steps in explicit wave propagation and mechanics calculations, J. Comput. Phys. 140, 421 (1998).
55. R. J. LeVeque and H. C. Yee, A study of numerical methods for hyperbolic conservation laws with stiff source terms, J. Comput. Phys. 86, 187 (1990).
56. D. Lanser and J. G. Verwer, Analysis of operator splitting for advection-diffusion-reaction problems from air pollution modelling, CWI Report NM-R9805; http://info4u.cwi.nl.
57. G. Strang, On the construction and comparison of difference schemes, SIAM J. Numer. Anal. 5(3), 506 (1968).
58. S. Z. Burstein and A. A. Mirin, Third order difference methods for hyperbolic equations, J. Comput. Phys. 5, 547 (1970).
59. H. Yoshida, Construction of higher order symplectic integrators, Phys. Lett. A 150(5-7), 262 (1990).
60. G. Goyal, P. J. Paul, H. S. Mukunda, and S. M. Deshpande, Time dependent operator-split and unsplit schemes for one dimensional premixed flames, Combust. Sci. Techol. 60, 167 (1988).
61. B. Yang and S. B. Pope, An investigation of the accuracy of manifold methods and splitting schemes in the computational implementation of combustion chemistry, Combust. Flame 112, 16 (1998).
62. A. Majda and J. Sethian, The derivation and numerical solution of the equations for zero Mach number combustion, Combust. Sci. Technol. 42, 185 (1985).
63. H. Schlichting, Boundary-Layer Theory (McGraw-Hill, New York, 1979), 7th ed.
64. F. A. Williams, Combustion Theory (Addison-Wesley, New York, 1985), 2nd ed.
65. R. J. Kee, F. M. Rupley, and J. A. Miller, Chemkin-II: A Fortran Chemical Kinetics Package for the Analysis of Gas Phase Chemical Kinetics, Sandia Report SAND89-8009B, Sandia National Laboratories, Livermore, CA, 1993.
66. A. J. Chorin, A numerical method for solving incompressible viscous flow problems, J. Comput. Phys. 2, 12 (1967).
67. J. B. Bell and D. L. Marcus, A second-order projection method for variable density flows, J. Comput. Phys. 101, 334 (1992).
68. A. S. Almgren, J. B. Bell, and M. L. Welcome, A conservative adaptive projection method for the variable density incompressible Navier-Stokes equations, J. Comput. Phys. 142, 1 (1998).
69. P. A. McMurtry, W.-H. Jou, J. J. Riley, and R. W. Metcalfe, Direct numerical simulations of a reacting mixing layer with chemical heat release, AIAA J. 24(6), 962 (1986).
70. C. Rutland, J. H. Ferziger, and B. J. Cantwell, Effects of Strain, Vorticity, and Turbulence on Premixed Flames, Report TF-44, Thermosciences Division, Mechanical Engineering, Stanford University, Stanford, CA, 1989.
71. C. J. Rutland and J. H. Ferziger, Simulations of flame-vortex interactions, Combust. Flame 84, 343 (1991).
72. S. Mahalingam, B. J. Cantwell, and J. H. Ferziger, Full numerical simulations of coflowing, axisymmetric jet diffusion flames, Phys. Fluids A 2, 720 (1990).
73. J. Kim and P. Moin, Application of a fractional-step method to incompressible Navier-Stokes equations, J. Comput. Phys. 59, 308 (1985).
74. H. N. Najm, A conservative low Mach number projection method for reacting flow modeling, in Transport Phenomena in Combustion, edited by S. Chan (Taylor & Francis, Washington DC, 1996), Vol. 2, pp. 921-932.
75. H. N. Najm and P. S. Wyckoff, Premixed flame response to unsteady strain-rate and curvature, Combust. Flame 110(1-2), 92 (1997).
76. H. N. Najm, O. M. Knio, and P. S. Wyckoff, Response of stoichiometric and rich methane-air flames to unsteady strain-rate and curvature, submitted for publication.
77. R. J. Kee, J. F. Grcar, M. D. Smooke, and J. A. Miller, A Fortran Program for Modeling Steady Laminar One-Dimensional Premixed Flames, Sandia Report SAND85-8240, Sandia National Laboratories, Livemore, CA, 1993.
78. W. T. Ashurst and P. A. McMurtry, Flame generation of vorticity: Vortex dipoles from monopoles, Combust. Sci. Technol. 66, 17 (1989).
79. T. Poinsot, D. Veynante, and S. Candel, Quenching processes and premixed turbulent combustion diagrams, J. Fluid Mech. 228, 561 (1991).
80. T. Mantel, Fundamental Mechanisms in Premixed Flame Propagation via Vortex-Flame Interactions -Numerical Simulations, Annual research briefs, Center for Turbulence Research, Stanford University/NASA Ames Research Center, 1994.
81. M. Hilka, D. Veynante, M. Baum, and T. J. Poinsot, Simulation of flame-vortex interactions using detailed and reduced chemical kinetics, in Proceedings, Tenth Symposium on Turbulent Shear Flows, Penn. State University, University Park, PA, 1995, Vol. 2, Section 19, pp. 19-24.
82. J. Jarosinski, J. H. Lee, and R. Knystautas, Interaction of a vortex ring and a laminar flame, in Proceedings, Twenty-Second Symposium (International) on Combustion (Combustion Institute, 1988), pp. 505-514.
83. W. L. Roberts, J. F. Driscoll, M. C. Drake, and L. P. Goss, Images of the quenching of a flame by a vortex - To quantify regimes of turbulent combustion, Combust. Flame 94, 58 (1993).
84. J.-M. Samaniego, Stretch-Induced Quenching in Flame-Vortex Interactions, Annual research briefs, Center for Turbulence Research, Stanford University/NASA Ames Research Center, 1993.
85. C. J. Mueller, J. F. Driscoll, D. J. Sutkus, W. L. Roberts, M. C. Drake, and M. D. Smooke, Effect of unsteady stretch rate on OH chemistry during a flame-vortex interaction: To assess flamelet models, Combust. Flame 100, 323 (1995).
86. Q.-V. Nguyen and P. H. Paul, The time evolution of a vortex-flame interaction observed via planar imaging of CH and OH, in Proceedings, Twenty-Sixth Symposium (International) on Combustion (Combustion Institute, 1996), pp. 357-364.
87. C. J. Mueller, J. F. Driscoll, D. L. Reuss, M. C. Drake, and M. E. Rosalik, Vorticity generation and attenuation as vortices convect through a premixed flame, Combust. Flame 112(3), 342 (1998).
88. H. N. Najm, P. H. Paul, C. J. Mueller, and P. S. Wyckoff, On the adequacy of certain experimental observables as measurements of flame burning rate, Combust. Flame 113(3), 312 (1998).
89. P. H. Paul and H. N. Najm, Planar laser-induced fluorescence imaging of flame heat release rate, in proceedings, 27th Symposium (International) on Combustion (Combustion Institute, 1998), pp. 43-50.
90. J. G. Verwer, Explicit Runge-Kutta methodsfor parabolic partial differential equations, Appl. Numer. Math. 22, 359 (1996).