Adaptive decomposition finite difference methods for solving singular problems—A review

Abstract

Decomposition, or splitting, finite difference methods have been playing an important role in the numerical solution of nonsingular differential equation problems due to their remarkable efficiency, simplicity, and flexibility in computations as compared with their peers. Although the numerical strategy is still in its infancy for solving singular differential equation problems arising from many applications, explorations of the next generation decomposition schemes associated with various kinds of adaptations can be found in many recent publications. The novel approaches have been proven to be highly effective and reliable in operations. In this article, we will focus on some of the latest developments in the area. Key comments and discussion will be devoted to two particularly interesting issues in the research, that is, direct solutions of degenerate singular reaction-diffusion equations and nonlinear sine-Gordon wave equations. Numerical experiments with simulated demonstrations will be given.

Keywords

Singularity / degeneracy / finite difference approximation / uniform and nonuniform grid / decomposition / adaptation / monotonicity and stability / large system of equations

References

[1] Ashyraliyev M., Blom J. G., Verwer J. G. On the numerical solution of diffusionreaction equations with singular source terms. J Comput Appl Math, 2008, 216: 20-38.

[2] Bandle C., Brunner H. Blowup in diffusion equations: a survey. Special Issue on Nonlinear Problems with Blow-Up Solutions: Applications & Numerical Analysis. J Comput Appl Math, 1998, 97: 3-22.

[3] Budd C. J., Huang W., Russell R. D. Moving mesh methods for problems with blowup. SIAM J Sci Comput, 1996, 17: 305-327.

[4] Budd C. J., Koomullil G. P., Stuart A. M. On the solution of convection-diffusion boundary value problems using equidistributed grids. SIAM J Sci Comput, 1998, 20: 591-618.

[5] Cao W., Huang W., Russell R. D. A study of monitor functions for two-dimensional adaptive mesh generation. SIAM J Sci Comput, 1999, 20: 1978-1994.

[6] Chan C. Y. Quenching and blow-up phenomena. Dynam Contin Discrete Impuls Systems (Ser A), 2003, 10: 403-409.

[7] Chan C. Y., Boonklurb R. A blow-up criterion for a degenerate parabolic problem due to a concentrated nonlinear source. Quart Appl Math, 2007, 65: 781-787.

[8] Chan C. Y., Tragoonsirisak P. A multi-dimensional quenching problem due to a concentrated nonlinear source in ℝ^N. Nonlinear Anal, 2008, 69: 1494-1514.

[9] Cheng H., Lin P., Sheng Q., Tan R. C. E. Solving degenerate reaction-diffusion equations via variable step Peaceman-Rachford splitting. SIAM J Sci Comput, 2003, 25: 1273-1292.

[10] Chin S. A. Quantum statistical calculations and symplectic corrector algorithms. Phys Rev E, 2004, 69: 046118/1-7.

[11] Christiansen P. L., Lomdahl P. S. Numerical study of 2+1 dimensional sine-Gordon solitons. Physica D, 1981, 2: 482-494.

[12] Cisneros L. A., Minzoni A. A. Asymptotics for kink propagation in the discrete sine-Gordon equation. Physica D, 2008, 237: 50-65.

[13] Di Y., Li R., Tang T. A general moving mesh framework in 3D and its application for simulating the mixture of multi-phase flows. Commun Comput Phys, 2008, 3: 582-603.

[14] Djidjeli K., Price W. G., Twizell E. H. Numerical solutions of a damped sine-Gordon equation in two space variables. J Engrg Math, 1995, 29: 347-369.

[15] Douglas J. Jr. On the numerical integration of u_xx + u_yy = u_t by implicit methods. J Soc Indust Appl Math, 1955, 3: 42-65.

[16] Douglas J. Jr. A note on the alternating direction implicit method for the numerical solution of heat flow problems. Proc Amer Math Soc, 1957, 8: 409-412.

[17] Faragó I., Geiser J. Iterative operator-splitting methods for linear problems. Inter J Comput Sci Eng, 2007, 3: 255-263.

[18] Geiser J. Iterative operator-splitting methods with higher-order time integration methods and applications for parabolic partial differential equations. J Comput Appl Math, 2008, 217: 227-242.

[19] Hairer E., Lubich C., Wanner G. Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations, 2006, Springer: New York and Berlin.

[20] Hong J., Liu Y., Munthe-Kaas H., Zanna A. Globally conservative properties and error estimation of a multi-symplectic scheme for Schrödinger equations with variable coefficients. Appl Numer Math, 2006, 56: 814-843.

[21] Huang W., Ma J., Russell R. D. A study of moving mesh PDE methods for numerical simulation of blowup in reaction diffusion equations. J Comput Phys, 2008, 227: 6532-6552.

[22] Iserles A, Sheng Q. Implementation of splitting methods. DAMTP Research Report, NA8/1987. Cambridge University, 1987

[23] Iserles A., Zanna A. Efficient computation of the matrix exponential by generalized polar decompositions. SIAM J Numer Anal, 2005, 42: 2218-2256.

[24] Jain B., Sheng A. D. An exploration of the approximation of derivative functions via finite differences. Rose-Hulman Undergrad Math J, 2007, 8: 172-188.

[25] Kawarada H. On solutions of initial-boundary problem for u_t = u_xx + 1/(1 − u). Pul Res Inst Math Sci, 1975, 10: 729-736.

[26] Khaliq A., Abukhodair B., Sheng Q. A predictor-corrector scheme for the sine-Gordon equation. Numer Methods Partial Differential Eqns, 2000, 16: 133-146.

[27] Khaliq A., Sheng Q. On monotonicity of an adaptive splitting scheme for twodimensional singular reaction-diffusion equations. Inter J Computer Math, 2007, 84: 795-806.

[28] Lang J., Verwer J. G. On global error estimation and control for initial value problems. SIAM J Sci Comput, 2007, 29: 1460-1475.

[29] LeVeque R. J. Numerical Methods for Conservation Laws, 2006, Basel: Birkhäuser.

[30] Levine H. A. Quenching, nonquenching and beyond quenching for solution of some parabolic equations. Ann Math Pura Appl, 1989, 155: 243-260.

[31] Liang K., Lin P., Tan R. C. E. Numerical solution of quenching problems using meshdependent variable temporal steps. Appl Numer Math, 2007, 57: 791-800.

[32] Liu F., Ji S., Liao G. An adaptive grid method and its application to steady Euler flow calculations. SIAM J Sci Comput, 1998, 20: 811-825.

[33] McLachlan R. I. Explicit symplectic splitting methods applied to PDEs. Lectures in Appl Math, 1993, 29: 325-337.

[34] McLachlan R. I., Quispel G. Splitting methods. Acta Numerica, 2002, 11: 341-434.

[35] Moler C., Van Loan C. Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later. SIAM Review, 2003, 45: 1-46.

[36] Oliger J., Zhu X. Stability and error estimation for component adaptive grid methods. Appl Numer Math, 1996, 20: 407-426.

[37] Peaceman D. W., Rachford H. H. Jr. The numerical solution of parabolic and elliptic differential equations. J Soc Indust Appl Math, 1955, 3: 28-41.

[38] Ren W., Wang X. -P. An iterative grid redistribution method for singular problems in multiple dimensions. J Comput Phys, 2000, 159: 246-273.

[39] Serbin S. M. Some cosine schemes for second-order systems of DFE’s with timevarying coefficients. SIAM J Sci Stat Comput, 1985, 6: 61-68.

[40] Sheng Q. Solving linear partial differential equations by exponential splitting. IMA J Numer Anal, 1989, 9: 199-212.

[41] Sheng Q. Global error estimate for exponential splitting. IMA J Numer Anal, 1993, 14: 27-56.

[42] Sheng Q. A note on asymptotic splitting and its applications. Math Comput Modelling, 1994, 20: 45-58.

[43] Sheng Q. Hybrid approximations via second order combined dynamic derivatives on time scales. Electr J Qualit Theory Diff Eqns, 2007, 17: 1-13.

[44] Sheng Q. A moving-mesh splitting scheme for 2-dimensional quenching problems. Proc Appl Math Mech, 2008, 7: 1023303-1023304.

[45] Sheng Q., Cheng H. An adaptive grid method for degenerate semilinear quenching problems. Computers Math Appl, 2000, 39: 57-71.

[46] Sheng Q., Haus J. W. On the n-th mode numerical solutions of a second order boundary value problem on semi-infinite domain. Comm Appl Nonlinear Anal, 2007, 14: 121-133.

[47] Sheng Q., Khaliq A., Voss D. Numerical simulation of two-dimensional sine-Gordon solitons via a split cosine scheme. Math Computers Simulations, 2005, 68: 355-373.

[48] Strang G. On the construction and comparison of difference schemes. SIAM J Numer Anal, 1968, 5: 507-517.

[49] Suzuki M. General theory of fractal path integrals with applications to many-body theories and statistical physics. J Math Phys, 1991, 32: 400-407.

[50] Tabor M. The Sine-Gordon Equation, 1989, New York: Wiley.

[51] Tan Z., Zhang Z., Huang Y., Tang T. Moving mesh methods with locally varying time steps. J Comput Phys, 2004, 200: 347-367.

[52] Tang H. Z., Tang T. Adaptive mesh methods for one- and two-dimensional hyperbolic conservation laws. SIAM J Numer Anal, 2003, 41: 487-515.

[53] Thalhammer M. High-order exponential operator splitting methods for timedependent Schrödinger equations. SIAM J Numer Anal, 2008, 46: 2022-2038.

[54] Wang H., Liu J., Espedal M. S., Ewing R. E. A characteristic nonoverlapping domain decomposition method for multidimensional convection-diffusion equations. Numer Methods Partial Differential Eqns, 2005, 21: 89-103.