Bound-preserving and entropy stable enriched Galerkin methods for nonlinear hyperbolic equations

IF 3.8 2区 物理与天体物理 Q2 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS
Dmitri Kuzmin , Sanghyun Lee , Yi-Yung Yang
{"title":"Bound-preserving and entropy stable enriched Galerkin methods for nonlinear hyperbolic equations","authors":"Dmitri Kuzmin ,&nbsp;Sanghyun Lee ,&nbsp;Yi-Yung Yang","doi":"10.1016/j.jcp.2025.114323","DOIUrl":null,"url":null,"abstract":"<div><div>In this paper, we develop monolithic limiting techniques for enforcing nonlinear stability constraints in enriched Galerkin (EG) discretizations of nonlinear scalar hyperbolic equations. To achieve local mass conservation and gain control over the cell averages, the space of continuous (multi-)linear finite element approximations is enriched with piecewise-constant functions. The resulting spatial semi-discretization has the structure of a variational multiscale method. For linear advection equations, it is inherently stable but generally not bound preserving. To satisfy discrete maximum principles and ensure entropy stability in the nonlinear case, we use limiters adapted to the structure of our locally conservative EG method. The cell averages are constrained using a flux limiter, while the nodal values of the continuous component are constrained using a clip-and-scale limiting strategy for antidiffusive element contributions. The design and analysis of our new algorithms build on recent advances in the fields of convex limiting and algebraic entropy fixes for finite element methods. In addition to proving the claimed properties of the proposed approach, we conduct numerical studies for two-dimensional nonlinear hyperbolic problems. The numerical results demonstrate the ability of our limiters to prevent violations of the imposed constraints, while preserving the optimal order of accuracy in experiments with smooth solutions.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"541 ","pages":"Article 114323"},"PeriodicalIF":3.8000,"publicationDate":"2025-08-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computational Physics","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0021999125006059","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0

Abstract

In this paper, we develop monolithic limiting techniques for enforcing nonlinear stability constraints in enriched Galerkin (EG) discretizations of nonlinear scalar hyperbolic equations. To achieve local mass conservation and gain control over the cell averages, the space of continuous (multi-)linear finite element approximations is enriched with piecewise-constant functions. The resulting spatial semi-discretization has the structure of a variational multiscale method. For linear advection equations, it is inherently stable but generally not bound preserving. To satisfy discrete maximum principles and ensure entropy stability in the nonlinear case, we use limiters adapted to the structure of our locally conservative EG method. The cell averages are constrained using a flux limiter, while the nodal values of the continuous component are constrained using a clip-and-scale limiting strategy for antidiffusive element contributions. The design and analysis of our new algorithms build on recent advances in the fields of convex limiting and algebraic entropy fixes for finite element methods. In addition to proving the claimed properties of the proposed approach, we conduct numerical studies for two-dimensional nonlinear hyperbolic problems. The numerical results demonstrate the ability of our limiters to prevent violations of the imposed constraints, while preserving the optimal order of accuracy in experiments with smooth solutions.
非线性双曲型方程的保界和熵稳定富Galerkin方法
在本文中,我们开发了用于非线性标量双曲方程的富伽辽金(EG)离散化的非线性稳定性约束的单片限制技术。为了实现局部质量守恒并获得对单元平均值的控制,连续(多)线性有限元近似的空间被分段常数函数丰富。所得到的空间半离散化具有变分多尺度方法的结构。对于线性平流方程,它是固有稳定的,但通常不保界。为了满足非线性情况下的离散极大值原则和保证熵的稳定性,我们使用了适应局部保守EG方法结构的约束器。单元平均使用通量限制器进行约束,而连续分量的节点值使用抗扩散单元贡献的剪切和尺度限制策略进行约束。我们的新算法的设计和分析建立在凸极限和有限元素方法的代数熵固定领域的最新进展之上。除了证明所提出方法的性质外,我们还对二维非线性双曲型问题进行了数值研究。数值结果表明,我们的限制器能够防止违反所施加的约束,同时在光滑解的实验中保持最优的精度顺序。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
求助全文
约1分钟内获得全文 求助全文
来源期刊
Journal of Computational Physics
Journal of Computational Physics 物理-计算机:跨学科应用
CiteScore
7.60
自引率
14.60%
发文量
763
审稿时长
5.8 months
期刊介绍: Journal of Computational Physics thoroughly treats the computational aspects of physical problems, presenting techniques for the numerical solution of mathematical equations arising in all areas of physics. The journal seeks to emphasize methods that cross disciplinary boundaries. The Journal of Computational Physics also publishes short notes of 4 pages or less (including figures, tables, and references but excluding title pages). Letters to the Editor commenting on articles already published in this Journal will also be considered. Neither notes nor letters should have an abstract.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
确定
请完成安全验证×
copy
已复制链接
快去分享给好友吧!
我知道了
右上角分享
点击右上角分享
0
联系我们:info@booksci.cn Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。 Copyright © 2023 布克学术 All rights reserved.
京ICP备2023020795号-1
ghs 京公网安备 11010802042870号
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术官方微信