Computers & Mathematics with Applications最新文献

{"title":"A hybrid reduced-order model for segregated fluid-structure interaction solvers in an ALE approach at high Reynolds number","authors":"Valentin Nkana Ngan ,&nbsp;Giovanni Stabile ,&nbsp;Andrea Mola ,&nbsp;Gianluigi Rozza","doi":"10.1016/j.camwa.2025.01.004","DOIUrl":"10.1016/j.camwa.2025.01.004","url":null,"abstract":"<div><div>This study introduces a first step for constructing a hybrid reduced-order models (ROMs) for segregated fluid-structure interaction in an Arbitrary Lagrangian-Eulerian (ALE) approach at a high Reynolds number using the Finite Volume Method (FVM). The ROM is driven by proper orthogonal decomposition (POD) with hybrid techniques that combines the classical Galerkin projection and two data-driven methods (radial basis networks, and neural networks/ long short term memory). Results demonstrate the ROM's ability to accurately capture the physics of fluid-structure interaction phenomena. This approach is validated through a case study focusing on flow-induced vibration (FIV) of a pitch-plunge airfoil at a high Reynolds number (<span><math><mi>R</mi><mi>e</mi><mo>=</mo><msup><mrow><mn>10</mn></mrow><mrow><mn>7</mn></mrow></msup></math></span>).</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"180 ","pages":"Pages 299-321"},"PeriodicalIF":2.9,"publicationDate":"2025-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143055382","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":"Block ω-circulant preconditioners for parabolic equations","authors":"Po Yin Fung,&nbsp;Sean Y. Hon","doi":"10.1016/j.camwa.2025.01.019","DOIUrl":"10.1016/j.camwa.2025.01.019","url":null,"abstract":"<div><div>In this study, a novel class of block <em>ω</em>-circulant preconditioners is developed for the all-at-once linear system that emerges from solving parabolic equations using first and second order discretization schemes for time. We establish a unifying preconditioning framework for <em>ω</em>-circulant preconditioners, extending and modifying the preconditioning approach recently proposed in (Zhang and Xu, 2024 <span><span>[27]</span></span>) and integrating some existing results in the literature. The proposed preconditioners leverage fast Fourier transforms for efficient diagonalization, facilitating parallel-in-time execution. Theoretically, these preconditioners ensure that eigenvalue clustering around ±1 is achieved, fostering fast convergence under the minimal residual method. Furthermore, when using the generalized minimal residual method, the effectiveness of these preconditioners is supported by the singular value clustering at unity. Numerical experiments validate the performance of the developed preconditioning strategies.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"182 ","pages":"Pages 122-138"},"PeriodicalIF":2.9,"publicationDate":"2025-01-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143020141","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 staggered discontinuous Galerkin method for solving SN transport equation on arbitrary polygonal grids","authors":"Deng Wang ,&nbsp;Zupeng Jia","doi":"10.1016/j.camwa.2025.01.018","DOIUrl":"10.1016/j.camwa.2025.01.018","url":null,"abstract":"<div><div>This paper proposes a staggered discontinuous Galerkin (SDG) method for solving the 2D <span><math><msub><mrow><mi>S</mi></mrow><mrow><mi>N</mi></mrow></msub></math></span> transport equation on arbitrary polygonal mesh. The new method allows rough grids such as highly distorted quadrilateral grids and general polygonal grids. More importantly, it is numerical flux free, different from the standard discontinuous Galerkin (DG) method using upwind flux, and thus gains the advantage of not demanding a sweeping algorithm to determine the computational ordering of all elements. The sweeping process not only requires a significant amount of computation, but also encounters deadlocks due to the presence of cycles in the corresponding directed graph on deformed three-dimensional polyhedral meshes. Additionally, there are challenges with sweeping stability on curved meshes. Our method naturally avoids these problems such that it can be generalized to high-order schemes on curved meshes. Convergence of the new method is analyzed in the linear case, and we have shown that it is optimally convergent for sufficiently smooth transport solution and appropriate total cross section. Numerical results are consistent with the theoretical analysis. The tests also show that our method employing linear elements can maintain second-order accuracy on rough grids mentioned earlier, demonstrating good robustness, and be able to model material interfaces with sharp changes of the angular flux. Moreover, the asymptotic-preserving property can be observed by scaling the cross section parameter in the thick diffusive limit problem. A test employing quadratic elements is also provided, which preliminarily demonstrates the feasibility and effectiveness of our method in higher-order scenarios.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"182 ","pages":"Pages 102-121"},"PeriodicalIF":2.9,"publicationDate":"2025-01-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143020191","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 handy tool for assessing tetrahedron-based finite-cell methods and for numerical simulations in spheroidal domains","authors":"Andrés León Baldelli ,&nbsp;Vitoriano Ruas ,&nbsp;Marco Antonio Silva Ramos","doi":"10.1016/j.camwa.2025.01.017","DOIUrl":"10.1016/j.camwa.2025.01.017","url":null,"abstract":"<div><div>A straightforward procedure is presented for the generation of finite-cell meshes consisting of tetrahedrons for curved domains, whose boundary can be expressed in spherical coordinates with origin at a suitable location in its interior. Besides the equation of the boundary, the generation of the mesh depends only on an integer parameter, whose value is associated with its degree of refinement. Several examples indicate that the meshes of a given domain form a quasi-uniform family of partitions, as the value of the integer parameter increases. Mesh quality is optimal in the case of a ball but it remains quite correct as the shape of the domain moves away from perfect sphericity, with a gradual but in all natural downgrade. The procedure is a handy tool for an a priori order-checking of a new finite-cell method, as applied to a given type of boundary value problem posed in curved domains. A MATLAB code was developed to implement this tetrahedrization procedure for domains with three symmetry planes.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"181 ","pages":"Pages 323-339"},"PeriodicalIF":2.9,"publicationDate":"2025-01-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143020142","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":"Simulation of cartesian cut-cell technique for modeling turbulent flow in asymmetric diffusers using various turbulence models","authors":"A.S. Dawood ,&nbsp;A.S. Amer ,&nbsp;R.M. Abumandour ,&nbsp;W.A. El-Askary","doi":"10.1016/j.camwa.2025.01.015","DOIUrl":"10.1016/j.camwa.2025.01.015","url":null,"abstract":"<div><div>The process of accurately predicting the behavior of the separated turbulent flow requires extraordinary efforts, whether it is choosing the appropriate computational mesh or using the appropriate turbulence model for that. The present study introduces a comparative numerical investigation for predicting the behavior of turbulent-separated flow in asymmetric diffusers. Numerical simulation using a finite volume approach of incompressible Reynolds Averaged Navier Stokes equations (RANS) with three turbulence models (Standard <span><math><mrow><mi>k</mi><mo>−</mo><mrow><mi>ε</mi></mrow></mrow></math></span>, Chen-kim, and modified Chen-kim) is here performed in a self-developed FORTRAN code. The treatment of asymmetric diffusers poses challenges due to the complex flow behavior and geometry. To address this, a developed cartesian cut-cell technique is employed, which provides compatibility with solid boundaries and efficiently handles complex geometries This developed cut-cell technique is checked to treat its ability to predict complex turbulent flow with the presence of strong pressure gradient for correctness and convergence, as well as testing the proposed turbulence-models performance. So, verifications are performed by comparing the present computational results of asymmetric-diffuser flow characteristics with available experimental and LES data. The proposed models reveal acceptable agreements in most cases, especially the modified Chen-kim model, which shows a great match with the experimental and LES results for all flow characteristics. The standard k-ε model fails to predict the flow-separation well in most comparisons. Extended computational studies are also introduced to investigate the effects of diffuser cant angles (4 to 15°) and area ratio (2.4 to 7) on the diffuser flow behavior using the successfully modified Chen-kim turbulence model. The parametric study reveals that these two factors strongly affect the diffuser performance, where the pressure recovery, skin friction coefficients, and separation bubble size are provided.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"182 ","pages":"Pages 84-101"},"PeriodicalIF":2.9,"publicationDate":"2025-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143020188","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 of the Picard-Newton finite element iteration for the stationary incompressible inductionless MHD equations","authors":"Xiaodi Zhang ,&nbsp;Meiying Zhang ,&nbsp;Xianghai Zhou","doi":"10.1016/j.camwa.2025.01.016","DOIUrl":"10.1016/j.camwa.2025.01.016","url":null,"abstract":"<div><div>In this paper, we propose and analyze the Picard-Newton finite element iteration for the stationary incompressible inductionless magnetohydrodynamics (MHD) equations. In finite element discretization, the hydrodynamic unknowns are approximated by stable finite element pairs, and the electromagnetic system is discretized by using the face-volume pairs. To solve the nonlinear discretized problem efficiently, our method consists of first applying the Picard iteration and then applying the Newton iteration. The Picard-Newton iteration is proved to be globally stable under the uniqueness condition and quadratically convergent under the stronger uniqueness condition. Thanks to the improved stability property, this solver has a larger convergence basin than the usual Newton iteration. Numerical tests confirm our theoretical analysis and show that the Picard-Newton iteration dramatically excels both the Picard and Newton iterations in several benchmark problems.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"182 ","pages":"Pages 46-65"},"PeriodicalIF":2.9,"publicationDate":"2025-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143020189","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":"Efficient reliability analysis method for non-linear truss structures using machine learning-based uncertainty quantification","authors":"Trung-Hieu Nguyen,&nbsp;Truong-Thang Nguyen,&nbsp;Duc-Minh Hoang,&nbsp;Viet-Hung Dang,&nbsp;Xuan-Dat Pham","doi":"10.1016/j.camwa.2025.01.014","DOIUrl":"10.1016/j.camwa.2025.01.014","url":null,"abstract":"<div><div>Truss structures typically involve a large number of similar elements; hence, it is necessary to employ reliability analysis algorithms that can handle high-dimensional problems to analyze the reliability of truss structures. Moreover, when considering non-linear behaviors in terms of both material properties and geometry, developing such an algorithm is challenging. For this purpose, this study proposes a novel method, named t-LQR that combines the advancements from three domains: i) a high-performance gradient boosting model from machine learning for a highly accurate prediction model, ii) an active learning process from reliability analysis for adaptively improving the prediction model, and iii) quantile regression for uncertainty quantification from probabilistic information to identify the relevant candidates used to refine the prediction model. The validity and robustness of the proposed method are verified through planar and spatial truss structures, showing that t-LQR significantly reduces the computational time of structural analysis-up to 25 times-compared to the conventional Monte Carlo methods. Furthermore, t-LQR outperforms competing Kirging-based models in terms of accuracy for non-linear problems.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"182 ","pages":"Pages 66-83"},"PeriodicalIF":2.9,"publicationDate":"2025-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143020190","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":"The mechanism analysis of the oscillatory flow of Oldroyd-B fluid in a tube with a right triangular cross-section","authors":"Lin Liu ,&nbsp;Sen Zhang ,&nbsp;Lang Liu ,&nbsp;Chiyu Xie ,&nbsp;Lixiu Wang ,&nbsp;Yijia Liu ,&nbsp;Jing Zhu ,&nbsp;Libo Feng","doi":"10.1016/j.camwa.2024.12.001","DOIUrl":"10.1016/j.camwa.2024.12.001","url":null,"abstract":"<div><div>This paper considers the oscillatory flow of the Oldroyd-B fluid in a tube with a right triangular cross-section. The partial differential equation for describing the unidirectional flow of Oldroyd-B fluid is derived. For treating the triangular region, the unstructured mesh finite element method is applied. For verifying the accuracy of the finite element method, the source term is introduced and the exact solution is constructed. We obtain an analytical solution for the isosceles triangular cross-section tube and analyze the effects of different parameters on the velocity amplitude and phase difference of the Oldroyd-B fluid flow. In addition, the dynamic characteristics for the unsteady-state Oldroyd-B fluid flow are discussed and the influence of the involved parameters, such as the pressure gradient, the relaxation and retardation time parameters on the velocity distribution are also investigated and the flow mechanism is analyzed graphically. At the same condition, we find that only the tube with an isosceles right triangle cross-section can produce the oscillatory flow. Finally, the fluid mechanism between the Newtonian fluid and the Oldroyd-B fluid for special cases is analyzed and discussed in detail.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"178 ","pages":"Pages 179-196"},"PeriodicalIF":2.9,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142804451","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 of tumor metastasis model in microenvironment based on coupled fractional reaction diffusion equation","authors":"Yating Huang,&nbsp;Zhenyou Wang","doi":"10.1016/j.camwa.2025.01.012","DOIUrl":"10.1016/j.camwa.2025.01.012","url":null,"abstract":"<div><div>This paper explores the licensing stage in tumor metastasis, focusing on the mathematical relationships between the tumor microenvironment and cells. Building on Academician Cao's insights, we develop a fractional order spatiotemporal model using a reaction-diffusion approach, solving it numerically to account for tumor metastasis memory and heritability. Our analysis identifies causal relationships during this stage, providing essential theoretical support and tools for advancing cancer treatment research. This study aims to innovate tumor therapy methods, offering new insights and strategies.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"182 ","pages":"Pages 24-45"},"PeriodicalIF":2.9,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143020192","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":"In-plane vibration analysis of elastically restrained FGM skew plates using variational differential quadrature method","authors":"C. Li ,&nbsp;S. Li ,&nbsp;Y.M. Zhang ,&nbsp;J.W. Cai ,&nbsp;S.K. Lai","doi":"10.1016/j.camwa.2024.11.026","DOIUrl":"10.1016/j.camwa.2024.11.026","url":null,"abstract":"<div><div>This work presents an accurate in-plane vibration analysis of functionally graded material (FGM) skew plates with elastically restrained boundaries using the variational differential quadrature method (VDQM). The weak form of the governing equations is derived by integrating two-dimensional elasticity theory with Hamilton's principle. The differential and integral operators are directly converted into matrix forms, thereby removing the necessity for higher-order derivative approximations in the displacement field. Transformation matrices are then developed for both interior and boundary nodes to link the governing equations with the boundary conditions, leading to the formulation of the vibration eigenvalue equations for FGM skew plates. Various factors, including aspect ratios, skew angles, boundary restraints, and gradient indices, are considered to investigate the in-plane vibration mode characteristics of FGM skew plates. Detailed solution procedures are provided, and numerical examples using the proposed solutions indicate that the VDQM exhibits superior numerical convergence and stability compared to other existing methods. The study also investigates the influence of highly skewed plates (75°), where stress singularities arise at the corners. This aspect is crucial for in-plane vibration analysis and has garnered limited attention in the existing literature. Furthermore, the dynamic analysis of FGM skew plates reveals a coupling between normal and tangential vibration modes.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":"178 ","pages":"Pages 136-153"},"PeriodicalIF":2.9,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142762994","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}