Journal of Nonlinear Science最新文献

For three specific singular perturbed three-component reaction–diffusion systems that admit N-spike solutions in one of the components on a finite domain, we present a detailed analysis for the dynamics of temporal oscillations in the spike positions. The onset of these oscillations is induced by N Hopf bifurcations with respect to the translation modes that are excited nearly simultaneously. To understand the dynamics of N spikes in the vicinity of Hopf bifurcations, we combine the center manifold reduction and the matched asymptotic method to derive a set of ordinary differential equations (ODEs) of dimension 2N describing the spikes’ locations and velocities, which can be recognized as normal forms of multiple Hopf bifurcations. The reduced ODE system then is represented in the form of linear oscillators with weakly nonlinear damping. By applying the multiple-time method, the leading order of the oscillation amplitudes is further characterized by an N-dimensional ODE system of the Stuart–Landau type. Although the leading order dynamics of these three systems are different, they have the same form after a suitable transformation. On the basis of the reduced systems for the oscillation amplitudes, we prove that there are at most (lfloor N/2 rfloor +1) stable equilibria, corresponding to (lfloor N /2 rfloor +1) types of different oscillations. This resolves an open problem proposed by Xie et al. (Nonlinearity 34(8):5708–5743, 2021) for a three-component Schnakenberg system and generalizes the results to two other classic systems. Numerical simulations are presented to verify the analytic results.

Coupled oscillator networks provide mathematical models for interacting periodic processes. If the coupling is weak, phase reduction—the reduction of the dynamics onto an invariant torus—captures the emergence of collective dynamical phenomena, such as synchronization. While a first-order approximation of the dynamics on the torus may be appropriate in some situations, higher-order phase reductions become necessary, for example, when the coupling strength increases. However, these are generally hard to compute and thus they have only been derived in special cases: This includes globally coupled Stuart–Landau oscillators, where the limit cycle of the uncoupled nonlinear oscillator is circular as the amplitude is independent of the phase. We go beyond this restriction and derive second-order phase reductions for coupled oscillators for arbitrary networks of coupled nonlinear oscillators with phase-dependent amplitude, a scenario more reminiscent of real-world oscillations. We analyze how the deformation of the limit cycle affects the stability of important dynamical states, such as full synchrony and splay states. By identifying higher-order phase interaction terms with hyperedges of a hypergraph, we obtain natural classes of coupled phase oscillator dynamics on hypergraphs that adequately capture the dynamics of coupled limit cycle oscillators.

Both smooth subsonic and transonic spiral flows to steady Euler–Poisson system with nonzero angular velocity and vorticity in a concentric cylinder are studied. On the one hand, we investigate the structural stability of smooth cylindrically symmetric subsonic flows under three-dimensional perturbations on the inner and outer cylinders. On the other hand, the structural stability of smooth transonic flows under the axi-symmetric perturbations is examined. There are no any restrictions on the background subsonic and transonic solutions. A deformation-curl-Poisson decomposition to the steady Euler–Poisson system is utilized to deal with the hyperbolic-elliptic mixed structure in the subsonic region. We emphasize that there is a special structure of the steady Euler–Poisson system which yields a priori estimates and uniqueness of the linearized elliptic system.

In this paper, the notion of the catenary curve in the sphere and in the hyperbolic plane is introduced. In both spaces, a catenary is defined as the shape of a hanging chain when its potential energy is determined by the distance to a given geodesic of the space. Several characterizations of the catenary are established in terms of the curvature of the curve and of the angle that its unit normal makes with a vector field of the ambient space. Furthermore, in the hyperbolic plane, we extend the concept of catenary substituting the reference geodesic by a horocycle or the hyperbolic distance by the horocycle distance.

Predicting the evolution of dynamics from a given trajectory history of an unknown system is an important and challenging problem. This paper presents a model-free method of forecasting unknown chaotic systems through reconstructing vector fields from noisy measured data via an adaptation of optimal control methods. This technique is also applicable to partially observed systems using a Takens delay embedding approach. The algorithms are validated on the Lorenz system and the four-dimensional hyperchaotic Rössler system, and demonstrate successful predictions well beyond the Lyapunov timescale. It is found that for small datasets or datasets with large levels of noise, the prediction accuracy of partially observed systems approaches that of fully observed systems. The presented approach also allows the model-free assessment of local predictability on the attractor by evolving initial condition density through the reconstructed vector fields via estimation of the transfer operator. The method is compared to predictions made by an imperfect model which highlights the utility of model-free approaches when the only available models have significant model error. The capability of this method for reconstruction of continuous and global vector fields may be applied to model validation, forecasting of initial conditions not in the training set, and model-free filtering.

The study of nonlinear Schrödinger-type equations with nonzero boundary conditions introduces challenging problems both for the continuous (partial differential equation) and the discrete (lattice) counterparts. They are associated with fascinating dynamics emerging by the ubiquitous phenomenon of modulation instability. In this work, we consider the discrete nonlinear Schrödinger equation with linear gain and nonlinear loss. For the infinite lattice supplemented with nonzero boundary conditions, which describe solutions decaying on the top of a finite background, we provide a rigorous proof that for the corresponding initial boundary value problem, solutions exist for any initial condition, if and only if, the amplitude of the background has a precise value (A_*) defined by the gain-loss parameters. We argue that this essential property of this infinite lattice cannot be captured by finite lattice approximations of the problem. Commonly, such approximations are provided by lattices with periodic boundary conditions or as it is shown herein, by a modified problem closed with Dirichlet boundary conditions. For the finite-dimensional dynamical system defined by the periodic lattice, the dynamics for all initial conditions are captured by a global attractor. Analytical arguments corroborated by numerical simulations show that the global attractor is trivial, defined by a plane wave of amplitude (A_*). Thus, any instability effects or localized phenomena simulated by the finite system can be only transient prior the convergence to this trivial attractor. Aiming to simulate the dynamics of the infinite lattice as accurately as possible, we study the dynamics of localized initial conditions on the constant background and investigate the potential impact of the global asymptotic stability of the background with amplitude (A_*) in the long-time evolution of the system.

The occurrence of finite time blow-up phenomenon in the Keller–Segel (KS) model has always been a significant area of interest for mathematicians. Despite extensive research on the blow-up phenomenon in KS models with signal production, Understanding of this phenomenon in models with signal consumption mechanisms has been scarce.This paper marks a preliminary investigation into this unexplored field. In this study, we employ a backward self-similar solution to demonstrate that the finite time blowup indeed occurs within this model. More precisely, in one-dimensional space, finite time blowing up corresponding to the chemotactic collapse phenomenon (the formation of Dirac (delta )-singularity ) happens; in high-dimensional space, the self-similar solution will blow up everywhere. Finally, we also consider the special cases where the diffusion coefficient of bacteria or oxygen is 0. For these cases, chemotactic collapse phenomenon occurs in both one-dimensional and two-dimensional spaces.

We consider piecewise smooth vector fields (Z=(Z_+, Z_-)) defined in ({mathbb {R}}^n) where both vector fields are tangent to the switching manifold (Sigma ) along a submanifold (Msubset Sigma ). We shall see that, under suitable assumptions, Filippov convention gives rise to a unique sliding mode on M, governed by what we call the tangential sliding vector field. Here, we will provide the necessary and sufficient conditions for characterizing such a vector field. Additionally, we prove that the tangential sliding vector field is conjugated to the reduced dynamics of a singular perturbation problem arising from the Sotomayor–Teixeira regularization of Z around M. Finally, we analyze several examples where tangential sliding vector fields can be observed, including a model for intermittent treatment of HIV.

In this paper, we study the existence of spread speeds and periodic traveling waves for a class of space–time periodic and degenerate cooperative systems with nonlocal dispersal and delays. In order to characterize the spread speeds, we first establish the principal eigenvalue theory for linear space–time periodic and degenerate systems with nonlocal dispersal and delays and give some sufficient conditions for the existence of the principal eigenvalue. Then, we prove the existence of a single spreading speed and give the computational formulae of it. Finally, we generalize the monotone iteration scheme combined with the method of sub–super-solutions to prove the existence of the space–time periodic traveling waves.

We develop a convex integration scheme for constructing nonunique weak solutions to the hydrostatic Euler equations (also known as the inviscid primitive equations of oceanic and atmospheric dynamics) in both two and three dimensions. We also develop such a scheme for the construction of nonunique weak solutions to the three-dimensional viscous primitive equations, as well as the two-dimensional Prandtl equations. While in Boutros et al. (Calc Var Partial Differ Equ 62(8):219, 2023) the classical notion of weak solution to the hydrostatic Euler equations was generalised, we introduce here a further generalisation. For such generalised weak solutions, we show the existence and nonuniqueness for a large class of initial data. Moreover, we construct infinitely many examples of generalised weak solutions which do not conserve energy. The barotropic and baroclinic modes of solutions to the hydrostatic Euler equations (which are the average and the fluctuation of the horizontal velocity in the z-coordinate, respectively) that are constructed have different regularities.