Numerische Mathematik最新文献_第2页

Fast solution of incompressible flow problems with two-level pressure approximation 用两级压力近似法快速解决不可压缩流动问题

IF 2.1 2区数学

Numerische Mathematik Pub Date : 2024-06-17 DOI: 10.1007/s00211-024-01420-z

Jennifer Pestana, David J. Silvester

引用次数: 0

A structure-preserving finite element method for the multi-phase Mullins–Sekerka problem with triple junctions 具有三重结点的多相穆林斯-塞克尔卡问题的结构保持有限元方法

IF 2.1 2区数学

Numerische Mathematik Pub Date : 2024-06-14 DOI: 10.1007/s00211-024-01414-x

Tokuhiro Eto, Harald Garcke, Robert Nürnberg

引用次数: 0

Moment-SoS methods for optimal transport problems 优化运输问题的矩-SoS 方法

IF 2.1 2区数学

Numerische Mathematik Pub Date : 2024-06-11 DOI: 10.1007/s00211-024-01422-x

Olga Mula, Anthony Nouy

{"title":"Moment-SoS methods for optimal transport problems","authors":"Olga Mula, Anthony Nouy","doi":"10.1007/s00211-024-01422-x","DOIUrl":"https://doi.org/10.1007/s00211-024-01422-x","url":null,"abstract":"Most common optimal transport (OT) solvers are currently based on an approximation of underlying measures by discrete measures. However, it is sometimes relevant to work only with moments of measures instead of the measure itself, and many common OT problems can be formulated as moment problems (the most relevant examples being (L^p)-Wasserstein distances, barycenters, and Gromov–Wasserstein discrepancies on Euclidean spaces). We leverage this fact to develop a generalized moment formulation that covers these classes of OT problems. The transport plan is represented through its moments on a given basis, and the marginal constraints are expressed in terms of moment constraints. A practical computation then consists in considering a truncation of the involved moment sequences up to a certain order, and using the polynomial sums-of-squares hierarchy for measures supported on semi-algebraic sets. We prove that the strategy converges to the solution of the OT problem as the order increases. We also show how to approximate linear quantities of interest, and how to estimate the support of the optimal transport map from the computed moments using Christoffel–Darboux kernels. Numerical experiments illustrate the good behavior of the approach.","PeriodicalId":49733,"journal":{"name":"Numerische Mathematik","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-06-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141503109","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 family of conforming finite element divdiv complexes on cuboid meshes 立方体网格上的保形有限元 divdiv 复合物系列

IF 2.1 2区数学

Numerische Mathematik Pub Date : 2024-06-06 DOI: 10.1007/s00211-024-01418-7

Jun Hu, Yizhou Liang, Rui Ma, Min Zhang

引用次数: 0

Validated integration of semilinear parabolic PDEs 半线性抛物线 PDE 的验证整合

IF 2.1 2区数学

Numerische Mathematik Pub Date : 2024-06-04 DOI: 10.1007/s00211-024-01415-w

Jan Bouwe van den Berg, Maxime Breden, Ray Sheombarsing

{"title":"Validated integration of semilinear parabolic PDEs","authors":"Jan Bouwe van den Berg, Maxime Breden, Ray Sheombarsing","doi":"10.1007/s00211-024-01415-w","DOIUrl":"https://doi.org/10.1007/s00211-024-01415-w","url":null,"abstract":"Integrating evolutionary partial differential equations (PDEs) is an essential ingredient for studying the dynamics of the solutions. Indeed, simulations are at the core of scientific computing, but their mathematical reliability is often difficult to quantify, especially when one is interested in the output of a given simulation, rather than in the asymptotic regime where the discretization parameter tends to zero. In this paper we present a computer-assisted proof methodology to perform rigorous time integration for scalar semilinear parabolic PDEs with periodic boundary conditions. We formulate an equivalent zero-finding problem based on a variation of constants formula in Fourier space. Using Chebyshev interpolation and domain decomposition, we then finish the proof with a Newton–Kantorovich type argument. The final output of this procedure is a proof of existence of an orbit, together with guaranteed error bounds between this orbit and a numerically computed approximation. We illustrate the versatility of the approach with results for the Fisher equation, the Swift–Hohenberg equation, the Ohta–Kawasaki equation and the Kuramoto–Sivashinsky equation. We expect that this rigorous integrator can form the basis for studying boundary value problems for connecting orbits in partial differential equations.","PeriodicalId":49733,"journal":{"name":"Numerische Mathematik","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-06-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141254147","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

Space-time CutFEM on overlapping meshes II: simple discontinuous mesh evolution 重叠网格上的时空 CutFEM II：简单不连续网格演化

IF 2.1 2区数学

Numerische Mathematik Pub Date : 2024-05-27 DOI: 10.1007/s00211-024-01413-y

Mats G. Larson, Carl Lundholm

{"title":"Space-time CutFEM on overlapping meshes II: simple discontinuous mesh evolution","authors":"Mats G. Larson, Carl Lundholm","doi":"10.1007/s00211-024-01413-y","DOIUrl":"https://doi.org/10.1007/s00211-024-01413-y","url":null,"abstract":"We present a cut finite element method for the heat equation on two overlapping meshes: a stationary background mesh and an overlapping mesh that evolves inside/“on top” of it. Here the overlapping mesh is prescribed by a simple discontinuous evolution, meaning that its location, size, and shape as functions of time are discontinuous and piecewise constant. For the discrete function space, we use continuous Galerkin in space and discontinuous Galerkin in time, with the addition of a discontinuity on the boundary between the two meshes. The finite element formulation is based on Nitsche’s method. The simple discontinuous mesh evolution results in a space-time discretization with a slabwise product structure between space and time which allows for existing analysis methodologies to be applied with only minor modifications. We follow the analysis methodology presented by Eriksson and Johnson (SIAM J Numer Anal 28(1):43–77, 1991; SIAM J Numer Anal 32(3):706–740, 1995). The greatest modification is the introduction of a Ritz-like “shift operator” that is used to obtain the discrete strong stability needed for the error analysis. The shift operator generalizes the original analysis to some methods for which the discrete subspace at one time does not lie in the space of the stiffness form at the subsequent time. The error analysis consists of an a priori error estimate that is of optimal order with respect to both time step and mesh size. We also present numerical results for a problem in one spatial dimension that verify the analytic error convergence orders.","PeriodicalId":49733,"journal":{"name":"Numerische Mathematik","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-05-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141167102","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

Space-time CutFEM on overlapping meshes I: simple continuous mesh motion 重叠网格上的时空 CutFEM I：简单连续网格运动

IF 2.1 2区数学

Numerische Mathematik Pub Date : 2024-05-25 DOI: 10.1007/s00211-024-01417-8

Mats G. Larson, Anders Logg, Carl Lundholm

{"title":"Space-time CutFEM on overlapping meshes I: simple continuous mesh motion","authors":"Mats G. Larson, Anders Logg, Carl Lundholm","doi":"10.1007/s00211-024-01417-8","DOIUrl":"https://doi.org/10.1007/s00211-024-01417-8","url":null,"abstract":"We present a cut finite element method for the heat equation on two overlapping meshes: a stationary background mesh and an overlapping mesh that moves around inside/“on top” of it. Here the overlapping mesh is prescribed by a simple continuous motion, meaning that its location as a function of time is continuous and piecewise linear. For the discrete function space, we use continuous Galerkin in space and discontinuous Galerkin in time, with the addition of a discontinuity on the boundary between the two meshes. The finite element formulation is based on Nitsche’s method and also includes an integral term over the space-time boundary between the two meshes that mimics the standard discontinuous Galerkin time-jump term. The simple continuous mesh motion results in a space-time discretization for which standard analysis methodologies either fail or are unsuitable. We therefore employ what seems to be a relatively uncommon energy analysis framework for finite element methods for parabolic problems that is general and robust enough to be applicable to the current setting. The energy analysis consists of a stability estimate that is slightly stronger than the standard basic one and an a priori error estimate that is of optimal order with respect to both time step and mesh size. We also present numerical results for a problem in one spatial dimension that verify the analytic error convergence orders.\u0000","PeriodicalId":49733,"journal":{"name":"Numerische Mathematik","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-05-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141148070","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

Drift approximation by the modified Boris algorithm of charged-particle dynamics in toroidal geometry 环状几何中带电粒子动力学的改进鲍里斯算法漂移近似

IF 2.1 2区数学

Numerische Mathematik Pub Date : 2024-05-24 DOI: 10.1007/s00211-024-01416-9

Yanyan Shi

引用次数: 0

Finite Elements with Switch Detection for direct optimal control of nonsmooth systems 有限元与开关检测，用于非平滑系统的直接优化控制

IF 2.1 2区数学

Numerische Mathematik Pub Date : 2024-05-24 DOI: 10.1007/s00211-024-01412-z

Armin Nurkanović, Mario Sperl, Sebastian Albrecht, Moritz Diehl

{"title":"Finite Elements with Switch Detection for direct optimal control of nonsmooth systems","authors":"Armin Nurkanović, Mario Sperl, Sebastian Albrecht, Moritz Diehl","doi":"10.1007/s00211-024-01412-z","DOIUrl":"https://doi.org/10.1007/s00211-024-01412-z","url":null,"abstract":"This paper introduces Finite Elements with Switch Detection (FESD), a numerical discretization method for nonsmooth differential equations. We consider the Filippov convexification of these systems and a transformation into dynamic complementarity systems introduced by Stewart (Numer Math 58(1):299–328, 1990). FESD is based on solving nonlinear complementarity problems and can automatically detect nonsmooth events in time. If standard time-stepping Runge–Kutta (RK) methods are naively applied to a nonsmooth ODE, the accuracy is at best of order one. In FESD, we let the integrator step size be a degree of freedom. Additional complementarity conditions, which we call cross complementarities, enable exact switch detection, hence FESD can recover the high order accuracy that the RK methods enjoy for smooth ODE. Additional conditions called step equilibration allow the step size to change only when switches occur and thus avoid spurious degrees of freedom. Convergence results for the FESD method are derived, local uniqueness of the solution and convergence of numerical sensitivities are proven. The efficacy of FESD is demonstrated in several simulation and optimal control examples. In an optimal control problem benchmark with FESD, we achieve up to five orders of magnitude more accurate solutions than a standard time-stepping approach for the same computational time.","PeriodicalId":49733,"journal":{"name":"Numerische Mathematik","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-05-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141148206","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

Super-localization of spatial network models 空间网络模型的超定位

IF 2.1 2区数学

Numerische Mathematik Pub Date : 2024-05-21 DOI: 10.1007/s00211-024-01410-1

Moritz Hauck, Axel Målqvist

引用次数: 0