{"title":"Hopf bifurcation in a generalized Goodwin model with delay","authors":"Eysan Sans , Melisa Akdemir , Ayse Tiryakioglu , Ayse Peker-Dobie , Cihangir Ozemir","doi":"10.1016/j.matcom.2025.04.012","DOIUrl":"10.1016/j.matcom.2025.04.012","url":null,"abstract":"<div><div>Goodwin’s model is a cornerstone in the study of dynamical systems within macroeconomics, explaining the interaction between employment ratio and wage share in a closed economy. Analogous to predator–prey dynamics in mathematical economics, the Goodwin model, despite its simplicity, effectively captures the periodic behavior of state variables over specific time intervals. By relaxing the initial assumptions, the model can be adapted to account for more complex economic scenarios. In this article, we study a higher-dimensional extension of the Goodwin model that incorporates variable capacity utilization and capital coefficient alongside employment ratio and wage share. In particular instances, the wage share and employment rate equations decouple from the overall system. For these cases, by incorporating a delay effect in the Phillips curve, we demonstrate that while the equilibrium of the generalized system remains stable within certain parameter domains in the absence of delay, the introduction of delay can induce a Hopf bifurcation, leading to periodic oscillations. We analytically derive the critical delay parameter value that destabilizes the equilibrium point via a Hopf bifurcation.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"237 ","pages":"Pages 86-106"},"PeriodicalIF":4.4,"publicationDate":"2025-04-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143906032","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"N-policy for redundant machining system with double retrial orbits using soft computing techniques","authors":"Vijay Pratap Singh , Madhu Jain , Richa Sharma","doi":"10.1016/j.matcom.2025.04.025","DOIUrl":"10.1016/j.matcom.2025.04.025","url":null,"abstract":"<div><div>The present study is concerned with the performance prediction of a double retrial orbit redundant repairable machining system. Both primary and secondary orbits are available as waiting/buffer space for the failed units. In these orbits, the failed units can reside and make re-attempts for the repair. As per N-policy, if there are no units in the orbits for the repairing job, the repairman goes on vacation and further starts the repair job when N-failed units are accumulated. The objective of this investigation is to evaluate the transient and steady-state distributions of the queue length of failed units under N-policy. The matrix analytic and matrix recursive methods are utilized for solution purpose while an adaptive neuro-fuzzy inference system (ANFIS) is employed for validating the feasibility of designing the AI-based controller. The harmonic search (HS) and particle swarm optimization (PSO) methods have been implemented for the cost optimization purpose so as to evaluate the optimal design parameters. The outputs of study provides critical insights into optimal system performance and improving the repair policy. Furthermore, a practical application of this investigation is demonstrated in a telecommunications network traffic system, where the proposed methods can be utilized to manage the maintenance issues of routers in the network traffic.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"237 ","pages":"Pages 42-69"},"PeriodicalIF":4.4,"publicationDate":"2025-04-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143906036","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Composite likelihood expectation–maximization algorithm for the first-order multivariate integer-valued autoregressive model with multivariate mixture distributions","authors":"Weiyang Yu, Haitao Zheng","doi":"10.1016/j.matcom.2025.04.017","DOIUrl":"10.1016/j.matcom.2025.04.017","url":null,"abstract":"<div><div>In this paper, we develops a estimation framework for first-order multivariate integer-valued autoregressive (MINAR(1)) models with multivariate Poisson-lognormal (MPL) or multivariate geometric-logitnormal (MGL) innovations. Owing to the structural complexity of the MINAR(1) framework and the latent dependence in MPL/MGL innovation processes, the likelihood functions involve summations and integrals. Traditional expectation–maximization (EM) algorithms face prohibitive computational demands as model dimensionality increases. To address this, we propose a composite likelihood expectation–maximization (CLEM) algorithm that strategically combines composite likelihood with the EM algorithm, reducing the computational burden. Furthermore, we implement a Cholesky parameterization for the covariance matrix to ensure positive definiteness. Through comprehensive Monte Carlo simulations, we demonstrate the performance of CLEM algorithm. Finally, we validate the method’s practical utility by analyzing a real-world dataset.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"237 ","pages":"Pages 167-187"},"PeriodicalIF":4.4,"publicationDate":"2025-04-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143912495","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Local convergence analysis of L1/finite element scheme for a constant delay reaction-subdiffusion equation with uniform time mesh","authors":"Weiping Bu, Xin Zheng","doi":"10.1016/j.matcom.2025.04.014","DOIUrl":"10.1016/j.matcom.2025.04.014","url":null,"abstract":"<div><div>The aim of this paper is to develop a refined error estimate of L1/finite element scheme for a reaction-subdiffusion equation with constant delay <span><math><mi>τ</mi></math></span> based on the uniform time mesh. Under the non-uniform multi-singularity assumption of exact solution in time, the local truncation error of L1 scheme is investigated. Then we introduce a fully discrete finite element scheme of the considered problem. In order to investigate the error estimate, a novel discrete fractional Grönwall inequality with delay term is proposed, which does not include the increasing Mittag-Leffler function. By applying this Grönwall inequality, we obtain a local error estimate of the above fully discrete scheme without including the Mittag-Leffler function. In particular, the convergence result implies that, for the considered time interval <span><math><mrow><mo>(</mo><mrow><mo>(</mo><mi>i</mi><mo>−</mo><mn>1</mn><mo>)</mo></mrow><mi>τ</mi><mo>,</mo><mi>i</mi><mi>τ</mi><mo>]</mo></mrow></math></span>, although the convergence rate in the sense of maximum time error is low for the first interval, i.e. <span><math><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow></math></span>, it will be improved as the increasing <span><math><mi>i</mi></math></span>, which is consistent with the factual assumption that the smoothness of the solution will be improved as the increasing <span><math><mi>i</mi></math></span>. Finally, we present some numerical tests to verify the developed theory.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"237 ","pages":"Pages 70-85"},"PeriodicalIF":4.4,"publicationDate":"2025-04-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143906037","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Lanndon Ocampo , Dharyll Prince Abellana , Julius Selle , Rhoda Namoco , Fatima Maturan , John Kevin Padro , Samantha Shane Evangelista
{"title":"Dynamical system modelling for picture fuzzy multi-criteria public opinion on transportation infrastructure","authors":"Lanndon Ocampo , Dharyll Prince Abellana , Julius Selle , Rhoda Namoco , Fatima Maturan , John Kevin Padro , Samantha Shane Evangelista","doi":"10.1016/j.matcom.2025.04.024","DOIUrl":"10.1016/j.matcom.2025.04.024","url":null,"abstract":"<div><div>This work presents a novel approach for modelling dynamic fuzzy multi-criteria decision-making problems. We consider a problem where members of a population evaluate alternatives under a given criterion by taking either of three options: agree, neutral, or disagree. The physical partitioning of the population gives rise to the notion of describing the performance of alternatives under multiple criteria as Picture fuzzy sets, with operations that can aggregate multi-criteria performance into real numbers. Then, the dynamic performance of the alternatives is described by considering the interactions of the partitions of the populations, modelled as a system of differential equations constructed for each criterion. We provided proof of the basic properties of each system, including nonnegativity and boundedness of solutions, and presented the points of convergence. Parameter values were estimated by performing time-series clustering. The proposed approach was applied numerically in modelling multi-criteria public opinion of a bus rapid transit project—a single alternative, where nine criteria were evaluated. Initial performance values under multiple criteria were obtained using a survey, and four scenarios representing possible government interventions were evaluated. The time-dependent performance profiles of the public about the project were obtained, with insights indicating the efficacy of the proposed model.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"236 ","pages":"Pages 455-474"},"PeriodicalIF":4.4,"publicationDate":"2025-04-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143899208","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Hierarchical Newton iterative identification methods for a class of input multi-piecewise Hammerstein models with autoregressive noise","authors":"Yamin Fan , Ximei Liu , Meihang Li","doi":"10.1016/j.matcom.2025.04.019","DOIUrl":"10.1016/j.matcom.2025.04.019","url":null,"abstract":"<div><div>Advancements in mathematical theories and applied technologies have driven research on modeling and identification of complex nonlinear systems, yet existing models still face challenges in terms of model structures, accuracy of parameter estimation, and efficiency. In this study, we present a class of input multi-piecewise Hammerstein models that utilize piecewise-linear function to capture arbitrary nonlinear characteristics. By employing the key term separation technique and the Newton iterative algorithm for parameter estimation, a generalized Newton iterative algorithm is proposed, which overcomes the limitations of conventional identification methods in handling complex nonlinearities. Additionally, considering the computational load caused by the numerous parameters in multi-piecewise linear function and the inversion of the Hessian matrix in the Newton iterative algorithm, the hierarchical identification principle is introduced and a three-stage generalized Newton iterative algorithm is derived for enhancing the computational efficiency. The feasibility of the presented methods are demonstrated through a simulation example.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"237 ","pages":"Pages 247-262"},"PeriodicalIF":4.4,"publicationDate":"2025-04-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143918095","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Linfeng Zhang , Guorong Wang , Lin Zhong , Qingping Li , Qiang Fu , Mengdi Wu , Lei Chen
{"title":"Nonlocal modeling for gas-fluid seepage coupled with heat transfer in gas hydrate by using peridynamic differential operator","authors":"Linfeng Zhang , Guorong Wang , Lin Zhong , Qingping Li , Qiang Fu , Mengdi Wu , Lei Chen","doi":"10.1016/j.matcom.2025.04.021","DOIUrl":"10.1016/j.matcom.2025.04.021","url":null,"abstract":"<div><div>A coupled non-local thermal-hydro numerical based on the peridynamic differential operator (PDDO) is presented for simulating gas-fluid flow coupled with heat transfer of natural gas hydrate (NGH). The non-local integral replaces the governing equations for nonlinear seepage motion heat conduction. The nonlocal governing equations are resolved by employing an explicit iterative technique. This method is applied to model heat transfer in saturated NGH porous media with accompanying gas-liquid seepage. The effects of boundary pressure and temperature changes on seepage behavior are also investigated. The efficacy and precision of the established coupled model were validated by comparing the temperature field and pressure field outcomes derived from Lagrange multiplier variation analysis with the finite element findings and analytical solutions documented in previous literature. This paper introduces a different approach to numerically simulate the flow of gas and fluid, while also considering heat transfer issues.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"237 ","pages":"Pages 188-212"},"PeriodicalIF":4.4,"publicationDate":"2025-04-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143912496","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Numerical steepest descent method for computing oscillatory-type Bessel integral transforms","authors":"Ruyun Chen, Yu Li, Yongxiong Zhou","doi":"10.1016/j.matcom.2025.04.016","DOIUrl":"10.1016/j.matcom.2025.04.016","url":null,"abstract":"<div><div>In this paper, numerical steepest descent method is implemented to approximate highly oscillatory Bessel-type integral transforms. We begin our analysis by utilizing an important relationship between Bessel function of the first kind and modified Bessel function of the second kind. Subsequently, we transform new integrals into the forms on the interval <span><math><mrow><mo>[</mo><mn>0</mn><mo>,</mo><mo>+</mo><mi>∞</mi><mo>)</mo></mrow></math></span>, where the integrands do not oscillate and decay exponentially fast. These integrals can then be efficiently computed using Gauss–Laguerre quadrature rule. Furthermore, we derive the theoretical error estimates that depend on the frequency <span><math><mi>ω</mi></math></span> and the number of nodes <span><math><mi>n</mi></math></span>. Numerical examples based on the theoretical results are provided to demonstrate the effectiveness of these methods.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"236 ","pages":"Pages 320-333"},"PeriodicalIF":4.4,"publicationDate":"2025-04-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143874477","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Numerical algorithms for the reconstruction of space-dependent sources in thermoelasticity","authors":"Frederick Maes , Karel Van Bockstal","doi":"10.1016/j.matcom.2025.04.004","DOIUrl":"10.1016/j.matcom.2025.04.004","url":null,"abstract":"<div><div>This paper investigates the inverse problems of determining a space-dependent source for thermoelastic systems of type-III under adequate time-averaged or final-in-time measurements and conditions on the time-dependent part of the sought source. Several numerical methods are proposed and examined, including a Landweber scheme and minimisation methods for the corresponding cost functionals, which are based on the gradient and conjugate gradient method. A shortcoming of these methods is that the values of the sought source are fixed ab initio and remain fixed during the iterations. The Sobolev gradient method is applied to overcome the possible inaccessibility of the source values at the boundary. Numerical examples are presented to discuss the different approaches and support our findings based on the implementation on the FEniCSx platform.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"236 ","pages":"Pages 426-454"},"PeriodicalIF":4.4,"publicationDate":"2025-04-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143882296","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Convergence of individual opinions during a bi-virus epidemic and its impact","authors":"Shidong Zhai , Zhenmei Zhang , Jun Ma","doi":"10.1016/j.matcom.2025.04.006","DOIUrl":"10.1016/j.matcom.2025.04.006","url":null,"abstract":"<div><div>This paper investigates the intricate interplay between individual opinions and bi-virus transmission in a city or country. The proposed model considers two competing viruses, where virus transmission is shaped by individuals’ perceptions of the epidemic, which are, in turn, influenced by the transmission dynamics. Building on the concept of an opinion-dependent reproduction number, we conduct a rigorous theoretical analysis of the proposed model, focusing on equilibrium stability and parameter sensitivity. From the analysis of simulated data, we derive the following insights: (1) Negative opinions during an epidemic, such as believing the epidemic is not severe or distrusting public health measures, may result in inaccurate community responses and potentially reignite outbreaks. (2) In the bi-virus-opinion model, while perceived epidemic severity can reduce infection rates, it does not alter the dominance of the prevailing virus under conditions of high epidemic severity. Furthermore, we examine the role of stubborn opinions to reshape public perception, which can facilitate epidemic eradication under mild epidemic conditions. In the case of dual-virus competition, stubborn opinions can alter the final competitive outcomes. Finally, numerical simulations are performed to validate the findings and corroborate them with real-world data from the SARS-CoV-2 pandemic.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"236 ","pages":"Pages 334-353"},"PeriodicalIF":4.4,"publicationDate":"2025-04-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143874387","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}