Normalized to: Umurhan, O.

**41** article(s) in total. **278** co-authors, from **1** to **14** common article(s). Median position in authors list is **2,0**.

[1] oai:arXiv.org:2003.05576 [pdf] - 2063346

The solar nebula origin of (486958) Arrokoth, a primordial contact
binary in the Kuiper belt

Comments: Published in Science 28 Feb 2020 (First release 13 Feb 2020)

Submitted: **2020-03-11**

The New Horizons spacecraft's encounter with the cold classical Kuiper belt
object (486958) Arrokoth (formerly 2014 MU69) revealed a contact-binary
planetesimal. We investigate how it formed, finding it is the product of a
gentle, low-speed merger in the early Solar System. Its two lenticular lobes
suggest low-velocity accumulation of numerous smaller planetesimals within a
gravitationally collapsing, solid particle cloud. The geometric alignment of
the lobes indicates the lobes were a co-orbiting binary that experienced
angular momentum loss and subsequent merger, possibly due to dynamical friction
and collisions within the cloud or later gas drag. Arrokoth's contact-binary
shape was preserved by the benign dynamical and collisional environment of the
cold classical Kuiper belt, and so informs the accretion processes that
operated in the early Solar System.

[2] oai:arXiv.org:1906.05371 [pdf] - 2055773

Streaming instability in turbulent protoplanetary disks

Comments: Second revision, under final review for ApJ. New sections added. Old
section pertaining to nonlinear saturation removed. New title

Submitted: **2019-06-12**, last modified: **2020-02-26**

The streaming instability for solid particles in protoplanetary disks is
re-examined assuming the familiar alpha ($\alpha$) model for isotropic
turbulence. Turbulence always reduces the growth rates of the streaming
instability relative to values calculated for globally laminar disks. While for
small values of the turbulence parameter, $\alpha < 10^{-5}$, the wavelengths
of the fastest-growing disturbances are small fractions of the local gas
vertical scale height $H$, we find that for moderate values of the turbulence
parameter, i.e., $\alpha \sim 10^{-5}-10^{-3}$, the lengthscales of maximally
growing disturbances shift toward larger scales, approaching $H$. At these
moderate turbulent intensities and for local particle to gas mass density
ratios $\epsilon < 0.5$, the vertical scales of the most unstable modes begin
to exceed the corresponding radial scales so that the instability appears in
the form of vertically oriented sheets extending well beyond the particle scale
height. We find that for hydrodynamical turbulent disk models reported in the
literature, with { $\alpha = 4\times 10^{-5} - 5\times 10^{-4}$, } together
with state of the art global evolution models of particle growth, the streaming
instability is predicted to be viable within a narrow triangular patch of
$\alpha$--$\tau_s$ parameter space centered on Stokes numbers, $\tau_s \sim
0.01$ and $\alpha \sim 4\times 10^{-5}$ and, further, exhibits growth rates on
the order of several hundred to thousands of orbit times for disks with cosmic
solids abundance or metallicity. Our results are consistent with, and place in
context, published numerical studies of streaming instabilities.

[3] oai:arXiv.org:2002.06720 [pdf] - 2049475

Color, Composition, and Thermal Environment of Kuiper Belt Object
(486958) Arrokoth

Comments: 31 pages, 8 figures

Submitted: **2020-02-16**

The outer Solar System object (486958) Arrokoth (provisional designation 2014
MU$_{69}$) has been largely undisturbed since its formation. We study its
surface composition using data collected by the New Horizons spacecraft.
Methanol ice is present along with organic material, which may have formed
through radiation of simple molecules. Water ice was not detected. This
composition indicates hydrogenation of carbon monoxide-rich ice and/ or
energetic processing of methane condensed on water ice grains in the cold,
outer edge of the early Solar System. There are only small regional variations
in color and spectra across the surface, suggesting Arrokoth formed from a
homogeneous or well-mixed reservoir of solids. Microwave thermal emission from
the winter night side is consistent with a mean brightness temperature of
29$\pm$5 K.

[4] oai:arXiv.org:1909.03853 [pdf] - 2025784

Migrating Scarps as a Significant Driver for Cometary Surface Evolution

Comments: 14 pages, 4 figures, accepted GRL

Submitted: **2019-09-09**

Rosetta observations of 67P/Churyumov-Gerasimenko (67P) reveal that most
changes occur in the fallback-generated smooth terrains, vast deposits of
granular material blanketing the comet's northern hemisphere. These changes
express themselves both morphologically and spectrally across the nucleus, yet
we lack a model that describes their formation and evolution. Here we present a
self-consistent model that thoroughly explains the activity and mass loss from
Hapi's smooth terrains. Our model predicts the removal of dust via re-radiated
solar insolation localized within depression scarps that are substantially more
ice-rich than previously expected. We couple our model with numerous Rosetta
observations to thoroughly capture the seasonal erosion of Hapi's smooth
terrains, where local scarp retreat gradually removes the uppermost dusty
mantle. As sublimation-regolith interactions occur on rocky planets, comets,
icy moons and KBOs, our coupled model and observations provide a foundation for
future understanding of the myriad of sublimation-carved worlds.

[5] oai:arXiv.org:1903.05571 [pdf] - 1848689

Convection in a volatile nitrogen-ice-rich layer drives Pluto's
geological vigor

Comments:

Submitted: **2019-03-13**

The vast, deep, volatile-ice-filled basin informally named Sputnik Planum is
central to Pluto's geological activity[1,2]. Composed of molecular nitrogen,
methane, and carbon monoxide ices[3], but dominated by N2-ice, this ice layer
is organized into cells or polygons, typically ~10-40 km across, that resemble
the surface manifestation of solid state convection[1,2]. Here we report, based
on available rheological measurements[4], that solid layers of N2 ice
approximately greater than 1 km thick should convect for estimated present-day
heat flow conditions on Pluto. More importantly, we show numerically that
convective overturn in a several-km-thick layer of solid nitrogen can explain
the great lateral width of the cells. The temperature dependence of N2-ice
viscosity implies that the SP ice layer convects in the so-called sluggish lid
regime[5], a unique convective mode heretofore not definitively observed in the
Solar System. Average surface horizontal velocities of a few cm/yr imply
surface transport or renewal times of ~500,000 years, well under the 10 Myr
upper limit crater retention age for Sputnik Planum[2]. Similar convective
surface renewal may also occur on other dwarf planets in the Kuiper belt, which
may help explain the high albedos of some of them.

[6] oai:arXiv.org:1903.05574 [pdf] - 1848691

Reorientation of Sputnik Planitia implies a Subsurface Ocean on Pluto

Comments:

Submitted: **2019-03-13**

The deep nitrogen-covered Sputnik Planitia (SP; informal name) basin on Pluto
is located very close to the longitude of Pluto's tidal axis[1] and may be an
impact feature [2], by analogy with other large basins in the solar
system[3,4]. Reorientation[5-7] due to tidal and rotational torques can explain
SP's location, but requires it to be a positive gravity anomaly[7], despite its
negative topography. Here we argue that if SP formed via impact and if Pluto
possesses a subsurface ocean, a positive gravity anomaly would naturally result
because of shell thinning and ocean uplift, followed by later modest N2
deposition. Without a subsurface ocean a positive gravity anomaly requires an
implausibly thick N2 layer (greater than 40 km). A rigid, conductive ice shell
is required to prolong such an ocean's lifetime to the present day[8] and
maintain ocean uplift. Because N2 deposition is latitude-dependent[9], nitrogen
loading and reorientation may have exhibited complex feedbacks[7].

[7] oai:arXiv.org:1903.04546 [pdf] - 1847026

Planet formation: The case for large efforts on the computational side

Comments: White paper submitted to the Astro2020 decadal survey

Submitted: **2019-03-11**

Modern astronomy has finally been able to observe protoplanetary disks in
reasonable resolution and detail, unveiling the processes happening during
planet formation. These observed processes are understood under the framework
of disk-planet interaction, a process studied analytically and modeled
numerically for over 40 years. Long a theoreticians' game, the wealth of
observational data has been allowing for increasingly stringent tests of the
theoretical models. Modeling efforts are crucial to support the interpretation
of direct imaging analyses, not just for potential detections but also to put
meaningful upper limits on mass accretion rates and other physical quantities
in current and future large-scale surveys. This white paper addresses the
questions of what efforts on the computational side are required in the next
decade to advance our theoretical understanding, explain the observational
data, and guide new observations. We identified the nature of accretion, ab
initio planet formation, early evolution, and circumplanetary disks as major
fields of interest in computational planet formation. We recommend that
modelers relax the approximations of alpha-viscosity and isothermal equations
of state, on the grounds that these models use flawed assumptions, even if they
give good visual qualitative agreement with observations. We similarly
recommend that population synthesis move away from 1D hydrodynamics. The
computational resources to reach these goals should be developed during the
next decade, through improvements in algorithms and the hardware for hybrid
CPU/GPU clusters. Coupled with high angular resolution and great line
sensitivity in ground based interferometers, ELTs and JWST, these advances in
computational efforts should allow for large strides in the field in the next
decade.

[8] oai:arXiv.org:1903.03724 [pdf] - 1845998

Formation of Charon's Red Poles From Seasonally Cold-Trapped Volatiles

Comments:

Submitted: **2019-03-08**

A unique feature of Pluto's large satellite Charon is its dark red northern
polar cap. Similar colours on Pluto's surface have been attributed to organic
macromolecules produced by energetic radiation processing of hydrocarbons. The
polar location of this material on Charon implicates the temperature extremes
that result from Charon's high obliquity and long seasons. The escape of
Pluto's atmosphere provides a potential feed stock for production of complex
chemistry. Gas from Pluto that is transiently cold-trapped and processed at
Charon's winter pole was proposed as an explanation on the basis of an image of
Charon's northern hemisphere, but not modelled quantitatively. Here we report
images of the southern hemisphere illuminated by Pluto-shine and also images
taken during the approach phase showing the northern polar cap over a range of
longitudes. We model the surface thermal environment on Charon, the supply and
temporary cold-trapping of material escaping from Pluto, and, while
cold-trapped, its photolytic processing into more complex and less volatile
molecules. The model results are consistent with the proposed mechanism
producing the observed colour pattern on Charon.

[9] oai:arXiv.org:1903.02096 [pdf] - 1871512

The CH4 cycles on Pluto over seasonal and astronomical timescales

Comments: Accepted in Icarus

Submitted: **2019-03-05**

New Horizons observations suggest that CH4 on Pluto has a complex history,
involving reservoirs of different composition, thickness and stability
controlled by volatile processes occurring on different timescales. In order to
interpret these observations, we use a Pluto volatile transport model able to
simulate the cycles of N2 and CH4 ices over millions of years. By assuming
fixed solid mixing ratios, we explore how changes in surface albedos,
emissivities and thermal inertias impact volatile transport. This work is
therefore a direct and natural continuation of the work by Bertrand et al.
(2018), which only explored the N2 cycles. Results show that bright CH4
deposits can create cold traps for N2 ice outside Sputnik Planitia, leading to
a strong coupling between the N2 and CH4 cycles. Depending on the assumed
albedo for CH4 ice, the model predicts CH4 ice accumulation (1) at the same
equatorial latitudes where the Bladed Terrain Deposits are observed, supporting
the idea that these CH4-rich deposits are massive and perennial, or (2) at
mid-latitudes (25{\deg}N-70{\deg}N), forming a thick mantle which is consistent
with New Horizons observations. In our simulations, both CH4 ice reservoirs are
not in an equilibrium state and either one can dominate the other over long
timescales, depending on the assumptions made for the CH4 albedo. This suggests
that long-term volatile transport exists between the observed reservoirs. The
model also reproduces the formation of N2 deposits at mid-latitudes and in the
equatorial depressions surrounding the Bladed Terrain, as observed by New
Horizons. At the poles, only seasonal CH4 and N2 deposits are obtained in
Pluto's current orbital configuration. Finally, we show that Pluto's atmosphere
always contained, over the last astronomical cycles, enough gaseous CH4 to
absorb most of the incoming Lyman-flux.

[10] oai:arXiv.org:1808.08681 [pdf] - 1901742

The initial conditions for planet formation: Turbulence driven by
hydrodynamical instabilities in disks around young stars

Comments: Review to appear in Publications of the Astronomical Society of the
Pacific (PASP). 32 pages, 16 figures. Accepted

Submitted: **2018-08-26**, last modified: **2018-11-24**

This review examines recent theoretical developments in our understanding of
turbulence in cold, non-magnetically active, planetesimal forming regions of
protoplanetary disks which we refer to throughout as "Ohmic zones". We give a
brief background introduction to the subject of disk turbulence followed by a
terse pedagogical review of the phenomenology of hydrodynamic turbulence. The
equations governing the dynamics of cold astrophysical disks are given and
basic flow states are described. We discuss the Solberg-Hoiland conditions
required for stability, and the three recently identified turbulence generating
mechanisms possibly active in protoplanetary disk Ohmic zones, namely, (i) the
Vertical Shear Instability, (ii) The Convective Overstability and (iii) the
Zombie Vortex Instability. We summarize the properties of these processes,
identify their limitations and discuss where and under what conditions these
processes are active in protoplanetary disk Ohmic zones.

[11] oai:arXiv.org:1808.02118 [pdf] - 1729244

Great Expectations: Plans and Predictions for New Horizons Encounter
with Kuiper Belt Object 2014 MU69 ('Ultima Thule')

Comments:

Submitted: **2018-08-06**

The New Horizons encounter with the cold classical Kuiper Belt object (KBO)
2014 MU69 (informally named 'Ultima Thule,' hereafter Ultima) on 1 January 2019
will be the first time a spacecraft has ever closely observed one of the
free-orbiting small denizens of the Kuiper Belt. Related to but not thought to
have formed in the same region of the Solar System as the comets that been
explored so far, it will also be the largest, most distant, and most primitive
body yet visited by spacecraft. In this letter we begin with a brief overview
of cold classical KBOs, of which Ultima is a prime example. We give a short
preview of our encounter plans. We note what is currently known about Ultima
from earth-based observations. We then review our expectations and capabilities
to evaluate Ultima's composition, surface geology, structure, near space
environment, small moons, rings, and the search for activity.

[12] oai:arXiv.org:1606.05734 [pdf] - 1722878

Modeling glacial flow on and onto Pluto's Sputnik Planitia

Comments: This definitive version has several typos corrected including
equation expressions in Appendix C.2 and C.3. Updated version with associated
errata in process at Icarus (July 26, 2018). All previous results appearing
in journal version --especially those pertaining to glacial flow and their
related timescales (Sections 4 & 5) -- are unchanged

Submitted: **2016-06-18**, last modified: **2018-07-26**

Observations of Pluto's surface made by the New Horizons spacecraft indicates
present-day nitrogen ice glaciation in and around the basin known as Sputnik
Planum. Motivated by these observations, we have developed an evolutionary
glacial flow model of solid nitrogen ice taking into account its published
thermophysical and rheologies properties. This model assumes that glacial ice
layers flow laminarly and have low aspect ratios which permits a vertically
integrated mathematical formulation. We assess the conditions for the validity
of laminar nitrogen ice motion by revisiting the problem of the onset of
solid-state buoyant convection of nitrogen ice for a variety of bottom thermal
boundary conditions. Subject to uncertainties in nitrogen ice rheology,
nitrogen ice layers are estimated to flow laminarly for thicknesses less than
400-1000 meters. The resulting mass-flux formulation for when the nitrogen ice
flows as a laminar dry glacier is characterized by an Arrhenius-Glen functional
form. The flow model developed is used here to qualitatively answer some
questions motivated by observed glacial flow features found on Sputnik Planum.
We find that the wavy transverse dark features found along the northern
shoreline of Sputnik Planum may be a transitory imprint of shallow topography
just beneath the ice surface suggesting the possibility that a major shoreward
flow event happened relatively recently within the last few hundred years.
Model results also support the interpretation that the prominent darkened
features resembling flow lobes observed along the eastern shoreline of the
Sputnik Planum basin may be a result of wet nitrogen glacial ice flowing into
the basin from the pitted highlands of eastern Tombaugh Regio.

[13] oai:arXiv.org:1804.02434 [pdf] - 1666894

The nitrogen cycles on Pluto over seasonal and astronomical timescales

Comments: 49 pages, 16 figures

Submitted: **2018-04-06**

Pluto's landscape is shaped by the cycles of the volatile ices covering its
surface. In particular, the Sputnik Planitia (SP) ice sheet displays a large
diversity of terrains, with bright and dark plains, pits, topographic
depressions and evidences of recent and past glacial flows. Outside SP, New
Horizons also revealed numerous N2 ice deposits, in Tombaugh Regio and at
mid-northern latitudes. These observations suggest a complex history involving
volatile and glacial processes on different timescales. We present numerical
simulations of volatile transport on Pluto performed with a model able to
simulate the N2 cycle over millions of years (Myrs), taking into account the
changes of obliquity and orbital parameters as experienced by Pluto. Results
show that over one obliquity cycle, the latitudes of SP between
25{\deg}S-30{\deg}N are dominated by N2 condensation, while the latitudes
between 30-50{\deg}N are dominated by N2 sublimation. We find that a net amount
of 1 km of ice has sublimed at the northern edge of SP during the last 2 Myrs.
By comparing these results with the observed geology of SP, we can relate the
formation of the pits and the brightness of the ice to the ice flux occurring
at the annual timescale, while the glacial flows at its eastern edge and the
erosion of the water ice mountains all around the ice sheet are related to the
astronomical timescale. We also perform simulations with a glacial flow scheme
which shows that SP is currently at its minimum extent. We also explore the
stability of N2 ice outside SP. Results show that it is not stable at the poles
but rather in the equatorial regions, in particular in depressions, where thick
deposits may persist over tens of Myrs, before being trapped in SP. Finally,
another key result is that the minimum and maximum surface pressures obtained
over the simulated Myrs remain in the range of mm-Pa and Pa, respectively.

[14] oai:arXiv.org:1601.00382 [pdf] - 1501895

Critical layers and protoplanetary disk turbulence

Comments: Accepted to ApJ (25 pages). Results of previously appearing versions
unchanged. Changes and additions include: updated references and
nomenclature, viscous solutions added, higher resolution solutions displayed,
expanded discussion, and shortened abstract

Submitted: **2016-01-04**, last modified: **2016-07-09**

A linear analysis of the zombie vortex instability is performed in a
stratified shearing sheet setting for three model barotropic shear flows. The
linear analysis is done by utilizing a Green's function formulation to resolve
the critical layers of the associated normal-mode problem. The instability is
the result of a resonant interaction between a Rossby wave and a gravity wave
which we refer to as Z-modes. The associated critical layer is the location
where the Doppler shifted frequency of a distant Rossby wave equals the local
Brunt-Vaisala frequency. The minimum required Rossby number for instability,
Ro= 0.2, is confirmed for parameter values reported in the literature. It is
also found that the shear layer supports the instability in the limit where
stratification vanishes. The zombie vortex instability is examined in a jet
model, finding that the instability can occur for Ro= 0.05. Nonlinear vorticity
forcing due to unstable Z-modes is shown to result in the creation of a jet
flow at the critical layer emerging as the result of the competition between
the vertical lifting of perturbation radial vorticity and the radial transport
of perturbation vertical vorticity. We find that the picture of this
instability leading to a form of nonlinearly driven self-replicating pattern of
creation and destruction is warranted: a parent jet spawns a growing child jet
at associated critical layers. A mature child jet creates a next generation of
child jets at associated critical layers while simultaneously destroying its
parent jet via the Rossby wave instability.

[15] oai:arXiv.org:1603.00821 [pdf] - 1542656

Mean radius and shape of Pluto and Charon from New Horizons images

Comments: revised version submitted to Icarus (Pluto special issue)

Submitted: **2016-03-02**, last modified: **2016-06-22**

Approach images taken by the LORRI imaging system during the New Horizons
spacecraft encounter have been used to determine the mean radii and shapes of
Pluto and Charon. The primary observations are limb locations derived using
three independent approaches. The resulting mean radii of Pluto and Charon are
1188.3 +/- 1.6 km and 606.0 +/- 1.0 km, respectively (2-sigma). The
corresponding densities are 1854 +/- 11 kg m-3 and 1701 +/- 33 kg m-3
(2-sigma). The Charon radius value is consistent with previous Earth-based
occultation estimates. The Pluto radius estimate is consistent with solar
occultation measurements performed by the ALICE and Fine Sun Sensor instruments
on New Horizons. Neither Pluto nor Charon show any evidence for
tidal/rotational distortions; upper bounds on the oblateness are <0.6% and
<0.5%, respectively.

[16] oai:arXiv.org:1604.05702 [pdf] - 1393089

The Geology of Pluto and Charon Through the Eyes of New Horizons

Comments:

Submitted: **2016-04-19**

NASA's New Horizons spacecraft has revealed the complex geology of Pluto and
Charon. Pluto's encounter hemisphere shows ongoing surface geological activity
centered on a vast basin containing a thick layer of volatile ices that appears
to be involved in convection and advection, with a crater retention age no
greater than $\approx$10 Ma. Surrounding terrains show active glacial flow,
apparent transport and rotation of large buoyant water-ice crustal blocks, and
pitting, likely by sublimation erosion and/or collapse. More enigmatic features
include tall mounds with central depressions that are conceivably cryovolcanic,
and ridges with complex bladed textures. Pluto also has ancient cratered
terrains up to ~4 Ga old that are extensionally fractured and extensively
mantled and perhaps eroded by glacial or other processes. Charon does not
appear to be currently active, but experienced major extensional tectonism and
resurfacing (probably cryovolcanic) nearly 4 billion years ago. Impact crater
populations on Pluto and Charon are not consistent with the steepest proposed
impactor size-frequency distributions proposed for the Kuiper belt.

[17] oai:arXiv.org:1503.08470 [pdf] - 1396662

On the mechanism of self gravitating Rossby interfacial waves in
proto-stellar accretion discs

Comments: 20 pages, 5 figures and in GAFD article sytle

Submitted: **2015-03-29**, last modified: **2016-01-12**

The dynamical response of edge waves under the influence of self-gravity is
examined in an idealized two-dimensional model of a proto-stellar disc,
characterized in steady state as a rotating vertically infinite cylinder of
fluid with constant density except for a single density interface at some
radius r0. The fluid in basic state is prescribed to rotate with a Keplerian
profile $\Omega_k(r)\sim r^{-3/2}$ modified by some additional azimuthal
sheared flow. A linear analysis shows that there are two azimuthally
propagating edge waves, kin to the familiar Rossby waves and surface gravity
waves in terrestrial studies, which move opposite to one another with respect
to the local basic state rotation rate at the interface. Instability only
occurs if the radial pressure gradient is opposite to that of the density jump
(unstably stratified) where self-gravity acts as a wave stabilizer irrespective
of the stratification of the system. The propagation properties of the waves
are discussed in detail in the language of vorticity edge waves. The roles of
both Boussinesq and non-Boussinesq effects upon the stability and propagation
of these waves with and without the inclusion of self-gravity are then
quantified. The dynamics involved with self-gravity non- Boussinesq effect is
shown to be a source of vorticity production where there is a jump in the basic
state density, in addition, self-gravity also alters the dynamics via the
radial main pressure gradient, which is a Boussinesq effect . Further
applications of these mechanical insights are presented in the conclusion
including the ways in which multiple density jumps or gaps may or may not be
stable.

[18] oai:arXiv.org:1601.01921 [pdf] - 1365574

Vortex formation in protoplanetary discs induced by the vertical shear
instability

Comments:

Submitted: **2016-01-08**

We present the results of 2D and 3D hydrodynamic simulations of idealized
protoplanetary discs that examine the formation and evolution of vortices by
the vertical shear instability (VSI). In agreement with recent work, we find
that discs with radially decreasing temperature profiles and short thermal
relaxation time-scales, are subject to the axisymmetric VSI. In three
dimensions, the resulting velocity perturbations give rise to
quasi-axisymmetric potential vorticity perturbations that break-up into
discrete vortices, in a manner that is reminiscent of the Rossby wave
instability. Discs with very short thermal evolution time-scales (i.e.
{\tau}<0.1 local orbit periods) develop strong vorticity perturbations that
roll up into vortices that have small aspect ratios ({\chi}<2) and short
lifetimes (~ a few orbits). Longer thermal time-scales give rise to vortices
with larger aspect ratios (6<{\chi}<10), and lifetimes that depend on the
entropy gradient. A steeply decreasing entropy profile leads to vortex
lifetimes that exceed the simulation run times of hundreds of orbital periods.
Vortex lifetimes in discs with positive or weakly decreasing entropy profiles
are much shorter, being 10s of orbits at most, suggesting that the subcritical
baroclinic instability plays an important role in sustaining vortices against
destruction through the elliptical instability. Applied to the outer regions of
protoplanetary discs, where the VSI is most likely to occur, our results
suggest that vortices formed by the VSI are likely to be short lived
structures.

[19] oai:arXiv.org:1510.07704 [pdf] - 1300465

The Pluto system: Initial results from its exploration by New Horizons

Comments: 8 pages - Initial Science paper from NASA's New Horizons Pluto
Encounter

Submitted: **2015-10-26**

The Pluto system was recently explored by NASA's New Horizons spacecraft,
making closest approach on 14 July 2015. Pluto's surface displays diverse
landforms, terrain ages, albedos, colors, and composition gradients. Evidence
is found for a water-ice crust, geologically young surface units, surface ice
convection, wind streaks, volatile transport, and glacial flow. Pluto's
atmosphere is highly extended, with trace hydrocarbons, a global haze layer,
and a surface pressure near 10 microbars. Pluto's diverse surface geology and
long-term activity raise fundamental questions about how small planets remain
active many billions of years after formation. Pluto's large moon Charon
displays tectonics and evidence for a heterogeneous crustal composition, its
north pole displays puzzling dark terrain. Small satellites Hydra and Nix have
higher albedos than expected.

[20] oai:arXiv.org:1505.01892 [pdf] - 1351423

Linear analysis of the vertical shear instability: outstanding issues
and improved solutions (Research Note)

Comments: 8 pages, 5 figures, submitted to A+A as a Research Note (Version 2
reflects typo corrections to some formulae - results unaltered from previous
version - still under review)

Submitted: **2015-05-07**, last modified: **2015-05-28**

The Vertical Shear Instability is one of two known mechanisms potentially
active in the so-called dead zones of protoplanetary accretion disks. A recent
analysis indicates that a subset of unstable modes shows unbounded growth -
both as resolution is increased and when the nominal lid of the atmosphere is
extended, possibly indicating ill-posedness in previous attempts of linear
analysis. The reduced equations governing the instability are revisited and the
generated solutions are examined using both the previously assumed separable
forms and an improved non-separable solution form that is herewith introduced.
Analyzing the reduced equations using the separable form shows that, while the
low-order body modes have converged eigenvalues and eigenfunctions (as both the
vertical boundaries of the atmosphere are extended and with increased radial
resolution), it is also confirmed that the corresponding high-order body modes
and the surface modes do indeed show unbounded growth rates. However, the
energy contained in both the higher-order body modes and surface modes
diminishes precipitously due to the disk's Gaussian density profile. Most of
the energy of the instability is contained in the low-order modes. An
inseparable solution form is introduced which filters out the inconsequential
surface modes leaving only body modes (both low and high-order ones). The
analysis predicts a fastest growing mode with a specific radial length scale.
The growth rates associated with the fundamental corrugation and breathing
modes matches the growth and length scales observed in previous nonlinear
studies of the instability.

[21] oai:arXiv.org:1311.5134 [pdf] - 1201760

Intercomparison of General Circulation Models for Hot Extrasolar Planets

Comments: 29 pages, 11 figures: Accepted for publication in Icarus

Submitted: **2013-11-20**

We compare five general circulation models (GCMs) which have been recently
used to study hot extrasolar planet atmospheres (BOB, CAM, IGCM, MITgcm, and
PEQMOD), under three test cases useful for assessing model convergence and
accuracy. Such a broad, detailed intercomparison has not been performed thus
far for extrasolar planets study. The models considered all solve the
traditional primitive equations, but employ different numerical algorithms or
grids (e.g., pseudospectral and finite volume, with the latter separately in
longitude-latitude and `cubed-sphere' grids). The test cases are chosen to
cleanly address specific aspects of the behaviors typically reported in hot
extrasolar planet simulations: 1) steady-state, 2) nonlinearly evolving
baroclinic wave, and 3) response to fast timescale thermal relaxation. When
initialized with a steady jet, all models maintain the steadiness, as they
should -- except MITgcm in cubed-sphere grid. A very good agreement is obtained
for a baroclinic wave evolving from an initial instability in pseudospectral
models (only). However, exact numerical convergence is still not achieved
across the pseudospectral models: amplitudes and phases are observably
different. When subject to a typical `hot-Jupiter'-like forcing, all five
models show quantitatively different behavior -- although qualitatively
similar, time-variable, quadrupole-dominated flows are produced. Hence, as have
been advocated in several past studies, specific quantitative predictions (such
as the location of large vortices and hot regions) by GCMs should be viewed
with caution. Overall, in the tests considered here, pseudospectral models in
pressure coordinate (PEBOB and PEQMOD) perform the best and MITgcm in
cubed-sphere grid performs the worst.

[22] oai:arXiv.org:1301.0285 [pdf] - 609612

The Equations of Magnetoquasigeostrophy

Comments: 13 pages. Under consideration for publication in Astronomy and
Astrophysics

Submitted: **2013-01-02**

The dynamics contained in magnetized layers of exoplanet atmospheres are
important to understand in order to characterize what observational signatures
they may provide for future observations. It is important to develop a
framework to begin studying and learning the physical processes possible under
those conditions and what, if any, features contained in them may be observed
in future observation missions. The aims of this study is to formally derive,
from scaling arguments, a manageable reduced set of equations for analysis,
i.e. a magnetic formulation of the equations of quasigeostrophy appropriate for
a multi-layer atmosphere. The main goal is to provide a simpler theoretical
platform to explore the dynamics possible within confined magnetized layers of
exoplanet atmospheres. We primarily use scaling arguments to derive the reduced
equations of "magnetoquasigeostrophy" which assumes dynamics to take place in
an atmospheric layer which is vertically thin compared to its horizontal
scales. The derived equation set retains features existing in standard
shallow-water magnetohydrodynamic equations but are absent in more classical
derivations of the quasi-geostrophic limit, namely, the non-divergence of the
in-plane components of the magnetic field. We liken this non-divergence of the
in-plane magnetic fields as indicative of a quantity whose behaviour mimics a
two-dimensional "pseudo"-magnetic monopole source. We also find, using the same
scaling argument procedures, appropriate limits of the fundamental parameters
of the system which yield reduced equations describing the flow dynamics
primarily characterized by magnetostrophic balance. The standard scaling
arguments employed here show how traditional magnetized quasigeostrophic
equations connect to their magnetized shallow water forms. The equations
derived are amenable to analysis using well-known techniques.

[23] oai:arXiv.org:1210.5343 [pdf] - 1152283

Non-Dissipative Saturation of the Magnetorotational Instability in Thin
Disks

Comments: 4 pages, 3 figures; Phys. Rev. Lett. accepted

Submitted: **2012-10-19**

A new non-dissipative mechanism is proposed for the saturation of the
axisymmetric magnetorotational (MRI) instability in thin Keplerian disks that
are subject to an axial magnetic field. That mechanism relies on the energy
transfer from the MRI to stable magnetosonic (MS) waves. Such mode interaction
is enabled due to the vertical stratification of the disk that results in the
discretization of its MRI spectrum, as well as by applying the appropriate
boundary conditions. A second order Duffing-like amplitude equation for the
initially unstable MRI modes is derived. The solutions of that equation exhibit
bursty nonlinear oscillations with a constant amplitude that signifies the
saturation level of the MRI. Those results are verified by a direct numerical
solution of the full nonlinear reduced set of thin disk magnetohydrodynamics
equations.

[24] oai:arXiv.org:1209.2753 [pdf] - 1151374

Linear and nonlinear evolution of the vertical shear instability in
accretion discs

Comments: 19 pages, 15 figures, submitted to MNRAS for publication. A version
with full resolution figures is available at
http://www.maths.qmul.ac.uk/~rpn/Corrugation/NGU.pdf

Submitted: **2012-09-12**

(Abridged) We analyse the stability and evolution of power-law accretion disc
models. These have midplane densities that follow radial power-laws, and have
either temperature or entropy distributions that are power-law functions of
cylindrical radius. We employ two different hydrodynamic codes to perform
2D-axisymmetric and 3D simulations that examine the long-term evolution of the
disc models as a function of the power-law indices of the temperature or
entropy, the thermal relaxation time of the fluid, and the viscosity. We
present a stability analysis of the problem that we use to interpret the
simulation results. We find that disc models whose temperature or entropy
profiles cause the equilibrium angular velocity to vary with height are
unstable to the growth of modes with wavenumber ratios |k_R/k_Z| >> 1 when the
thermodynamic response of the fluid is isothermal, or the thermal evolution
time is comparable to or shorter than the local dynamical time scale. These
discs are subject to the Goldreich-Schubert-Fricke (GSF) or `vertical shear'
linear instability. Development of the instability involves excitation of
vertical breathing and corrugation modes in the disc, with the corrugation
modes in particular being a feature of the nonlinear saturated state.
Instability operates when the dimensionless disc kinematic viscosity nu <
10^{-6} (Reynolds numbers Re>H c_s/nu > 2500). In 3D the instability generates
a quasi-turbulent flow, and the Reynolds stress produces a fluctuating
effective viscosity coefficient whose mean value reaches alpha ~ 6 x 10^{-4} by
the end of the simulation. The vertical shear instability in disc models which
include realistic thermal physics has yet to be examined. Should it occur,
however, our results suggest that it will have significant consequences for
their internal dynamics, transport properties, and observational appearance.

[25] oai:arXiv.org:1203.1686 [pdf] - 1117152

Potential vorticity dynamics in the framework of disk shallow-water
theory: II. Mixed Barotropic-Baroclinic Instability

Comments: Under consideration for publication in Astronomy and Astrophysics.
Comments to the author welcome

Submitted: **2012-03-07**

We extend exploration of potential vorticity instabilities in cold
astrophysical disks whose mean states are baroclinic. In particular, we seek to
demonstrate the potential existence of traditional baroclinic instabilities of
meteorological studies in a simplified two-layer Philips Disk Model. Each disk
layer is of constant but differing densities. The resulting mean azimuthal
velocity profile shows a variation in the vertical direction implying that the
system is baroclinic in the mean state. The stability of the system is treated
in the context of disk shallow water theory wherein azimuthal disturbances are
much longer than the corresponding radial or vertical scales. The normal-mode
problem is solved numerically using two different methods. The results of a
symmetric single layer barotropic model is considered and it is found that
instability persists for models in which the potential vorticity profiles are
not symmetric, consistent with previous results. The instaiblity is interpreted
in terms of interacting Rossby waves. For a two layer system in which the flow
is fundamentally baroclinic we report here that instability takes on the form
of mixed barotropic-baroclinic type: instability occurs but it qualitatively
follows the pattern of instability found in the barotropic models. Instability
arises because of the phase locking and interaction of the Rossby waves between
the two layers. The strength of the instability weakens as the density contrast
between layers increases. (For full abstract see article.)

[26] oai:arXiv.org:1008.2073 [pdf] - 1034237

Potential vorticity dynamics in the framework of disk shallow-water
theory: I. The Rossby wave instability

Comments: To appear in Astronomy and Astrophysics

Submitted: **2010-08-12**

The Rossby wave instability in astrophysical disks is as a potentially
important mechanism for driving angular momentum transport in disks. We aim to
understand this instability in an approximate three-dimensional disk model
environment which we assume to be a single homentropic annular layer we analyze
using disk shallow-water theory. We consider the normal mode stability analysis
of two kinds of radial profiles of the mean potential vorticity: The first type
is a single step and the second kind is a symmetrical step of finite width
describing either a localized depression or peak of the mean potential
vorticity. For single potential vorticity steps we find there is no
instability. There is no instability when the symmetric step is a localized
peak. However, the Rossby wave instability occurs when the symmetrical step
profile is a depression, which, in turn, corresponds to localized peaks in the
mean enthalpy profile. This is in qualitative agreement with previous
two-dimensional investigations of the instability. For all potential vorticity
depressions, instability occurs for regions narrower than some maximum radial
length scale. We interpret the instability as resulting from the interaction of
at least two Rossby edgewaves. We identify the Rossby wave instability in the
restricted three-dimensional framework of disk shallow water theory. Additional
examinations of generalized barotropic flows are needed. Viewing disk vortical
instabilities from the conceptual perspective of interacting edgewaves can be
useful.

[27] oai:arXiv.org:1002.1719 [pdf] - 1025060

Non-exponential hydrodynamical growth in density-stratified thin
Keplerian discs

Comments:

Submitted: **2010-02-08**, last modified: **2010-03-12**

The short time evolution of three dimensional small perturbations is studied.
Exhibiting spectral asymptotic stability, thin discs are nonetheless shown to
host intensive hydrodynamical activity in the shape of non modal growth of
initial small perturbations. Two mechanisms that lead to such behavior are
identified and studied, namely, non-resonant excitation of vertically confined
sound waves by stable planar inertia-coriolis modes that results in linear
growth with time, as well as resonant coupling of those two modes that leads to
a quadratic growth of the initial perturbations. It is further speculated that
the non modal growth can give rise to secondary strato-rotational instabilities
and thus lead to a new route to turbulence generation in thin discs.

[28] oai:arXiv.org:1002.0506 [pdf] - 1024944

Low magnetic-Prandtl number flow configurations for cold astrophysical
disk models: speculation and analysis

Comments: To Appear in Astronomy and Astrophysics

Submitted: **2010-02-02**

Simulations of astrophysical disks in the shearing box that are subject to
the magnetorotational instability (MRI) show that activity appears to be
reduced as the magnetic Prandtl number (Pm) is lowered. On the other hand,
calculations for laboratory experiments show that saturation is achieved
through modification of the background shear for Pm << 1. Guided by the results
of calculations appropriate for laboratory experiments when Pm is very low, the
axisymmetric stability of inviscid disturbances in a shearing box model
immersed in a constant vertical background magnetic field is considered under a
variety of shear profiles and boundary conditions in order to evaluate the
hypothesis that modifications of the shear bring about saturation of the
instability. It is found that the emergence and stability of the MRI is
sensitive to the boundary conditions adopted. Channel modes do not appear to be
stabilized through any modification of the background shear whose average
remains Keplerian. However, systems that have non-penetrative boundaries can
saturate the MRI through modification of the background shear. Conceptually
equating the qualitative results from laboratory experiments to the conditions
in a disk may therefore be misleading.

[29] oai:arXiv.org:0906.0004 [pdf] - 1002403

Global transient dynamics of three-dimensional hydrodynamical
disturbances in a thin viscous accretion disk

Comments: Accepted in Physics of Fluids

Submitted: **2009-06-01**, last modified: **2009-06-11**

Thin viscous Keplerian accretion disks are considered asymptotically stable,
even though they can show significant dynamic activity on short timescales. In
this paper the dynamics of non-axisymmetric hydrodynamical disturbances of
disks are investigated analytically building upon the steady state
three-dimensional structure and evolution of axisymmetric perturbations
explored in previous work. Assuming a polytropic equation of state solutions
are found by means of an asymptotic expansion in the small parameter measuring
the ratio of the disk thickness to characteristic radius. In-depth analysis
shows that every perturbation that disturbs the radial velocity induces
significant transient growth in the (acoustic) energy of the evolving
disturbance. This effect is most evident in the density and vertical velocity.
The transient growth observed is tied to the non-separable nature of the
solutions where, in particular, pattern evolution is controlled by a similarity
variable composed of the radial coordinate and time. This leads to growing
winding perturbations that display successive radial peaks and troughs. We
argue that these transient non-axisymmetric structures may precipitate
secondary instabilities which, consequently, may be a critical element for a
new alternative picture of turbulence arousal in non-magnetized astrophysical
disks.

[30] oai:arXiv.org:0810.5498 [pdf] - 18022

Linear dynamics of weakly viscous accretion disks: A disk analog of
Tollmien-Schlichting waves

Comments: Accepted for publication in Astronomy and Astrophysics. 18 pages.
This version 3 with corrected style file

Submitted: **2008-10-30**, last modified: **2008-11-04**

This paper discusses new perspectives and approaches to the problem of disk
dynamics where, in this study, we focus on the effects of viscous instabilities
influenced by boundary effects. The Boussinesq approximation of the viscous
large shearing box equations is analyzed in which the azimuthal length scale of
the disturbance is much larger than the radial and vertical scales. We examine
the stability of a non-axisymmetric potential vorticity mode, i.e. a
PV-anomaly. in a configuration in which buoyant convection and the
strato-rotational instability do not to operate. We consider a series of
boundary conditions which show the PV-anomaly to be unstable both on a finite
and semi-infinite radial domains. We find these conditions leading to an
instability which is the disk analog of Tollmien-Schlichting waves. When the
viscosity is weak, evidence of the instability is most pronounced by the
emergence of a vortex sheet at the critical layer located away from the
boundary where the instability is generated. For some boundary conditions a
necessary criterion for the onset of instability for vertical wavelengths that
are a sizable fraction of the layer's thickness and when the viscosity is small
is that the appropriate Froude number of the flow be greater than one. This
instability persists if more realistic boundary conditions are applied,
although the criterion on the Froude number is more complicated. The unstable
waves studied here share qualitative features to the instability seen in
rotating Blasius boundary layers. The implications of these results are
discussed. An overall new strategy for exploring and interpreting disk
instability mechanisms is also suggested.

[31] oai:arXiv.org:0802.3486 [pdf] - 10361

A shallow-water theory for annular sections of Keplerian Disks

Comments: Accepted (July 21, 2008), now in final form

Submitted: **2008-02-23**, last modified: **2008-09-22**

A scaling argument is presented that leads to a shallow water theory of
non-axisymmetric disturbances in annular sections of thin Keplerian disks. To
develop a theoretical construction that will aid in physically understanding
the relationship of known two-dimensional vortex dynamics to their
three-dimensional counterparts in Keplerian disks. Using asymptotic scaling
arguments varicose disturbances of a Keplerian disk are considered on radial
and vertical scales consistent with the height of the disk while the azimuthal
scales are the full $2\pi$ angular extent of the disk. The scalings lead to
dynamics which are radially geostrophic and vertically hydrostatic. It follows
that a potential vorticity quantity emerges and is shown to be conserved in a
Lagrangian sense. Uniform potential vorticity linear solutions are explored and
the theory is shown to contain an incarnation of the strato-rotational
instability under channel flow conditions. Linearized solutions of a single
defect on an infinite domain is developed and is shown to support a propagating
Rossby edgewave. Linear non-uniform potential vorticity solutions are also
developed and are shown to be similar in some respects to the dynamics of
strictly two-dimensional inviscid flows. Based on the framework of this theory,
arguments based on geophysical notions are presented to support the assertion
that the strato-rotational instability is in a generic class of
barotropic/baroclinic potential vorticity instabilities. Extensions of this
formalism are also proposed. The shallow water formulation achieved by the
asymptotic theory developed here opens a new approach to studying disk
dynamics.

[32] oai:arXiv.org:astro-ph/0703358 [pdf] - 90167

Hydrodynamic response of rotationally supported flows in the Small
Shearing Box model

Comments: 6 pages, resubmitted to Astronomy and Astrophysics, shortened with
references added

Submitted: **2007-03-14**, last modified: **2008-03-04**

The hydrodynamic response of the inviscid small shearing box model of a
midplane section of a rotationally supported astrophysical disk is examined. An
energy functional ${\cal E}$ is formulated for the general nonlinear problem.
It is found that the fate of disturbances is related to the conservation of
this quantity which, in turn, depends on the boundary conditions utilized:
${\cal E}$ is conserved for channel boundary conditions while it is not
conserved in general for shearing box conditions. Linearized disturbances
subject to channel boundary conditions have normal-modes described by Bessel
Functions and are qualitatively governed by a quantity $\Sigma$ which is a
measure of the ratio between the azimuthal and vertical wavelengths. Inertial
oscillations ensue if $\Sigma >1$ - otherwise disturbances must in general be
treated as an initial value problem. We reflect upon these results and offer a
speculation.

[33] oai:arXiv.org:0711.0794 [pdf] - 6743

On the viability of the shearing box approximation for numerical studies
of MHD turbulence in accretion disks

Comments: Accepted for publication in Astronomy and Astrophysics

Submitted: **2007-11-05**, last modified: **2008-02-27**

Most of our knowledge on the nonlinear development of the magneto-rotational
instability (MRI) relies on the results of numerical simulations employing the
shearing box (SB) approximation. A number of difficulties arising from this
approach have recently been pointed out in the literature. We thoroughly
examine the effects of the assumptions made and numerical techniques employed
in SB simulations. This is done in order to clarify and gain better
understanding of those difficulties as well as of a number of additional
serious problems, raised here for the first time, and of their impact on the
results. Analytical derivations and estimates as well as comparative analysis
to methods used in the numerical study of turbulence are used. Numerical
experiments are performed to support some of our claims and conjectures. The
following problems, arising from the (virtually exclusive) use of the SB
simulations as a tool for the understanding and quantification of the nonlinear
MRI development in disks, are analyzed and discussed: (i) inconsistencies in
the application of the SB approximation itself; (ii) the limited spatial scale
of the SB; (iii) the lack of convergence of most ideal MHD simulations; (iv)
side-effects of the SB symmetry and the non-trivial nature of the linear MRI;
(v) physical artifacts arising on the too small box scale due to periodic
boundary conditions. The computational and theoretical challenge posed by the
MHD turbulence problem in accretion disks cannot be met by the SB
approximation, as it has been used to date. A new strategy to confront this
challenge is proposed, based on techniques widely used in numerical studies of
turbulent flows - developing (e.g., with the help of local numerical studies) a
sub-grid turbulence model and implementing it in global calculations.

[34] oai:arXiv.org:astro-ph/0702747 [pdf] - 89801

On the nonlinear saturation of the magnetorotational instability near
threshold in a thin-gap Taylor-Couette setup

Comments: Accepted with minor changes to Physical Review E

Submitted: **2007-02-28**, last modified: **2007-07-01**

We study the saturation near threshold of the axisymmetric magnetorotational
instability (MRI) of a viscous, resistive, incompressible fluid in a thin-gap
Taylor-Couette configuration. A vertical magnetic field, Keplerian shear and
no-slip, conducting radial boundary conditions are adopted. The weakly
non-linear theory leads to a real Ginzburg-Landau equation for the disturbance
amplitude, like in our previous idealized analysis. For small magnetic Prandtl
number (P \ll 1), the saturation amplitude scales as P^{2/3} while the
magnitude of angular momentum transport scales as P^{4/3}. The difference with
the previous scalings (~P^{1/2} and P respectively) is attributed to the
emergence of radial boundary layers. Away from those, steady-state non-linear
saturation is achieved through a modest reduction in the destabilizing shear.
These results will be useful to understand MRI laboratory experiments and
associated numerical simulations.

[35] oai:arXiv.org:astro-ph/0701217 [pdf] - 88346

An idealized model for dust-gas interaction in a rotating channel

Comments: Updated with more current references

Submitted: **2007-01-09**, last modified: **2007-01-11**

A 2D model representing the dynamical interaction of dust and gas in a
planetary channel is explored. The two components are treated as
interpenetrating fluids in which the gas is treated as a Boussinesq fluid while
the dust is treated as pressureless. The only coupling between both fluid
states is kinematic drag. The channel gas experiences a temperature gradient in
the spanwise direction and it is adverse the constant force of gravity. The
latter effects only the gas and not the dust component which is considered to
free float in the fluid. The channel is also considered on an f-plane so that
the background vorticity gradient can cause any emerging vortex structure to
drift like a Rossby wave. A linear theory analysis is explored and a nonlinear
amplitude theory is developed for disturbances of this arrangement. It is found
that the presence of the dust can help generate and shape emerging convection
patterns and dynamics in the gas so long as the state of the gas exceeds a
suitably defined Rayleigh number appropriate for describing drag effects. In
the linear stage the dust particles collect quickly onto sites in the gas where
the vorticity is minimal, i.e. where the disturbance vorticity is anticylonic
which is consistent with previous studies. The nonlinear theory shows that, in
turn, the local enhancement of dust concentration in the gas effects the vigor
of the emerging convective roll by modifying the effective local Rayleigh
number of the fluid. It is also found that without the f-plane approximation
built into the model the dynamics there is an algebraic runaway caused by
unrestrained growth in the dust concentration. The background vorticity
gradient forces the convective roll to drift like a Rossby wave and this causes
the dust concentration enhancements to not runaway.

[36] oai:arXiv.org:astro-ph/0608374 [pdf] - 84251

A weakly nonlinear analysis of the magnetorotational instability in a
model channel flow

Comments: Accepted to Physical Review Letters - Nov. 30, 2006. In final form

Submitted: **2006-08-17**, last modified: **2007-01-08**

We show by means of a perturbative weakly nonlinear analysis that the
axisymmetric magnetorotational instability (MRI) of a viscous, resistive,
incompressible rotating shear flow in a thin channel gives rise to a real
Ginzburg-Landau equation for the disturbance amplitude. For small magnetic
Prandtl number (${\cal P}_{\rm m}$), the saturation amplitude is $\propto
\sqrt{{\cal P}_{\rm m}}$ and the resulting momentum transport scales as ${\cal
R}^{-1}$, where $\cal R$ is the {\em hydrodynamic} Reynolds number. Simplifying
assumptions, such as linear shear base flow, mathematically expedient boundary
conditions and continuous spectrum of the vertical linear modes, are used to
facilitate this analysis. The asymptotic results are shown to comply with
numerical calculations using a spectral code. They suggest that the transport
due to the nonlinearly developed MRI may be very small in experimental setups
with ${\cal P}_{\rm m} \ll 1$.

[37] oai:arXiv.org:astro-ph/0506016 [pdf] - 73416

On the stratorotational instability in the quasi-hydrostatic
semi-geostrophic limit

Comments: 16 pages. Accepted in MNRAS (09/05). New sections added and abstract
changed

Submitted: **2005-06-01**, last modified: **2005-09-13**

The linear normal-mode stratorotational instability (SRI) is analytically
reexamined in the inviscid limit where the length scales of horizontal
disturbances are large compared their vertical and radial counterparts.
Boundary conditions different than channel walls are also considered. This
quasi-hydrostatic, semi-geostrophic (QHSG) approximation allows one to examine
the effect of a vertically varying Brunt-Vaisaila frequency, $N^2$. It is found
that the normal-mode instability persists when $N^2$ increases quadratically
with respect to the disc vertical coordinate. However we also find that the SRI
seems to exist in this inviscid QHSG extreme only for channel wall conditions:
when one or both of the reflecting walls are removed there is no instability in
the asymptotic limit explored here. It is also found that only exponential-type
SRI modes (as defined by Dubrulle et al. 2005) exist under these conditions.
These equations also admit non-normal mode behaviour. Fixed Lagrangian pressure
conditions on both radial boundaries predicts there to be no normal mode
behaviour in the QHSG limit. The mathematical relationship between the results
obtained here and that of the classic Eady (1949) problem for baroclinic
instability is drawn. We conjecture as to the mathematical/physical nature of
the SRI.
The general linear problem, analyzed without approximation in the context of
the Boussinesq equations, admits a potential vorticity-like quantity that is
advectively conserved by the shear. Its existence means that a continuous
spectrum \emph{is a generic feature of this system}. It also implies that in
places where the Brunt-Vaisaila frequency becomes dominant the linearized flow
may two-dimensionalize by advectively conserving its vertical vorticity.

[38] oai:arXiv.org:astro-ph/0502389 [pdf] - 71180

Hydrodynamical stability of thin accretion discs: transient growth of
global axisymmetric perturbations

Comments: 20 pages, 4 figues (last figure in jpg format). Submitted to MNRAS

Submitted: **2005-02-20**

The purpose of this paper is to explore how accretion discs manifest the
phenomenon of transient growth on a global scale. We investigate analytically
the time response of a thin accretion disc to particular axisymmetric
perturbations. To facilitate an analytical treatment we replace the energy
equation with a general polytropic assumption. The asymptotic expansion of
Klu\'zniak & Kita (2000), which extended the method of Regev (1983) to a full
steady polytropic disc (with $n=3/2$), is further developed and implemented for
both the steady (for any polytropic index) and time-dependent problems. The
spatial form and temporal behaviour of selected dynamical disturbances are
studied in detail. We identify the perturbation space which leads to transient
growth and provide analytical solutions which manifest this expected transient
growth behaviour. Three terms (physical causes) responsible for the appearance
of transient growth are identified. Two depend explicitly on the viscosity
while the third one is relevant also for inviscid discs. The main conclusion we
draw is that the phenomenon of transient growth exists in discs on a global
scale.

[39] oai:arXiv.org:astro-ph/0411393 [pdf] - 68992

On the nature of the hydrodynamic stability of accretion disks

Comments: Submitted to A+A Letters on November 15, 2004. 4 pages, no figures.
Accepted January 20, 2005. Section 3 rewritten and references updated

Submitted: **2004-11-15**, last modified: **2005-01-30**

The linear stability of accretion disks is revisited. The governing equations
are expanded asymptotically and solved to first order in the expansion
parameter $\epsilon$ defined by the ratio of the disk's vertical thickness to
its radial extent. Algebraically growing solutions are found for global
perturbations on the radial accretion flow of thin inviscid compressible
Keplerian disks. The algebraic temporal behavior is exhibited in the vertical
velocities and the thermodynamic variables and has the form $t\sin\Omega_0 t$
locally in the disk where $\Omega_0$ is the Keplerian rotation rate. The
physical implications and relations to the Solberg-Hoiland stability criteria
are discussed.

[40] oai:arXiv.org:astro-ph/0404020 [pdf] - 63928

Hydrodynamic stability of rotationally supported flows: Linear and
nonlinear 2D shearing box results

Comments: Submitted to A+A on April 1, 2004, Accepted August 3, 2004

Submitted: **2004-04-01**, last modified: **2004-08-31**

We present here both analytical and numerical results of hydrodynamic
stability investigations of rotationally supported circumstellar flows using
the shearing box formalism. Asymptotic scaling arguments justifying the
shearing box approximation are systematically derived, showing that there exist
two limits which we call small shearing box (SSB) and large shearing box (LSB).
The physical meaning of these two limits and their relationship to model
equations implemented by previous investigators are discussed briefly. Two
dimensional (2D) dynamics of the SSB are explored and shown to contain
transiently growing (TG) linear modes, whose nature is first discussed within
the context of linear theory. The fully nonlinear regime in 2D is investigated
numerically for very high Reynolds (Re) numbers. Solutions exhibiting sustained
dynamics are found and manifest episodic but recurrent TG behavior and these
are associated with the formation and long-term survival of coherent vortices.
The life-time of this spatio-temporal complexity depends on the Re number and
the strength and nature of the initial disturbance. We show results for a case
in which the dynamical activity persists for the entire duration of the
simulation (hundreds of box orbits). Some combinations of the Re number and the
initial perturbation spectrum (and strength) show the transiently growing
phenomenon to ultimately fade away or, in some instances, to be absent
altogether. Because the SSB approximation used here is equivalent to a 2D
incompressible flow, the dynamics can not depend on the Coriolis force.
Therefore, three dimensional (3D) simulations are needed in order to decide if
this force indeed suppresses nonlinear hydrodynamical instability in
rotationally supported disks in the shearing box.

[41] oai:arXiv.org:astro-ph/9806209 [pdf] - 101802

Stellar Oscillons

Comments: 12 Latex pages using aasms4.sty, 2 postscript figures, submitted to
the proceedings of the Florida Workshop in Nonlinear Astrophysics and Physics

Submitted: **1998-06-15**

We study the weakly nonlinear evolution of acoustic instability of a plane-
parallel polytrope with thermal dissipation in the form of Newton's law of
cooling. The most unstable horizontal wavenumbers form a band around zero and
this permits the development of a nonlinear pattern theory leading to a complex
Ginzburg-Landau equation (CGLE). Numerical solutions for a subcritical, quintic
CGLE produce vertically oscillating, localized structures that resemble the
oscillons observed in recent experiments of vibrated granular material.