Nonlinear Sciences

Vortex coherent structures on arrays of nonlinear oscillators joined by weak links into topologically nontrivial two-dimensional discrete manifolds have been theoretically studied. A circuit of nonlinear electric oscillators coupled by relatively weak capacitances has been considered as a possible physical implementation of such objects. Numerical experiments have shown that a time-monochromatic external force applied to several oscillators leads to the formation of long-lived and nontrivially interacting vortices in the system against a quasistationary background in a wide range of parameters. The dynamics of vortices depends on the method of "coupling" of the opposite sides of a rectangular array by links, which determines the topology of the resulting manifold (torus, Klein bottle, projective plane, Möbius strip, ring, or disk).

We touch upon the wide topic of discrete breather formation with a special emphasis on the the ϕ 4 model. We start by introducing the model and discussing some of the application areas/motivational aspects of exploring time periodic, spatially localized structures, such as the discrete breathers. Our main emphasis is on the existence, and especially on the stability features of such solutions. We explore their spectral stability numerically, as well as in special limits (such as the vicinity of the so-called anti-continuum limit of vanishing coupling) analytically. We also provide and explore a simple, yet powerful stability criterion involving the sign of the derivative of the energy vs. frequency dependence of such solutions. We then turn our attention to nonlinear stability, bringing forth the importance of a topological notion, namely the Krein signature. Furthermore, we briefly touch upon linearly and nonlinearly unstable dynamics of such states. Some special aspects/extensions of such structures are only touched upon, including moving breathers and dissipative variations of the model and some possibilities for future work are highlighted.

We consider an array of double oligomers in an optical waveguide device. A mathematical model for the system is the coupled discrete nonlinear Schrödinger (NLS) equations, where the gain-and-loss parameter contributes to the complex-valued linear coupling. The array caters to an optical simulation of the parity-time ( PT )-symmetry property between the coupled arms. The system admits fundamental bright discrete soliton solutions. We investigate their existence and spectral stability using perturbation theory analysis. These analytical findings are verified further numerically using the Newton-Raphson method and a standard eigenvalue-problem solver. Our study focuses on two natural discrete modes of the solitons: single- and double-excited-sites, also known as onsite and intersite modes, respectively. Each of these modes acquires three distinct configurations between the dimer arms, i.e., symmetric, asymmetric, and antisymmetric. Although both intersite and onsite discrete solitons are generally unstable, the latter can be stable, depending on the combined values of the propagation constant, horizontal linear coupling coefficient, and gain-loss parameter.

We study discrete solitons in zigzag discrete waveguide arrays with different types of linear mixing between nearest-neighbor and next-nearest-neighbor couplings. The waveguide array is constructed from two layers of one-dimensional (1D) waveguide arrays arranged in zigzag form. If we alternately label the number of waveguides between the two layers, the cross-layer couplings (which couple one waveguide in one layer with two adjacent waveguides in the other layer) construct the nearestneighbor couplings, while the couplings that couple this waveguide with the two nearest-neighbor waveguides in the same layer, i.e., self-layer couplings, contribute the next-nearest-neighbor couplings. Two families of discrete solitons are found when these couplings feature different types of linear mixing. As the total power is increased, a phase transition of the second kind occurs for discrete solitons in one type of setting, which is formed when the nearest-neighbor coupling and next-nearest-neighbor coupling feature positive and negative linear mixing, respectively. The mobilities and collisions of these two families of solitons are discussed systematically throughout the paper, revealing that the width of the soliton plays an important role in its

Two-dimensional arrays of nonlinear electric oscillators are considered theoretically, where nearest neighbors are coupled by relatively small, constant, but non-equal capacitors. The dynamics is approximately reduced to a weakly dissipative defocusing discrete nonlinear Schrödinger equation with translationally non-invariant linear dispersive coefficients. Behavior of quantized discrete vortices in such systems is shown to depend strongly on the spatial profile of the inter-node coupling as well as on the ratio between time-increasing healing length and lattice spacing. In particular, vortex clusters can be stably trapped for a some initial period of time by a circular barrier in the coupling profile, but then, due to gradual dissipative broadening of vortex cores, they lose stability and suddenly start to move.

We consider the one-dimensional dynamics of nonlinear non-dispersive waves. The problem can be mapped onto a linear one by means of the hodograph transform. We propose an approximate scheme for solving the corresponding Euler-Poisson equation which is valid for any kind of nonlinearity. The approach is exact for monoatomic classical gas and agrees very well with exact results and numerical simulations for other systems. We also provide a simple and accurate determination of the wave breaking time for typical initial conditions.

We investigate, both analytically and numerically, dispersive fractalization and quantization of solutions to periodic linear and nonlinear Fermi-Pasta-Ulam-Tsingou systems. When subject to periodic boundary conditions and discontinuous initial conditions, e.g., a step function, both the linearized and nonlinear continuum models for FPUT exhibit fractal solution profiles at irrational times (as determined by the coefficients and the length of the interval) and quantized profiles (piecewise constant or perturbations thereof) at rational times. We observe a similar effect in the linearized FPUT chain at times t where these models have validity, namely t=O( h −2 ) , where h is proportional to the intermass spacing or, equivalently, the reciprocal of the number of masses. For nonlinear periodic FPUT systems, our numerical results suggest a somewhat similar behavior in the presence of small nonlinearities, which disappears as the nonlinear force increases in magnitude. However, these phenomena are manifested on very long time intervals, posing a severe challenge for numerical integration as the number of masses increases. Even with the high-order splitting methods used here, our numerical investigations are limited to nonlinear FPUT chains with a smaller number of masses than would be needed to resolve this question unambiguously.

Long time dynamics of the smoothed step initial value problem or dispersive Riemann problem for the Benjamin-Bona-Mahony (BBM) equation u t +u u x = u xxt are studied using asymptotic methods and numerical simulations. The catalog of solutions of the dispersive Riemann problem for the BBM equation is much richer than for the related, integrable, Korteweg-de Vries equation u t +u u x + u xxx =0. The transition width of the initial smoothed step is found to significantly impact the dynamics. Narrow width gives rise to rarefaction and dispersive shock wave (DSW) solutions that are accompanied by the generation of two-phase linear wavetrains, solitary wave shedding, and expansion shocks. Both narrow and broad initial widths give rise to two-phase nonlinear wavetrains or DSW implosion and a new kind of dispersive Lax shock for symmetric data. The dispersive Lax shock is described by an approximate self-similar solution of the BBM equation whose limit as t?��? is a stationary, discontinuous weak solution. By introducing a slight asymmetry in the data for the dispersive Lax shock, the generation of an incoherent solitary wavetrain is observed. Further asymmetry leads to the DSW implosion regime that is effectively described by a pair of coupled nonlinear Schrödinger equations. The complex interplay between nonlocality, nonlinearity and dispersion in the BBM equation underlies the rich variety of nonclassical dispersive hydrodynamic solutions to the dispersive Riemann problem.

Dispersive shock waves (DSWs) in the three dimensional Benjamin- Ono (3DBO) equation is studied with step-like initial condition along a paraboloid front. By using a similarity reduction, problem of studying DSWs in three space one time (3+1) dimensions reduces to finding DSW solution of a (1+1) dimensional equation. By using a special ansatz, the 3DBO equation exactly reduces to the spherical Benjamin-Ono (sBO) equation. Whitham modulation equations are derived which describes DSW evolution in the sBO equation by using a perturbation method and these equations are written in terms of appropriate Riemmann type variables to obtain the sBO- Whitham system. DSW solution which obtained from the numerical solutions of the Whitham system and the direct numerical solution of the sBO equation are compared. In this comparison, a good agreement is found between these solutions. Also, some physical qualitative results about DSWs in sBO equation are presented. It is concluded that DSW solutions in the reduced sBO equation provide some information about DSW behaviour along the paraboloid fronts in the 3DBO equation.

Wave amplification in nonlinear dispersive wave equations may be caused by nonlinear focussing of waves from a certain background. In the model of nonlinear Schrödinger equation we will introduce a transformation to displaced phase-amplitude variables with respect to a background of monochromatic waves. The potential energy in the Hamiltonian then depends essentially on the phase. Looking as a special case to phases that are time independent, the oscillator equation for the signal at each position becomes autonomous, with the change of phase with position as only driving force for a spatial evolution towards extreme waves. This is observed to be the governing process of wave amplification in classes of already known solutions of NLS, namely the Akhemediev-, Ma- and Peregrine-solitons. We investigate the case of the soliton on finite background in detail in this Letter as the solution that describes the complete spatial evolution of modulational instability from background to extreme waves.