On the numerical solution to a nonlinear wave equation associated with the first painlevé equation: an operator-splitting approach

Abstract

The main goal of this article is to discuss the numerical solution to a nonlinear wave equation associated with the first of the celebrated Painlevé transcendent ordinary differential equations. In order to solve numerically the above equation, whose solutions blow up in finite time, the authors advocate a numerical methodology based on the Strang’s symmetrized operator-splitting scheme. With this approach, one can decouple nonlinearity and differential operators, leading to the alternate solution at every time step of the equation as follows: (i) The first Painlevé ordinary differential equation, (ii) a linear wave equation with a constant coefficient. Assuming that the space dimension is two, the authors consider a fully discrete variant of the above scheme, where the space-time discretization of the linear wave equation sub-steps is achieved via a Galerkin/finite element space approximation combined with a second order accurate centered time discretization scheme. To handle the nonlinear sub-steps, a second order accurate centered explicit time discretization scheme with adaptively variable time step is used, in order to follow accurately the fast dynamic of the solution before it blows up. The results of numerical experiments are presented for different coefficients and boundary conditions. They show that the above methodology is robust and describes fairly accurately the evolution of a rather “violent” phenomenon.

Keywords

Painlevé equation / Nonlinear wave equation / Blow-up solution / Operator-Splitting

References

[1] Bornemann F, Clarkson P, Deift P A request: The Painlevé Project. Notices of the American Mathematical Society, 2010, 57(11): 1938

[2] Jimbo M. Monodromy problem and the boundary condition for some Painlevé equations. Publ. Res. Inst. Sci., 1982, 18(3): 1137-1161

[3] Wong R, Zhang H Y. On the connection formulas of the fourth Painlevé transcendent. Analysis and Applications, 2009, 7(4): 419-448

[4] Clarkson, P. A., Painlevé transcendents, NIST Handbook of Mathematical Functions, F. W. J. Olver, D. W. Lozier, R. F. Boisvert and C. W. Clark (eds)., Cambridge University Press, Cambridge, UK, 2010, 723–740.

[5] Fornberg B, Weideman J A C. A numerical methodology for the Painlevé equations. J. Comp. Phys., 2011, 230(15): 5957-5973

[6] Strang G. On the construction and comparison of difference schemes. SIAM J. Numer. Anal., 1968, 5(3): 506-517

[7] Glowinski, R., Finite element methods for incompressible viscous flow, Handbook of Numerical Analysis, Vol. IX, P. G. Ciarlet and J. L. Lions (eds.), North-Holland, Amsterdam, 2003, 3–1176.

[8] Bokil V A, Glowinski R. An operator-splitting scheme with a distributed Lagrange multiplier based fictitious domain method for wave propagation problems. J. Comput. Phys., 2005, 205(1): 242-268

[9] Glowinski R, Shiau L, Sheppard M. Numerical methods for a class of nonlinear integro-differential equations. Calcolo, 2012

[10] Lions J L. Quelques Méthodes de Résolution des Problèmes aux Limites Non Linéaires, 1969, Paris: Dunod

[11] Glowinski R, Dean E J, Guidoboni G Application of operator-splitting methods to the direct numerical simulation of particulate and free-surface flows and to the numerical solution to the two-dimensional elliptic Monge-Ampère equation. Japan J. Indust. Appl. Math., 2008, 25(1): 1-63

[12] Chorin A J, Hughes T J R, McCracken M F Product formulas and numerical algorithms. Comm. Pure Appl. Math., 1978, 31: 205-256

[13] Beale J T, Majda A. Rates of convergence for viscous splitting of the Navier-Stokes equations. Math. Comp., 1981, 37(156): 243-259

[14] Leveque R J, Oliger J. Numerical methods based on additive splittings for hyperbolic partial differential equations. Math. Comp., 1983, 40(162): 469-497

[15] Marchuk, G. I., Splitting and alternating direction method, Handbook of Numerical Analysis, Vol. I, P. G. Ciarlet and J. L. Lions (eds.), North-Holland, Amsterdam, 1990, 197–462.

[16] Temam R. Navier-Stokes Equations: Theory and Numerical Analysis, 2001, Providence, RI: AMS Chelsea Publish

[17] Glowinski R. Numerical Methods for Nonlinear Variational Problems, 1984, New York: Springer-Verlag

[18] http://en.wikipedia.org/wiki/Painleve_transcendents

[19] Keller J B. On solutions of nonlinear wave equations. Comm. Pure Appl. Math., 1957, 10(4): 523-530