Source author record

Daisuke Namekata

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

4works
4topics
4close collaborators

Published work

4 publications

preprint2019arXiv

Accelerated FDPS --- Algorithms to Use Accelerators with FDPS

In this paper, we describe the algorithms we implemented in FDPS to make efficient use of accelerator hardware such as GPGPUs. We have developed FDPS to make it possible for many researchers to develop their own high-performance parallel particle-based simulation programs without spending large amount of time for parallelization and performance tuning. The basic idea of FDPS is to provide a high-performance implementation of parallel algorithms for particle-based simulations in a "generic" form, so that researchers can define their own particle data structure and interparticle interaction functions and supply them to FDPS. FDPS compiled with user-supplied data type and interaction function provides all necessary functions for parallelization, and using those functions researchers can write their programs as though they are writing simple non-parallel program. It has been possible to use accelerators with FDPS, by writing the interaction function that uses the accelerator. However, the efficiency was limited by the latency and bandwidth of communication between the CPU and the accelerator and also by the mismatch between the available degree of parallelism of the interaction function and that of the hardware parallelism. We have modified the interface of user-provided interaction function so that accelerators are more efficiently used. We also implemented new techniques which reduce the amount of work on the side of CPU and amount of communication between CPU and accelerators. We have measured the performance of N-body simulations on a systems with NVIDIA Volta GPGPU using FDPS and the achieved performance is around 27 \% of the theoretical peak limit. We have constructed a detailed performance model, and found that the current implementation can achieve good performance on systems with much smaller memory and communication bandwidth.

preprint2016arXiv

Subparsec-scale dynamics of a dusty gas disk exposed to anisotropic AGN radiation with frequency-dependent radiative transfer

We explore the gas dynamics near the dust sublimation radius of active galactic nucleus (AGN). For the purpose, we perform axisymmetric radiation hydrodynamic simulations of a dusty gas disk of radius ≈1 pc\approx 1\,\mathrm{pc} around a supermassive black hole of mass 107 M⊙10^{7}\,\mathrm{M_{\odot}} taking into account (1) anisotropic radiation of accretion disk, (2) X-ray heating by corona, (3) radiative transfer of infrared (IR) photons reemitted by dust, (4) frequency dependency of direct and IR radiations, and (5) separate temperatures for gas and dust. As a result, we find that for Eddington ratio ≈0.77\approx 0.77, a nearly neutral, dense (≈106-⁡8  cm−3\approx 10^{6\operatorname{-}8}\;\mathrm{cm^{-3}}), geometrically-thin (h/r<0.06h/r<0.06) disk forms with a high velocity (≈200∼3000  km/s\approx 200 \sim 3000\;\mathrm{km/s}) dusty outflow launched from the disk surface. The disk temperature is determined by the balance between X-ray heating and various cooling, and the disk is almost supported by thermal pressure. Contrary to \citet{krolik07:_agn}, the radiation pressure by IR photons is not effective to thicken the disk, but rather compresses it. Thus, it seems difficult for a radiation-supported, geometrically-thick, obscuring torus to form near the dust sublimation radius as far as the Eddington ratio is high (∼1\sim 1). The mass outflow rate is 0.050.05-0.1  M⊙/yr0.1\;\mathrm{M_{\odot}}/\mathrm{yr} and the column density of the outflow is NH≲1021  cm−2N_{\mathrm{H}}\lesssim 10^{21}\;\mathrm{cm^{-2}}. To explain observed type-II AGN fraction, it is required that outflow gas is extended to larger radii (r≳10  pcr\gtrsim 10\;\mathrm{pc}) or that a denser dusty wind is launched from smaller radii (r∼104  Rgr\sim 10^{4}\;R_{g}).

preprint2014arXiv

On the evolution of gas clouds exposed to AGN radiation. I. Three-dimensional radiation hydrodynamic simulations

We perform three-dimensional radiation hydrodynamic simulations of uniform dusty gas clouds irradiated by an active galactic nucleus (AGN) to investigate the dependence of evolution of clouds on the ionization parameter U\mathcal{U} and the Str{ö}mgren number NS\mathcal{N}_{S}. We find that the evolution can be classified into two cases depending on U\mathcal{U}. In low U\mathcal{U} cases (U≈10−2\mathcal{U}\approx 10^{-2}), the evolution is mainly driven by photo-evaporation. A approximately spherically-symmetric evaporation flow with velocity of 100-⁡150  km  s−1100\operatorname{-}150\;\mathrm{km\;s^{-1}} is launched from the irradiated face. The cloud is compressed by a D-type shock with losing its mass due to photo-evaporation and is finally turned into a dense filament by t≲1.5tsct\lesssim 1.5t_{\mathrm{sc}}. In high U\mathcal{U} cases (U≈5×10−2\mathcal{U}\approx 5\times 10^{-2}), radiation pressure suppresses photo-evaporation from the central part of the irradiated face, reducing photo-evaporation rate. A evaporation flow from the outskirts of the irradiated face is turned into a high velocity (≲500  km  s−1\lesssim 500\;\mathrm{km\;s^{-1}}) gas wind because of radiation pressure on dust. The cloud is swept by a radiation pressure-driven shock and becomes a dense gas disk by t≈tsweept\approx t_{\mathrm{sweep}}. Star formation is expected in these dense regions for both cases of U\mathcal{U}. We discuss the influences of the AGN radiation on the clumpy torus. A simple estimate suggests that the clumps are destroyed in timescales shorter than their orbital periods. For the clumpy structure to be maintained over long period, the incident radiation field needs to be sufficiently weaken for most of the clumps, or, some mechanism that creates the clumps continuously is needed.

preprint2008arXiv

Mass Supply to Galactic Center due to Nested Bars in the Galaxy

We investigate rapid mass supply process by nested bars in the Galaxy by numerical simulation. We simulate gas flow in the whole galaxy disk with nested bars, which are the outer bar and the inner bar, especially with highly spatial resolution in the galactic central region. We assume two cases of inner bar size which are a smaller one and a larger one than the radius of the 200 pc gas ring which is corresponds to the Central Molecular Zone. From our numerical results, in the large size bar cases, the inner bars with large elongation induce sufficient mass inflow and destroy the 200 pc gas ring. On the other hand, in the small size bar cases, the inner bars with large elongation induce large mass inflow and do not destroy the 200 pc gas ring. This mass inflow is caused by straight shocks excited by the inner bar. In this case, nuclear gas disks of ~ 15 pc radius are formed. The nuclear gas disks are self-gravitationally unstable and we expect formation of compact star clusters under strong tidal force in the nuclear gas disks. We discuss evolution of the nuclear gas disk.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.