On the positivity, monotonicity, and stability of a semi-adaptive LOD method for solving three-dimensional degenerate Kawarada equations

Abstract This paper concerns the numerical solution of three-dimensional degenerate Kawarada equations. These partial differential equations possess highly nonlinear source terms, and exhibit strong quenching singularities which pose severe challenges to the design and analysis of highly reliable schemes. Arbitrary fixed nonuniform spatial grids, which are not necessarily symmetric, are considered throughout this study. The numerical solution is advanced through a semi-adaptive Local One-Dimensional (LOD) integrator. The temporal adaptation is achieved via a suitable arc-length monitoring mechanism. Criteria for preserving the positivity and monotonicity are investigated and acquired. The numerical stability of the splitting method is proven in the von Neumann sense under the spectral norm. Extended stability expectations are proposed and investigated.

[1]  Hong Cheng,et al.  Solving Degenerate Reaction-Diffusion Equations via Variable Step Peaceman-Rachford Splitting , 2004, SIAM J. Sci. Comput..

[2]  Robert D. Russell,et al.  A Study of Monitor Functions for Two-Dimensional Adaptive Mesh Generation , 1999, SIAM J. Sci. Comput..

[3]  ADI, LOD and Modern Decomposition Methods for Certain Multiphysics Applications , 2015 .

[4]  Andrew Acker,et al.  The quenching problem for nonlinear parabolic differential equations , 1976 .

[5]  Hideo Kawarada,et al.  On Solutions of Initial-Boundary Problem for ut=uxx+\frac{1}{1−u} , 1973 .

[6]  Jens Lang,et al.  An adaptive Rothe method for nonlinear reaction-diffusion systems , 1993 .

[7]  Edward H. Twizell,et al.  Chaos-free numerical solutions of reaction-diffusion equations , 1990, Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences.

[8]  Howard A. Levine,et al.  Quenching, nonquenching, and beyond quenching for solution of some parabolic equations , 1989 .

[9]  Tasso J. Kaper,et al.  N th-order operator splitting schemes and nonreversible systems , 1996 .

[10]  J. Verwer,et al.  A numerical study of three moving-grid methods for one-dimensional partial differential equations which are based on the method of lines , 1990 .

[11]  Hideo Kawarada,et al.  On Solutions of Initial-Bound ary Problem 1 for ut=uxx+- l-u , 1975 .

[12]  Michelle Schatzman,et al.  Stability of the Peaceman–Rachford Approximation , 1999 .

[13]  Qin Sheng,et al.  Adaptive decomposition finite difference methods for solving singular problems—A review , 2009 .

[14]  Q. Sheng and A.Q.M. Khaliq.  Linearly Implicit Adaptive Schemes for Singular Reaction-Diffusion Equations , 2001 .

[15]  Qin Sheng,et al.  A Revisit of the Semi-Adaptive Method for Singular Degenerate Reaction-Diffusion Equations , 2012 .

[16]  Jerrold Bebernes,et al.  Mathematical Problems from Combustion Theory , 1989 .

[17]  C. Y. Chan,et al.  Parabolic Quenching for Nonsmooth Convex Domains , 1994 .

[18]  N. N. Yanenko.  Splitting Methods for Partial Differential Equations , 1971, IFIP Congress.

[19]  Matthew A. Beauregard,et al.  An adaptive splitting approach for the quenching solution of reaction-diffusion equations over nonuniform grids , 2013, J. Comput. Appl. Math..

[20]  H. Ockendon.  Channel flow with temperature-dependent viscosity and internal viscous dissipation , 1979 .

[21]  Joseph E. Flaherty,et al.  On the stability of mesh equidistribution strategies for time-dependent partial differential equations , 1986 .

[22]  Gene H. Golub,et al.  Matrix computations , 1983 .

[23]  Andrew Acker,et al.  Remarks on quenching , 1988, Documenta Mathematica.