Higher-order Terms Research Articles

Study of Soliton Interaction in Optical Fibers with Third Order Dispersion and Higher Order Nonlinear Effects

In this paper, we present a numerical approach to solve the GNLSE and analyze soliton interaction phenomena using COMSOL environment. By leveraging the capabilities of COMSOL's PDE module, we can accurately capture the dynamics of solitons and investigate their interactions. We analyze the impact of different parameters such as soliton power, initial separation distance, and dispersion characteristics on the soliton dynamics. Furthermore, we examine the role of higher-order dispersion terms in shaping the soliton interactions. Our findings demonstrate the effectiveness of the proposed numerical approach in accurately simulating and analyzing soliton interaction phenomena. The COMSOL-based methodology provides a flexible and efficient framework for studying complex nonlinear optical systems, enabling researchers to gain insights into the behavior of solitons in different media and design optimized communication systems. This paper contributes to the understanding of soliton dynamics and provides a practical tool for investigating the behavior of solitons in nonlinear dispersive media. The presented numerical approach using COMSOL opens avenues for further research in nonlinear optics and fiber optic communication systems.

On linear stabilization of a class of nonlinear systems in a critical case

In this paper, we address the stabilization problem for nonlinear systems in a critical case. Namely, we study the class of canonical nonlinear systems. Canonical nonlinear systems or chain of power integrators is an important subject of research. Studying such systems is complicated by the fact that they cannot be mapped onto linear systems. Moreover, they have the uncontrollable first approximation. Previous results on smooth stabilization of such systems were obtained under the assumption that the powers in the right-hand side are strictly decreasing. In this work, we consider a case of non-increasing powers in the right-hand side for a three-dimensional system. A popular approach for studying such systems is the backstepping method, which is a method of step-wise stabilization. This method requires a sequential investigation of lower-dimensional subsystems. Backstepping enables the study of a wide range of nonlinear triangular systems but requires technically complex and cumbersome computations. Therefore, a natural question arises about constructing stabilizing controls of a simple form. Polynomial controls can serve as an example of such controls. In the paper, we demonstrate that linear controls can be considered as stabilizing controls. We derive sufficient conditions for the coefficients of the linear control that ensure the asymptotic stability of the zero equilibrium point of the corresponding closed-loop system. The asymptotic stability is proven using the Lyapunov function method, which is found as the sum of squares. The negative definiteness of the Lyapunov function derivative in a neighborhood of the origin guarantees asymptotic stability. In contrast to the case of strictly decreasing powers, additional conditions on the control coefficients, apart from their negativity, emerge. The obtained result extends to a broader class of nonlinear systems through stabilization by nonlinear approximation. This allows the consideration of systems with higher-order terms in the right-hand side. The effectiveness of the applied approach is illustrated by several model examples. The method used in this work to investigate the case of non-increasing powers can be applied to systems of higher dimensions.

Cyclops states in repulsive theta-neuron networks

Networks of phase oscillators have become a widely established paradigmatic model for studying emergent collective behavior across several real-world systems, including neuronal networks, populations of chemical oscillators, and power grids. The Kuramoto model, involving one-dimensional or two-dimensional phase oscillators, demonstrates the potential for networks to showcase exceptional collective dynamics. This encompasses various outcomes such as full, partial, explosive, and asymmetry-induced synchronization, clusters, chimeras, solitary states, and generalized splay states. Notably, increasing all-to-all coupling in the classical Kuramoto model induces full synchronization as the most probable outcome and dominant rhythm. Kuramoto networks with repulsive coupling usually display splay, generalized, and cluster splay states, but the conditions under which a certain rhythm can arise and prevail are not entirely understood. Equally important for connecting Kuramoto networks to practical physical systems is understanding the function of higher-order coupling terms. These terms display a Fourier decomposition of a general 2π-periodic interaction function [1]. Previous studies have demonstrated that the inclusion of higher-order terms in the classical Kuramoto model of oscillators with all-to-all attractive coupling can lead to multiple synchronous states and switching between synchronization clusters. However, the impact of higher-order coupling modes on rhythm generation in repulsive networks remains unexplored. In this work, we present significant progress in addressing the critical issue related to repulsive Kuramoto–Sakaguchi networks of phase oscillators with phase-lagged first-order and higher-order coupling. We demonstrate that weakly repulsive networks of even and odd numbers of oscillators with first-order coupling are dominated by two-cluster and three-cluster splay states, respectively. The three-cluster splay states consist of two distinct coherent clusters and one solitary oscillator. These tripod states can be considered a fusion of a two-body chimera and a solitary state. We have dubbed these patterns of three oscillators as “Cyclops states” in reference to the Greek mythological giant with a single eye. The solitary oscillator and synchronous clusters respectively represent the Cyclops’ eye and shoulders. We present a remarkable discovery that the inclusion of higher-order coupling modes leads to worldwide stability of cyclops states across almost the entire range of the phase-lag parameter controlling repulsion [2]. Beyond the Kuramoto oscillators, we demonstrate the robust presence of this effect in networks of canonical theta-neurons with adaptive coupling. Furthermore, our results provide insight into identifying dominant rhythms within repulsive physical and biological networks.

Boundary-value problems for differential-difference equations with nite and in nite orbits of boundaries

We consider boundary-value problems for differential-difference equations containing incommensurable shifts of arguments in the higher-order terms. We show that for the case when the orbits of the domain boundary generated by the set of shifts of the difference operator are nite, the original problem is similar to the boundary-value problem for differential-difference equations with integer shifts of arguments. The case of an in nite boundary orbit is also studied.

A decoupled, linearly implicit and high-order structure-preserving scheme for Euler–Poincaré equations

It is challenging to develop high-order structure-preserving finite difference schemes for the modified two-component Euler–Poincaré equations due to the nonlinear terms and high-order derivative terms. To overcome the difficulties, we introduce a bi-variate function and carefully choose the intermediate average variable in the temporal discretization. Then, we obtain a decoupled and linearly implicit scheme. It is shown that the fully-discrete scheme can keep both the discrete mass and energy conserved. And the fully-discrete scheme has fourth-order accuracy in the spatial direction and second-order accuracy in the temporal direction. Several numerical examples are given to confirm the theoretical results.

The effect of strong nonlinearity on wave-induced vertical mixing

Abstract. A semi-analytical solution to an advection–diffusion equation is coupled with a nonlinear wavemaker model to investigate the effect of strong nonlinearity on wave-induced mixing. The comparisons with weakly nonlinear model predictions reveal that in the case of waves of higher steepness, enhanced mixing affects the subsurface layer of the water column. A fully nonlinear model captures the neglected higher-order terms from a weakly nonlinear solution and provides a reliable estimation of the time-mean velocity field. The corrected wave-induced mass-transport velocity leads to improved estimates of subsurface mixing intensity and ocean surface temperature.

Temporal discretization for improving kinetic-energy and entropy preservation properties in KEEP schemes

A temporal discretization for improving kinetic-energy and entropy preservation on compressible flow simulations is proposed. The key of the proposed scheme is to introduce the auxiliary velocity, density, and internal energy to the temporal discretization and derive the quantities that satisfy the analytical relations of temporal-derivative terms in the kinetic-energy and entropy equations at the discrete level by utilizing the Maclaurin expansion and truncating high-order terms. The numerical fluxes in the kinetic-energy and entropy preserving (KEEP) schemes proposed by Kuya et al. (2018) are also evaluated by the auxiliary quantities to enhance kinetic-energy and entropy preservation in both spatial and temporal discretizations. The numerical experiments of the inviscid Taylor–Green vortex demonstrate that the proposed scheme significantly improves the kinetic energy and entropy preservation and achieves stable and non-dissipative compressible flow simulations at CFL ≈1.0, at which the temporal-discretization errors for the kinetic-energy and entropy preservation are not negligible when the widely used fourth-order four-step Runge–Kutta scheme is employed.

Delay differential equation modeling of social contagion with higher-order interactions

In this paper, we propose a social contagion model with group interactions on heterogeneous network, in which the group interactions are represented by incorporating higher-order terms. The dynamics of the proposed model are analyzed, revealing the effect of group interactions on the system dynamics. We derive a global asymptotically stability condition of the zero equilibrium point. If the parameters enhancement factor and transform probability meet the condition, the social contagion will eventually disappear. The bifurcation behavior arising from group interactions is investigated, when the enhancement factor exceeds a threshold, the system undergoes a backward bifurcation. This implies that R0<1 does not guarantee the disappearance of social contagion, we also need to control the initial values of social contagion at a lower level. The optimal control strategies for the model are provided. Moreover, the numerical simulations validate the accuracy of the theoretical analysis. It is worth noting that the group interactions lead to the emergence of a bistable phenomenon, in which the social contagion will either fade away or persist eventually, depending on the initial values.

Euler Continuity Equation with High-Order Terms in Time

Euler Continuity Equation with High-Order Terms in Time

Improved mean-field dynamical equations are able to detect the two-step relaxation in glassy dynamics at low temperatures

We study the stochastic relaxation dynamics of the Ising p-spin model on a random graph, which is a well-known model with glassy dynamics at low temperatures. We introduce and discuss a new closure scheme for the master equation governing the continuous-time relaxation of the system, which translates into a set of differential equations for the evolution of local probabilities. The solution to these dynamical mean-field equations describes the out-of-equilibrium dynamics at high temperatures very well, notwithstanding the key observation that the off-equilibrium probability measure contains higher-order interaction terms not present in the equilibrium measure. In the low-temperature regime, the solution to the dynamical mean-field equations shows the correct two-step relaxation (a typical feature of glassy dynamics), but with a too-short relaxation timescale. We propose a solution to this problem by identifying the range of energies where entropic barriers play a key role and defining a renormalized microscopic timescale for the dynamical mean-field solution. The final result perfectly matches the complex out-of-equilibrium dynamics computed through extensive Monte Carlo simulations.

Evaluation of a high-order central-difference solver for highly compressible flows out of thermochemical equilibrium

A high-order shock-capturing central finite-difference scheme is evaluated for numerical simulations of hypersonic high-enthalpy flows out of thermochemical equilibrium. The scheme is an extension to thermochemical out-of-equilibrium flows of the technique presented in Sciacovelli et al. (2021) for high-speed flows in chemical nonequilibrium. It relies on a standard tenth-order accurate central-difference approximation of the inviscid fluxes, supplemented with a high-order accurate nonlinear artificial dissipation term of ninth-order accuracy in smooth flow regions. Close to flow discontinuities, a shock-capturing low-order term is activated based on a highly selective shock sensor. To enable robust and non oscillatory solutions in regions of strong discontinuities of the thermodynamic variables, including the vibrational temperature, a shock detector consisting of a combination of a pressure-based term and Ducros’ vorticity/dilatation sensor is applied to all conservation equations except that of vibrational energy. For the latter, a shock sensor based on the vibrational temperature itself is adopted instead, to account for the loose coupling between vibrational energy and pressure. The accuracy and robustness of the proposed approach is demonstrated for a variety of thermochemical non-equilibrium configurations, ranging from one-dimensional benchmarks to three-dimensional turbulent flows, for which the ILES capabilities of the selective high-order numerical dissipation are also showcased.

Modeling of Nonlinear Sea Wave Modulation in the Presence of Ice Coverage

A model accounting for the influence of ice coverage on the propagation of surface sea waves is suggested. The model includes higher-order linear and nonlinear terms in the equation of wave motion. The asymptotic solution is obtained to account for nonlinear modulated wave propagation and attenuation. Two kinds of attenuation are revealed. The influence of the higher-order nonlinear, dispersion, and dissipative terms on the shape and velocity of the modulated nonlinear wave is studied. Despite the presence of higher-order terms in the original equation, the modulated solitary wave solution contains free parameters, which is important for the possible generation of such waves.

Linear and nonlinear optical absorption coefficients in InGaN/GaN quantum wells: Interplay between intense laser field and higher-order anharmonic potentials

This computational investigation delves into the electronic and optical attributes of InGaN/GaN nanostructures subjected to both harmonic and anharmonic confinement potentials, coupled with the influence of a nonresonant intense laser field (ILF). The theoretical framework incorporates higher-order anharmonic terms, specifically quartic and sextic terms. The solutions to the Schrödinger equation have been computed employing the finite element method and the effective mass theory. Moreover, linear and third-order nonlinear optical absorption coefficients are derived via a density matrix expansion. Our analysis reveals the feasibility of manipulating electronic and optical properties by adjusting confinement potential parameters, system attributes, and laser field intensity. In addition, the ILF induces remarkable modifications, characterized by reduced resonance peak amplitudes and a blue shift in absorption coefficients. Intriguingly, regardless of potential harmonicity, the impact of incident electromagnetic intensity is notably more pronounced in the absence of the ILF. These findings hold significant promise for advancing theoretical predictions, providing valuable insights into the intricate interplay between confinement potentials, laser fields, and their effects on electronic and optical behaviors within nanostructures.

Inferring connectivity of an oscillatory network via the phase dynamics reconstruction

We review an approach for reconstructing oscillatory networks’ undirected and directed connectivity from data. The technique relies on inferring the phase dynamics model. The central assumption is that we observe the outputs of all network nodes. We distinguish between two cases. In the first one, the observed signals represent smooth oscillations, while in the second one, the data are pulse-like and can be viewed as point processes. For the first case, we discuss estimating the true phase from a scalar signal, exploiting the protophase-to-phase transformation. With the phases at hand, pairwise and triplet synchronization indices can characterize the undirected connectivity. Next, we demonstrate how to infer the general form of the coupling functions for two or three oscillators and how to use these functions to quantify the directional links. We proceed with a different treatment of networks with more than three nodes. We discuss the difference between the structural and effective phase connectivity that emerges due to high-order terms in the coupling functions. For the second case of point-process data, we use the instants of spikes to infer the phase dynamics model in the Winfree form directly. This way, we obtain the network’s coupling matrix in the first approximation in the coupling strength.

Potential vision tester using adaptive optics, Maxwellian view, and small pupil.

We demonstrate a free-space, trolley-mounted potential vision tester (PVT), designed to study and improve the accuracy of visual acuity (VA) measurements in the aging eye. Key features include a high-resolution visual display presented in Maxwellian view, a 3 mm pupil to limit wavefront (WF) aberrations, and a moderate cost deformable mirror to induce or correct higher order optical aberrations. The visual display supported accurate measurement of visual acuities down to 20/5. The moderate cost, piezo deformable mirror induced seven nominal aberrations, calibrated as 0, -0.32, -0.23, + 0.27, and +0.39 microns spherical aberration; + 0.49 microns Y coma; and -0.51 microns X coma. A custom Hartmann Shack (HS) calibration (HSc) system demonstrated that induced aberrations were repeatable and stable. A Badal optometer provided the coarse focus. WF aberrations were measured for five normal subjects with a commercially available HS device (HSP) (OCULUS Pentacam AXL Wave), providing estimates of WF errors for 3 mm and other pupil sizes. VA was measured using four alternative forced-choice for a single black on white E stimulus in each trial. Using the method of constant stimuli yielded robust standard deviation measurements. The 50% fit for VA plotted against induced aberration resulted in linear functions for each subject for the range of our positive and negative spherical aberration data. Subjects differed, but higher order terms were unnecessary to describe data across spherical aberrations.

A dimensionless analytical analysis for buckling and lateral buckling interaction of thin-walled beams with open cross sections

This paper investigates the interaction between buckling and lateral buckling of thin-walled beams with arbitrary open cross-section using nonlinear modeling in large rotation and with warping. First, the equilibrium equations have been transformed into dimensionless ones, and within a nonlinear stability model, various dimensionless parameters are introduced to control the interaction between buckling and lateral buckling. In comparison to a reference of literature, this work focuses on the case of arbitrary sections and includes a greater number of higher-order terms in our analysis. These additions have an impact on the coefficients of the stability formula by expanding its applicability to non- symmetric sections. Then, the effect of the different dimensionless parameters, thus introduced on the stability curve, has been studied. The proposed solutions are validated and compared to numerical solutions obtained by a general finite element package including the warping developed by the authors using Asymptotic Numerical Method (ANM) to solve the initial nonlinear equations.

Parametric Updating of Steel Strand Stay Cable Force Calculation Model Based on Response Surface Method

In order to solve the problem that the calculation result of cable force deviates from the measured value because the parameters of steel strand stay cable are not exactly equal to the designed value in practical engineering, a model updating method selecting the steel strand stay cable parameters as the updating parameters and using more measured data as the objective function was proposed. This method selects the central composite design in the response surface method. In the fine-fitting model affected by the cross-term and high-order terms between parameters, the value of the parameter to be corrected is optimized according to the measured data. Using the above method, taking a cable-stayed bridge without a backstay as the engineering background, the model updating is carried out by using the cable force monitoring data during the construction monitoring process. The results show that using this method, the updated values of the parameters of the cable can be obtained according to the measured data, and the updated values of the parameters are consistent with the reality. The deviation between the calculated value of cable force and the actual cable tension obtained according to the parameter updating decreases significantly, and most of the deviation is reduced from 30% to 50% before the updating to less than 5%. The updating effect is obvious and meets the engineering requirements.

A New Regression Model for the Prediction of the Stress-Strain Relations of Different Materials.

Experimental flow stress-strain data under different stress states are often used to calibrate the plastic constitutive model of anisotropic metal materials or identify the appropriate model that is able to reproduce their plastic deformation behavior. Since the experimental stress-strain data are discrete, they need to be mathematically returned to a continuous function to be used to describe an equivalent hardening increment. However, the regression results obtained using existing regression models are not always accurate, especially for stress-strain curves under biaxial stress loading conditions. Therefore, a new regression model is proposed in this paper. The highest-order term in the recommended form of the new model is quadratic, so the functional relationships between stress-strain components can be organized into explicit expressions. All the experimental data of the uniform deformation stage can be substituted into the new model to reasonably reproduce the biaxial experimental stress-strain data. The regression results of experimental data show that the regression accuracy of the new model is greatly improved, and the residual square sum SSE of the regression curves of the new model reduced to less than 50% of the existing three models. The regression results of stress-strain curves show significant differences in describing the yield and plastic flow characteristics of anisotropic metal materials, indicating that accurate regression results are crucial for accurately describing the anisotropic yielding and plastic flow behaviors of anisotropic metal materials.

Analysis of probability of inserting a hard spherical particle with small diameter in hard-sphere fluid

The probability of inserting, without overlap, a hard spherical particle of diameter σ in a hard-sphere fluid of diameter σ0 and packing fraction η determines its excess chemical potential at infinite dilution, μex(σ, η). In our previous work [R. L. Davidchack and B. B. Laird, J. Chem. Phys. 157, 074701 (2022)], we used Widom’s particle insertion method within molecular dynamics simulations to obtain high precision results for μex(σ, η) with σ/σ0 ≤ 4 and η ≤ 0.5. In the current work, we investigate the behavior of this quantity at small σ. In particular, using the inclusion-exclusion principle, we relate the insertion probability to the hard-sphere fluid distribution functions and thus derive the higher-order terms in the Taylor expansion of μex(σ, η) at σ = 0. We also use direct evaluation of the excluded volume for pairs and triplets of hard spheres to obtain simulation results for μex(σ, η) at σ/σ0 ≤ 0.2247 that are of much higher precision than those obtained earlier with Widom’s method. These results allow us to improve the quality of the small-σ correction in the empirical expression for μex(σ, η) presented in our previous work.

New simple local smoothness indicators for fifth-order WENO schemes simulating compressible flows

For the WENO-Z scheme, the structure of the local smoothness indicators has a significant impact on the computational performance of the weighted essentially non-oscillatory scheme with a given global smoothness indicator. In this paper, a new class of two-parameter local smoothness indicators is constructed by combining the classical fifth-order WENO-JS and WENO-UD5 schemes and locally appending the coefficients of the higher-order terms of the new local smoothness indicators. A new WENO scheme, WENO-NSLI, is constructed using the global smoothness indicators of WENO-UD5. Numerical experiments show that the new scheme can achieve optimal accuracy at the critical point when the parameter (λ,μ)=(0.5,10) is taken and has a higher resolution than WENO-JS, WENO-Z, and WENO-UD5.