GREENY: A full-F 2D gyrofluid reconnection code

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

Computer Physics Communications Pub Date : 2025-08-07 DOI:10.1016/j.cpc.2025.109807

F.F. Locker , M. Held , T.M. Stocker-Waldhuber , A. Stürz , M. Rinner , A. Kendl

{"title":"GREENY: A full-F 2D gyrofluid reconnection code","authors":"F.F. Locker , M. Held , T.M. Stocker-Waldhuber , A. Stürz , M. Rinner , A. Kendl","doi":"10.1016/j.cpc.2025.109807","DOIUrl":null,"url":null,"abstract":"<div><div>We present a novel Full-<em>F</em> gyrofluid model and its implementation, the 2D gyrofluid magnetic reconnection code GREENY (Gyrofluid Reconnection with Extended Electromagnetic Nonlinearity). After a brief introduction to gyrofluids, magnetic reconnection, and the implemented models, we discuss the numerical framework and the algorithmic treatment of the quasi-neutrality condition and Amperè's law with special focus on arbitrary wavelength polarisation and induction. Next, we present solver tests, conservation laws, and the influence of artificial subgrid dissipation on Harris-sheet magnetic reconnection. Finally, we show different applications, initial conditions and present example simulations.</div></div><div><h3>Program summary</h3><div><em>Program Title:</em> GREENY</div><div><em>CPC Library link to program files:</em> <span><span>https://doi.org/10.17632/rc32ngbxft.1</span><svg><path></path></svg></span></div><div><em>Developer's repository link:</em> <span><span>https://git.uibk.ac.at/c7441315/greeny</span><svg><path></path></svg></span></div><div><em>Licensing provisions:</em> MIT</div><div><em>Programming language:</em> C/C++</div><div><em>Supplementary material:</em> Granalysis, <span>vortex_experiments.pdf</span>, <span>gyromod_derivation.pdf</span></div><div><em>Nature of problem:</em> Solves 2D isothermal electromagnetic gyrofluid equations of magnetic reconnection with self-consistent finite Larmor radius effects. The simulations can be done, using Full-<em>F</em>, <em>δ</em>F models with arbitrary wavelength polarisation or long-wavelength limit. To invert Ampère's law, multiple solvers, with and without Boussinesq-Oberbeck approximation, are available.</div><div><em>Solution method:</em> Finite difference solver for the dynamical gyrofluid density and momentum equations (Adams-Bashforth scheme, Arakawa scheme) with spectral and iterative solvers for evaluation of the gyrofluid polarisation equation, gyroaveraging operators and Ampère's law.</div><div><em>Additional comments including restrictions and unusual features:</em> Requires OpenMP, FFTW3 and NetCDF</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"316 ","pages":"Article 109807"},"PeriodicalIF":3.4000,"publicationDate":"2025-08-07","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/S0010465525003091","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

We present a novel Full-F gyrofluid model and its implementation, the 2D gyrofluid magnetic reconnection code GREENY (Gyrofluid Reconnection with Extended Electromagnetic Nonlinearity). After a brief introduction to gyrofluids, magnetic reconnection, and the implemented models, we discuss the numerical framework and the algorithmic treatment of the quasi-neutrality condition and Amperè's law with special focus on arbitrary wavelength polarisation and induction. Next, we present solver tests, conservation laws, and the influence of artificial subgrid dissipation on Harris-sheet magnetic reconnection. Finally, we show different applications, initial conditions and present example simulations.

Program summary

Program Title: GREENY

CPC Library link to program files: https://doi.org/10.17632/rc32ngbxft.1

Developer's repository link: https://git.uibk.ac.at/c7441315/greeny

Licensing provisions: MIT

Programming language: C/C++

Supplementary material: Granalysis, vortex_experiments.pdf, gyromod_derivation.pdf

Nature of problem: Solves 2D isothermal electromagnetic gyrofluid equations of magnetic reconnection with self-consistent finite Larmor radius effects. The simulations can be done, using Full-F, δF models with arbitrary wavelength polarisation or long-wavelength limit. To invert Ampère's law, multiple solvers, with and without Boussinesq-Oberbeck approximation, are available.

Solution method: Finite difference solver for the dynamical gyrofluid density and momentum equations (Adams-Bashforth scheme, Arakawa scheme) with spectral and iterative solvers for evaluation of the gyrofluid polarisation equation, gyroaveraging operators and Ampère's law.

Additional comments including restrictions and unusual features: Requires OpenMP, FFTW3 and NetCDF

查看原文本刊更多论文

绿色：全f二维回旋流体重连接代码

我们提出了一种新的全f陀螺仪流体模型及其实现，即二维陀螺仪流体磁重联代码GREENY（扩展电磁非线性的陀螺仪流体重联）。在简要介绍了回旋流体、磁重联和实现模型之后，我们讨论了数值框架和准中性条件的算法处理和Amperè定律，特别关注任意波长的极化和感应。接下来，我们给出了求解器测试、守恒定律以及人工亚网格耗散对哈里斯片磁重联的影响。最后给出了不同的应用、初始条件和实例仿真。程序摘要程序标题：GREENYCPC库链接到程序文件：https://doi.org/10.17632/rc32ngbxft.1Developer's存储库链接：https://git.uibk.ac.at/c7441315/greenyLicensing规定：mit编程语言：C/ c++补充材料：Granalysis， vortex_experiments.pdf, gyromod_derivation。pdf问题性质：求解具有自洽有限拉莫尔半径效应的二维等温电磁回旋流体重联方程。模拟可以使用Full-F， δF模型，具有任意波长偏振或长波长限制。为了反演安培特雷定律，有或没有Boussinesq-Oberbeck近似的多个求解器是可用的。求解方法：动态回旋流体密度和动量方程（Adams-Bashforth格式，Arakawa格式）的有限差分求解器，以及回旋流体极化方程、回旋平均算子和安培特雷定律的谱和迭代求解器。附加说明，包括限制和不寻常的功能：需要OpenMP， FFTW3和NetCDF

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

求助全文

约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.