Mathematics and Computers in Simulation最新文献

{"title":"An innovative methodology in scrutinizing the nonlinear instability of two immiscible MHD viscoelastic liquids","authors":"Galal M. Moatimid,&nbsp;Yasmeen M. Mohamed","doi":"10.1016/j.matcom.2025.03.017","DOIUrl":"10.1016/j.matcom.2025.03.017","url":null,"abstract":"<div><div>This study examines the nonlinear stability of two distinct viscoelastic magneto-rheological planner fluids that are immersed in porous media. The lower zone is filled with the Reiner-Rivlin fluid (RRF); meanwhile, the upper one is occupied by the Eyring-Powell fluid (EPF). An unchanged magnetic field (MF) is applied to the whole structure, and the effects of surface tension (ST) and mass and heat transfer (MHT) are also documented. This approach offers insights into stability thresholds and flow dynamics crucial for applications in energy systems, medicinal devices, and industrial processes that involve multi-layered magneto-rheological fluids in porous settings. One use of these fluids is the real-time adjustment of damping qualities in adaptive vehicle suspension systems. Improvements in vehicle dynamic performance, comfort, and safety are directly impacted by this research. The calculations are shortened by making use of viscous potential theory (VPT). Therefore, the viscoelastic influences are considered in order to show how the nonlinear boundary conditions (BCs) produce their contributions. Consequently, the impacts of the viscoelasticity parameters are removed from the solution of the equations of motion. Merging the fundamental linear hydrodynamic equations with Maxwell's equations over the quasi-static approximations, the boundary-value problem is raised. A popular nonlinear ordinary differential equation (ODE) can be transformed into a linear via the He’s frequency formula (HFF), which forms the basis of the non-perturbative approach (NPA). The non-dimensional analysis reveals a set of physical dimensionless numerals. Additionally, they help to reduce the amount of variables that are needed to comprehend the framework. The stability constraints are numerically tested in the complex scenario, and the stability mechanism remains consistent for both real and imaginary coefficients within the nonlinear characteristic equation governing interface displacement. Polar plots of unstable solutions are omitted, since these solutions are not desired.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"234 ","pages":"Pages 472-496"},"PeriodicalIF":4.4,"publicationDate":"2025-03-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143726019","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":"Three-dimensional D3Q27 multiple-relaxation-time lattice Boltzmann simulation of Herschel–Bulkley viscoelastic fluids in a cubic cavity with top lid driven diagonally","authors":"Md. Mamun Molla,&nbsp;Md. Mahadul Islam","doi":"10.1016/j.matcom.2025.03.015","DOIUrl":"10.1016/j.matcom.2025.03.015","url":null,"abstract":"<div><div>Graphics Processing Unit (GPU) accelerated multiple-relaxation-time (MRT) lattice Boltzmann method (LBM) is used for the simulation of Herschel–Bulkley non-Newtonian fluids in a three-dimensional (3D) cubic cavity with the top lid-driven diagonally. For the 3D simulation, a D3Q27 lattices model, which is more stable and well-accepted in the LBM community, is used in the present MRT-LBM. Simulations using numerical models are run for a variety of dimensionless variables, including the Reynolds numbers (<span><math><mrow><mi>R</mi><mi>e</mi><mo>=</mo><mn>300</mn><mo>,</mo><mn>600</mn><mo>,</mo><mn>1000</mn><mo>,</mo><mn>1200</mn></mrow></math></span>), Bingham number (<span><math><mrow><mi>B</mi><mi>n</mi><mo>=</mo><mn>0</mn><mo>.</mo><mn>0</mn><mo>,</mo><mn>0</mn><mo>.</mo><mn>5</mn><mo>,</mo><mn>1</mn><mo>.</mo><mn>0</mn><mo>,</mo><mn>2</mn><mo>.</mo><mn>0</mn></mrow></math></span>), Power-law index, (<span><math><mrow><mi>n</mi><mo>=</mo><mn>0</mn><mo>.</mo><mn>8</mn></mrow></math></span>). In the present numerical simulation, the GPU has used a parallel computing technique based on the Compute Unified Device Architecture (CUDA) C++ programming. MRT-LBM code is validated for the Newtonian and non-Newtonian power law fluid with a lid-driven cubic cavity. The numerical results obtained regarding the streamlines, velocity, viscosity distributions, and the iso-surfaces of the non-Newtonian viscosity are presented. The current numerical findings could potentially function as benchmark results for validating 3D codes validation for the non-Newtonian fluids.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"234 ","pages":"Pages 419-437"},"PeriodicalIF":4.4,"publicationDate":"2025-03-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143685635","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":"Dynamical behavior and bifurcation analysis for a theoretical model of dengue fever transmission with incubation period and delayed recovery","authors":"Burcu Gürbüz ,&nbsp;Aytül Gökçe ,&nbsp;Segun I. Oke ,&nbsp;Michael O. Adeniyi ,&nbsp;Mayowa M. Ojo","doi":"10.1016/j.matcom.2025.03.008","DOIUrl":"10.1016/j.matcom.2025.03.008","url":null,"abstract":"<div><div>As offered by the World Health Organisation (WHO), close to half of the population in the world’s resides in dengue-risk zones. Dengue viruses are transmitted to individuals by Aedes mosquito species infected bite (Ae. Albopictus of Ae. aegypti). These mosquitoes can transmit other viruses, including Zika and Chikungunya. In this research, a mathematical model is formulated to reflect different time delays considered in both extrinsic and intrinsic incubation periods, as well as in the recovery periods of infectious individuals. Preliminary results for the non-delayed model including positivity and boundedness of solutions, non-dimensionalization and equalibria analysis are presented. The threshold parameter (reproduction number) of the model is obtained via next generation matrix schemes. The stability analysis of the model revealed that various dynamical behavior can be observed depending on delay parameters, where in particular the effect of delay in the recovery time of infectious individuals may lead to substantial changes in the dynamics. The ideas presented in this paper can be applied to other infectious diseases, providing qualitative evaluations for understanding time delays influencing the transmission dynamics.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"234 ","pages":"Pages 497-513"},"PeriodicalIF":4.4,"publicationDate":"2025-03-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143739609","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Global dynamics of a chemotaxis system with toxicity in invasive species","authors":"Xiaoyue Yuan ,&nbsp;Xuebing Zhang ,&nbsp;Wenjun Liu ,&nbsp;Ali Moussaoui ,&nbsp;Pierre Auger","doi":"10.1016/j.matcom.2025.03.009","DOIUrl":"10.1016/j.matcom.2025.03.009","url":null,"abstract":"<div><div>Invasive species threaten the integrity of ecosystems by altering the structure and function of natural systems. It is crucial to predict the mode of biological invasion and control intruders. In this paper, we establish a diffusion biological invasion model with toxicant-taxis and conduct research through theoretical analysis and numerical simulation methods. We investigate the local stability of the system and find that it undergoes saddle–node bifurcation and transcritical bifurcation. We further prove the boundedness and global existence of the classical solutions of the system. By constructing appropriate Lyapunov functionals, the global stability of the positive steady state is analyzed, and the decay rate of the solution is provided. In addition, to investigate the effects of toxins and competition intensity on the survival of invasive species, we demonstrate the existence conditions of steady-state solutions through numerical simulations. By comparison, it is found that invasive species can only survive in new environments if they possess at least one advantage.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"235 ","pages":"Pages 16-36"},"PeriodicalIF":4.4,"publicationDate":"2025-03-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143734913","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":"Nonsingularity of unsymmetric Kansa matrices: Random collocation by MultiQuadrics and Inverse MultiQuadrics","authors":"R. Cavoretto ,&nbsp;A. De Rossi ,&nbsp;F. Dell’Accio ,&nbsp;A. Sommariva ,&nbsp;M. Vianello","doi":"10.1016/j.matcom.2025.03.005","DOIUrl":"10.1016/j.matcom.2025.03.005","url":null,"abstract":"<div><div>Unisolvence of unsymmetric Kansa collocation is still a substantially open problem. We prove that Kansa matrices with MultiQuadrics and Inverse MultiQuadrics for the Dirichlet problem of the Poisson equation are almost surely nonsingular, when the collocation points are chosen by any continuous random distribution in the domain interior and arbitrarily on its boundary.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"234 ","pages":"Pages 390-395"},"PeriodicalIF":4.4,"publicationDate":"2025-03-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143685633","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A family of C1 Clough–Tocher spline spaces on C0 piecewise quadratic domain partitions","authors":"Jan Grošelj,&nbsp;Marjeta Knez","doi":"10.1016/j.matcom.2025.03.006","DOIUrl":"10.1016/j.matcom.2025.03.006","url":null,"abstract":"<div><div>The paper addresses the construction of <span><math><msup><mrow><mi>C</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span> splines on a curved domain that is parametrized by a <span><math><msup><mrow><mi>C</mi></mrow><mrow><mn>0</mn></mrow></msup></math></span> piecewise geometry mapping composed of quadratic Bézier triangles. The <span><math><msup><mrow><mi>C</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span> splines are assembled from polynomials of a chosen total degree greater than or equal to four, and their construction is based on the Clough–Tocher splitting technique that ensures locality. In particular, the splines are locally characterized by an interpolation problem described by Hermite data, which resembles the standard macro-element concepts developed for <span><math><msup><mrow><mi>C</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span> splines on triangulations.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"234 ","pages":"Pages 368-389"},"PeriodicalIF":4.4,"publicationDate":"2025-03-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143685114","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Stability and bifurcation analysis of a time-order fractional model for water-plants: Implications for vegetation pattern formation","authors":"Shanwei Li,&nbsp;Yimamu Maimaiti","doi":"10.1016/j.matcom.2025.03.007","DOIUrl":"10.1016/j.matcom.2025.03.007","url":null,"abstract":"<div><div>The water-plant model is a significant tool for studying vegetation patterns. It helps researchers understand the interactions between water availability and plant growth, which are crucial for analyzing ecological dynamics and predicting changes in vegetation due to environmental factors. However, there has been limited research on the memory effect associated with the water-plant model. This paper investigates a fractional-order water-plant model with cross-diffusion, in which the fractional order signifies the memory effect. First, we examine the conditions for the equilibrium point in a spatially homogeneous model, followed by an analysis of the model’s linear stability and the existence of Hopf bifurcation. Subsequently, we analyze the stability of spatiotemporal models incorporating cross-diffusion, along with the presence of Turing bifurcation, Hopf bifurcation, and Turing–Hopf bifurcation. Finally, we present several numerical simulations to validate the theoretical results. The results indicate that the Hopf bifurcation parameters increase with the fractional order <span><math><mi>τ</mi></math></span>, leading to a larger parameter space for Hopf instability. As the fractional order <span><math><mi>τ</mi></math></span> increases, it results in a smaller parameter space for Turing instability and a reduced parameter space for stability. This indicates that an increase in the fractional order <span><math><mi>τ</mi></math></span> accelerates the transition of vegetation patterns, thereby affecting the stability of the system.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"234 ","pages":"Pages 342-358"},"PeriodicalIF":4.4,"publicationDate":"2025-03-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143654553","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":"An asymptotic approximation of the solution for nearly tridiagonal quasi-Toeplitz linear systems","authors":"Philsu Kim ,&nbsp;Sangbeom Park ,&nbsp;Seonghak Kim,&nbsp;Soyoon Bak","doi":"10.1016/j.matcom.2025.02.024","DOIUrl":"10.1016/j.matcom.2025.02.024","url":null,"abstract":"<div><div>We introduce an asymptotic approximate algorithm for solving nearly tridiagonal quasi-Toeplitz linear systems. When addressing low-rank perturbations of a tridiagonal Toeplitz matrix system based on the Sherman–Morrison–Woodbury formula (or Woodbury identity), conventional methods require solving at least two simpler systems. The proposed algorithm overcomes this limitation by providing an explicit asymptotic formula for one of these systems. This asymptotic approximation enables a rapid resolution of the original system with minimal additional computation. To validate the accuracy and efficiency of the proposed algorithm, we conduct numerical experiments on two cases, comparing the results with those of existing methods. The results demonstrate that the proposed algorithm significantly reduces computation time while maintaining accuracy compared to the existing methods.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"234 ","pages":"Pages 359-367"},"PeriodicalIF":4.4,"publicationDate":"2025-03-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143685634","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":"Dynamical analysis and optimal control strategy of seasonal brucellosis","authors":"Huidi Chu ,&nbsp;Xinmiao Rong ,&nbsp;Liu Yang ,&nbsp;Meng Fan","doi":"10.1016/j.matcom.2025.03.003","DOIUrl":"10.1016/j.matcom.2025.03.003","url":null,"abstract":"<div><div>Brucellosis exhibits typical seasonal patterns and shows a notable rising trend in recent years, posing a serious threat to public health and economic development. Experimental research indicates that increased tick activity may elevate brucellosis transmission risk although the quantitative impact of ticks remains insufficiently explored. To investigate the seasonal transmission mechanisms of Brucella, identify the key factors, and assess ticks’ potential role, a multi-population non-autonomous periodic dynamical model is developed. The global dynamics of the model such as extinction, uniform persistence, disease-free periodic solution, and endemic periodic solution are well explored in terms of the basic reproduction number. Theoretical and numerical analyses demonstrate that, while tick control helps mitigate transmission risks, it is insufficient to eliminate periodic transmission. Effective control of brucellosis requires a comprehensive approach, especially culling infected sheep and improving vaccination coverage to curb the overall rising trend. Additionally, adjusting sheep reproductive schedules within the sheep’s life cycle, such as delaying the peak time of birth and advancing the peak time of abortion, is crucial for managing seasonal transmission. Numerical simulations of the optimal control strategies reveal that adjusting interventions based on seasonal fluctuations in infections balances the cost and effectiveness while highlighting the importance of effective tick control.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"234 ","pages":"Pages 299-324"},"PeriodicalIF":4.4,"publicationDate":"2025-03-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143628609","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":"Estimating the region of attraction on fractional-order complex networks with time-varying delay","authors":"Feifei Du ,&nbsp;Jun-Guo Lu ,&nbsp;Qing-Hao Zhang","doi":"10.1016/j.matcom.2025.02.030","DOIUrl":"10.1016/j.matcom.2025.02.030","url":null,"abstract":"<div><div>Recently, there has been increasing attention towards reckoning the region of attraction (ROA) for networks. However, applying existing theory to networks with fractional-order and delays presents significant challenges. This article addresses the estimation of ROA for fractional-order complex networks with time-varying delay. Initially, two generalized fractional-order Halanay inequalities are formulated. Subsequently, leveraging the first Halanay inequality, a method for ROA estimation is developed, which is unaffected by both delay and fractional-order. However, this method tends to be conservative. To mitigate this conservatism, a delay-dependent and order-dependent ROA estimation technique is proposed based on our two developed Halanay inequalities. Additionally, numerical examples are presented to validate the proposed methodologies.</div></div>","PeriodicalId":49856,"journal":{"name":"Mathematics and Computers in Simulation","volume":"234 ","pages":"Pages 438-458"},"PeriodicalIF":4.4,"publicationDate":"2025-03-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143697557","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}