BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Natsuki Hosono

Natsuki Hosono contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.EPastro-ph.IMastro-ph.GAastro-ph.COastro-ph.SRphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - Emerging
6works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

6 published item(s)

preprint2020arXiv

N-body Simulations of Ring Formation Process around the Dwarf Planet Haumea

Haumea is the only known trans-Neptunian object with a ring. The ring is in a position that causes a 3:1 spin-orbit resonance with the rotational period of Haumea, which has a triaxial shape. The non-axisymmetric gravitational field around Haumea is thought to affect the dynamics of the ring; however, the process of ring formation has not been elucidated. In this study, we analyze in some detail a potential ring formation scenario for Haumea. We first calculated the gravitational field around the triaxial ellipsoid and estimated the distance at which an object revolving around Haumea can exist stably using simulation that incorporated the time-varying gravitational field. The results of this simulation showed that the trajectory of the object became unstable just inside its current ring position. Next, we analytically derived the Roche radius for a rigid body revolving around a triaxial ellipsoid and showed that the Roche radius could be near the current ring position. Furthermore, as a parameter study, we performed N-body simulations using the coefficient of rubble pile restitution as a variable. Results demonstrated that, according to numerous parameters, the position of the Roche radius was near the current position of Haumea's ring. Based on these findings, we can assume that there is a high possibility that the ring formed in the region between the boundary of the unstable region of the orbit and the Roche radius. The scenario presented in this study could help explain the process by which Haumea's ring formed.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Can the Uranian Satellites Form from a Debris Disk Generated by a Giant Impact?

Hydrodynamic simulations of a giant impact to proto-Uranus indicated that such an impact could tilt its rotational axis and produce a circumplanetary debris disk beyond the corotation radius of Uranus. However, whether Uranian satellites can actually be formed from such a wide disk remains unclear. Herein, we modeled a wide debris disk of solids with several initial conditions inferred from the hydrodynamic simulations, and performed $N$-body simulations to investigate in-situ satellite formation from the debris disk. We also took account of orbital evolutions of satellites due to the planetary tides after the growth of satellites. We found that, in any case, the orbital distribution of the five major satellites could not be reproduced from the disk as long as the power index of its surface density is similar to that of the disk generated just after the giant impact: Satellites in the middle region obtained much larger masses than Ariel or Umbriel, while the outermost satellite did not grow to the mass of Oberon. Our results indicate that we should consider the thermal and viscous evolution of the evaporated disk after the giant impact to form the five major satellites through the in-situ formation scenario. On the other hand, the small inner satellites would be formed from the rings produced by the disrupted satellites which migrated from around the corotation radius of Uranus due to the planetary tides.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A comparison of SPH artificial viscosities and their impact on the Keplerian disk

Hydrodynamical simulations of rotating disk play important roles in the field of astrophysical and planetary science. Smoothed Particle Hydrodynamics (SPH) has been widely used for such simulations. It, however, has been known that with SPH, a cold and thin Kepler disk breaks up due to the unwanted angular momentum transfer. Two possible reasons have been suggested for this breaking up of the disk; the artificial viscosity (AV) and the numerical error in the evaluation of pressure gradient in SPH. Which one is dominant has been still unclear. In this paper, we investigate the reason for this rapid breaking up of the disk. We implemented most of popular formulations of AV and switches and measured the angular momentum transfer due to both AV and the error of SPH estimate of pressure gradient. We found that the angular momentum transfer due to AV at the inner edge triggers the breaking up of the disk. We also found that the classical von-Neumann-Richtmyer-Landshoff type AV with a high order estimate for $\nabla \cdot \vec{v}$ can maintain the disk for $\sim 100$ orbits even when used with the standard formulation of SPH.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Implementation and performance of FDPS: A Framework Developing Parallel Particle Simulation Codes

We present the basic idea, implementation, measured performance and performance model of FDPS (Framework for developing particle simulators). FDPS is an application-development framework which helps the researchers to develop particle-based simulation programs for large-scale distributed-memory parallel supercomputers. A particle-based simulation program for distributed-memory parallel computers needs to perform domain decomposition, redistribution of particles, and gathering of particle information for interaction calculation. Also, even if distributed-memory parallel computers are not used, in order to reduce the amount of computation, algorithms such as Barnes-Hut tree method should be used for long-range interactions. For short-range interactions, some methods to limit the calculation to neighbor particles are necessary. FDPS provides all of these necessary functions for efficient parallel execution of particle-based simulations as "templates", which are independent of the actual data structure of particles and the functional form of the interaction. By using FDPS, researchers can write their programs with the amount of work necessary to write a simple, sequential and unoptimized program of O(N^2) calculation cost, and yet the program, once compiled with FDPS, will run efficiently on large-scale parallel supercomputers. A simple gravitational N-body program can be written in around 120 lines. We report the actual performance of these programs and the performance model. The weak scaling performance is very good, and almost linear speedup was obtained for up to the full system of K computer. The minimum calculation time per timestep is in the range of 30 ms (N=10^7) to 300 ms (N=10^9). These are currently limited by the time for the calculation of the domain decomposition and communication necessary for the interaction calculation. We discuss how we can overcome these bottlenecks.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Giant Impact Simulations with Density Independent Smoothed Particle Hydrodynamics

At present, the giant impact (GI) is the most widely accepted model for the origin of the Moon. Most of the numerical simulations of GI have been carried out with the smoothed particle hydrodynamics (SPH) method. Recently, however, it has been pointed out that standard formulation of SPH (SSPH) has difficulties in the treatment of a contact discontinuity such as a core-mantle boundary and a free surface such as a planetary surface. This difficulty comes from the assumption of differentiability of density in SSPH. We have developed an alternative formulation of SPH, density independent SPH (DISPH), which is based on differentiability of pressure instead of density to solve the problem of a contact discontinuity. In this paper, we report the results of the GI simulations with DISPH and compare them with those obtained with SSPH. We found that the disk properties, such as mass and angular momentum produced by DISPH is different from that of SSPH. In general, the disks formed by DISPH are more compact: while formation of a smaller mass moon for low-oblique impacts is expected with DISPH, inhibition of ejection would promote formation of a larger mass moon for high-oblique impacts. Since only the improvement of core-mantle boundary significantly affects the properties of circumplanetary disks generated by GI and DISPH has not been significantly improved from SSPH for a free surface, we should be very careful when some conclusions are drawn from the numerical simulations for GI. And it is necessary to develop the numerical hydrodynamical scheme for GI that can properly treat the free surface as well as the contact discontinuity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Density Independent Smoothed Particle Hydrodynamics for Non-Ideal Equation of State

The smoothed particle hydrodynamics (SPH) method is a useful numerical tool for the study of a variety of astrophysical and planetlogical problems. However, it turned out that the standard SPH algorithm has problems in dealing with hydrodynamical instabilities. This problem is due to the assumption that the local density distribution is differentiable. In order to solve this problem, a new SPH formulation, which does not require the differentiability of the density, have been proposed. This new SPH method improved the treatment of hydrodynamical instabilities. This method, however, is applicable only to the equation of state (EOS) of the ideal gas. In this paper, we describe how to extend the new SPH method to non-ideal EOS. We present the results of various standard numerical tests for non-ideal EOS. Our new method works well for non-ideal EOS. We conclude that our new SPH can handle hydrodynamical instabilities for an arbitrary EOS and that it is an attractive alternative to the standard SPH.

0 citations0 reviewsNo code linkedOpen access