A new two-level implicit discretization of O ( k 2 + kh 2 + h 4 ) for the solution of singularly perturbed two-space dimensional non-linear parabolic equations

• Published in: Journal of computational and applied mathematics Vol. 208; no. 2; pp. 391 - 403
• Main Authors: Mohanty, R.K., Singh, Swarn
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier B.V 15.11.2007; Elsevier
• Subjects: Burgers’ equation; Exact sciences and technology; Functional analysis; Heat equation; High order compact scheme; Implicit scheme; Mathematical analysis; Mathematics; Non-linear parabolic equation; Normal derivatives; Numerical analysis; Numerical analysis. Scientific computation; Ordinary differential equations; Partial differential equations; Polar cylindrical coordinates; Sciences and techniques of general use; Singularly perturbed; High order compact scheme; Heat equation; Polar cylindrical coordinates; 65M06; Non-linear parabolic equation; Burgers’ equation; G.1.8; Implicit scheme; Singularly perturbed; Normal derivatives; Alternating direction method; Two dimensional space; Grid pattern; Initial condition; Differential equation; Perturbed equation; Boundary condition; Burgers equation; Algorithm; Discretization method; Exact solution; Non linear equation; 65M06;G.1.8; Parabolic equation; Numerical analysis; Linear equation; Applied mathematics; Implicit scheme; High order compact scheme; Singularly perturbed; Non-linear parabolic equation; Normal derivatives; Heat equation; Polar cylindrical coordinates; Burgers' equation; Two-dimensional calculations; Diffusion equation
• ISSN: 0377-0427; 1879-1778
• DOI: 10.1016/j.cam.2006.10.023

We propose a new two-level implicit difference method of O ( k 2 + kh 2 + h 4 ) for the solution of singularly perturbed non-linear parabolic differential equation ε ( u xx + u yy ) = f ( x , y , t , u , u x , u y , u t ) , 0 < x , y < 1 , t > 0 subject to appropriate initial and Dirichlet boundary conditions, where k > 0 and h > 0 are grid sizes in time and space directions, respectively, and ε > 0 is a small parameter. We also develop new methods of O ( kh 2 + h 4 ) for the estimates of ( ∂ u / ∂ x ) and ( ∂ u / ∂ y ) . In all cases, we use 9-spatial grid points and a single computational cell. The proposed methods are directly applicable to singular problems. We do not require any special scheme to solve singular problems. We also discuss alternating direction implicit (ADI) method for solving diffusion equation in polar cylindrical coordinates. This method permits multiple use of the one-dimensional tri-diagonal algorithm with a considerable saving in computing time, and produces a very efficient solver. It is shown that the ADI method is unconditionally stable. Numerical experiments are conducted to test the high accuracy of the proposed methods and compared with the exact solutions.