当前位置:首页>python>AGATE 程序:一个用Python实现高性能磁流体动力学模拟的程序包

AGATE 程序:一个用Python实现高性能磁流体动力学模拟的程序包

  • 2026-08-22 05:58:13
AGATE 程序:一个用Python实现高性能磁流体动力学模拟的程序包


分享一个AGATE 程序:一个用Python实现高性能磁流体动力学模拟的程序包。

感谢论文的原作者!

================================

    🧠 为什么需要AGATE?

    磁流体动力学(MHD)模拟是研究空间物理、天体物理、聚变物理等领域的重要工具。传统HPC代码多采用C++或Fortran编写,如:

    • GAMERA、SWMF、OpenGGCM:Fortran

    • COCONUT、EUHFORIA:C++

    • Kokkos、Parthenon、FLASH、AMREX:C++

    这些代码性能卓越,但对普通科研人员而言,学习曲线陡峭、调试困难、扩展性差。

    Python近年来已成为最受欢迎的编程语言之一,但执行速度慢一直是其在HPC领域的短板。AGATE的出现,正是为了打破这一局面!

    🔬 支持的MHD方程

    1️⃣ 理想MHD & Hall MHD

    • 归一化:密度、磁场、长度、速度、压力、时间

    • 支持GLM散度清理

    • Hall速度 vH=−δiJ/ρ,δi 为归一化离子惯性长度

    • 波速:快磁声速 cfcf + 哨声波速 cwcw

    2️⃣ CGL(Chew-Goldberger-Low)各向异性MHD

    • 压力张量:

    • 守恒变量:能量密度 EE + CGL不变量 A

    • 压力各向异性限制:firehose、mirror、ion cyclotron不稳定性

    • 点隐式松弛:将各向异性拉回稳定区域

    图1. AGATE顶层类结构。

    图2. 对于CGL磁流体力学,AGATE成功再现了接触间断区域的物理量变化,以及慢激波处平行压强的选择性增强现象

    图3. 理想、霍尔及CGL磁流体力学下的奥尔扎格-唐涡旋演化。

    图4. 前三幅子图展示了CGL驻立磁声波测试中,三种各向异性比率下vₓ随位置x(横轴)和时间t(纵轴)变化的时空图。颜色表示vₓ的幅值(蓝色:负值,红色:正值)。

    ================================

    以上是我们分享的一些经验或者文章的搬运,或有不足,欢迎大家指出。若留言未回复,重要的消息可以留言再提醒一下。

    如有侵权,请联系我们立马删除!

    👇

    文章题目:

    High-Performance computational magnetohydrodynamics with Python

    https://doi.org/10.1016/j.cpc.2026.110077

    👇

    1. 【程序分享】TransOpt 程序:在恒定电子-声子耦合近似下求解半导体电输运性质的程序
    2. 【程序分享】VTAnDeM 程序:用于同时可视化材料的相稳定性、缺陷能量学和载流子浓度的 python 工具包
    3. 【程序分享】PyCDT 程序:用于对半导体和绝缘体中的点缺陷进行建模的 Python 工具包

    最新文章

    随机文章