Low Mach number modeling of Type I X-ray burst deflagrations

David J. Lin,Alvin Bayliss,Ronald E. Taam

DOI: https://doi.org/10.1086/508863

2006-09-05

Abstract:The Low Mach Number Approximation (LMNA) is applied to 2D hydrodynamical modeling of Type I X-ray bursts on a rectangular patch on the surface of a non-rotating neutron star. Because such phenomena involve decidedly subsonic flows, the timestep increase offered by the LMNA makes routine simulations of these deflagrations feasible in an environment where strong gravity produces significant stratification, while allowing for potentially significant lateral differences in temperature and density. The model is employed to simulate the heating, peak, and initial cooling stages in the deep envelope layers of a burst. During the deflagration, Benard-like cells naturally fill up a vertically expanding convective layer. The Mach number is always less than 0.15 throughout the simulation, thus justifying the low Mach number approximation. While the convective layer is superadiabatic on average, significant fluctuations in adiabaticity occur within it on subconvective timescales. Due to convective layer expansion, significant compositional mixing naturally occurs, but tracer particle penetration through the convective layer boundaries on convective timescales is temporary and spatially limited. Thus, mixing occurs on the relatively slow burst timescale through thermal expansion of the convective layer rather than from mass penetration of the convective layer boundary through particle convection. At the convective layer boundaries where mixing is less efficient, the actual temperature gradient more closely follows the Ledoux criteria.

Astrophysics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the modeling of the low - Mach - number combustion process in Type I X - ray bursts occurring on the surface of neutron stars. Specifically, the author applies the Low Mach Number Approximation (LMNA) to conduct two - dimensional hydrodynamic simulations to study the combustion and convection processes in these outburst phenomena.

### Main problems
1. **Simulation of the low - Mach - number combustion process**: Since Type I X - ray bursts involve subsonic flows, traditional fully compressible simulation methods will lead to an extremely small time step and a prohibitively high computational cost. Therefore, using LMNA can significantly increase the time step, making conventional simulations possible.
2. **Stratified structure in a strong - gravity environment**: On the surface of a neutron star, strong gravity will produce a significant stratified structure, and at the same time, the lateral temperature and density differences are also relatively obvious. How to conduct accurate simulations in such an environment is a challenge.
3. **Super - adiabaticity and mixing processes in the convective layer**: During the combustion process, whether the convective layer always remains in a super - adiabatic state and what the specific mechanism of component mixing is are all issues that need to be explored.

### Main contributions of the paper
- **Introduction of the low - Mach - number approximation**: By introducing LMNA, the influence of sound waves is excluded, the problems of sound - wave time - step limitation and boundary reflection are avoided, and the computational time is greatly reduced.
- **Consideration of lateral fluctuations**: Allowing for significant lateral temperature and density fluctuations more realistically simulates the actual physical process.
- **Complete physical inputs**: Including thermal diffusion, a real equation of state, and a combustion network, making the model more comprehensive and accurate.

### Simulation results
- **Height and characteristics of the combustion layer**: The height of the combustion layer is roughly maintained at about 500 cm from the bottom of the distance domain, and the maximum Energy Generation Rate (EGR) reaches \(2\times10^{19}\text{erg g}^{-1}\text{s}^{-1}\).
- **Characteristics of the convective layer**: During the combustion process, the convective layer naturally forms a structure similar to Bénard cells, and the adiabaticity within the convective layer fluctuates significantly on different time scales.
- **Component - mixing mechanism**: Component mixing mainly occurs on a slower outburst time scale rather than through particle convection penetrating the convective - layer boundary.

In conclusion, by introducing the low - Mach - number approximation, this paper successfully solves several key problems in the simulation of Type I X - ray bursts, providing a new perspective and tool for understanding this complex astrophysical phenomenon.