Stability Of Equilibrium Point Research Articles

Analysis of a Resistive-Capacitive Shunted Josephson Junction with Topologically Nontrivial Barrier Coupledto a RLC Resonator

This paper presents the resistive-capacitive shunted Josephson junction (RCSJJ) with a topologically nontrivial barrier (TNB) coupled to a linear RLC resonator. The rate equations describing RCSJJ with TNB coupled to the linear RLC resonator are established via Kirchhoff’s current and voltage laws. The model exhibits four, two, or no equilibrium points depending on the external direct current (DC) source and the fractional parameter. The stability analysis of the equilibrium points with credit to the Routh-Hurwitz stability criterion reveals that the stability of equilibrium points depends on the DC source and the fractional parameter. Current-voltage characteristic reveals the presence of a birhythmicity zone which is sensitive to the fractional parameter m. As the fractional parameter increases, the coexistence of the resonant state is destroyed, which is followed simultaneously by the appearance of a new resonance state. Depending on initial conditions, birhythmic behaviour is characterized by the existence of a limit cycle. The projection of the phase space in the specific plane and the time evolution of charge is predicted in which the amplitude of attractors reported is sensitive to the parameter m. Lastly, with a defined fractional parameter, the amplitude of the branch locked to the resonator is greater than the unlocked branch.

When and why microbial-explicit soil organic carbon models can be unstable

Abstract. Microbial-explicit soil organic carbon (SOC) cycling models are increasingly being recognized for their advantages over linear models in describing SOC dynamics. These models are known to exhibit oscillations, but it is not clear when they yield stable vs. unstable equilibrium points (EPs) – i.e., EPs that exist analytically but are not stable in relation to small perturbations and cannot be reached by transient simulations. The occurrence of such unstable EPs can lead to unexpected model behavior in transient simulations or unrealistic predictions of steady-state soil organic carbon (SOC) stocks. Here, we ask when and why unstable EPs can occur in an archetypal microbial-explicit model (representing SOC, dissolved OC (DOC), microbial biomass, and extracellular enzymes) and some simplified versions of it. Further, if a model formulation allows for physically meaningful but unstable EPs, can we find constraints in the model parameters (i.e., environmental conditions and microbial traits) that ensure stability of the EPs? We use analytical, numerical, and descriptive tools to answer these questions. We found that instability can occur when the resupply of a growth substrate (DOC) is (via a positive feedback loop) dependent on its abundance. We identified a conservative, sufficient condition in terms of model parameters to ensure the stability of EPs. Principally, three distinct strategies can avoid instability: (1) neglecting explicit DOC dynamics, (2) biomass-independent uptake rate, or (3) correlation between parameter values to obey the stability criterion. While the first two approaches simplify some mechanistic processes, the third approach points to the interactive effects of environmental conditions and parameters describing microbial physiology, highlighting the relevance of basic ecological principles for the avoidance of unrealistic (i.e., unstable) simulation outcomes. These insights can help to improve the applicability of microbial-explicit models, aid our understanding of the dynamics of these models, and highlight the relation between mathematical requirements and (in silico) microbial ecology.

Enhancing Water Pollution Management Through a Comprehensive Fractional Modeling Framework and Optimal Control Techniques

Water pollution poses a significant threat to public health, particularly in developing nations. This study aims to enhance our understanding of water pollution transmission dynamics by introducing a novel mathematical model within a general fractional framework. The model offers a comprehensive analysis, including assessing solution feasibility and the stability of equilibrium points. To effectively implement this model, we develop an efficient numerical scheme based on the trapezoidal method. Furthermore, we conduct a thorough error discussion and convergence analysis. Notably, we devise and evaluate a new strategy to effectively tackle the associated fractional optimal control problem concerning water pollutant transmission. Applying this approach to a real case demonstrates its substantial impact, notably in reducing non-soluble contaminants within the water transmission system.

Generating grid double-scroll attractors from magnetized SC-CNN and its application in image encryption

Cellular neural network has attracted widespread attention due to its easy hardware implementation and general application in information processing. In this paper, a new memristive State-Controlled-CNN (MSC-CNN) model is established by introducing two novel flux-controlled memristors into a three-cell network, where one of which simulates the electromagnetic radiation and the other characterizes the coupling weight. The present MSC-CNN can generate multi-double-scroll attractors (MDSAs) in one direction or grid multi-double-scroll attractors (GMDSAs). Stability analysis shows that each asymptotically stable equilibrium point of the memristors serves as index-1 saddle-focus in the MSC-CNN, which induces the corresponding connecting bond orbits between two scrolls. Thus, the number of GMDSAs can be flexibly adjusted by varying the number of the asymptotically stable equilibrium points of the memristors. Moreover, the memristive coupling strength also controls the number of MDSAs. Through numerical simulations, complex dynamic behaviors such as extreme multistability and amplitude control are also observed. Then, an image encryption scheme based on the proposed MSC-CNN is designed from the application perspective. Finally, a microcontroller-based hardware system is implemented to verify the effectiveness of numerical simulations.

Hysteresis behavior and generalized Hopf bifurcation in a three-degrees-of-freedom aeroelastic system with concentrated nonlinearities

A three-degrees-of-freedom aeroelastic system with concentrated structural nonlinearities is considered. The aerodynamic loads in the model are simulated by using unsteady aerodynamic forces. First, we evaluate the linear stability region of the aerodynamic system employing the Lienard-Chipart criterion so that the stability boundaries of the center of gravity position and the uncoupled pitch natural frequency are obtained. Next, we present the distribution of the first Lyapunov coefficient l1 in the parameter space by using the Poincaré projection method to determine the type of Hopf bifurcation (subcritical or supercritical). We identify three types of hysteresis loops associated with Hopf bifurcations: the first type contains one stable equilibrium point and one stable limit cycle, the second type consists of one stable equilibrium point and two stable limit cycles, and the third type comprises two stable limit cycles. Nonlinear analysis shows that the variations in hysteresis loop types are mainly caused by the nonlinear flutter critical points, namely the saddle–node bifurcation points of limit cycles. Furthermore, we determine a degenerate Hopf bifurcation type known as generalized Hopf bifurcation (Bautin bifurcation) and its corresponding second Lyapunov coefficient, l2. Subsequently, we conduct nonlinear flutter analyses for positive and negative values of l2 by using bifurcation diagrams and phase portraits. The results show that when the uncoupled pitch natural frequency is used as the control parameter, the occurring generalized Hopf bifurcations are all supercritical (l2<0). However, when the center of gravity position is used as the control parameter, the generalized Hopf bifurcations that occur have both supercritical and subcritical (l2>0) modes.

Study on the interior equilibrium point of a special class of 2 × 2 × 2 asymmetric evolutionary games.

Many behavioural interactions in real life involve three individuals. When each individual has two alternative strategies, they can be abstracted into mathematical models by means of asymmetric games. In this paper, we explore a special class of asymmetric games satisfying fixed conditions. Firstly, we analyse two solitary interior equilibrium points and provide the judgement condition for their instability based on the Jacobi matrix local stability analysis method. Secondly, we analyse the interior equilibrium points that are continuously distributed within a line and probe into their stability conditions based on generalized Hamiltonian systems theory. Under the circumstances, the stable interior equilibrium point is surrounded by closed orbits in phase space, which presents an observable stable state where two strategies coexist and fluctuate in each of the three game populations. This work enriches the study of asymmetric games' evolutionary dynamics.

Modeling and analysis of a human papilloma virus transmission model with impact of media

The human papillomavirus (HPV) is threatening human health as it spreads globally in varying degrees. On the other hand, the speed and scope of information transmission continues to increase, as well as the significant increase in the number of HPV-related news reports, it has never been more important to explore the role of media news coverage in the spread and control of the virus. Using a decreasing factor that captures the impact of media on the actions of people, this paper develops a model that characterizes the dynamics of HPV transmission with media impact, vaccination and recovery. We obtain global stability of equilibrium points employing geometric method, and further yield effective methods to contain the HPV pandemic by sensitivity analysis. With the center manifold theory, we show that there is a forward bifurcation when R0=1. Our study suggested that, besides controlling contact between infected and susceptible populations and improving effective vaccine coverage, a better intervention would be to strengthen media coverage. In addition, we demonstrated that contact rate and the effect of media coverage result in multiple epidemics of infection when certain conditions are met, implying that interventions need to be tailored to specific situations.

Dengue dynamics in Nepal: A Caputo fractional model with optimal control strategies

An infectious disease called dengue is a significant health concern nowadays. The dengue outbreak occurred with a single serotype all over Nepal in 2023. In the tropical and subtropical regions, dengue fever is a leading cause of sickness and death. Currently, there is no specified treatment for dengue fever. Avoiding mosquito bites is strongly advised to reduce the likelihood of controlling this disease. In underdeveloped countries like Nepal, the implementation of appropriate control measures is the most important factor in preventing and controlling the spread of dengue illness. The Caputo fractional dengue model with optimum control variables, including mosquito repellent and insecticide use, investigates the impact of alternative control strategies to minimize dengue prevalence. Using the fixed point theorem, the existence and uniqueness of a solution will be demonstrated for the problem. Ulam-Hyers stability, disease-free equilibrium point stability, and basic reproduction number are studied for the proposed model. The model is simulated using a two-step Lagrange interpolation technique, and the least squares method is used to estimate parameter values using real monthly infected data. We then analyze the sensitivity analysis to determine influencing parameters and the control measure effects on the basic reproduction number. The Pontryagin Maximum Principle is used to determine the optimal control variable in the dengue model for control strategies. The present study suggests that the deployment of control measures is extremely successful in lowering infectious disease incidences. Which facilitates the decision-makers to practice rigorous evaluation of such an epidemiological scenario while implementing appropriate control measures to prevent dengue disease transmission in Nepal.

An evolutionary game for analysis of implementation strategies in inter-regional diffusion of clean technology

Green transition driven by clean technology diffusion attaches importance to sustainable development. Therefore, choosing the optimal strategy in the inter-regional diffusion of clean technology has attracted the attention of local governments. In this regard, developed regions leading the R&D activities of clean technology can choose to fund or not fund diffusion, and developing regions can adopt absorption strategy or substitution strategy as well. An evolutionary game model was developed in this paper to analyze the stable strategies of local governments in the inter-regional diffusion of clean technology. In general, there are four possible stable equilibrium points. The result indicates that a collaboration relationship aiming at technology transfer among regions is likely to be established when developed regions and developing regions gain more from funding diffusion and subsidizing absorption, respectively. Especially, an increase of incentives from the central government, as well as an increase of cost performing independent clean R&D activities, improves the inter-regional diffusion network. This paper, therefore, provides a mathematical approach to examine the interaction of clean technology diffusion, and sheds light on actor behavior in games of inter-regional coordination.

Bifurcation and Chaos in a Fractional-Order Cournot Duopoly Game Model on Scale-Free Networks

In this study, a Cournot duopoly model describing Caputo fractional-order differential equations with piecewise constant arguments is discussed. We have obtained two-dimensional discrete dynamical system as a result of applying the discretization process to the model. By using the center manifold theory and the bifurcation theory, it is shown that the discrete dynamical system undergoes flip bifurcation about the Nash equilibrium point. Phase portraits, bifurcation diagrams, and Lyapunov exponents show the existence of many complex dynamical behaviors in the model such as the stable equilibrium point, period-2 orbit, period-4 orbit, period-8 orbit, period-16 orbit, and chaos according to changing the speed of the adjustment parameter [Formula: see text]. The discrete Cournot duopoly game model is also considered on two scale-free networks with different numbers of nodes. It is observed that the complex dynamical networks exhibit similar dynamical behaviors such as the stable equilibrium point, flip bifurcation, and chaos depending on changing the coupling strength parameter [Formula: see text]. Moreover, flip bifurcation and transition chaos take place earlier in more heterogeneous networks. Calculating the largest Lyapunov exponents guarantees the transition from nonchaotic to chaotic states in complex dynamical networks.

Complexity analysis of technology spillover supply chain with dual business objectives

ABSTRACT Based on the existence of technology spillover effect of manufacturers, considering the manufacturer's commercial goals of profit and market share, a dynamic price game model of technology spillover supply chain with dual commercial goals is established. Based on game theory and nonlinear dynamics theory, the equilibrium point and local stability of the system are analyzed. In addition, the effects of price adjustment rate, balance coefficient of business objectives, and technology spillover on model complexity, the change of manufacturer price and profit were analyzed through numerical simulation. The results show that the price adjustment rate, the balance coefficient of business objectives, and the technology spillover have an impact on the stability and complexity of the manufacturer's price competition model. Rationally controlling the price adjustment rate, and enhancing technological spillovers are beneficial to the healthy development of manufacturers and the healthy competition among manufacturers, maintaining the stability of product market prices, and effectively avoiding market chaos. Moreover, when considering market share, manufacturers can expand the system stability domain.

Study of chronic myeloid leukemia with T-cell under fractal-fractional order model

Abstract This research work is devoted to investigate myeloid leukemia mathematical model. We give some details about the existence of trivial and nontrivial equilibrium points and their stability. Also, local asymptotical stability of disease-free and endemic equilibrium points is discussed. Also, positivity of the solution has been discussed. Some sufficient results are achieved to study the local existence and uniqueness of solution to the considered model for Mittag–Leffler kernel using the Banach contraction theorem. Three numerical algorithms are derived in obtaining the numerical solution of suggested model under three different kernels using Adams–Basforth technique. Numerical results have been presented for different fractals and fractional orders to show the behavior of the proposed model.

Nonlinear dynamics and optimal control of a fractional order cotton leaf curl virus model incorporating climate change influences

This work investigates the transmission dynamics of a proposed fractional-order cotton leaf curl virus (CLCuV) model considering the effects of climate change. High temperatures resulting from climate change and global warming significantly impact the cotton crop and CLCuV transmission, causing substantial economic losses for farmers and agricultural industries. The study aims to address these impacts and contribute to the Sustainable Development Goals (SDGs) by promoting well-being and healthy lives. The characteristics of the CLCuV model’s solution are examined, including equilibrium points and the basic reproduction number. The stability of equilibrium points is analyzed in relation to key parameters, and sensitivity analysis of controllable parameters is performed. Additionally, optimal management strategies for controlling CLCuV spread are explored, aiming for minimal cost, improved health, and reduced climate change hazards. Numerical simulations are conducted to validate the theoretical findings.

Evolutionary Game Analysis of Governments’ and Enterprises’ Carbon-Emission Reduction

With the increasingly serious problem of global climate change, many countries are positively promoting carbon-emission-reduction actions. In order to deeply explore the interaction between enterprises’ carbon-emission reduction and governments’ regulation, this paper builds evolutionary game models between governments and enterprises under the reward-and-punishment mechanism. The peer-incentive mechanism is introduced to incentivize enterprises to reduce carbon emissions and coordinate governments and enterprises. The evolutionary-stability strategies are obtained by solving the evolutionary game models. The stability of equilibrium points under different situations is theoretically and numerically studied. The results show that the existence of peer incentives makes enterprises more inclined to positively reduce carbon emissions and governments more inclined to positively regulate. A sufficiently large peer fund can always encourage enterprises to choose positive carbon-reduction emission strategies, while governments choose positive regulation strategies. Not only the increasing rewards and fines but also lowering regulatory costs will promote carbon-emission-reduction behaviors of enterprises. Peer incentives are more effective in promoting positive emission reduction of enterprises compared with rewards and punishments. This study can provide important guidance for governments to formulate regulatory strategies and for enterprises to formulate emission-reduction strategies.

Investigating the Dynamic Behavior of Integer and Noninteger Order System of Predation with Holling’s Response

The paper’s primary objective is to examine the dynamic behavior of an integer and noninteger predator–prey system with a Holling type IV functional response in the Caputo sense. Our focus is on understanding how harvesting influences the stability, equilibria, bifurcations, and limit cycles within this system. We employ qualitative and quantitative analysis methods rooted in bifurcation theory, dynamical theory, and numerical simulation. We also delve into studying the boundedness of solutions and investigating the stability and existence of equilibrium points within the system. Leveraging Sotomayor’s theorem, we establish the presence of both the saddle-node and transcritical bifurcations. The analysis of the Hopf bifurcation is carried out using the normal form theorem. The model under consideration is extended to the fractional reaction–diffusion model which captures non-local and long-range effects more accurately than integer-order derivatives. This makes fractional reaction–diffusion systems suitable for modeling phenomena with anomalous diffusion or memory effects, improving the fidelity of simulations in turn. An adaptable numerical technique for solving this class of differential equations is also suggested. Through simulation results, we observe that one of the Lyapunov exponents has a negative value, indicating the potential for the emergence of a stable-limit cycle via bifurcation as well as chaotic and complex spatiotemporal distributions. We supplement our analytical investigations with numerical simulations to provide a comprehensive understanding of the system’s behavior. It was discovered that both the prey and predator populations will continue to coexist and be permanent, regardless of the choice of fractional parameter.

Wiener and Lévy processes to prevent disease outbreaks: Predictable vs stochastic analysis

The study considers stochastic behavior along with modeling, mathematical analysis, theory development, and numerical simulation of the COVID-19 virus. To evaluate current trends and make future projections regarding the basic reproduction number and infection case, we have taken the modified five-compartment SEIRD mathematical model. Since the basic reproduction number R0 is not sufficient to predict the outbreak, we applied the Itô stochastic differential equations (SDEs) along with the Weiner process and Lévy jump to investigate the nature of the disease outbreak. Since COVID-19 is an infectious disease caused by the severe acute respiratory syndrome coronavirus 2 “(SARS-CoV-2)”, that has common symptoms including fever, cough, and difficulty breathing, we have illustrated the boundedness and positivity of solutions and investigated the stability of endemic and disease-free equilibrium points. The findings show that the dynamics of disease epidemics are influenced by contact patterns. For all the models, the threshold quantity, R0 is calculated which is the key reason to prove the global and local stability analysis of disease-free equilibrium and endemic equilibrium points. Using the least squares method for data fitting, we have performed a case study on COVID-19 data in Italy for this article. A sensitivity analysis is part of our investigation to find important factors. This paper investigates a stochastic epidemic model (SDE) using Lévy jumps and Weiner processes. Local and global stability in a stochastic setting are examined in the analysis. To view the analytical study, a multitude of numerical results are achieved. The study is important to comprehend the stochastic behavior of the virus in the contagious model in order to prevent similar disease outbreaks in the future.

Data-driven economic MPC with asymptotic stability and strong duality verification using Hankel matrix

We consider the problem of dynamic regulation with an economic cost function to control unknown linear systems, in which improving the economic performance and guaranteeing the stability of economical optimal equilibrium point are control objectives. A data-driven economic MPC scheme is presented using measured input-output trajectories without a prior system identification step. Our method uses Hankel matrices which include one input-output data trajectory for prediction in economic MPC, while persistently exciting of the input generating the data is needed. One of the novelties of the presented framework is directly verifying the strong duality property from input-output trajectory with the general cost function, considered as the supply rate. This is used to find a Lyapunov function for data-driven economic MPC. Under the strong duality assumption, asymptotic stability of the economical optimal equilibrium point for the closed-loop system with terminal equality constraint is guaranteed. The proposed data-driven economic MPC approach needs only persistently exciting data trajectory along with an upper bound on the system order and need no model description and no online parameter estimation. The proposed scheme applicability compared to the existing model-based economic MPC and data-driven MPC is illustrated for continuous stirred tank reactor (CSTR) and a numerical example and the robustness of the proposed scheme is evaluated in the case of measurement noise, as well as nonlinear model for CSTR system.

Energy flow investigations of Rayleigh-Plesset equation for cavitation simulations

Nonlinear Rayleigh-Plesset Equation (RPE) for cavitation simulations is investigated using energy-flow theory. Nondimensional RPE affected by Reynolds-number, surface-tension, bubble-pressure, and liquid static-dynamic pressure ratio is derived to examine its equilibrium points with stability, possible periodical/chaotic motions by the energy-flow criteria. An example is numerically analyzed to illustrate the developed method. Five cases of the example reveal: 1) It is a damped system, where initial disturbances are gradually reduced with its phase point tends to its stable equilibrium-point; 2) Reynolds-number affects the damping of system, large one corresponds small damping; 3) Bubble-pressure, surface-tension and liquid static-dynamic pressure ratio affect the position of equilibrium point; 4) Periodical orbit appears only in forced vibrations, in which free vibration is reduced with time forward and the system finally shows a stable periodical oscillation; 5) Energy flow criteria for chaotic motions is not reached, and there are no chaotic motions for the cases of example. Numerical simulations confirm the developed energy-flow means with available computer code is effective to investigate generalized RPEs to reveal their inherent characteristics affected by required parameters in engineering cavitation analysis and designs, such as considering mass transports across the boundary of bubble by evaporation or condensation of liquids.

Estimating Markov Chain Mixing Times: Convergence Rate Towards Equilibrium of a Stochastic Process Traffic Assignment Model

Network equilibrium models have been extensively used for decades. The rationale for using equilibrium as a predictor is essentially that (i) a unique and globally stable equilibrium point is guaranteed to exist and (ii) the transient period over which a system adapts to a change is sufficiently short in time that it can be neglected. However, we find transport problems without a unique and stable equilibrium in the literature. Even if it exists, it is not certain how long it takes for the system to reach an equilibrium point after an external shock onto the transport system, such as infrastructure improvement and damage by a disaster. The day-to-day adjustment process must be analysed to answer these questions. Among several models, the Markov chain approach has been claimed to be the most general and flexible. It is also advantageous as a unique stationary distribution is guaranteed in mild conditions, even when a unique and stable equilibrium does not exist. In the present paper, we first aim to develop a methodology for estimating the Markov chain mixing time (MCMT), a worst-case assessment of the convergence time of a Markov chain to its stationary distribution. The main tools are coupling and aggregation, which enable us to analyse MCMTs in large-scale transport systems. Our second aim is to conduct a preliminary examination of the relationships between MCMTs and critical properties of the system, such as travellers’ sensitivity to differences in travel cost and the frequency of travellers’ revisions of their choices. Through analytical and numerical analyses, we found key relationships in a few transport problems, including those without a unique and stable equilibrium. We also showed that the proposed method, combined with coupling and aggregation, can be applied to larger transport models. History: This paper has been accepted for the Transportation Science Special Issue on the 25th International Symposium on Transportation and Traffic Theory. Funding: This study was financially supported by the Japan Society for the Promotion of Science [Grant-in-Aid 20H00265]. Supplemental Material: The online appendices are available at https://doi.org/10.1287/trsc.2024.0523 .

Dynamical behavior of fractal-fractional order monkeypox virus model

A virus that infects both humans and animals, monkeypox is common in many West African countries and has sporadically spread to other parts of the world. There are serious public health concerns around the globe as a result of the recent spike in cases of monkeypox among both endemic and non-endemic populations. This paper proposes to use the Caputo sense of fractal fractional-order derivatives to examine the dynamics of monkeypox transmission. The study uses Schauder’s fixed point theorem to evaluate the solutions qualitatively and establish their uniqueness inside the model. The next generation matrix technique is used to compute the essential reproduction number. The stability of equilibrium points is further investigated in the research, and a sensitivity analysis of model parameters is carried out. When the basic reproduction number R0 is less than 1, the equilibrium without infections is locally stable. Also, when R0 exceeds 1, this equilibrium becomes unstable. The proposed model incorporates Ulam–Hyers stability through nonlinear functional analysis. To estimate solutions for the fractal-fractional order monkeypox model, the Lagrange’s interpolation method is utilized. In addition, numerical simulations are carried out to examine the influence of some parameters on the overall dynamics of the model. Numerical simulations are performed using MATLAB software to exemplify the model behavior in the context of the Nigerian case study. The graphical representations suggest that the fractal fractional order affects the dynamics of the monkeypox. The findings indicate that isolation of infected individuals in the human population helps to reduce disease transmission.