仿星器MHD模拟的非轴对称柔性网格有限差分格式

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

Computer Physics Communications Pub Date : 2025-07-18 DOI:10.1016/j.cpc.2025.109776

J. Wang , Y. Zhou , J.M. Duan , Z.W. Ma , W. Zhang

{"title":"仿星器MHD模拟的非轴对称柔性网格有限差分格式","authors":"J. Wang , Y. Zhou , J.M. Duan , Z.W. Ma , W. Zhang","doi":"10.1016/j.cpc.2025.109776","DOIUrl":null,"url":null,"abstract":"<div><div>The adaptive moving mesh CLT code is extended to be applicable for the stellarator magneto-hydrodynamic (MHD) simulations. Compared with the tokamak version, the mesh can not only be non-uniform, but can also be non-axisymmetric and strongly shaped with a concave boundary. The extra toroidal transformation from the physical domain to the computational domain has to be taken into consideration. In the computational domain, the fourth-order finite difference scheme is constructed on a Cartesian computational mesh. To verify the code, we simulate the internal kink mode in a W7-X configuration, and benchmark the linear and non-linear results with the M3D-C<sup>1</sup> code. The method can also be used for different 3D equilibria, and a calculation of the resistive ballooning mode in a NCSX equilibrium is given.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"316 ","pages":"Article 109776"},"PeriodicalIF":3.4000,"publicationDate":"2025-07-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"CLT-fx: Non-axisymmetric flexible mesh finite difference scheme for stellarator MHD simulations\",\"authors\":\"J. Wang , Y. Zhou , J.M. Duan , Z.W. Ma , W. Zhang\",\"doi\":\"10.1016/j.cpc.2025.109776\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<div><div>The adaptive moving mesh CLT code is extended to be applicable for the stellarator magneto-hydrodynamic (MHD) simulations. Compared with the tokamak version, the mesh can not only be non-uniform, but can also be non-axisymmetric and strongly shaped with a concave boundary. The extra toroidal transformation from the physical domain to the computational domain has to be taken into consideration. In the computational domain, the fourth-order finite difference scheme is constructed on a Cartesian computational mesh. To verify the code, we simulate the internal kink mode in a W7-X configuration, and benchmark the linear and non-linear results with the M3D-C<sup>1</sup> code. The method can also be used for different 3D equilibria, and a calculation of the resistive ballooning mode in a NCSX equilibrium is given.</div></div>\",\"PeriodicalId\":285,\"journal\":{\"name\":\"Computer Physics Communications\",\"volume\":\"316 \",\"pages\":\"Article 109776\"},\"PeriodicalIF\":3.4000,\"publicationDate\":\"2025-07-18\",\"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/S0010465525002784\",\"RegionNum\":2,\"RegionCategory\":\"物理与天体物理\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q1\",\"JCRName\":\"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computer Physics Communications","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0010465525002784","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}

引用次数: 0

摘要

将自适应运动网格CLT代码扩展到仿星器磁流体动力学（MHD）仿真中。与托卡马克版本相比，网格不仅可以是不均匀的，而且可以是非轴对称的，并且具有凹边界的强烈形状。必须考虑从物理域到计算域的额外环面变换。在计算域中，在笛卡尔计算网格上构造了四阶有限差分格式。为了验证代码，我们在W7-X配置中模拟了内部扭结模式，并使用M3D-C1代码对线性和非线性结果进行了基准测试。该方法也可用于不同的三维平衡，并给出了NCSX平衡中阻力气胀模式的计算。

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

查看原文本刊更多论文

CLT-fx: Non-axisymmetric flexible mesh finite difference scheme for stellarator MHD simulations

The adaptive moving mesh CLT code is extended to be applicable for the stellarator magneto-hydrodynamic (MHD) simulations. Compared with the tokamak version, the mesh can not only be non-uniform, but can also be non-axisymmetric and strongly shaped with a concave boundary. The extra toroidal transformation from the physical domain to the computational domain has to be taken into consideration. In the computational domain, the fourth-order finite difference scheme is constructed on a Cartesian computational mesh. To verify the code, we simulate the internal kink mode in a W7-X configuration, and benchmark the linear and non-linear results with the M3D-C¹ code. The method can also be used for different 3D equilibria, and a calculation of the resistive ballooning mode in a NCSX equilibrium is given.

求助全文

通过发布文献求助，成功后即可免费获取论文全文。去求助

来源期刊

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.