8 research outputs found
A GPU-enabled implicit Finite Volume solver for the ideal two-fluid plasma model on unstructured grids
This paper describes the main features of a pioneering unsteady solver for
simulating ideal two-fluid plasmas on unstructured grids, taking profit of
GPGPU (General-purpose computing on graphics processing units). The code, which
has been implemented within the open source COOLFluiD platform, is implicit,
second-order in time and space, relying upon a Finite Volume method for the
spatial discretization and a three-point backward Euler for the time
integration. In particular, the convective fluxes are computed by a multi-fluid
version of the AUSM+up scheme for the plasma equations, in combination with a
modified Rusanov scheme with tunable dissipation for the Maxwell equations.
Source terms are integrated with a one-point rule, using the cell-centered
value. Some critical aspects of the porting to GPU's are discussed, as well as
the performance of two open source linear system solvers (i.e. PETSc,
PARALUTION). The code design allows for computing both flux and source terms on
the GPU along with their Jacobian, giving a noticeable decrease in the
computational time in comparison with the original CPU-based solver. The code
has been tested in a wide range of mesh sizes and in three different systems,
each one with a different GPU. The increased performance (up to 14x) is
demonstrated in two representative 2D benchmarks: propagation of circularly
polarized waves and the more challenging Geospace Environmental Modeling (GEM)
magnetic reconnection challenge.Comment: 22 pages, 7 figure
衝突性および無衝突性の降着円盤における角運動量輸送メカニズムに関する研究
学位の種別: 課程博士審査委員会委員 : (主査)東京大学准教授 横山 央明, 東京大学教授 今村 剛, 東京大学教授 鈴木 建, 東京大学教授 藤本 正樹, 千葉大学教授 松元 亮治University of Tokyo(東京大学