BOUT++: A framework for parallel plasma fluid simulations
BOUT++: A framework for parallel plasma fluid simulations
复制标题
DOI：
10.1016/j.cpc.2009.03.008
复制
发表时间：
2009-09-01
影响因子：
6.3
通讯作者：
Wilson, H. R.
中科院分区：
文献类型：
--
作者：
Dudson, B. D.;Umansky, M. V.;Wilson, H. R.
提出了一种名为BOUT++的新型模块化代码，它在曲线坐标系中模拟三维流体方程。尽管该代码旨在模拟托卡马克X点几何中的边缘局域模（ELMs），但它能够在多种几何形状下模拟涉及任意数量标量场和矢量场的多种流体模型（磁化的和未磁化的）。时间演化是完全隐式的，并且实现了三阶加权本质非振荡（WENO）格式。给出了线性和非线性问题（奥萨格 - 唐涡旋）的基准测试，结果显示吻合良好。通过随问题规模和处理器数量的扩展来测试代码的性能，结果表明可高效扩展到数千个处理器。给出了使用简化理想磁流体动力学（MHD）对ELMs进行的线性初值模拟，并将结果与ELITE线性MHD特征值代码进行了比较。发现所得的模式结构和增长率非常吻合（BOUT++中的γ = 0.245ω(A)，ELITE中的γ = 0.239ω(A)，其中阿尔芬时间尺度1/ω(A)=R/V - A）。据我们所知，这是首次成功展示无耗散的ELMs初值模拟。（C）2009爱思唯尔有限公司保留所有权利。
A new modular code called BOUT++ is presented, which simulates 3D fluid equations in curvilinear coordinates. Although aimed at simulating Edge Localised Modes (ELMS) in tokamak x-point geometry, the code is able to simulate a wide range Of fluid models (magnetised and unmagnetised) involving an arbitrary number of scalar and vector fields, in a wide range of geometries. Time evolution is fully implicit, and 3rd-order WENO schemes are implemented. Benchmarks are presented for linear and nonlinear problems (the Orszag-Tang vortex) showing good agreement. Performance of the code is tested by scaling with problem size and processor number, showing efficient scaling to thousands of processors. Linear initial-value simulations of ELMS using reduced ideal MHD are presented, and the results compared to the ELITE linear MHD eigenvalue code. The resulting mode-structures and growth-rate are found to be in good agreement (gamma(BOUT++) = 0.245 omega(A), gamma(ELITE) = 0.239 omega(A), with Alfvenic timescale 1/omega(A) = R/V-A). To our knowledge, this is the first time dissipationless, initial-value Simulations of ELMS have been successfully demonstrated. (C) 2009 Elsevier B.V. All rights reserved.