Recognition, decoding and tracking of vortex patterns is of increasing importance in many fields reaching from astronomical observation of distant galaxies to turbulence phenomena in liquids or gases. Currently, coherent light beams with orbital angular momentum (OAM) are of particular interest for optical communication, metrology, micro-machining or particle manipulation. One common task is to identify characteristic spiral patterns in pixelated intensity maps at real-world signal-to-noise ratios. Recently, Fast Fourier Transform (FFT) was applied to the quantitative analysis of spiral interference patterns of OAM beams. Our study considers this approach in more detail. It is demonstrated that specific information can be extracted by combining conformal mapping by polar transform with determining one- or two-dimensional angular frequency spectra either in spatially fixed or scanning mode. Results of proof-of-principle experiments are presented.

The application of non-uniform spiral gratings to control structure, topological parameters and propagation of orbital angular momentum (OAM) beams was studied experimentally with coherent near-infrared light. Adapted digital spiral grating structures were programmed into the phase map of a high-resolution liquid-crystal-on-silicon spatial light modulator (LCoS-SLM). It is shown that characteristic spatio-spectral anomalies related to Gouy phase shift can be used as pointers to quantify rotational beam properties. Depending on sign and gradient of spatially variable periods of chirped spiral gratings (CSG), variations of rotation angle and angular velocity as a function of propagation distance were measured. Propagation-dependent self-torque is introduced in analogy to known local self-torque phenomena of OAM beams as obtained by the superposition of temporally chirped or phase modulated wavepackets.

In this study, we focus on an abstract form of summation equations, which can be regarded as a particular class of nonlinear difference equations. We investigate sufficient conditions for the existence of discrete almost automorphic solutions by appealing to the fixed point theory and establish a linkage between the existence of discrete almost automorphic solutions and the existence of bounded solutions for the nonlinear difference equations in the spotlight. Consequentially, we propose a discrete counterpart of the Bohr-Neugebauer theorem for the qualitative theory of difference equations.

In the paper, by virtue of techniques in combinatorial analysis, the authors simplify two families of nonlinear ordinary differential equations in terms of the Stirling numbers of the first kind.

The study of different forms of preconditioners for solving a system of nonlinear equations, by using Newton’s method, is presented. The preconditioners provide numerical stability and rapid convergence with reasonable computation cost, whenever chosen accurately. Different families of iterative methods can be constructed by using a different kind of preconditioners. The multi-step iterative method consists of a base method and multi-step part. The convergence order of base method is quadratic and each multi-step add an additive factor of one in the previously achieved convergence order. Hence the convergence of order of an m-step iterative method is m + 1. Numerical simulations confirm the claimed convergence order by calculating the computational order of convergence. Finally, the numerical results clearly show the benefit of preconditioning for solving system of nonlinear equations.

We present three reduced integrable hierarchies of nonlocal integrable NLS type equations from a vector integrable hierarchy associated with a matrix Lie algebra, not being A type. Three similarity transformations are taken to keep the invariance of the transformed zero curvature equations. The key step is to formulate a solution to a reduced stationary zero curvature equation so that the zero curvature formulation works for a reduced case.

We propose a new hierarchy of the vector derivative nonlinear Schr{\"o}dinger equations and consider the simplest multiphase solutions of this hierarchy. The study of the simplest solutions of these equations led to the following results. First, the three-leaf spectral curves $\Gamma=\{(\mu,\lambda)\}$ of the simplest multiphase solutions have a quite simple symmetry. They are invariant with respect to holomorphic involution $\tau$. The type of this involution depends on the genus of the spectral curve. Or the involution has the form $\tau:(\mu,\lambda)\to(\mu,-\lambda)$, or $\tau:(\mu,\lambda)\to(-\mu,-\lambda)$. The presence of symmetry leads to the fact that the dynamics of the solution is determined not by the entire spectral curve $\Gamma$, but by its factor $\Gamma/\tau$, which has a smaller genus. Secondly, it turned out that the dynamics of the two-component vector $\bp=(p_1,p_2)^t$ is determined, first of all, by the dynamics of its length $\abs{\bp}$. Independent equations determine the dependence of the direction of the vector $\bp$ from its length. In cases where the direction of the vector $\bp$ is fixed, the corresponding spectral curve splits into separate components. In conclusion, we note that, as in the case of the Manakov system, the equation of the spectral curve is invariant with respect to the orthogonal transformation of the vector solutions. I.e., the solution can be found from the spectral curve up to the orthogonal transformation. This fact indicates that the spectral curve does not depend on the individual components of the solution, but on their symmetric functions. Thus, the spectral data of multiphase solutions have two symmetries. These symmetries make it difficult to reconstruct signals from their spectral data. The work contains examples illustrating these statements.

In the last years, a recent growing line of research has been developed with fruitful results by using fractional and conformable derivatives in the iterative procedures of classical methods for solving nonlinear equations. In that sense, the use of conformable derivatives has shown better behavior than fractional ones, not only in the theory, but also in the practice. In this work, we adapt the approximation of conformable derivatives in order to design the first conformable derivative-free iterative schemes to solve nonlinear equations: a Steffensen’s type method and a Secant type method; the latter with memory. Convergence analysis is made, preserving the order of classical cases, and the numerical performance is studied in order to confirm the results in the theory. It is shown that these methods can present some numerical advantages versus their classical partners, with wide sets of converging initial estimations.

The necessary and sufficient conditions of existence of the nonlinear operator equations' branches of solutions in the neighbourhood of branching points are derived. The approach is based on reduction of the nonlinear operator equations to finite-dimensional problems. Methods of nonlinear functional analysis, integral equations, spectral theory based on index of Kronecker-Poincare, Morse-Conley index, power geometry and other methods are employed. Proposed methodology enables justification of the theorems on existence of bifurcation points and bifurcation sets in the nonstandard models. Formulated theorems are constructive. For a certain smoothness of the nonlinear operator, the asymptotic behaviour of the solutions is analysed in the neighbourhood of the branch points and uniformly converging iterative schemes with a choice of the uniformization parameter enables the comprehensive analysis of the problems details. General theorems are illustrated on the nonlinear integral equations.

In the present work we study the existence of solutions for hybrid nonlinear fractional integrodifferential equations. We developed an algorithm by using the operator theoretical techniques in order to obtain the approximate solutions. The main results depend on the Dhage iteration method that were incorporated with the modern hybrid fixed point theorems. The approximate solutions were obtained by using Lipschitz conditions and weaker form of mixed partial continuity. Further, we provide some examples to explain the hypotheses and the related results.

Traditional population-based metaheuristic algorithms are effective in solving complex real-world problems but require careful strategy selection and parameter tuning. Parameter-less population-based optimization algorithms have gained importance due to their simplicity and efficiency. However, research on their applicability for solving large systems of nonlinear equations is still incipient. In this paper, a systematic review and detailed description of metaphor-less optimization algorithms is made, including the Jaya algorithm, the main Jaya variants for continuous optimization problems, with a focus on the Enhanced Jaya (EJAYA) algorithm, the three Rao algorithms, the best-worst-play (BWP) algorithm, and the very recent max-min-greedy-interaction (MaGI) algorithm. This paper also discusses a few previous parallelizations of the Jaya optimization algorithm, as well as recent GPU-based massively parallel implementations of the Jaya, Enhanced Jaya (EJAYA), Rao, and BWP algorithms developed by the authors. In addition, a novel GPU-accelerated version of the MaGI algorithm is proposed. The GPU-accelerated versions of the metaphor-less algorithms developed by the authors were implemented using the Julia programming language and tested primarily on GeForce RTX 3090 and NVIDIA A100 GPUs, as well as on Tesla V100S and Tesla T4 GPUs, using a set of difficult, large-scale nonlinear equation system problems. The computational experiments carried out produced quite significant speedups, which highlights the efficiency of the GPU-based versions of the metaphor-less algorithms developed for solving large-scale systems of nonlinear equations.

In this paper, we investigate the equilibrium points, stability of two equilibrium points, convergences of negative equilibrium point, periodic solutions, and existence of bounded or unbounded solutions of a system of nonlinear difference equations x_n+1 =x_n-1y_n - 1, y_n+1 = y_n-1x_n - 1 n = 0,1,..., where the initial values are real numbers. Additionally we present some numerical examples to verify our theoretical results.

This paper deals with the form, the stability character, the periodicity and the global behavior of solutions of the following four rational difference equations x_{n+1} = ((±1)/(x_{n}(x_{n-1}±1)-1)) x_{n+1} = ((±1)/(x_{n}(x_{n-1}∓1)+1)).

We investigate the exact travelling wave solutions of the Tzitzéica-Dodd–Bullough–Mikhailov (TDBM) and Tzitzéica-type (TT) equations, appear in nonlinear optics, using simple ansatz approach. We use the Painlevé and traveling wave transformations to convert the mentioned nonlinear equations into nonlinear ordinary differential equations. Then, the simple ansatz approach has been employed to procure wave solutions, in terms of bright, singular soliton and trigonometric functions, for the Tzitzéica-Dodd–Bullough–Mikhailov (TDBM) and Tzitzéica-type (TT) equations. The study shows that the simple ansatz approach is straightforward and easy to implement for solving nonlinear problems.

This study presents the solution of the second type of a two-dimensional nonlinear integral equation in Banach space. Also, the existence and uniqueness of this equation’s solution are discussed. We utilize a numerical approach involving hybrid and block-pulse functions to obtain the approximate solution of a two-dimensional nonlinear integral equation. Nonlinear integral equation in two dimensions is reduced numerically to a system of nonlinear algebraic equations that can be solved using numerical methods. This study focuses on showing the convergence analysis for the numerical approach and obtaining an error estimate. Some numerical examples have been provided to demonstrate the approach’s viability and efficacy

This paper discusses approaches to the numerical integration of the coupled nonlinear Schrödinger equations system in case of few-mode wave propagation. The wave propagation assumes the propagation of up to nine modes of light in an optical fiber. In this case, the light propagation is described by the non-linear coupled Schrödinger equation system, where propagation of each mode is described by own Schrödinger equation with other modes interactions. In this case, the non-linear coupled Schrödinger equation system solving becomes increasingly complex, because each mode affects the propagation of other modes. The suggested solution is based on the direct numerical integration approach, which is based on a finite-difference integration scheme. The well-known explicit finite-difference integration scheme approach fails, due to the non-stability of the computing scheme. Due to this fact, the combined explicit/implicit finite-difference integration scheme, based on the implicit Crank–Nicolson finite-difference scheme, is used. It allows ensuring the stability of the computing scheme. Moreover, this approach allows separating the whole equation system on the independent equation system for each wave mode at each integration step. Additionally, the algorithm of numerical solution refining at each step and the integration method with automatic integration step selection are used. The suggested approach has performance gains (or resolutions) up to three or more orders of magnitude in comparison with the split-step Fourier method due to the fact that there is no need to produce direct and inverse Fourier transforms at each integration step. The main advantage of the proposed method is the ability to calculate the propagation of an arbitrary number of modes in the fiber.

In this paper, a class of nonlinear ordinary differential equations with impulses at variable times is considered. The existence and uniqueness of solution are given. At the same time, modifying the classical definitions of continuous dependence and Ga^teaux differentiability, some results on continuous dependence and Ga^teaux differentiable of solution relative to the initial value also are presented in new topology sense. For the autonomous impulsive system, the periodicity of solution is given. As an application, properties of solution for a type of controlled nonlinear ordinary differential equation with impulses at variable times is obtained. These results are foundation to study optimal control problems of systems governed by the differential equations with impulses at variable times.

In this work, the analytical treatment of Fourth order pd equations has been discussed with the variational iteration algorithm with an auxiliary parameter. The technique has the capability of reducing the size of computational work and easily overcomes the difficulty of the perturbation method or Adomian polynomials. Comparison with the classic variational iteration algorithm- is carried out, showing that the modification is more efficient and reliable.

At present paper, we aim to study the long-time behavior of a stochastic semi-linear degenerate parabolic equation on bounded or unbounded domain and driven by a nonlinear noise and defined. Since the theory of pathwise random dynamical systems can not be applied directly to the equation with nonlinear noise, first, we establish the existence of weak pullback mean random attractors for the equation by applying the theory of mean-square random dynamical systems; then, we prove the existence of (pathwise) pullback random attractors for the Wong-Zakai approximate system of the equation. In addition, we establish the upper semicontinuity of pullback random attractors for the Wong-Zakai approximate system of the equation under consideration driven by a linear multiplicative noise.

The aim of this paper is to deep into the dynamic study of the well known Chebyshev-Halley family of iterative methods for solving nonlinear equations. Our objetives are twofold: one the one hand, we are interested in characterize the existence of extraneous attracting fixed points when the methods in the family are applied to polynomial equations. On the other hand, we are also interested in study the free critical points of the methods in the family, as a previous step to determine the existence of attracting cycles. In both cases, we want to identify situations where the methods in the family have a bad behavior from the root-finding point of view. Finally, and joining these two studies, we look for polynomials for which there are methods in the family where this two situations happen simultaneously. That is, the rational map obtained by applying a method in the Chebyshev-Halley family to a polynomial has both super-attracting extraneous fixed point and super-attracting cycles different to the roots of the polynomial.

The splitting method developed in a previous study for linear, multi – component, convection-diffusion mass transfer equations is extended to non-linear, multi-component, convection-diffusion mass transfer equations. The diffusion term is non-linear (the diffusion coefficients depend on chemical species concentrations). Theoretical results about the splitting error are presented, as well as numerical simulations for mixtures with n ≤ 5. The numerical experiments show second-order accuracy and stability for the present splitting algorithm.

In this work, a mixture of Elzaki transform and projected differential transform method is applied to solve a nonlinear wave-like equations with variable coefficients. Nonlinear terms can be easily manipulated by using the projected differential transformation method. The method gives the results show that the proposed method is very efficient, simple and can be applied to other applications.

This paper uses the $\exp(-\Phi(\xi))$-Expansion method to investigate solitons to the M-fractional nonlinear Schrödingers equation with cubic nonlinearity. The results obtained are dark solitons, trigonometric function solutions, hyperbolic solutions and rational solutions. Thus, the constraint relations between the model coefficients and the traveling wave frequency coefficient for the existence of solitons solutions are also derived.

In wireless sensor networks (WSNs), each sensor node can estimate the global parameter from the local data in distributed manner. This paper proposed a robust diffusion estimation algorithm based on minimum error entropy criterion with self-adjusting step-size, which are referred to as diffusion MEE-SAS (DMEE-SAS) algorithm. The DMEE-SAS algorithm has fast speed of convergence and is robust against non-Gaussian noise in the measurements. The detailed performance analysis of the DMEE-SAS algorithm is performed. By combining the DMEE-SAS with diffusion minimum error entropy (DMEE) algorithms, an Improving DMEE-SAS algorithm is proposed, in non-stationary environment where tracking is very important. The Improving DMEE-SAS algorithm can avoid insensitivity of the DMEE-SAS algorithm due to the small effective step-size near the optimal estimator, and obtain a fast convergence speed. Numerical simulations are given to verify the effectiveness and advantages of these proposed algorithms.

The chiral optical Tamm state (COTS) is a special localized state at the interface of a handedness-preserving mirror and a structurally chiral medium such as a cholesteric liquid crystal or a chiral sculptured thin film. The spectral behavior of COTS, observed as reflection resonances, is described by the temporal coupled-mode theory. Mode coupling is different for two circular light polarizations because COTS has a helix structure replicating that of the cholesteric. The mode coupling for co-handed circularly polarized light exponentially attenuates with the cholesteric layer thickness since the COTS frequency falls into the stop band. Cross-handed circularly polarized light freely goes through the cholesteric layer and can excite COTS when reflected from the handedness-preserving mirror. The coupling in this case is proportional to anisotropy of the cholesteric and theoretically it is only anisotropy of magnetic permittivity that can ultimately cancel this coupling. These two couplings being equal results in a polarization crossover (the Kopp--Genack effect) for which a linear polarization is optimal to excite COTS. The corresponding cholesteric thickness and scattering matrix for COTS are generally described by simple expressions.

In this paper, we present a new parametric class of optimal fourth-order iterative methods to estimate the solutions of nonlinear equations. After the convergence analysis, we study the stability of this class by using the tools of complex discrete dynamics. This allows us to select those elements of the class with lower dependence on initial estimations. We verify by means of numerical tests that the stable members on quadratic polynomials perform better than the unstable ones, when applied to other non-polynomial functions. We also compare the performance of the best elements of the family with known methods, such as the schemes by Newton or Chun, among others, showing robust and stable behaviour.

The goal of the current study is to analyse several nonlinear two-dimensional time-fractional Rosenau Hymanequations. The two-dimensional fractional Rosenau-Hyman equation has extensive use in engineering and applied sciences. The fractional view analysis of two-dimensional time-fractional Rosenau-Hyman equations is discussed using the homotopy perturbation approach, adomian decomposition method, and Yang transformation. Some examples involving two-dimensional time-fractional Rosenau-Hyman equations are provided in order to better understand the suggested approaches. The solutions appear as infinite series. We offer a comparison between the accurate solutions and those that are generated employing the proposed approaches in order to demonstrate the effectiveness and applicability of the proposed techniques. The results are graphically illustrated using 2D and 3D graphs. It has been noted that the obtained results and the targeted problems real solutions are quite similar. Calculated solutions at various fractional levels describe some of the problems useful dynamics. Other fractional problems that arise in other fields of science and engineering can be solved using a modified version of the current techniques.

The numerical implementation of the controllers allows the use of very complex algorithms with ease. However, in practice, due to its proven advantages, the PID controller (and its variants) is widely used in industrial control systems as well as in many other applications that require continuous control. Most of the methods for tuning the parameters of PID controllers are based on time-invariant linear models of the processes, which in practice can lead to poor performance of the control system. The paper presents an application of reinforcement learning algorithms in tuning of PID controllers for the control of some classes of continuous nonlinear systems. Tuning the parameters of PID controllers is done with the help of Twin Delayed Deep Deterministic Policy Gradients (TD3) algorithm which presents a series of advantages compared to other similar methods from Machine Learning dedicated to continuous state and action spaces. TD3 algorithm is an off-policy Actor-Critic based method and was used as it does not require a system model. The presented technique is applied for control of a biotechnological system which has a strongly nonlinear dynamic. The proposed tuning method is compared to the classical tuning methods of PID controllers. The performance of the tuning method based on the TD3 algorithm is demonstrated through simulation illustrating the effectiveness of the proposed methodology.

This paper applied the space-time generalized finite difference scheme for solving the nonlinear dispersive shallow water waves as the modified Camassa−Holm equation, modified Degasperis-Procesi equation, Fornberg-Whitham equation, and its modified form. The proposed meshless numerical scheme was composed of the space-time generalized finite difference method, the two-step Newton-Raphson method, and the time-marching method. The space-time approach can treat the temporal derivative as one of the spatial derivatives. This numerical technique enables all the partial derivatives in the governing equation can be discretized by a spatial discretization method, and the mixed derivative can efficiently deal with using the proposed meshless numerical scheme. The space-time generalized finite difference method is advanced from the Taylor series expansion and the moving-least square method. The numerical discretization process is only related to functional data and weighting coefficients on the central node and its nearby nodes. Thus, the matrix system composed of nonlinear algebraic equations will be a sparse matrix and can be efficiently solved by the two-step Newton-Raphson method. Furthermore, the time-marching method was utilized to proceed with the space-time domain along the time axis. In this paper, several numerical examples were tested to verify the capability of the proposed space-time generalized finite difference scheme.

While numerous studies have explored the correlation between the built environment and ride-hailing demand, few have assessed their nonlinear interplay. Utilizing ride-hailing order data and multi-source built environment data from Nanjing, China, this paper uses the machine learning method, eXtreme Gradient Boosting (XGBoost), combined with SHapley Additive exPlanations (SHAP) and Partial Dependence Plots (PDP) to investigate the impact of built environment factors on ride-hailing travel demand, including their non-linear and threshold effects. The findings reveal that dining facilities have the most significant impact, with a contribution rate of 30.75%, on predicting ride-hailing travel demand. Additionally, financial, corporate, and medical facilities also exert considerable influence. The built environment factors need to reach a certain threshold or within a certain range to maximize the impact of ride-hailing travel demand. population density, land use mix and distance to the subway station collectively influence ride-hailing demand. The results are helpful for TNCs to allocate network ride-hailing resources reasonably and effectively.

This paper proposed an optimal tracking control scheme through adaptive dynamic programming(ADP) for a class of partially unknown discrete-time nonlinear systems based on radial basis function neural network(RBF-NN). In order to acquire the unknown system dynamics, we use two RBF-NNs, the one is used to construct the identifier, and the another is used to directly approximate the steady-state control input, where a novel adaptive law is proposed to update neural network weights. While the optimal feedback control and the cost function are derived via feedforward neural networks approximating, it is proposed to regulate the tracking error, the critic network and the actor network are then trained online to obtain the solution of the associated Hamilton–Jacobi–Bellman (HJB) equation being built under the ADP framework. Simulations verify the effectiveness of the optimal tracking control technique using the neural networks.

Partial Differential Equations (PDEs) constitute a fundamental framework for modeling various physical phenomena across diverse fields such as physics, engineering, and finance. Solving PDEs accurately and efficiently remains a significant challenge, particularly when dealing with complex systems or high-dimensional data. In this paper, we present a novel approach based on deep learning techniques to investigate the effective factors in solving PDEs. Our methodology leverages the expressive power of deep neural networks to learn and represent the underlying relationships between input parameters and the solutions of PDEs. By incorporating domain knowledge into the network architecture and training process, our approach aims to provide insights into the key factors influencing the solution of PDEs. We demonstrate the effectiveness of our method through numerical experiments on various PDEs, showcasing its capability to identify important features and improve predictive accuracy. Our findings suggest that deep learning techniques offer a promising avenue for understanding and analyzing the complex dynamics inherent in PDEs, paving the way for enhanced computational methods in scientific and engineering applications.

The fractional advection-dispersion equation FADE plays an important role, especially in explaining the process of solute-transport in natural porous media with different degrees of heterogeneity. However, studies of the flow and transport of contaminants in the underground environment are based on the knowledge of many functional parameters of the subsoil. The delay factor caused by nonlinear absorption is one of these parameters inheriting significant ambiguities in its estimation in heterogeneous media. In the study, a fractional advection-dispersion FADE model was developed by the implicit Euler method, based on a Grunwald and Caputo-Liousville approximation to explicitly evaluate and simulate the transport mechanism of contaminants in media with unsaturated porous materials subjected to a nonlinear sorption model. Better solute mass retention was observed for fractional advection-dispersion FADE model. Moreover, the absorption parameter of the retardation coefficient has a better impact on the motion of the solute in a saturated porous medium. In other words, it has a significant effect on the dispersion of pollutants in aquifers. This study highlights the need to characterize the aquifer to predict solute transport in porous media.

This paper considers the Hermitian solutions of a new system of commutative quaternion matrix equations, we establish both necessary and sufficient conditions for the existence of solutions. Furthermore, we derive an explicit general expression when it is solvable. To illustrate our main findings, we present two numerical algorithms and examples in this paper.

Local polynomial smoother (LPS) is a weighted local least-squares nonparametric method. It provides a local Taylor series fit of the data at any location and can be directly used in a differential equation to provide a numerical scheme. In this article, we introduce this new nonparametric idea based on local polynomial smoother, for acquiring the numerical solution of the Bagley-Torvik fractional-order differential equations. Furthermore, this paper will present a numerical comparison with some methods, such as legendre operational matrix and pseudo-spectral method. The efficiency and accuracy of the LPS method are demonstrated by two numerical examples.

The stability properties of the Hill equation are discussed, and especially those of the Mathieu equation that characterize ion motion in electrodynamic traps. The solutions of the Mathieu equation for a trapped ion are characterized by using the Floquet theory and Hill's Method solution, which yields an infinite system of linear and homogeneous equations whose coefficients are recursively determined. Stability is discussed for parameters $a$ and $q$ that are real. Characteristic curves are introduced naturally by the Sturm-Liouville problem for the well known even and odd Mathieu equations $ce_m(z,q)$ and $se_m(z,q)$. We illustrate the stability diagram for a combined (Paul and Penning) trap and represent the frontiers of the stability domains for axial and radial motion. In case of a Paul trap the stable solution corresponds to a superposition of harmonic motions. The problem of evaluating the maximum amplitudes of stable oscillations for the ideal conditions (taken into consideration) is also approached. Anharmonic corrections are discussed within the frame of the perturbation theory, while the frontiers of the modified stability domains are determined as a function of the chosen perturbation parameter. The results apply to 2D and 3D ion traps used for different applications in quantum engineering, among which optical clocks, quantum logic and quantum metrology, but not restricted only to these.

In this paper, we consider a fourth-order suspension bridge equation with nonlinear damping term |u_t|^m-2u_t and source term |u|^p-2u. We give necessary and sufficient condition for global existence and energy decay results without considering the relation between m and p. Moreover, when p>m, we give sufficient condition for finite time blow-up of solutions. The lower bound of the blow-up time T_max is also established. It worth to mention that our obtained results extend the recent results of Wang (J. Math. Anal. Appl., 2014) to the nonlinear damping case.

In current transportation hubs, it is a common phenomenon that passengers travel with necessary wheeled-luggage or suitcases. Due to the fact that most luggage occupies certain space within dense passenger crowds with giant mass inertia, its abnormal motions, such as drifting, frequently trig unavoidable local disturbances and turbulence in surrounding pedestrian flows, further increasing security risk. While, current researches primarily concentrate on examining the impact of luggage on crowd evacuation efficiency through scenario-based experiments and analyzing the spatial characteristics of crowd distribution in surveillance videos. There still lacks of the dynamic disturbance mechanisms and crowd stability associated with the combined motion of passengers and luggage. Therefore, this study considers the luggage-laden passenger as a deformable particle, the resulting disturbance on surrounding non-luggage passengers is analyzed and quantified into a nonlinear pressure term. Subsequently, the disturbance propagation model of passenger-luggage is developed by adapting the classical Aw-Rascle traffic flow model with incorporation of the pressure term. Simulation experiments of disturbances caused by luggage drifting and retrograding validate that the proposed propagation model aligns with real-world situations in transportation hubs. Quantitively, when luggage-laden passenger drifted left, the disturbing force of left side crowd can reach a peak of 238N in the passenger density of 3.0 p/m2, and the maximum difference between the left and right sides disturbing force can reach 153N. Furthermore, it is recommended that the proposed model can be applied to crowd flow analysis and risk assessment.

In economics, especially econometrics, fixed-effects are widely recognized as powerful tools for heterogeneity analysis of longitudinal data sets. However, the existing research pays more attention to such as the variable selection, parameter estimation and dimension reduction on linear relationships with too many assumptions while nonlinear relationships are ubiquitous. In this paper, an extension fixed-effects model is proposed on the curve_fit nonlinear least squares. First, with the help of expert knowledge, the series regressions between independent variables including time and dependent variable is estimated on each series data, and the parameters are estimated by curve_fit. Similarly, the coeffection between independent variables including time and dependent variable is also initially estimated on the whole panel data. Second, the distance between regressions is defined, and then the coeffect function between the series regressions is obtained on the initial co-effect with curve_fit. Third, the time-effects and individual fixed-effects for each series data in the panel data are estimated on the differences of series regressions and co-effect function. Additionally, the effectiveness of the proposed method is varified on the synthetic data.

This paper aims to address the exponential stability and stabilization problems for a class of delayed nonlinear Markov jump systems under randomly occurring Denial-of-Service (DoS) attack and packet loss. Firstly, the stochastic characteristics of DoS attack and packet loss are depicted by the attack success rate and packet loss rate. Secondly, a Period Observation Window (POW) method and a hybrid-input strategy are proposed to compensate for the impact of DoS attack and packet loss on the system. Thirdly, A Dynamic Event-triggered Mechanism (DETM) is introduced to save more network resources and ensure the security and reliability of systems. Then, by constructing a general common Lyapunov functional and combining with the DETM and other inequality analysis techniques, the less conservative stability and stabilization criteria for the underlying systems are derived. In the end, the effectiveness of our result is verified through a numerical example.

In this paper, we study the semilinear pseudo-parabolic equations $\displaystyle u_{t} - \triangle_{\mathbb{B}}u - \triangle_{\mathbb{B}}u_{t} = \left|u\right|^{p-1}u$ on a manifold with conical singularity, where $\triangle _{\mathbb{B}}$ is Fuchsian type Laplace operator investigated with totally characteristic degeneracy on the boundary $x_{1} = 0$. Firstly, we discuss the invariant sets and the vacuum isolating behavior of solutions with the help of a family of potential wells. Then, we derive a threshold result of existence and nonexistence of global weak solution: for the low initial energy $J(u_{0})<d$, the solution is global in time with $I(u_{0}) >0$ or $\displaystyle\Vert \nabla_{\mathbb{B}}u_{0}\Vert_{L_{2}^{\frac{n}{2}}(\mathbb{B})} = 0$ and blows up in finite time with $I(u_{0}) < 0$; for the critical initial energy $J(u_{0}) = d$, the solution is global in time with $I(u_{0}) \geq0$ and blows up in finite time with $I(u_{0}) < 0$. The decay estimate of the energy functional for the global solution and the estimates of the lifespan of local solution are also given.

In this paper, we present an innovative technique that improves the convergence order of iterative schemes that do not require the evaluation of Jacobian matrices. Using this procedure, we achieve a remarkable increase in the order of convergence, raising it from p to p + 3 units, which results in a remarkable improvement in the overall performance. We have conducted comprehensive numerical tests in various scenarios to validate the theoretical results, showing the efficiency and effectiveness of the new Jacobian-free schemes.

This work considers the existence of heteroclinic solutions for second order differential equations with ϕ-Laplacians and generalized impulsive conditions on the real line. For the construction of the results, it was only imposed that ϕ is a homeomorphism and the Schauder’s fixed point theorem allied with the concepts of L1-Carathéodory functions and sequences together with the equiconvergence on infinity and on each impulsive moments. Finally, a practical part illustrates the main theorem.

We primarily investigate the existence of solutions for fractional neutral integro-differential equations subjected to Neumann-type boundary conditions, which is crucial for understanding natural phenomena. Taking into account factors such as neutral type, fractional-order integrals, and fractional-order derivatives, we employ probability density functions, Laplace transforms, and resolvent operators to formulate a well-defined concept of a mild solution for the specified equation. Following this, by integrating fixed point theorems, we establish the existence of mild solutions under more relaxed conditions.

The study aims to explain the different effects between the multiplication, addition, division and subtraction equations. In this sense, multiplication and addition equations develop a new element or fact, meanwhile, division equations express information within the element or fact, subtraction is the selection of element or fact. Furthermore, for a multiplication and addition equation results in the new element or fact, it is necessary to depend on specific physical concepts, without extreme alteration, as it would result in the alteration of the element or fact. It is worth highlighting the importance of fixed values ​​in the equation to describe the specific nature of the element or fact. The study makes it possible to understand open equations (without solutions) in mathematics and physics, in addition to understanding existing equations. The study makes reference to the theory of obligatory necessity and the theory of differences between elements.

This paper addresses the analysis for the nonlinear vibration response of a rectangular tube with a flexible end and non-rigid acoustic boundaries. This is a further work of the linear structural acoustic problem in a well-known acoustic book. In fact, the acoustic boundaries of an enclosed space sometimes are non-rigid and the structural vibration responses are nonlinear. These two points are the focuses of this paper. The multi-level residue harmonic balance method is applied to this nonlinear structural acoustic problem. The results obtained from the multi-level residue harmonic balance method and numerical method are generally in good agreement. The effects of excitation magnitude, tube length, and phase shift parameter etc. are examined.

In the framework of the self-consistent Maxwell-Pauli theory, the non-linear Pauli equation is obtained. Stationary and nonstationary solutions of the nonlinear Pauli equation for the hydrogen atom are studied. We show that spontaneous emission and the related rearrangement of the internal structure of an atom, which is traditionally called a spontaneous transition, have a simple and natural description in the framework of classical field theory without any quantization and additional hypotheses. The behavior of the intrinsic magnetic moment (spin) of an electron wave in an external magnetic field is considered. We show that, according to the self-consistent Maxwell-Pauli theory, in a weak magnetic field, the intrinsic magnetic moment of an electron wave is always oriented parallel to the magnetic field strength vector, while in a strong magnetic field, depending on the initial orientation of the intrinsic magnetic moment, two orientations are realized: either parallel or antiparallel to the magnetic field strength vector.

In this paper, we compile the fractional power series method and the Laplace transform to design a new algorithm for solving the fractional Volterra integro-differential equation. For that, we assume the Laplace power series (LPS) solution in terms of power q=1m,m∈Z+, where the fractional derivative of order α=qγ for which γ∈Z+. This assumption will help us to write the integral, the kernel, and the nonhomogeneous terms as a LPS with the same power. The recurrence relations for finding the series coefficients can be constructed using this form. To demonstrate the algorithm's accuracy, the residual error is defined and calculated for several values of the fractional derivative. Two strongly nonlinear examples are discussed to provide the efficiency of the algorithm. The algorithm gains powerful results for this kind of problem. Under Caputo meaning the obtained results are illustrated numerically, and graphically. Geometrically, the behavior of the solution declares that the changing of the fractional derivative parameter values' in their domain alters the style of the attained solution in a symmetrical meaning and fully coinciding to the ordinary derivative value'. From these simulations, the results report that the recommended novel algorithm is a straightforward, accurate, and superb tool to generate analytic-approximate solutions for Integral, and integro-differental equations of fractional order.

The well-known equations for the process of powder compaction (PCP) in a rigid die published over the period from the beginning of the last century to the present are considered. Most of the considered equations are converted into the dependence of the densification pressure on the relative density. The equations were analyzed and their ability to describe PCP was assessed by determining the coefficient of determination when approximating experimental data on the compaction of various powders. It is shown that most of the equations contain two constants, the values of which are determined by fitting the mathematical dependence to the experimental curve. Such equations are able to describe PCP with high accuracy in the case of compaction of powders up to a relative density of 0.9-0.95. It is also shown that different equations can describe PCP in the range of density change from the initial to 0.9 with the same and sufficiently high accuracy, but when the process of compaction is extrapolated to higher values of density, the curves diverge. This indicates the importance of equations that can unambiguously describe PCP to a relative density equal to or close to 1.0. For an adequate description of PCP to a relative density greater than 0.95, equations containing three or four constants have proven themselves well.

In the paper, the authors view some ordinary differential equations and their solutions from the angle of (the generalized) derivative polynomials and simplify some known identities for the Bernoulli numbers and polynomials, the Frobenius-Euler polynomials, the Euler numbers and polynomials, in terms of the Stirling numbers of the first and second kinds.

The paper is considering of the basic differential equations of hydrodynamics: the equation of continuity and motion. On the simplest example, it is shown that the equation of continuity in a system with the equation of motion leads to contradictions and erroneous results of modeling. A more correct form of the continuity equation is described. It is shown that the equations of motion can be written in the form of complete differentials. Three possible integral forms of the equations of motion are presented. As a conclusion, the existence and smoothness of the solution of the Navier-Stokes equations are considered.

Boundary behavior of a given important function or its limit values are essential in the whole spectrum of mathematics and science. We consider some tractable cases of limit values in which either a difference of two ingredients or a difference equation is used coupled with the relevant functional equations to give rise to unexpected results. This involves the expression for the Laurent coefficients including the residue, the Kronecker limit formulas and higher order coefficients as well as the difference formed to cancel the inaccessible part, typically the Clausen functions. We also state Abelian results which yield asymptotic formulas for weighted summatory function from that for the original summatory function

Maxwell’s equations are valid only for a stationary observation point, therefore, to adequately describe real processes so far we have had to move to a moving reference frame. This paper presents the equations of electrodynamics for the moving observation point, it is shown that plane and spherical electromagnetic waves are their solutions, while the spherical wave propagates only outward, which cannot be said about Maxwell’s equations. The fields of uniformly moving charges are also solutions of the equations. Now there is no need to move to a moving reference frame, to use four-dimensional space and covariant form of equations. The question of finding a universal form of the equations that allows a solution in the form of the field of an arbitrarily moving charge remains open. This raises the question of the existence of a two-parameter group of transformations of electromagnetic fields along with the known one-parameter group has been posed. The phenomena derived from the equations, which make an additional contribution to the phase overrun in the Aharonov-Bohm effect are considered. The equation of motion of a charged particle in an electromagnetic field without simplifying approximations is considered, which allows us to take into account the radiation effects. It is shown that the fields in a moving observation point depend on its velocity and acceleration. In particular, although in a constant uniform electric field a force qE acts on a motionless charged particle, but on the same motionless but not fixed particle the force 4/3qE acts already, because it has a nonzero acceleration and the electric field at this point is larger. As the speed increases, the field decreases, and when it reaches the speed of light, when the particle stops accelerating, the force again becomes equal to qE The principle of operation of an unconventional alternator in a constant electric field and its corresponding engine, as well as new types of direct and impulse current generators, predicted by the equations, are described.

Time evolution of the quantum correlation of the hybrid qubit-qutrit system in multilocal and global dephasing under the effect of a classical noises environment is presented. We considered the qutrit in z direction and the qubit x direction. In addition, our system presents two configurations: the different and the common environment and our discussion involves three types of noises: A Static Noise (SN), Random Telegraph Noise (RTN) and Ornstein- Uhlenbeck noise (OU). By using the density matrix approach to determine the evolution of the state of the system, the found that the static noise is fatal for the quantum entanglement, the quantum correlation and for RTN and OU noises the quantum entanglement, the quantum correlation decreases monotonically with a sufficiently long time. Therefore, we can say that, entanglement and the quantum correlation are more robust to the noise precisely those Random telegraph noise and Ornstein-Uhlenbeck noise.

In a bid to resolve lingering problems in cosmology, more focus is being tilted towards cosmological models in which physical constants of nature are not necessarily real constants, but varying with cosmic time. In this paper we have study cosmology in nonlinear electrodynamics with the Newton's gravitational constant $G$ not a constant but vary in form of power-law of the scale factor of the universe. The evolution of the scale factor $a (t)$ is studied in this model which depends on nonlinear electrodynamics fine tuning term of $\alpha$. Then we check the stability of the model using the speed of sound.

We discuss methods for calculation of critical indices and amplitudes from the perturbative expansions. Several important examples of the Stokes flow through 2D and 3D channels are brought up. Power series for the permeability derived for small values of amplitude are employed to calculation of various critical exponents in the regime of large amplitudes. Special nonlinear approximations valid for arbitrary values of the wave amplitude are derived from the expansions. The technique developed for critical phenomena is applied then for relaxation phenomena. The concept of time-translation invariance is discussed, its spontaneous violation and restoration considered. Emerging probabilistic patterns correspond to a local breakdown of time-translation invariance. Their evolution leads to the timetranslation symmetry complete (or partial) restoration. We estimate typical time extent, amplitude and direction for such restorative process. The new technique is based on explicit introduction of origin in time. After some transformations we come to the exponential and generalized, exponential-type solution with explicit finite time scale, which was only implicit in initial parametrization with polynomial approximation. The concept of crash as a relaxation phenomenon, consisting of time-translation invariance breaking and restoration, is put forward. COVID-19 related mini-crash in the time series for Shanghai Composite is discussed as an illustration.

In this paper, we report orthogonal fuzzy versions of some celebrated iterative mappings. We provide various concrete conditions on the real valued functions J,S:(0,1]→(−∞,∞) for the existence of fixed-points of (J,S)-fuzzy iterative mappings. We obtain many fixed point theorems in orthogonal fuzzy metric spaces. We apply (J,S)-fuzzy version of Banach fixed point theorem to show the existence and uniqueness of the solution. We support these results with several non-trivial examples and applications to Volterra-type integral equations and fractional differential equations.

This paper considers a nonlinear model for population dynamics with age structure. The fertility rate with respect to age is non constant and has the form proposed by . Moreover, its multiplicative structure and the multiplicative structure of mortality makes the model separable. In this setting it is shown that the number of births in unit time is given by a system of nonlinear ordinary differential equations. The steady state solution together with the equilibrium solution is found explicitly.

Forest resources are renewable,rational exploitation and utilization of forest resources are not only conducive to sustainable development of population, but also can obtain higher economic benefits. Based on actual timber harvesting problems,a combined harvesting model of timber and non-timber with nonlinear harvesting term was established.Firstly, using the qualitative theory of ordinary differential equations, the dynamic behavior of the model is studied, and the existence and stability of the equilibrium point of the model are proved .Secondly, the optimal control solution is obtained by using the optimal control theory.Finally, the optimal harvesting strategy of timber and non-timber products is given based on the numerical simulation results, which provides a certain theoretical basis for the sustainable development of forestry ecological economy.

The paper examines the impact of inward remittances flows on per capita GDP growth in Bangladesh during 1976-2012. We find that the growth effect of remittances is negative at first but becomes positive at a later stage, an evidence of a non-linear. Unproductive use of remittances was rampant in the beginning when they were received by migrant families but better social and economic investments led to more productive utilisation of remittances receipts at later periods. This was the possible mechanism behind the U-shaped relationship. Unlike what is suggested in the literature that the effect of remittances is more pronounced in a less financially developed economy, our evidence do not show that the effect of remittances on per capita GDP growth in Bangladesh is conditional on the level of financial development.

The ordinary continuity equation relating the current and density of a system is extended to incorporate systems with dual (longitudinal and transverse) currents. Such a system of equations is found to have the same mathematical structure as that of Maxwell equations. The horizontal and transverse currents and the densities associated with them are found to be coupled to each other. Each of these quantities are found to obey a wave equation traveling at the velocity of light in vacuum. London's equations of super-conductivity are shown to emerge from some sort of continuity equations. The new London's equations are symmetric and are shown to be dual to each other. It is shown that London's equations are Maxwell's equations with massive electromagnetic field (photon). These equations preserve the gauge invariance that is broken in other massive electrodynamics. The duality invariance may allow magnetic monopoles to be present inside superconductors. The new duality is called the comprehensive duality transformation.

The second article in the series presents the application of the Smoothed Particle Hydrodynamics (SPH) method to modelling of batch crystallisation in stirred tanks. A methodology to integrate the population balance equations (PBE) in parallel and independently from the Navier-Stokes equations is demonstrated. The benefits of the proposed methodology in terms of computational requirements, accuracy and availability of the crystal size distribution are discussed. The specific formulation of the SPH equations where the resulting system of ordinary differential equations is solved using the weighted contributions rather than numerically by solving a linear system of equations allows for massive parallelisation and a very loose coupling of the population balance and the fluid dynamics. It has been demonstrated, that the population balance equations can be solved on a Shared Memory Architecture (SMA) system using the OpenMP interface while the fluid dynamics equations being computed independently on a General Purpose Graphics Processing Unit (GPGPU) using the NVidia CUDA technology. This way, a significant portion of the computational overhead due to the large number of additional transport equations resulting from the discretisation of the population balance was removed: the SPH simulation coupled with 200 population balance equations was only 40% slower compared to SPH-only simulation. Two methods for the solution of population balance equations that preserve full crystal size distribution were implemented: discretised population balance (DPB) and method of characteristics (MOCH). The DPB equations are solved using the high-resolution finite-volume method with flux limiter and the effect of a large number of different flux limiters have been investigated. Both methods were validated using the case studies from the literature where an analytical solution can be derived. The developed models were applied to a numerical solution of coupled computational fluid dynamics and population balance equations to model a batch crystallization process. The effect of the hydrodynamics on the local temperature/supersaturation and the resulting crystal size distribution was captured and compared to the ideal mixing case. The simulation results from the DPB and MOCH methods were compared in terms of computational requirements and accuracy and MOCH selected as computationally more efficient and accurate.

Imperfection sensitivity of large amplitude vibration of curved single-walled carbon nanotubes (SWCNTs) is considered in this study. The SWCNT is modeled as a Timoshenko nano-beam and its curved shape is included as an initial geometric imperfection term in the displacement field. Geometric nonlinearities of von Kármán type and nonlocal elasticity theory of Eringen are employed to derive governing equations of motion of nano-beam. Spatial discretization of governing equations and associated boundary conditions is performed using differential quadrature (DQ) method and the corresponding nonlinear eigenvalue problem is iteratively solved. Effects of value and location of the geometric imperfection, and the nonlocal small-scale parameter on the nonlinear frequency ratio and imperfection sensitivity of a curved SWCNT for various boundary conditions are investigated. The results show that the geometric imperfection plays a significant role in the nonlinear vibration characteristics of curved SWCNTs.

In this paper, we study the problem of controllability of impulsive neutral evolution integrodifferential equations with state dependent delay in Banach spaces. The main results are completely new and are obtained by using Sadovskii's fixed point theorem, theory of resolvent operators, and an abstract phase space. An example is given to illustrate the theory.

The objective of this study is to evaluate the effect of constitutive equations on the prediction accuracy for springback in cold stamping with various deformation modes. In this study, two types of yield functions—Hill’48 and Yld2000-2d—were considered to describe yield behavior. Isotropic and kinematic hardening models based on the Yoshida–Uemori model were also adopted to describe hardening behavior. Various material tests (such as uniaxial tension, tension- compression, loading-unloading, and hydraulic bulging tests) were carried out to determine the material parameters of the models. The obtained parameters were implemented in the finite element (FE) simulation to predict springback, and the results were compared with experimental data. U-bending and T-shape drawing were employed to evaluate the springback prediction accuracy. Obviously, the springback prediction accuracy was greatly influenced by constitutive equations. Therefore, it is important to choose appropriate constitutive equations for accurate description of material behaviors in FE simulation.

We show that one of the two important uncertainty principles derived by Maccone and Pati \textit{[Phys. Rev. Lett., 2014]} can be derived for arbitrary maps defined on subsets of $\mathcal{L}^p$ spaces for $1< p<\infty$. Our main tool is the Clarkson inequalities. We also derive a nonlinear uncertainty principle for weak parallelogram spaces and Type-p Banach spaces.

According to our current understanding of quantum mechanics, a `measurement' violates unitarity. In other words as the act of measurement `resets' the probabilities described by the Schrödinger equation, in the very `moment' of the collapse of the wavefunction, conservation of probability does not hold. To make progress in our understanding of the measurement problem it is thus expected that one must encompass in a single equation both otherwise mutually-exclusive processes of measurement and unitary evolution. In this light, turning to the continuity equation, we realize the possibility that if we allow for existence of sources and sinks in the `fluid of probability' we get closer to such a theory by arriving at nonlinear generalizations of Schrödinger and Klein-Gordon equations. The generalized equations derived are then shown to satisfy all conditions that are expected from a proper generalization: simplification to their linear counterparts by a well-defined dynamical condition, and Galilean and Lorentz invariance.

Continuation methods are presented that are capable of treating frequency-domain flutter equations including multiple nonlinearities represented by describing functions. A small problem demonstrates how a series of continuation processes can find all limit-cycle oscillations within a specified region with a reasonable degree of confidence. Curves of the limit-cycle amplitude variation with velocity, indicating regions of stability and instability with colors give a compact view of the nonlinear behavior throughout the flight regime. A continuation technique for reducing limit-cycle amplitudes by adjusting various system parameters is presented. These processes are economical enough to be a routine part of aircraft design and certification.

We investigated the dynamics of a simple master equation describing the evolution of a probability distribution between two states, \( p_1 \) and \( p_2 \) (with \( p_1 + p_2 = 1 \)). Our focus is on the behavior of the entropy \( S \), the distance \( D \) from the uniform distribution (\( p_1 = p_2 = 1/2 \)), and the free energy \( F \). We define two ratios, \( \eta_S = \frac{dS/dt}{dF/dt} \) and \( \eta_D = \frac{dD/dt}{dF/dt} \), and find that both diverge to plus-minus infinity as the system evolves towards the uniform distribution. This divergence marks a critical transition point where the master equation temporarily loses physical meaning. However, after this divergence, the system settles into a new, well-behaved regime where the solutions stabilize at finite values, representing a new equilibrium or steady state. This two-regime behavior underscores the intricate dynamics of simple probabilistic systems as they approach and surpass equilibrium, providing insights into the complex interplay of entropy, distance, and free energy in statistical mechanics.

The research paper provides a mathematical framework for understanding phase shift in wave phenomena, bridging theoretical foundations with real-world applications. It emphasizes the importance of phase shift in physics and engineering, particularly in fields like telecommunications and acoustics. Key equations are introduced to explain phase angle, time delay, frequency, and wavelength relationships. The study also introduces the concept of time distortion due to a 1° phase shift, crucial for precise time measurements in precision instruments. The research also addresses infinitesimal wave energy loss related to phase shift, enriching our understanding of wave behavior, and impacting scientific and engineering disciplines.

Electrically charged particles such as Electrons and Protons carry electric, E, and magnetic, B, fields. In addition to these fields, Quantum Mechanics (QM) endows these particles with an `arcane and spooky' field --- the wavefunction. This wavefunction of QM is not only assumed to be separate but distinct from the electromagnetic field. We herein upend this view by demonstrating otherwise. That is, we demonstrate that the four components of the Dirac wavefunction, can be shown to not only be an intimate, but, a direct function of the electromagnetic field carried by the particle in question. Insofar as unity, depth in our understanding and insight into both Dirac and Maxwell's equations as major pillars of Modern Physics, we believe that this work may very well inch us one-step-closer to the truth.

The principles of causality and unitarity are studied within rational nonlinear electrodynamics proposed earlier. We investigate dyonic and magnetized black holes and show that in the self-dual case, when the electric charge equals the magnetic charge, corrections to Coulomb's law and Reissner-Nordstrom solutions are absent. In the case of the magnetic black hole, the Hawking temperature, the heat capacity and the Helmholtz free energy are calculated. It is shown that there are second-order phase transitions and it was demonstrated that at some range of parameters the black holes are stable.

We explored an extension of the Cauchy-Riemann equations to the algebra of treons, recently described by Alejandro Bermejo, whose elements are ordered 3-tuples. We leveraged the isomorphism between the algebra of treons and algebra B, and deduced the Cauchy-Riemann equations for the algebra of treons, establishing the necessary conditions for analyticity in this algebraic structure. This work significantly broadened our horizons in complex analysis and introduced new possibilities for applications across various fields of advanced mathematics.

The paper considers the physical consequences of bimetric geodesically (pro jectively) invariant equations of gravitation. The correct expression for the energy- impulse tensor of the gravitational field is obtained and with its help the key role of the gravitational defect of mass-energy for the energy of the newborn universe is shown..

Most fully developed galaxies have vivid spiral structure, but the formation and evolution of spiral structure is still a mystery that is not fully understood in astrophysics. We find that the currently used equations of galactic dynamics contain some unreasonable components. In this paper, the following three working assumptions are introduced to simplify the galactic structural equations. 1. In the research of large-scale structure, the retarded potential of the gravitational field should be taken into account. The propagating time of the gravitational field from center to border is longer than the revolution periods of the stars near the center of galaxy. Newton's gravitational potential is unreasonable for such case, and the weak field and low velocity approximation of Einstein's field equation should be adopted. 2. The stars in a fully developed galaxy should be zero-pressure and inviscid fluid, and the equation of motion is different from that of ordinary continuum mechanics. Stars move along geodesics. 3. The structure of the galaxy is only related to the total mass density distribution. The equation of state of dark halo is different from that of ordinary luminous interstellar matter, so their trajectories are also very different. Dark halo and ordinary matter in galaxy are automatically separated. The total mass density distribution can be presupposed according to the observation data, and then it can be determined by comparing the solution of the equations with the observed data. These assumptions and treatments are supported by theory and observation. The variables of the equations of simplified galactic dynamics are separated from each other, and the equations are well-posed and can be solved according to a definite procedure. Therefore, this simplified dynamic equation system provides a more reasonable and practical framework for the further study of galactic structure, and can solve many practical problems. Besides, it is closely related to the study of dark energy and dark matter.

In the paper, the authors present an explicit form for a family of inhomogeneous linear ordinary differential equations, find a more significant expression for all derivatives of a function related to the solution to the family of inhomogeneous linear ordinary differential equations in terms of the Lerch transcendent, establish an explicit formula for computing all derivatives of the solution to the family of inhomogeneous linear ordinary differential equations, acquire the absolute monotonicity, complete monotonicity, the Bernstein function property of several functions related to the solution to the family of inhomogeneous linear ordinary differential equations, discover a diagonal recurrence relation of the Stirling numbers of the first kind, and derive an inequality for bounding the logarithmic function.

We show that Universe inflation occurs when Einstein's gravity couples to nonlinear electromagnetic fields with cosmic stochastic magnetic fields background. Nonlinear electrodynamics with two parameters is used. The strength of magnetic fields to have the Universe inflation is obtained. It is shown that singularities of the energy density and pressure are absent for any scale factors. At large scale factors one arrives at equation of state for ultra-relativistic case. It is demonstrated that the curvature invariants do not have singularities. By computing the deceleration parameter we show that graceful exit takes place with reasonable e-folding number. The duration of Universe inflation as a function of model parameters is obtained. The classical stability and causality are analyzed by calculating the speed of the sound.

We derive an uncertainty principle for Lipschitz maps acting on subsets of Banach spaces. We show that this nonlinear uncertainty principle reduces to the Heisenberg-Robertson-Schrodinger uncertainty principle for linear operators acting on Hilbert spaces.

The article analyzes existing materials and structures with quadratic-nonlinear optical properties that can be used to generate a difference frequency in the terahertz and sub-terahertz frequency ranges. The principle of constructing a nonlinear optical-radio converter, based on an optical focon (a focusing cone), is proposed. Based on the assumption that this focon can be implemented from the metal-organic framework (MOF), we propose a technique for modeling its parameters. The mathematical model of the process of propagation and nonlinear interaction of waves inside the focon is based on a simplification of the nonlinear wave equation. Within the framework of the developed model, the following parameters are approximately determined: the 3D gradient of the linear refractive index and the function determining the geometric profile of the focon, which provide a few-mode-based generation of the difference frequency.

The nonlinear deformation and stress analysis of a circular torus is a difficult undertaking due to its complicated topology and the variation of the Gauss curvature. A nonlinear deformation (only one term in strain is omitted) of Mindlin torus was formulated in terms of the generalized displacement, and a general Maple code was written for numerical simulations. Numerical investigations show that the results obtained by nonlinear Mindlin, linear Mindlin, nonlinear Kirchhoff-Love, and linear Kirchhoff-Love models are close to each other. The study further reveals that the linear Kirchhoff-Love modeling of the circular torus gives good accuracy and provides assurance that the nonlinear deformation and stress analysis (not dynamics) of a Mindlin torus can be replaced by a simpler formulation, such as a linear Kirchhoff-Love theory of the torus, which has not been reported in the literature.

Based on the nonlinear failure criterion and modified tangential technique, the upper bound solutions of supporting pressure on the deep tunnel face were obtained under pore water pressure conditions. The influence of parameters on the supporting pressure and collapse range was investigated according to the unlimited block failure mechanism. It is found that the upper bound solutions of supporting pressure increase with the growth of the nonlinear coefficient and pore water pressure coefficient. The collapse range of the tunnel face scales out with the increase of the nonlinear coefficient and shrinks when increasing the pore water pressure coefficient. Moreover, with the increase of the nonlinear coefficient, the impact strength on supporting pressure and collapse range declines gradually. According to the calculating results, it is manifest that both the pore water pressure and nonlinear criterion factors have negative impacts on the stability of the tunnel face. Thus, more attention should be paid to these parameters to ensure face stability in deep tunnel construction.

We consider a Cauchy problem for differential equations in a Hilbert space $X$. The problem is stated in a time interval $I$, which can be finite or infinite. Weuse a fixed point argument for history-dependent operators to prove the unique solvability of the problem. Then, we state and prove convergence criteria for both a general fixed point problem and thecorresponding Cauchy problem. These criteria provide necessary and sufficient conditions on a sequence $\{u_n\}$ which guarantee its convergence to the solution of the corresponding problem, in the space of both continuous and continuously differentiable functions. We then specify our results in the study of a particular differential equation governed by two nonlinear operators. Finally, we provide an application in viscoelasticity and give a mechanical interpretation of the corresponding convergence result.

A careful examination of the fundamentals of electromagnetic theory shows that due to the underlying mathematical assumptions required for Stokes' Theorem, global charge conservation cannot be guaranteed in topologically non-trivial spacetimes. However, in order to break the charge conservation mechanism we must also allow the electromagnetic excitation fields D, H to possess a gauge freedom, just as the electromagnetic scalar and vector potentials $\varphi$ and A do. This has implications for the treatment of electromagnetism in spacetimes where black holes both form and then evaporate, as well as extending the possibilities for treating vacuum polarisation. Using this gauge freedom of D, H we also propose an alternative to the accepted notion that a charge passing through a wormhole necessarily leads to an additional (effective) charge on the wormhole's mouth.

We explore the field equations in a 4-d complex space-time, in the same way, that general relativity does for our usual 4-d real space-time, forming this way, a new "general relativity" in C⁴ space-time, free of sources. Afterwards, by embedding our usual 4-d real space-time in C⁴ space-time, we describe geometrically the energy-momentum tensor T_μν as the lost geometric information of this embedding. We further give possible explanation of dark eld and dark energy.

This research aims to investigate psychological anxiety factors, particularly Music Performance Anxiety (MPA), that influence the sustainability of music learning systems. This study employs an innovative integration of qualitative and quantitative methodologies, marking the first use of a nonlinear system paired with a comprehensive data framework for analyzing questionnaire responses. This approach allows for a detailed examination of the effects of these factors on academic outcomes. Previous research has primarily focused on the fragmented music learning system from the perspective of teachers' instructional strategies, emphasizing phased teaching methods and learning objectives which focus often neglected a significant subset of music learners. Addressing this oversight, the current study specifically centers on the emotional aspects of anxiety within music learning systems, with a particular emphasis on students, who play a crucial role in music education. This perspective enables a deep exploration of how MPA and other related emotional factors affect music learning systems. A notable innovation of this research is the development of a stochastic nonlinear system designed to model the psychological factors impacting student music learning, complemented by a quantitative model that assesses actual learning outcomes. By analyzing this nonlinear system, the research identifies weight values for various negative emotions, particularly those related to MPA, that significantly influence the music learning process. The integration of MPA and other anxiety related emotional factors with practical learning outcomes provides a comprehensive understanding of their combined impact on the music learning system. These insights are invaluable for educators and policymakers aiming to enhance both the effectiveness and emotional well-being within music education. This comprehensive approach offers a novel perspective in understanding and improving the dynamics of music education through a more empathetic and scientifically grounded lens.

When discussing the Diophantine quintic equation $p(a^5 + b^5) = q(c^5 + d^5)$ where p is a prime and q is an integer, there is a clear gap between the mathematical literature and online forums. Because of the intrinsic complexity of parameterizing fifth-degree equations, this equation re- mains largely unexplored. In this work, we approach this quintic problem through algebraic methods with the goal of providing numerical solutions that illuminate its properties and behaviors.Our research uncovers in- triguing patterns and connections that shed light on the enigmatic nature of Diophantine quintic equations in this particular form.

Coupling the Maxwell tensor to the Riemann-Christoffel curvature tensor is shown to lead to a geometricized theory of electrodynamics. While this geometricized theory leads directly to the classical Maxwell equations, it also extends their physical interpretation by giving charge density, mass density and the four-velocity that describes their motions geometric definitions. Assuming the existence of a conserved energy-momentum tensor, all solutions to the geometricized theory of electrodynamics developed here are shown to be consistent with the emergence of gravity obeying the General Relativity field equation augmented by a term that mimics the properties of dark matter and/or dark energy. Finally, due to the symmetries of the theory, the properties and phenomenology of antimatter emerge in solutions.

In this paper, we investigate the equilibrium points of following a system of difference equations x_n+1 = x_n 2y_n−1, y _n+1 = y_n 2x_n−1. We also study the asymptotic stability of related system of difference equations. Further we examine the periodic solutions of related system with period two. Additionally, we find out the invariant interval and periodic cycles of related system of difference equations.

We have shown that the use of evolutionary gravitational constants G, the speed of light c, and cosmological constant Λ in the Einstein’s field equations leads to a very simple model that fits the supernovae Ia data with a single parameter as well as the standard ΛCDM model with two parameters, and has the predictive capability superior to the latter. We have shown that at the current time G and c both increase as dG/dt = 5.4GH₀ and dc/dt = 1.8cH₀ with H₀ as the Hubble constant, and Λ decreases as dΛ/dt = −1.2ΛH₀. Negative finding on the variation of G with time derived from the analysis of lunar laser ranging data in several studies has been attributed to the fact that these studies did not consider the variation of c simultaneously with G.

The widespread integration and growing dependence among currency, stock and commodity markets make these markets often very vulnerable to shocks and at risk of collapse at the same time. As a result, these trends threaten the sustainability of the entire financial system. In this study, we explore the spillovers and nonlinear dependences between the seven major foreign ex-change rates, crude oil and gold prices, a global stock price index, and oil and stock implied volatility indices as proxy variables for global risk factors by em-ploying directional spillover network approach. We also use multi-scale de-composition method and nonlinear causality test between these variables to capture multi-level relationships at short and long horizons. Major findings are summarized as follows. First, from the multi-scale decomposition analysis, we identify that Granger causality test results and the direction and strength of return spillovers change with the level of decomposition. Second, the results of nonlinear causality tests show variation in both the significance and direction of Granger causality relationships between the decomposed currency and other series at different timescales, especially for the decomposed oil, gold, and OVX series. Third, the measured directional spillover indices identify the EUR as the largest contributor of connectedness to the other series. The central role of the EUR is a net transmitter of connectedness to gold, oil, the GBP, JPY, and CHF.

The article is devoted to issues related to the propagation and transformation of vortexes in the optical range of frequency. Within the framework of the traditional and modified model of slowly varying envelope approximation (SVEA), the process of converting vortex beams of the optical domain into vortex beams of the terahertz radio range based on nonlinear generation of a difference frequency in a medium with a second-order susceptibility is considered. The modified SVEA splits a slowly varying amplitude into two factors, which makes it possible to more accurately describe the three-wave mixing process. The theoretical substantiation of the rule of vortex beams topological charges conversion is given – the topological charge of the output radio-vortex beam is equal to the difference between the topological charges of the input optical vortex beams. A numerical simulation model of the processes under consideration has been implemented and analyzed.

We prove general stability theorems for $n$-dimensional quartic-cubic-quadratic-additive type functional equations of the form \begin{eqnarray*} \sum_{i=1}^\ell c_i f \big( a_{i1}x_1 + a_{i2}x_2 + \cdots + a_{in}x_n \big) = 0 \end{eqnarray*} by applying the direct method. These stability theorems can save us much trouble of proving the stability of relevant solutions repeatedly appearing in the stability problems for various functional equations.

The particulars of stimulus-response experiments done on dynamic biosystems clearly limit what one can learn and validate about their structural interconnectivity (topology), even when collected kinetic output data are perfect (noise-free). As always, available access ports and other data limitations rule. For linear systems, exponential modes, visible and hidden, play an important role in understanding data limitations, embodied in what we call dynamical signatures in the data. We show here how to circumscribe and analyze modal response data in compartmentalizing model structures – so that modal analysis can be used constructively in systems biology model building – for nonlinear (NL) as well as linear biosystems. We do this by developing and exploiting the modal basis for dynamical signatures in hypothetical (perfect) input-output data associated with a structural model – one that includes inputs and outputs explicitly – and for NL as well as linear biosystems. The methodology establishes model dimensionality (size, complexity) from particular data sets; helps select among multiple candidate models (model distinguishability); helps in designing new experiments to extract “hidden” structure; and helps to simplify (reduce) models to their essentials. For NL biosystems, results are not as comprehensive, similarly informative about their dominant dynamical properties, and unified with linear models on invariant 2-dimensional manifolds in phase space. Some automation of these highly technical aspects of biomodeling is also introduced.

An algorithm for the treatment of images affected by both blurring and salt&pepper noise is proposed with a cost only proportional to the number of pixels. The methodology uses a discretization scheme for the Laplace operator multiplied by a suitable nonlinear term depending on the gradient. Even if this approach resembles a diffusion type algorithm, only one step of the procedure is applied, leading to significant time savings. The procedure is successfully tested on some standard black&white natural images.

DGJ (Daftardar-Gejii-Jafaris) method is used to obtain numerical solution of the third order fractional differential equation. Providing the DGJ method converges, the approximate solution is a good and effective numerical result which is close to the exact solution or the exact solution. For this,the examples of the explaning the method are presented. The proposed method is implemented for the approximation solution of the third order nonlinear fractional partial differential equations. The method was shown to be unsuitable and inconsistent for an example of a nonlinear fractional partial differential equation depend on initial-boundary value conditions. The fact that these numerical results are not consistent can be explained by the fact that the method is not convergent.

Reliable fatigue design rules affect the proactive identification of safety parameters in aerospace industry. Numerous fatigue crack initiation and propagation models, linear and nonlinear, have been developed for designing purposes or estimation of the remaining life of aging airplanes. Depending on the adopted assumptions, the accuracy varies for different loading histories, loading types, and materials. Semi empirical models are simple but yield significant inaccuracies. Models with better theoretical basis provide better accuracy, but implementation in real conditions is problematic. In the present work, a review of author’s recent fatigue crack initiation and propagation models based on physical mechanisms is presented and improvements are proposed. Verification of the models on test results is provided and discussed.

Our starting point is the covariant coordinate transformation equation of a relativistic positioning system in Minkowski space-time, which maps the receiver’s emission coordinates (proper times broadcast by the emitters) to its coordinates in some inertial reference frame. Bancroft’s analytical (closed-form) solution to the basic pseudorange navigation equations with four emitters is recovered and the subjacent elements are geometrically interpreted. The case of four static beacons is analized as a clarifying situation.

It is universally accepted that Maxwell equations do not remain invariant under the Galilean transformation. This conflicts the principle of relativity which states that the physical law must remain invariant in the mathematical form in all inertial frames of reference. For this reason, the Lorentz transformation is invented, and the Galilean transformation is nowadays superseded by the Lorentz transformation. However, this paper challenges this widely held belief that the Maxwell equations are not invariant under the Galilean transformation. By applying the Galilean transformation to Lienard-Wiechert electromagnetic fields it is mathematically proven that the Maxwell equations indeed remain invariant under the Galilean transformation. In addition, the critical error in Lorentz's proof of Galilean non-invariance of Maxwell equations is pointed out, and it turns out that Lorentz's conclusion that Maxwell equations are not invariant under the Galilean transformation is the result of a mathematical error.

We establish certain modular equations of degree three based on some results of Ramanujan's theta functions.

Considerable efforts have been devoted to modifying gravity, which aim to elucidate the possible existence or nature of dark matter and energy, describe observational data more effectively, and formulate quantum gravity. In addition, despite the immense success of the quantum field theory, the framework requires renormalization techniques and breaks down at high energies. Recently, the Planck legacy 2018 release has confirmed the existence of an enhanced lensing amplitude in the cosmic microwave background power spectra, which prefers a positively curved early Universe with a confidence level higher than 99%. This study considers the implied curvature of the early Universe as the curvature of “the background or 4D conformal bulk as a manifestation of vacuum energy” and distinguishes it from the localized curvature that is induced in the bulk by the presence of celestial objects that are regarded as ”4D relativistic cloud-worlds.” Analogously, because gravity appears to emerge owing to spacetime curvature and does not exhibit critical characteristics shared by other fields, it has been incorporated as the local curvature of the bulk affecting the embedded quantum fields that are regarded as propagating “4D relativistic quantum clouds.” To consider the effects of the bulk on embedded clouds, this paper presents interaction field equations in terms of brane-world modified gravity and the perspective of geometrization of the quantum mechanics wherein gravity is manifested by the curvature of 4D conformal bulk as an indicator of the field strength of vacuum energy on the embedded 4D relativistic clouds in addition to the boundary interactions, which could remove the singularities and satisfy a conformal invariance theory. A visualization of the evolution of the 4D relativistic cloud-worlds over the conformal spacetime of the 4D bulk is presented, whereas the standard said theories can be recovered from the interaction field equations.