Nonlinear Solutions Research Articles

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.

Multi-fidelity wavelet neural operator surrogate model for time-independent and time-dependent reliability analysis

Operator learning frameworks have recently emerged as an effective scientific machine learning tool for learning complex nonlinear operators of differential equations. Since neural operators learn an infinite-dimensional functional mapping, it is useful in applications requiring rapid prediction of solutions for a wide range of input functions. A task of a similar nature arises in many applications of uncertainty quantification, including reliability estimation and design under uncertainty, each of which demands thousands of samples subjected to a wide range of possible input conditions, an aspect to which neural operators are specialized. Although the neural operators are capable of learning complex nonlinear solution operators, they require an extensive amount of data for successful training. Unlike the applications in computer vision, the computational complexity of the numerical simulations and the cost of physical experiments contributing to the synthetic and real training data compromise the performance of the trained neural operator model, thereby directly impacting the accuracy of uncertainty quantification results. We aim to alleviate the data bottleneck by using multi-fidelity learning in neural operators, where a neural operator is trained by using a large amount of inexpensive low-fidelity data along with a small amount of expensive high-fidelity data. We propose the multi-fidelity wavelet neural operator, capable of learning solution operators from a multi-fidelity dataset, for efficient and effective data-driven reliability analysis of dynamical systems. We illustrate the performance of the proposed framework on bi-fidelity data simulated on coarse and refined grids for spatial and spatiotemporal systems.

Vibration Correlation Technique: Methodology Applied to Buckling of Large-Scale Sandwich Cylindrical Shell

This paper assesses the Vibration Correlation Technique's (VCT) applicability within the context of sandwich cylindrical shells through detailed numerical investigation. The study utilizes a sandwich shell from NASA's Shell Buckling Knockdown Factor project. Finite element models are implemented, incorporating initial mid-surface and thickness imperfections to replicate NASA's buckling test campaign. The numerical results are compared to the experimental data, validating the nonlinear solution with an approximate deviation of 1%. Subsequently, the VCT experimental campaign is numerically conducted through free vibration analyses at different load levels, enabling a comprehensive evaluation of VCT's applicability. The results verified the effectiveness of the VCT, demonstrating a deviation of less than 10% when estimating the buckling load for load levels below 65% of the nonlinear buckling load. Overall, the findings confirm the non-destructive nature of the VCT when employed on such structures, supporting future practical applications.

A Vertical Surface Flow Analysis of MHD Nanofluids Influenced by Radiation, Viscous Dissipation, and Joule Heating with Soret and Dufour Effects

The purpose of this study is to investigate the effects of radiation and activation energy on mass transfer nanofluids flowing across an expanding sheet in a vertical direction, as well as joule heating, viscous dissipation, Soret and Dufour, and magnetohydrodynamic (MHD) flow. The boundary layer (BL) equations for nonlinear momentum, energy, solute, and nanoparticle volume fraction are reduced by using similarity transformations. After that, the shooting technique and bvp4c are used to numerically solve the boundary layer equations. Schmidt and Eckert’s numbers, along with Soret effects on velocity, temperature, and the volume percentage of nanoparticles, cause an increase in skin friction. Fluid temperature rises significantly with Eckert number. As the Dufour and Schmidt numbers rise, so does the local Nusselt number. A greater thermophoresis parameter, enhanced skin friction, and an increased Schmidt number.

CFD stability improvement using dynamic mode decomposition of solution update vectors

We present a novel method for improving the numerical stability of finite volume simulations by optimizing the mesh through dynamic mode decomposition of solution update vectors. Our approach leverages dynamic mode decomposition to approximate the non-linear solution evolution as a linear mapping, enabling us to represent the dynamic system with fewer degrees of freedom. By conducting eigenanalysis on the reduced system, we gain insights into the growth rate and oscillation frequency of dominant solution modes. We then identify the control volumes and vertices that have a significant influence on each dynamic solution mode. We compute the gradients of the Jacobian diagonal elements with respect to the movement of the selected vertices. Based on these gradients, we adjust the positions of the vertices to improve the diagonal dominance of the corresponding Jacobian rows. We utilize this approach in conjunction with our in-house flow solver as well as with Ansys Fluent to test its effectiveness when applied to different types of CFD software architecture. We show that the presented methodology is fully non-invasive to the host flow solver and does not require any modifications or access to the source code. Our approach offers a substantial computational cost reduction compared to existing methods for numerical stability improvement. Various results demonstrate the strength and efficacy of this state-of-the-art approach for improving numerical stability through mesh optimization.

Thermomechanical responses of steel frames under aftershock after fires

A novel advanced fiber beam-column element, which incorporates stability functions and the fiber distributed plasticity model to capture geometric and material nonlinearities, respectively, has been successfully developed to predict the nonlinear inelastic thermo-mechanical behaviors of steel frames subjected to aftershock after fires. To address the nonlinear dynamic problems resulting from aftershocks after fires, a sequential nonlinear thermal incremental-iterative solution scheme has been developed, which is based on the Newton-Raphson algorithm and Newmark-β method. In addition, a transcendental probability-based method is employed to express the statistical characteristics of aftershock intensity relative to the mainshock. The proposed method is applied to analyze three framed structures: a steel portal frame, a two-story plane steel frame, and a four-story asymmetric space steel frame. This analysis shows the accuracy, computational performance, and applicability of the proposed method. Remarkably, the results obtained from the code work, using just one or two elements per member, are in close agreement with all the results analyzed by Abaqus. This emphasizes that the proposed method effectively predicts the response of framed structures under aftershock after fires with excellent computational efficiency, significantly surpassing traditional finite element approaches in commercial software like Abaqus. The successful results of the dynamic analysis of the frames under aftershocks following fires underscore the importance of considering aftershocks after a fire incident. It is demonstrated that the occurrence of an aftershock accelerates the failure of the frame compared to the case when the frame is subjected to fire alone. This highlights the critical need to account for aftershocks when assessing the safety and structural integrity of framed structures after a fire event.

Lump, Breather, Ma-Breather, Kuznetsov–Ma-Breather, Periodic Cross-Kink and Multi-Waves Soliton Solutions for Benney–Luke Equation

The goal of this research is to utilize some ansatz forms of solutions to obtain novel forms of soliton solutions for the Benney–Luke equation. It is a mathematically valid approximation that describes the propagation of two-way water waves in the presence of surface tension. By using ansatz forms of solutions, with an appropriate set of parameters, the lump soliton, periodic cross-kink waves, multi-waves, breather waves, Ma-breather, Kuznetsov–Ma-breather, periodic waves and rogue waves solutions can be obtained. Breather waves are confined, periodic, nonlinear wave solutions that preserve their amplitude and shape despite alternating between compression and expansion. For some integrable nonlinear partial differential equations, a lump soliton is a confined, stable solitary wave solution. Rogue waves are unusually powerful and sharp ocean surface waves that deviate significantly from the surrounding wave pattern. They pose a threat to maritime safety. They typically show up in solitary, seemingly random circumstances. Periodic cross-kink waves are a particular type of wave pattern that has frequent bends or oscillations that cross at right angles. These waves provide insights into complicated wave dynamics and arise spontaneously in a variety of settings. In order to predict the wave dynamics, certain 2D, 3D and contour profiles are also analyzed. Since these recently discovered solutions contain certain arbitrary constants, they can be used to describe the variation in the qualitative characteristics of wave phenomena.