sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Carroll-Nellenback, Jonathan

Normalized to: Carroll-Nellenback, J.

26 article(s) in total. 35 co-authors, from 1 to 25 common article(s). Median position in authors list is 3,5.

[1]  oai:arXiv.org:1912.01647  [pdf] - 2107803
Bipolar Planetary Nebulae from Outflow Collimation by Common Envelope Evolution
Comments: 17 pages, 16 figures
Submitted: 2019-12-03, last modified: 2020-06-03
The morphology of bipolar planetary nebulae (PNe) can be attributed to interactions between a fast wind from the central engine and dense toroidal shaped ejecta left over from common envelope (CE) evolution. Here we use the 3-D hydrodynamic AMR code AstroBEAR to study the possibility that bipolar PN outflows can emerge collimated even from an uncollimated spherical wind in the aftermath of a CE event. The output of a single CE simulation via the SPH code PHANTOM serves as the initial conditions. Four cases of winds, all with high enough momenta to account for observed high momenta preplanetary nebula outflows, are injected spherically from the region of the CE binary remnant into the ejecta. We compare cases with two different momenta and cases with no radiative cooling versus application of optically thin emission via a cooling curve to the outflow. Our simulations show that in all cases highly collimated bipolar outflows result from deflection of the spherical wind via the interaction with the CE ejecta. Significant asymmetries between the top and bottom lobes are seen in all cases. The asymmetry is strongest for the lower momentum case with radiative cooling. While real post CE winds may be aspherical, our models show that collimation via "inertial confinement" will be strong enough to create jet-like outflows even beginning with maximally uncollimated drivers. Our simulations reveal detailed shock structures in the shock focused inertial confinement (SFIC) model and develop a lens-shaped inner shock that is a new feature of SFIC driven bipolar lobes.
[2]  oai:arXiv.org:2004.06829  [pdf] - 2103438
Common Envelope Evolution on the Asymptotic Giant Branch: Unbinding within a Decade?
Comments: 12 pages, 12 figures, published in MNRAS
Submitted: 2020-04-14, last modified: 2020-05-27
Common envelope (CE) evolution is a critical but still poorly understood progenitor phase of many high-energy astrophysical phenomena. Although 3D global hydrodynamic CE simulations have become more common in recent years, those involving an asymptotic giant branch (AGB) primary are scarce, due to the high computational cost from the larger dynamical range compared to red giant branch (RGB) primaries. But CE evolution with AGB progenitors is desirable to simulate because such events are the likely progenitors of most bi-polar planetary nebulae (PNe), and prominent observational testing grounds for CE physics. Here we present a high resolution global simulation of CE evolution involving an AGB primary and $1\,\mathrm{M}_\odot$ secondary, evolved for $20$ orbital revolutions. During the last $16$ of these orbits, the envelope unbinds at an almost constant rate of about $0.1$-$0.2\,\mathrm{M}_\odot\,\mathrm{yr}^{-1}$. If this rate were maintained, the envelope would be unbound in less than $10\,\mathrm{yr}$. The dominant source of this unbinding is consistent with inspiral; we assess the influence of the ambient medium to be subdominant. We compare this run with a previous run that used an RGB phase primary evolved from the same $2\,\mathrm{M}_\odot$ main sequence star to assess the influence of the evolutionary state of the primary. When scaled appropriately, the two runs are quite similar, but with some important differences.
[3]  oai:arXiv.org:1910.08027  [pdf] - 2076687
A 3D radiation-hydrodynamic AGB binary model
Comments: 21 pages, 8 figures, and 6 tables
Submitted: 2019-10-17, last modified: 2020-02-28
The origin of chemically peculiar stars and non-zero eccentricity in evolved close binaries have been long-standing problems in binary stellar evolution. Answers to these questions may trace back to an intense mass transfer during AGB binary phase. In this work, we use AstroBEAR to solve the 3D radiation-hydrodynamic equations and calculate the mass transfer rate in asymptotic-giant-branch (AGB) binaries that undergo the wind-Roche-lobe-overflow or Bondi-Hoyle-Lyttleton (BHL) accretion. MESA produces the density and temperature of the boundary condition of the AGB star. To improve the resolution of the dynamics of a circumbinary disk, we implement an azimuthal angle-dependent 3D radiation transfer. We consider optically thin cooling and obtain the number density of the coolants by solving Saha equations. One of the goals of this work is to illustrate the transition from the wind-Roche-lobe-overflow to BHL accretion. Both circumbinary disks and spiral structure outflows can appear in the simulations. Circumbinary disks may form when the optical thickness in the equatorial region increases. The increase of the optical thickness is due to the deflected wind. The resulting mass transfer efficiency in our models is up to a factor of 8 times higher than what the standard BHL accretion scenario predicts, and the outflow gains up to 91% of its initial angular momentum when it reaches 1.3 binary separations. Consequently, some AGB binaries may undergo orbit shrinkage, and some will expand. The high mass transfer efficiency is closely related to the presence of the circumbinary disks.
[4]  oai:arXiv.org:1902.04450  [pdf] - 2043498
The Fermi Paradox and the Aurora Effect: Exo-civilization Settlement, Expansion and Steady States
Comments:
Submitted: 2019-02-12, last modified: 2020-02-05
We model the settlement of the galaxy by space-faring civilizations in order to address issues related to the Fermi Paradox. We explore the problem in a way that avoids assumptions about the intent and motivation of any exo-civilization seeking to settle other planetary systems. We first consider the speed of an advancing settlement via probes of finite velocity and range to determine if the galaxy can become inhabited with space-faring civilizations on timescales shorter than its age. We also include the effect of stellar motions on the long term behavior of the settlement front which adds a diffusive component to its advance. The results of these models demonstrate that the Milky Way can be readily 'filled-in' with settled stellar systems under conservative assumptions about interstellar spacecraft velocities and launch rates. We then consider the question of the galactic steady-state achieved in terms of the fraction of settled planets. We do this by considering the effect of finite settlement civilization lifetimes on the steady states. We find a range of parameters for which the galaxy supports a population of interstellar space-faring civilizations even though some settleable systems are uninhabited. Both results point to ways in which Earth might remain unvisited in the midst of an inhabited galaxy. Finally we consider how our results can be combined with the finite horizon for evidence of previous settlements in Earth's geologic record. Our steady-state model can constrain the probabilities for an Earth visit by a settling civilization before a given time horizon. These results break the link between Hart's famous "Fact A" (no interstellar visitors on Earth now) and the conclusion that humans must, therefore, be the only technological civilization in the galaxy.
[5]  oai:arXiv.org:1906.00075  [pdf] - 2046213
Effects of Radiation Pressure on the Evaporative Wind of HD 209458b
Comments:
Submitted: 2019-05-31, last modified: 2020-02-03
The role of radiation pressure in shaping exoplanet photoevaporation remains a topic of contention. Radiation pressure from the exoplanet's host star has been proposed as a mechanism to drive the escaping atmosphere into a "cometary" tail and explain the high velocities observed in systems where mass loss is occurring. In this paper we present results from high-resolution 3-D hydrodynamic simulations of a planet similar to HD 209458b. We self-consistently launch a wind flowing outward from the planet by calculating the ionization and heating resulting from incident high-energy radiation, and account for radiation pressure. We first present a simplified calculation, setting a limit on the Lyman-$\alpha$ flux required to drive the photo-evaporated planetary material to larger radii and line-of-sight velocities. We then present the results of our simulations, which confirm the limits determined by our analytic calculation. We thus demonstrate that, within the limits of our hydrodynamic simulation and for the Lyman-$\alpha$ fluxes expected for HD 209458, radiation pressure is unlikely to significantly affect photoevaporative winds or to explain the high velocities at which wind material is observed, though further possibilities remain to be investigated.
[6]  oai:arXiv.org:1908.06195  [pdf] - 1979681
How Drag Force Evolves in Global Common Envelope Simulations
Comments: 14 pages, 8 figures, MNRAS in press
Submitted: 2019-08-16, last modified: 2019-10-03
We compute the forces, torque and rate of work on the companion-core binary due to drag in global simulations of common envelope (CE) evolution for three different companion masses. Our simulations help to delineate regimes when conventional analytic drag force approximations are applicable. During and just prior to the first periastron passage of the in-spiral phase, the drag force is reasonably approximated by conventional analytic theory and peaks at values proportional to the companion mass. Good agreement between global and local 3D "wind tunnel" simulations, including similar net drag force and flow pattern, is obtained for comparable regions of parameter space. However, subsequent to the first periastron passage, the drag force is up to an order of magnitude smaller than theoretical predictions, quasi-steady, and depends only weakly on companion mass. The discrepancy is exacerbated for larger companion mass and when the inter-particle separation reduces to the Bondi-Hoyle-Lyttleton accretion radius, creating a turbulent thermalized region. Greater flow symmetry during this phase leads to near balance of opposing gravitational forces in front of and behind the companion, hence a small net drag. The reduced drag force at late times helps explain why companion-core separations necessary for envelope ejection are not reached by the end of limited duration CE simulations.
[7]  oai:arXiv.org:1908.00157  [pdf] - 1971328
Hydrodynamic Simulations of Disrupted Planetary Accretion Discs Inside the Core of an AGB Star
Comments: Submitted to MNRAS
Submitted: 2019-07-31
Volume complete sky surveys provide evidence for a binary origin for the formation of isolated white dwarfs with magnetic fields in excess of a MegaGauss. Interestingly, not a single high-field magnetic white dwarf has been found in a detached system suggesting that if the progenitors are indeed binaries, the companion must be removed or merge during formation. An origin scenario consistent with observations involves the engulfment, inspiral, and subsequent tidal disruption of a low-mass companion in the interior of a giant star during a common envelope phase. Material from the shredded companion forms a cold accretion disc embedded in the hot ambient around the proto-white dwarf. Entrainment of hot material may evaporate the disc before it can sufficiently amplify the magnetic field, which typically requires at least a few orbits of the disc.Using three-dimensional hydrodynamic simulations of accretion discs with masses between 1 and 10 times the mass of Jupiter inside the core of an Asymptotic Giant Branch star, we find that the discs survive for at least 10 orbits (and likely for 100 orbits), sufficient for strong magnetic fields to develop.
[8]  oai:arXiv.org:1902.05063  [pdf] - 1953410
Hydrodynamic and Magnetohydrodynamic Simulations of Wire Turbulence
Comments: Accepted to High Energy Density Physics, in press, 10 pages, 7 figures
Submitted: 2019-02-12, last modified: 2019-07-10
We report on simulations of laboratory experiments in which magnetized supersonic flows are driven through a wire mesh. The goal of the study was to investigate the ability of such a configuration to generate supersonic, MHD turbulence. We first report on the morphological structures that develop in both magnetized and non-magnetized cases. We then analyze the flow using a variety of statistical measures, including power spectra and probability distribution functions of the density. Using these results we estimate the sonic mach number in the flows downstream of the wire mesh. We find the initially hypersonic (M=20) planar shock through the wire mesh does lead to downstream turbulent conditions. However, in both magnetized and non-magnetized cases, the resultant turbulence was marginally supersonic to transonic (M~1), and highly anisotropic in structure.
[9]  oai:arXiv.org:1812.11196  [pdf] - 1863939
Energy Budget and Core-Envelope Motion in Common Envelope Evolution
Comments: 17 pages, 10 figures, 5 tables, accepted for publication in MNRAS
Submitted: 2018-12-28, last modified: 2019-03-25
We analyze a 3D hydrodynamic simulation of common envelope evolution to understand how energy is transferred between various forms and whether theory and simulation are mutually consistent given the setup. Virtually all of the envelope unbinding in the simulation occurs before the end of the rapid plunge-in phase, here defined to coincide with the first periastron passage. In contrast, the total envelope energy is nearly constant during this time because positive energy transferred to the gas from the core particles is counterbalanced by the negative binding energy from the closer proximity of the inner layers to the plunged-in secondary. During the subsequent slow spiral-in phase, energy continues to transfer to the envelope from the red giant core and secondary core particles. We also propose that relative motion between the centre of mass of the envelope and the centre of mass of the particles could account for the offsets of planetary nebula central stars from the nebula's geometric centre.
[10]  oai:arXiv.org:1811.09314  [pdf] - 1800001
Photoevaporative Flows From Exoplanet Atmospheres: A 3-D Radiative Hydrodynamic Parameter Study
Comments:
Submitted: 2018-11-22
The photoionization-driven evaporation of planetary atmospheres has emerged as a potentially fundamental process for planets on short period orbits. While 1-D studies have proven the effectiveness of stellar fluxes at altering the atmospheric mass and composition for sub-Jupiter mass planets, there remains much that is uncertain with regard to the larger-scale, multidimensional nature of such "planetary wind" flows. In this paper we use a new radiation-hydrodynamic platform to simulate atmospheric evaporative flows. Using the AstroBEAR AMR multiphysics code in a co-rotating frame centered on the planet, we model the transfer of ionizing photons into the atmosphere, the subsequent launch of the wind and the wind's large scale evolution subject to tidal and non-inertial forces. We run simulations for planets of 0.263 and 0.07 Jupiter masses and stellar fluxes of $2 \times 10^{13}$ and $2 \times 10^{14}$ photons/cm^2/s. Our results reveal new, potentially observable planetary wind flow patterns, including the development, in some cases, of an extended neutral tail lagging behind the planet in its orbit.
[11]  oai:arXiv.org:1810.04757  [pdf] - 2050118
Accretion in Common Envelope Evolution
Comments: 4 pages, 3 figures, to appear in the proceedings of IAU Symposium 343: Why Galaxies Care About AGB Stars, (eds.) F. Kerschbaum, M. Groenewegen and H. Olofsson
Submitted: 2018-10-10
Common envelope evolution (CEE) occurs in some binary systems involving asymptotic giant branch (AGB) or red giant branch (RGB) stars, and understanding this process is crucial for understanding the origins of various transient phenomena. CEE has been shown to be highly asymmetrical and global 3D simulations are needed to help understand the dynamics. We perform and analyze hydrodynamic CEE simulations with the adaptive mesh refinement (AMR) code AstroBEAR, and focus on the role of accretion onto the companion star. We bracket the range of accretion rates by comparing a model that removes mass and pressure using a subgrid accretion prescription with one that does not. Provided a pressure-release valve, such as a bipolar jet, is available, super-Eddington accretion could be common. Finally, we summarize new results pertaining to the energy budget, and discuss the overall implications relating to the feasibility of unbinding the envelope in CEE simulations.
[12]  oai:arXiv.org:1805.03607  [pdf] - 1724963
Accretion in common envelope evolution
Comments: 14 pages, submitted to MNRAS
Submitted: 2018-05-09
Common envelope evolution (CEE) is presently a poorly understood, yet critical, process in binary stellar evolution. Characterizing the full 3D dynamics of CEE is difficult in part because simulating CEE is so computationally demanding. Numerical studies have yet to conclusively determine how the envelope ejects and a tight binary results, if only the binary potential energy is used to propel the envelope. Additional power sources might be necessary and accretion onto the inspiraling companion is one such source. Accretion is likely common in post-asymptotic giant branch (AGB) binary interactions but how it operates and how its consequences depend on binary separation remain open questions. Here we use high resolution global 3D hydrodynamic simulations of CEE with the adaptive mesh refinement (AMR) code AstroBEAR, to bracket the range of CEE companion accretion rates by comparing runs that remove mass and pressure via a subgrid accretion model with those that do not. The results show that if a pressure release valve is available, super-Eddington accretion may be common. Jets are a plausible release valve in these environments, and they could also help unbind and shape the envelopes.
[13]  oai:arXiv.org:1805.00596  [pdf] - 1682727
Generation of a Circumstellar Gas Disk by Hot Jupiter WASP-12b
Comments:
Submitted: 2018-05-01
Observations of transiting extra-solar planets provide rich sources of data for probing the in-system environment. In the WASP-12 system, a broad depression in the usually-bright MgII h&k lines has been observed, in addition to atmospheric escape from the extremely hot Jupiter WASP-12b. It has been hypothesized that a translucent circumstellar cloud is formed by the outflow from the planet, causing the observed signatures. We perform 3D hydrodynamic simulations of the full system environment of WASP-12, injecting a planetary wind and stellar wind from their respective surfaces. We find that a torus of density high enough to account for the lack of MgII h&k line core emission in WASP-12 can be formed in approximately 13 years. We also perform synthetic observations of the Lyman-alpha spectrum at different points in the planet's orbit, which demonstrate that significant absorption occurs at all points in the orbit, not just during transits, as suggested by the observations.
[14]  oai:arXiv.org:1705.01998  [pdf] - 1583131
Wind-accelerated orbital evolution in binary systems with giant stars
Comments: 10 pages, 10 figures and 6 tables submitted to MNRAS
Submitted: 2017-05-04, last modified: 2017-08-31
Using 3D radiation-hydrodynamic simulations and analytic theory, we study the orbital evolution of asymptotic-giant-branch (AGB) binary systems for various initial orbital separations and mass ratios, and thus different initial accretion modes. The time evolution of binary separations and orbital periods are calculated directly from the averaged mass loss rate, accretion rate and angular momentum loss rate. We separately consider spin-orbit synchronized and zero spin AGB cases. We find that the the angular momentum carried away by the mass loss together with the mass transfer can effectively shrink the orbit when accretion occurs via wind-Roche-lobe overflow. In contrast, the larger fraction of mass lost in Bondi-Hoyle-Lyttleton accreting systems acts to enlarge the orbit. Synchronized binaries tend to experience stronger orbital period decay in close binaries. We also find that orbital period decay is faster when we account for the nonlinear evolution of the accretion mode as the binary starts to tighten. This can increase the fraction of binaries that result in common envelope, luminous red novae, Type Ia supernovae and planetary nebulae with tight central binaries. The results also imply that planets in the the habitable zone around white dwarfs are unlikely to be found.
[15]  oai:arXiv.org:1609.02918  [pdf] - 1580226
Reorienting MHD Colliding Flows: A Shock Physics Mechanism for Generating Filaments Normal to Magnetic Fields
Comments: In press, MNRAS, 11 pages, 13 figures
Submitted: 2016-09-09, last modified: 2017-06-07
We present numerical simulations of reorienting oblique shocks that form in the collision layer between magnetized colliding flows. Reorientation aligns parsec-scale post-shock filaments normal to the background magnetic field. We find that reorientation begins with pressure gradients between the collision region and the ambient medium. This drives a lateral expansion of post-shock gas, which reorients the growing filament from the outside-in (i.e. from the flow-ambient boundary, toward the colliding flows axis). The final structures of our simulations resemble polarization observations of filaments in Taurus and Serpens South, as well as the integral-shaped filament in Orion A. Given the ubiquity of colliding flows in the interstellar medium, shock reorientation may be relevant to the formation of filaments normal to magnetic fields.
[16]  oai:arXiv.org:1702.06160  [pdf] - 1562773
Mass Transfer and Disc Formation in AGB Binary Systems
Comments: 14 pages, 31 figures and 4 tables. accepted by MNRAS
Submitted: 2017-02-20, last modified: 2017-04-13
We investigate mass transfer and the formation of disc in binary systems using a combination of numerical simulations and theory. We consider six models distinguished by binary separation, secondary mass and outflow mechanisms. Each system consists of an asymptotic-giant-branch (AGB) star and an accreting secondary. The AGB star loses its mass via a wind. In one of our six models, the AGB star incurs a short period of outburst. In all cases, the secondary accretes part of the ejected mass and also influences the mass-loss rate of the AGB star. The ejected mass may remain gravitationally bound to the binary system and form a circumbinary disk, or contribute to an accretion disk around the secondary. In other cases, the ejecta will escape the binary system. The accretion rate on to the secondary changes non-linearly with binary separation. In our closest binary simulations, our models exemplify the wind Roche lobe overflow while in our wide binary cases, the mass transfer exhibits Bondi-Hoyle accretion. The morphologies of the outflows in the binary systems are varied. The variety may provide clues to how the late AGB phase influences planetary nebulae shaping. We employ the adaptive-mesh-refinement code ASTROBEAR for our simulations and include ray-tracing, radiation transfer, cooling and dust formation. To attain the highest computational efficiency and the most stable results, all simulations are run in the corotating frame.
[17]  oai:arXiv.org:1701.08764  [pdf] - 1581315
Mass transfer in asymptotic-giant-branch binary systems
Comments: 10 figures, 2 pages, revised for IAU323 proceedings
Submitted: 2017-01-30
Binary stars can interact via mass transfer when one member (the primary) ascends onto a giant branch. The amount of gas ejected by the binary and the amount of gas accreted by the secondary over the lifetime of the primary influence the subsequent binary phenomenology. Some of the gas ejected by the binary will remain gravitationally bound and its distribution will be closely related to the formation of planetary nebulae. We investigate the nature of mass transfer in binary systems containing an AGB star by adding radiative transfer to the AstroBEAR AMR Hydro/MHD code.
[18]  oai:arXiv.org:1602.01417  [pdf] - 1414144
Molecular cloud formation in high-shear, magnetized colliding flows
Comments: MNRAS, in press, Monthly Notices of the Royal Astronomical Society 2016
Submitted: 2016-02-01, last modified: 2016-05-30
The colliding flows (CF) model is a well-supported mechanism for generating molecular clouds. However, to-date most CF simulations have focused on the formation of clouds in the normal-shock layer between head-on colliding flows. We performed simulations of magnetized colliding flows that instead meet at an oblique-shock layer. Oblique shocks generate shear in the post-shock environment, and this shear creates inhospitable environments for star formation. As the degree of shear increases (i.e. the obliquity of the shock increases), we find that it takes longer for sink particles to form, they form in lower numbers, and they tend to be less massive. With regard to magnetic fields, we find that even a weak field stalls gravitational collapse within forming clouds. Additionally, an initially oblique collision interface tends to reorient over time in the presence of a magnetic field, so that it becomes normal to the oncoming flows. This was demonstrated by our most oblique shock interface, which became fully normal by the end of the simulation.
[19]  oai:arXiv.org:1604.08532  [pdf] - 1530645
Hot planetary winds near a star: dynamics, wind-wind interactions, and observational signatures
Comments: 17 pages, 14 figures
Submitted: 2016-04-28
Signatures of "evaporative" winds from exo-planets on short (hot) orbits around their host star have been observed in a number of systems. In this paper we present global AMR simulations that track the launching of the winds, their expansion through the circumstellar environment, and their interaction with a stellar wind. We focus on purely hydrodynamic flows including the anisotropy of the wind launching and explore the orbital/fluid dynamics of the resulting flows in detail. In particular we find that a combination of the tidal and Coriolis forces strongly distorts the planetary "Parker" wind creating "up-orbit" and "down-orbit" streams. We characterize the flows in terms of their orbital elements which change depending on their launch position on the planet. We find that the anisotropy in the atmospheric temperature leads to significant backflow on to the planet. The planetary wind interacts strongly with the stellar wind creating instabilities that cause eventual deposition of planetary gas onto the star. We present synthetic observations of both transit and absorption line-structure for our simulations. For our initial conditions, we find that the orbiting wind material produces absorption signatures at significant distances from the planet and substantial orbit to orbit variability. Ly-{\alpha} absorption shows red and blueshifted features out to 70 km/s. Finally, using semi-analytic models we constrain the effect of radiation pressure, given the approximation of uniform stellar absorption.
[20]  oai:arXiv.org:1603.09525  [pdf] - 1392974
The inner cavity of the circumnuclear disc
Comments: 18 pages, 11 figures, accepted by MNRAS
Submitted: 2016-03-31
The circumnuclear disc (CND) orbiting the Galaxy's central black hole is a reservoir of material that can ultimately provide energy through accretion, or form stars in the presence of the black hole, as evidenced by the stellar cluster that is presently located at the CND's centre. In this paper, we report the results of a computational study of the dynamics of the CND. The results lead us to question two paradigms that are prevalent in previous research on the Galactic Centre. The first is that the disc's inner cavity is maintained by the interaction of the central stellar cluster's strong winds with the disc's inner rim, and second, that the presence of unstable clumps in the disc implies that the CND is a transient feature. Our simulations show that, in the absence of a magnetic field, the interaction of the wind with the inner disc rim actually leads to a filling of the inner cavity within a few orbital time-scales, contrary to previous expectations. However, including the effects of magnetic fields stabilizes the inner disc rim against rapid inward migration. Furthermore, this interaction causes instabilities that continuously create clumps that are individually unstable against tidal shearing. Thus the occurrence of such unstable clumps does not necessarily mean that the disc is itself a transient phenomenon. The next steps in this investigation are to explore the effect of the magnetorotational instability on the disc evolution and to test whether the results presented here persist for longer time-scales than those considered here.
[21]  oai:arXiv.org:1508.00880  [pdf] - 1351464
Planetary Evaporation and the Dynamics of Planet Wind/Stellar Wind Bow Shocks
Comments: To appear in the proceedings of IAU Symposium 314 "Young Stars and Planets Near the Sun"
Submitted: 2015-08-04
We present initial results of a new campaign of simulations focusing on the interaction of planetary winds with stellar environments using Adaptive Mesh Refinement methods. We have confirmed the results of Stone & Proga 2009 that an azimuthal flow structure is created in the planetary wind due to day/night temperatures differences. We show that a backflow towards the planet will occur with a strength that depends on the escape parameter. When a stellar outflow is included, we see unstable bow waves forming through the outflow's interaction with the planetary wind.
[22]  oai:arXiv.org:1308.4198  [pdf] - 1173500
Size of discs formed by wind accretion in binaries can be underestimated if the role of wind-driving force is ignored
Comments: 8 pages, 2 figures, accepted by MNRAS; in press
Submitted: 2013-08-19, last modified: 2013-09-11
Binary systems consisting of a secondary accreting form a wind-emitting primary are ubiquitous in astrophysics. The phenomenology of such Bondi-Hoyle-Lyttleton (BHL) accretors is particularly rich when an accretion disc forms around the secondary. The outer radius of such discs is commonly estimated from the net angular momentum produced by a density variation of material across the BHL or Bondi accretion cylinder, as the latter is tilted with respect to the direction to the primary due to orbital motion. But this approach has ignored the fact that the wind experiences an outward driving force that the secondary does not. In actuality, the accretion stream falls toward a retarded point in the secondary's orbit as the secondary is pulled toward the primary relative to the stream. The result is a finite separation or "accretion stream impact parameter" (ASIP) separating the secondary and stream. When the orbital radius a_o exceeds the BHL radius r_b, the ratio of outer disc radius estimated as the ASIP to the conventional estimate a_o^{1/2}/r_b^{1/2}>1. We therefore predict that discs will form at larger radii from the secondary than traditional estimates. This agrees with the importance of the ASIP emphasized by Huarte-Espinosa et al. (2013) and the practical consequence that resolving the initial outer radius of such an accretion disc in numerical simulations can be less demanding than what earlier estimates would suggest.
[23]  oai:arXiv.org:1211.1672  [pdf] - 1157608
The Formation and Evolution of Wind-Capture Disks In Binary Systems
Comments: 13 pages, 9 figures, in press, MNRAS
Submitted: 2012-11-07, last modified: 2013-05-03
We study the formation, evolution and physical properties of accretion disks formed via wind capture in binary systems. Using the AMR code AstroBEAR, we have carried out high resolution 3D simulations that follow a stellar mass secondary in the co-rotating frame as it orbits a wind producing AGB primary. We first derive a resolution criteria, based on considerations of Bondi-Hoyle flows, that must be met in order to properly resolve the formation of accretion disks around the secondary. We then compare simulations of binaries with three different orbital radii (10, 15, 20 AU). Disks are formed in all three cases, however the size of the disk and, most importantly, its accretion rate decreases with orbital radii. In addition, the shape of the orbital motions of material within the disk becomes increasingly elliptical with increasing binary separation. The flow is mildly unsteady with "fluttering" around the bow shock observed. The disks are generally well aligned with the orbital plane after a few binary orbits. We do not observe the presence of any large scale, violent instabilities (such as the flip-flop mode). For the first time, moreover, it is observed that the wind component that is accreted towards the secondary has a vortex tube-like structure, rather than a column-like one as it was previously thought. In the context of AGB binary systems that might be precursors to Pre-Planetary and Planetary Nebula, we find that the wind accretion rates at the chosen orbital separations are generally too small to produce the most powerful outflows observed in these systems if the companions are main sequence stars but marginally capable if the companions are white dwarfs. It is likely that many of the more powerful PPN and PN involve closer binaries than the ones considered here. The results also demonstrate principles of broad relevance to all wind-capture binary systems.
[24]  oai:arXiv.org:1304.1367  [pdf] - 1055645
The Effects of Inhomogeneities within Colliding Flows on the Formation and Evolution of Molecular Clouds
Comments: Submitted to ApJ
Submitted: 2013-04-04, last modified: 2013-04-09
Observational evidence from local star-forming regions mandates that star formation occurs shortly after, or even during, molecular cloud formation. Models of the formation of molecular clouds in large-scale converging flows have identified the physical mechanisms driving the necessary rapid fragmentation. They also point to global gravitational collapse driving supersonic turbulence in molecular clouds. Previous cloud formation models have focused on turbulence generation, gravitational collapse, magnetic fields, and feedback. Here, we explore the effect of structure in the flow on the resulting clouds and the ensuing gravitational collapse. We compare two extreme cases, one with a collision between two smooth streams, and one with streams containing small clumps. We find that structured converging flows lead to a delay of local gravitational collapse ("star formation"). Thus, more gas has time to accumulate, eventually leading to a strong global collapse, and thus to a high star formation rate. Uniform converging flows fragment hydrodynamically early on, leading to the rapid onset of local gravitational collapse and an overall low sink formation rate.
[25]  oai:arXiv.org:1112.1710  [pdf] - 967235
Efficient Parallelization for AMR MHD Multiphysics Calculations; Implementation in AstroBEAR
Comments: Updated version of paper with improved scaling results
Submitted: 2011-12-07, last modified: 2013-02-21
Current Adaptive Mesh Refinement (AMR) simulations require algorithms that are highly parallelized and manage memory efficiently. As compute engines grow larger, AMR simulations will require algorithms that achieve new levels of efficient parallelization and memory management. We have attempted to employ new techniques to achieve both of these goals. Patch or grid based AMR often employs ghost cells to decouple the hyperbolic advances of each grid on a given refinement level. This decoupling allows each grid to be advanced independently. In AstroBEAR we utilize this independence by threading the grid advances on each level with preference going to the finer level grids. This allows for global load balancing instead of level by level load balancing and allows for greater parallelization across both physical space and AMR level. Threading of level advances can also improve performance by interleaving communication with computation, especially in deep simulations with many levels of refinement. While we see improvements of up to 30% on deep simulations run on a few cores, the speedup is typically more modest (5-20%) for larger scale simulations. To improve memory management we have employed a distributed tree algorithm that requires processors to only store and communicate local sections of the AMR tree structure with neighboring processors. Using this distributed approach we are able to get reasonable scaling efficiency (> 80%) out to 12288 cores and up to 8 levels of AMR - independent of the use of threading.
[26]  oai:arXiv.org:1110.1616  [pdf] - 422694
Efficient Parallelization for AMR MHD Multiphysics Calculations; Implementation in AstroBEAR
Comments: Proceedings From 6th Annual International Conference on Numerical Modeling of Space Plasma Flows
Submitted: 2011-10-07
Current AMR simulations require algorithms that are highly parallelized and manage memory efficiently. As compute engines grow larger, AMR simulations will require algorithms that achieve new levels of efficient parallelization and memory management. We have attempted to employ new techniques to achieve both of these goals. Patch or grid based AMR often employs ghost cells to decouple the hyperbolic advances of each grid on a given refinement level. This decoupling allows each grid to be advanced independently. In AstroBEAR we utilize this independence by threading the grid advances on each level with preference going to the finer level grids. This allows for global load balancing instead of level by level load balancing and allows for greater parallelization across both physical space and AMR level. Threading of level advances can also improve performance by interleaving communication with computation, especially in deep simulations with many levels of refinement. To improve memory management we have employed a distributed tree algorithm that requires processors to only store and communicate local sections of the AMR tree structure with neighboring processors.