Numerische Mathematik最新文献

筛选
英文 中文
Physics-preserving enriched Galerkin method for a fully-coupled thermo-poroelasticity model 针对全耦合热-弹塑性模型的物理保留富集伽勒金方法
IF 2.1 2区 数学
Numerische Mathematik Pub Date : 2024-05-03 DOI: 10.1007/s00211-024-01406-x
Son-Young Yi, Sanghyun Lee
{"title":"Physics-preserving enriched Galerkin method for a fully-coupled thermo-poroelasticity model","authors":"Son-Young Yi, Sanghyun Lee","doi":"10.1007/s00211-024-01406-x","DOIUrl":"https://doi.org/10.1007/s00211-024-01406-x","url":null,"abstract":"<p>This paper proposes a new numerical method for a fully-coupled, quasi-static thermo-poroelasticity model in a unified enriched Galerkin (EG) method framework. In our method, the mechanics sub-problem is solved using a locking-free EG method, and the flow and heat sub-problems are solved using a locally-conservative EG method. The proposed method offers mass and energy conservation properties with much lower costs than other methods with the same properties, including discontinuous Galerkin methods and mixed finite element methods. The well-posedness and optimal a priori error estimates are carefully derived. Several numerical tests confirm the theoretical optimal convergence rates and the mass and energy conservation properties of the new method.\u0000</p>","PeriodicalId":49733,"journal":{"name":"Numerische Mathematik","volume":"51 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-05-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140882596","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
High-order bounds-satisfying approximation of partial differential equations via finite element variational inequalities 通过有限元变分不等式实现偏微分方程的高阶边界满足逼近
IF 2.1 2区 数学
Numerische Mathematik Pub Date : 2024-04-30 DOI: 10.1007/s00211-024-01405-y
Robert C. Kirby, Daniel Shapero
{"title":"High-order bounds-satisfying approximation of partial differential equations via finite element variational inequalities","authors":"Robert C. Kirby, Daniel Shapero","doi":"10.1007/s00211-024-01405-y","DOIUrl":"https://doi.org/10.1007/s00211-024-01405-y","url":null,"abstract":"<p>Solutions to many important partial differential equations satisfy bounds constraints, but approximations computed by finite element or finite difference methods typically fail to respect the same conditions. Chang and Nakshatrala (Comput Methods Appl Mech Eng 320:287–334, 2017) enforce such bounds in finite element methods through the solution of variational inequalities rather than linear variational problems. Here, we provide a theoretical justification for this method, including higher-order discretizations. We prove an abstract best approximation result for the linear variational inequality and estimates showing that bounds-constrained polynomials provide comparable approximation power to standard spaces. For any unconstrained approximation to a function, there exists a constrained approximation which is comparable in the <span>(W^{1,p})</span> norm. In practice, one cannot efficiently represent and manipulate the entire family of bounds-constrained polynomials, but applying bounds constraints to the coefficients of a polynomial in the Bernstein basis guarantees those constraints on the polynomial. Although our theoretical results do not guaruntee high accuracy for this subset of bounds-constrained polynomials, numerical results indicate optimal orders of accuracy for smooth solutions and sharp resolution of features in convection–diffusion problems, all subject to bounds constraints.</p>","PeriodicalId":49733,"journal":{"name":"Numerische Mathematik","volume":"22 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-04-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140838974","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
Trefftz discontinuous Galerkin discretization for the Stokes problem 斯托克斯问题的 Trefftz 非连续伽勒金离散法
IF 2.1 2区 数学
Numerische Mathematik Pub Date : 2024-04-10 DOI: 10.1007/s00211-024-01404-z
Philip L. Lederer, Christoph Lehrenfeld, Paul Stocker
{"title":"Trefftz discontinuous Galerkin discretization for the Stokes problem","authors":"Philip L. Lederer, Christoph Lehrenfeld, Paul Stocker","doi":"10.1007/s00211-024-01404-z","DOIUrl":"https://doi.org/10.1007/s00211-024-01404-z","url":null,"abstract":"<p>We introduce a new discretization based on a polynomial Trefftz-DG method for solving the Stokes equations. Discrete solutions of this method fulfill the Stokes equations pointwise within each element and yield element-wise divergence-free solutions. Compared to standard DG methods, a strong reduction of the degrees of freedom is achieved, especially for higher polynomial degrees. In addition, in contrast to many other Trefftz-DG methods, our approach allows us to easily incorporate inhomogeneous right-hand sides (driving forces) by using the concept of the embedded Trefftz-DG method. On top of a detailed a priori error analysis, we further compare our approach to other (hybrid) discontinuous Galerkin Stokes discretizations and present numerical examples.\u0000</p>","PeriodicalId":49733,"journal":{"name":"Numerische Mathematik","volume":"6 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-04-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140593245","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 iterative method for the solution of Laplace-like equations in high and very high space dimensions 高维和超高维空间拉普拉斯方程的迭代求解方法
IF 2.1 2区 数学
Numerische Mathematik Pub Date : 2024-04-01 DOI: 10.1007/s00211-024-01401-2
{"title":"An iterative method for the solution of Laplace-like equations in high and very high space dimensions","authors":"","doi":"10.1007/s00211-024-01401-2","DOIUrl":"https://doi.org/10.1007/s00211-024-01401-2","url":null,"abstract":"<h3>Abstract</h3> <p>This paper deals with the equation <span> <span>(-varDelta u+mu u=f)</span> </span> on high-dimensional spaces <span> <span>({mathbb {R}}^m)</span> </span>, where the right-hand side <span> <span>(f(x)=F(Tx))</span> </span> is composed of a separable function <em>F</em> with an integrable Fourier transform on a space of a dimension <span> <span>(n&gt;m)</span> </span> and a linear mapping given by a matrix <em>T</em> of full rank and <span> <span>(mu ge 0)</span> </span> is a constant. For example, the right-hand side can explicitly depend on differences <span> <span>(x_i-x_j)</span> </span> of components of <em>x</em>. Following our publication (Yserentant in Numer Math 146:219–238, 2020), we show that the solution of this equation can be expanded into sums of functions of the same structure and develop in this framework an equally simple and fast iterative method for its computation. The method is based on the observation that in almost all cases and for large problem classes the expression <span> <span>(Vert T^tyVert ^2)</span> </span> deviates on the unit sphere <span> <span>(Vert yVert =1)</span> </span> the less from its mean value the higher the dimension <em>m</em> is, a concentration of measure effect. The higher the dimension <em>m</em>, the faster the iteration converges.</p>","PeriodicalId":49733,"journal":{"name":"Numerische Mathematik","volume":"5 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140593521","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
Locality of the windowed local density of states 加窗局部态密度的位置性
IF 2.1 2区 数学
Numerische Mathematik Pub Date : 2024-03-27 DOI: 10.1007/s00211-024-01400-3
Terry A. Loring, Jianfeng Lu, Alexander B. Watson
{"title":"Locality of the windowed local density of states","authors":"Terry A. Loring, Jianfeng Lu, Alexander B. Watson","doi":"10.1007/s00211-024-01400-3","DOIUrl":"https://doi.org/10.1007/s00211-024-01400-3","url":null,"abstract":"<p>We consider a generalization of local density of states which is “windowed” with respect to position and energy, called the windowed local density of states (wLDOS). This definition generalizes the usual LDOS in the sense that the usual LDOS is recovered in the limit where the position window captures individual sites and the energy window is a delta distribution. We prove that the wLDOS is local in the sense that it can be computed up to arbitrarily small error using spatial truncations of the system Hamiltonian. Using this result we prove that the wLDOS is well-defined and computable for infinite systems satisfying some natural assumptions. We finally present numerical computations of the wLDOS at the edge and in the bulk of a “Fibonacci SSH model”, a one-dimensional non-periodic model with topological edge states.\u0000</p>","PeriodicalId":49733,"journal":{"name":"Numerische Mathematik","volume":"60 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-03-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140313129","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 the number of terms in the COS method for European option pricing 关于欧式期权定价 COS 方法中的条款数
IF 2.1 2区 数学
Numerische Mathematik Pub Date : 2024-03-25 DOI: 10.1007/s00211-024-01402-1
Gero Junike
{"title":"On the number of terms in the COS method for European option pricing","authors":"Gero Junike","doi":"10.1007/s00211-024-01402-1","DOIUrl":"https://doi.org/10.1007/s00211-024-01402-1","url":null,"abstract":"<p>The Fourier-cosine expansion (COS) method is used to price European options numerically in a very efficient way. To apply the COS method, one has to specify two parameters: a truncation range for the density of the log-returns and a number of terms <i>N</i> to approximate the truncated density by a cosine series. How to choose the truncation range is already known. Here, we are able to find an explicit and useful bound for <i>N</i> as well for pricing and for the sensitivities, i.e., the Greeks Delta and Gamma, provided the density of the log-returns is smooth. We further show that the COS method has an exponential order of convergence when the density is smooth and decays exponentially. However, when the density is smooth and has heavy tails, as in the Finite Moment Log Stable model, the COS method does not have exponential order of convergence. Numerical experiments confirm the theoretical results.</p>","PeriodicalId":49733,"journal":{"name":"Numerische Mathematik","volume":"273 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-03-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140298701","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 structure-preserving parametric finite element method for geometric flows with anisotropic surface energy 具有各向异性表面能的几何流动的结构保持参数有限元方法
IF 2.1 2区 数学
Numerische Mathematik Pub Date : 2024-03-11 DOI: 10.1007/s00211-024-01398-8
Weizhu Bao, Yifei Li
{"title":"A structure-preserving parametric finite element method for geometric flows with anisotropic surface energy","authors":"Weizhu Bao, Yifei Li","doi":"10.1007/s00211-024-01398-8","DOIUrl":"https://doi.org/10.1007/s00211-024-01398-8","url":null,"abstract":"<p>We propose and analyze a structure-preserving parametric finite element method (SP-PFEM) for the evolution of a closed curve under different geometric flows with arbitrary anisotropic surface energy density <span>(gamma (varvec{n}))</span>, where <span>(varvec{n}in mathbb {S}^1)</span> represents the outward unit normal vector. We begin with the anisotropic surface diffusion which possesses two well-known geometric structures—area conservation and energy dissipation—during the evolution of the closed curve. By introducing a novel surface energy matrix <span>(varvec{G}_k(varvec{n}))</span> depending on <span>(gamma (varvec{n}))</span> and the Cahn-Hoffman <span>(varvec{xi })</span>-vector as well as a nonnegative stabilizing function <span>(k(varvec{n}))</span>, we obtain a new conservative geometric partial differential equation and its corresponding variational formulation for the anisotropic surface diffusion. Based on the new weak formulation, we propose a full discretization by adopting the parametric finite element method for spatial discretization and a semi-implicit temporal discretization with a proper and clever approximation for the outward normal vector. Under a mild and natural condition on <span>(gamma (varvec{n}))</span>, we can prove that the proposed full discretization is structure-preserving, i.e. it preserves the area conservation and energy dissipation at the discretized level, and thus it is unconditionally energy stable. The proposed SP-PFEM is then extended to simulate the evolution of a close curve under other anisotropic geometric flows including anisotropic curvature flow and area-conserved anisotropic curvature flow. Extensive numerical results are reported to demonstrate the efficiency and unconditional energy stability as well as good mesh quality (and thus no need to re-mesh during the evolution) of the proposed SP-PFEM for simulating anisotropic geometric flows.</p>","PeriodicalId":49733,"journal":{"name":"Numerische Mathematik","volume":"33 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-03-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140097371","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
Parabolic PDE-constrained optimal control under uncertainty with entropic risk measure using quasi-Monte Carlo integration 利用准蒙特卡洛积分实现具有熵风险度量的不确定性下抛物线 PDE 受限最优控制
IF 2.1 2区 数学
Numerische Mathematik Pub Date : 2024-03-11 DOI: 10.1007/s00211-024-01397-9
Philipp A. Guth, Vesa Kaarnioja, Frances Y. Kuo, Claudia Schillings, Ian H. Sloan
{"title":"Parabolic PDE-constrained optimal control under uncertainty with entropic risk measure using quasi-Monte Carlo integration","authors":"Philipp A. Guth, Vesa Kaarnioja, Frances Y. Kuo, Claudia Schillings, Ian H. Sloan","doi":"10.1007/s00211-024-01397-9","DOIUrl":"https://doi.org/10.1007/s00211-024-01397-9","url":null,"abstract":"<p>We study the application of a tailored quasi-Monte Carlo (QMC) method to a class of optimal control problems subject to parabolic partial differential equation (PDE) constraints under uncertainty: the state in our setting is the solution of a parabolic PDE with a random thermal diffusion coefficient, steered by a control function. To account for the presence of uncertainty in the optimal control problem, the objective function is composed with a risk measure. We focus on two risk measures, both involving high-dimensional integrals over the stochastic variables: the expected value and the (nonlinear) entropic risk measure. The high-dimensional integrals are computed numerically using specially designed QMC methods and, under moderate assumptions on the input random field, the error rate is shown to be essentially linear, independently of the stochastic dimension of the problem—and thereby superior to ordinary Monte Carlo methods. Numerical results demonstrate the effectiveness of our method.</p>","PeriodicalId":49733,"journal":{"name":"Numerische Mathematik","volume":"87 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-03-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140097473","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
Wave scattering problems in exterior domains with the method of fundamental solutions 用基本解法解决外部域中的波散射问题
IF 2.1 2区 数学
Numerische Mathematik Pub Date : 2024-02-28 DOI: 10.1007/s00211-024-01395-x
Carlos J. S. Alves, Pedro R. S. Antunes
{"title":"Wave scattering problems in exterior domains with the method of fundamental solutions","authors":"Carlos J. S. Alves, Pedro R. S. Antunes","doi":"10.1007/s00211-024-01395-x","DOIUrl":"https://doi.org/10.1007/s00211-024-01395-x","url":null,"abstract":"<p>The method of fundamental solutions has been mainly applied to wave scattering problems in bounded domains and to our knowledge there have not been works addressing density results for general shapes, or addressing the calculation of the complex resonance frequencies that occur in exterior problems. We prove density and convergence of the fundamental solutions approximation in the context of wave scattering problems, with and without a priori knowledge of the frequency, which is of particular importance to detect resonance frequencies for trapping domains. We also present several numerical results that illustrate the good performance of the method in the calculation of complex resonance frequencies for trapping and non trapping domains in 2D and 3D.</p>","PeriodicalId":49733,"journal":{"name":"Numerische Mathematik","volume":"2 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-02-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140008655","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 randomized operator splitting scheme inspired by stochastic optimization methods 受随机优化方法启发的随机算子分割方案
IF 2.1 2区 数学
Numerische Mathematik Pub Date : 2024-02-26 DOI: 10.1007/s00211-024-01396-w
Monika Eisenmann, Tony Stillfjord
{"title":"A randomized operator splitting scheme inspired by stochastic optimization methods","authors":"Monika Eisenmann, Tony Stillfjord","doi":"10.1007/s00211-024-01396-w","DOIUrl":"https://doi.org/10.1007/s00211-024-01396-w","url":null,"abstract":"<p>In this paper, we combine the operator splitting methodology for abstract evolution equations with that of stochastic methods for large-scale optimization problems. The combination results in a randomized splitting scheme, which in a given time step does not necessarily use all the parts of the split operator. This is in contrast to deterministic splitting schemes which always use every part at least once, and often several times. As a result, the computational cost can be significantly decreased in comparison to such methods. We rigorously define a randomized operator splitting scheme in an abstract setting and provide an error analysis where we prove that the temporal convergence order of the scheme is at least 1/2. We illustrate the theory by numerical experiments on both linear and quasilinear diffusion problems, using a randomized domain decomposition approach. We conclude that choosing the randomization in certain ways may improve the order to 1. This is as accurate as applying e.g. backward (implicit) Euler to the full problem, without splitting.</p>","PeriodicalId":49733,"journal":{"name":"Numerische Mathematik","volume":"30 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-02-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139978704","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学术官方微信