发布求助
文献互助 智能选刊 最新文献

IMA Journal of Numerical Analysis最新文献

筛选
英文 中文
Inf-sup stable discretization of the quasi-static Biot’s equations in poroelasticity 孔隙弹性准静态Biot方程的非稳定离散化
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2025-07-07 DOI: 10.1093/imanum/draf032
Christian Kreuzer, Pietro Zanotti
{"title":"Inf-sup stable discretization of the quasi-static Biot’s equations in poroelasticity","authors":"Christian Kreuzer, Pietro Zanotti","doi":"10.1093/imanum/draf032","DOIUrl":"https://doi.org/10.1093/imanum/draf032","url":null,"abstract":"We propose a new full discretization of the Biot’s equations in poroelasticity. The construction is driven by the inf-sup theory, which we recently developed. It builds upon the four-field formulation of the equations obtained by introducing the total pressure and the total fluid content. We discretize in space with Lagrange finite elements and in time with backward Euler. We establish inf-sup stability and quasi-optimality of the proposed discretization, with robust constants with respect to all material parameters and the time horizon. We further construct an interpolant showing how the error decays for smooth solutions.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"189 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2025-07-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144577905","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
Finite element approximation of penalized elastoplastic torsion problem with nonconstant source term 非常源项惩罚弹塑性扭转问题的有限元逼近
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2025-07-03 DOI: 10.1093/imanum/draf052
Franz Chouly, Tom Gustafsson, Patrick Hild
{"title":"Finite element approximation of penalized elastoplastic torsion problem with nonconstant source term","authors":"Franz Chouly, Tom Gustafsson, Patrick Hild","doi":"10.1093/imanum/draf052","DOIUrl":"https://doi.org/10.1093/imanum/draf052","url":null,"abstract":"This study is concerned with the finite element approximation of the elastoplastic torsion problem. We focus on the case of a nonconstant source term, which cannot be easily recast into an obstacle problem as can be done in the case of a constant source term. We present a simple formulation that penalizes the constraint directly on the gradient norm of the solution. We study its well-posedness, derive error estimates and present numerical results to illustrate the theory.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"29 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2025-07-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144565936","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 Riemannian inexact Newton method for solving the orthogonal INDSCAL problem in multidimensional scaling 求解多维标度中正交INDSCAL问题的黎曼非精确牛顿法
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2025-07-03 DOI: 10.1093/imanum/draf047
Xue-lin Zhou, Chao-qian Li, Jiao-fen Li, Xue-feng Duan
{"title":"A Riemannian inexact Newton method for solving the orthogonal INDSCAL problem in multidimensional scaling","authors":"Xue-lin Zhou, Chao-qian Li, Jiao-fen Li, Xue-feng Duan","doi":"10.1093/imanum/draf047","DOIUrl":"https://doi.org/10.1093/imanum/draf047","url":null,"abstract":"The well-known individual differences scaling (INDSCAL) model is intended for simultaneous metric multidimensional scaling (MDS) of several doubly centered matrices of squared dissimilarities. In this work the problem of fitting the orthogonal INDSCAL model to the data is reformulated and studied as a matrix optimization problem on the product manifold of orthonormal and diagonal matrices. A Riemannian inexact Newton method is proposed to address the underlying problem, with the global and quadratic convergence of the proposed method established under some mild assumptions. Furthermore, the positive definiteness condition of the Riemannian Hessian of the objective function at a solution is derived. Some numerical experiments are provided to illustrate the efficiency of the proposed method.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"152 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2025-07-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144566072","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
Reflection coupling for unadjusted generalized Hamiltonian Monte Carlo in the nonconvex stochastic gradient case 非凸随机梯度情况下非调整广义哈密顿蒙特卡罗的反射耦合
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2025-06-29 DOI: 10.1093/imanum/draf045
Martin Chak, Pierre Monmarché
{"title":"Reflection coupling for unadjusted generalized Hamiltonian Monte Carlo in the nonconvex stochastic gradient case","authors":"Martin Chak, Pierre Monmarché","doi":"10.1093/imanum/draf045","DOIUrl":"https://doi.org/10.1093/imanum/draf045","url":null,"abstract":"Contraction in Wasserstein 1-distance with explicit rates is established for generalized Hamiltonian Monte Carlo with stochastic gradients under possibly nonconvex conditions. The algorithms considered include splitting schemes of kinetic Langevin diffusion, which are commonly used in molecular dynamics simulations. To accommodate the degenerate noise structure corresponding to inertia existing in the chain, a characteristically discrete-in-time coupling and contraction proof is devised. As consequence, quantitative Gaussian concentration bounds are provided for empirical averages. Convergence in Wasserstein 2-distance and total variation are also given, together with numerical bias estimates.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"32 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2025-06-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144513176","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
Stochastic theta methods for free stochastic differential equations 自由随机微分方程的随机θ方法
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2025-06-29 DOI: 10.1093/imanum/draf044
Yuanling Niu, Jiaxin Wei, Zhi Yin, Dan Zeng
{"title":"Stochastic theta methods for free stochastic differential equations","authors":"Yuanling Niu, Jiaxin Wei, Zhi Yin, Dan Zeng","doi":"10.1093/imanum/draf044","DOIUrl":"https://doi.org/10.1093/imanum/draf044","url":null,"abstract":"We introduce free probability analogues of the stochastic theta methods for free stochastic differential equations in this work. Assuming that the drift coefficient of the free stochastic differential equations is operator Lipschitz and the diffusion coefficients are locally operator Lipschitz we prove the strong convergence of the numerical methods. Moreover, we investigate the exponential stability in mean square of the equations and the numerical methods. In particular, the free stochastic theta methods with $theta in [1/2, 1]$ can inherit the exponential stability of original equations for any given step size. Our methods offer better stability than the free Euler–Maruyama method. Numerical results are reported to confirm these theoretical findings and show the efficiency of our methods compared with the free Euler–Maruyama method.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"18 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2025-06-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144513174","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
Adaptive approximations of inclusions in a semilinear elliptic problem related to cardiac electrophysiology 与心脏电生理相关的半线性椭圆问题中内含物的自适应逼近
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2025-06-28 DOI: 10.1093/imanum/draf041
Bangti Jin, Fengru Wang, Yifeng Xu
{"title":"Adaptive approximations of inclusions in a semilinear elliptic problem related to cardiac electrophysiology","authors":"Bangti Jin, Fengru Wang, Yifeng Xu","doi":"10.1093/imanum/draf041","DOIUrl":"https://doi.org/10.1093/imanum/draf041","url":null,"abstract":"In this work, we investigate the numerical reconstruction of inclusions in a semilinear elliptic equation arising in the mathematical modeling of cardiac ischemia. We propose an adaptive finite-element method for the resulting constrained minimization problem that is relaxed by a phase-field approach. The a posteriori error estimators of the adaptive algorithm consist of three components, i.e., the state variable, the adjoint variable and the complementary relation. Moreover, using tools from adaptive finite-element analysis and nonlinear optimization, we establish the strong convergence for a subsequence of adaptively generated discrete solutions to a solution of the continuous optimality system. Several numerical examples are presented to illustrate the convergence and efficiency of the adaptive algorithm.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"87 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2025-06-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144503630","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
Uniformly higher order accurate schemes for dynamics of charged particles under fast oscillating magnetic fields 快速振荡磁场下带电粒子动力学的均匀高阶精确格式
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2025-06-28 DOI: 10.1093/imanum/draf048
Megala Anandan, Benjamin Boutin, Nicolas Crouseilles
{"title":"Uniformly higher order accurate schemes for dynamics of charged particles under fast oscillating magnetic fields","authors":"Megala Anandan, Benjamin Boutin, Nicolas Crouseilles","doi":"10.1093/imanum/draf048","DOIUrl":"https://doi.org/10.1093/imanum/draf048","url":null,"abstract":"This work deals with the numerical approximation of plasmas that are confined by the effect of a fast oscillating magnetic field [Bostan, M. (2012), Transport of charged particles under fast oscillating magnetic fields. SIAM J. Math. Anal., 44, 1415–1447] in the Vlasov model. The presence of this magnetic field induces oscillations (in time) to the solution of the characteristic equations. Due to its multiscale character, a standard time discretization would lead to an inefficient solver. In this work, time integrators are derived and analyzed for a class of highly oscillatory differential systems. We prove the uniform accuracy property of these time integrators, meaning that the accuracy does not depend on the small parameter $varepsilon $. Moreover, we construct an extension of the scheme, which degenerates towards an energy preserving numerical scheme for the averaged model, when $varepsilon to 0$. Several numerical results illustrate the capabilities of the method.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"87 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2025-06-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144503359","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
Riemannian optimization methods for ground states of multicomponent Bose–Einstein condensates 多组分玻色-爱因斯坦凝聚基态的黎曼优化方法
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2025-06-28 DOI: 10.1093/imanum/draf046
Robert Altmann, Martin Hermann, Daniel Peterseim, Tatjana Stykel
{"title":"Riemannian optimization methods for ground states of multicomponent Bose–Einstein condensates","authors":"Robert Altmann, Martin Hermann, Daniel Peterseim, Tatjana Stykel","doi":"10.1093/imanum/draf046","DOIUrl":"https://doi.org/10.1093/imanum/draf046","url":null,"abstract":"This paper addresses the computation of ground states of multicomponent Bose–Einstein condensates, defined as the global minimizer of an energy functional on an infinite-dimensional generalised oblique manifold. We establish the existence of the ground state, prove its uniqueness up to scaling and characterize it as the solution to a coupled nonlinear eigenvector problem. By equipping the manifold with several Riemannian metrics, we introduce a suite of Riemannian gradient descent and Riemannian Newton methods. Metrics that incorporate first- or second-order information about the energy are particularly advantageous, effectively preconditioning the resulting methods. For a Riemannian gradient descent method with an energy-adaptive metric, we provide a qualitative global and quantitative local convergence analysis, confirming its reliability and robustness with respect to the choice of the spatial discretization. Numerical experiments highlight the computational efficiency of both the Riemannian gradient descent and Newton methods.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"46 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2025-06-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144503629","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
Using generalized simplex methods to approximate derivatives 用广义单纯形法逼近导数
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2025-06-28 DOI: 10.1093/imanum/draf053
Gabriel Jarry-Bolduc, Chayne Planiden
{"title":"Using generalized simplex methods to approximate derivatives","authors":"Gabriel Jarry-Bolduc, Chayne Planiden","doi":"10.1093/imanum/draf053","DOIUrl":"https://doi.org/10.1093/imanum/draf053","url":null,"abstract":"This paper presents two methods for approximating a proper subset of the entries of a Hessian using only function evaluations. It is also shown how to approximate a Hessian-vector product with a minimal number of function evaluations. These approximations are obtained using the techniques called generalized simplex Hessian and generalized centred simplex Hessian. We show how to choose the matrices of directions involved in the computation of these two techniques, depending on the entries of the Hessian of interest. We discuss the number of function evaluations required in each case and develop a general formula to approximate all order-$P$ partial derivatives. Since only function evaluations are required to compute the methods discussed in this paper they are suitable for use in derivative-free optimization methods.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"643 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2025-06-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144503628","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
Exponential convergence of hp-ILGFEM for semilinear elliptic boundary value problems with monomial reaction 半线性椭圆型单反应边值问题hp-ILGFEM的指数收敛性
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2025-06-26 DOI: 10.1093/imanum/draf030
Yanchen He, Paul Houston, Christoph Schwab, Thomas P Wihler
{"title":"Exponential convergence of hp-ILGFEM for semilinear elliptic boundary value problems with monomial reaction","authors":"Yanchen He, Paul Houston, Christoph Schwab, Thomas P Wihler","doi":"10.1093/imanum/draf030","DOIUrl":"https://doi.org/10.1093/imanum/draf030","url":null,"abstract":"We study the fully explicit numerical approximation of a semilinear elliptic boundary value model problem, which features a monomial reaction and analytic forcing, in a bounded polygon $varOmega subset{mathbb{R}}^{2}$ with a finite number of straight edges. In particular, we analyse the convergence of $hp$-type iterative linearized Galerkin ($hp$-ILG) solvers. Our convergence analysis is carried out for conforming $hp$-finite element (FE) Galerkin discretizations on sequences of regular, simplicial partitions of $varOmega $, with geometric corner refinement, with polynomial degrees increasing in tandem with the geometric mesh refinement towards the corners of $varOmega $. For a sequence of discrete solutions generated by the ILG solver with a stopping criterion that is consistent with the exponential convergence of the exact $hp$-FE Galerkin solution we prove exponential convergence in $text{H}^{1}(varOmega )$ to the unique weak solution of the boundary value problem. Numerical experiments illustrate the exponential convergence of the numerical approximations obtained from the proposed scheme in terms of the number of degrees of freedom as well as of the computational complexity involved.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"26 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2025-06-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144500473","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学术文献互助群
群 号：604180095
Book学术官方微信