-
Leptonic and semi-leptonic neutrino interactions with muons in the proto-neutron star cooling
Authors:
Ken'ichi Sugiura,
Shun Furusawa,
Kohsuke Sumiyoshi,
Shoichi Yamada
Abstract:
It is known that muons are scarce just after the birth of a proto-neutron star via a supernova explosion but get more abundant as the proto-neutron star cools via neutrino emissions on the Kelvin-Helmholtz timescale.
In this paper we evaluate all the relevant rates of the neutrino interactions with muons at different times in the proto-neutron star cooling.
We are particularly interested in th…
▽ More
It is known that muons are scarce just after the birth of a proto-neutron star via a supernova explosion but get more abundant as the proto-neutron star cools via neutrino emissions on the Kelvin-Helmholtz timescale.
In this paper we evaluate all the relevant rates of the neutrino interactions with muons at different times in the proto-neutron star cooling.
We are particularly interested in the late phase ($ t \gtrsim 10 \operatorname{s}$), which will be accessible in the next Galactic supernova but has not been studied well so far.
We calculate both leptonic and semi-leptonic processes, for the latter of which we pay attention also to the form factors with their dependence on the transferred momentum as well as to the modification of the dispersion relations for nucleons on the mean field level.
We find that the flavor-exchange reactions $ν_e + μ^- \rightarrow ν_μ + e^-$ and $\barν_μ + μ^- \rightarrow \barν_e + e^-$ can be dominant, particularly at low energies, over the capture of $ν_e$ on neutron and the scatterings of $\barν_μ$ on nucleons as the opacity sources for these species and that the inverse muon decay $ \barν_e + ν_μ + e^- \leftrightarrows μ^- $ can overwhelm the scatterings of $\barν_e$ and $ν_μ$ on nucleons again at low energies.
At high energies, on the other hand, the corrections in the semi-leptonic processes mentioned above are more important.
We also show the non-trivial energy- and angular dependences of the flavor-exchange reactions and the inverse muon decay.
In the study of the diffusion coefficients from these reactions, we find that $\barν_μ$ is most affected.
These pieces of information are indispensable for numerical computations and the interpretation of results thereof for the proto-neutron star cooling particularly at the very late phase.
△ Less
Submitted 13 November, 2022;
originally announced November 2022.
-
Formation of moons and equatorial ridge around top-shaped asteroids after surface landslide
Authors:
Ryuki Hyodo,
Keisuke Sugiura
Abstract:
Top-shaped asteroids have been observed among near-Earth asteroids. About half of them are reported to have moons (on the order of $\sim 1$wt.\% of the top-shaped primary) and many of them have an equatorial ridge. A recent study has shown that the enigmatic top-shaped figure of asteroids (e.g., Ryugu, Bennu, and Didymos) could result from an axisymmetric landslide of the primary during a fast spi…
▽ More
Top-shaped asteroids have been observed among near-Earth asteroids. About half of them are reported to have moons (on the order of $\sim 1$wt.\% of the top-shaped primary) and many of them have an equatorial ridge. A recent study has shown that the enigmatic top-shaped figure of asteroids (e.g., Ryugu, Bennu, and Didymos) could result from an axisymmetric landslide of the primary during a fast spin-up near the breakup rotation period. Such a landslide would inevitably form a particulate disk around an asteroid with a short timescale ($\sim 3$ hours). However, the long-term full dynamical evolution is not investigated. Here, we perform a continuous simulation ($\sim 700$ hours) that investigates the sequence of events from the surface landslide that forms a top-shaped asteroid and a particulate disk to disk evolution. We show that the disk quickly spreads and produces moons (within $\sim 300$ hours). The mass of the formed moon is consistent with what is observed around the top-shaped asteroids. We also demonstrate that an equatorial ridge is naturally formed because a fraction of the disk particles re-accretes selectively onto the equatorial region of the primary. We envision that Ryugu and Bennu could once have an ancient moon that was later lost due to a successive moon's orbital evolution. Alternatively, at the top-shaped asteroid that has a moon, such as Didymos, no significant orbital evolution of the moon has occurred that would result in its loss. Our study would also be qualitatively applicable to any rubble-pile asteroids near the breakup rotation period.
△ Less
Submitted 15 September, 2022;
originally announced September 2022.
-
Giant Impact onto a Vesta-Like Asteroid and Formation of Mesosiderites through Mixing of Metallic Core and Surface Crust
Authors:
Keisuke Sugiura,
Makiko K. Haba,
Hidenori Genda
Abstract:
Mesosiderites are a type of stony-iron meteorites composed of a mixture of silicates and Fe-Ni metals. The mesosiderite silicates and metals are considered to have originated from the crust and metal core, respectively, of a differentiated asteroid. In contrast, mesosiderites rarely contain the olivine that is mainly included in a mantle. Although a giant impact onto a differentiated asteroid is c…
▽ More
Mesosiderites are a type of stony-iron meteorites composed of a mixture of silicates and Fe-Ni metals. The mesosiderite silicates and metals are considered to have originated from the crust and metal core, respectively, of a differentiated asteroid. In contrast, mesosiderites rarely contain the olivine that is mainly included in a mantle. Although a giant impact onto a differentiated asteroid is considered to be a probable mechanism to mix crust and metal materials to form mesosiderites, it is not obvious how such a giant impact can form mesosiderite-like materials without including mantle materials. We conducted numerical simulations of giant impacts onto differentiated asteroids, using the smoothed particle hydrodynamics method, to investigate the detailed distribution of mixed materials on the resultant bodies. For the internal structure of a target body, we used a thin-crust model derived from the magma ocean crystallization model of the asteroid Vesta and a thick-crust and a large-core model suggested from the proximity observation of Vesta by the Dawn probe. In the simulations with the former model, excavation of the metal core requires nearly catastrophic impacts and mantle is exposed over large surface areas. Thus, stony-iron materials produced on its surface are likely to include mantle materials and it is difficult to produce mesosiderite-like materials. Conversely, in the simulations with the latter model, mantle materials are exposed only at impact sites, even when the impacts excavate the metal core, and the formation of a surface with little mantle material and the formation of mesosiderite-like materials are possible. Therefore, our simulations suggest that an internal structure with a thick crust and a large core is more likely as a mesosiderite parent body rather than the thin-crust internal structure inferred from the conventional magma ocean model.
△ Less
Submitted 14 February, 2022;
originally announced February 2022.
-
Operational solar flare prediction model using Deep Flare Net
Authors:
Naoto Nishizuka,
Yuki Kubo,
Komei Sugiura,
Mitsue Den,
Mamoru Ishii
Abstract:
We developed an operational solar flare prediction model using deep neural networks, named Deep Flare Net (DeFN). DeFN can issue probabilistic forecasts of solar flares in two categories, such as >=M-class and <M-class events or >=C-class and <C-class events, occurring in the next 24 h after observations and the maximum class of flares occurring in the next 24 h. DeFN is set to run every 6 h and h…
▽ More
We developed an operational solar flare prediction model using deep neural networks, named Deep Flare Net (DeFN). DeFN can issue probabilistic forecasts of solar flares in two categories, such as >=M-class and <M-class events or >=C-class and <C-class events, occurring in the next 24 h after observations and the maximum class of flares occurring in the next 24 h. DeFN is set to run every 6 h and has been operated since January 2019. The input database of solar observation images taken by the Solar Dynamic Observatory (SDO) is downloaded from the data archive operated by the Joint Science Operations Center (JSOC) of Stanford University. Active regions are automatically detected from magnetograms, and 79 features are extracted from each region nearly in real time using multiwavelength observation data. Flare labels are attached to the feature database, and then, the database is standardized and input into DeFN for prediction. DeFN was pretrained using the datasets obtained from 2010 to 2015. The model was evaluated with the skill score of the true skill statistics (TSS) and achieved predictions with TSS = 0.80 for >=M-class flares and TSS = 0.63 for >=C-class flares. For comparison, we evaluated the operationally forecast results from January 2019 to June 2020. We found that operational DeFN forecasts achieved TSS = 0.70 (0.84) for >=C-class flares with the probability threshold of 50 (40)%, although there were very few M-class flares during this period and we should continue monitoring the results for a longer time. Here, we adopted a chronological split to divide the database into two for training and testing. The chronological split appears suitable for evaluating operational models. Furthermore, we proposed the use of time-series cross-validation. The procedure achieved TSS = 0.70 for >=M-class flares and 0.59 for >=C-class flares using the datasets obtained from 2010 to 2017.
△ Less
Submitted 1 December, 2021;
originally announced December 2021.
-
SPH Simulations for Shape Deformation of Rubble-Pile Asteroids Through Spinup: The Challenge for Making Top-Shaped Asteroids Ryugu and Bennu
Authors:
Keisuke Sugiura,
Hiroshi Kobayashi,
Sei-ichiro Watanabe,
Hidenori Genda,
Ryuki Hyodo,
Shu-ichiro Inutsuka
Abstract:
Asteroid Ryugu and asteroid Bennu, which were recently visited by spacecraft Hayabusa2 and OSIRIS-REx, respectively, are spinning top-shaped rubble piles. Other axisymmetric top-shaped near-Earth asteroids have been observed with ground-based radar, most of which rotate near breakup rotation periods of ~ 3 hours. This suggests that rotation-induced deformation of asteroids through rotational spinu…
▽ More
Asteroid Ryugu and asteroid Bennu, which were recently visited by spacecraft Hayabusa2 and OSIRIS-REx, respectively, are spinning top-shaped rubble piles. Other axisymmetric top-shaped near-Earth asteroids have been observed with ground-based radar, most of which rotate near breakup rotation periods of ~ 3 hours. This suggests that rotation-induced deformation of asteroids through rotational spinup produces top shapes. Although some previous simulations using the Discrete Element Method showed that spinup of rubble piles may produce oblate top shapes, it is still unclear what kinds of conditions such as friction angles of constituent materials and spinup timescales are required for top-shape formation. Here we show, through Smoothed Particle Hydrodynamics simulations of granular bodies spinning-up at different rates, that the rotation-induced deformation of spherical rubble piles before breakup can be classified into three modes according to the friction angle φ_{d}: quasi-static and internal deformation for φ_{d} < 40 degrees, dynamical and internal deformation for 50 degrees < φ_{d} < 60 degrees, and surface landslides for φ_{d} > 70 degrees. Note that these apparent large values of friction angle can be acceptable if we consider the effect of cohesion among blocks of a rubble pile under weak gravity. Bodies with φ_{d} < 60 degrees evolve into oblate spheroids through internal deformation, but never form pronounced equators defining a top shape. In contrast, bodies with φ_{d} > 70 degrees deform into axisymmetric top shapes through an axisymmetric surface landslides if spinup timescales are < a few days. In addition, through slow spinups with timescales > 1 month, bodies with φ_{d} > 70 degrees deform into non-axisymmetric shapes via localized landslides. We suggest that rapid spinup mechanisms are preferable for the formation of axisymmetric top shapes.
△ Less
Submitted 27 April, 2021;
originally announced April 2021.
-
Reliable Probability Forecast of Solar Flares: Deep Flare Net-Reliable (DeFN-R)
Authors:
Naoto Nishizuka,
Yûki Kubo,
Komei Sugiura,
Mitsue Den,
Mamoru Ishii
Abstract:
We developed a reliable probabilistic solar flare forecasting model using a deep neural network, named Deep Flare Net-Reliable (DeFN-R). The model can predict the maximum classes of flares that occur in the following 24 h after observing images, along with the event occurrence probability. We detected active regions from 3x10^5 solar images taken during 2010-2015 by Solar Dynamic Observatory and e…
▽ More
We developed a reliable probabilistic solar flare forecasting model using a deep neural network, named Deep Flare Net-Reliable (DeFN-R). The model can predict the maximum classes of flares that occur in the following 24 h after observing images, along with the event occurrence probability. We detected active regions from 3x10^5 solar images taken during 2010-2015 by Solar Dynamic Observatory and extracted 79 features for each region, which we annotated with flare occurrence labels of X-, M-, and C-classes. The extracted features are the same as used by Nishizuka et al. (2018); for example, line-of-sight/vector magnetograms in the photosphere, brightening in the corona, and the X-ray emissivity 1 and 2 h before an image. We adopted a chronological split of the database into two for training and testing in an operational setting: the dataset in 2010-2014 for training and the one in 2015 for testing. DeFN-R is composed of multilayer perceptrons formed by batch normalizations and skip connections. By tuning optimization methods, DeFN-R was trained to optimize the Brier skill score (BSS). As a result, we achieved BSS = 0.41 for >=C-class flare predictions and 0.30 for >=M-class flare predictions by improving the reliability diagram while keeping the relative operating characteristic curve almost the same. Note that DeFN is optimized for deterministic prediction, which is determined with a normalized threshold of 50%. On the other hand, DeFN-R is optimized for a probability forecast based on the observation event rate, whose probability threshold can be selected according to users' purposes.
△ Less
Submitted 6 July, 2020;
originally announced July 2020.
-
High-Resolution Simulations of Catastrophic Disruptions: Resultant Shape Distributions
Authors:
Keisuke Sugiura,
Hiroshi Kobayashi,
Shu-ichiro Inutsuka
Abstract:
The members of asteroid families have various shapes. We investigate the origin of their shapes by high-resolution impact simulations for catastrophic disruptions using a Smoothed Particle Hydrodynamics code. Collisional remnants produced through our simulations of the catastrophic disruptions mainly have spherical or bilobed shapes. However, no flat remnants with the ratio of minor to major axis…
▽ More
The members of asteroid families have various shapes. We investigate the origin of their shapes by high-resolution impact simulations for catastrophic disruptions using a Smoothed Particle Hydrodynamics code. Collisional remnants produced through our simulations of the catastrophic disruptions mainly have spherical or bilobed shapes. However, no flat remnants with the ratio of minor to major axis lengths c/a < 0.5 are formed. The results of the simulations provide various shapes of asteroids and explain most of the shapes in asteroid families that are supposed to be produced through catastrophic disruptions. However, the present simulations do not explain significantly flat asteroids. We suggest that these flat asteroids may be interlopers or formed through low-velocity collisions between member asteroids.
△ Less
Submitted 21 November, 2019; v1 submitted 19 November, 2019;
originally announced November 2019.
-
Collisional Elongation: Possible Origin of Extremely Elongated Shape of 1I/`Oumuamua
Authors:
Keisuke Sugiura,
Hiroshi Kobayashi,
Shu-ichiro Inutsuka
Abstract:
Light curve observations of a recently discovered interstellar object 1I/`Oumuamua suggest that this object has an extremely elongated shape with the axis ratio 0.3 or smaller. Planetesimal collisions can produce irregular shapes including elongated shapes. In this paper, we suggest that the extremely elongated shape of 1I/`Oumuamua may be the result of such an impact. To find detailed impact cond…
▽ More
Light curve observations of a recently discovered interstellar object 1I/`Oumuamua suggest that this object has an extremely elongated shape with the axis ratio 0.3 or smaller. Planetesimal collisions can produce irregular shapes including elongated shapes. In this paper, we suggest that the extremely elongated shape of 1I/`Oumuamua may be the result of such an impact. To find detailed impact conditions to form the extremely elongated objects, we conduct numerical simulations of planetesimal collisions using Smoothed Particle Hydrodynamics method for elastic dynamics with self-gravity and interparticle friction. Impacts into strengthless target planetesimals with radius 50 m are conducted with various ratios of impactor mass to target mass q, friction angles phi_d, impact velocities v_imp, and impact angles theta_imp. We find that impacts with q \geq 0.5, phi_d \geq 40 degrees, v_imp \leq 40 degrees, and theta_imp \leq 30 degrees produce remnants with the ratio of intermediate to major axis length less than 0.3. This impact condition suggests that the parent protoplanetary disk in the planetesimal collision stage was weakly turbulent (alpha < 10^{-4} for the inner disk) and composed of planetesimals smaller than ~ 7 km to ensure small impact velocity.
△ Less
Submitted 15 March, 2019;
originally announced March 2019.
-
Linear Analysis of the Shock Instability in Core-collapse Supernovae: Influences of Acoustic Power and Fluctuations of Neutrino Luminosity
Authors:
Ken'ichi Sugiura,
Kazuya Takahashi,
Shoichi Yamada
Abstract:
This paper is a sequel to Takahashi et al. (2016), in which the authors investigated the influences of fluctuations in pre-shock accreting matter on the linear stability of the standing accretion shock in core-collapse supernovae (CCSNe). If one understands that this was concerning the effect of the outer boundary condition for the post-shock accretion flows, the present paper should be regarded a…
▽ More
This paper is a sequel to Takahashi et al. (2016), in which the authors investigated the influences of fluctuations in pre-shock accreting matter on the linear stability of the standing accretion shock in core-collapse supernovae (CCSNe). If one understands that this was concerning the effect of the outer boundary condition for the post-shock accretion flows, the present paper should be regarded as an investigation on the possible influences of the inner boundary conditions. More specifically, we impose a time-dependent, oscillating condition at the inner boundary, having in mind the injection of acoustic power by an oscillating proto-neutron star. We also consider possible correlations between the inner and outer boundary conditions as invoked in the argument for Lepton-number Emission Self-sustained Asymmetry, or LESA. In this paper, we conduct the linear stability analysis of the standing accretion shock commonly encountered in CCSNe based on Laplace transform. We find that the acoustic power enhances the standing accretion shock instability, or SASI, especially when the luminosity is low. On the other hand, the correlation between the fluctuations of neutrino luminosity at the neutrino sphere has little influences on the instability, changing the amplitudes of eigenmodes only slightly. We further investigate steady solution of perturbation equations, being motivated by LESA, and conclude that not the difference but the sum of the fluxes of electron-type neutrinos and anti-neutrino is the key ingredient to production of the self-sustained steady perturbed configuration.
△ Less
Submitted 1 March, 2019;
originally announced March 2019.
-
Deep Flare Net (DeFN) model for solar flare prediction
Authors:
Naoto Nishizuka,
Komei Sugiura,
Yuki Kubo,
Mitsue Den,
Mamoru Ishii
Abstract:
We developed a solar flare prediction model using a deep neural network (DNN), named Deep Flare Net (DeFN). The model can calculate the probability of flares occurring in the following 24 h in each active region, which is used to determine the most likely maximum classes of flares via a binary classification (e.g., >=M class versus <M class or >=C class versus <C class). From 3x10^5 observation im…
▽ More
We developed a solar flare prediction model using a deep neural network (DNN), named Deep Flare Net (DeFN). The model can calculate the probability of flares occurring in the following 24 h in each active region, which is used to determine the most likely maximum classes of flares via a binary classification (e.g., >=M class versus <M class or >=C class versus <C class). From 3x10^5 observation images taken during 2010-2015 by Solar Dynamic Observatory, we automatically detected sunspots and calculated 79 features for each region, to which flare occurrence labels of X-, M-, and C-class were attached. We adopted the features used in Nishizuka et al. (2017) and added some features for operational prediction: coronal hot brightening at 131 A (T>=10^7 K) and the X-ray and 131 A intensity data 1 and 2 h before an image. For operational evaluation, we divided the database into two for training and testing: the dataset in 2010-2014 for training and the one in 2015 for testing. The DeFN model consists of deep multilayer neural networks, formed by adapting skip connections and batch normalizations. To statistically predict flares, the DeFN model was trained to optimize the skill score, i.e., the true skill statistic (TSS). As a result, we succeeded in predicting flares with TSS=0.80 for >=M-class flares and TSS=0.63 for >=C-class flares. Note that in usual DNN models, the prediction process is a black box. However, in the DeFN model, the features are manually selected, and it is possible to analyze which features are effective for prediction after evaluation.
△ Less
Submitted 9 May, 2018;
originally announced May 2018.
-
Toward Understanding the Origin of Asteroid Geometries: Variety in Shapes Produced by Equal-Mass Impacts
Authors:
Keisuke Sugiura,
Hiroshi Kobayashi,
Shu-ichiro Inutsuka
Abstract:
More than a half of asteroids in the main belt have irregular shapes with the ratios of the minor to major axis lengths less than 0.6. One of the mechanisms to create such shapes is collisions between asteroids. The relationship between shapes of collisional outcomes and impact conditions such as impact velocities may provide information on the collisional environments and its evolutionary stages…
▽ More
More than a half of asteroids in the main belt have irregular shapes with the ratios of the minor to major axis lengths less than 0.6. One of the mechanisms to create such shapes is collisions between asteroids. The relationship between shapes of collisional outcomes and impact conditions such as impact velocities may provide information on the collisional environments and its evolutionary stages when those asteroids are created. In this study, we perform numerical simulations of collisional destruction of asteroids with radii 50 km and subsequent gravitational reaccumulation using Smoothed Particle Hydrodynamics for elastic dynamics with self-gravity, a model of fracture of rock, and a model of friction of completely damaged rock. We systematically vary the impact velocity from 50 m/s to 400 m/s and the impact angle from 5 degrees to 45 degrees. We investigate shapes of the largest remnants resulting from collisional simulations. As a result, various shapes (bilobed, spherical, flat, elongated, and hemispherical shapes) are formed through equal-mass and low-velocity (50 - 400 m/s) impacts. We clarify a range of the impact angle and velocity to form each shape. Our results indicate that irregular shapes, especially flat shapes, of asteroids with diameters larger than 80 km are likely to be formed through similar-mass and low-velocity impacts, which are likely to occur in the primordial environment prior to the formation of Jupiter.
△ Less
Submitted 10 October, 2018; v1 submitted 29 April, 2018;
originally announced April 2018.
-
Solar Flare Prediction Model with Three Machine-Learning Algorithms Using Ultraviolet Brightening and Vector Magnetogram
Authors:
N. Nishizuka,
K. Sugiura,
Y. Kubo,
M. Den,
S. Watari,
M. Ishii
Abstract:
We developed a flare prediction model using machine learning, which is optimized to predict the maximum class of flares occurring in the following 24 h. Machine learning is used to devise algorithms that can learn from and make decisions on a huge amount of data. We used solar observation data during the period 2010-2015, such as vector magnetogram, ultraviolet (UV) emission, and soft X-ray emissi…
▽ More
We developed a flare prediction model using machine learning, which is optimized to predict the maximum class of flares occurring in the following 24 h. Machine learning is used to devise algorithms that can learn from and make decisions on a huge amount of data. We used solar observation data during the period 2010-2015, such as vector magnetogram, ultraviolet (UV) emission, and soft X-ray emission taken by the Solar Dynamics Observatory and the Geostationary Operational Environmental Satellite. We detected active regions from the full-disk magnetogram, from which 60 features were extracted with their time differentials, including magnetic neutral lines, the current helicity, the UV brightening, and the flare history. After standardizing the feature database, we fully shuffled and randomly separated it into two for training and testing. To investigate which algorithm is best for flare prediction, we compared three machine learning algorithms: the support vector machine (SVM), k-nearest neighbors (k-NN), and extremely randomized trees (ERT). The prediction score, the true skill statistic (TSS), was higher than 0.9 with a fully shuffled dataset, which is higher than that for human forecasts. It was found that k-NN has the highest performance among the three algorithms. The ranking of the feature importance showed that the previous flare activity is most effective, followed by the length of magnetic neutral lines, the unsigned magnetic flux, the area of UV brightening, and the time differentials of features over 24 h, all of which are strongly correlated with the flux emergence dynamics in an active region.
△ Less
Submitted 6 November, 2016;
originally announced November 2016.
-
An Extension of Godunov SPH II: Application to Elastic Dynamics
Authors:
Keisuke Sugiura,
Shu-ichiro Inutsuka
Abstract:
Godunov Smoothed Particle Hydrodynamics (Godunov SPH) method is a computational fluid dynamics method that utilizes a Riemann solver and achieves the second-order accuracy in space. In this paper, we extend the Godunov SPH method to elastic dynamics by incorporating deviatoric stress tensor that represents the stress for shear deformation or anisotropic compression. Analogously to the formulation…
▽ More
Godunov Smoothed Particle Hydrodynamics (Godunov SPH) method is a computational fluid dynamics method that utilizes a Riemann solver and achieves the second-order accuracy in space. In this paper, we extend the Godunov SPH method to elastic dynamics by incorporating deviatoric stress tensor that represents the stress for shear deformation or anisotropic compression. Analogously to the formulation of the original Godunov SPH method, we formulate the equation of motion, the equation of energy, and the time evolution equation of deviatoric stress tensor so that the resulting discretized system achieves the second-order accuracy in space. The standard SPH method tends to suffer from the tensile instability that results in unphysical clustering of particles especially in tension-dominated region. We find that the tensile instability can be suppressed by selecting appropriate interpolation for density distribution in the equation of motion for the Godunov SPH method even in the case of elastic dynamics. Several test calculations for elastic dynamics are performed, and the accuracy and versatility of the present method are shown.
△ Less
Submitted 13 December, 2016; v1 submitted 24 May, 2016;
originally announced May 2016.
-
An Extension of Godunov SPH: Application to Negative Pressure Media
Authors:
Keisuke Sugiura,
Shu-ichiro Inutsuka
Abstract:
The modification of Smoothed Particle Hydrodynamics (SPH) method with Riemann Solver is called Godunov SPH. We further extend the Godunov SPH to the description of a medium with negative pressure. Under certain circumstances, the SPH method shows an unphysical instability that results in particle clustering. This instability is called the tensile instability. The tensile instability occurs in posi…
▽ More
The modification of Smoothed Particle Hydrodynamics (SPH) method with Riemann Solver is called Godunov SPH. We further extend the Godunov SPH to the description of a medium with negative pressure. Under certain circumstances, the SPH method shows an unphysical instability that results in particle clustering. This instability is called the tensile instability. The tensile instability occurs in positive pressure regions in a regular fluid if a very large number of neighbor particles are used with certain shapes of kernel functions, and it is significant in negative pressure regions that emerge in stretched elastic bodies. We must suppress the tensile instability in SPH for calculations of elastic bodies. In this study, we develop a new technique to remove the tensile instability by extending the Godunov SPH method and conducting a linear stability analysis of the equation of motion for the extended method. We find that the tensile instability can be suppressed by choosing an appropriate order of interpolation in the equation of motion of the Godunov SPH method. We also derive an analytic solution for a Riemann solver for a simple equation of state of an elastic body, and construct a Godunov SPH method for the equation of state that allows negative pressure.
△ Less
Submitted 14 December, 2015; v1 submitted 19 May, 2015;
originally announced May 2015.