{"title":"An alternating shifted higher order power method based algorithm for rank-R Hermitian approximation and solving Hermitian CP-decomposition problems","authors":"Xiaofang Xin, Guyan Ni, Ying Li","doi":"10.1016/j.cam.2024.116385","DOIUrl":"10.1016/j.cam.2024.116385","url":null,"abstract":"<div><div>The Hermitian tensor is a higher order extension of the Hermitian matrix that can be used to represent quantum mixed states and solve problems such as entanglement and separability of quantum mixed states. In this paper, we propose a novel numerical algorithm, an alternating shifted higher order power method (AS-HOPM), for rank-<span><math><mi>R</mi></math></span> Hermitian approximation, which can also be used to compute Hermitian Candecomp/Parafac (CP) decomposition. At the same time, for the choice of initial points, we give a Broyden–Fletcher–Goldfarb–Shanno (BFGS) method based on unconstrained optimization, and propose a BFGS-AS-HOPM algorithm for rank-<span><math><mi>R</mi></math></span> Hermitian approximation. For solving the Hermitian CP-decomposition problem, numerical experiments show that using the BFGS-AS-HOPM algorithm has a higher success rate than using the AS-HOPM algorithm alone.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"459 ","pages":"Article 116385"},"PeriodicalIF":2.1,"publicationDate":"2024-11-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142699826","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":"Decoupled and energy stable schemes for phase-field surfactant model based on mobility operator splitting technique","authors":"Nan Lu , Chenxi Wang , Lun Zhang , Zhen Zhang","doi":"10.1016/j.cam.2024.116365","DOIUrl":"10.1016/j.cam.2024.116365","url":null,"abstract":"<div><div>In this paper, we investigate numerical methods for the phase-field surfactant (PFS) model, which is a gradient flow system consisting of two nonlinearly coupled Cahn-Hilliard type equations. The main challenge in developing high-order efficient energy stable methods for this system results from the nonlinearity and the strong coupling in the two variables in the free energy functional. We propose two fully decoupled, linear and energy stable schemes based on a linear stabilization approach and an operator splitting technique. We rigorously prove that both schemes can preserve the original energy dissipation law. The techniques employed in these schemes are then summarized into an innovative approach, which we call the mobility operator splitting (MOS), to design high-order decoupled energy stable schemes for a wide class of gradient flow systems. As a particular case, MOS allows different time steps for updating respective variables, leading to a multiple time-stepping strategy for fast-slow dynamics and thus serious improvement of computational efficiency. Various numerical experiments are presented to validate the accuracy, efficiency and other desired properties of the proposed schemes. In particular, detailed phenomena in thin-film pinch-off dynamics can be clearly captured by using the proposed schemes.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"459 ","pages":"Article 116365"},"PeriodicalIF":2.1,"publicationDate":"2024-11-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142699970","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}
Shipra Mahata , Samala Rathan , Juan Ruiz-Álvarez , Dionisio F. Yáñez
{"title":"A general correction for numerical integration rules over piece-wise continuous functions","authors":"Shipra Mahata , Samala Rathan , Juan Ruiz-Álvarez , Dionisio F. Yáñez","doi":"10.1016/j.cam.2024.116378","DOIUrl":"10.1016/j.cam.2024.116378","url":null,"abstract":"<div><div>This article presents a novel approach to enhance the accuracy of classical quadrature rules by incorporating correction terms. The proposed method is particularly effective when the position of an isolated discontinuity in the function and the jump in the function and its derivatives at that position are known. Traditional numerical integration rules are exact for polynomials of certain degree. However, they may not provide accurate results for piece-wise polynomials or functions with discontinuities without modifying the location and number of data points in the formula. Our proposed correction terms address this limitation, enabling the integration rule to conserve its accuracy even in the presence of a jump discontinuity. The numerical experiments that we present support the theoretical results obtained.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"459 ","pages":"Article 116378"},"PeriodicalIF":2.1,"publicationDate":"2024-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142699823","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":"An efficient second-order predictor–corrector infeasible primal–dual IPM algorithm with large iteration path updates for solving well-known SDO problems","authors":"Hadis Abedi, Behrouz Kheirfam","doi":"10.1016/j.cam.2024.116379","DOIUrl":"10.1016/j.cam.2024.116379","url":null,"abstract":"<div><div>In this paper, we propose a second-order predictor–corrector infeasible interior-point algorithm for semidefinite optimization in a new large neighborhood. The new large neighborhood, which is based on the spectral norm, is wider than the popular large neighborhoods based on the negative pseudo-infinity norm and the Frobenius norm. In each iteration, our algorithm calculates a new predictor direction using two modified systems and Yang et al. strategy. Then, this algorithm calculates a second-order corrector direction using the directions obtained in the predictor step. The iterates are determined by taking the largest possible step lengths along the search directions within the new large neighborhood. We prove that the algorithm is globally convergent and has <span><math><mrow><mi>O</mi><mrow><mo>(</mo><msup><mrow><mi>n</mi></mrow><mrow><mfrac><mrow><mn>5</mn></mrow><mrow><mn>4</mn></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>1</mn></mrow><mrow><mi>q</mi></mrow></mfrac></mrow></msup><mo>log</mo><msup><mrow><mi>ɛ</mi></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow></mrow></math></span> iteration complexity bound. Finally, the numerical experiments of the proposed algorithm confirm the efficiency and reliability of this approach.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"459 ","pages":"Article 116379"},"PeriodicalIF":2.1,"publicationDate":"2024-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660937","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 simulation of the generalized modified Benjamin–Bona–Mahony equation using SBP-SAT in time","authors":"Vilma Kjelldahl, Ken Mattsson","doi":"10.1016/j.cam.2024.116377","DOIUrl":"10.1016/j.cam.2024.116377","url":null,"abstract":"<div><div>In this paper we present high-order accurate finite difference approximations for solving the generalized modified Benjamin–Bona–Mahony (BBM) equation, a non-linear soliton model. The spatial discretization uses high-order accurate summation-by-parts (SBP) finite difference operators combined with both weak and strong enforcement of boundary conditions. For time integration we compare the explicit RK4 method against an implicit SBP time integrator. These time-marching methods are evaluated and compared in terms of accuracy and efficiency. It is shown that the implicit SBP time-integrator is more efficient than the explicit RK4 method for non-linear soliton models.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"459 ","pages":"Article 116377"},"PeriodicalIF":2.1,"publicationDate":"2024-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142661499","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":"Greedy Randomized Kaczmarz with momentum method for nonlinear equation","authors":"Li Liu, Wei-Guo Li, Li-Li Xing, Wen-Di Bao","doi":"10.1016/j.cam.2024.116359","DOIUrl":"10.1016/j.cam.2024.116359","url":null,"abstract":"<div><div>Based on the nonlinear greedy randomized Kaczmarz (NGRK) method and the heavy ball method, a nonlinear greedy randomized Kaczmarz with momentum (NGRKm) method is presented for solving large nonlinear systems. Furthermore, the convergence of the NGRKm method is proved and the corresponding numerical experiments show the effectiveness of the NGRKm method.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"459 ","pages":"Article 116359"},"PeriodicalIF":2.1,"publicationDate":"2024-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142699822","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}
Djulustan Nikiforov , Leonardo A. Poveda , Dmitry Ammosov , Yesy Sarmiento , Juan Galvis
{"title":"Meshfree generalized multiscale exponential integration method for parabolic problems","authors":"Djulustan Nikiforov , Leonardo A. Poveda , Dmitry Ammosov , Yesy Sarmiento , Juan Galvis","doi":"10.1016/j.cam.2024.116367","DOIUrl":"10.1016/j.cam.2024.116367","url":null,"abstract":"<div><div>This paper considers flow problems in multiscale heterogeneous porous media. The multiscale nature of the modeled process significantly complicates numerical simulations due to the need to compute huge and ill-conditioned sparse matrices, which negatively affect both the computational cost and the stability of the numerical solution. We propose a novel combined approach of the meshfree Generalized Multiscale Finite Element Method (MFGMsFEM) and exponential time integration for solving such problems. MFGMsFEM provides a robust and efficient spatial approximation, allowing us to consider complex heterogeneities without constructing a coarse computational grid. At the same time, using the cost-effective MFGMsFEM matrix, exponential integration provides a robust temporal approximation for stiff multiscale problems, allowing larger time steps. For the proposed multiscale approach, we provide a rigorous convergence analysis, including the new analysis of the MFGMsFEM spatial approximation. We conduct numerical experiments to computationally verify the proposed approach by solving linear and semi-linear flow problems in multiscale media. Numerical results demonstrate that the proposed multiscale method achieves significant reductions in computational cost and improved stability, even with larger time steps, confirming the theoretical analysis.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"459 ","pages":"Article 116367"},"PeriodicalIF":2.1,"publicationDate":"2024-11-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142661497","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 combined mixed finite element method and discontinuous Galerkin method for hybrid-dimensional fracture models of two-phase flow","authors":"Shuangshuang Chen","doi":"10.1016/j.cam.2024.116373","DOIUrl":"10.1016/j.cam.2024.116373","url":null,"abstract":"<div><div>In this paper, a combined numerical method consisting of the mixed finite element method (MFE) for the pressure equation and the discontinuous Galerkin (DG) method for the saturation equation is proposed to solve hybrid-dimensional fracture models of incompressible two-phase flow in porous media. The hybrid-dimensional fracture models treat fractures as <span><math><mrow><mo>(</mo><mi>d</mi><mo>−</mo><mn>1</mn><mo>)</mo></mrow></math></span>-dimensional interfaces immersed in <span><math><mi>d</mi></math></span>-dimensional matrix domains and take fluid exchange between fractures and surrounding matrix into account. Fully implicit approximation schemes combining the MFE-DG method with the backward Euler time discretization for the models with both a single fracture and an intersecting fractures network are all formulated successfully. The stability of the discrete solution is analyzed, and optimal error estimates in <span><math><mrow><mi>H</mi><mrow><mo>(</mo><mi>d</mi><mi>i</mi><mi>v</mi><mo>)</mo></mrow></mrow></math></span>-norm for the velocity and in <span><math><msup><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span>-norm for the pressure are derived, as well as in the discrete <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span>-norm for the saturation. Numerical experiments with a single fracture and a T-junction intersecting fractures network are conducted to verify the accuracy of our theoretical analysis.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"459 ","pages":"Article 116373"},"PeriodicalIF":2.1,"publicationDate":"2024-11-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660935","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 analysis of discontinuous Galerkin method for Electrical Impedance Tomography with partial data","authors":"Xiaosheng Li, Wei Wang","doi":"10.1016/j.cam.2024.116376","DOIUrl":"10.1016/j.cam.2024.116376","url":null,"abstract":"<div><div>In this work, we extend our previous work (Li and Wang, 2023) of the discontinuous Galerkin (DG) method for Electrical Impedance Tomography (EIT) with full data to partial data where the current and voltage measurements are taken only on part of the boundary. Additionally, we provide the convergence analysis of the DG approximation for EIT for partial data based on an iterative method with Tikhonov regularization. We prove that the minimizers of the discrete optimization problems converge to a minimizer of the continuous optimization problem as the mesh sizes of the discretization approach zero. Numerical results for the recovery of conductivities are tested with different types of partial data. The partial data results are demonstrated to be comparable to the full data results.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"459 ","pages":"Article 116376"},"PeriodicalIF":2.1,"publicationDate":"2024-11-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142661495","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":"Determination of the modified exterior Steklov eigenvalues via the reciprocity gap method","authors":"Wensong Qiu, Hongyan Wang, Yuan Li, Lixin Feng","doi":"10.1016/j.cam.2024.116360","DOIUrl":"10.1016/j.cam.2024.116360","url":null,"abstract":"<div><div>The modified exterior Steklov eigenvalues (MESEs) arise from the inverse scattering problem for inhomogeneous media with a cavity and may serve as potential target signatures in nondestructive testing. In this paper we are interested in the determination of the MESEs from the measured Cauchy data of the total field on some manifold inside the cavity due to interior point sources. To this end, the reciprocity gap (RG) method based on a linear integral equation is employed. We provide the related theory and show that the blow-up property of the approximate solution to the integral equation can be used to characterize the MESEs. Numerical examples are presented to demonstrate the viability of our method.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"459 ","pages":"Article 116360"},"PeriodicalIF":2.1,"publicationDate":"2024-11-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142661493","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}