IMA Journal of Numerical Analysis最新文献

筛选
英文 中文
Linearly implicit energy-preserving integrating factor methods and convergence analysis for the 2D nonlinear Schrödinger equation with wave operator 带波算子的二维非线性Schrödinger方程的线性隐式保能积分因子方法及其收敛性分析
2区 数学
IMA Journal of Numerical Analysis Pub Date : 2023-09-25 DOI: 10.1093/imanum/drad067
Xuelong Gu, Wenjun Cai, Yushun Wang, Chaolong Jiang
{"title":"Linearly implicit energy-preserving integrating factor methods and convergence analysis for the 2D nonlinear Schrödinger equation with wave operator","authors":"Xuelong Gu, Wenjun Cai, Yushun Wang, Chaolong Jiang","doi":"10.1093/imanum/drad067","DOIUrl":"https://doi.org/10.1093/imanum/drad067","url":null,"abstract":"Abstract In this paper, we develop a novel class of linearly implicit and energy-preserving integrating factor methods for the 2D nonlinear Schrödinger equation with wave operator (NLSW), combining the scalar auxiliary variable approach and the integrating factor methods. To begin, a second-order scheme is proposed, which is rigorously proved to be energy-preserving. By using the energy methods, we analyze its optimal convergence without any restrictions on the grid ratio, where a novel technique and an improved induction argument are proposed to circumvent the difficulty arising from the unavailability of a priori$L^{infty }$ estimates of numerical solutions. Based on the integrating factor Runge–Kutta methods, we extend the proposed scheme to arbitrarily high order, which is also linearly implicit and conservative. Numerical experiments are presented to confirm the theoretical analysis and demonstrate the advantages of the proposed methods.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"20 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-09-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135859525","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}
引用次数: 0
A general collocation analysis for weakly singular Volterra integral equations with variable exponent 变指数弱奇异Volterra积分方程的一般配位分析
2区 数学
IMA Journal of Numerical Analysis Pub Date : 2023-09-19 DOI: 10.1093/imanum/drad072
Hui Liang, Martin Stynes
{"title":"A general collocation analysis for weakly singular Volterra integral equations with variable exponent","authors":"Hui Liang, Martin Stynes","doi":"10.1093/imanum/drad072","DOIUrl":"https://doi.org/10.1093/imanum/drad072","url":null,"abstract":"Abstract Piecewise polynomial collocation of weakly singular Volterra integral equations (VIEs) of the second kind has been extensively studied in the literature, where integral kernels of the form $(t-s)^{-alpha }$ for some constant $alpha in (0,1)$ are considered. Variable-order fractional-derivative differential equations currently attract much research interest, and in Zheng & Wang (2020, An optimal-order numerical approximation to variable-order space-fractional diffusion equations on uniform or graded meshes. SIAM J. Numer. Anal., 58, 330–352), such a problem is transformed to a weakly singular VIE whose kernel has the above form with variable $alpha = alpha (t)$, then solved numerically by piecewise linear collocation, but it is unclear whether this analysis could be extended to more general problems or to polynomials of higher degree. In the present paper, the general theory (existence, uniqueness, regularity of solutions) of variable-exponent weakly singular VIEs is developed using novel techniques. These results then underpin an error analysis of collocation methods where piecewise polynomials of any degree can be used. This error analysis is also novel—it makes no use of the usual resolvent representation, which is a key technique in the error analysis of collocation methods for VIEs in the current research literature. Furthermore, all the above analysis for a scalar VIE can be extended to certain nonlinear VIEs and to systems of VIEs. The sharpness of the theoretical error bounds obtained for the collocation methods is demonstrated by numerical examples.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"111 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135108645","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}
引用次数: 0
Hybrid neural-network FEM approximation of diffusion coefficient in elliptic and parabolic Problems 椭圆型和抛物型问题扩散系数的混合神经网络有限元逼近
2区 数学
IMA Journal of Numerical Analysis Pub Date : 2023-09-18 DOI: 10.1093/imanum/drad073
Siyu Cen, Bangti Jin, Qimeng Quan, Zhi Zhou
{"title":"Hybrid neural-network FEM approximation of diffusion coefficient in elliptic and parabolic Problems","authors":"Siyu Cen, Bangti Jin, Qimeng Quan, Zhi Zhou","doi":"10.1093/imanum/drad073","DOIUrl":"https://doi.org/10.1093/imanum/drad073","url":null,"abstract":"Abstract In this work we investigate the numerical identification of the diffusion coefficient in elliptic and parabolic problems using neural networks (NNs). The numerical scheme is based on the standard output least-squares formulation where the Galerkin finite element method (FEM) is employed to approximate the state and NNs act as a smoothness prior to approximate the unknown diffusion coefficient. A projection operation is applied to the NN approximation in order to preserve the physical box constraint on the unknown coefficient. The hybrid approach enjoys both rigorous mathematical foundation of the FEM and inductive bias/approximation properties of NNs. We derive a priori error estimates in the standard $L^2(varOmega )$ norm for the numerical reconstruction, under a positivity condition which can be verified for a large class of problem data. The error bounds depend explicitly on the noise level, regularization parameter and discretization parameters (e.g., spatial mesh size, time step size and depth, upper bound and number of nonzero parameters of NNs). We also provide extensive numerical experiments, indicating that the hybrid method is very robust for large noise when compared with the pure FEM approximation.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"19 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135208123","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}
引用次数: 0
A posteriori error estimates for fully coupled McKean–Vlasov forward-backward SDEs 完全耦合McKean-Vlasov正向向后SDEs的后验误差估计
2区 数学
IMA Journal of Numerical Analysis Pub Date : 2023-09-15 DOI: 10.1093/imanum/drad060
Christoph Reisinger, Wolfgang Stockinger, Yufei Zhang
{"title":"<i>A posteriori</i> error estimates for fully coupled McKean–Vlasov forward-backward SDEs","authors":"Christoph Reisinger, Wolfgang Stockinger, Yufei Zhang","doi":"10.1093/imanum/drad060","DOIUrl":"https://doi.org/10.1093/imanum/drad060","url":null,"abstract":"Abstract Fully coupled McKean–Vlasov forward-backward stochastic differential equations (MV-FBSDEs) arise naturally from large population optimization problems. Judging the quality of given numerical solutions for MV-FBSDEs, which usually require Picard iterations and approximations of nested conditional expectations, is typically difficult. This paper proposes an a posteriori error estimator to quantify the $L^2$-approximation error of an arbitrarily generated approximation on a time grid. We establish that the error estimator is equivalent to the global approximation error between the given numerical solution and the solution of a forward Euler discretized MV-FBSDE. A crucial and challenging step in the analysis is the proof of stability of this Euler approximation to the MV-FBSDE, which is of independent interest. We further demonstrate that, for sufficiently fine time grids, the accuracy of numerical solutions for solving the continuous MV-FBSDE can also be measured by the error estimator. The error estimates justify the use of residual-based algorithms for solving MV-FBSDEs. Numerical experiments for MV-FBSDEs arising from mean field control and games confirm the effectiveness and practical applicability of the error estimator.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"41 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-09-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135438201","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}
引用次数: 0
Group sparse optimization for inpainting of random fields on the sphere 球上随机场绘制的群稀疏优化
2区 数学
IMA Journal of Numerical Analysis Pub Date : 2023-09-15 DOI: 10.1093/imanum/drad071
Chao Li, Xiaojun Chen
{"title":"Group sparse optimization for inpainting of random fields on the sphere","authors":"Chao Li, Xiaojun Chen","doi":"10.1093/imanum/drad071","DOIUrl":"https://doi.org/10.1093/imanum/drad071","url":null,"abstract":"Abstract We propose a group sparse optimization model for inpainting of a square-integrable isotropic random field on the unit sphere, where the field is represented by spherical harmonics with random complex coefficients. In the proposed optimization model, the variable is an infinite-dimensional complex vector and the objective function is a real-valued function defined by a hybrid of the $ell _2$ norm and non-Lipschitz $ell _p (0&amp;lt;p&amp;lt;1)$ norm that preserves rotational invariance property and group structure of the random complex coefficients. We show that the infinite-dimensional optimization problem is equivalent to a convexly-constrained finite-dimensional optimization problem. Moreover, we propose a smoothing penalty algorithm to solve the finite-dimensional problem via unconstrained optimization problems. We provide an approximation error bound of the inpainted random field defined by a scaled Karush–Kuhn–Tucker (KKT) point of the constrained optimization problem in the square-integrable space on the sphere with probability measure. Finally, we conduct numerical experiments on band-limited random fields on the sphere and images from Cosmic Microwave Background (CMB) data to show the promising performance of the smoothing penalty algorithm for inpainting of random fields on the sphere.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"205 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-09-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135396059","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}
引用次数: 0
Maximum norm error bounds for the full discretization of nonautonomous wave equations 非自治波动方程完全离散化的最大范数误差界
2区 数学
IMA Journal of Numerical Analysis Pub Date : 2023-09-09 DOI: 10.1093/imanum/drad065
Benjamin Dörich, Jan Leibold, Bernhard Maier
{"title":"Maximum norm error bounds for the full discretization of nonautonomous wave equations","authors":"Benjamin Dörich, Jan Leibold, Bernhard Maier","doi":"10.1093/imanum/drad065","DOIUrl":"https://doi.org/10.1093/imanum/drad065","url":null,"abstract":"Abstract In the present paper, we consider a specific class of nonautonomous wave equations on a smooth, bounded domain and their discretization in space by isoparametric finite elements and in time by the implicit Euler method. Building upon the work of Baker and Dougalis (1980, On the ${L}^{infty }$-convergence of Galerkin approximations for second-order hyperbolic equations. Math. Comp., 34, 401–424), we prove optimal error bounds in the $W^{1,infty } times L^infty $-norm for the semidiscretization in space and the full discretization. The key tool is the gain of integrability coming from the inverse of the discretized differential operator. For this, we have to pay with (discrete) time derivatives on the error in the $H^{1} times L^2$-norm, which are reduced to estimates of the differentiated initial errors. To confirm our theoretical findings, we also present numerical experiments.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"40 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"136192193","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}
引用次数: 0
Local absorbing boundary conditions on fixed domains give order-one errors for high-frequency waves 固定域上的局部吸收边界条件对高频波产生一阶误差
2区 数学
IMA Journal of Numerical Analysis Pub Date : 2023-09-09 DOI: 10.1093/imanum/drad058
Jeffrey Galkowski, David Lafontaine, Euan A. Spence
{"title":"Local absorbing boundary conditions on fixed domains give order-one errors for high-frequency waves","authors":"Jeffrey Galkowski, David Lafontaine, Euan A. Spence","doi":"10.1093/imanum/drad058","DOIUrl":"https://doi.org/10.1093/imanum/drad058","url":null,"abstract":"Abstract We consider approximating the solution of the Helmholtz exterior Dirichlet problem for a nontrapping obstacle, with boundary data coming from plane-wave incidence, by the solution of the corresponding boundary value problem where the exterior domain is truncated and a local absorbing boundary condition coming from a Padé approximation (of arbitrary order) of the Dirichlet-to-Neumann map is imposed on the artificial boundary (recall that the simplest such boundary condition is the impedance boundary condition). We prove upper- and lower-bounds on the relative error incurred by this approximation, both in the whole domain and in a fixed neighbourhood of the obstacle (i.e., away from the artificial boundary). Our bounds are valid for arbitrarily-high frequency, with the artificial boundary fixed, and show that the relative error is bounded away from zero, independent of the frequency, and regardless of the geometry of the artificial boundary.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"54 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"136107442","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}
引用次数: 3
Error estimates of high-order compact finite difference schemes for the nonlinear abcd Boussinesq systems 非线性abcd Boussinesq系统高阶紧致有限差分格式的误差估计
2区 数学
IMA Journal of Numerical Analysis Pub Date : 2023-09-08 DOI: 10.1093/imanum/drad069
Su-Cheol Yi, Kai Fu, Shusen Xie
{"title":"Error estimates of high-order compact finite difference schemes for the nonlinear <i>abcd</i> Boussinesq systems","authors":"Su-Cheol Yi, Kai Fu, Shusen Xie","doi":"10.1093/imanum/drad069","DOIUrl":"https://doi.org/10.1093/imanum/drad069","url":null,"abstract":"Abstract In this paper, some fourth-order compact finite difference schemes are derived and analyzed for the nonlinear $abcd$ Boussinesq systems. The optimal order error estimates for the semidiscrete compact finite difference schemes with different cases of dispersion coefficients $a, b, c, d$, are presented. The third-order and fourth-order linearized implicit multistep schemes are adopted for time discretization, and numerical experiments are conducted on the model problems. Numerical results show that the proposed schemes have high accuracy and are consistent with the theoretical analysis.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"40 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-09-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"136299906","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}
引用次数: 0
Discrete Gagliardo–Nirenberg inequality and application to the finite volume approximation of a convection–diffusion equation with a Joule effect term 离散gagliado - nirenberg不等式及其在具有焦耳效应项的对流扩散方程有限体积近似中的应用
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2023-08-30 DOI: 10.1093/imanum/drad063
C. Calgaro, C. Cancès, E. Creusé
{"title":"Discrete Gagliardo–Nirenberg inequality and application to the finite volume approximation of a convection–diffusion equation with a Joule effect term","authors":"C. Calgaro, C. Cancès, E. Creusé","doi":"10.1093/imanum/drad063","DOIUrl":"https://doi.org/10.1093/imanum/drad063","url":null,"abstract":"A discrete order-two Gagliardo–Nirenberg inequality is established for piecewise constant functions defined on a two-dimensional structured mesh composed of rectangular cells. As in the continuous framework, this discrete Gagliardo–Nirenberg inequality allows to control in particular the $L^4$ norm of the discrete gradient of the numerical solution by the $L^2$ norm of its discrete Hessian times its $L^infty $ norm. This result is crucial for the convergence analysis of a finite volume method for the approximation of a convection–diffusion equation involving a Joule effect term on a uniform mesh in each direction. The convergence proof relies on compactness arguments and on a priori estimates under a smallness assumption on the data, which is essential also in the continuous framework.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":" ","pages":""},"PeriodicalIF":2.1,"publicationDate":"2023-08-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48381868","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}
引用次数: 0
An arbitrary-order discrete rot-rot complex on polygonal meshes with application to a quad-rot problem 多边形网格上任意阶离散rot-rot复形及其在四次rot问题中的应用
2区 数学
IMA Journal of Numerical Analysis Pub Date : 2023-08-28 DOI: 10.1093/imanum/drad045
Daniele Antonio Di Pietro
{"title":"An arbitrary-order discrete rot-rot complex on polygonal meshes with application to a quad-rot problem","authors":"Daniele Antonio Di Pietro","doi":"10.1093/imanum/drad045","DOIUrl":"https://doi.org/10.1093/imanum/drad045","url":null,"abstract":"Abstract In this work, following the discrete de Rham approach, we develop a discrete counterpart of a two-dimensional de Rham complex with enhanced regularity. The proposed construction supports general polygonal meshes and arbitrary approximation orders. We establish exactness on a contractible domain for both the versions of the complex with and without boundary conditions and, for the former, prove a complete set of Poincaré-type inequalities. The discrete complex is then used to derive a novel discretization method for a quad-rot problem, which, unlike other schemes in the literature, does not require the forcing term to be prepared. We carry out complete stability and convergence analyses for the proposed scheme and provide numerical validation of the results.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"42 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-08-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"136348995","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}
引用次数: 0
0
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
确定
请完成安全验证×
相关产品
×
本文献相关产品
联系我们:info@booksci.cn Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。 Copyright © 2023 布克学术 All rights reserved.
京ICP备2023020795号-1
ghs 京公网安备 11010802042870号
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术官方微信