-
Modeling Solids in Nuclear Astrophysics with Smoothed Particle Hydrodynamics
Authors:
Irina Sagert,
Oleg Korobkin,
Ingo Tews,
Bing-Jyun Tsao,
Hyun Lim,
Michael J. Falato,
Julien Loiseau
Abstract:
Smoothed Particle Hydrodynamics (SPH) is a frequently applied tool in computational astrophysics to solve the fluid dynamics equations governing the systems under study. For some problems, for example when involving asteroids and asteroid impacts, the additional inclusion of material strength is necessary in order to accurately describe the dynamics. In compact stars, that is white dwarfs and neut…
▽ More
Smoothed Particle Hydrodynamics (SPH) is a frequently applied tool in computational astrophysics to solve the fluid dynamics equations governing the systems under study. For some problems, for example when involving asteroids and asteroid impacts, the additional inclusion of material strength is necessary in order to accurately describe the dynamics. In compact stars, that is white dwarfs and neutron stars, solid components are also present. Neutron stars have a solid crust which is the strongest material known in nature. However, their dynamical evolution, when modeled via SPH or other computational fluid dynamics codes, is usually described as a purely fluid dynamics problem. Here, we present the first 3D simulations of neutron-star crustal toroidal oscillations including material strength with the Los Alamos National Laboratory SPH code FleCSPH. In the first half of the paper, we present the numerical implementation of solid material modeling together with standard tests. The second half is on the simulation of crustal oscillations in the fundamental toroidal mode. Here, we dedicate a large fraction of the paper to approaches which can suppress numerical noise in the solid. If not minimized, the latter can dominate the crustal motion in the simulations.
△ Less
Submitted 17 November, 2023; v1 submitted 3 November, 2022;
originally announced November 2022.
-
An Explosively Driven Launcher Capable of 10 km/s Projectile Velocities
Authors:
Justin Huneault,
Jason Loiseau,
Myles T. Hildebrand,
Andrew J. Higgins
Abstract:
Launching large (> 1 g) well-characterized projectiles to velocities beyond 10 km/s is of interest for a number of scientific fields but is beyond the reach of current hypervelocity launcher technology. This paper reports the development of an explosively driven light-gas gun that has demonstrated the ability to launch 8 mm diameter, 0.36 g magnesium projectiles to 10.4 km/s. The implosion-driven…
▽ More
Launching large (> 1 g) well-characterized projectiles to velocities beyond 10 km/s is of interest for a number of scientific fields but is beyond the reach of current hypervelocity launcher technology. This paper reports the development of an explosively driven light-gas gun that has demonstrated the ability to launch 8 mm diameter, 0.36 g magnesium projectiles to 10.4 km/s. The implosion-driven launcher (IDL) uses the linear implosion of a pressurized tube to shock-compress helium gas to a pressure of \SI{5}{\giga\pascal}, which then expands to propel a projectile to hypervelocity. The launch cycle of the IDL is explored with the use of down-bore velocimetry experiments and a quasi-one-dimensional internal ballistics solver. A detailed overview of the design of the 8 mm launcher is presented, with an emphasis on the unique considerations which arise from the explosively driven propellant compression and the resulting extreme pressures and temperatures. The high average driving pressure results in a launcher that is compact, with a total length typically less than a meter. The possibility to scale the design to larger projectile sizes (25 mm diameter) is demonstrated. Finally, concepts for a modified launch cycle which may allow the IDL to reach significantly greater projectile velocities are explored conceptually and with preliminary experiments.
△ Less
Submitted 29 June, 2022;
originally announced June 2022.
-
Realistic Kilonova Up Close
Authors:
Alexandra Ruth Stewart,
Li-Ta Lo,
Oleg Korobkin,
Irina Sagert,
Julien Loiseau,
Hyun Lim,
Mark Alexander Kaltenborn,
Christopher Michael Mauney,
Jonah Maxwell Miller
Abstract:
Neutron star mergers are cosmic catastrophes that produce some of the most energetic observed phenomena: short gamma-ray bursts, gravitational wave signals, and kilonovae. The latter are optical transients, powered by radioactive nuclides which are synthesized when the neutron-rich ejecta of a disrupted neutron star undergoes decompression. We model this decompression phase using data from simulat…
▽ More
Neutron star mergers are cosmic catastrophes that produce some of the most energetic observed phenomena: short gamma-ray bursts, gravitational wave signals, and kilonovae. The latter are optical transients, powered by radioactive nuclides which are synthesized when the neutron-rich ejecta of a disrupted neutron star undergoes decompression. We model this decompression phase using data from simulations of post-merger accretion disk winds. We use smoothed particle hydrodynamics with realistic nuclear heating to model the expansion over multiple scales, from initially several thousand km to billions of km. We then render a realistic image of a kilonova ejecta as it would appear for a nearby observer. This is the first time such a visualization is performed using input from state-of-the-art accretion disk simulations, nuclear physics and atomic physics. The volume rendering in our model computes an opacity transfer function on the basis of the physical opacity, varying significantly with the inhomogeneity of the neutron richness in the ejecta. Other physical quantities such as temperature or electron fraction can be visualized using an independent color transfer function. We discuss several difficulties with the ParaView application that we encountered during the visualization process, and give descriptions of our solutions and workarounds which could be used for future improvements.
△ Less
Submitted 5 January, 2022;
originally announced January 2022.
-
Modeling Neutron Star Oscillations in a Fixed General Relativistic Background Including Solid Crust Dynamics
Authors:
Bing-Jyun Tsao,
Irina Sagert,
Oleg Korobkin,
Ingo Tews,
Hyun Lim,
Gary Dilts,
Julien Loiseau
Abstract:
Measurements of the gravitational-wave signals from neutron star mergers allow scientists to learn about the interior of neutron stars and the properties of dense nuclear matter. The study of neutron star mergers is usually performed with computational fluid dynamics codes, mostly in Eulerian but also in Lagrangian formulation such as smoothed particle hydrodynamics (SPH). Codes include our best k…
▽ More
Measurements of the gravitational-wave signals from neutron star mergers allow scientists to learn about the interior of neutron stars and the properties of dense nuclear matter. The study of neutron star mergers is usually performed with computational fluid dynamics codes, mostly in Eulerian but also in Lagrangian formulation such as smoothed particle hydrodynamics (SPH). Codes include our best knowledge of nuclear matter in the form of an equation of state as well as effects of general relativity (GR). However, one important aspect of neutron stars is usually ignored: the solid nature of their crust. The solid matter in the crust is the strongest material known in nature which could lead to a multitude of possible observational effects that have not been studied in dynamical simulations yet. The crust could change the way a neutron star deforms during a merger, leaving an imprint in the gravitational wave signal. It could even shatter during the inspiral, producing a potentially observable electromagnetic signal. Here, we present a first study of the dynamical behavior of neutron stars with a solid crust and fixed GR background with FleCSPH. FleCSPH is a general-purpose SPH code, developed at Los Alamos National Laboratory. It features an efficient algorithm for gravitational interactions via the Fast Multipole Method, which, together with the implemented nuclear equation of state, makes it appropriate for astrophysical applications. The solid material dynamics is described via the elastic-perfectly plastic model with maximum-strain breaking. Despite its simplicity, the model reproduces the stress-strain curve of crustal material as extracted from microphysical simulations very well. We present first tests of our implementation via simulations of neutron star oscillations and give an outlook on our study of the dynamical behavior of the solid crust in neutron star merger events.
△ Less
Submitted 15 July, 2021;
originally announced July 2021.
-
Conservation of Angular Momentum in the Fast Multipole Method
Authors:
Oleg Korobkin,
Hyun Lim,
Irina Sagert,
Julien Loiseau,
Christopher Mauney,
M. Alexander R. Kaltenborn,
Bing-Jyun Tsao,
Wesley P. Even
Abstract:
Smoothed particle hydrodynamics (SPH) is positioned as having ideal conservation properties. When properly implemented, conservation of total mass, energy, and both linear and angular momentum is guaranteed exactly, up to machine precision. This is particularly important for some applications in computational astrophysics, such as binary dynamics, mergers, and accretion of compact objects (neutron…
▽ More
Smoothed particle hydrodynamics (SPH) is positioned as having ideal conservation properties. When properly implemented, conservation of total mass, energy, and both linear and angular momentum is guaranteed exactly, up to machine precision. This is particularly important for some applications in computational astrophysics, such as binary dynamics, mergers, and accretion of compact objects (neutron stars, black holes, and white dwarfs). However, in astrophysical applications that require the inclusion of gravity, calculating pairwise particle interactions becomes prohibitively expensive. In the Fast Multipole Method (FMM), they are, therefore, replaced with symmetric interactions between distant clusters of particles (contained in the tree nodes) Although such an algorithm is linear momentum-conserving, it introduces spurious torques that violate conservation of angular momentum. We present a modification of FMM that is free of spurious torques and conserves angular momentum explicitly. The new method has practically no computational overhead compared to the standard FMM.
△ Less
Submitted 15 July, 2021;
originally announced July 2021.
-
FleCSPH: The Next Generation FleCSIble Parallel Computational Infrastructure for Smoothed Particle Hydrodynamics
Authors:
Julien Loiseau,
Hyun Lim,
Mark Alexander Kaltenborn,
Oleg Korobkin,
Christopher M. Mauney,
Irina Sagert,
Wesley P. Even,
Benjamin K. Bergen
Abstract:
FleCSPH is a smoothed particle hydrodynamics simulation tool, based on the compile-time configurable framework FleCSI. The asynchronous distributed tree topology combined with a fast multipole method allows FleCSPH to efficiently compute hydrodynamics and long range particle-particle interactions. FleCSPH provides initial data generators, particle relaxation techniques, and standard evolution driv…
▽ More
FleCSPH is a smoothed particle hydrodynamics simulation tool, based on the compile-time configurable framework FleCSI. The asynchronous distributed tree topology combined with a fast multipole method allows FleCSPH to efficiently compute hydrodynamics and long range particle-particle interactions. FleCSPH provides initial data generators, particle relaxation techniques, and standard evolution drivers, which can be easily modified and extended to user-specific setups. Data input/output uses the H5part format, compatible with modern visualization software.
△ Less
Submitted 6 July, 2020;
originally announced July 2020.