首页 > 最新文献

Journal of Numerical Mathematics最新文献

英文 中文
New non-augmented mixed finite element methods for the Navier–Stokes–Brinkman equations using Banach spaces 利用Banach空间求解Navier-Stokes-Brinkman方程的非增广混合有限元新方法
IF 3 2区 数学 Q1 Mathematics Pub Date : 2023-03-24 DOI: 10.1515/jnma-2022-0073
G. Gatica, Nicolás Núñez, R. Ruiz-Baier
Abstract In this paper we consider the Navier–Stokes–Brinkman equations, which constitute one of the most common nonlinear models utilized to simulate viscous fluids through porous media, and propose and analyze a Banach spaces-based approach yielding new mixed finite element methods for its numerical solution. In addition to the velocity and pressure, the strain rate tensor, the vorticity, and the stress tensor are introduced as auxiliary unknowns, and then the incompressibility condition is used to eliminate the pressure, which is computed afterwards by a postprocessing formula depending on the stress and the velocity. The resulting continuous formulation becomes a nonlinear perturbation of, in turn, a perturbed saddle point linear system, which is then rewritten as an equivalent fixed-point equation whose operator involved maps the velocity space into itself. The well-posedness of it is then analyzed by applying the classical Banach fixed point theorem, along with a smallness assumption on the data, the Babuška–Brezzi theory in Banach spaces, and a slight variant of a recently obtained solvability result for perturbed saddle point formulations in Banach spaces as well. The resulting Galerkin scheme is momentum-conservative. Its unique solvability is analyzed, under suitable hypotheses on the finite element subspaces, using a similar fixed-point strategy as in the continuous problem. A priori error estimates are rigorously derived, including also that for the pressure. We show that PEERS and AFW elements for the stress, the velocity and the rotation, together with piecewise polynomials of a proper degree for the strain rate tensor, yield stable discrete schemes. Then, the approximation properties of these subspaces and the Céa estimate imply the respective rates of convergence. Finally, we include two and three dimensional numerical experiments that serve to corroborate the theoretical findings, and these tests illustrate the performance of the proposed mixed finite element methods.
本文考虑了用于模拟多孔介质中粘性流体的最常用非线性模型之一的Navier-Stokes-Brinkman方程,提出并分析了一种基于Banach空间的方法,并对其数值解给出了新的混合有限元方法。除了速度和压力外,还引入应变率张量、涡量和应力张量作为辅助未知量,然后利用不可压缩条件消除压力,然后根据应力和速度的后处理公式计算压力。由此产生的连续公式变成一个非线性扰动,反过来,一个扰动鞍点线性系统,然后将其重写为一个等价的不动点方程,其算子涉及将速度空间映射到自身。然后通过应用经典的Banach不动点定理,以及对数据的小假设,Banach空间中的Babuška-Brezzi理论,以及最近获得的Banach空间中摄动鞍点公式的可解性结果的轻微变化来分析它的适定性。得到的伽辽金格式是动量保守的。在适当的有限元子空间假设下,采用与连续问题类似的不动点策略,分析了该问题的唯一可解性。严格推导了先验误差估计,包括对压力的估计。我们证明了应力、速度和旋转的PEERS和AFW单元,以及应变率张量的适当程度的分段多项式,可以产生稳定的离散格式。然后,这些子空间的近似性质和csama估计暗示了各自的收敛速度。最后,我们包括二维和三维数值实验,以证实理论结果,这些试验说明了所提出的混合有限元方法的性能。
{"title":"New non-augmented mixed finite element methods for the Navier–Stokes–Brinkman equations using Banach spaces","authors":"G. Gatica, Nicolás Núñez, R. Ruiz-Baier","doi":"10.1515/jnma-2022-0073","DOIUrl":"https://doi.org/10.1515/jnma-2022-0073","url":null,"abstract":"Abstract In this paper we consider the Navier–Stokes–Brinkman equations, which constitute one of the most common nonlinear models utilized to simulate viscous fluids through porous media, and propose and analyze a Banach spaces-based approach yielding new mixed finite element methods for its numerical solution. In addition to the velocity and pressure, the strain rate tensor, the vorticity, and the stress tensor are introduced as auxiliary unknowns, and then the incompressibility condition is used to eliminate the pressure, which is computed afterwards by a postprocessing formula depending on the stress and the velocity. The resulting continuous formulation becomes a nonlinear perturbation of, in turn, a perturbed saddle point linear system, which is then rewritten as an equivalent fixed-point equation whose operator involved maps the velocity space into itself. The well-posedness of it is then analyzed by applying the classical Banach fixed point theorem, along with a smallness assumption on the data, the Babuška–Brezzi theory in Banach spaces, and a slight variant of a recently obtained solvability result for perturbed saddle point formulations in Banach spaces as well. The resulting Galerkin scheme is momentum-conservative. Its unique solvability is analyzed, under suitable hypotheses on the finite element subspaces, using a similar fixed-point strategy as in the continuous problem. A priori error estimates are rigorously derived, including also that for the pressure. We show that PEERS and AFW elements for the stress, the velocity and the rotation, together with piecewise polynomials of a proper degree for the strain rate tensor, yield stable discrete schemes. Then, the approximation properties of these subspaces and the Céa estimate imply the respective rates of convergence. Finally, we include two and three dimensional numerical experiments that serve to corroborate the theoretical findings, and these tests illustrate the performance of the proposed mixed finite element methods.","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":3.0,"publicationDate":"2023-03-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"87723083","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}
引用次数: 4
Exploring numerical blow-up phenomena for the Keller–Segel–Navier–Stokes equations 探索Keller-Segel-Navier-Stokes方程的数值爆破现象
IF 3 2区 数学 Q1 Mathematics Pub Date : 2023-01-31 DOI: 10.48550/arXiv.2302.00139
J. Bonilla, J. V. Guti'errez-Santacreu
Abstract The Keller-Segel-Navier-Stokes system governs chemotaxis in liquid environments. This system is to be solved for the organism and chemoattractant densities and for the fluid velocity and pressure. It is known that if the total initial organism density mass is below 2π there exist globally defined generalised solutions, but what is less understood is whether there are blow-up solutions beyond such a threshold and its optimality. Motivated by this issue, a numerical blow-up scenario is investigated. Approximate solutions computed via a stabilised finite element method founded on a shock capturing technique are such that they satisfy a priori bounds as well as lower and L1(Ω) bounds for the organism and chemoattractant densities. In particular, these latter properties are essential in detecting numerical blow-up configurations, since the non-satisfaction of these two requirements might trigger numerical oscillations leading to non-realistic finite-time collapses into persistent Dirac-type measures. Our findings show that the existence threshold value 2π encountered for the organism density mass may not be optimal and hence it is conjectured that the critical threshold value 4π may be inherited from the fluid-free Keller-Segel equations. Additionally it is observed that the formation of singular points can be neglected if the fluid flow is intensified.
Keller-Segel-Navier-Stokes系统控制着液体环境中的趋化性。该系统需要求解生物体和化学引诱剂密度以及流体速度和压力。已知,如果总初始生物密度质量低于2π,则存在全局定义的广义解,但不太了解的是,是否存在超出此阈值的爆破解及其最优性。基于这一问题,本文研究了一个数值爆破场景。通过建立在冲击捕获技术基础上的稳定有限元方法计算的近似解是这样的,它们满足生物体和化学引诱剂密度的先验边界以及下限和L1(Ω)边界。特别是,后两种性质对于探测数值爆破构型至关重要,因为不满足这两种要求可能引发数值振荡,导致非现实的有限时间坍缩为持久的狄拉克型测量。我们的研究结果表明,生物密度质量遇到的存在阈值2π可能不是最优的,因此推测临界阈值4π可能继承自无流体的Keller-Segel方程。此外,还观察到,如果流体流动加剧,奇点的形成可以忽略不计。
{"title":"Exploring numerical blow-up phenomena for the Keller–Segel–Navier–Stokes equations","authors":"J. Bonilla, J. V. Guti'errez-Santacreu","doi":"10.48550/arXiv.2302.00139","DOIUrl":"https://doi.org/10.48550/arXiv.2302.00139","url":null,"abstract":"Abstract The Keller-Segel-Navier-Stokes system governs chemotaxis in liquid environments. This system is to be solved for the organism and chemoattractant densities and for the fluid velocity and pressure. It is known that if the total initial organism density mass is below 2π there exist globally defined generalised solutions, but what is less understood is whether there are blow-up solutions beyond such a threshold and its optimality. Motivated by this issue, a numerical blow-up scenario is investigated. Approximate solutions computed via a stabilised finite element method founded on a shock capturing technique are such that they satisfy a priori bounds as well as lower and L1(Ω) bounds for the organism and chemoattractant densities. In particular, these latter properties are essential in detecting numerical blow-up configurations, since the non-satisfaction of these two requirements might trigger numerical oscillations leading to non-realistic finite-time collapses into persistent Dirac-type measures. Our findings show that the existence threshold value 2π encountered for the organism density mass may not be optimal and hence it is conjectured that the critical threshold value 4π may be inherited from the fluid-free Keller-Segel equations. Additionally it is observed that the formation of singular points can be neglected if the fluid flow is intensified.","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":3.0,"publicationDate":"2023-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"75007947","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
Unifying a posteriori error analysis of five piecewise quadratic discretisations for the biharmonic equation 统一双调和方程五个分段二次离散的后验误差分析
IF 3 2区 数学 Q1 Mathematics Pub Date : 2023-01-31 DOI: 10.1515/jnma-2022-0092
C. Carstensen, Benedikt Gräßle, N. Nataraj
Abstract An abstract property (H) is the key to a complete a priori error analysis in the (discrete) energy norm for several nonstandard finite element methods in the recent work [Lowest-order equivalent nonstandard finite element methods for biharmonic plates, Carstensen and Nataraj, M2AN, 2022]. This paper investigates the impact of (H) to the a posteriori error analysis and establishes known and novel explicit residualbased a posteriori error estimates. The abstract framework applies to Morley, two versions of discontinuous Galerkin, C0 interior penalty, as well as weakly overpenalized symmetric interior penalty schemes for the biharmonic equation with a general source term in H−2(Ω).
摘要在最近的研究中[双谐板的最低阶等效非标准有限元方法,Carstensen和Nataraj, M2AN, 2022],一个抽象性质(H)是完成几种非标准有限元方法(离散)能量范数的先验误差分析的关键。本文研究了(H)对后验误差分析的影响,并建立了已知的和新的基于后验误差估计的显式残差。抽象框架适用于H−2中具有一般源项的双调和方程的Morley、两个版本的不连续Galerkin、C0内罚以及弱过罚对称内罚格式(Ω)。
{"title":"Unifying a posteriori error analysis of five piecewise quadratic discretisations for the biharmonic equation","authors":"C. Carstensen, Benedikt Gräßle, N. Nataraj","doi":"10.1515/jnma-2022-0092","DOIUrl":"https://doi.org/10.1515/jnma-2022-0092","url":null,"abstract":"Abstract An abstract property (H) is the key to a complete a priori error analysis in the (discrete) energy norm for several nonstandard finite element methods in the recent work [Lowest-order equivalent nonstandard finite element methods for biharmonic plates, Carstensen and Nataraj, M2AN, 2022]. This paper investigates the impact of (H) to the a posteriori error analysis and establishes known and novel explicit residualbased a posteriori error estimates. The abstract framework applies to Morley, two versions of discontinuous Galerkin, C0 interior penalty, as well as weakly overpenalized symmetric interior penalty schemes for the biharmonic equation with a general source term in H−2(Ω).","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":3.0,"publicationDate":"2023-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"83762618","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 POD-DEIM correction for Turing pattern approximation in reaction–diffusion PDE systems 反应扩散PDE系统图灵模式逼近的自适应POD-DEIM校正
2区 数学 Q1 Mathematics Pub Date : 2023-01-20 DOI: 10.1515/jnma-2022-0025
Alessandro Alla, Angela Monti, Ivonne Sgura
Abstract We investigate a suitable application of Model Order Reduction (MOR) techniques for the numerical approximation of Turing patterns, that are stationary solutions of reaction–diffusion PDE (RD-PDE) systems. We show that solutions of surrogate models built by classical Proper Orthogonal Decomposition (POD) exhibit an unstable error behaviour over the dimension of the reduced space. To overcome this drawback, first of all, we propose a POD-DEIM technique with a correction term that includes missing information in the reduced models. To improve the computational efficiency, we propose an adaptive version of this algorithm in time that accounts for the peculiar dynamics of the RD-PDE in presence of Turing instability. We show the effectiveness of the proposed methods in terms of accuracy and computational cost for a selection of RD systems, i.e., FitzHugh–Nagumo, Schnakenberg and the morphochemical DIB models, with increasing degree of nonlinearity and more structured patterns.
摘要研究了模型阶降简(MOR)技术在反应-扩散PDE (RD-PDE)系统稳态解图灵模式数值逼近中的合适应用。我们证明了由经典固有正交分解(POD)建立的代理模型的解在约简空间的维度上表现出不稳定的误差行为。为了克服这一缺点,首先,我们提出了一种POD-DEIM技术,该技术具有包含简化模型中缺失信息的校正项。为了提高计算效率,我们提出了该算法的自适应版本,该算法考虑了RD-PDE在存在图灵不稳定性时的特殊动力学。我们在准确性和计算成本方面展示了所提出方法的有效性,以选择RD系统,即FitzHugh-Nagumo, Schnakenberg和形态化学DIB模型,随着非线性程度的增加和更结构化的模式。
{"title":"Adaptive POD-DEIM correction for Turing pattern approximation in reaction–diffusion PDE systems","authors":"Alessandro Alla, Angela Monti, Ivonne Sgura","doi":"10.1515/jnma-2022-0025","DOIUrl":"https://doi.org/10.1515/jnma-2022-0025","url":null,"abstract":"Abstract We investigate a suitable application of Model Order Reduction (MOR) techniques for the numerical approximation of Turing patterns, that are stationary solutions of reaction–diffusion PDE (RD-PDE) systems. We show that solutions of surrogate models built by classical Proper Orthogonal Decomposition (POD) exhibit an unstable error behaviour over the dimension of the reduced space. To overcome this drawback, first of all, we propose a POD-DEIM technique with a correction term that includes missing information in the reduced models. To improve the computational efficiency, we propose an adaptive version of this algorithm in time that accounts for the peculiar dynamics of the RD-PDE in presence of Turing instability. We show the effectiveness of the proposed methods in terms of accuracy and computational cost for a selection of RD systems, i.e., FitzHugh–Nagumo, Schnakenberg and the morphochemical DIB models, with increasing degree of nonlinearity and more structured patterns.","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-01-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"134916455","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}
引用次数: 1
Transformed primal-dual methods for nonlinear saddle point systems 非线性鞍点系统的变换原对偶方法
2区 数学 Q1 Mathematics Pub Date : 2023-01-15 DOI: 10.1515/jnma-2022-0056
Long Chen, Jingrong Wei
Abstract A transformed primal-dual (TPD) flow is developed for a class of nonlinear smooth saddle point system. The flow for the dual variable contains a Schur complement which is strongly convex. Exponential stability of the saddle point is obtained by showing the strong Lyapunov property. Several TPD iterations are derived by implicit Euler, explicit Euler, implicit-explicit and Gauss-Seidel methods with accelerated overrelaxation of the TPD flow. Generalized to the symmetric TPD iterations, linear convergence rate is preserved for convex-concave saddle point systems under assumptions that the regularized functions are strongly convex. The effectiveness of augmented Lagrangian methods can be explained as a regularization of the non-strongly convexity and a preconditioning for the Schur complement. The algorithm and convergence analysis depends crucially on appropriate inner products of the spaces for the primal variable and dual variable. A clear convergence analysis with nonlinear inexact inner solvers is also developed.
摘要针对一类非线性光滑鞍点系统,建立了一种变换的原对偶流。对偶变量流包含一个强凸的Schur补。通过证明强李雅普诺夫性质,得到了鞍点的指数稳定性。采用隐式欧拉法、显式欧拉法、隐式-显式法和高斯-赛德尔法推导了TPD流的多次迭代,并加速了TPD流的超松弛。推广到对称TPD迭代,在正则函数为强凸的假设下,凸凹鞍点系统保持线性收敛速率。增广拉格朗日方法的有效性可以解释为非强凸性的正则化和舒尔补的先决条件。该算法和收敛性分析关键取决于原变量和对偶变量空间的适当内积。并给出了非线性非精确内解的清晰收敛分析。
{"title":"Transformed primal-dual methods for nonlinear saddle point systems","authors":"Long Chen, Jingrong Wei","doi":"10.1515/jnma-2022-0056","DOIUrl":"https://doi.org/10.1515/jnma-2022-0056","url":null,"abstract":"Abstract A transformed primal-dual (TPD) flow is developed for a class of nonlinear smooth saddle point system. The flow for the dual variable contains a Schur complement which is strongly convex. Exponential stability of the saddle point is obtained by showing the strong Lyapunov property. Several TPD iterations are derived by implicit Euler, explicit Euler, implicit-explicit and Gauss-Seidel methods with accelerated overrelaxation of the TPD flow. Generalized to the symmetric TPD iterations, linear convergence rate is preserved for convex-concave saddle point systems under assumptions that the regularized functions are strongly convex. The effectiveness of augmented Lagrangian methods can be explained as a regularization of the non-strongly convexity and a preconditioning for the Schur complement. The algorithm and convergence analysis depends crucially on appropriate inner products of the spaces for the primal variable and dual variable. A clear convergence analysis with nonlinear inexact inner solvers is also developed.","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":0.0,"publicationDate":"2023-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135743985","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}
引用次数: 1
Diagonally implicit Runge-Kutta schemes: Discrete energy-balance laws and compactness properties 对角隐式龙格-库塔格式:离散能量平衡定律和紧性
IF 3 2区 数学 Q1 Mathematics Pub Date : 2022-12-24 DOI: 10.1515/jnma-2022-0069
Abner J. Salgado, Ignacio Tomas
We study diagonally implicit Runge-Kutta (DIRK) schemes when applied to abstract evolution problems that fit into the Gelfand-triple framework. We introduce novel stability notions that are well-suited to this setting and provide simple, necessary and sufficient, conditions to verify that a DIRK scheme is stable in our sense and in Bochner-type norms. We use several popular DIRK schemes in order to illustrate cases that satisfy the required structural stability properties and cases that do not. In addition, under some mild structural conditions on the problem we can guarantee compactness of families of discrete solutions with respect to time discretization.
我们研究了对角隐式龙格-库塔(DIRK)格式在适合Gelfand-triple框架的抽象进化问题中的应用。我们引入了新的稳定性概念,非常适合于这种情况,并提供了简单的、必要的和充分的条件来验证DIRK方案在我们的意义上和在bochner型范数下是稳定的。我们使用几种流行的DIRK方案来说明满足结构稳定性要求的情况和不满足结构稳定性要求的情况。此外,在一些温和的结构条件下,我们可以保证离散解族相对于时间离散的紧性。
{"title":"Diagonally implicit Runge-Kutta schemes: Discrete energy-balance laws and compactness properties","authors":"Abner J. Salgado, Ignacio Tomas","doi":"10.1515/jnma-2022-0069","DOIUrl":"https://doi.org/10.1515/jnma-2022-0069","url":null,"abstract":"We study diagonally implicit Runge-Kutta (DIRK) schemes when applied to abstract evolution problems that fit into the Gelfand-triple framework. We introduce novel stability notions that are well-suited to this setting and provide simple, necessary and sufficient, conditions to verify that a DIRK scheme is stable in our sense and in Bochner-type norms. We use several popular DIRK schemes in order to illustrate cases that satisfy the required structural stability properties and cases that do not. In addition, under some mild structural conditions on the problem we can guarantee compactness of families of discrete solutions with respect to time discretization.","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":3.0,"publicationDate":"2022-12-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138530262","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
Error analysis for a Crouzeix–Raviart approximation of the p-Dirichlet problem p-Dirichlet问题的Crouzeix-Raviart近似的误差分析
IF 3 2区 数学 Q1 Mathematics Pub Date : 2022-10-21 DOI: 10.48550/arXiv.2210.12116
A. Kaltenbach
Abstract In the present paper, we examine a Crouzeix–Raviart approximation for non-linear partial differential equations having a (p, δ)-structure for some p ∈ (1, ∞) and δ⩾0. We establish a priori error estimates, which are optimal for all p ∈ (1, ∞) and δ⩾0, medius error estimates, i.e., best-approximation results, and a primal-dual a posteriori error estimate, which is both reliable and efficient. The theoretical findings are supported by numerical experiments.
在本文中,我们研究了对于某些p∈(1,∞)和δ大于或等于0具有(p, δ)结构的非线性偏微分方程的Crouzeix-Raviart近似。我们建立了先验误差估计,这对于所有p∈(1,∞)和δ小于或等于0是最优的,中等误差估计,即最佳近似结果,以及原始-对偶后验误差估计,这既可靠又有效。理论结果得到数值实验的支持。
{"title":"Error analysis for a Crouzeix–Raviart approximation of the p-Dirichlet problem","authors":"A. Kaltenbach","doi":"10.48550/arXiv.2210.12116","DOIUrl":"https://doi.org/10.48550/arXiv.2210.12116","url":null,"abstract":"Abstract In the present paper, we examine a Crouzeix–Raviart approximation for non-linear partial differential equations having a (p, δ)-structure for some p ∈ (1, ∞) and δ⩾0. We establish a priori error estimates, which are optimal for all p ∈ (1, ∞) and δ⩾0, medius error estimates, i.e., best-approximation results, and a primal-dual a posteriori error estimate, which is both reliable and efficient. The theoretical findings are supported by numerical experiments.","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":3.0,"publicationDate":"2022-10-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"85113245","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}
引用次数: 4
A subspace of linear nonconforming finite element for nearly incompressible elasticity and Stokes flow 近似不可压缩弹性和Stokes流的线性非协调有限元子空间
IF 3 2区 数学 Q1 Mathematics Pub Date : 2022-09-14 DOI: 10.1515/jnma-2022-0010
Shangyou Zhang
Abstract The linear nonconforming finite element, combined with constant finite element for pressure, is stable for the Stokes problem. But it does not satisfy the discrete Korn inequality. The linear conforming finite element satisfies the discrete Korn inequality, but is not stable for the Stokes problem and fails for the nearly incompressible elasticity problems. We enrich the linear conforming finite element by some nonconforming P1 bubbles, i.e., select a subspace of the linear nonconforming finite element space, so that the resulting linear nonconforming element is both stable and conforming enough to satisfy the Korn inequality, on HTC-type triangular and tetrahedral grids. Numerical tests in 2D and 3D are presented, confirming the analysis.
对于Stokes问题,线性非协调有限元与压力的恒定有限元相结合是稳定的。但它不满足离散Korn不等式。线性拟合有限元满足离散Korn不等式,但对于Stokes问题不稳定,对于近不可压缩弹性问题不稳定。我们通过一些不协调P1泡来丰富线性不协调有限元,即在线性不协调有限元空间中选择一个子空间,使得到的线性不协调单元在htc型三角形和四面体网格上既稳定又协调,足以满足Korn不等式。给出了二维和三维数值试验,验证了分析结果。
{"title":"A subspace of linear nonconforming finite element for nearly incompressible elasticity and Stokes flow","authors":"Shangyou Zhang","doi":"10.1515/jnma-2022-0010","DOIUrl":"https://doi.org/10.1515/jnma-2022-0010","url":null,"abstract":"Abstract The linear nonconforming finite element, combined with constant finite element for pressure, is stable for the Stokes problem. But it does not satisfy the discrete Korn inequality. The linear conforming finite element satisfies the discrete Korn inequality, but is not stable for the Stokes problem and fails for the nearly incompressible elasticity problems. We enrich the linear conforming finite element by some nonconforming P1 bubbles, i.e., select a subspace of the linear nonconforming finite element space, so that the resulting linear nonconforming element is both stable and conforming enough to satisfy the Korn inequality, on HTC-type triangular and tetrahedral grids. Numerical tests in 2D and 3D are presented, confirming the analysis.","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":3.0,"publicationDate":"2022-09-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"75820637","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
Transformed primal-dual methods for nonlinear saddle point systems 非线性鞍点系统的变换原对偶方法
IF 3 2区 数学 Q1 Mathematics Pub Date : 2022-08-04 DOI: 10.48550/arXiv.2208.02444
Long Chen, Jingrong Wei
Abstract A transformed primal–dual (TPD) flow is developed for a class of nonlinear smooth saddle point systemThe flow for the dual variable contains a Schur complement which is strongly convex. Exponential stability of the saddle point is obtained by showing the strong Lyapunov property. Several TPD iterations are derived by implicit Euler, explicit Euler, implicit–explicit, and Gauss–Seidel methods with accelerated overrelaxation of the TPD flow. Generalized to the symmetric TPD iterations, linear convergence rate is preserved for convex–concave saddle point systems under assumptions that the regularized functions are strongly convex. The effectiveness of augmented Lagrangian methods can be explained as a regularization of the non-strongly convexity and a preconditioning for the Schur complement. The algorithm and convergence analysis depends crucially on appropriate inner products of the spaces for the primal variable and dual variable. A clear convergence analysis with nonlinear inexact inner solvers is also developed.
摘要针对一类非线性光滑鞍点系统,建立了一种变换的原对偶流。对偶变量流包含一个强凸的Schur补。通过证明强李雅普诺夫性质,得到了鞍点的指数稳定性。采用隐式欧拉法、显式欧拉法、隐式-显式法和高斯-赛德尔法推导了TPD流的多次迭代,并加速了TPD流的超松弛。推广到对称TPD迭代,在正则函数为强凸的假设下,凸凹鞍点系统保持线性收敛速率。增广拉格朗日方法的有效性可以解释为非强凸性的正则化和舒尔补的先决条件。该算法和收敛性分析关键取决于原变量和对偶变量空间的适当内积。并给出了非线性非精确内解的清晰收敛分析。
{"title":"Transformed primal-dual methods for nonlinear saddle point systems","authors":"Long Chen, Jingrong Wei","doi":"10.48550/arXiv.2208.02444","DOIUrl":"https://doi.org/10.48550/arXiv.2208.02444","url":null,"abstract":"Abstract A transformed primal–dual (TPD) flow is developed for a class of nonlinear smooth saddle point systemThe flow for the dual variable contains a Schur complement which is strongly convex. Exponential stability of the saddle point is obtained by showing the strong Lyapunov property. Several TPD iterations are derived by implicit Euler, explicit Euler, implicit–explicit, and Gauss–Seidel methods with accelerated overrelaxation of the TPD flow. Generalized to the symmetric TPD iterations, linear convergence rate is preserved for convex–concave saddle point systems under assumptions that the regularized functions are strongly convex. The effectiveness of augmented Lagrangian methods can be explained as a regularization of the non-strongly convexity and a preconditioning for the Schur complement. The algorithm and convergence analysis depends crucially on appropriate inner products of the spaces for the primal variable and dual variable. A clear convergence analysis with nonlinear inexact inner solvers is also developed.","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":3.0,"publicationDate":"2022-08-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"74265900","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}
引用次数: 4
The deal.II library, Version 9.4 这笔交易。II库,版本9.4
IF 3 2区 数学 Q1 Mathematics Pub Date : 2022-07-17 DOI: 10.1515/jnma-2022-0054
D. Arndt, W. Bangerth, Marco Feder, M. Fehling, Rene Gassmöller, T. Heister, L. Heltai, M. Kronbichler, Matthias Maier, Peter Munch, Jean-Paul Pelteret, S. Sticko, Bruno Turcksin, David R. Wells
Abstract This paper provides an overview of the new features of the finite element library deal.II, version 9.4.
摘要本文概述了有限元库协议的新特点。II,版本9.4。
{"title":"The deal.II library, Version 9.4","authors":"D. Arndt, W. Bangerth, Marco Feder, M. Fehling, Rene Gassmöller, T. Heister, L. Heltai, M. Kronbichler, Matthias Maier, Peter Munch, Jean-Paul Pelteret, S. Sticko, Bruno Turcksin, David R. Wells","doi":"10.1515/jnma-2022-0054","DOIUrl":"https://doi.org/10.1515/jnma-2022-0054","url":null,"abstract":"Abstract This paper provides an overview of the new features of the finite element library deal.II, version 9.4.","PeriodicalId":50109,"journal":{"name":"Journal of Numerical Mathematics","volume":null,"pages":null},"PeriodicalIF":3.0,"publicationDate":"2022-07-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"88943081","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}
引用次数: 113
期刊
Journal of Numerical Mathematics
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1