We update the capabilities of the software instrument Modules for Experiments in Stellar Astrophysics (MESA) and enhance its ease of use and availability. Our new approach to locating convective boundaries is consistent with the physics of convection, and yields reliable values of the convective core mass during both hydrogen and helium burning phases. Stars with M < 8 M become white dwarfs and cool to the point where the electrons are degenerate and the ions are strongly coupled, a realm now available to study with MESA due to improved treatments of element diffusion, latent heat release, and blending of equations of state. Studies of the final fates of massive stars are extended in MESA by our addition of an approximate Riemann solver that captures shocks and conserves energy to high accuracy during dynamic epochs. We also introduce a 1D capability for modeling the effects of Rayleigh-Taylor instabilities that, in combination with the coupling to a public version of the STELLA radiation transfer instrument, creates new avenues for exploring Type II supernovae properties. These capabilities are exhibited with exploratory models of pair-instability supernova, pulsational pair-instability supernova, and the formation of stellar mass black holes. The applicability of MESA is now widened by the capability of importing multi-dimensional hydrodynamic models into MESA. We close by introducing software modules for handling floating point exceptions and stellar model optimization, and four new software tools − MESA-Web, MESA-Docker, pyMESA, and mesastar.org − to enhance MESA's education and research impact.
We present the results of 2D hydrodynamical simulations of circumbinary disk accretion using the finite-volume code DISCO. This code solves the 2D viscous Navier-Stokes equations on a highresolution moving mesh which shears with the fluid flow, greatly reducing advection errors in comparison with a fixed grid. We perform a series of simulations for binary mass ratios in the range 0.026 ≤ q ≤ 1.0, each lasting longer than a viscous time so that we reach a quasi-steady accretion state. In each case, we find that gas is efficiently stripped from the inner edge of the circumbinary disk and enters the cavity along accretion streams, which feed persistent "mini-disks" surrounding each black hole. We find that for q 0.1, the binary excites eccentricity in the inner region of the circumbinary disk, creating an overdense lump which gives rise to enhanced periodicity in the accretion rate. The dependence of the periodicity on mass ratio may provide a method for observationally inferring mass ratios from measurements of the accretion rate. We also find that for all mass ratios studied, the magnitude of the accretion onto the secondary is sufficient to drive the binary toward larger mass ratio. This suggests a mechanism for biasing mass ratio distributions toward equal mass.
By numerically integrating the compressible Navier-Stokes equations in two dimensions, we calculate the criterion for gap formation by a very low mass (q ∼ 10 −4 ) protoplanet on a fixed orbit in a thin viscous disk. In contrast with some previously proposed gap-opening criteria, we find that a planet can open a gap even if the Hill radius is smaller than the disk scale height. Moreover, in the lowviscosity limit, we find no minimum mass necessary to open a gap for a planet held on a fixed orbit. In particular, a Neptune-mass planet will open a gap in a minimum mass solar nebula with suitably low viscosity (α 10 −4 ). We find that the mass threshold scales as the square root of viscosity in the low mass regime. This is because the gap width for critical planet masses in this regime is a fixed multiple of the scale height, not of the Hill radius of the planet.
We present the results of 2D, moving mesh, viscous hydrodynamical simulations of accretion onto merging supermassive black hole (SMBH) binaries. We include viscous heating, shock heating, and radiative cooling, and simulate the transition from the "pre-decoupling" epoch, where the inspiral timescale is longer than the viscous timescale, to the "post-decoupling" epoch, where the inspiral timescale is shorter than the viscous timescale. We find that there is no abrupt halt to the accretion at decoupling, but rather the accretion shows a slow decay, with significant accretion well after the expected decoupling. Moreover, we find that the luminosity in X-rays is significantly higher prior to the merger, as orbital energy from the SMBH binary is converted to heat via strong shocks inside the cavity, and radiated away. Following the merger, the cavity refills viscously and the accretion rate relaxes to the Shakura-Sunyaev value, while the X-ray luminosity drops as the shocks quickly dissipate.
Using numerical hydrodynamics calculations and a novel method for densely sampling parameter space, we measure the accretion and torque on a binary system from a circumbinary disk. In agreement with some earlier studies, we find that the net torque on the binary is positive for mass ratios close to unity, and that accretion always drives the binary toward equal mass. Accretion variability depends sensitively on the numerical sink prescription, but the torque and relative accretion onto each component do not depend on the sink timescale. Positive torque and highly variable accretion occurs only for mass ratios greater than around 0.05. This means that for mass ratios below 0.05, the binary would migrate inward until the secondary accreted sufficient mass, after which it would execute a U-turn and migrate outward. We explore a range of viscosities, from α = 0.03 to α = 0.15, and find that this outward torque is proportional to the viscous torque, so that torque per unit accreted mass is independent of α. Dependence of accretion and torque on mass ratio is explored in detail, densely sampling mass ratios between 0.01 and unity. For mass ratio q > 0.2, accretion variability is found to exhibit a distinct sawtooth pattern, typically with a five-orbit cycle that provides a smoking gun prediction for variable quasars observed over long periods, as a potential means to confirm the presence of a binary.
We have generalized a method for the numerical solution of hyperbolic systems of equations using a dynamic Voronoi tessellation of the computational domain. The Voronoi tessellation is used to generate moving computational meshes for the solution of multi-dimensional systems of conservation laws in finite-volume form. The mesh generating points are free to move with arbitrary velocity, with the choice of zero velocity resulting in an Eulerian formulation. Moving the points at the local fluid velocity makes the formulation effectively Lagrangian. We have written the TESS code to solve the equations of compressible hydrodynamics and magnetohydrodynamics for both relativistic and nonrelativistic fluids on a dynamic Voronoi mesh. When run in Lagrangian mode, TESS is significantly less diffusive than fixed mesh codes and thus preserves contact discontinuities to high precision while also accurately capturing strong shock waves. TESS is written for Cartesian, spherical and cylindrical coordinates and is modular so that auxilliary physics solvers are readily integrated into the TESS framework and so that the TESS framework can be readily adapted to solve general systems of equations. We present results from a series of test problems to demonstrate the performance of TESS and to highlight some of the advantages of the dynamic tessellation method for solving challenging problems in astrophysical fluid dynamics.
Most standard descriptions of Type II migration state that massive, gap-opening planets must migrate at the viscous drift rate. This is based on the idea that the disk is separated into an inner and outer region and gas is considered unable to cross the gap. In fact, gas easily crosses the gap on horseshoe orbits, nullifying this necessary premise which would set the migration rate. In this work, it is demonstrated using highly accurate numerical calculations that the actual migration rate is dependent on disk and planet parameters, and can be significantly larger or smaller than the viscous drift rate. In the limiting case of a disk much more massive than the secondary, the migration rate saturates to a constant which is sensitive to disk parameters and is not necessarily of order viscous rate. In the opposite limit of a low-mass disk, the migration rate decreases linearly with disk mass. Steady-state solutions in the low disk mass limit show no pile-up outside the secondary's orbit, and no corresponding drainage of the inner disk.
The merger of binary neutron stars (BNSs) can lead to large amplifications of the magnetic field due to the development of turbulence and instabilities in the fluid, such as the Kelvin-Helmholtz shear instability, which drive small-scale dynamo activity. In order to properly resolve such instabilities and obtain the correct magnetic field amplification, one would need to employ resolutions that are currently unfeasible in global general relativistic magnetohydrodynamic (GRMHD) simulations of BNS mergers. Here, we present a subgrid model that allows global simulations to take into account the small-scale amplification of the magnetic field which is caused by the development of turbulence during BNS mergers. Assuming dynamo saturation, we show that magnetar-level fields (∼ 10 16 G) can be easily reached, and should therefore be expected from the merger of magnetized BNSs. The total magnetic energy can reach values up to ∼ 10 51 erg and the post-merger remnant can therefore emit strong electromagnetic signals and possibly produce short gamma-ray bursts.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.