An efficient numerical treatment for the asymptotic behaviour of the nonlinear Airy-type problems

Abstract This study focuses on symplectic integrators for numerical evaluation of the asymptotic solutions of the nonlinear Airy-type equations obtained by reducing the nonlinear dispersive equations. Since the nature of Airy-type equations has both highly oscillatory slow decay and exponential fast decay, most of classical integrators are not able to correctly exhibit challenging physical behaviour. We use specially designed symplectic integrators combining splitting methods with Magnus integrators to catch asymptotic behaviour of nonlinear Airy-type equations efficiently, even for large step sizes. Efficiency of the proposed methods for given problems is discussed. Moreover, numerical results obtained by the proposed methods are compared with the existing results in the literature.

[1]  Tasawar Hayat,et al.  On comparison of series and numerical solutions for second Painlevé equation , 2010 .

[2]  Mechthild Thalhammer,et al.  A fourth-order commutator-free exponential integrator for nonautonomous differential equations , 2006, SIAM J. Numer. Anal..

[3]  Yuhao Cong,et al.  A sixth order diagonally implicit symmetric and symplectic Runge-Kutta method for solving Hamiltonian systems , 2015 .

[4]  Hüseyin Koçak,et al.  Similarity solutions of nonlinear third-order dispersive PDEs: The first critical exponent , 2017, Appl. Math. Lett..

[5]  S. Blanes,et al.  The Magnus expansion and some of its applications , 2008, 0810.5488.

[6]  Sergio Blanes,et al.  Fourier methods for the perturbed harmonic oscillator in linear and nonlinear Schrödinger equations. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[7]  Fernando Casas,et al.  Symplectic integrators for second-order linear non-autonomous equations , 2018, J. Comput. Appl. Math..

[8]  C. Klein,et al.  Numerical study of a multiscale expansion of the Korteweg–de Vries equation and Painlevé-II equation , 2007, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[9]  Sergio Blanes,et al.  Symplectic integrators for the matrix Hill equation , 2017, J. Comput. Appl. Math..

[10]  Lun Zhang,et al.  Location of Poles for the Hastings–McLeod Solution to the Second Painlevé Equation , 2014, 1410.3338.

[11]  Fernando Casas,et al.  Optimized high-order splitting methods for some classes of parabolic equations , 2011, Math. Comput..

[12]  W. Magnus.  On the exponential solution of differential equations for a linear operator , 1954 .

[13]  S. P. Hastings,et al.  A boundary value problem associated with the second painlevé transcendent and the Korteweg-de Vries equation , 1980 .

[14]  Very Singular Similarity Solutions and Hermitian Spectral Theory for Semilinear Odd-Order PDEs , 2009, 0910.4916.

[15]  Alexander Ostermann,et al.  Exponential splitting for unbounded operators , 2009, Math. Comput..

[16]  Arieh Iserles,et al.  On the Global Error of Discretization Methods for Highly-Oscillatory Ordinary Differential Equations , 2002 .

[17]  A. Dragt,et al.  Lie methods for nonlinear dynamics with applications to accelerator physics , 2011 .

[18]  Fasma Diele,et al.  Splitting and composition methods for explicit time dependence in separable dynamical systems , 2010, J. Comput. Appl. Math..

[19]  H. Koçak.  Blow-up and Global Similarity Solutions for Semilinear Third-Order Dispersive PDEs , 2015 .

[20]  Bruce R. Fabijonas,et al.  Algorithm 838: Airy Functions , 2004, TOMS.

[21]  Holger Fehske,et al.  High-order commutator-free exponential time-propagation of driven quantum systems , 2011, J. Comput. Phys..

[22]  Bengt Fornberg,et al.  A numerical methodology for the Painlevé equations , 2011, J. Comput. Phys..

[23]  Jacques Laskar,et al.  New families of symplectic splitting methods for numerical integration in dynamical astronomy , 2012, 1208.0689.

[24]  R. Rosales.  The similarity solution for the Korteweg-de Vries equation and the related Painlevé transcendent , 1978, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[25]  Bengt Fornberg,et al.  A Computational Exploration of the Second Painlevé Equation , 2014, Found. Comput. Math..

[26]  H. Koçak,et al.  Numerical study of the asymptotics of the second Painlevé equation by a functional fitting method , 2013 .

[27]  A. Iserles,et al.  Explicit Magnus expansions for nonlinear equations , 2006 .

[28]  S. Blanes,et al.  Fourth-and sixth-order commutator-free Magnus integrators for linear and non-linear dynamical systems , 2006 .

[29]  S. Blanes,et al.  High-order splitting methods for separable non-autonomous parabolic equations , 2013, 1311.1023.

[30]  Alexander Ostermann,et al.  Magnus integrators on multicore CPUs and GPUs , 2017, Comput. Phys. Commun..

[31]  Stéphane Descombes,et al.  Splitting methods with complex times for parabolic equations , 2009 .

[32]  J. Miles.  The Second Painlevé Transcendent: A Nonlinear Airy Function , 1980 .

[33]  Uniform semiclassical approximations of the nonlinear Schrödinger equation by a Painlevé mapping , 2006, quant-ph/0608099.

[34]  S. Blanes,et al.  Practical symplectic partitioned Runge--Kutta and Runge--Kutta--Nyström methods , 2002 .