Abstract
A new scheme that tightly couples kinetic turbulence codes across a spatial interface is introduced. This scheme evolves from considerations of competing strategies and down-selection. It is found that the use of a composite kinetic distribution function and fields with global boundary conditions as if the coupled code were one makes the coupling problem tractable. In contrast, coupling the two solutions from each code across the overlap region is found to be more difficult due to numerical dephasing of the turbulent solutions between two solvers. Another advantage of the new scheme is that the data movement can be limited to the 3D fluid quantities, instead of higher dimensional kinetic information, which is computationally more efficient for large scale simulations on leadership class computers.
Original language | English |
---|---|
Article number | 072308 |
Journal | Physics of Plasmas |
Volume | 25 |
Issue number | 7 |
DOIs | |
State | Published - Jul 1 2018 |
Funding
This manuscript is based upon work supported by the U.S. Department of Energy, Office of Science, Office of Fusion Energy Sciences, and has been authored by Princeton University under Contract No. DE-AC02-09CH11466 with the U.S. Department of Energy. This research used resources of the Oak Ridge Leadership Computing Facility at the Oak Ridge National Laboratory (ORNL) and resources of the National Energy Research Scientific Computing Center (NERSC), which are supported by the Office of Science of the U.S. Department of Energy under Contract Nos. DE-AC05-00OR22725 and DE-AC02-05CH11231, respectively.
Funders | Funder number |
---|---|
U.S. Department of Energy | |
Office of Science | |
Fusion Energy Sciences | |
Oak Ridge National Laboratory | |
Princeton University |