Nonlinear Solutions Research Articles

Acceleration of ion phase-space holes due to interactions with ion solitons in a wave-guided plasma

Ion phase-space holes are solitary kinetic structures found in the ion phase-space of collision-less plasmas, and are nonlinear solutions to the Vlasov-Poisson equations, identified as Bernstein-Greene-Kruskal (BGK) modes. In this study, interactions between an ion phase-space hole and a travelling ion KdV soliton is presented. This interaction, which is simulated in a fully ionised highly magnetised plasma within a cylindrical wave-guide, exhibits acceleration and deacceleration of the ion hole, depending on its mode of collision with the travelling ion soliton. We present these interactions and discuss the mechanism of this interaction between the two solitary waves.

Testing a Nonlinear Solution of the Israel–Stewart Theory

Testing a Nonlinear Solution of the Israel–Stewart Theory

Modon solutions in an N-layer quasi-geostrophic model

Modons, or dipolar vortices, are common and long-lived features of the upper ocean, consisting of a pair of counter-rotating monopolar vortices moving through self-advection. Such structures remain stable over long times and may be important for fluid transport over large distances. Here, we present a semi-analytical method for finding fully nonlinear modon solutions in a multi-layer quasi-geostrophic model with arbitrarily many layers. Our approach is to reduce the problem to a multi-parameter linear eigenvalue problem which can be solved using numerical techniques from linear algebra. The method is shown to replicate previous results for one- and two-layer models and is applied to a three-layer model to find a solution describing a mid-depth propagating, topographic vortex.

The localized excitation on the Weierstrass elliptic function periodic background for the (3+1)-dimensional generalized variable-coefficient Kadomtsev-Petviashvili equation

In this paper, the linear spectral problem associated with the (3+1)-dimensional generalized variable-coefficient Kadomtsev-Petviashvili (gvcKP) equation with the Weierstrass function as the external potential is investigated based on the Lamé function, from which some new localized nonlinear wave solutions on the Weierstrass elliptic ℘-function periodic background are obtained by the Darboux transformation. The degenerate solutions on the ℘-function periodic background for the gvcKP equation can be derived by taking the limits of the half-periods ω 1, ω 2 of ℘(x), whose evolution and corresponding dynamics are also discussed. The findings show that nonlinear waves on the ℘-function periodic background behave as different types of nonlinear waves in different spaces, including periodic waves, vortex solitons and interaction solutions, aiding in elucidating some physical phenomena in the related fields, such as the physical ocean and nonlinear optics.

An analytical thermal model for the optimization of EU DEMO feeder thermal shields

The present work aims at defining a first design of the DEMO EU’s feeder thermal shield for the CFT, by exploring the temperature distribution resulting from radiation and conduction coming from different sources including the vacuum duct, the containment duct and the cryogenic He lines on the thermal shield. In order to achieve these results, at first a numerical FEM representation was produced and parameterized. To ease the burden of such a complex and highly non-linear solution, an analytical model however was created and validated against the numerical one with very good results. The sensitivity of several parameters was computed with respect to the temperature distribution of the thermal shield and the heat flow to the containment duct, finding a candidate configuration for the thermal shield. To complete the study, an evolutive optimization based on gradient descent was implemented to understand the influence of the circumferential positioning of the cooling He lines.

Nonlinear dynamic analysis of an inclined micro-beam under a moving mass

Nonlinear vibration of an inclined simply supported micro-beam under a moving mass is investigated for Euler–Bernoulli beam theory (EBT) and Timoshenko beam theory (TBT) respectively. Based on a modified couple stress theory (MCST) and the von-Karman geometric nonlinearity, the nonlinear coupled dynamic equations of the system are established through the Hamilton’s principle with the assumed mode method. A wide range of numerical examples are employed to study the influence of slenderness ratio, cross-section height, inclined angle, the size and velocity of the moving mass and the scale factor of the material on the solutions of nonlinear and linear, the solutions of EBT and TBT and the solutions of moving mass and moving load. By comparing the differences between the nonlinear and linear solutions under different parameters and beam theories, the importance and significance of nonlinear dynamic analysis of the inclined micro-beam are revealed.

Robust recursive sigma point Kalman filtering for Huber-based generalized M-estimation

For nonlinear state estimation driven by non-Gaussian noise, the estimator is required to be updated iteratively. Since the iterative update approximates a linear process, it fails to capture the nonlinearity of observation models, and this further degrades filtering accuracy and consistency. Given the flaws of nonlinear iteration, this work incorporates a recursive strategy into generalized M-estimation rather than the iterative strategy. The proposed algorithm extends nonlinear recursion to nonlinear systems using the statistical linear regression method. The recursion allows for the gradual release of observation information and consequently enables the update to proceed along the nonlinear direction. Considering the correlated state and observation noise induced by recursions, a separately reweighting strategy is adopted to build a robust nonlinear system. Analogous to the nonlinear recursion, a robust nonlinear recursive update strategy is proposed, where the associated covariances and the observation noise statistics are updated recursively to ensure the consistency of observation noise statistics, thereby completing the nonlinear solution of the robust system. Compared with the iterative update strategies under non-Gaussian observation noise, the recursive update strategy can facilitate the estimator to achieve higher filtering accuracy, stronger robustness, and better consistency. Therefore, the proposed strategy is more suitable for the robust nonlinear filtering framework.

Envelope vector solitons in nonlinear flexible mechanical metamaterials

In this paper, we employ a combination of analytical and numerical techniques to investigate the dynamics of lattice envelope vector soliton solutions propagating within a one-dimensional chain of flexible mechanical metamaterial. To model the system, we formulate discrete equations that describe the longitudinal and rotational displacements of each individual rigid unit mass using a lump element approach. By applying the multiple-scales method in the context of a semi-discrete approximation, we derive an effective nonlinear Schrödinger equation that characterizes the evolution of rotational and slowly varying envelope waves from the aforementioned discrete motion equations. We thus show that this flexible mechanical metamaterial chain supports envelope vector solitons where the rotational component has the form of either a bright or a dark soliton. In addition, due to nonlinear coupling, the longitudinal displacement displays kink-like profiles thus forming the 2-components vector soliton. These findings, which include specific vector envelope solutions, enrich our knowledge on the nonlinear wave solutions supported by flexible mechanical metamaterials and open new possibilities for the control of nonlinear waves and vibrations.

On the stability of fully nonlinear hydraulic-fall solutions to the forced water wave problem

Two-dimensional free-surface flow over localised topography is examined, with the emphasis on the stability of hydraulic-fall solutions. A Gaussian topography profile is assumed with a positive or negative amplitude modelling a bump or a dip, respectively. Steady hydraulic-fall solutions to the full incompressible, irrotational Euler equations are computed, and their linear and nonlinear stability is analysed by computing eigenspectra of the pertinent linearised operator and by solving an initial value problem. The computations are carried out numerically using a specially developed computational framework based on the finite-element method. The Hamiltonian structure of the problem is demonstrated, and stability is determined by computing eigenspectra of the pertinent linearised operator. It is found that a hydraulic-fall flow over a bump is spectrally stable. The corresponding flow over a dip is found to be linearly unstable. In the latter case, time-dependent simulations show that ultimately, the flow settles into a time-periodic motion that corresponds to an invariant solution in an appropriately defined phase space. Physically, the solution consists of a localised large-amplitude wave that pulsates above the dip while simultaneously emitting nonlinear cnoidal waves in the upstream direction and multi-harmonic linear waves in the downstream direction.

A scalar field theory of 1+1-dimensional laser wakefield accelerators

Abstract A relativistic non-linear scalar field theory is developed from a 2+2-dimensional decomposition of the cold plasma field equations, and the theory is used to investigate a 1+1-dimensional description of a laser wakefield accelerator. The relationship between the properties of a compact laser pulse and its wake is explored. Non-linear solutions are sought describing a regular (i.e. unbroken) wake driven by a prescribed rectangular circularly-polarised laser pulse. An upper bound on the dimensionless amplitude a 0 of the laser pulse is determined as a function of the phase speed v of the wake. The asymptotic behaviour of the upper bound on a 0 as v → c is shown to agree with well-established, but approximate, results obtained using the conventional encoding of the plasma degrees of freedom. Our approach leads to a closed-form expression for the upper bound on a 0 which is exact for all values of the phase speed of the wake, unlike conventional results that are applicable only when v is sufficiently close to c.

A Novel Approach for Material Handling-Driven Facility Layout

Material handling is a widely used process in manufacturing and is generally considered a non-value-added process. The Dynamic Facility Layout Problem (DFLP) considered in this paper minimizes the total material handling and re-arrangement cost. In this study, an integrated DFLP model with unequal facility areas, assignment of material handling devices (MHD), and flexible bay structure (FBS) is considered, and it is aimed to propose fast solution approaches. Two different solution methods are proposed for the problem, which are the genetic algorithm and the simulated annealing algorithm, respectively. In both methods, a non-linear mathematical model solution was used to calculate the fitness values. Thus, the solutions in the feasible solution space are utilized. The proposed solution approaches were applied to solve four problems published in the literature. The computational experiments have validated the effectiveness of the algorithms and the quality of solutions produced.

Nonlinear Poisson-Boltzmann solutions for charged parallel plates: When opposite charges repel.

I present an exact solution of the Poisson-Boltzmann equation for two parallel plates and discuss the solution properties. I discuss in more detail plates with opposite charges: In this case, there are two critical separations, Lc,1 < Lc,2. For separations less than Lc,1, the force between plates is repulsive. It switches to attractive at Lc,1, but with the electric potential having the same signon both plates. For L > Lc,2, the force remains attractive, and the potential at the plates has the same signas the charge on each plate. I also describe charge regulation, determined by pKa, and provide formulas for both the critical distance where oppositely charged plates repel and their charging process. The implications of these results for the nanoparticle assembly, as driven by electrostatic interactions, are also discussed.

Local spherical collapsing box in Athena++ : Numerical implementation and benchmark tests

We implement a local model for a spherical collapsing or expanding gas cloud in the Athena++ magnetohydrodynamic code. This local model consists of a Cartesian periodic box with time-dependent geometry. We present a series of benchmark test problems, including nonlinear solutions and linear perturbations of the local model, confirming the code's desired performance. During a spherical collapse, a horizontal shear flow is amplified, corresponding to angular momentum conservation of zonal flows in the global problem; wave speed and the amplitude of sound waves increase in the local frame, due to the reduction in the characteristic length scale of the box, which can lead to an anisotropic effective sound speed in the local box. Our code conserves both mass and momentum-to-machine precision. This numerical implementation of the local model has potential applications to the study of local physics and hydrodynamic instabilities during protostellar collapse, providing a powerful framework for better understanding the earliest stages of star and planet formation.

Theoretical solution of vortex splitting due to islands effect on the vortex

Islands, acting as transitional areas between land and sea, significantly influence their surrounding environments. Studying the changes in vortex structure resulting from the interaction between islands and mesoscale vortices is crucial for understanding the dynamic characteristics and ecological processes of the marine environment. Previous theoretical studies have shown that in a domain without boundaries, due to the conservation of angular momentum, a vortex cannot split by itself. This paper establishes the conditions for the splitting of an anticyclonic vortex when it collides with two square islands. By linking the initial and final states of islands effect on the vortex and utilizing conserved quantities such as angular momentum and mass, along with the slow variation approximation, a nonlinear theoretical solution is constructed. The analysis shows that when the interaction between the vortex and the two islands leads to vortex splitting, the length of the islands, the vortex radius, and the distance between the islands must satisfy a certain condition O(LR)∼O(wR)∼O(1). These results provide support for subsequent analyses of the impact of various parameters on vortex structure when the North Brazil Current (NBC) ring encounters the Lesser Antilles in the tropical western Atlantic.

Wronskian solution, Bäcklund transformation and Painlevé analysis to a (2 + 1)-dimensional Konopelchenko–Dubrovsky equation

Abstract The main work of this paper is to construct the Wronskian solution and investigate the integrability characteristics of the (2 + 1)-dimensional Konopelchenko–Dubrovsky equation. Firstly, the Wronskian technique is used to acquire a sufficient condition of the Wronskian solution. According to the Wronskian form, the soliton solution is obtained by selecting the elements in the determinant that satisfy the linear partial differential systems. Secondly, the bilinear Bäcklund transformation and Bell-polynomial-typed Bäcklund transformation are derived directly via the Hirota bilinear method and the Bell polynomial theory, respectively. Finally, Painlevé analysis proves that this equation possesses the Painlevé property, and a Painlevé-typed Bäcklund transformation is constructed to solve a family of exact solutions by selecting appropriate seed solution. It shows that the Wronskian technique, Bäcklund transformation, Bell polynomial and Painlevé analysis are applicable to obtain the exact solutions of the nonlinear evolution equations, e.g., soliton solution, single-wave solution and two-wave solution.

Complex Ginzburg–Landau equation for time‐varying anisotropic media

AbstractWhen extending the complex Ginzburg–Landau equation (CGLE) to more than one spatial dimension, there is an underlying question of whether one is capturing all the interesting physics inherent in these higher dimensions. Although spatial anisotropy is far less studied than its isotropic counterpart, anisotropy is fundamental in applications to superconductors, plasma physics, and geology, to name just a few examples. We first formulate the CGLE on anisotropic, time‐varying media, with this time variation permitting a degree of control of the anisotropy over time, focusing on how time‐varying anisotropy influences diffusion and dispersion within both bounded and unbounded space domains. From here, we construct a variety of exact dissipative nonlinear wave solutions, including analogs of wavetrains, solitons, breathers, and rogue waves, before outlining the construction of more general solutions via a dissipative, nonautonomous generalization of the variational method. We finally consider the problem of modulational instability within anisotropic, time‐varying media, obtaining generalizations to the Benjamin–Feir instability mechanism. We apply this framework to study the emergence and control of anisotropic spatiotemporal chaos in rectangular and curved domains. Our theoretical framework and specific solutions all point to time‐varying anisotropy being a potentially valuable feature for the manipulation and control of waves in anisotropic media.

The localized excitation on the Jacobi elliptic function periodic background for the Gross–Pitaevskii equation

In this paper, the nonlinear wave solutions for Gross–Pitaevskii equation on the periodic wave background are investigated by Darboux-Bäcklund transformation, from which the soliton and breather wave solutions on the Jacobi elliptic cn and dn functions backgrounds are derived. The corresponding evolutions and dynamical properties of nonlinear wave solutions under different parameters are discussed. These results reported in this paper may raise the possibility of related experiments and potential applications in nonlinear science fields, such as nonlinear optics, oceanography and so on.

Analyzing bifurcation, stability, and wave solutions in nonlinear telecommunications models using transmission lines, Hamiltonian and Jacobian techniques

This study presents a comprehensive analysis of a nonlinear telecommunications model, exploring bifurcation, stability, and wave solutions using Hamiltonian and Jacobian techniques. The investigation begins with a thorough examination of bifurcation behavior, identifying critical points and their stability characteristics, leading to the discovery of diverse bifurcation scenarios. The stability of critical points is further assessed through graphical and numerical methods, highlighting the sensitivity to parameter variations. The study delves into the derivation of both numerical and analytical wave solutions, aligning them with energy orbits depicted in phase portraits, revealing a spectrum of wave behaviors. Additionally, the analysis extends to traveling wave solutions, providing insights into wave propagation dynamics. Notably, the study underscores the efficacy of the planar dynamical approach in capturing system behavior in harmony with phase portrait orbits. The findings have significant implications for telecommunications engineers and researchers, offering insights into system behavior, stability, and signal propagation, ultimately advancing our understanding of complex nonlinear dynamics in telecommunications networks.

Nonlinear Analysis of Electric Potential in a Piezoelectric Thin Plate as Pressure Sensor

Piezoelectric materials are widely applied in electronic components due to their ability to couple electric and mechanical fields. For the piezoelectric thin plate as a pressure sensor, we extend the modified first-order piezoelectric plate theory to explore the nonlinear solution of electric potential and discuss the relevant regulation mechanisms. Based on the piezoelectric effect, the electric potential generated by deformation is not linear with thickness. A quadratic function of electric potential across thickness is introduced. However, the electric potential on upper and lower electrodes is an unknown constant, which can be determined by the nonlinear complementary equation based on the Gauss theorem. It is shown that the nonlinear results of electric potential are consistent with outcomes by the three-dimensional finite element method. Furthermore, relevant regulation mechanisms are discussed on electric potential, including supported conditions, plate thickness, and load area. In addition, the influence of load locations on the electric potential response is analyzed. It is expected that these findings will be instructive to the structural design of piezoelectric pressure sensors.

Frequency Domain Thermal Analysis of Future Spacecraft and Systems

Future space telescopes require tight thermomechanical stability. The thermal analysis may be based on a nonlinear steady solution with a superimposed linear transient perturbation. This linearity allows frequency domain (FD) thermal solutions that retain high resolution when scaled. Linearity allows generality through superposition. Linear superelements can represent subsystems and thereby simplify interfaces. The (rare) prior spacecraft FD thermal analyses have been with N×N system matrices, converting lumped parameter (LP) time domain models to FD. This paper includes continuously distributed resistance and capacitance elements. Exact distributed element results differ greatly from LP renditions at higher frequency. Two FD methods are used together: a 2×2 matrix method with a bridge to an N×N matrix method. The 2×2 method allows the combination of elements in series/parallel configurations into overall elements. The 2×2 matrices can then be used in an N×N system matrix formulation. This paper includes FD analysis and programming methods, FD testing methods, and model to data comparisons. A thermal signal generator was developed for FD testing.