A high-order finite-difference solver for direct numerical simulations of magnetohydrodynamic turbulence

IF 7.2 2区物理与天体物理 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computer Physics Communications Pub Date : 2024-10-15 DOI:10.1016/j.cpc.2024.109400

Jian Fang , Sylvain Laizet , Alex Skillen

{"title":"A high-order finite-difference solver for direct numerical simulations of magnetohydrodynamic turbulence","authors":"Jian Fang , Sylvain Laizet , Alex Skillen","doi":"10.1016/j.cpc.2024.109400","DOIUrl":null,"url":null,"abstract":"<div><div>This paper presents the development and validation of a Magnetohydrodynamics (MHD) module integrated into the Xcompact3d framework, an open-source high-order finite-difference suite of solvers designed to study turbulent flows on supercomputers. Leveraging the Fast Fourier Transform library already implemented in Xcompact3d, alongside sixth-order compact finite-difference schemes and a direct spectral Poisson solver, both the induction and potential-based MHD equations can be efficiently solved at scale on CPU-based supercomputers for fluids with strong and weak magnetic field, respectively. Validation of the MHD solver is conducted against established benchmarks, including Orszag-Tang vortex and MHD channel flows, demonstrating the module's capability to accurately capture complex MHD phenomena, providing a powerful tool for research in both engineering and astrophysics. The scalability of the Xcompact3d framework remains intact with the incorporation of the MHD module, ensuring efficient performance on modern high-performance clusters. This paper also presents new findings on the evolution of the Taylor-Green vortex under an external magnetic field for different flow regimes.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"307 ","pages":"Article 109400"},"PeriodicalIF":7.2000,"publicationDate":"2024-10-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computer Physics Communications","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0010465524003230","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}

引用次数: 0

Abstract

This paper presents the development and validation of a Magnetohydrodynamics (MHD) module integrated into the Xcompact3d framework, an open-source high-order finite-difference suite of solvers designed to study turbulent flows on supercomputers. Leveraging the Fast Fourier Transform library already implemented in Xcompact3d, alongside sixth-order compact finite-difference schemes and a direct spectral Poisson solver, both the induction and potential-based MHD equations can be efficiently solved at scale on CPU-based supercomputers for fluids with strong and weak magnetic field, respectively. Validation of the MHD solver is conducted against established benchmarks, including Orszag-Tang vortex and MHD channel flows, demonstrating the module's capability to accurately capture complex MHD phenomena, providing a powerful tool for research in both engineering and astrophysics. The scalability of the Xcompact3d framework remains intact with the incorporation of the MHD module, ensuring efficient performance on modern high-performance clusters. This paper also presents new findings on the evolution of the Taylor-Green vortex under an external magnetic field for different flow regimes.

查看原文本刊更多论文

用于磁流体动力湍流直接数值模拟的高阶有限差分求解器

本文介绍了集成到 Xcompact3d 框架中的磁流体动力学（MHD）模块的开发和验证情况，Xcompact3d 框架是一个开源的高阶有限差分求解器套件，旨在研究超级计算机上的湍流。利用 Xcompact3d 中已经实现的快速傅立叶变换库、六阶紧凑有限差分方案和直接谱泊松求解器，可以在基于 CPU 的超级计算机上高效求解基于感应和电势的 MHD 方程，分别用于具有强磁场和弱磁场的流体。MHD 求解器根据既定基准（包括 Orszag-Tang 涡流和 MHD 通道流）进行了验证，证明该模块能够准确捕捉复杂的 MHD 现象，为工程和天体物理学研究提供了强大的工具。加入 MHD 模块后，Xcompact3d 框架的可扩展性保持不变，确保了在现代高性能集群上的高效性能。本文还介绍了泰勒-格林涡旋在外部磁场作用下不同流动状态下的演变过程。

本文章由计算机程序翻译，如有差异，请以英文原文为准。

求助全文

约1分钟内获得全文求助全文

来源期刊

Computer Physics Communications 物理-计算机：跨学科应用

CiteScore

12.10

自引率

3.20%

发文量

287

审稿时长

5.3 months

期刊介绍： The focus of CPC is on contemporary computational methods and techniques and their implementation, the effectiveness of which will normally be evidenced by the author(s) within the context of a substantive problem in physics. Within this setting CPC publishes two types of paper. Computer Programs in Physics (CPiP) These papers describe significant computer programs to be archived in the CPC Program Library which is held in the Mendeley Data repository. The submitted software must be covered by an approved open source licence. Papers and associated computer programs that address a problem of contemporary interest in physics that cannot be solved by current software are particularly encouraged. Computational Physics Papers (CP) These are research papers in, but are not limited to, the following themes across computational physics and related disciplines. mathematical and numerical methods and algorithms; computational models including those associated with the design, control and analysis of experiments; and algebraic computation. Each will normally include software implementation and performance details. The software implementation should, ideally, be available via GitHub, Zenodo or an institutional repository.In addition, research papers on the impact of advanced computer architecture and special purpose computers on computing in the physical sciences and software topics related to, and of importance in, the physical sciences may be considered.