{"title":"Error estimates of exponential wave integrators for the Dirac equation in the massless and nonrelativistic regime","authors":"Ying Ma, Lizhen Chen","doi":"10.1002/num.23128","DOIUrl":"https://doi.org/10.1002/num.23128","url":null,"abstract":"We present exponential wave integrator Fourier pseudospectral (EWI‐FP) methods and establish their error estimates of the fully discrete schemes for the Dirac equation in the massless and nonrelativistic regime. This regime involves a small dimensionless parameter where , and is inversely proportional to the speed of light. The solution exhibits highly oscillatory behavior in time and rapid wave propagation in space in this regime. Specifically, the time oscillations have a wavelength of , while the spatial oscillations have a wavelength of , with a wave speed of . We employ (symmetric) exponential wave integrators for temporal derivatives and Fourier spectral discretization for spatial derivatives. We rigorously derive the error bounds which explicitly depend on the mesh size , the time step and the small dimensionless parameter . The error estimates for the EWI‐FP methods demonstrate that their meshing strategy requirement (‐scalability) necessitates setting and when . Finally, some numerical examples are provided to validate the error bounds.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-08-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141925068","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Determining a time‐varying potential in time‐fractional diffusion from observation at a single point","authors":"Siyu Cen, Kwancheol Shin, Zhi Zhou","doi":"10.1002/num.23136","DOIUrl":"https://doi.org/10.1002/num.23136","url":null,"abstract":"We discuss the identification of a time‐dependent potential in a time‐fractional diffusion model from a boundary measurement taken at a single point. Theoretically, we establish a conditional Lipschitz stability for this inverse problem. Numerically, we develop an easily implementable iterative algorithm to recover the unknown coefficient, and also derive rigorous error bounds for the discrete reconstruction. These results are attained by leveraging the (discrete) solution theory of direct problems, and applying error estimates that are optimal with respect to problem data regularity. Numerical simulations are provided to demonstrate the theoretical results.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-08-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141930782","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"On the optimality of voronoi‐based column selection","authors":"Maria Emelianenko, Guy B. Oldaker","doi":"10.1002/num.23137","DOIUrl":"https://doi.org/10.1002/num.23137","url":null,"abstract":"While there exists a rich array of matrix column subset selection problem (CSSP) algorithms for use with interpolative and CUR‐type decompositions, their use can often become prohibitive as the size of the input matrix increases. In an effort to address these issues, in earlier work we developed a general framework that pairs a column‐partitioning routine with a column‐selection algorithm. Two of the four algorithms presented in that work paired the Centroidal Voronoi Orthogonal Decomposition (<jats:styled-content>CVOD</jats:styled-content>; <jats:italic>Centroidal Voronoi Tessellation Based Proper Orthogonal Decomposition Analysis</jats:italic>, 2003, 137–150) and an adaptive variant (<jats:styled-content>adaptCVOD</jats:styled-content>) with the Discrete Empirical Interpolation Method (<jats:styled-content>DEIM; <jats:italic>SIAM J. Sci. Computer</jats:italic>. 38 (2016), no. 3, A1454–A1482</jats:styled-content>). In this work, we extend this framework and pair the <jats:styled-content>CVOD</jats:styled-content>‐type algorithms with any CSSP algorithm that returns linearly independent columns. Our results include detailed error bounds for the solutions provided by these paired algorithms, as well as expressions that explicitly characterize how the quality of the selected column partition affects the resulting CSSP solution. In addition to examples involving matrix approximation, we test several of our partition‐based constructions on tasks commonly encountered in model order reduction (MOR).","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-08-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141930772","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Semiexplicit K‐symplectic‐like methods with energy conservation for noncanonical Hamiltonian systems","authors":"Beibei Zhu, Ran Gu","doi":"10.1002/num.23138","DOIUrl":"https://doi.org/10.1002/num.23138","url":null,"abstract":"For the nonseparable noncanonical Hamiltonian systems, we propose efficient K‐symplectic‐like methods which are semiexplicit and energy‐preserving. By introducing two copies of the phase space and constructing an augmented Hamiltonian, we can separate the noncanonical Hamiltonian system into two explicitly integrable parts. Subsequently, explicit K‐symplectic methods can be constructed by using the splitting and composing method. To enforce constraints on the two copies of the phase space, we provide two transformations with energy conservation property. This enables us to obtain semiexplicit K‐symplectic‐like methods that preserve energy. Two algorithms are provided to implement the semiexplicit K‐symplectic‐like methods with energy conservation and their convergence has been proved. Numerical results on two noncanonical Hamiltonian systems demonstrate that the energy errors of our proposed methods remain bounded within machine precision over long time without exhibiting energy drift. Furthermore, the proposed methods exhibit superior computational efficiency compared to the canonicalized symplectic methods of the same order.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-08-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141930774","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Ernesto Pimentel‐García, Manuel J. Castro, Christophe Chalons, Carlos Parés
{"title":"High‐order in‐cell discontinuous reconstruction path‐conservative methods for nonconservative hyperbolic systems–DR.MOOD method","authors":"Ernesto Pimentel‐García, Manuel J. Castro, Christophe Chalons, Carlos Parés","doi":"10.1002/num.23133","DOIUrl":"https://doi.org/10.1002/num.23133","url":null,"abstract":"In this work, we develop a new framework to deal numerically with discontinuous solutions in nonconservative hyperbolic systems. First an extension of the MOOD methodology to nonconservative systems based on Taylor expansions is presented. This extension combined with an in‐cell discontinuous reconstruction operator are the key points to develop a new family of high‐order methods that are able to capture exactly isolated shocks. Several test cases are proposed to validate these methods for the Modified Shallow Water equations and the Two‐Layer Shallow Water system.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-08-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141884916","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Uniform and optimal error estimates of a nested Picard integrator for the nonlinear Schrödinger equation with wave operator","authors":"Yongyong Cai, Yue Feng, Yichen Guo, Zhiguo Xu","doi":"10.1002/num.23135","DOIUrl":"https://doi.org/10.1002/num.23135","url":null,"abstract":"We propose a second‐order nested Picard iterative integrator sine pseudospectral (NPI‐SP) method for the nonlinear Schrödinger equation with wave operator (NLSW) involving a parameter and carry out rigorous error estimates. Actually, the equation propagates wave with wavelength in time, while the amplitude of the leading oscillation is for well‐prepared initial data, and for ill‐prepared initial data, respectively. Based on the exponential integrator and nested Picard iteration, the uniformly accurate (w.r.t. ) NPI‐SP scheme is proposed with the optimal uniform error bounds at in time and spectral accuracy in space for both well‐prepared and ill‐prepared data in ‐norm. This result significantly improves the error bounds of the finite difference methods and exponential wave integrator for the NLSW. Error estimates are rigorously carried out and numerical examples are provided to confirm the theoretical analysis.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-08-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141884757","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Partitioning method for the finite element approximation of a 3D fluid‐2D plate interaction system","authors":"Pelin G. Geredeli, Hemanta Kunwar, Hyesuk Lee","doi":"10.1002/num.23132","DOIUrl":"https://doi.org/10.1002/num.23132","url":null,"abstract":"We consider the finite element approximation of a coupled fluid‐structure interaction (FSI) system, which comprises a three‐dimensional (3D) Stokes flow and a two‐dimensional (2D) fourth‐order Euler–Bernoulli or Kirchhoff plate. The interaction of these parabolic and hyperbolic partial differential equations (PDE) occurs at the boundary interface which is assumed to be fixed. The vertical displacement of the plate dynamics evolves on the flat portion of the boundary where the coupling conditions are implemented via the matching velocities of the plate and fluid flow, as well as the Dirichlet boundary trace of the pressure. This pressure term also acts as a coupling agent, since it appears as a forcing term on the flat, elastic plate domain. Our main focus in this work is to generate some numerical results concerning the approximate solutions to the FSI model. For this, we propose a numerical algorithm that sequentially solves the fluid and plate subsystems through an effective decoupling approach. Numerical results of test problems are presented to illustrate the performance of the proposed method.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-08-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141884918","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A type of multigrid method for semilinear elliptic problems based on symmetric interior penalty discontinuous Galerkin method","authors":"Fan Chen, Ming Cui, Chenguang Zhou","doi":"10.1002/num.23130","DOIUrl":"https://doi.org/10.1002/num.23130","url":null,"abstract":"This article introduces a new kind of multigrid approach for semilinear elliptic problems, which is based on the symmetric interior penalty discontinuous Galerkin (SIPDG) method. We first give an optimal error estimate of the SIPDG method for the problem. Then, we design a type of multigrid method, which is called the multilevel correction method, and derive a priori error estimates. The primary idea of this method is to take the solution of the semilinear problem and utilize it to establish a sequence of solutions for associated linear boundary value problem on discontinuous finite element spaces and a newly defined low dimensional augmented subspace. Lastly, numerical experiments are offered to confirm the suggested method's precision and effectiveness.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-08-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141884835","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Optimal control of variably distributed‐order time‐fractional diffusion equation: Analysis and computation","authors":"Xiangcheng Zheng, Huan Liu, Hong Wang, Xu Guo","doi":"10.1002/num.23134","DOIUrl":"https://doi.org/10.1002/num.23134","url":null,"abstract":"Fractional diffusion equations exhibit competitive capabilities in modeling many challenging phenomena such as the anomalously diffusive transport and memory effects. We prove the well‐posedness and regularity of an optimal control of a variably distributed‐order fractional diffusion equation with pointwise constraints, where the distributed‐order operator accounts for, for example, the effect of uncertainties. We accordingly develop and analyze a fully‐discretized finite element approximation to the optimal control without any artificial regularity assumption of the true solution. Numerical experiments are also performed to substantiate the theoretical findings.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-08-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141884917","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Numerical solution of convected wave equation in free field using artificial boundary method","authors":"Xin Wang, Jihong Wang, Yana Di, Jiwei Zhang","doi":"10.1002/num.23131","DOIUrl":"https://doi.org/10.1002/num.23131","url":null,"abstract":"In this article, we propose two procedures focusing on the computation of the time‐dependent convected wave equation in a free field with a uniform background flow. Both procedures are based on a framework, expended from Du et al. (SIAM J. Sci. Comput. 40 (2018), A1430–A1445.), of constructing the Dirichlet‐to‐Dirichlet (DtD)‐type discrete absorbing boundary conditions (ABCs). The first procedure is dedicated to reducing the infinite problem into a finite problem by a direct application of the framework on the finite difference discretization of the convected wave equation. However, the presence of convection terms makes the stability analysis hard to implement, which motivates us to develop the second procedure. First, the convected wave equation is transformed into a standard wave equation by using the Prandtl‐Glauert‐Lorentz transformation. After that, we obtain the DtD‐type ABC by using the above framework, and on this basis, derive an equivalent Dirichlet‐to‐Neumann‐type ABCs, which makes stability and convergence analysis easy to be obtained by the classical energy method. The effectiveness and comparison of these two procedures are investigated through numerical experiments.","PeriodicalId":19443,"journal":{"name":"Numerical Methods for Partial Differential Equations","volume":null,"pages":null},"PeriodicalIF":3.9,"publicationDate":"2024-07-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141863952","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}