-
Do not forget the electrons: Extending moderately-sized nuclear networks for multidimensional hydrodynamic codes
Authors:
Domingo García-Senz,
Rubén M. Cabezón,
Moritz Reichert,
Axel S. Lechuga,
José A. Escartín,
Athanasios Psaltis,
Almudena Arcones,
Friedrich-Karl Thielemann
Abstract:
We present here an extended nuclear network, with 90 species, designed for being coupled with hydrodynamic simulations, which includes neutrons, protons, electrons, positrons, and the corresponding neutrino and anti-neutrino emission. This network is also coupled with temperature, making it extremely robust and, together with its size, unique of its kind. The inclusion of electron captures on free…
▽ More
We present here an extended nuclear network, with 90 species, designed for being coupled with hydrodynamic simulations, which includes neutrons, protons, electrons, positrons, and the corresponding neutrino and anti-neutrino emission. This network is also coupled with temperature, making it extremely robust and, together with its size, unique of its kind. The inclusion of electron captures on free protons makes the network very appropriate for multidimensional studies of Type Ia supernova explosions, especially when the exploding object is a massive white dwarf.
The results obtained with the proposed medium-sized network compare fairly well, to a few percent, with those computed with the extended network WinNet (> 2000 isotopes) in scenarios reproducing the gross physical conditions of current Type Ia supernova explosion models. In those cases where the carbon and oxygen fuel ignites at high density, the high-temperature plateau typical of the nuclear statistical equilibrium regime is well defined and stable, allowing large integration time steps. We show that the inclusion of electron captures on free protons substantially improves the estimation of the electron fraction of the mixture. Therefore, the pressure is better determined than in networks where electron captures are excluded, which will ultimately lead to more reliable hydrodynamic models. Explosive combustion of helium at low density, occurring near the surface layer of a white dwarf, is also better described with the proposed network, which gives nuclear energy generation rates much closer to WinNet than typical reduced alpha networks.
△ Less
Submitted 11 July, 2024; v1 submitted 6 March, 2024;
originally announced March 2024.
-
Mixing Sinc kernels to improve interpolations in smoothed particle hydrodynamics without pairing instability
Authors:
Rubén M. Cabezón,
Domingo García-Senz
Abstract:
The smoothed particle hydrodynamic technique is strongly based on the proper choice of interpolation functions. This statement is particularly relevant for the study of subsonic fluxes and turbulence, where inherent small errors in the averaging procedures introduce excessive damping on the smallest scales. To mitigate these errors we can increase both the number of interpolating points and the or…
▽ More
The smoothed particle hydrodynamic technique is strongly based on the proper choice of interpolation functions. This statement is particularly relevant for the study of subsonic fluxes and turbulence, where inherent small errors in the averaging procedures introduce excessive damping on the smallest scales. To mitigate these errors we can increase both the number of interpolating points and the order of the interpolating kernel function. However, this approach leads to a higher computational burden across all fluid regions. Ideally, the development of a single kernel function capable of effectively accommodating varying numbers of interpolating points in different fluid regions, providing good resolution and minimal errors would be highly desirable. In this work, we revisit and extend the main properties of a family of interpolators called $Sinc~kernels$ and compare them with the widely used family of Wendland kernels. We show that a linear combination of low- and high-order Sinc kernels generates good-quality interpolators, which are resistant to pairing instability while maintaining good sampling properties in a wide range of neighbor interpolating points, $60\le n_b\le 400$. We show that a particular case of this linear mix of Sincs produces a well-balanced and robust kernel that improves previous results in the Gresho-Chan vortex experiment even when the number of neighbors is not large, while yielding a good convergence rate. Although such a mixing technique is ideally suited for Sinc kernels owing to their excellent flexibility, it can be easily applied to other interpolating families such as the B-splines and Wendland kernels.
△ Less
Submitted 13 February, 2024; v1 submitted 20 July, 2023;
originally announced July 2023.
-
Axisymmetric smoothed particle magneto-hydrodynamics
Authors:
Domingo García-Senz,
Robert Wissing,
Rubén M. Cabezón,
Eda Vurgun,
Manuel Linares
Abstract:
Many astrophysical and terrestrial scenarios involving magnetic fields can be approached in axial geometry. Although the smoothed particle hydrodynamics (SPH) technique has been successfully extended to magneto-hydrodynamics (MHD), a well-verified, axisymmetric MHD scheme based on such technique does not exist yet. In this work we fill that gap in the scientific literature and propose and check a…
▽ More
Many astrophysical and terrestrial scenarios involving magnetic fields can be approached in axial geometry. Although the smoothed particle hydrodynamics (SPH) technique has been successfully extended to magneto-hydrodynamics (MHD), a well-verified, axisymmetric MHD scheme based on such technique does not exist yet. In this work we fill that gap in the scientific literature and propose and check a novel axisymmetric MHD hydrodynamic code, that can be applied to physical problems which display the adequate geometry. We show that the hydrodynamic code built following these axisymmetric hypothesis is able to produce similar results than standard 3D-SPMHD codes with equivalent resolution but with much lesser computational load.
△ Less
Submitted 11 November, 2022;
originally announced November 2022.
-
Sub-stellar engulfment by a main sequence star: where is the lithium?
Authors:
Rubén M. Cabezón,
Carlos Abia,
Inma Domínguez,
Domingo García-Senz
Abstract:
In this work, we study whether the engulfment of a brown dwarf (BD) by a solar-like main-sequence (MS) star can significantly alter the structure of the star and the Li content on its surface. We perform 3D Smoothed Particle Hydrodynamics simulations of the engulfment of a BD with masses 0.01 and 0.019 Msun, onto an MS star of 1 Msun and solar composition, in three different scenarios: a head-on c…
▽ More
In this work, we study whether the engulfment of a brown dwarf (BD) by a solar-like main-sequence (MS) star can significantly alter the structure of the star and the Li content on its surface. We perform 3D Smoothed Particle Hydrodynamics simulations of the engulfment of a BD with masses 0.01 and 0.019 Msun, onto an MS star of 1 Msun and solar composition, in three different scenarios: a head-on collision, a grazing collision, and a merger. We study the dynamics of the interaction in detail, and the relevance of the type of interaction and the mass of the BD on the final fate of the sub-stellar object and the host star in terms of mass loss of the system, angular momentum transfer, and changes in the Li abundance in the surface of the host star. We found that most of the BD mass is diluted in the denser region of the MS star. Only in the merger scenario a significant fraction (40%) of the BD material would remain in the outer layers. We find a clear increase in the surface rotational velocity of the host star after the interaction, ranging between 25 km/s (grazing collision) to 50 km/s (merger). We also find a significant mass loss from the system (1e-4 - 1e-3 Msun) due to the engulfment, which in the case of the merger, may form a circumstellar disk-like structure. Assuming that neither the depth of the convective envelope of the host star nor its mass content are modified during the interaction, a small change in the surface Li abundance in the head-on and grazing collisions is found. However, in the merger we find large Li enhancements, by factors 20-30, depending on the BD mass. Some of these features could be detected observationally in the host star provide they remain long enough time.
△ Less
Submitted 2 November, 2022; v1 submitted 31 October, 2022;
originally announced October 2022.
-
Axisymmetric magneto-hydrodynamics with SPH
Authors:
Domingo García-Senz,
Robert Wissing,
Rubén M. Cabezón
Abstract:
Many interesting terrestrial and astrophysical scenarios involving magnetic fields can be approached in axial geometry. Even though the Lagrangian smoothed particle hydrodynamics (SPH) technique has been successfully extended to handle magneto-hydrodynamic (MHD) problems, a well-verified, axisymmetric MHD scheme based on the SPH technique does not exist. In this work, we propose and check a new ax…
▽ More
Many interesting terrestrial and astrophysical scenarios involving magnetic fields can be approached in axial geometry. Even though the Lagrangian smoothed particle hydrodynamics (SPH) technique has been successfully extended to handle magneto-hydrodynamic (MHD) problems, a well-verified, axisymmetric MHD scheme based on the SPH technique does not exist. In this work, we propose and check a new axisymmetric MHD hydrodynamic code that can be applied to astrophysical and engineering problems which display an adequate geometry. We show that a hydrodynamic code built on these axisymmetric premises is able to produce similar results to standard 3D-SPHMHD codes but with much lesser computational effort.
△ Less
Submitted 3 July, 2022; v1 submitted 10 June, 2022;
originally announced June 2022.
-
Conservative, density-based smoothed particle hydrodynamics with improved partition of the unity and better estimation of gradients
Authors:
Domingo García-Senz,
Rubén M. Cabezón,
José A. Escartín
Abstract:
The correct evaluation of gradients is at the cornerstone of the smoothed particle hydrodynamics (SPH) technique. Using an integral approach to estimate gradients has proven to enhance accuracy substantially. Such approach retains the Lagrangian structure of SPH equations and is fully conservative. But, in practice, it is difficult to make the Lagrangian formulation totally consistent to an exact…
▽ More
The correct evaluation of gradients is at the cornerstone of the smoothed particle hydrodynamics (SPH) technique. Using an integral approach to estimate gradients has proven to enhance accuracy substantially. Such approach retains the Lagrangian structure of SPH equations and is fully conservative. But, in practice, it is difficult to make the Lagrangian formulation totally consistent to an exact partition of the unity.
In this paper we study, among other things, the connection between the choice of the volume elements (VEs), which enters in the SPH summations, and the accuracy in the gradient estimation within the integral approach scheme (ISPH). A new variant of VEs are proposed which improve the partition of the unity and is fully compatible with the Lagrangian formulation of SPH, including the grad-h corrections. Using analytic considerations, simple static toy models in 1D, and a few full 3D test cases, we show that any improvement in the partition of the unity also leads to a better calculation of gradients when the integral approach is used jointly. Additionally, we propose an easy-to-implement modification of the ISPH scheme, which makes it more flexible and better suited to handle sharp density contrasts.
The ISPH code built with the proposed scheme has been validated with a good number of standard tests, some of them involving contact discontinuities. The performance of the code was excellent in all of them, showing that an improvement in the partition of the unity is not detrimental of the good conservation of energy, momentum, and entropy typical of Lagrangian schemes.
△ Less
Submitted 30 March, 2022; v1 submitted 18 January, 2021;
originally announced January 2021.
-
Self-gravitating barotropic equilibrium configurations of rotating bodies with SPH
Authors:
Domingo García-Senz,
Rubén M. Cabezón,
José M. Blanco Iglesias,
Pablo Lorén-Aguilar
Abstract:
We present a novel relaxation method to build three-dimensional rotating structures of barotropic bodies using the SPH technique. The method is able to relax gaseous structures in rigid as well as differential rotation. The relaxation procedure strongly relies on the excellent conservation of angular momentum that characterizes the SPH technique. The method has been successfully applied to a varie…
▽ More
We present a novel relaxation method to build three-dimensional rotating structures of barotropic bodies using the SPH technique. The method is able to relax gaseous structures in rigid as well as differential rotation. The relaxation procedure strongly relies on the excellent conservation of angular momentum that characterizes the SPH technique. The method has been successfully applied to a variety of zero-temperature white dwarfs and polytropic self-gravitating structures. Our SPH results have been validated by comparing the main features (energies, central densities and the polar to equatorial radius ratio) to those obtained with independent, albeit grid-based methods, as for example, the self-consistent field method, showing that both methods agree within few percents.
△ Less
Submitted 21 March, 2020;
originally announced March 2020.
-
3D Hydrodynamical Simulations of a Brown Dwarf Accretion by a Main-Sequence Star and its Impact on the Surface Li Abundance
Authors:
C. Abia,
R. M. Cabezón,
I. Domí nguez
Abstract:
Li-depleted (enhanced) stars in the main-sequence (MS) and (or) the RGB, pose a puzzling mystery. Presently, there is still no clear answer to the mechanism(s) that enables such Li depletion (enhancement). One possible explanation comes from the, still controversial, observational evidence of Li underabundances in MS stars hosting planets, and of a positive correlation between the Li abundance and…
▽ More
Li-depleted (enhanced) stars in the main-sequence (MS) and (or) the RGB, pose a puzzling mystery. Presently, there is still no clear answer to the mechanism(s) that enables such Li depletion (enhancement). One possible explanation comes from the, still controversial, observational evidence of Li underabundances in MS stars hosting planets, and of a positive correlation between the Li abundance and rotational velocity in some RGB stars, which suggests a stellar collision with a planet-like object as a possible solution. In this study we explore this scenario, performing for first time 3D-hydrodynamical simulations of a 0.019 Mo brown dwarf collision with a MS star under different initial conditions. This enables us to gather information about the impact on the physical structure and final Li content in the hosting star.
△ Less
Submitted 14 February, 2020;
originally announced February 2020.
-
Two-level Dynamic Load Balancing for High Performance Scientific Applications
Authors:
Ali Mohammed,
Aurelien Cavelan,
Florina M. Ciorba,
Ruben M. Cabezon,
Ioana Banicesu
Abstract:
Scientific applications are often complex, irregular, and computationally-intensive. To accommodate the ever-increasing computational demands of scientific applications, high-performance computing (HPC) systems have become larger and more complex, offering parallelism at multiple levels (e.g., nodes, cores per node, threads per core). Scientific applications need to exploit all the available multi…
▽ More
Scientific applications are often complex, irregular, and computationally-intensive. To accommodate the ever-increasing computational demands of scientific applications, high-performance computing (HPC) systems have become larger and more complex, offering parallelism at multiple levels (e.g., nodes, cores per node, threads per core). Scientific applications need to exploit all the available multilevel hardware parallelism to harness the available computational power. The performance of applications executing on such HPC systems may adversely be affected by load imbalance at multiple levels, caused by problem, algorithmic, and systemic characteristics. Nevertheless, most existing load balancing methods do not simultaneously address load imbalance at multiple levels. This work investigates the impact of load imbalance on the performance of three scientific applications at the thread and process levels. We jointly apply and evaluate selected dynamic loop self-scheduling (DLS) techniques to both levels. Specifically, we employ the extended LaPeSD OpenMP runtime library at the thread level and extend the DLS4LB MPI-based dynamic load balancing library at the process level. This approach is generic and applicable to any multiprocess-multithreaded computationally-intensive application (programmed using MPI and OpenMP). We conduct an exhaustive set of experiments to assess and compare six DLS techniques at the thread level and eleven at the process level. The results show that improved application performance, by up to 21%, can only be achieved by jointly addressing load imbalance at the two levels. We offer insights into the performance of the selected DLS techniques and discuss the interplay of load balancing at the thread level and process level.
△ Less
Submitted 15 November, 2019;
originally announced November 2019.
-
Core-collapse supernovae in the hall of mirrors. A three-dimensional code-comparison project
Authors:
Rubén M. Cabezón,
Kuo-Chuan Pan,
Matthias Liebendörfer,
Takami Kuroda,
Kevin Ebinger,
Oliver Heinimann,
Friedrich-Karl Thielemann,
Albino Perego
Abstract:
Modeling core-collapse supernovae (CCSNe) with neutrino transport in three dimensions (3D) requires tremendous computing resources and some level of approximation. We present a first comparison study of CCSNe in 3D with different physics approximations and hydrodynamics codes. We aim to assess the impact of the hydrodynamics code, approximations for the neutrino and gravity treatments, and rotatio…
▽ More
Modeling core-collapse supernovae (CCSNe) with neutrino transport in three dimensions (3D) requires tremendous computing resources and some level of approximation. We present a first comparison study of CCSNe in 3D with different physics approximations and hydrodynamics codes. We aim to assess the impact of the hydrodynamics code, approximations for the neutrino and gravity treatments, and rotation on the simulation of CCSNe in 3D. We use four different hydrodynamics codes in this work (ELEPHANT, FLASH, fGR1, and SPHYNX) in combination with two different neutrino treatments, the isotropic diffusion source approximation (IDSA) and two-moment M1, and three different gravity treatments: Newtonian, 1D General Relativity (GR) correction, and full GR). Additional parameters discussed in this study are the inclusion of neutrino-electron scattering via a parametrized deleptonization (PD) and the influence of rotation. The four codes compared in this work include Eulerian and fully Lagrangian (smoothed particle hydrodynamics) codes for the first time. They show agreement in the overall evolution of the collapse phase and early post-bounce within the range of 10% (20% in some cases). The comparison of the different neutrino treatments highlights the need to further investigate the antineutrino luminosities in IDSA, which tend to be relatively high. We also demonstrate the requirement for a more detailed heavy-lepton neutrino leakage. When comparing with a full GR code, including an M1 transport method, we confirm the influence of neutrino-electron scattering during the collapse phase, which is adequately captured by the PD scheme. Also, the effective GR potential reproduces the overall dynamic evolution correctly in all Newtonian codes. Additionally, we verify that rotation aids the shock expansion and estimate the overall angular momentum losses for each code in rotating scenarios.
△ Less
Submitted 18 December, 2018; v1 submitted 24 June, 2018;
originally announced June 2018.
-
Surface and core detonations in rotating white dwarfs
Authors:
Domingo García-Senz,
Rubén M. Cabezón,
Inmaculada Domínguez
Abstract:
The feasibility of the Double Detonation mechanism, -a surface Helium-detonation followed by the complete carbon detonation of the core-, in a rotating white dwarf with a mass $\simeq 1 M_{\odot}$ is studied using three-dimensional hydrodynamic simulations. Assuming rigid rotation, the rotational speed is taken high enough as to considerably distort the initial spherical geometry of the white dwar…
▽ More
The feasibility of the Double Detonation mechanism, -a surface Helium-detonation followed by the complete carbon detonation of the core-, in a rotating white dwarf with a mass $\simeq 1 M_{\odot}$ is studied using three-dimensional hydrodynamic simulations. Assuming rigid rotation, the rotational speed is taken high enough as to considerably distort the initial spherical geometry of the white dwarf. Unlike spherically symmetric models, we found that when helium ignition is located far from the spinning axis the detonation fronts converge asynchronically at the antipodes of the igniting point. Nevertheless, the detonation of the carbon core still remains as the most probable outcome. The detonation of the core gives rise to a strong explosion, matching many of the basic observational constraints of Type Ia Supernova. We conclude that the Double Detonation mechanism also works when the white dwarf is spinning fast. This confirms the sub-Chandrasekhar-mass models and, maybe some Double Degenerate models (those having some helium fuel at the merging moment), as appealing channels to produce Type Ia Supernova events.
△ Less
Submitted 31 January, 2018;
originally announced January 2018.
-
SPHYNX: an accurate density-based SPH method for astrophysical applications
Authors:
Ruben M. Cabezon,
Domingo Garcia-Senz,
Joana Figueira
Abstract:
Hydrodynamical instabilities and shocks are ubiquitous in astrophysical scenarios. Therefore, an accurate numerical simulation of these phenomena is mandatory to correctly model and understand many astrophysical events, such as Supernovas, stellar collisions, or planetary formation. In this work, we attempt to address many of the problems that the smoothed particle hydrodynamics (SPH) technique ha…
▽ More
Hydrodynamical instabilities and shocks are ubiquitous in astrophysical scenarios. Therefore, an accurate numerical simulation of these phenomena is mandatory to correctly model and understand many astrophysical events, such as Supernovas, stellar collisions, or planetary formation. In this work, we attempt to address many of the problems that the smoothed particle hydrodynamics (SPH) technique has when dealing with subsonic hydrodynamical instabilities or shocks. To that aim we built a new SPH code named SPHYNX, that includes many of the recent advances in the SPH technique and some other new ones, which we present here. SPHYNX is of Newtonian type and grounded in the Euler-Lagrange formulation of the SPH technique. Its distinctive features are: the use of an integral approach to estimating the gradients; the use of a flexible family of interpolators called sinc kernels, which suppress pairing instability; and the incorporation of a new type of volume element which provides a better partition of the unity. Unlike other modern formulations, which consider volume elements linked to pressure, our volume element choice relies on density. SPHYNX conserves mass, linear and angular momentum, energy, entropy, and preserves kernel normalization even in strong shocks. The coupling between the integral approach to calculate gradients and the new family of volume elements reduces the so-called tensile instability. Both features help to suppress the damp which often prevents the growth of hydrodynamic instabilities in regular SPH codes. On the whole, SPHYNX has passed the verification tests described below. For identical particle setting and initial conditions the results were similar (or better in some particular cases) than those obtained with other SPH schemes such as GADGET-2, PSPH or with the recent density-independent formulation (DISPH) and conservative reproducing kernel (CRKSPH) techniques.
△ Less
Submitted 7 July, 2017; v1 submitted 6 July, 2016;
originally announced July 2016.
-
Type Ia Supernovae: Can Coriolis force break the symmetry of the gravitational confined detonation explosion mechanism?
Authors:
D. García-Senz,
R. M. Cabezón,
I. Domínguez,
F. K. Thielemann
Abstract:
Nowadays the number of models aimed at explaining the Type Ia supernova phenomenon is high and discriminating between them is a must-do. In this work we explore the influence of rotation in the evolution of the nuclear flame which drives the explosion in the so called gravitational confined detonation models. Assuming that the flame starts in a point-like region slightly above the center of the wh…
▽ More
Nowadays the number of models aimed at explaining the Type Ia supernova phenomenon is high and discriminating between them is a must-do. In this work we explore the influence of rotation in the evolution of the nuclear flame which drives the explosion in the so called gravitational confined detonation models. Assuming that the flame starts in a point-like region slightly above the center of the white dwarf (WD) and adding a moderate amount of angular velocity to the star we follow the evolution of the deflagration using a smoothed particle hydrodynamics code. We find that the results are very dependent on the angle between the rotational axis and the line connecting the initial bubble of burned material with the center of the white dwarf at the moment of the ignition. The impact of rotation is larger for angles close to 90° because the Coriolis force on a floating element of fluid is maximum, and its principal effect is to break the symmetry of the deflagration. Such symmetry breaking weakens the convergence of the nuclear flame at the antipodes of the initial ignition volume, changing the environmental conditions around the convergence region with respect to non-rotating models. These changes seem to disfavor the emergence of a detonation in the compressed volume at the antipodes, thus compromising the viability of the so called gravitational confined detonation mechanism.
△ Less
Submitted 21 July, 2015;
originally announced July 2015.
-
Equalizing resolution in smoothed-particle hydrodynamics calculations using self-adaptive sinc kernels
Authors:
Domingo García-Senz,
Rubén M. Cabezón,
José A. Escartín,
Kevin Ebinger
Abstract:
The smoothed-particle hydrodynamics (SPH) technique is a numerical method for solving gas-dynamical problems. It has been applied to simulate the evolution of a wide variety of astrophysical systems. The method has a second-order accuracy, with a resolution that is usually much higher in the compressed regions than in the diluted zones of the fluid. In this work, we propose and check a scheme to…
▽ More
The smoothed-particle hydrodynamics (SPH) technique is a numerical method for solving gas-dynamical problems. It has been applied to simulate the evolution of a wide variety of astrophysical systems. The method has a second-order accuracy, with a resolution that is usually much higher in the compressed regions than in the diluted zones of the fluid. In this work, we propose and check a scheme to balance and equalize the resolution of SPH between high- and low-density regions. This method relies on the versatility of a family of interpolators called Sinc kernels, which allows increasing the interpolation quality by varying only a single parameter (the exponent of the Sinc function). The scheme is checked and validated through a number of numerical tests, from standard one-dimensional Riemann problems in shock tubes, to multidimensional simulations of explosions, hydrodynamic instabilities and the collapse of a sun-like polytrope. The analysis of the hydrodynamical simulations suggests that the scheme devised to equalizing accuracy improves the treatment of the post-shock regions and, in general, of the rarefacted zones of fluids while causing no harm to the growth of hydrodynamic instabilities. The method is robust and easy to implement with a low computational overload. It conserves mass, energy, and momentum and reduces to the standard SPH scheme in regions of the fluid that have smooth density gradients.
△ Less
Submitted 25 August, 2014; v1 submitted 28 May, 2014;
originally announced May 2014.
-
High resolution simulations of the head-on collision of white dwarfs
Authors:
D. García-Senz,
R. M. Cabezón,
A. Arcones,
A. Relaño,
F. K. Thielemann
Abstract:
The direct impact of white dwarfs has been suggested as a plausible channel for type Ia supernovae. In spite of their (a priori) rareness, in highly populated globular clusters and in galactic centers, where the amount of white dwarfs is considerable, the rate of violent collisions between two of them might be non-negligible. Even more, there are indications that binary white dwarf systems orbited…
▽ More
The direct impact of white dwarfs has been suggested as a plausible channel for type Ia supernovae. In spite of their (a priori) rareness, in highly populated globular clusters and in galactic centers, where the amount of white dwarfs is considerable, the rate of violent collisions between two of them might be non-negligible. Even more, there are indications that binary white dwarf systems orbited by a third stellar-mass body have an important chance to induce a clean head-on collision. Therefore, this scenario represents a source of contamination for the supernova light-curves sample that it is used as standard candles in cosmology, and it deserves further investigation. Some groups have conducted numerical simulations of this scenario, but their results show several differences. In this paper we address some of the possible sources of these differences, presenting the results of high resolution hydrodynamical simulations jointly with a detailed nuclear post-processing of the nuclear abundances, to check the viability of white dwarf collisions to produce significant amounts of 56Ni. To that purpose, we use a 2D-axial symmetric smoothed particle hydrodynamic code to obtain a resolution considerably higher than in previous studies. In this work, we also study how the initial mass and nuclear composition affect the results. The gravitational wave emission is also calculated, as this is a unique signature of this kind of events. All calculated models produce a significant amount of 56Ni, ranging from 0.1 Msun to 1.1 Msun, compatible not only with normal-Branch type Ia supernova but also with the subluminous and super-Chandrasekhar subset. Nevertheless, the distribution mass-function of white dwarfs favors collisions among 0.6-0.7 Msun objects, leading to subluminous events.
△ Less
Submitted 26 September, 2013;
originally announced September 2013.
-
Testing the concept of integral approach to derivatives within the smoothed particle hydrodynamics technique in astrophysical scenarios
Authors:
Ruben M. Cabezon,
Domingo Garcia-Senz,
Jose Antonio Escartin
Abstract:
The behavior of IAD_0 scheme, a fully conservative SPH scheme based on a tensor formulation, is analyzed in connection with several astrophysical scenarios, and compared to the same simulations carried out with the standard SPH technique. The proposed hydrodynamic scheme is validated using a variety of numerical tests that cover important topics in astrophysics, such as the evolution of supernova…
▽ More
The behavior of IAD_0 scheme, a fully conservative SPH scheme based on a tensor formulation, is analyzed in connection with several astrophysical scenarios, and compared to the same simulations carried out with the standard SPH technique. The proposed hydrodynamic scheme is validated using a variety of numerical tests that cover important topics in astrophysics, such as the evolution of supernova remnants, the stability of self-gravitating bodies and the coalescence of compact objects. The results suggest that the SPH scheme built with the integral approach to the derivatives premise improves the results of the standard SPH technique. In particular, it is observed a better development of hydrodynamic instabilities, an improved description of self-gravitant structures in equilibrium and a reasonable description of the process of coalescence of two white dwarfs. A good energy, and linear and angular momentum conservation, generally better than that of standard SPH, was also obtained. In addition the new scheme is less susceptible to suffer pairing instability.
△ Less
Submitted 28 August, 2012; v1 submitted 23 July, 2012;
originally announced July 2012.
-
Improving smoothed particle hydrodynamics with an integral approach to calculating gradients
Authors:
Domingo Garcia-Senz,
Ruben M. Cabezon,
Jose Antonio Escartin
Abstract:
In this paper we develop and test a fully conservative SPH scheme based on a tensor formulation that can be applied to simulate astrophysical systems. In the proposed scheme, derivatives are calculated from an integral expression that leads to a tensor (instead of a vectorial) estimation of gradients and reduces to the standard formulation in the continuum limit. The new formulation improves the i…
▽ More
In this paper we develop and test a fully conservative SPH scheme based on a tensor formulation that can be applied to simulate astrophysical systems. In the proposed scheme, derivatives are calculated from an integral expression that leads to a tensor (instead of a vectorial) estimation of gradients and reduces to the standard formulation in the continuum limit. The new formulation improves the interpolation of physical magnitudes, leading to a set of conservative equations that resembles those of standard SPH. The resulting scheme is verified using a variety of well-known tests, all of them simulated in two dimensions. We also discuss an application of the proposed tensor method to astrophysics by simulating the stability of a Sun-like polytrope calculated in three dimensions.
△ Less
Submitted 5 December, 2011; v1 submitted 14 November, 2011;
originally announced November 2011.
-
Pulsating reverse detonation models of Type Ia supernovae. II: Explosion
Authors:
Eduardo Bravo,
Domingo Garcia-Senz,
Ruben M. Cabezon,
Inmaculada Dominguez
Abstract:
Observational evidences point to a common explosion mechanism of Type Ia supernovae based on a delayed detonation of a white dwarf. However, all attempts to find a convincing ignition mechanism based on a delayed detonation in a destabilized, expanding, white dwarf have been elusive so far. One of the possibilities that has been invoked is that an inefficient deflagration leads to pulsation of a…
▽ More
Observational evidences point to a common explosion mechanism of Type Ia supernovae based on a delayed detonation of a white dwarf. However, all attempts to find a convincing ignition mechanism based on a delayed detonation in a destabilized, expanding, white dwarf have been elusive so far. One of the possibilities that has been invoked is that an inefficient deflagration leads to pulsation of a Chandrasekhar-mass white dwarf, followed by formation of an accretion shock that confines a carbon-oxygen rich core, while transforming the kinetic energy of the collapsing halo into thermal energy of the core, until an inward moving detonation is formed. This chain of events has been termed Pulsating Reverse Detonation (PRD). In this work we present three dimensional numerical simulations of PRD models from the time of detonation initiation up to homologous expansion. Different models characterized by the amount of mass burned during the deflagration phase, M_defl, give explosions spanning a range of kinetic energies, K ~ (1.0-1.2) foes, and 56Ni masses, M(56Ni) ~ 0.6-0.8 M_sun, which are compatible with what is expected for typical Type Ia supernovae. Spectra and light curves of angle-averaged spherically symmetric versions of the PRD models are discussed. Type Ia supernova spectra pose the most stringent requirements on PRD models.
△ Less
Submitted 20 January, 2009;
originally announced January 2009.
-
Axisymmetric smoothed particle hydrodynamics with self-gravity
Authors:
D. Garcia-Senz,
A. Relano,
R. M. Cabezon,
E. Bravo
Abstract:
The axisymmetric form of the hydrodynamic equations within the smoothed particle hydrodynamics (SPH) formalism is presented and checked using idealized scenarios taken from astrophysics (free fall collapse, implosion and further pulsation of a sun-like star), gas dynamics (wall heating problem, collision of two streams of gas) and inertial confinement fusion (ICF, -ablative implosion of a small…
▽ More
The axisymmetric form of the hydrodynamic equations within the smoothed particle hydrodynamics (SPH) formalism is presented and checked using idealized scenarios taken from astrophysics (free fall collapse, implosion and further pulsation of a sun-like star), gas dynamics (wall heating problem, collision of two streams of gas) and inertial confinement fusion (ICF, -ablative implosion of a small capsule-). New material concerning the standard SPH formalism is given. That includes the numerical handling of those mass points which move close to the singularity axis, more accurate expressions for the artificial viscosity and the heat conduction term and an easy way to incorporate self-gravity in the simulations. The algorithm developed to compute gravity does not rely in any sort of grid, leading to a numerical scheme totally compatible with the lagrangian nature of the SPH equations.
△ Less
Submitted 11 October, 2008; v1 submitted 10 October, 2008;
originally announced October 2008.
-
A one-parameter family of interpolating kernels for Smoothed Particle Hydrodynamics studies
Authors:
Ruben M. Cabezon,
Domingo Garcia-Senz,
Antonio Relaño
Abstract:
A set of interpolating functions of the type f(v)={(sin[v pi/2])/(v pi/2)}^n is analyzed in the context of the smoothed-particle hydrodynamics (SPH) technique. The behaviour of these kernels for several values of the parameter n has been studied either analytically as well as numerically in connection with several tests carried out in two dimensions. The main advantage of this kernel relies in i…
▽ More
A set of interpolating functions of the type f(v)={(sin[v pi/2])/(v pi/2)}^n is analyzed in the context of the smoothed-particle hydrodynamics (SPH) technique. The behaviour of these kernels for several values of the parameter n has been studied either analytically as well as numerically in connection with several tests carried out in two dimensions. The main advantage of this kernel relies in its flexibility because for n=3 it is similar to the standard widely used cubic-spline, whereas for n>3 the interpolating function becomes more centrally condensed, being well suited to track discontinuities such as shock fronts and thermal waves.
△ Less
Submitted 16 September, 2008;
originally announced September 2008.
-
Constraining deflagration models of Type Ia supernovae through intermediate-mass elements
Authors:
D. Garcia-Senz,
E. Bravo,
R. M. Cabezon,
S. E. Woosley
Abstract:
The physical structure of a nuclear flame is a basic ingredient of the theory of Type Ia supernovae (SNIa). Assuming an exponential density reduction with several characteristic times we have followed the evolution of a planar nuclear flame in an expanding background from an initial density 6.6 10^7 g/cm3 down to 2 10^6 g/cm3. The total amount of synthesized intermediate-mass elements (IME), fro…
▽ More
The physical structure of a nuclear flame is a basic ingredient of the theory of Type Ia supernovae (SNIa). Assuming an exponential density reduction with several characteristic times we have followed the evolution of a planar nuclear flame in an expanding background from an initial density 6.6 10^7 g/cm3 down to 2 10^6 g/cm3. The total amount of synthesized intermediate-mass elements (IME), from silicon to calcium, was monitored during the calculation. We have made use of the computed mass fractions, X_IME, of these elements to give an estimation of the total amount of IME synthesized during the deflagration of a massive white dwarf. Using X_IME and adopting the usual hypothesis that turbulence decouples the effective burning velocity from the laminar flame speed, so that the relevant flame speed is actually the turbulent speed on the integral length-scale, we have built a simple geometrical approach to model the region where IME are thought to be produced. It turns out that a healthy production of IME involves the combination of not too short expansion times, t_c > 0.2 s, and high turbulent intensities. According to our results it could be difficult to produce much more than 0.2 solar masses of intermediate-mass elements within the deflagrative paradigma. The calculations also suggest that the mass of IME scales with the mass of Fe-peak elements, making it difficult to conciliate energetic explosions with low ejected nickel masses, as in the well observed SN1991bg or in SN1998de. Thus a large production of Si-peak elements, especially in combination with a low or a moderate production of iron, could be better addressed by either the delayed detonation route in standard Chandrasekhar-mass models or, perhaps, by the off-center helium detonation in the sub Chandrasekhar-mass scenario.
△ Less
Submitted 15 September, 2006;
originally announced September 2006.