zbMATH — the first resource for mathematics

Numerical simulation of two-dimensional combustion using mesh-free methods. (English) Zbl 1244.76075
Summary: The purpose of this research was to develop tools for numerical simulations of flame propagation with mesh-free radial basis functions (RBFs). Mesh-free methods offer many distinct advantages over traditional finite difference, finite element, and finite volume methods. Traditional Lagrangian methods with significant swirl require mesh stiffeners and periodic remeshing to avoid excessive mesh distortion; such codes often require user interaction to repair the meshes before the simulation can proceed again.
A propagating flame of infinite extent is simulated as a collection of normalized cells with periodic boundary conditions. Rather than capturing the flame front, it is tracked as a discontinuity. The flame front is approximated as a product of a Heaviside function in the normal propagation direction and a piece-wise continuous function represented by RBFs in the tangential direction. The cells are subdivided into the burned and unburned sub-domains approximated by two-dimensional periodic RBFs that are constrained to be strictly conservative. The underlying steady flow is vortical with an input turbulent intensity. The governing equations are rotationally and translationally transformed to produce exact differentials that are integrated exactly in time.
In the present paper, the previous results of Aldredge who used a finite-difference level-set method were compared. The physical behavior was remarkably similar, whereas the finite-difference level-set method required 14 h of CPU time, the RBF approach required only 120 CPU seconds on a desktop computer for the case with the largest turbulent intensity. Although there are no other papers that tried to duplicate the original results of Aldredge, the results that are reported here are consistent with the physics observed in other experimental and numerical investigations.

76M25 Other numerical methods (fluid mechanics) (MSC2010)
76V05 Reaction effects in flows
80A25 Combustion
Full Text: DOI
[1] Akkerman, V.; Bychkov, V.; Eriksson, L.-E., Combust flame, 151, 452-471, (2007)
[2] Aldredge, R.C., Premixed flame propagation in a high-intensity, large-scale vortical flow, Combust flame, 106, 29-40, (1996)
[3] Aldredge, R.C., An analytical model for flame propagation in low-Mach-number, variable-density flow, Int commun heat mass transfer, 24, 1163-1169, (1997)
[4] Buhmann, M.D.; Micchelli, C.A., Multiquadric interpolation improved, Comput math appl, 24, 12, 21-25, (1992) · Zbl 0764.41001
[5] Cheng, A.H.D.; Golberg, M.A.; Kansa, E.J.T.; Zammito, T., Exponential convergence and h – c multiquadric collocation method for partial differential equations, Numer methods pdes, 19, 571-594, (2003) · Zbl 1031.65121
[6] de Boor, C.; Kreiss, H.O., On the condition of linear systems associated with discretized BVPs of odes, SIAM J numer anal, 23, 936-939, (1986) · Zbl 0613.65085
[7] Emdadi, A.; Kansa, E.J.; Libre, N.A.; Rahimian, M.; Shekarchi, M., Stable PDE solution methods for large multiquadric shape parameters, Comput methods eng sci, 25, 23-41, (2008) · Zbl 1232.65154
[8] Fasshauer, G.E., Meshfree approximation methods with MATLAB, (2007), World Scientific Singapore · Zbl 1123.65001
[9] Fedoseyev, A.I.; Friedman, M.J.; Kansa, E.J., Improved multiquadric method for elliptic partial differential equations via PDE collocation on the boundary, Comput math appl, 43, 3-5, 491-500, (2002) · Zbl 0999.65137
[10] Fornberg, B.; Driscoll, T.A., Interpolation in the limit of increasingly flat radial basis functions, Comput math appl, 43, 3-5, 413-421, (2002) · Zbl 1006.65013
[11] Fornberg, B.; Fornberg, B., Theoretical and computational aspects of multivariate interpolation with increasingly flat radial basis functions, Comp math appl, 49, 1, 103-130, (2005) · Zbl 1074.41012
[12] Fornberg, B.; Wright, G., Stable computation of multiquadric interpolants for all values of the shape parameter, Comput math appl, 48, 5-6, 853-867, (2004) · Zbl 1072.41001
[13] Fornberg, B.; Wright, G.; Larsson, E., Some observations regarding interpolants in the limit of flat radial basis functions, Comput math appl, 47, 37-55, (2004) · Zbl 1048.41017
[14] Franke, R., Scattered data interpolation: tests of some methods, Math comput, 38, 181-200, (1982) · Zbl 0476.65005
[15] Fureby, C., A fractal flame-wrinkling large eddy model for premixed turbulent simulation combustion, Proc combust inst, 30, 593-601, (2005)
[16] Hon, Y.C.; Schaback, R.; Zhou, X., An adaptive greedy algorithm for solving large RBF collocation problems, Numer algor, 24, 13-25, (2003) · Zbl 1019.65093
[17] Huang, C.-S.; Lee, C.-F.; Cheng, A.H.D., Error estimate, optimal shape parameter, and high precision computation of multiquadric collocation method, Eng anal bound elem, 31, 615-623, (2007) · Zbl 1195.65176
[18] Ingber, M.S.; Chen, C.S.; Tanski, J.A., A mesh free approach using radial basis functions and parallel domain decomposition for solving three dimensional diffusion equations, Int J numer methods eng, 60, 2183-2201, (2004) · Zbl 1178.76276
[19] Kagan, L.; Sivashinsky, G.; Makhviladze, G., On flame extinction by a spatially periodic shear flow, Combus theory modeling, 2, 399-409, (1998) · Zbl 0944.76093
[20] Kansa, E.J., Multiquadrics—a scattered data approximation scheme with applications to computational fluid dynamics: I. surface approximations and partial derivative estimates, Comput math appl, 19, 6-8, 127-145, (1990) · Zbl 0692.76003
[21] Kansa, E.J., Multiquadrics—a scattered data approximation scheme with applications to computational fluid dynamics: II. solutions to parabolic, hyperbolic, and elliptic partial differential equations, Comput math appl, 19, 6-8, 147-161, (1990) · Zbl 0850.76048
[22] Kansa, E.J.; Hon, Y.C., Circumventing the ill-conditioning problem with multiquadric radial basis functions: applications to elliptic partial differential equations, Comput math appl, 39, 7/8, 123-137, (2000) · Zbl 0955.65086
[23] Kansa, E.J., Exact explicit time integration of hyperbolic partial differential equations with radial basis functions, Eng anal bound elem, 31, 577-585, (2007) · Zbl 1195.65144
[24] Lee C-F, Ling L, Schaback R. On convergent numerical algorithms for unsymmetric collocation. Adv Comput Math 2008, to appear. · Zbl 1168.65419
[25] Libre NA, Emdadi A, Kansa EJ, Shekarchi M, Rahimian M. A fast adaptive wavelet scheme in RBF collocation for nearly singular potential PDEs. CMES, 2009, to appear. · Zbl 1357.65272
[26] Ling, L.; Opfer, R.; Schaback, R., Results on meshless collocation techniques, Eng anal bound elem, 30, 147-253, (2006) · Zbl 1195.65177
[27] Ling, L.; Schaback, R., Stable and convergent unsymmetric meshless collocation methods, SIAM J numer anal, 46, 3, 1097-1115, (2008) · Zbl 1167.65059
[28] Ling L, Schaback R. An improved subspace selection algorithm for meshless collocation method. Int J Numer Methods Eng, 2009, to appear. · Zbl 1183.65153
[29] Liu, Y.; Lau, K.S.; Chan, C.K.1.; Guo, Y.C.; Lin, W.Y., Structures of scalar transport in 2D transitional jet diffusion flames by LES, Int J heat mass transfer, 46, 3841-3851, (2003)
[30] Madych, W.R.; Nelson, S.A., Multivariate interpolation and conditionally positive definite functions II, Math comput, 54, 211-230, (1990) · Zbl 0859.41004
[31] Madych, W.R., Miscellaneous error bounds for multiquadric and related interpolators, Comput math appl, 24, 12, 121-138, (1992) · Zbl 0766.41003
[32] Markstein, G.H., Nonsteady flame propagation, (1964), Macmillian New York
[33] Marzouk, Y.M.; Ghoniem, A.F.; Najm, H.N., Toward a flame embedding model for turbulent combustion simulation, Aiaa j, 41, 641-652, (2003)
[34] Osher, S.; Sethian, J.A., Fronts propagating with curvature-dependent speed-algorithms based on hamilton – jacobi formulations, J comput phys, 79, 12-49, (1988) · Zbl 0659.65132
[35] Qian, H.; Tryggvason, G.; Law, C.K., A front tracking method for the motion of premixed flames, J comput phys, 144, 52-69, (1998) · Zbl 1392.76050
[36] Schaback R. 1995. Multivariate interpolation and approximation by translates of a basis function. In: Approximation theory VIII, College Station, TX, vol. 1. · Zbl 1139.41301
[37] Schaback, R., On the efficiency of interpolation by radial basis functions, (), 309-318 · Zbl 0937.65013
[38] Vaezi, V.; Aldredge, R.C., Laminar-flame instabilities in a taylor – couette combustor, Combust flame, 121, 356-366, (2000)
[39] Volokh, K.Y.; Vilnay, O., Pin-pointing solution of ill-conditioned square systems of linear equations, Appl math lett, 13, 119-124, (2000) · Zbl 0955.65031
[40] Warnatz, J.; Maas, U.; Dibble, R.W., Combustion: physical and chemical fundamentals, modeling and simulation, experiments, and pollution formation, (2001), Springer Berlin, New York · Zbl 0976.80005
[41] Wicksall, D.M.; Agrawal, A.K.; Schefer, R.W.; Keller, J.O., The interaction of flame and flow field in a Lean premixed swirl-stabilized combustor operated on H-2/CH4/air, Proc combust inst, 30, 2875-2883, (2005)
[42] Yakhot, V., Propagation velocity of premixed turbulent flames, Combust sci technol, 60, 191-214, (1988)
[43] Brown D, Ling L, Kansa EJ, Levesley J. On approximate cardinal preconditioning methods for solving PDEs with radial basis functions. Eng Anal Bound Elem 2005;29:343-53. · Zbl 1182.65174
[44] Driscoll TA, Heryudono ARH. Adaptive residual subsampling for radial basis function interpolation collocations problems. Comput Math Appl 2007;53:927-39.
[45] Fasshauer GE. Solving differential equations with radial basis functions: multilevel methods and smoothing. Adv Comput Math 1999;11:139-59. · Zbl 0940.65122
[46] Hardy RL. Multiquadric equations of topography and other irregular surfaces. J Geophys Rev 1971;76:1905-915.
[47] Hardy RL. Theory and applications of the multiquadric-biharmonic method: 20 years of discovery 1968-1988. Comput Math Appl 1990;19(8/9):163-208. · Zbl 0692.65003
[48] Ling L, Kansa EJ. A least-squares preconditioner for radial basis function collocation methods. Adv Comput Math 2005;23:31-54. · Zbl 1067.65136
[49] Ling L, Hon YC. Improved numerical solver for Kansa’s method based on affine space decomposition. Eng Anal Bound Elem 2005;29:1077-1085. · Zbl 1182.65176
[50] Wertz J, Kansa EJ, Ling L. The role of the multiquadric shape parameters in solving elliptic partial differential equations. Comput Math Appl 2006;51:1335-48. · Zbl 1146.65078
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. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.