Muhammad Usman , Muhammad Hamid , Dianchen Lu , Zhengdi Zhang , Wojciech Sumelka
{"title":"Robust iterative spectral algorithms for smooth solutions of time-fractional nonlinear diffusion problems and convergence analysis","authors":"Muhammad Usman , Muhammad Hamid , Dianchen Lu , Zhengdi Zhang , Wojciech Sumelka","doi":"10.1016/j.camwa.2024.10.015","DOIUrl":"10.1016/j.camwa.2024.10.015","url":null,"abstract":"<div><div>Nonlinear time-fractional diffusion problems, a significant class of parabolic-type problems, appear in various diffusion phenomena that seem extensively in nature. Such physical problems arise in numerous fields, such as phase transition, filtration, biochemistry, and dynamics of biological groups. Because of its massive involvement, its accurate solutions have become a challenging task among researchers. In this framework, this article proposed two operational-based robust iterative spectral schemes for accurate solutions of the nonlinear time-fractional diffusion problems. Temporal and spatial variables are approximated using Vieta-Lucas polynomials, and derivative operators are approximated using novel operational matrices. The approximated solution, novel operational matrices, and uniform collection points convert the problem into a system of nonlinear equations. Here, two robust methods, namely Picard's iterative and Newton's, are incorporated to tackle a nonlinear system of equations. Some problems are considered in authenticating the present methods' accuracy, credibility, and reliability. An inclusive comparative study demonstrates that the proposed computational schemes are effective, accurate, and well-matched to find the numerical solutions to the problems mentioned above. The proposed methods improve the accuracy of numerical solutions from 27 % to 100 % when <span><math><mrow><mi>M</mi><mo>></mo><mn>2</mn></mrow></math></span> as compared to the existing results. The suggested methods' convergence, error bound, and stability are investigated theoretically and numerically.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2024-10-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142538691","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":"Generalised adaptive cross approximation for convolution quadrature based boundary element formulation","authors":"A.M. Haider , S. Rjasanow , M. Schanz","doi":"10.1016/j.camwa.2024.10.025","DOIUrl":"10.1016/j.camwa.2024.10.025","url":null,"abstract":"<div><div>The acoustic wave equation is solved in time domain with a boundary element formulation. The time discretisation is performed with the generalised convolution quadrature method and for the spatial approximation standard lowest order elements are used. Collocation and Galerkin methods are applied. In the interest of increasing the efficiency of the boundary element method, a low-rank approximation such as the adaptive cross approximation (ACA) is carried out. We discuss a generalisation of the ACA to approximate a three-dimensional array of data, i.e., usual boundary element matrices at several complex frequencies. This method is used within the generalised convolution quadrature (gCQ) method to obtain a real time domain formulation. The behaviour of the proposed method is studied with three examples, a unit cube, a unit cube with a reentrant corner, and a unit ball. The properties of the method are preserved in the data sparse representation and a significant reduction in storage is obtained.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2024-10-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142530310","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":"Innovative discretizations of PDEs: Towards an accurate representation of the reality","authors":"Fleurianne Bertrand, Daniele Boffi, Alexander Düster, Jean-Luc Guermond, Norbert Heuer, Jichun Li, Waldemar Rachowicz","doi":"10.1016/j.camwa.2024.10.013","DOIUrl":"10.1016/j.camwa.2024.10.013","url":null,"abstract":"","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2024-10-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142528145","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-priori and a-posteriori error estimates for discontinuous Galerkin method of the Maxwell eigenvalue problem","authors":"Jun Zhang , Zijiang Luo , Jiayu Han , Hu Chen","doi":"10.1016/j.camwa.2024.10.026","DOIUrl":"10.1016/j.camwa.2024.10.026","url":null,"abstract":"<div><div>This paper is devoted to a-priori and a-posteriori error analysis of discontinuous Galerkin (DG) method for the Maxwell eigenvalue problem. The discrete compactness of DG space is proved so that the Babuška and Osborn spectral approximation theory can be applicable in the a-priori error analysis. Then we prove the optimal error estimates for DG eigenfunctions in mesh-dependent norm and DG eigenvalues. A special contribution of this work is to prove that the error in <span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span>-norm for smooth eigenfunctions is of higher order than that in mesh-dependent norm, so that the DG eigenvalues can approximate the true eigenvalues from upper. Another contribution of this work is to provide a-posteriori error analysis for the DG method. A reliable a-posteriori error estimator is analyzed. The upper bound property of DG eigenvalues and the robustness of adaptive methods are verified through numerical experiments.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2024-10-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142528061","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}
Yujie Wang , Shengquan Wang , Yong Cai , Guidong Wang , Guangyao Li
{"title":"Fully parallel and pipelined sparse direct solver for large symmetric indefinite finite element problems","authors":"Yujie Wang , Shengquan Wang , Yong Cai , Guidong Wang , Guangyao Li","doi":"10.1016/j.camwa.2024.10.017","DOIUrl":"10.1016/j.camwa.2024.10.017","url":null,"abstract":"<div><div>Sparse linear system solving is a primary computational cost in large-scale finite element analysis, and improving its performance is a key technological challenge in this field. Real-world engineering problems involve diverse materials, elements, and connectivity relationships, making it difficult for iterative methods to handle their global stiffness matrices. Direct methods, owing to their robustness, emerge as the preferred choice. In this paper, a novel block-based supernodal LDL<sup>T</sup> numerical factorization method is introduced. The computational process is disassembled into distinct tasks, and the dependency relationships between these tasks are expressed via a directed acyclic graph to guide the calculation sequence. Based on this approach, a global task pool and local task stack are established to store task queues, enhancing data reuse and multicore collaboration efficiency. Additionally, an effective task dispatch and work-stealing mechanism is implemented to prevent performance degradation caused by load imbalances. Numerical experiments, including a publicly available matrix test set and real-world engineering finite element problems, are conducted to compare the parallel performances of the Pardiso, MUMPS, and proposed solver. The results illustrate that the proposed solver performs significantly better than the other solvers when handling various types of sparse matrices and diverse architectures of multicore processors.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2024-10-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142529732","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 least-squares Fourier frame method for nonlocal diffusion models on arbitrary domains","authors":"Mengxia Shen , Haiyong Wang","doi":"10.1016/j.camwa.2024.10.024","DOIUrl":"10.1016/j.camwa.2024.10.024","url":null,"abstract":"<div><div>We introduce a least-squares Fourier frame method for solving nonlocal diffusion models with Dirichlet volume constraint on arbitrary domains. The mathematical structure of a frame rather than a basis allows using a discrete least-squares approximation on irregular domains and imposing non-periodic boundary conditions. The method has inherited the one-dimensional integral expression of Fourier symbols of the nonlocal diffusion operator from Fourier spectral methods for any <em>d</em> spatial dimensions. High precision of its solution can be achieved via a direct solver such as pivoted QR decomposition even though the corresponding system is extremely ill-conditioned, due to the redundancy in the frame. The extension of AZ algorithm improves the complexity of solving the rectangular linear system to <span><math><mi>O</mi><mo>(</mo><mi>N</mi><msup><mrow><mi>log</mi></mrow><mrow><mn>2</mn></mrow></msup><mo></mo><mi>N</mi><mo>)</mo></math></span> for 1<em>d</em> problems and <span><math><mi>O</mi><mo>(</mo><msup><mrow><mi>N</mi></mrow><mrow><mn>2</mn></mrow></msup><msup><mrow><mi>log</mi></mrow><mrow><mn>2</mn></mrow></msup><mo></mo><mi>N</mi><mo>)</mo></math></span> for 2<em>d</em> problems, compared with <span><math><mi>O</mi><mo>(</mo><msup><mrow><mi>N</mi></mrow><mrow><mn>3</mn></mrow></msup><mo>)</mo></math></span> of the direct solvers, where <em>N</em> is the number of degrees of freedom. We present ample numerical experiments to show the flexibility, fast convergence and asymptotical compatibility of the proposed method.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2024-10-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142530295","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 local projection stabilised HHO method for the Oseen problem","authors":"Gouranga Mallik , Rahul Biswas , Thirupathi Gudi","doi":"10.1016/j.camwa.2024.10.030","DOIUrl":"10.1016/j.camwa.2024.10.030","url":null,"abstract":"<div><div>In this article, we consider a local projection stabilisation for a Hybrid High-Order (HHO) approximation of the Oseen problem. We prove an existence-uniqueness result under a stronger SUPG-like norm. We improve the stability and provide error estimation in stronger norm for convection dominated Oseen problem. We also derive an optimal order error estimate under the SUPG-like norm for equal-order polynomial discretisation of velocity and pressure spaces. Numerical experiments are performed to validate the theoretical results.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2024-10-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142528144","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":"Accurate and parallel simulation of the anisotropic dendrite crystal growth by Lagrangian data assimilation with directional operator splitting","authors":"Fenglian Zheng , Yan Wang , Xufeng Xiao","doi":"10.1016/j.camwa.2024.10.020","DOIUrl":"10.1016/j.camwa.2024.10.020","url":null,"abstract":"<div><div>Dendritic crystal growth is a prevalent natural phenomenon that generates a crystalline structure resembling a tree during a phase transition. In practical computations, inaccuracies in model parameters and initial conditions can introduce observation errors, even leading to inaccurate results. To enhance the precision and efficiency of the numerical simulation, the data assimilation method with parallel simulation are considered in this study. Firstly, based on the phase-field dendritic crystal growth model, a Lagrangian data assimilation method, which adds Lagrange multiplier terms into the phase-field partial differential equations (PDEs), is presented to integrate the observed data of physical information to modify the numerical solution, thereby improving simulation accuracy. Secondly, to achieve efficient data assimilation, a parallel directional operator splitting method is presented to solve the modified data assimilation PDEs. Thirdly, in the section of numerical experiments, we investigate the validity of the method and assess the impact of various factors such as the Lagrange multiplier parameter, spatio-temporal sampling rate and parameter perturbation ratio on the effectiveness of data assimilation. The evaluation is conducted for two distinct problem categories: initial observation errors and model parameter errors. Experimental results demonstrate that our method can effectively assimilate experimental observations in simulations, thereby enhancing more accurate dendritic crystal growth processes.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2024-10-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142530294","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 peridynamic modeling of damage and penetration in composite plates exposed to localized explosive blasts","authors":"D.A. Abdoh","doi":"10.1016/j.camwa.2024.10.022","DOIUrl":"10.1016/j.camwa.2024.10.022","url":null,"abstract":"<div><div>Localized failures in composites can lead to catastrophic consequences in some vehicles and systems, such as airplanes and submarines. The paper presents a new three-dimensional (3D) peridynamic model to explore the localized explosive effects on composite plates. The study utilizes the peridynamic method to simulate fractures and deformations in composite plates when exposed to localized explosive blasts. Accurate prediction of composites' performance during explosive blast events is crucial in the design process to avoid the deadly effects of their failures. This study highlights the following novelties: (1) We present, for the first time, a novel 3D mesh-free model to simulate the fracture and damage behavior of composite plates when exposed to explosive blasts; (2) The adopted numerical modeling technique enables highly efficient simulations of fractures and failures in composites when compared with other mesh-based numerical models; (3) We introduce a new mathematical framework to reflect the explosive pressure loads through different composite layers; (4) The study provides an accessible and efficient tool for engineers and researchers to enhance the design of composites in several industries instead relying on limited-access commercial software packages. In addition to the previous novelties, the paper presents a new parametric study that investigates the performance of several composite plates, such as titanium-aramid and aluminum-aramid composites, in explosive blast scenarios. Moreover, the roles of explosive mass and plate thickness in the failure mechanisms of composites are examined.</div></div>","PeriodicalId":55218,"journal":{"name":"Computers & Mathematics with Applications","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2024-10-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142530293","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}