IMA Journal of Numerical Analysis最新文献

筛选
英文 中文
Gamma-convergent LDG method for large bending deformations of bilayer plates 双层板大弯曲变形的伽马收敛 LDG 方法
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2024-01-19 DOI: 10.1093/imanum/drad100
Andrea Bonito, Ricardo H Nochetto, Shuo Yang
{"title":"Gamma-convergent LDG method for large bending deformations of bilayer plates","authors":"Andrea Bonito, Ricardo H Nochetto, Shuo Yang","doi":"10.1093/imanum/drad100","DOIUrl":"https://doi.org/10.1093/imanum/drad100","url":null,"abstract":"Bilayer plates are slender structures made of two thin layers of different materials. They react to environmental stimuli and undergo large bending deformations with relatively small actuation. The reduced model is a constrained minimization problem for the second fundamental form, with a given spontaneous curvature that encodes material properties, subject to an isometry constraint. We design a local discontinuous Galerkin (LDG) method, which imposes a relaxed discrete isometry constraint and controls deformation gradients at barycenters of elements. We prove $varGamma $-convergence of LDG, design a fully practical gradient flow, which gives rise to a linear scheme at every step, and show energy stability and control of the isometry defect. We extend the $varGamma $-convergence analysis to piecewise quadratic creases. We also illustrate the performance of the LDG method with several insightful simulations of large deformations, one including a curved crease.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-01-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139505921","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
Darcy’s problem coupled with the heat equation under singular forcing: analysis and discretization 奇异强迫下与热方程耦合的达西问题:分析与离散化
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2024-01-11 DOI: 10.1093/imanum/drad094
Alejandro Allendes, Gilberto Campaña, Francisco Fuica, Enrique Otárola
{"title":"Darcy’s problem coupled with the heat equation under singular forcing: analysis and discretization","authors":"Alejandro Allendes, Gilberto Campaña, Francisco Fuica, Enrique Otárola","doi":"10.1093/imanum/drad094","DOIUrl":"https://doi.org/10.1093/imanum/drad094","url":null,"abstract":"We study the existence of solutions for Darcy’s problem coupled with the heat equation under singular forcing; the right-hand side of the heat equation corresponds to a Dirac measure. The model studied involves thermal diffusion and viscosity depending on the temperature. We propose a finite element solution technique and analyze its convergence properties. In the case where thermal diffusion is independent of temperature, we propose an a posteriori error estimator and study its reliability and efficiency properties. We illustrate the theory with numerical examples.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-01-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139431258","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
Compactness estimates for difference schemes for conservation laws with discontinuous flux 非连续通量守恒定律差分方案的紧凑性估计
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2024-01-04 DOI: 10.1093/imanum/drad096
Kenneth H Karlsen, John D Towers
{"title":"Compactness estimates for difference schemes for conservation laws with discontinuous flux","authors":"Kenneth H Karlsen, John D Towers","doi":"10.1093/imanum/drad096","DOIUrl":"https://doi.org/10.1093/imanum/drad096","url":null,"abstract":"We establish quantitative compactness estimates for finite difference schemes used to solve nonlinear conservation laws. These equations involve a flux function $f(k(x,t),u)$, where the coefficient $k(x,t)$ is $BV$-regular and may exhibit discontinuities along curves in the $(x,t)$ plane. Our approach, which is technically elementary, relies on a discrete interaction estimate and one entropy function. While the details are specifically outlined for the Lax-Friedrichs scheme, the same framework can be applied to other difference schemes. Notably, our compactness estimates are new even in the homogeneous case ($kequiv 1$).","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-01-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139101268","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 constraint dissolving approach for nonsmooth optimization over the Stiefel manifold 在 Stiefel 流形上进行非平滑优化的约束消解方法
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2023-12-23 DOI: 10.1093/imanum/drad098
Xiaoyin Hu, Nachuan Xiao, Xin Liu, Kim-Chuan Toh
{"title":"A constraint dissolving approach for nonsmooth optimization over the Stiefel manifold","authors":"Xiaoyin Hu, Nachuan Xiao, Xin Liu, Kim-Chuan Toh","doi":"10.1093/imanum/drad098","DOIUrl":"https://doi.org/10.1093/imanum/drad098","url":null,"abstract":"This paper focuses on the minimization of a possibly nonsmooth objective function over the Stiefel manifold. The existing approaches either lack efficiency or can only tackle prox-friendly objective functions. We propose a constraint dissolving function named NCDF and show that it has the same first-order stationary points and local minimizers as the original problem in a neighborhood of the Stiefel manifold. Furthermore, we show that the Clarke subdifferential of NCDF is easy to achieve from the Clarke subdifferential of the objective function. Therefore, various existing approaches for unconstrained nonsmooth optimization can be directly applied to nonsmooth optimization problems over the Stiefel manifold. We propose a framework for developing subgradient-based methods and establishing their convergence properties based on prior works. Furthermore, based on our proposed framework, we can develop efficient approaches for optimization over the Stiefel manifold. Preliminary numerical experiments further highlight that the proposed constraint dissolving approach yields efficient and direct implementations of various unconstrained approaches to nonsmooth optimization problems over the Stiefel manifold.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2023-12-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139041468","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
On best p-norm approximation of discrete data by polynomials 论多项式对离散数据的最佳 p-norm 近似值
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2023-12-09 DOI: 10.1093/imanum/drad086
Michael S Floater
{"title":"On best p-norm approximation of discrete data by polynomials","authors":"Michael S Floater","doi":"10.1093/imanum/drad086","DOIUrl":"https://doi.org/10.1093/imanum/drad086","url":null,"abstract":"In this note, we derive a solution to the problem of finding a polynomial of degree at most $n$ that best approximates data at $n+2$ points in the $l_{p}$ norm. Analogous to a result of de la Vallée Poussin, one can express the solution as a convex combination of the Lagrange interpolants over subsets of $n+1$ points, and the error oscillates in sign.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2023-12-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138559326","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
Semilinear optimal control with Dirac measures 具有狄拉克测度的半线性最优控制
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2023-11-29 DOI: 10.1093/imanum/drad091
Enrique Otárola
{"title":"Semilinear optimal control with Dirac measures","authors":"Enrique Otárola","doi":"10.1093/imanum/drad091","DOIUrl":"https://doi.org/10.1093/imanum/drad091","url":null,"abstract":"The purpose of this work is to study an optimal control problem for a semilinear elliptic partial differential equation with a linear combination of Dirac measures as a forcing term; the control variable corresponds to the amplitude of such singular sources. We analyze the existence of optimal solutions and derive first- and, necessary and sufficient, second-order optimality conditions. We develop a solution technique that discretizes the state and adjoint equations with continuous piecewise linear finite elements; the control variable is already discrete. We analyze the convergence properties of discretizations and obtain, in two dimensions, an a priori error estimate for the underlying approximation of an optimal control variable.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2023-11-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138455384","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
Unconditionally stable small stencil enriched multiple point flux approximations of heterogeneous diffusion problems on general meshes 一般网格上非均匀扩散问题的无条件稳定小模板富多点通量近似
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2023-11-25 DOI: 10.1093/imanum/drad087
Julien Coatléven
{"title":"Unconditionally stable small stencil enriched multiple point flux approximations of heterogeneous diffusion problems on general meshes","authors":"Julien Coatléven","doi":"10.1093/imanum/drad087","DOIUrl":"https://doi.org/10.1093/imanum/drad087","url":null,"abstract":"We derive new multiple point flux approximations (MPFA) for the finite volume approximation of heterogeneous and anisotropic diffusion problems on general meshes, in dimensions 2 and 3. The resulting methods are unconditionally stable while preserving the small stencil typical of MPFA finite volumes. The key idea is to solve local variational problems with a well-designed stabilization term from which we deduce conservative flux instead of directly prescribing a flux formula and solving the usual flux continuity equations. The boundary conditions of our local variational problems are handled through additional cell-centered unknowns, leading to an overall scheme with the same number of unknowns than first-order discontinuous Galerkin methods. Convergence results follow from well-established frameworks, while numerical experiments illustrate the good behavior of the method.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2023-11-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138449679","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 VEM for variable data: convergence and optimality 可变数据的自适应VEM:收敛性和最优性
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2023-11-16 DOI: 10.1093/imanum/drad085
L Beirão da Veiga, C Canuto, R H Nochetto, G Vacca, M Verani
{"title":"Adaptive VEM for variable data: convergence and optimality","authors":"L Beirão da Veiga, C Canuto, R H Nochetto, G Vacca, M Verani","doi":"10.1093/imanum/drad085","DOIUrl":"https://doi.org/10.1093/imanum/drad085","url":null,"abstract":"We design an adaptive virtual element method (AVEM) of lowest order over triangular meshes with hanging nodes in 2d, which are treated as polygons. AVEM hinges on the stabilization-free a posteriori error estimators recently derived in Beirão da Veiga et al. (2023, Adaptive VEM: stabilization-free a posteriori error analysis and contraction property. SIAM J. Numer. Anal., 61, 457–494). The crucial property, which also plays a central role in this paper, is that the stabilization term can be made arbitrarily small relative to the a posteriori error estimators upon increasing the stabilization parameter. Our AVEM concatenates two modules, GALERKIN and DATA. The former deals with piecewise constant data and is shown in the above article to be a contraction between consecutive iterates. The latter approximates general data by piecewise constants to a desired accuracy. AVEM is shown to be convergent and quasi-optimal, in terms of error decay versus degrees of freedom, for solutions and data belonging to appropriate approximation classes. Numerical experiments illustrate the interplay between these two modules and provide computational evidence of optimality.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2023-11-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138293500","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 two-dimensional boundary value problem of elliptic type with nonlocal conjugation conditions 具有非局部共轭条件的椭圆型二维边值问题
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2023-11-14 DOI: 10.1093/imanum/drad084
Zorica Milovanović Jeknić, Aleksandra Delić, Sandra Živanović
{"title":"A two-dimensional boundary value problem of elliptic type with nonlocal conjugation conditions","authors":"Zorica Milovanović Jeknić, Aleksandra Delić, Sandra Živanović","doi":"10.1093/imanum/drad084","DOIUrl":"https://doi.org/10.1093/imanum/drad084","url":null,"abstract":"We consider an elliptic boundary value problem with nonlocal conjugation conditions. An a priori estimate for its weak solution in an appropriate Sobolev-like space is proved. A finite difference scheme approximating this problem is proposed and analyzed. An estimate of the convergence rate, compatible with the smoothness of the input data, is obtained.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2023-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"109126944","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
Convergent finite element methods for the perfect conductivity problem with close-to-touching inclusions 近接触夹杂物完美电导率问题的收敛有限元方法
IF 2.1 2区 数学
IMA Journal of Numerical Analysis Pub Date : 2023-11-14 DOI: 10.1093/imanum/drad088
Buyang Li, Haigang Li, Zongze Yang
{"title":"Convergent finite element methods for the perfect conductivity problem with close-to-touching inclusions","authors":"Buyang Li, Haigang Li, Zongze Yang","doi":"10.1093/imanum/drad088","DOIUrl":"https://doi.org/10.1093/imanum/drad088","url":null,"abstract":"In the perfect conductivity problem (i.e., the conductivity problem with perfectly conducting inclusions), the gradient of the electric field is often very large in a narrow region between two inclusions and blows up as the distance between the inclusions tends to zero. The rigorous error analysis for the computation of such perfect conductivity problems with close-to-touching inclusions of general geometry still remains open in three dimensions. We address this problem by establishing new asymptotic estimates for the second-order partial derivatives of the solution with explicit dependence on the distance $varepsilon $ between the inclusions, and use the asymptotic estimates to design a class of graded meshes and finite element spaces to solve the perfect conductivity problem with possibly close-to-touching inclusions. In particular, we propose a special finite element basis function that resolves the asymptotic singularity of the solution by making the interpolation error bounded in $W^{1,infty }$ in a neighborhood of the close-to-touching point, even though the solution itself is blowing up in $W^{1,infty }$. This is crucial in the error analysis for the numerical approximations. We prove that the proposed method yields optimal-order convergence in the $H^1$ norm, uniformly with respect to the distance $varepsilon $ between the inclusions, in both two and three dimensions for general convex smooth inclusions, which are possibly close-to-touching. Numerical experiments are presented to support the theoretical analysis and to illustrate the convergence of the proposed method for different shapes of inclusions in both two- and three-dimensional domains.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2023-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"109126945","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学术官方微信