Chaotic scattering in solitary wave interactions: A singular iterated-map description

We derive a family of singular iterated maps--closely related to Poincare maps--that describe chaotic interactions between colliding solitary waves. The chaotic behavior of such solitary wave collisions depends on the transfer of energy to a secondary mode of oscillation, often an internal mode of the pulse. Unlike previous analyses, this map allows one to understand the interactions in the case when this mode is excited prior to the first collision. The map is derived using Melnikov integrals and matched asymptotic expansions and generalizes a ``multi-pulse'' Melnikov integral and allows one to find not only multipulse heteroclinic orbits, but exotic periodic orbits. The family of maps derived exhibits singular behavior, including regions of infinite winding. This problem is shown to be a singular version of the conservative Ikeda map from laser physics and connections are made with problems from celestial mechanics and fluid mechanics.Comment: 29 pages, 17 figures, submitted to Chaos, higher-resolution figures available at author's website: http://m.njit.edu/goodman/publication

R-adaptive multisymplectic and variational integrators

Moving mesh methods (also called r-adaptive methods) are space-adaptive strategies used for the numerical simulation of time-dependent partial differential equations. These methods keep the total number of mesh points fixed during the simulation, but redistribute them over time to follow the areas where a higher mesh point density is required. There are a very limited number of moving mesh methods designed for solving field-theoretic partial differential equations, and the numerical analysis of the resulting schemes is challenging. In this paper we present two ways to construct r-adaptive variational and multisymplectic integrators for (1+1)-dimensional Lagrangian field theories. The first method uses a variational discretization of the physical equations and the mesh equations are then coupled in a way typical of the existing r-adaptive schemes. The second method treats the mesh points as pseudo-particles and incorporates their dynamics directly into the variational principle. A user-specified adaptation strategy is then enforced through Lagrange multipliers as a constraint on the dynamics of both the physical field and the mesh points. We discuss the advantages and limitations of our methods. Numerical results for the Sine-Gordon equation are also presented.Comment: 65 pages, 13 figure

Decomposition Method for Kdv Boussinesq and Coupled Kdv Boussinesq Equations

This paper obtains the solitary wave solutions of two different forms of Boussinesq equations that model the study of shallow water waves in lakes and ocean beaches. The decomposition method using He’s polynomials is applied to solve the governing equations. The travelling wave hypothesis is also utilized to solve the generalized case of coupled Boussinesq equations, and, thus, an exact soliton solution is obtained. The results are also supported by numerical simulations. Keywords: Decomposition Method, He’s polynomials, cubic Boussinesq equation, Coupled Boussinesq equation

Some Problems on Variational Iteration Method

In this research project paper, I introduce some basic idea of Variational iteration method and its algorithm to solve the equations ODE & PDE, fractional differential equation, fractal differential equation and differential-difference equations. Also, some linear and nonlinear differential equations like Burger’s equation, Fisher’s equation,Wave equation and Schrodinger equation are solved by using Variational iteration method. Then I compare this method with Adomian decomposition method (ADM) and modified Variational iteration method (MVIM). The advantage of VIM, it does not require a small parameter in an equation as perturbation technique needs. The VIM is used to solve effectively, easily, and accurately a large class of non-linear problems with approximations which converge rapidly to accurate solutions. For linear problems, its exact solution can be obtained by only one iteration step due to the fact that the Lagrange multiplier can be exactly identified

Wavelet Methods for the Solutions of Partial and Fractional Differential Equations Arising in Physical Problems

The subject of fractional calculus has gained considerable popularity and importance during the past three decades or so, mainly due to its demonstrated applications in numerous seemingly diverse and widespread fields of science and engineering. It deals with derivatives and integrals of arbitrary orders. The fractional derivative has been occurring in many physical problems, such as frequency-dependent damping behavior of materials, motion of a large thin plate in a Newtonian fluid, creep and relaxation functions for viscoelastic materials, the PI D controller for the control of dynamical systems etc. Phenomena in electromagnetics, acoustics, viscoelasticity, electrochemistry, control theory, neutron point kinetic model, anomalous diffusion, Brownian motion, signal and image processing, fluid dynamics and material science are well described by differential equations of fractional order. Generally, nonlinear partial differential equations of fractional order are difficult to solve. So for the last few decades, a great deal of attention has been directed towards the solution (both exact and numerical) of these problems. The aim of this dissertation is to present an extensive study of different wavelet methods for obtaining numerical solutions of mathematical problems occurring in disciplines of science and engineering. This present work also provides a comprehensive foundation of different wavelet methods comprising Haar wavelet method, Legendre wavelet method, Legendre multi-wavelet methods, Chebyshev wavelet method, Hermite wavelet method and Petrov-Galerkin method. The intension is to examine the accuracy of various wavelet methods and their efficiency for solving nonlinear fractional differential equations. With the widespread applications of wavelet methods for solving difficult problems in diverse fields of science and engineering such as wave propagation, data compression, image processing, pattern recognition, computer graphics and in medical technology, these methods have been implemented to develop accurate and fast algorithms for solving integral, differential and integro-differential equations, especially those whose solutions are highly localized in position and scale. The main feature of wavelets is its ability to convert the given differential and integral equations to a system of linear or nonlinear algebraic equations, which can be solved by numerical methods. Therefore, our main focus in the present work is to analyze the application of wavelet based transform methods for solving the problem of fractional order partial differential equations. The introductory concept of wavelet, wavelet transform and multi-resolution analysis (MRA) have been discussed in the preliminary chapter. The basic idea of various analytical and numerical methods viz. Variational Iteration Method (VIM), Homotopy Perturbation Method (HPM), Homotopy Analysis Method (HAM), First Integral Method (FIM), Optimal Homotopy Asymptotic Method (OHAM), Haar Wavelet Method, Legendre Wavelet Method, Chebyshev Wavelet Method and Hermite Wavelet Method have been presented in chapter 1. In chapter 2, we have considered both analytical and numerical approach for solving some particular nonlinear partial differential equations like Burgers’ equation, modified Burgers’ equation, Huxley equation, Burgers-Huxley equation and modified KdV equation, which have a wide variety of applications in physical models. Variational Iteration Method and Haar wavelet Method are applied to obtain the analytical and numerical approximate solution of Huxley and Burgers-Huxley equations. Comparisons between analytical solution and numerical solution have been cited in tables and also graphically. The Haar wavelet method has also been applied to solve Burgers’, modified Burgers’, and modified KdV equations numerically. The results thus obtained are compared with exact solutions as well as solutions available in open literature. Error of collocation method has been presented in this chapter. Methods like Homotopy Perturbation Method (HPM) and Optimal Homotopy Asymptotic Method (OHAM) are very powerful and efficient techniques for solving nonlinear PDEs. Using these methods, many functional equations such as ordinary, partial differential equations and integral equations have been solved. We have implemented HPM and OHAM in chapter 3, in order to obtain the analytical approximate solutions of system of nonlinear partial differential equation viz. the Boussinesq-Burgers’ equations. Also, the Haar wavelet method has been applied to obtain the numerical solution of BoussinesqBurgers’ equations. Also, the convergence of HPM and OHAM has been discussed in this chapter. The mathematical modeling and simulation of systems and processes, based on the description of their properties in terms of fractional derivatives, naturally leads to differential equations of fractional order and the necessity to solve such equations. The mathematical preliminaries of fractional calculus, definitions and theorems have been presented in chapter 4. Next, in this chapter, the Haar wavelet method has been analyzed for solving fractional differential equations. The time-fractional Burgers-Fisher, generalized Fisher type equations, nonlinear time- and space-fractional Fokker-Planck equations have been solved by using two-dimensional Haar wavelet method. The obtained results are compared with the Optimal Homotopy Asymptotic Method (OHAM), the exact solutions and the results available in open literature. Comparison of obtained results with OHAM, Adomian Decomposition Method (ADM), VIM and Operational Tau Method (OTM) has been demonstrated in order to justify the accuracy and efficiency of the proposed schemes. The convergence of two-dimensional Haar wavelet technique has been provided at the end of this chapter. In chapter 5, the fractional differential equations such as KdV-Burger-Kuramoto (KBK) equation, seventh order KdV (sKdV) equation and Kaup-Kupershmidt (KK) equation have been solved by using two-dimensional Legendre wavelet and Legendre multi-wavelet methods. The main focus of this chapter is the application of two-dimensional Legendre wavelet technique for solving nonlinear fractional differential equations like timefractional KBK equation, time-fractional sKdV equation in order to demonstrate the efficiency and accuracy of the proposed wavelet method. Similarly in chapter 6, twodimensional Chebyshev wavelet method has been implemented to obtain the numerical solutions of the time-fractional Sawada-Kotera equation, fractional order Camassa-Holm equation and Riesz space-fractional sine-Gordon equations. The convergence analysis has been done for these wavelet methods. In chapter 7, the solitary wave solution of fractional modified Fornberg-Whitham equation has been attained by using first integral method and also the approximate solutions obtained by optimal homotopy asymptotic method (OHAM) are compared with the exact solutions acquired by first integral method. Also, the Hermite wavelet method has been implemented to obtain approximate solutions of fractional modified Fornberg-Whitham equation. The Hermite wavelet method is implemented to system of nonlinear fractional differential equations viz. the fractional Jaulent-Miodek equations. Convergence of this wavelet methods has been discussed in this chapter. Chapter 8 emphasizes on the application of Petrov-Galerkin method for solving the fractional differential equations such as the fractional KdV-Burgers’ (KdVB) equation and the fractional Sharma-TassoOlver equation with a view to exhibit the capabilities of this method in handling nonlinear equation. The main objective of this chapter is to establish the efficiency and accuracy of Petrov-Galerkin method in solving fractional differential equtaions numerically by implementing a linear hat function as the trial function and a quintic B-spline function as the test function. Various wavelet methods have been successfully employed to numerous partial and fractional differential equations in order to demonstrate the validity and accuracy of these procedures. Analyzing the numerical results, it can be concluded that the wavelet methods provide worthy numerical solutions for both classical and fractional order partial differential equations. Finally, it is worthwhile to mention that the proposed wavelet methods are promising and powerful methods for solving fractional differential equations in mathematical physics. This work also aimed at, to make this subject popular and acceptable to engineering and science community to appreciate the universe of wonderful mathematics, which is in between classical integer order differentiation and integration, which till now is not much acknowledged, and is hidden from scientists and engineers. Therefore, our goal is to encourage the reader to appreciate the beauty as well as the usefulness of these numerical wavelet based techniques in the study of nonlinear physical system