Numerical method satisfying the first two conservation laws for the Korteweg–de Vries equation

In this paper, we develop a finite-volume scheme for the KdV equation which conserves both the momentum and energy. The main ingredient of the method is a numerical device we developed in recent years that enables us to construct numerical method for a PDE that also simulates its related equations. In the method, numerical approximations to both the momentum and energy are conservatively computed. The operator splitting approach is adopted in constructing the method in which the conservation and dispersion parts of the equation are alternatively solved; our numerical device is applied in solving the conservation part of the equation. The feasibility and stability of the method is discussed, which involves an important property of the method, the so-called Jensen condition. The truncation error of the method is analyzed, which shows that the method is second-order accurate. Finally, several numerical examples, including the Zabusky–Kruskal’s example, are presented to show the good stability property of the method for long-time numerical integration.     

Soliton and kink jams in traffic flow with open boundaries

Soliton density wave is investigated numerically and analytically in the optimal velocity model (a car-following model) of a one-dimensional traffic flow with open boundaries. Soliton density wave is distinguished from the kink density wave. It is shown that the soliton density wave appears only at the threshold of occurrence of traffic jams. The Korteweg-de Vries (KdV) equation is derived from the optimal velocity model by the use of the nonlinear analysis. It is found that the traffic soliton appears only near the neutral stability line. The soliton solution is analytically obtained from the perturbed KdV equation. It is shown that the soliton solution obtained from the nonlinear analysis is consistent with that of the numerical simulation.     

The “backward looking” effect in the lattice hydrodynamic model

The novel lattice hydrodynamic model is presented by incorporating the “backward looking” effect. The stability condition for the the model is obtained using the linear stability theory. The result shows that considering one following site in vehicle motion leads to the stabilization of the system compared with the original lattice hydrodynamic model and the cooperative driving lattice hydrodynamic model. The Korteweg-de Vries (KdV, for short) equation near the neutral stability line is derived by using the reductive perturbation method to show the traffic jam which is proved to be described by KdV soliton solution obtained from the KdV equation. The simulation result is consistent with the nonlinear analysis.     

Analysis of car-following model considering driver’s physical delay in sensing headway

An extended car-following model is proposed by taking into account the delay of the driver’s response in sensing headway. The stability condition of this model is obtained by using the linear stability theory. The results show that the stability region decreases when the driver’s physical delay in sensing headway increases. The KdV equation and mKdV equation near the neutral stability line and the critical point are respectively derived by applying the reductive perturbation method. The traffic jams could be thus described by soliton solution and kink-antikink soliton solution for the KdV equation and mKdV equation respectively. The numerical results in the form of the space-time evolution of headway show that the stabilization effect is weakened when the driver’s physical delay increases. It confirms the fact that the delay of driver’s response in sensing headway plays an important role in jamming transition, and the numerical results are in good agreement with the theoretical analysis.     

Fermi-Pasta-Ulam, solitons and the fabric of nonlinear and computational science: history, synergetics, and visiometrics

This paper is mostly a history of the early years of nonlinear and computational physics and mathematics. I trace how the counterintuitive result of near-recurrence to an initial condition in the first scientific digital computer simulation led to the discovery of the soliton in a later computer simulation. The 1955 report by Fermi, Pasta, and Ulam (FPU) described their simulation of a one-dimensional nonlinear lattice which did not show energy equipartition. The 1965 paper by Zabusky and Kruskalshowed that the Korteweg-de Vries (KdV) nonlinear partial differential equation, a long wavelength model of the alpha-lattice (or cubic nonlinearity), derived by Kruskal, gave quantitatively the same results obtained by FPU. In 1967, Zabusky and Deem showed that a localized short wavelength initial excitation (then called an "optical" and now a "zone-boundary mode" excitation ) of the alpha-lattice revealed "n-curve" coherent states. If the initial amplitude was sufficiently large energy equipartition followed in a short time. The work of Kruskal and Miura (KM), Gardner and Greene (GG), and myself led to the appreciation of the infinity of denumerable invariants (conservation laws) for Hamiltonian systems and to a procedure by GGKM in 1967 for solving KdV exactly. The nonlinear science field exponentiated in diversity of linkages (as described in Appendix A). Included were pure and applied mathematics and all branches of basic and applied physics, including the first nonhydrodynamic application to optical solitons, as described in a brief essay (Appendix B) by Hasegawa. The growth was also manifest in the number of meetings held and institutes founded, as described briefly in Appendix D. Physicists and mathematicians in Japan, USA, and USSR (in the latter two, people associated with plasma physics) contributed to the diversification of the nonlinear paradigm which continues worldwide to the present. The last part of the paper (and Appendix C) discuss visiometrics: the visualization and quantification of simulation data, e.g., projection to lower dimensions, to facilitate understanding of nonlinear phenomena for modeling and prediction (or design). Finally, I present some recent developments that are linked to my early work by: Dritschel (vortex dynamics via contour dynamics/surgery in two and three dimensions); Friedland (pattern formation by synchronization in Hamiltonian nonlinear wave, vortex, plasma, systems, etc.); and the author ("n-curve" states and energy equipartition in a FPU lattice).     

A simple and direct method for generating travelling wave solutions for nonlinear equations

We propose a simple and direct method for generating travelling wave solutions for nonlinear integrable equations. We illustrate how nontrivial solutions for the KdV, the mKdV and the Boussinesq equations can be obtained from simple solutions of linear equations. We describe how using this method, a soliton solution of the KdV equation can yield soliton solutions for the mKdV as well as the Boussinesq equations. Similarly, starting with cnoidal solutions of the KdV equation, we can obtain the corresponding solutions for the mKdV as well as the Boussinesq equations. Simple solutions of linear equations can also lead to cnoidal solutions of nonlinear systems. Finally, we propose and solve some new families of KdV equations and show how soliton solutions are also obtained for the higher order equations of the KdV hierarchy using this method.     

Density waves in traffic flow model with relative velocity

The car-following model of traffic flow is extended to take into account the relative velocity. The stability condition of this model is obtained by using linear stability theory. It is shown that the stability of uniform traffic flow is improved by considering the relative velocity. From nonlinear analysis, it is shown that three different density waves, that is, the triangular shock wave, soliton wave and kink-antikink wave, appear in the stable, metastable and unstable regions of traffic flow respectively. The three different density waves are described by the nonlinear wave equations: the Burgers equation, Korteweg-de Vries (KdV) equation and modified Korteweg-de Vries (mKdV) equation, respectively.     

Nonliner Analysis of a Synthesized Optimal Velocity Model for Traffic Flow

We analyze a new car-following model described by a differential-difference equation with a synthesized optimal velocity function (SOVF), which depends on the front interactions between every two adjacent vehicles instead of the weighted average headway. The model is analyzed with the use of the linear stability theory and nonlinear analysis method. The stability and neutral stability condition are obtained. We also derive the modified KdV (Korteweg-de Vries) equation and the kink-antikink soliton solution near the critical point. A simulation is conducted with integrating the differential-difference equation by the Euler scheme. The results of the numerical simulation verify the validity of the new model.     

Two-component description of dynamical systems that can be approximated by solitons: The case of the ion acoustic wave equations of plasma physics

A new approach to the perturbative analysis of dynamical systems, which can be described approximately by soliton solutions of integrable non-linear wave equations, is employed in the case of small-amplitude solutions of the ion acoustic wave equations of plasma physics. Instead of pursuing the traditional derivation of a perturbed KdV equation, the ion velocity is written as a sum of two components: elastic and inelastic. In the single-soliton case, the elastic component is the full solution. In the multiple-soliton case, it is complemented by the inelastic component. The original system is transformed into two evolution equations: An asymptotically integrable Normal Form for ordinary KdV solitons, and an equation for the inelastic component. The zero-order term of the elastic component is a single-soliton or multiple-soliton solution of the Normal Form. The inelastic component asymptotes into a linear combination of single-soliton solutions of the Normal Form, with amplitudes determined by soliton interactions, plus a second-order decaying dispersive wave. Satisfaction of a conservation law by the inelastic component and of mass conservation by the disturbance to the ion density is determined solely by the initial data and/or boundary conditions imposed on the inelastic component. The electrostatic potential is a first-order quantity. It is affected by the inelastic component only in second order. The charge density displays a triple-layer structure. The analysis is carried out through the third order.     

Two velocity difference model for a car following theory

In the light of the optimal velocity model, a two velocity difference model for a car-following theory is put forward considering navigation in modern traffic. To our knowledge, the model is an improvement over the previous ones theoretically, because it considers more aspects in the car-following process than others. Then we investigate the property of the model using linear and nonlinear analyses. The Korteweg-de Vries equation (for short, the KdV equation) near the neutral stability line and the modified Korteweg-de Vries equation (for short, the mKdV equation) around the critical point are derived by applying the reductive perturbation method. The traffic jam could be thus described by the KdV soliton and the kink-anti-kink soliton for the KdV equation and mKdV equation, respectively. Numerical simulations are made to verify the model, and good results are obtained with the new model.     

The Korteweg-de Vires equation for the bidirectional pedestrian flow model considering the next-nearest-neighbor effect

This paper focuses on a two-dimensional bidirectional pedestrian flow model which involves the next-nearest-neighbor effect. The stability condition and the Korteweg-de Vries （KdV） equation are derived to describe the density wave of pedestrian congestion by linear stability and nonlinear analysis. Through theoretical analysis, the soliton solution is obtained.     

Multi-hump bright solitons in a Schrödinger–mKdV system

We consider the problem of energy transport in a Davydov model along an anharmonic crystal medium obeying quartic longitudinal interactions corresponding to rigid interacting particles. The Zabusky and Kruskal unidirectional continuum limit of the original discrete equations reduces, in the long wave approximation, to a coupled system between the linear Schrödinger (LS) equation and the modified Korteweg–de Vries (mKdV) equation. Single- and two-hump bright soliton solutions for this LS–mKdV system are predicted to exist by variational means and numerically confirmed. The one-hump bright solitons are found to be the anharmonic supersonic analogue of the Davydov's solitons while the two-hump (in both components) bright solitons are found to be a novel type of soliton consisting of a two-soliton solution of mKdV trapped by the wave function associated to the LS equation. This two-hump soliton solution, as a two component solution, represents a new class of polaron solution to be contrasted with the two-soliton interaction phenomena from soliton theory, as revealed by a variational approach and direct numerical results for the two-soliton solution.     

Engineering solitons and breathers in a deformed ferromagnet: Effect of localised inhomogeneities

We investigate the soliton dynamics of the electromagnetic wave propagating in an inhomogeneous or deformed ferromagnet. The dynamics of magnetization and the propagation of electromagnetic waves are governed by the Landau–Lifshitz–Maxwell (LLM) equation, a certain coupling between the Landau–Lifshitz and Maxwell's equations. In the framework of multiscale analysis, we obtain the perturbed integral modified KdV (PIMKdV) equation. Since the dynamic is governed by the nonlinear integro-differential equation, we rely on numerical simulations to study the interaction of its mKdV solitons with various types of inhomogeneities. Apart from simple one soliton experiments with periodic or localised inhomogeneities, the numerical simulations revealed an interesting dynamical scenario where the collision of two solitons on a localised inhomogeneity create a bound state which then produces either two separated solitons or a mKdV breather.     

Painlevé analysis,conservation laws,and symmetry of perturbed nonlinear equations

We consider the Lie-Backlund symmetries and conservation laws of a perturbed KdV equation and NLS equation. The arbitrary coefficients of the perturbing terms can be related to the condition of existence of nontrivial LB symmetry generator. When the perturbed KdV equation is subjected to Painlevé analysisa la Weiss, it is found that the resonance position changes compared to the unperturbed one. We prove the compatibility of the overdetermined set of equations obtained at the different stages of recursion relations, at least for one branch. All other branches are also indicated and difficulties associated them are discussed considering the perturbation parameter to be small. We determine the Lax pair for the aforesaid branch through the use of Schwarzian derivative. For the perturbed NLS equation we determine the conservation laws following the approach of Chen and Liu. From the recurrence of these conservation laws a Lax pair is constructed. But the Painlevé analysis does not produce a positive answer for the perturbed NLS equation. So here we have two contrasting examples of perturbed nonlinear equations: one passes the Painlevé test and its Lax pair can be found from the analysis itself, but the other equation does not meet the criterion of the Painlevé test, though its Lax pair is found in another way.     

拉曼增益对孤子传输特性的影响

利用考虑拉曼增益效应的非线性薛定谔方程, 在忽略光纤损耗的情况下, 采用基于MATLAB的分步傅里叶数值算法, 得出线性算符和非线性算符具体的表达式, 分步作用于光孤子脉冲传输方程, 仿真模拟了光孤子在光纤中传输时的演变. 与不考虑拉曼增益的光孤子在光纤中传输相对比, 探析了拉曼增益对孤子传输特性的影响.拉曼增益会破坏孤子的传输周期, 导致孤子在光纤中传输时快速衰减, 并且影响程度和输入孤子的脉冲峰值功率大小有关, 拉曼增益对基态孤子和高阶孤子的影响也不相同. 关键词： 拉曼增益 孤子 对称分步傅里叶法 非线性薛定谔方程     

Complete and orthonormal eigenfunctions of excitations to a one-dimensional dark soliton

We study linear excitations to a one-dimensional dark soliton described by a defocusing nonlinear Schödinger equation. By solving an eigenvalue problem for the excitations we obtain all eigenvalues and eigenfunctions and prove rigorously that these eigenfunctions are orthonormal and form a complete set. We then use the eigenfunctions to obtain the exact form of linear excitations for any given initial condition and to investigate the transverse stability of the dark soliton. The rigorous results reported in the present work can be applied to study the dynamics of dark solitons in various nonlinear optical media and Bose-Einstein condensates.     

A New Car Following Model： Comprehensive Optimal Velocity Model

In this paper, we present a new car-following model, i.e. comprehensive optimal velocity model (COVM), whose optimal velocity function not only depends on the following distance of the preceding vehicle, but also depends on the velocity difference with preceding vehicle. Simulation results show that COVM is an improvement over the previous ones theoretically. Then, the stability condition of the model is obtained by the linear stability analysis, which has shown that the model could obtain a bigger stable region thanprevious models in the phase diagram. Through the nonlinear analysis, the Burgers, Korteweg-de Vries (KdV) and modified KdV (mKdV) equations are derived for the triangular shock wave, the soliton wave, and the kink-antikink soliton wave. At the same time, numerical simulations are also carried out to show that the model could simulate these density waves.     

Conservation laws,modulation instability and solitons interactions for a nonlinear Schrödinger equation with the sextic operators in an optical fiber

Under investigation in this paper is a sextic nonlinear Schrödinger equation, which describes the pulses propagating along an optical fiber. Based on the symbolic computation, Lax pair and infinitely-many conservation laws are derived. Via the modiied Hirota method, bilinear forms and multi-soliton solutions are obtained. Propagation and interactions of the solitons are illustrated graphically: Initial position and velocity of the soliton are related to the coefficient of the sixth-order dispersion, while the amplitude of the soliton is not affected by it. Head-on, overtaking and oscillating interactions between the two solitons are displayed. Through the asymptotic analysis, interaction between the two solitons is proved to be elastic. Based on the linear stability analysis, the modulation instability condition for the soliton solutions is obtained.     

The Korteweg-de Vries soliton in the lattice hydrodynamic model

The lattice hydrodynamic model is not only a simplified version of the macroscopic hydrodynamic model, but is also closely connected with the microscopic car following model. The modified Korteweg-de Vries (mKdV) equation about the density wave in congested traffic has been derived near the critical point since Nagatani first proposed it. But the Korteweg-de Vries (KdV) equation near the neutral stability line has not been studied, which has been investigated in detail in the car following model. So we devote ourselves to obtaining the KdV equation from the lattice hydrodynamic model and obtaining the KdV soliton solution describing the traffic jam. Numerical simulation is conducted, to demonstrate the nonlinear analysis result.     

Nonlinear analysis of traffic flow on a gradient highway

We investigate the slope effects upon traffic flow on a single lane gradient (uphill/downhill) highway analytically and numerically. The stability condition, neutral stability condition and instability condition are obtained by the use of linear stability theory. It is found that stability of traffic flow on the gradient varies with the slopes. The Burgers, Korteweg-de Vries (KdV) and modified Korteweg-de Vries (mKdV) equations are derived to describe the triangular shock waves, soliton waves and kink-antikink waves in the stable, meta-stable and unstable region respectively. A series of simulations are carried out to reproduce the triangular shock waves, kink-antikink waves and soliton waves. Results show that amplitudes of the triangular shock waves and kink-antikink waves vary with the slopes, the soliton wave appears in an upward form when the average headway is less than the safety distance and a downward form when the average headway is more than the safety distance. Moreover both the kink-antikink waves and the solitary waves propagate backwards. The numerical simulation shows a good agreement with the analytical result.