Improved tanh-function method and the new exact solutions for the general variab le coefficient KdV equation and MKdV equation

An infinite number of form-invariant symmetries is obtained and a one-to-one correspondence between symmetries and conservation laws is established for the non-isospectral and variable coefficient generalizations of both the KdV and the MKdV equations. Two families of symmetries and their Lie algebraic structures are constructed. Some interesting facts about their Hamiltonian structures are presented.

First-order ordinary differential equations have some rather special properties, which result for the most part because they can only contain a limited number of terms. In fact, all linear first-order ordinary differential equations can easily be solved. This is in contrast to higher-order ordinary differential equations that become much more difficult to solve when, for example, they contain variable coefficients. Furthermore, some methods exist to solve a pretty large class of nonlinear first-order ordinary differential equations, which is generally not the case for higherorder nonlinear equations.

Based on a Riccati equation and a symbolic computation system-Maple, a generalized Riccati equation expansion method is presented for constructing soliton-like solutions and periodic form solutions for some nonlinear evolution equations (NEEs) or NEEs with variable coefficients. Compared with most of the existing tanh methods, the extended tanh-function method, the modified extended tanh-function method and the generalized hyperbolic-function method, the proposed method is more powerful. We study a (2+1)-dimensional general nonintegrable KdV equation, a KdV equation with variable coefficients. As a result, rich new families of exact solutions, including the non-travelling wave's and coefficient functions' soliton-like solutions, singular soliton-like solutions, periodic form solutions, are obtained. When setting the arbitrary functions in some solutions be equal to special constants or special functions, the solitary wave solutions can be recovered.

Fractal fractional analysis of modified KdV equation under three different kernels

Chapter 2 - Ordinary differential equations and recurrence relations

We consider the equivalence problem for first-order ordinary differential equations (ODEs) and recover the well-known result that any two first-order ODEs are equivalent under a point transformation by Cartan’s equivalence method. Curvature properties of some first-order and second-order ODEs motivate us to define the equivalence problem for such equations as the metric equivalence problem which defines the [Formula: see text]-structure on the corresponding manifolds. Accordingly, we show that there are first-order ODEs and linear second-order ODEs which may not be equivalent under a class of fiber-preserving diffeomorphisms.

Compound KdV equation and MKdV equation are investigated in the presented work, with the aid of Matlab software, using the trigonometric function transform method and the Wu elimination method, and more exact traveling wave solutions are obtained for compound KdV equation and MKdV equation.

SUMMARY The global asymptotic nonlinear behavior of 11 explicit and implicit time discretizations for four 2 × 2 systems of first-order autonomous nonlinear ordinary differential equations (ODEs) is analyzed. The objectives are to gain a basic understanding of the difference in the dynamics of numerics between the scalars and systems of nonlinear autonomous ODEs and to set a baseline global asymptotic solution behavior of these schemes for practical computations in computational fluid dynamics. We show how "numerical" basins of attraction can complement the bifurcation diagrams in gaining more detailed global asymptotic behavior of time discretizations for nonlinear differential equations (DCs). We show how in the presence of spurious asymptotes the basins of the true stable steady states can be segmented by the basins of the spurious stable and unstable asymptotes. One major consequence of this phenomenon which is not commonly known is that this spurious behavior can result in a dramatic distortion and, in most cases, a dramatic shrinkage and segmentation of the basin of attraction of the true solution for finite time steps. Such distortion, shrinkage and segmentation of the numerical basins of attraction will occur regardless ofthe stability ofthe spurious asymptotes, and will occur for unconditionally stable implicit linear multistep methods. In other words, for the same (common) steady-state solution the associated basin of attraction of the DE might be very different from the discretized counterparts and the numerical basin of attraction can be very different from numerical method to numerical method. The results can be used as an explanation for possible causes of error, and slow convergence and nonconvergence of steady-state numerical solutions when using the time-dependent approach for nonlinear hyperbolic or parabolic PDEs.

The effect of changing the direction of motion of a defect (a soliton of small amplitude) in soliton lattices described by the Korteweg–de Vries and modified Korteweg–de Vries integrable equations (KdV and mKdV) was studied. Manifestation of this effect is possible as a result of the negative phase shift of a small soliton at the moment of nonlinear interaction with large solitons, as noted in [1], within the KdV equation. In the recent paper [2], an expression for the mean soliton velocity in a “cold” KdV-soliton gas has been found using kinetic theory, from which this effect also follows, but this fact has not been mentioned. In the present paper, we will show that the criterion of negative velocity is the same for both the KdV and mKdV equations and it can be obtained using simple kinematic considerations without applying kinetic theory. The averaged dynamics of the “smallest” soliton (defect) in a soliton gas consisting of solitons with random amplitudes has been investigated and the average criterion of changing the sign of the velocity has been derived and confirmed by numerical solutions of the KdV and mKdV equations.

An approximate method for solving the Cauchy problem for nonlinear first-order ordinary differential equations is considered. The method is based on using the shifted Chebyshev series and a Markov quadrature formula. Some approaches are given to estimate the error of an approximate solution expressed by a partial sum of a certain order series. The error is estimated using the second approximation of the solution expressed by a partial sum of a higher order series. An algorithm for partitioning the integration interval into elementary subintervals to ensure the computation of the solution with a prescribed accuracy is constructed on the basis of the proposed approaches to error estimation.

In this paper, we present an integrable semi-discretization of the modified Korteweg-deVries (mKdV) equation. We discretize the “time” variable of the mKdV equation and get an integrable differential-difference system. Under a standard limit, the differential-difference system converges to the continuous mKdV equation. By Hirota’s bilinear method, we find some explicit solutions including solitons and breather solutions. From the semi-discrete system, we design a numerical scheme to the mKdV equation and carry out several numerical experiments with the 3-soliton solution and breather solution.

Damped and driven oscillators are generally modeled with a nonautonomous second-order nonlinear ordinary differential equation including a sinusoidal driving forcing term, such as the forced Duffing equation and the forced Holmes–Rand equation. These equations have been extensively studied during the last century and the last two decades. In the early 1990s, Abarbanel, Rabinovich and Sushchik proposed replacing the sinusoidal forcing term with a “force controlled by the movements of the oscillator itself”, i.e. by the product of two variables: the first being the solution of the oscillator itself, while the second is the solution of a first-order nonlinear ordinary differential equation. They referred to the resulting autonomous dynamical system of two coupled nonlinear ordinary differential equations as a “controlled nonlinear oscillator”. To that end, they introduced a change of variables and parameters to transform the “controlled nonlinear oscillator” that corresponds to a particular case of the forced Duffing equation into the Lorenz system. The aim of this work is to show that their idea can be further generalized and applied to many other dynamical systems, including the forced Holmes–Rand equation, Chua’s cubic circuit, Chen’s system and the forced Helmholtz oscillator. It is proved that a certain class of three-dimensional dynamical systems can be rewritten into the form of “generalized controlled nonlinear oscillators”, which can then be transformed into various Lorenz-like systems. Such a transformation could be very useful for the study of intermittent chaos.

Post-buckling behaviour of slender structures with a bi-linear bending moment–curvature relationship

In this paper, a mathematical model of a prey-predator system with infectious disease in the prey population is proposed and studied. It is assumed that there is a constant refuge in prey as a defensive property against predation and harvesting from the predator. The proposed mathematical model is consisting of three first-order nonlinear ordinary differential equations, which describe the interaction among the healthy prey, infected prey, and predator. The existence, uniqueness, and boundedness of the system’ solution are investigated. The system's equilibrium points are calculated with studying their local and global stability. The persistence conditions of the proposed system are established. Finally the obtained analytical results are justified by a numerical simulation.

Einstein's equations for a rotating pressure-free space-time are reduced to a system of four first-order non-linear ordinary differential equations in one self-similar dimensionless variable. Numerical results are given for the vacuum solution. A compatible thin disk can be specified by a surface density σ and an angular velocity ω. Self-similarity as a statement of the absence of scales implies that σ and ω can be written asσ=αc2/4πGr, ω=βc/r, and demands that α and β be pure numbers.