Normalized to: Escartín, J.
[1]
oai:arXiv.org:1405.7207 [pdf] - 1450520
Equalizing resolution in smoothed-particle hydrodynamics calculations
using self-adaptive sinc kernels
Submitted: 2014-05-28, last modified: 2014-08-25
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.
[2]
oai:arXiv.org:1110.1949 [pdf] - 442299
Type Ia supernovae and the ^{12}C+^{12}C reaction rate
Submitted: 2011-10-10
The experimental determination of the cross-section of the ^{12}C+^{12}C
reaction has never been made at astrophysically relevant energies (E<2 MeV).
The profusion of resonances throughout the measured energy range has led to
speculation that there is an unknown resonance at E\sim1.5 MeV possibly as
strong as the one measured for the resonance at 2.14 MeV. We study the
implications that such a resonance would have for the physics of SNIa, paying
special attention to the phases that go from the crossing of the ignition curve
to the dynamical event. We use one-dimensional hydrostatic and hydrodynamic
codes to follow the evolution of accreting white dwarfs until they grow close
to the Chandrasekhar mass and explode as SNIa. In our simulations, we account
for a low-energy resonance by exploring the parameter space allowed by
experimental data. A change in the ^{12}C+^{12}C rate similar to the one
explored here would have profound consequences for the physical conditions in
the SNIa explosion, namely the central density, neutronization, thermal
profile, mass of the convective core, location of the runaway hot spot, or time
elapsed since crossing the ignition curve. For instance, with the largest
resonance strength we use, the time elapsed since crossing the ignition curve
to the supernova event is shorter by a factor ten than for models using the
standard rate of ^{12}C+^{12}C, and the runaway temperature is reduced from
\sim8.14\times10^{8} K to \sim4.26\times10^{8} K. On the other hand, a
resonance at 1.5 MeV, with a strength ten thousand times smaller than the one
measured at 2.14 MeV, but with an {\alpha}/p yield ratio substantially
different from 1 would have a sizeable impact on the degree of neutronization
of matter during carbon simmering. We conclude that a robust understanding of
the links between SNIa properties and their progenitors will not be attained
until the ^{12}C+^{12}C reaction rate is measured at energies \sim1.5 MeV.