Applied Mathematical Modelling最新文献

{"title":"A hydrodynamic model for simulating particles with generalised cross-sectional shapes in two-phase flow systems","authors":"Wenhui Zhang ,&nbsp;Yanli Chen ,&nbsp;Xiaodong Wang ,&nbsp;Shuke Li","doi":"10.1016/j.apm.2025.116132","DOIUrl":"10.1016/j.apm.2025.116132","url":null,"abstract":"<div><div>Modeling the motion and interactions of particles with complex shapes in two-phase flow systems presents significant challenges. In this study, we aim to develop an efficient and stable hydrodynamic model for simulating two-phase flow systems involving particles with generalised cross-sectional shapes. Firstly, a geometric representation model (GRM) is developed to precisely characterize arbitrary particle geometries while maintaining uniform Lagrangian points spacing. Next, to overcome traditional limitations, a boundary thickening-based direct forcing Immersed Boundary Method (BTDF-based IBM) and a Characteristic-Based Split (CBS) method have been integrated within the finite element framework for the first time. Unlike iterative IBM, BTDF-based IBM enforces no-slip/no-penetration conditions without iteration, significantly improving efficiency and accuracy. The CBS method can ensure the stability of the flow field using equal-order velocity-pressure interpolation, providing reliable results for accurate particle force and motion calculations. The validity and robustness of our proposed method are demonstrated by the benchmark problems. The method has been applied to simulate two-phase flow systems involving particles with generalised cross-sectional shapes. The results offer new insights into the influence of particle shape and initial orientation on particle-fluid interactions, highlighting the applicability of the method as a powerful tool for simulating particle-fluid interactions in engineering applications.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"145 ","pages":"Article 116132"},"PeriodicalIF":4.4,"publicationDate":"2025-04-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143843382","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":"A three-dimensional numerical model for the simulation of electrohydrodynamic varicose and whipping instabilities","authors":"Bin He,&nbsp;Yin Guan,&nbsp;Yanxiu Sha,&nbsp;Jingze Zheng,&nbsp;Yang Liu,&nbsp;Yihang Lei","doi":"10.1016/j.apm.2025.116133","DOIUrl":"10.1016/j.apm.2025.116133","url":null,"abstract":"<div><div>Electrohydrodynamic (EHD) cone-jet is a widely studied research topic owing to its variable outcomes and broad applications in micro/nanoscopic additive manufacturing. Among all the EHD cone-jet outcomes, the varicose and whipping instabilities attract much attention because of their rich underlying physics involved in the unstable downward stretching jet. However, the mechanisms behind these two instabilities are still not adequately understood, partly due to the lack of a full three-dimensional numerical model for the simulation of these two instabilities. In this paper, we propose a three-dimensional numerical model and perform a numerical study on EHD varicose and whipping instabilities under variations of electric voltage, liquid electrical conductivity, and liquid flow rate. The numerical model is devised with a two-phase finite-volume-method based flow solver in the open-source CFD program OpenFOAM. From the numerical results, we obtain the distributions of electric charge density and liquid velocity at the liquid surface and around the liquid jet at critical instants and locations. This enables us to analyze the differences in Taylor cone formation with the variation of the three parameters and compare the differences between liquid jet emission of the varicose and whipping instabilities. Finally, we compare our numerical results with the experimental findings in recently published papers, with very good agreement achieved regarding the scaling laws for the transition between these two instabilities. To the best of our knowledge, this is the first three-dimensional numerical work on the transition between EHD varicose and whipping instabilities, from which we intend to provide a new approach to analyzing these interesting but complicated phenomena.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"145 ","pages":"Article 116133"},"PeriodicalIF":4.4,"publicationDate":"2025-04-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143843323","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":"A novel Hybrid PSO-Heuristic Algorithm with Combinatorial Benders' Cuts for maximal evacuation planning in wildfire disasters","authors":"Saeide Bigdellou ,&nbsp;Qian Chen ,&nbsp;Saeed Beheshti","doi":"10.1016/j.apm.2025.116131","DOIUrl":"10.1016/j.apm.2025.116131","url":null,"abstract":"<div><div>Natural disasters, such as floods and fires, affect various regions of the world every year. One of the most critical aspects of disaster management and planning is facilitating the evacuation of people. Therefore, this study develops a mathematical model for emergency evacuation, taking into account the constraints and limitations of transferring individuals to shelters. The model, called the maximal bus evacuation planning, considers constraints including road blockages, vehicle fuel, and passenger capacities. Given the Nondeterministic Polynomial-time hard complexity of this problem, a Hybrid Particle Swarm Optimization–Heuristic Algorithm, which combines the Particle Swarm Optimization Algorithm with a Heuristic Algorithm, is used to solve the nonlinear model. Furthermore, we linearize the model and solve it using the developed Combinatorial Benders' Cuts approach in three versions. Numerical computations are examined under various scenarios, and the outputs of the Hybrid Particle Swarm Optimization – Heuristic algorithm reach suitable solutions within a short timeframe. Additionally, the results from solving the linearized model using the Combinatorial Benders' Cuts show that, in most cases, the execution time of the second version is better than the other versions.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"145 ","pages":"Article 116131"},"PeriodicalIF":4.4,"publicationDate":"2025-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143823551","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":"Refraction and reflection patterns of an inclined turbulent conductive jet induced by a local transverse magnetic field","authors":"Alexander Vakhrushev ,&nbsp;Ebrahim Karimi-Sibaki ,&nbsp;Menghuai Wu ,&nbsp;Abdellah Kharicha","doi":"10.1016/j.apm.2025.116124","DOIUrl":"10.1016/j.apm.2025.116124","url":null,"abstract":"<div><div>This study investigates the interaction between turbulent conductive jets and transverse magnetic fields, focusing on the effects of jet inclination angle and magnetic field strength. This phenomenon is relevant to both metallurgy and astrophysics. In metallurgy, complex localized magnetic fields are frequently used to control molten metal flows. Similarly, in astrophysics, a comparable effect occurs with relativistic jets (RJs) emitted by active galactic nuclei (AGN), where these jets can be deflected or refracted when encountering magnetic fields. In this study, it was observed that the jet undergoes refraction, meaning its initial direction changes when passing through a direct current (DC) magnetic field. This refraction is accompanied by expansion and fluctuations in the jet's mean flow. The degree of refraction was measured using a refractive index (RI), denoted as <span><math><mi>n</mi></math></span>, which was related to the Stuart (magnetic interaction) number <span><math><mi>N</mi></math></span>. For moderate Stuart numbers (<span><math><mrow><mi>N</mi><mo>≤</mo><mn>0.5</mn></mrow></math></span>), the refraction index followed a simple relationship, <span><math><mrow><mi>n</mi><mo>≈</mo><mn>1</mn><mo>−</mo><mi>N</mi></mrow></math></span>, regardless of the jet's inclination angle. At higher values (<span><math><mrow><mi>N</mi><mo>&gt;</mo><mn>0.5</mn></mrow></math></span>), refraction index varies with inclination due to the development of reverse flow zones adjacent to the main jet, a behavior previously reported by the authors. At certain critical <span><math><mi>N</mi></math></span> values, full reflection occurs, where the jet cannot penetrate the magnetic field. When the jet is confined by conductive walls, refraction intensifies, leading to earlier full reflection. These findings provide insight into the mechanisms behind the bending and expansion of astrophysical relativistic jets, suggesting that similar magnetic interactions may influence the RJs behavior as they encounter magnetic fields in space.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"145 ","pages":"Article 116124"},"PeriodicalIF":4.4,"publicationDate":"2025-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143829132","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":"Modelling for assessing the compulsory vaccination strategy against highly pathogenic avian influenza in a poultry farm of China","authors":"Juan Li ,&nbsp;Bingxin Nie ,&nbsp;Lu Gao ,&nbsp;Junhui Zhang ,&nbsp;Lijun Liu ,&nbsp;Hui Xu ,&nbsp;Youming Wang ,&nbsp;Renjie Sun ,&nbsp;Hongli Zhang ,&nbsp;Huaiping Zhu","doi":"10.1016/j.apm.2025.116130","DOIUrl":"10.1016/j.apm.2025.116130","url":null,"abstract":"<div><div>The prevalence of highly pathogenic avian influenza (HPAI) has become substantial obstacles for the poultry industry and international trade with potential threat to human health. We developed a data-driven transmission model for quantitatively assessing the risk of HPAI spread in a henhouse under the China national mandatory vaccination plan for prevention and control of HPAI in the poultry industry. With data from farm investigations and literature reviews, we calibrated our model to estimate some key model parameters. We employed the Latin Hypercubic Sampling (LHS) technique for multiple sampling during simulations and subsequently conducted statistical analysis to enhance the accuracy and reliability of our results. We found that within a designated single henhouse, regardless of whether the initial virus level is 5 or 500 times the minimum viral load required to trigger highly fatal avian influenza, when the virus transmission capacity is low (the basic reproductive number <span><math><mrow><msub><mi>R</mi><mn>0</mn></msub><mi>∈</mi><mrow><mo>(</mo><mrow><mn>2</mn><mo>,</mo><mn>3</mn></mrow><mo>)</mo></mrow></mrow></math></span>), implementing mandatory vaccination measures according to the policy-stipulated timings will keep the average overall mortality rate of chickens with maternal antibodies within the range of 0.08 %-0.42 %. However, if the virus transmission capacity doubles, the loss of chickens will increase by 3–5 times, but it will still be significantly lower than the normal mortality rate (5 %). Therefore, to optimize post-vaccination management and control, we recommend monitoring flock antibodies regularly, evaluating vaccine effectiveness promptly, strengthening pathogen surveillance, and upgrading farm biosecurity level. These measures are essential for safeguarding flock immunity, facilitating early virus detection, and mitigating the risk of HPAI transmission.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"145 ","pages":"Article 116130"},"PeriodicalIF":4.4,"publicationDate":"2025-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143826236","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":"Determining flux terms in a time fractional model","authors":"Mohamed BenSalah ,&nbsp;Salih Tatar ,&nbsp;Süleyman Ulusoy ,&nbsp;Masahiro Yamamoto","doi":"10.1016/j.apm.2025.116125","DOIUrl":"10.1016/j.apm.2025.116125","url":null,"abstract":"<div><div>In this paper, we investigate the flux identification problem for a nonlinear time-fractional viscoelastic equation with a general source function, using boundary measurements. We first establish the well-posedness of the direct problem by proving the existence, uniqueness, and continuous dependence of the solution on the heat flux. Next, we demonstrate the Fréchet differentiability of the cost functional, providing a theoretical foundation for solving the inverse problem. To efficiently reconstruct the unknown flux, we develop a Conjugate Gradient Algorithm based on the derived gradient formula. Finally, we validate the effectiveness and robustness of our approach through numerical experiments, including both noise-free and noisy data, confirming its accuracy and practical applicability.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"145 ","pages":"Article 116125"},"PeriodicalIF":4.4,"publicationDate":"2025-04-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143814997","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":"Adaptive mollified prescribed performance controller for waverider vehicle subjected to mismatched disturbances","authors":"Rongyi Guo ,&nbsp;Yibo Ding ,&nbsp;Ran Tao ,&nbsp;Xiaokui Yue","doi":"10.1016/j.apm.2025.116122","DOIUrl":"10.1016/j.apm.2025.116122","url":null,"abstract":"<div><div>Considering gliding waverider vehicle is often subjected to mismatched and matched disturbances, designing a strongly robust controller to accurately limit transient performance and ensure high-precision convergence of tracking error is of great significance for flight safety. Therefore, this paper investigates an adaptive mollified prescribed performance controller for waverider vehicle, including an adaptive anti-saturation prescribed performance function, and a mollified fixed-time sliding mode surface. Firstly, a novel adaptive anti-saturation prescribed performance function is designed for waverider vehicle to exactly restrict transient and steady-state performance of attitude, and reduce saturated duration of actuator without anti-saturation compensator. Meanwhile, compared with traditional function, the performance function can adaptively adjust boundary to ensure that tracking error is always enveloped within the performance function, thus avoiding severe singularity caused by error exceeding boundary. Then, a mollified fixed-time sliding mode surface is proposed to effectively suppress mismatched disturbances acting on gliding waverider vehicle, by introducing a novel designed mollified function. The mollified function can smooth sign function, making sliding mode surface continuously differentiable to reduce chattering. Based on the novel performance function, sliding mode surface, and mollified reaching law, the adaptive mollified prescribed performance controller can accurately restrict transient and steady-state performance of waverider vehicle, reduce saturation duration of actuator, and exhibit strong robustness against mismatched and matched disturbances. Finally, numerical simulations demonstrate the effectiveness of the proposed method.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"145 ","pages":"Article 116122"},"PeriodicalIF":4.4,"publicationDate":"2025-04-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143820740","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":"Analysis for face stability of the underwater shield tunnel traversing longitudinally heterogeneous strata","authors":"Jijie Du ,&nbsp;Liansheng Tang ,&nbsp;Lujia Niu ,&nbsp;Xianzhou Lyu ,&nbsp;Jianing Huang","doi":"10.1016/j.apm.2025.116116","DOIUrl":"10.1016/j.apm.2025.116116","url":null,"abstract":"<div><div>Maintaining tunnel face stability during underwater shield tunneling through longitudinally heterogeneous strata remains a critical challenge, with limited theoretical research addressing this issue. This study presents an improved limit equilibrium model to assess tunnel face stability, specifically considering the failure zones on both sides of the interface. The model's validity is demonstrated through comparative verification with prior studies and numerical simulations, and a comprehensive parametric analysis reveals critical insights into tunnel face stability. The findings reveal the limit support pressure (LSP) increases when the shield machine approaches a low-strength stratum from a high-strength one and decreases during the reverse process. Neither water depth nor burial depth do not affect the evolution of the LSP but contribute to its increase as these parameters increase. The support ratio (<em>K</em>) rises with water depth increases while decreases with deeper burial depth. Furthermore, the tunnel face-interface distance at which LSP initiates modification (<em>L_m</em>) demonstrates an inverse correlation with stratum strength, while being positively influenced by burial depth and water depth, with stratum strength dominating over the latter two parameters.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"145 ","pages":"Article 116116"},"PeriodicalIF":4.4,"publicationDate":"2025-04-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143823550","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":"Nonlinear static and dynamic buckling analysis of GPL-reinforced spherical caps and circular plates with porous core and stepped spiderweb stiffeners","authors":"Hoai Nam Vu ,&nbsp;Quang Minh Tran ,&nbsp;Van Tien Nguyen ,&nbsp;Thanh Hieu Pham ,&nbsp;Thi Phuong Nguyen","doi":"10.1016/j.apm.2025.116123","DOIUrl":"10.1016/j.apm.2025.116123","url":null,"abstract":"<div><div>This paper analyzes the nonlinear static and dynamic buckling responses of graphene platelet (GPL) reinforced shallow spherical caps and circular plates with porous core and stepped spiderweb stiffeners. The new design of stepped spiderweb stiffeners is proposed by adding meridian stiffeners near the edge region, and three regions are created as the edge region, the middle region, and the top region. Different GPL distribution laws, including U, X, O, A, and V laws are considered for GPL-reinforced face sheets, and the U, X, O, V, and A GPL distribution laws for spiderweb stiffeners are respectively chosen, while the trigonometric distribution law of porosity is designed for the porous core. The smeared stiffener technique is improved for stepped spiderweb stiffeners, combining with nonlinear Donnell shell theory, and nonlinear foundation model to derive the static and dynamic responses of the structures. The energy method is employed to perform the equilibrium equations for static problems, and the motion equation for dynamic problems. The static postbuckling responses are explicitly obtained, and the Runge-Kutta method is applied to investigate the nonlinear dynamic buckling behavior of the considered plates and shell caps. The numerical examples of nonlinear static and dynamic buckling of spherical caps and circular plates can predict the significant effects of spiderweb stiffeners, geometrical and material parameters, and nonlinear foundation. The results show that the spiderweb stiffeners largely enhance the buckling strength of spherical caps, while a higher porosity coefficient reduces mechanical but increases thermal postbuckling strength.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"145 ","pages":"Article 116123"},"PeriodicalIF":4.4,"publicationDate":"2025-04-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143820739","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":"A reduced model for phase-change problems with radiation using simplified PN approximations","authors":"Youssef Belhamadia ,&nbsp;Mohammed Seaid","doi":"10.1016/j.apm.2025.116108","DOIUrl":"10.1016/j.apm.2025.116108","url":null,"abstract":"<div><div>Radiative heat transfer in phase-change media is of great interest in many thermal applications in sciences and engineering involving internal melting or solidification. In these problems at high temperature, a mathematical model used to describe the heat transfer and phase change should also include equations accounting for thermal radiation. Using the integro-differential equation for the radiative intensity in these models results in a system of coupled equations for which its numerical solution is computationally very demanding. In the present study, we develop a class of efficient reduced models for phase-change problems accounting for grey thermal radiation. The novelty in these models lies in the fact that effects of thermal radiation are well captured in phase-change materials without solving the computationally demanding radiative transfer equation. The model is derived from the enthalpy formulation and the simplified <span><math><msub><mrow><mi>P</mi></mrow><mrow><mi>N</mi></mrow></msub></math></span> approximations of spherical harmonics. The integro-differential equation for the full radiative transfer is replaced by a set of differential equations which are independent of the angle variable and easy to solve using conventional computational methods. To solve the coupled equations, we implement a second-order implicit scheme for the time integration and a mixed finite element method for the space discretization. A Newton-based algorithm is also adopted for solving the nonlinear systems resulting from the considered monolithic approach. The performance of the proposed reduced models is analyzed on several test examples for coupled radiative heat transfer and phase-change problems in two and three space dimensions. The results presented in this study demonstrate that the proposed models can accurately predict the temperature distributions and capture the phase-change interfaces in melting and solidification examples, all while maintaining a very low computational cost.</div></div>","PeriodicalId":50980,"journal":{"name":"Applied Mathematical Modelling","volume":"145 ","pages":"Article 116108"},"PeriodicalIF":4.4,"publicationDate":"2025-04-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143792189","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}