We present three-dimensional local hydrodynamic simulations of flows around objects embedded within stellar envelopes using a "wind tunnel" formalism. Our simulations model the common envelope dynamical inspiral phase in binary star systems in terms of dimensionless flow characteristics. We present suites of simulations that study the effects of varying the binary mass ratio, stellar structure, equation of state, relative Mach number of the object's motion through the gas, and density gradients across the gravitational focusing scale. For each model, we measure coefficients of accretion and drag experienced by the embedded object. These coefficients regulate the coupled evolution of the object's masses and orbital tightening during the dynamical inspiral phase of the common envelope. We extrapolate our simulation results to accreting black holes with masses comparable to that of the population of LIGO black holes. We demonstrate that the mass and spin accrued by these black holes per unit orbital tightening are directly related to the ratio of accretion to drag coefficients. We thus infer that the mass and dimensionless spin of initially non-rotating black holes change by of order $1\%$ and 0.05, respectively, in a typical example scenario. Our prediction that the masses and spins of black holes remain largely unmodified by a common envelope phase aids in the interpretation of the properties of the growing observed population of merging binary black holes. Even if these black holes passed through a common envelope phase during their assembly, features of mass and spin imparted by previous evolutionary epochs should be preserved.

Binary stars are common. While only those with small separations may exchange gas with one another, even the widest binaries interact with their gaseous surroundings. Drag forces and accretion rates dictate how these systems are transformed by these interactions. We perform three-dimensional hydrodynamic simulations of Bondi-Hoyle-Lyttleton flows, in which a binary moves supersonically relative to a homogeneous medium, using the adaptive mesh refinement code FLASH. We simulate a range of values of the initial semi-major axis of the orbit relative to the gravitational focusing impact parameter of the pair. When the binary separation is less than the gravitational focusing impact parameter, the pair orbits within a shared bow shock. When the pair is wider, each object has an individual bow-shock structure. The long-term evolution of the binary is determined by the timescales for accretion, slowing of the center of mass, and orbital inspiral. We find a clear hierarchy of these timescales; a binary's center-of-mass motion is slowed over a shorter timescale than the pair inspirals or accretes. In contrast to previous analytic predictions, which assume an unperturbed background medium, we find that the timescale for orbital inspiral is proportional to the semi-major axis to the $0.19 \pm 0.01$ power. This positive scaling indicates that gaseous drag forces can drive binaries either to coalescence or to the critical separation at which gravitational radiation dominates their further evolution. We discuss the implications of our results for binaries embedded in the interstellar medium, active galactic nuclei disks, and common envelope phases.

An isolated star moving supersonically through a uniform gas accretes material from its gravitationally-induced wake. The rate of accretion is set by the accretion radius of the star and is well-described by classical Bondi-Hoyle-Lyttleton theory. Stars, however, are not born in isolation. They form in clusters where they accrete material that is influenced by all the stars in the cluster. We perform three-dimensional hydrodynamic simulations of clusters of individual accretors embedded in a uniform-density wind in order to study how the accretion rates experienced by individual cluster members are altered by the properties of the ambient gas and the cluster itself. We study accretion as a function of number of cluster members, mean separation between them, and size of their individual accretion radii. We determine the effect of these key parameters on the aggregate and individual accretion rates, which we compare to analytic predictions. We show that when the accretion radii of the individual objects in the cluster substantially overlap, the surrounding gas is effectively accreted into the collective potential of the cluster prior to being accreted onto the individual stars. We find that individual cluster members can accrete drastically more than they would in isolation, in particular when the flow is able to cool efficiently. This effect could potentially modify the luminosity of accreting compact objects in star clusters and could lead to the rejuvenation of young star clusters as well as globular clusters with low-inclination and low-eccentricity.

During a common envelope episode in a binary system, the engulfed companion spirals to tighter orbital separations under the influence of drag from the surrounding envelope material. As this object sweeps through material with a steep radial gradient of density, net angular momentum is introduced into the flow, potentially leading to the formation of an accretion disk. The presence of a disk would have dramatic consequences for the outcome of the interaction because accretion might be accompanied by strong, polar outflows with enough energy to unbind the entire envelope. Without a detailed understanding of the necessary conditions for disk formation during common envelope, therefore, it is difficult to accurately predict the population of merging compact binaries. This paper examines the conditions for disk formation around objects embedded within common envelopes using the `wind tunnel' formalism developed by MacLeod et al. (2017). We find that the formation of disks is highly dependent on the compressibility of the envelope material. Disks form only in the most compressible of stellar envelope gas, found in envelopes' outer layers in zones of partial ionization. These zones are largest in low-mass stellar envelopes, but comprise small portions of the envelope mass and radius in all cases. We conclude that disk formation and associated accretion feedback in common envelope is rare, and if it occurs, transitory. The implication for LIGO black hole binary assembly is that by avoiding strong accretion feedback, common envelope interactions should still result in the substantial orbital tightening needed to produce merging binaries.

This paper examines the properties of flows around objects embedded within common envelopes in the simplified context of a "wind tunnel." We establish characteristic relationships between key common envelope flow parameters like the Mach number and density scale height. Our wind tunnel is a three-dimensional, cartesian geometry hydrodynamic simulation setup that includes the gravity of the primary and secondary stars and allows us to study the coefficients of drag and accretion experienced by the embedded object. Accretion and drag lead to a transformation of an embedded object and its orbit during a common envelope phase. We present two suites of simulations spanning a range of density gradients and Mach numbers -- relevant for flow near the limb of a stellar envelope to the deep interior. In one suite, we adopt an ideal gas adiabatic exponent of $\gamma=5/3$, in the other, $\gamma=4/3$. We find that coefficients of drag rise in flows with steeper density gradients and that coefficients of drag and accretion are consistently higher in the more compressible, $\gamma=4/3$ flow. We illustrate the impact of these newly derived coefficients by integrating the inspiral of a secondary object through the envelopes of $3M_\odot$ ($\gamma\approx5/3$) and $80M_\odot$ ($\gamma\approx4/3$) giants. In these examples, we find a relatively rapid initial inspiral because, near the stellar limb, dynamical friction drag is generated mainly from dense gas focussed from deeper within the primary-star's envelope. This rapid initial inspiral timescale carries potential implications for the timescale of transients from early common envelope interaction.