This study investigates a detailed analytical and numerical investigation of a nonlinear two-degree-of-freedom (2DOF) mechanical oscillator subjected to parametric excitation, magnetic stiffness nonlinearities, and dry friction. The considered system consists of two coupled oscillators, both of which are connected to a rotating rectangular beam that induces a time-periodic stiffness variation. The Complex Averaging (CxA) method is employed to derive approximate analytical solutions, which are thoroughly validated through time-domain simulations and bifurcation analyses. The dynamic analysis reveals a rich spectrum of nonlinear behaviors, including periodic, quasi-periodic, and chaotic responses. Detailed bifurcation diagrams, Lyapunov exponent analysis, and Poincaré maps demonstrate the influence of nonlinear stiffness degree, mass symmetry, and frictional effects on system stability and response amplitude. The obtained results give a significant understanding of the dynamic behavior of coupled nonlinear systems and establish a conceptual framework for the development of complex vibration abatement strategies, energy harvesting devices, and advanced mechanical systems.
Gas turbine engines are complex and highly nonlinear dynamical systems. Deriving their physics-based models can be challenging because it requires performance characteristics that are not always available, often leading to many simplifying assumptions. This paper discusses the limitations of conventional experimental methods used to derive component-level and locally linear parameter-varying models, and addresses these issues by employing identification techniques based on data collected from standard engine operation under closed-loop control. The rotor dynamics are estimated using the sparse identification of nonlinear dynamics. Subsequently, the autonomous part of the dynamics is mapped into an optimally constructed Koopman eigenfunction space. This process involves eigenvalue optimization using metaheuristic algorithms and temporal projection, followed by gradient-based eigenfunction identification. The resulting Koopman model is validated against an in-house reference component-level model. A globally optimal nonlinear feedback controller and a Kalman estimator are then designed within the eigenfunction space and compared to traditional and gain-scheduled proportional-integral
This paper investigates density driven flow in porous media, focusing on the roles of viscosity contrast, density contrast, and linear adsorption. In this setup, the fluid on top is heavier and more viscous than the fluid below. Under the effect of gravity, this system becomes unstable, and finger-like structures appear. The phenomenon is described mathematically by coupling Darcy's law with a convection-diffusion reaction equation. The nonlinearity in this model arises mainly from the concentration dependence of viscosity and the convective transport term. The existence of a unique pair of weak solutions is shown using the Galerkin approximation method and truncation technique. Moreover, an application of the maximum principle shows non-negativity of the concentration. Additionally, we analyze the long-time behavior of the solution and prove that the concentration converges exponentially to zero in the $L^p$-norm for all $1 \le p \le \infty$ as $t \to \infty.$ To complement the theoretical analysis, we perform numerical simulations based on a pressure formulation. By tracking total kinetic energy and mixing measures over time, we discuss the instability and the mixing efficiency,
We investigate the asymmetric integrable turbulence and rogue waves (RWs) emerging from the modulation instability (MI) of plane waves for the DNLS equation. The \(n\)-th moments and ensemble-averaged kinetic and potential energy exhibit oscillatory convergence towards their steady-state values. Specifically, the amplitudes of oscillations for these indexes decay asymptotically with time as \(t^{-1.36}\), while the phase shifts demonstrate a nonlinear decay with a rate of \(t^{-0.78}\). The frequency of these oscillations is observed to be twice the maximum growth rate of MI. These oscillations can be classified into two distinct types: one is in phase with ensemble-averaged potential energy modulus $|\langle H_4\rangle|$, and the other is anti-phase. At the same time, this unity is also reflected in the wave-action spectrum \( S_k(t) \) for a given \( k \), the auto-correlation function \( g(x,t) \) for a given \( x \), as well as the PDF \( P(I,t) \). The critical feature of the turbulence is the wave-action spectrum, which follows a power-law distribution of \( |k+3|^{-α} \) expect for $k=-3$. Unlike the NLS equation, the turbulence in the DNLS setting is asymmetric, primarily d
We demonstrate transmission stabilization against radiation emission by frequency-dependent linear gain-loss and perturbation-induced frequency shifting for solitons of the cubic-quintic nonlinear Schrödinger (CQNLS) equation. We consider soliton propagation in a nonlinear optical waveguide with focusing cubic nonlinearity, defocusing quintic nonlinearity, and dissipative perturbations due to weak frequency-dependent linear gain-loss, cubic loss, and delayed Raman response. The frequency shifting is induced by delayed Raman response. Our perturbation analysis and numerical simulations with the perturbed CQNLS equation show that transmission stabilization with CQNLS solitons is indeed possible, and in this way provide the first demonstration of the stabilization method for solitons of a nonintegrable nonlinear wave model. Moreover, we find that transmission stabilization with energetic CQNLS solitons is realized with significantly smaller frequency shifts and pulse distortion compared with stabilization with energetic solitons of the cubic nonlinear Schrödinger equation. Therefore, our study also demonstrates that stabilization of energetic solitons by the method is significantly im
Conservation laws are computed for various nonlinear partial differential equations that arise in elasticity and acoustics. Using a scaling homogeneity approach, conservation laws are established for two models describing shear wave propagation in a circular cylinder and a cylindrical annulus. Next, using the multiplier method, conservation laws are derived for a parameterized system of constitutive equations in cylindrical coordinates involving a general expression for the Cauchy stress. Conservation laws for the Khokhlov-Zabolotskaya-Kuznetsov equation and Westervelt-type equations in various coordinate systems are also presented.
This article extends the study of the dynamical properties of the symmetric McMillan map, emphasizing its utility in understanding and modeling complex nonlinear systems. Although the map features six parameters, we demonstrate that only two are irreducible: the linearized rotation number at the fixed point and a nonlinear parameter representing the ratio of terms in the biquadratic invariant. Through a detailed analysis, we classify regimes of stable motion, provide exact solutions to the mapping equations, and derive a canonical set of action-angle variables, offering analytical expressions for the rotation number and nonlinear tune shift. We further establish connections between general standard-form mappings and the symmetric McMillan map, using the area-preserving Hénon map and accelerator lattices with thin sextupole magnet as representative case studies. Our results show that, despite being a second-order approximation, the symmetric McMillan map provides a highly accurate depiction of dynamics across a wide range of system parameters, demonstrating its practical relevance in both theoretical and applied contexts.
Nonlinear science has evolved significantly over the 35 years since the launch of the journal Chaos. This Focus Issue, dedicated to the 80th Birthday of its founding editor-in-chief, David K. Campbell, brings together a selection of contributions on influential topics, many of which were advanced by Campbell's own research program and leadership role. The topics include new phenomena and method development in the realms of network dynamics, machine learning, quantum and material systems, chaos and fractals, localized states, and living systems, with a good balance of literature review, original contributions, and perspectives for future research.
We establish the necessary and sufficient conditions for constructing gradient Einstein-type warped metrics. One of these conditions leads us to a general Lichnerowicz equation with analytic and geometric coefficients for this class of metrics on the space of warping functions. In this way, we prove gradient estimates for positive solutions of a nonlinear elliptic differential equation on a complete Riemannian manifold with associated Bakry-Émery Ricci tensor bounded from below. As an application, we provide nonexistence and rigidity results for a large class of gradient Einstein-type warped metrics. Furthermore, we show how to construct gradient Einstein-type warped metrics, and then we give explicit examples which are not only meaningful in their own right, but also help to justify our results.
Quantifying the entanglement between two continuous bosonic modes, such as magnons and photons, is not trivial. The logarithmic negativity, calculated through the quantum Langevin equations is subjected to thermal noise. However this method requires further approximation. The phase space of a generic nonlinear system contains topologically different regions, and the steady state may correspond to the different types of fixed points, such as Saddle Points, Stable or unstable Spirals, and Nodes. In the present work, we propose a new procedure. Namely, we derived the complete set of nonlinear equations, which includes equations for the magnon and photon number operators and phases. We show that not only number operators but also phases are important for exploring the character of the fixed point, and magnon-photon entanglement. We showed that methods of the qualitative theory of nonlinear differential equations are also relevant for photon-magnon entanglement problems. Our main finding is that entanglement is not defined in the Saddle Point region. On the other hand, the maximum of the entanglement corresponds to the region near the border between the Stable node and Stable spiral reg
We consider reduced-order modeling of nonlinear dispersive waves described by a class of nonlinear Schrodinger (NLS) equations. We compare two nonlinear reduced-order modeling methods: (i) The reduced Lagrangian approach which relies on the variational formulation of NLS and (ii) The recently developed method of reduced-order nonlinear solutions (RONS). First, we prove the surprising result that, although the two methods are seemingly quite different, they can be obtained from the real and imaginary parts of a single complex-valued master equation. Furthermore, for the NLS equation in a stationary frame, we show that the reduced Lagrangian method fails to predict the correct group velocity of the waves whereas RONS predicts the correct group velocity. Finally, for the modified NLS equation, where the reduced Lagrangian approach is inapplicable, the RONS reduced-order model accurately approximates the true solutions.
It is shown that nonlinear electrodynamics of the Born--Infeld theory type may be exploited to shed insight into a few fundamental problems in theoretical physics, including rendering electromagnetic asymmetry to energetically exclude magnetic monopoles, achieving finite electromagnetic energy to relegate curvature singularities of charged black holes, and providing theoretical interpretation of equations of state of cosmic fluids via k-essence cosmology. Also discussed are some nonlinear differential equation problems.
We consider the Cauchy problem for nonlinear Schrödinger equations in a general domain $Ω\subset\mathbb{R}^N$. Construction of solutions has been only done by classical compactness method in previous results. Here, we construct solutions by a simple alternative approach. More precisely, solutions are constructed by proving that approximate solutions form a Cauchy sequence in some Banach space. We discuss three different types of nonlinearities: power type nonlinearities, logarithmic nonlinearities and damping nonlinearities.
The suggestion of writing, for some problems, nonlinear state equations not as dx/dt = F(x,u,t), but as dx/dt = [A(t,x)]x + [B(t,x)]u(t), which is more "constructive", is considered supported by arguments related to: the axiomatization of system theory, the classification of switched circuits as linear and nonlinear, the use of nonlinear sampling for measuring frequency without clock, the consideration of an ensemble of colliding particles in which the establishing of the chaotic movements is seen as a result of the same kind of nonlinearity as in electronic switching circuits, the modeling of a liquid medium as a "system" whose structure is directly influenced by the "input", and some other problems.
Overcontact binary stars are systems of two stars where the component stars are in contact with each other. This implies that they share a common envelope of gas. In this work we seek signatures of nonlinearity and chaos in these stars by using time series analysis techniques. We use three main techniques, namely the correlation dimension,f (α) spectrum and the bicoherence. The former two are calculated from the reconstructed dynamics, while the latter is calculated from the Fourier transforms of the time series of intensity variations(light curves) of these stars. Our dataset consists of data from 463 overcontact binary stars in the Kepler field of view [1]. Our analysis indicates nonlinearity and signatures of chaos in almost all the light curves. We also explore whether the underlying nonlinear properties of the stars are related to their physical properties like fill-out-factor, a measure of the extend of contact between the components of an overcontact binary system . We observe that significant correlations exist between the fill out factor and the nonlinear quantifiers. This correlation is more pronounced in specific subcategories constructed based on the mass ratios and eff
In this set of papers we formulate a stand alone method to derive maximal number of linearizing transformations for nonlinear ordinary differential equations (ODEs) of any order including coupled ones from a knowledge of fewer number of integrals of motion. The proposed algorithm is simple, straightforward and efficient and helps to unearth several new types of linearizing transformations besides the known ones in the literature. To make our studies systematic we divide our analysis into two parts. In the first part we confine our investigations to the scalar ODEs and in the second part we focuss our attention on a system of two coupled second order ODEs. In the case of scalar ODEs, we consider second and third order nonlinear ODEs in detail and discuss the method of deriving maximal number of linearizing transformations irrespective of whether it is local or nonlocal type and illustrate the underlying theory with suitable examples. As a by-product of this investigation we unearth a new type of linearizing transformation in third order nonlinear ODEs. Finally the study is extended to the case of general scalar ODEs. We then move on to the study of two coupled second order nonlinear
Let $0<q<1<p$. In this study, we investigate positive solutions of the logistic elliptic equation $-Δu = u(1-u^{p-1})$ in a smooth bounded domain $Ω$ of $\mathbb{R}^N$, $N\geq1$, with the nonlinear boundary condition $\frac{\partial u}{\partial ν}=-λu^q$ on $\partialΩ$. This nonlinear boundary condition arises from coastal fishery harvesting. When $p>1$ is subcritical, we prove that in the case of $λ_Ω>1$, there exist at least two positive solutions for $λ>0$ sufficiently small but no positive solutions for $λ>0$ large enough. In the case of $λ_Ω<1$, there exists at least one positive solution for every $λ>0$. Here, $λ_Ω>0$ is the smallest eigenvalue of $-Δ$ under the Dirichlet boundary condition. An interpretation of our main results from an ecological viewpoint is presented.
A bistable nonlinear energy sink conceived to mitigate the vibrations of host structural systems is considered in this paper. The hosting structure consists of two coupled symmetric linear oscillators (LOs) and the nonlinear energy sink (NES) is connected to one of them. The peculiar nonlinear dynamics of the resulting three-degree-of-freedom system is analytically described by means of its slow invariant manifold derived from a suitable rescaling, coupled with a harmonic balance procedure, applied to the governing equations transformed in modal coordinates. On the basis of the first-order reduced model, the absorber is tuned and optimized to mitigate both modes for a broad range of impulsive load magnitudes applied to the LOs. On the one hand, for low-amplitude, in-well, oscillations, the parameters governing the bistable NES are tuned in order to make it functioning as a linear tuned mass damper (TMD); on the other, for high-amplitude, cross-well, oscillations, the absorber is optimized on the basis of the invariant manifolds features. The analytically predicted performance of the resulting tuned bistable nonlinear energy sink (TBNES) are numerically validated in terms of dissipa
We present a method of deriving linearizing transformations for a class of second order nonlinear ordinary differential equations. We construct a general form of a nonlinear ordinary differential equation that admits Bernoulli equation as its first integral. We extract conditions for this integral to yield three different linearizing transformations, namely point, Sundman and generalized linearizing transformations. The explicit forms of these linearizing transformations are given. The exact forms and the general solution of the nonlinear ODE for these three linearizables cases are also enumerated. We illustrate the procedure with three different examples.
We describe a method to model nonlinear dynamical systems using periodic solutions of delay-differential equations. We show that any finite-time trajectory of a nonlinear dynamical system can be loaded approximately into the initial condition of a linear delay-differential system. It is further shown that the initial condition can be extended to a periodic solution of the delay-differential system if an appropriate choice of its parameters is made. As a result, any finite set of trajectories of a nonlinear dynamical system can be modeled with arbitrarily small error via a set of periodic solutions of a linear delay-differential equation. These results can be extended to some non-linear delay differential systems. One application of the method is for modeling memory and perception.