Normalized to: Carroll-Nellenback, J.
[1]
oai:arXiv.org:1912.01647 [pdf] - 2107803
Bipolar Planetary Nebulae from Outflow Collimation by Common Envelope
Evolution
Zou, Yangyuxin;
Frank, Adam;
Chen, Zhuo;
Reichardt, Thomas;
De Marco, Orsola;
Blackman, Eric G.;
Nordhaus, Jason;
Balick, Bruce;
Carroll-Nellenback, Jonathan;
Chamandy, Luke;
Liu, Baowei
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?
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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.