A finite difference approach for the numerical solution of non-smooth problems for Boundary Value ODEs

This paper concerns the numerical approximation of Boundary Value ODEs (BVPs) with non-smooth coefficients and solutions. Different strategies are presented to tackle the cases of known and unknown singularity locations. In the former, the original problem is transformed in a multipoint BVP and high order Extended Central Difference Formulas (ECDFs) are used to approximate the smooth branches of the solution and the Neumann boundary conditions (BCs) with same accuracy. In the latter, an iterative Hybrid method coupling ECDFs and the shooting technique has been introduced to approximate both the discontinuity point and the solution. Convergence analysis and numerical comparisons with other approaches from literature are also presented. Good performances in terms of errors and convergence order are reported by applying ECDFs and the Hybrid method to linear test BVPs and to a nonlinear bio-mechanical model, in both cases of mixed and Dirichlet BCs.

[1]  Eusebius J. Doedel,et al.  Finite Difference Methods for Singular Two-Point Boundary Value Problems , 1984 .

[2]  Manoj Kumar,et al.  Methods for solving singular boundary value problems using splines: a review , 2010 .

[3]  R. LeVeque,et al.  A comparison of the extended finite element method with the immersed interface method for elliptic equations with discontinuous coefficients and singular sources , 2006 .

[4]  Andreas Wiegmann,et al.  The Explicit-Jump Immersed Interface Method: Finite Difference Methods for PDEs with Piecewise Smooth Solutions , 2000, SIAM J. Numer. Anal..

[5]  Ole Østerby,et al.  Solving Ordinary Differential Equations with Discontinuities , 1984, TOMS.

[6]  M. K. Kadalbajoo,et al.  Numerical solution of singular boundary value problems via Chebyshev polynomial and B-spline , 2005, Appl. Math. Comput..

[7]  Pierluigi Amodio,et al.  High order generalized upwind schemes and numerical solution of singular perturbation problems , 2007 .

[8]  Pierluigi Amodio,et al.  High-order finite difference schemes for the solution of second-order BVPs , 2005 .

[9]  M. K. Kadalbajoo,et al.  Cubic spline for solving singular two-point boundary value problems , 2004, Appl. Math. Comput..

[10]  Xiaolin Zhong,et al.  A new high-order immersed interface method for solving elliptic equations with imbedded interface of discontinuity , 2007, J. Comput. Phys..

[11]  Vincent Acary,et al.  Numerical Methods for Nonsmooth Dynamical Systems: Applications in Mechanics and Electronics , 2008 .

[12]  Winfried Auzinger,et al.  Analysis of a New Error Estimate for Collocation Methods Applied to Singular Boundary Value Problems , 2005, SIAM J. Numer. Anal..

[13]  Shan Zhao,et al.  High order matched interface and boundary method for elliptic equations with discontinuous coefficients and singular sources , 2006, J. Comput. Phys..

[14]  Ivonne Sgura,et al.  Inhomogeneous shear of orthotropic incompressible non-linearly elastic solids: Singular solutions and biomechanical interpretation , 2009 .

[15]  Lawrence F. Shampine,et al.  Singular boundary value problems for ODEs , 2003, Appl. Math. Comput..

[16]  W. H. Enright,et al.  Effective solution of discontinuous IVPs using a Runge-Kutta formula pair with interpolants , 1988 .

[17]  L. Shampine,et al.  Event location for ordinary differential equations , 2000 .

[18]  Ivonne Sgura,et al.  A Mathematical Model for the Corrosion of Metallic Bipolar Plates in PEM Fuel Cells: Numerical and Experimental Issues , 2009, SIAM J. Appl. Math..

[19]  Deborah Lacitignola,et al.  Travelling waves in a reaction-diffusion model for electrodeposition , 2011, Math. Comput. Simul..

[20]  Reinhold Mannshardt.  One-step methods of any order for ordinary differential equations with discontinuous right-hand sides , 1978 .