×

A detailed verification procedure for compressible reactive multicomponent Navier-Stokes solvers. (English) Zbl 1391.76234

Summary: The present manuscript reports a numerical verification study based on a series of tests that allows to evaluate the numerical performance of a compressible reactive multicomponent Navier-Stokes solver. The verification procedure is applied to a density-based finite difference numerical scheme suited to compressible reactive flows representative of either combustion in high speed flows or detonation. The numerical algorithm is based on a third-order accurate total variation diminishing (TVD) Runge Kutta time integration scheme. It employs a seventh-order accurate weighted essentially non-oscillatory (WENO) scheme to discretize the non-linear advective terms while an eighth-order accurate centered finite difference scheme is retained for the molecular viscous and diffusive terms. These molecular contributions are evaluated with the library EGLIB that accounts for detailed multicomponent transport including Soret and Dufour effects. The developed numerical solver thus offers an interesting combination of existing methods suited to the present purpose of studying combustion in high speed flows and/or detonations. The numerical solver is verified by considering a complete procedure that gathers eight elementary verification subsets including, among others, the classical Sod’s shock tube problem, the ignition sequence of a multi-species mixture in a shock tube, the unsteady diffusion of a smoothed concentration profile and a one-dimensional laminar premixed flame. Such verification analyses are seldom reported in the literature but constitute an important part of computational research activities. It is presently completed with the application of the verified finite difference scheme to the numerical simulation of (i) shock (reactive) mixing layer interaction and (ii) combustion ignition downstream of a highly under-expanded jet. The corresponding results shed some light onto the robustness (stability) and performance of the numerical scheme, and also provide some very valuable insights onto the complex physics that prevails in the development of chemical reactions in such situations, which are considered as representative of the discharge or accidental release of high pressure flammable mixtures into the atmosphere.

MSC:

76F65 Direct numerical and large eddy simulation of turbulence
76N15 Gas dynamics (general theory)
76V05 Reaction effects in flows
80A25 Combustion

Software:

PREMIX; EGlib; VODE; CHEMKIN
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Ferri, A., Mixing controlled supersonic combustion, Annu Rev Fluid Mech, 5, 301-338, (1973)
[2] Barrère, M.; Mestré, A., Stabilisation de la flamme en combustion supersonique, La Recherche Aérospatiale, 1, 1-13, (1988)
[3] Billig, F. S., Research on supersonic combustion, J Propul Power, 9, 499-514, (1993)
[4] Im, H. G.; Helenbrook, B. T.; Lee, S. R.; Law, C. K., Ignition in the supersonic hydrogen/air mixing layer with reduced reaction mechanisms, J Fluid Mech, 322, 275-296, (1996) · Zbl 0888.76088
[5] Izard, J. F.; Mura, A., Stabilisation of non-premixed flames in supersonic reactive flows, Comptes Rendus Mecanique, 337, 362-372, (2009) · Zbl 1432.76137
[6] Mura, A.; Izard, J. F., Numerical simulation of supersonic non-premixed turbulent combustion in a scramjet combustor model, J Propul Power, 26, 858-868, (2010)
[7] Birch, A. D.; Brown, D. R.; Dodson, M. G.; Swaffield, F., The structure and concentration decay if high pressure jets of natural gas, Combust Sci Technol, 36, 249-261, (1984)
[8] Birch, A. D.; Brown, D. R.; Cook, D. K.; Hargrave, G. K., Flame stability in under-expanded natural gas jets, Combust Sci Technol, 58, 267-280, (1988)
[9] Ern, A.; Giovangigli, V., Fast and accurate multicomponent transport property evaluation, J Comput Phys, 120, 105-116, (1995) · Zbl 0938.76097
[10] Ern A, Giovangigli V. EGLIB: a general-purpose fortran library for multicomponent transport property evaluation. CERMICS Internal Report, 96-51; 1996.
[11] Ern, A.; Giovangigli, V., Multicomponent transport algorithms, Lecture notes in physics monographs, vol. M24, (1994), Springer-Verlag Heidelberg · Zbl 0820.76002
[12] Hirschfelder, J. O.; Curtiss, C. F., Molecular theory of gases and liquids, (1969), John Wiley & Sons New York
[13] Kee RJ, Rupley FM, Meeks E. CHEMKIN-III: a fortran chemical kinetics package for the analysis of gas-phase chemical and plasma kinetics. Sandia National Laboratories; 1996.
[14] Shahab, M. F.; Lehnasch, G.; Gatski, T. B.; Comte, P., Statistical evaluation of the shock wave boundary layer interaction phenomenon, Flow Turbul Combust, 86, 369-397, (2011) · Zbl 1432.76127
[15] Balsara, D. S.; Shu, C. W., Monotonicity preserving weighted essentially non-oscillatory schemes with increasingly high-order of accuracy, J Comput Phys, 106, 405-452, (2000) · Zbl 0961.65078
[16] Shu CW. Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws. ICASE Report, 97-65; 1997.
[17] Jiang, G. S.; Shu, C. W., Efficient implementation of weighted ENO schemes, J Comput Phys, 126, 202-228, (1996) · Zbl 0877.65065
[18] Strikwerda, J. C., Finite difference schemes and partial differential equations, (1989), Wadsworth Belmont · Zbl 0681.65064
[19] Fedkiw, R. P.; Merriman, B.; Osher, S., High accuracy numerical methods for thermally perfect gas flows with chemistry, J Comput Phys, 132, 175-190, (1997) · Zbl 0888.76053
[20] Gottlieb, S.; Shu, C. W., Total variation diminishing Runge-Kutta schemes, Math Comput, 67, 73-85, (1998) · Zbl 0897.65058
[21] Brown, P. N.; Byrne, G. D.; Hindmarsh, A. C., VODE: a variable coefficient ODE solver, J Sci Stat Comput, 10, 1038-1051, (1989) · Zbl 0677.65075
[22] Sod, G. A., Survey of several finite difference methods for systems of non-linear hyperbolic conservation laws, J Comput Phys, 27, 1-31, (1978) · Zbl 0387.76063
[23] Liska, R.; Wendroff, B., Comparison of several difference schemes on 1D and 2D test problems for the Euler equations, J Sci Comput, 25, 995-1017, (2003) · Zbl 1096.65089
[24] Harten, A., ENO schemes with subcell resolution, J Comput Phys, 83, 148-184, (1989) · Zbl 0696.65078
[25] Henrick, A. K.; Aslam, T. D.; Powers, J. M., Mapped weighted essentially non-oscillatory schemes: achieving optimal order near critical points, J Comput Phys, 207, 542-567, (2005) · Zbl 1072.65114
[26] Harten, A., The artificial compression method for computation of shocks and contact discontinuities. I. single conservation laws, Commun Pure Appl Math, 30, 611-638, (1977) · Zbl 0343.76023
[27] Shen Y, Zha G. Improved seventh-order WENO scheme. In: AIAA paper 2010-1451; 2010.
[28] Borges, R.; Carmona, M.; Costa, B.; Don, W. S., An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws, J Comput Phys, 227, 3191-3211, (2008) · Zbl 1136.65076
[29] Salari K, Knupp P. Code verification by the method of manufactured solutions. Technical report SAND2000-1444; 2000.
[30] Vedovoto, J. M.; Silveira Neto, A.; Mura, A.; Figueira da Silva, L. F., Application of the method of manufactured solutions to the verification of a pressure-based finite-volume numerical scheme, Comput Fluids, 51, 89-99, (2011) · Zbl 1271.76199
[31] Davoudzadeh, F.; McDonald, H.; Thompson, B. E., Accuracy evaluation of unsteady CFD numerical schemes by vortex preservation, Comput Fluids, 24, 883-895, (1995) · Zbl 0873.76050
[32] Yee, H. C.; Vinokur, M.; Djomehri, M. J., Entropy splitting and numerical dissipation, J Comput Phys, 162, 33-81, (2000) · Zbl 0987.65086
[33] Vicquelin R. PhD thesis, Ecole Centrale Paris; 2010.
[34] Vicquelin, R.; Fiorina, B.; Payet, S.; Darabiha, N.; Gicquel, O., Coupling tabulated chemistry with compressible CFD solvers, Proc Combust Inst, 33, 1481-1488, (2011)
[35] Goodwin D. Cantera: object-oriented software for reacting flows. <http://cantera.github.com/docs/sphinx/html/index.html>.
[36] Oran, E.; Young, T.; Boris, J. P.; Cohen, A., Weak and strong ignition. I. numerical simulations of shock tube experiments, Combust Flame, 48, 135-148, (1982)
[37] Deiterding, R., A parallel adaptive method for simulating shock-induced combustion with detailed chemical kinetics in complex domains, Comput Struct, 87, 769-783, (2009)
[38] Ziegler, J. L.; Deiterding, R.; Shepherd, J. E.; Pullin, D. I., An adaptive high-order hybrid scheme for compressive, viscous flows with detailed chemistry, J Comput Phys, 230, 7598-7630, (2011) · Zbl 1433.76144
[39] Kee RJ, Grcar JF, Smooke MD, Miller JA. A fortran program for modelling steady laminar one dimensional premixed flames. Sandia National Laboratories; 1988.
[40] Bray, K. N.C.; Libby, P. A.; Moss, J. B., Flamelet crossing frequencies and mean reaction rates in premixed turbulent combustion, Combust Sci Technol, 41, 143-172, (1984)
[41] Billet, G.; Giovangigli, V.; de Gassowski, G., Impact of volume viscosity on a shock/hydrogen bubble interaction, Combust Theory Modell, 12, 221-248, (2008) · Zbl 1148.80374
[42] Chaudhuri, A.; Hadjadj, A.; Chinnayya, A.; Palerm, S., Numerical study of compressible mixing layers using high-order WENO schemes, J Sci Comput, 47, 170-197, (2011) · Zbl 1213.76120
[43] Miller, M. F.; Bowman, C. T.; Mungal, M. G., An experimental investigation of the effects of compressibility on a turbulent reacting mixing layer, J Fluid Mech, 356, 25-64, (1998)
[44] Cheng, T. S.; Wehrmeyer, J. A.; Pitz, R. W., Raman measurement of mixing and finite-rate chemistry in a supersonic hydrogen-air diffusion flame, Combust Flame, 99, 157-173, (1994)
[45] Papamoschou, D.; Roshko, A., The compressible turbulent shear layer: an experimental study, J Fluid Mech, 197, 453-477, (1988)
[46] Marinov NM, Westbrook CK, Pitz WJ. Detailed and global chemical kinetics model for hydrogen. In: 8th International symposium on transport properties, San Francisco, CA; 1995.
[47] Ó Conaire, M.; Curran, H. J.; Simmie, J. M.; Pitz, W. J.; Westbrook, C. K., A comprehensive modeling study of hydrogen oxidation, Int J Chem Kinet, 36, 603-622, (2004)
[48] Mitchell RE, Kee RJ. A general-purpose computer code for predicting chemical kinetic behavior behind incident and reflected shocks. Sandia National Laboratories; 1989.
[49] Lehnasch, G.; Bruel, P., A robust methodology for RANS simulations of highly under-expanded jets, Int J Numer Methods Fluids, 56, 2079-2205, (2008) · Zbl 1388.76141
[50] Sheeran, W. J.; Dosanjh, D. S., Observations on jet flows from a two-dimensional, under-expanded, sonic nozzle, AIAA J, 6, 540-542, (1968)
[51] Werle, M. J.; Schaffer, D. G.; Driftmyer, R. T., Freejet terminal shocks, AIAA J, 8, 2295-2297, (1970)
[52] Vinson PW. Prediction of reaction control effectiveness at supersonic and hypersonic speeds. OR 6487, Research Division of Martin-Orlando; 1965.
[53] Roe, P. L., Approximate Riemann solvers, parameter vectors, and difference schemes, J Comput Phys, 43, 357-372, (1981) · Zbl 0474.65066
[54] Abgrall, R., Généralization du schéma de roe pour le calcul d’écoulements de mélanges de gaz à concentrations variables, La Recherche Aérospatiale, 3, 31-43, (1988) · Zbl 0662.76097
[55] Grossmann, B.; Cinella, P., Flux-split algorithms for flows with non-equilibrium chemistry and vibrational relaxation, J Comput Phys, 88, 131-168, (1990) · Zbl 0703.76047
[56] Glaister, P., An approximate linearised Riemann solver for the Euler equations of gas dynamics, J Comput Phys, 74, 382-408, (1988) · Zbl 0632.76079
[57] Montagné, J. L.; Yee, H. C.; Vinokur, M., Comparative study of high-resolution shock-capturing schemes for a real gas, AIAA J, 27, 1332-1346, (1989), see also report NASA TM 100004; 1987
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.