Structure preserving hybrid Finite Volume Finite Element method for compressible MHD

IF 3.8 2区 物理与天体物理 Q2 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Journal of Computational Physics Pub Date : 2024-12-20 DOI:10.1016/j.jcp.2024.113691
Francesco Fambri , Eric Sonnendrücker
{"title":"Structure preserving hybrid Finite Volume Finite Element method for compressible MHD","authors":"Francesco Fambri ,&nbsp;Eric Sonnendrücker","doi":"10.1016/j.jcp.2024.113691","DOIUrl":null,"url":null,"abstract":"<div><div>In this manuscript we present a novel and efficient numerical method for the compressible viscous and resistive MHD equations for all Mach number regimes. The time-integration strategy is a semi-implicit splitting, combined with a hybrid Finite Volume and Finite Element (FE) discretization in space. The nonlinear convection is solved by a robust explicit FV scheme, while the magneto-acoustic terms are treated implicitly in time. The resulting CFL stability condition depends only on the fluid velocity, and not on the Alfvénic and acoustic modes. The magneto-acoustic terms are discretized by compatible FE based on a continuous and a discrete de Rham complexes designed using Finite Element Exterior Calculus (FEEC). Thanks to the use of FEEC, energy stability, magnetic-helicity conservation and the divergence-free conditions can be preserved also at the discrete level. A very efficient splitting approach is used to separate the acoustic and the Alfvénic modes in such a fashion that the original symmetries of the PDE governing equations are preserved. In this way, the algorithm relies on the solution of linear, symmetric and positive-definite algebraic systems, that are very efficiently handled by the simple matrix-free conjugate-gradient method. The resulting algorithm showed to be robust and accurate in low and high Mach regimes even at large Courant numbers. Non-trivial tests are solved in one-, two- and three-space dimensions to confirm the robustness, accuracy, and the low-dissipative and conserving properties of the final algorithm. While the formulation of the method is very general, numerical results for a second-order accurate FV-FE scheme will be presented.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"523 ","pages":"Article 113691"},"PeriodicalIF":3.8000,"publicationDate":"2024-12-20","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/S0021999124009392","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 manuscript we present a novel and efficient numerical method for the compressible viscous and resistive MHD equations for all Mach number regimes. The time-integration strategy is a semi-implicit splitting, combined with a hybrid Finite Volume and Finite Element (FE) discretization in space. The nonlinear convection is solved by a robust explicit FV scheme, while the magneto-acoustic terms are treated implicitly in time. The resulting CFL stability condition depends only on the fluid velocity, and not on the Alfvénic and acoustic modes. The magneto-acoustic terms are discretized by compatible FE based on a continuous and a discrete de Rham complexes designed using Finite Element Exterior Calculus (FEEC). Thanks to the use of FEEC, energy stability, magnetic-helicity conservation and the divergence-free conditions can be preserved also at the discrete level. A very efficient splitting approach is used to separate the acoustic and the Alfvénic modes in such a fashion that the original symmetries of the PDE governing equations are preserved. In this way, the algorithm relies on the solution of linear, symmetric and positive-definite algebraic systems, that are very efficiently handled by the simple matrix-free conjugate-gradient method. The resulting algorithm showed to be robust and accurate in low and high Mach regimes even at large Courant numbers. Non-trivial tests are solved in one-, two- and three-space dimensions to confirm the robustness, accuracy, and the low-dissipative and conserving properties of the final algorithm. While the formulation of the method is very general, numerical results for a second-order accurate FV-FE scheme will be presented.
查看原文
分享 分享
微信好友 朋友圈 QQ好友 复制链接
本刊更多论文
求助全文
约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.
期刊最新文献
A two-stage two-derivative fourth order positivity-preserving discontinuous Galerkin method for hyperbolic conservation laws Unified gas-kinetic wave-particle method for multiscale flow simulation of partially ionized plasma Learning to solve PDEs with finite volume-informed neural networks in a data-free approach A discrete de Rham discretization of interface diffusion problems with application to the Leaky Dielectric Model An extension of the hybrid transport diffusion Monte Carlo in thermal radiative transfer: The difference formulation and the global optimal reference field
×
引用
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