BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Henrik N. Latter

Henrik N. Latter appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.HEastro-ph.SRphysics.flu-dynastro-ph.COastro-ph.GA

Catalog footprint

What is connected

24works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2024arXiv

A physical picture for the acoustic resonant drag instability

Mixtures of gas and dust are pervasive in the universe, from AGN and molecular clouds to proto-planetary discs. When the two species drift relative to each other, a large class of instabilities can arise, called resonant drag instabilities (RDIs). The most famous RDI is the streaming instability, which plays an important role in planet formation. On the other hand, acoustic RDIs, the simplest kind, feature in the winds of cool stars, AGN, or starburst regions. Unfortunately, owing to the complicated dynamics of two coupled fluids (gas and dust), the underlying physics of most RDIs is mysterious. In this paper, we develop a clear physical picture of how the acoustic RDI arises and support this explanation with transparent mathematics. We find that the acoustic RDI is built on two coupled mechanisms. In the first, the converging flows of a sound wave concentrate dust. In the second, a drifting dust clump excites sound waves. These processes feed into each other at resonance, thereby closing an unstable feedback loop. This physical picture helps decide where and when RDIs are most likely to happen, and what can suppress them. Additionally, we find that the acoustic RDI remains strong far from resonance. This second result suggests that one can simulate RDIs without having to fine-tune the dimensions of the numerical domain.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The vertical shear instability in poorly ionised, magnetized protoplanetary discs

Protoplanetary discs should exhibit a weak vertical variation in their rotation profiles. Typically this `vertical shear' issues from a baroclinic effect driven by the central star's radiation field, but it might also arise during the launching of a magnetocentrifugal wind. As a consequence, protoplanetary discs are subject to a hydrodynamical instability, the `vertical shear instability' (VSI), whose breakdown into turbulence could transport a moderate amount of angular momentum and facilitate, or interfere with, the process of planet formation. Magnetic fields may suppress the VSI, however, either directly via magnetic tension or indirectly through magnetorotational turbulence. On the other hand, protoplanetary discs exhibit notoriously low ionisation fractions, and non-ideal effects, if sufficiently dominant, may come to the VSI's rescue. In this paper we develop a local linear theory that explores how non-ideal MHD influences the VSI, while also launching additional diffusive shear instabilities. We derive a set of analytical criteria that establish when the VSI prevails, and then show how it can be applied to a realistic global model of a protoplanetary disc. Our calculations suggest that within ~10au the VSI should have little trouble emerging in the main body of the disk, but beyond that, and in the upper regions of the disc, its onset depends sensitively on the size of the preponderant dust grains.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Wavelike nature of the vertical shear instability in global protoplanetary disks

The vertical shear instability (VSI) is a robust phenomenon in irradiated protoplanetary disks (PPDs). The majority of previous numerical simulations have focused on the turbulent properties of its saturated state. However, the saturation of the VSI manifests as large-scale coherent radially travelling inertial waves. In this paper, we study inertial-wave-disk interactions and their impact on VSI saturation. Inertial-wave linear theory is developed and applied to a representative global 2D simulation using the Athena++ code. It is found that the VSI saturates by separating the disk into several radial wave zones roughly demarcated by corotation resonances (turning points); this structure also manifests in modest radial variations in the vertical turbulence strength. Future numerical work should employ large radial domains to accommodate this radial structure of the VSI, while concurrently adopting sufficiently fine resolutions to resolve the parametric instability that attacks the saturated VSI inertial waves.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

HFQPOs and discoseismic mode excitation in eccentric, relativistic discs. I. Hydrodynamic simulations

High-frequency quasi-periodic oscillations (HFQPOs) observed in the emission of black-hole X-ray binary systems promise insight into strongly curved spacetime. `Discoseismic' oscillations with frequencies set by the intrinsic properties of the central black hole, in particular `trapped inertial waves' (r-modes), offer an attractive explanation for HFQPOs. To produce an observable signature, however, such oscillations must be excited to sufficiently large amplitudes. Turbulence driven by the magnetorotational instability (MRI) fails to provide the necessary amplification, but r-modes may still be excited via interaction with accretion disc warps or eccentricities. We present 3D global hydrodynamic simulations of relativistic accretion discs, which demonstrate for the first time the excitation of trapped inertial waves by an imposed eccentricity in the flow. While the r-modes' saturated state depends on the vertical boundary conditions used in our unstratified, cylindrical framework, their excitation is unambiguous in all runs with eccentricity >0.005 near the ISCO. These simulations provide a proof of concept, demonstrating the robustness of trapped inertial wave excitation in a non-magnetized context. In a companion paper, we explore the competition between this excitation, and damping by magnetohydrodynamic turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

HFQPOs and discoseismic mode excitation in eccentric, relativistic discs. II. Magnetohydrodynamic simulations

Trapped inertial oscillations (r-modes) provide a promising explanation for high-frequency quasi-periodic oscillations (HFQPOs) observed in the emission from black hole X-ray binary systems. An eccentricity (or warp) can excite r-modes to large amplitudes, but concurrently the oscillations are likely damped by magnetohydrodynamic (MHD) turbulence driven by the magnetorotational instability (MRI). We force eccentricity in global, unstratified, zero-net flux MHD simulations of relativistic accretion discs, and find that a sufficiently strong disc distortion generates trapped inertial waves despite this damping. In our simulations, eccentricities above ~ 0.03 in the inner disc excite trapped waves. In addition to the competition between r-mode damping and driving, we observe that larger amplitude eccentric structures modify and in some cases suppress MRI turbulence. Given the variety of distortions (warps as well as eccentricities) capable of amplifying r-modes, the robustness of trapped inertial wave excitation in the face of MRI turbulence in our simulations provides support for a discoseismic explanation for HFQPOs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Local and global aspects of the linear MRI in accretion disks

We revisit the linear MRI in a cylindrical model of an accretion disk and uncover a number of attractive results overlooked in previous treatments. In particular, we elucidate the connection between local axisymmetric modes and global modes, and show that a local channel flow corresponds to the evanescent part of a global mode. In addition, we find that the global problem reproduces the local dispersion relation without approximation, a result that helps explain the success the local analysis enjoys in predicting global growth rates. MRI channel flows are nonlinear solutions to the governing equations in the local shearing box. However, only a small subset of MRI modes share the same property in global disk models, providing further evidence that the prominence of channels in local boxes is artificial. Finally, we verify our results via direct numerical simulations with the Godunov code RAMSES.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the vertical-shear instability in astrophysical discs

We explore the linear stability of astrophysical discs exhibiting vertical shear, which arises when there is a radial variation in the temperature or entropy. Such discs are subject to a "vertical-shear instability", which recent nonlinear simulations have shown to drive hydrodynamic activity in the MRI-stable regions of protoplanetary discs. We first revisit locally isothermal discs using the quasi-global reduced model derived by Nelson et al. (2013). This analysis is then extended to global axisymmetric perturbations in a cylindrical domain. We also derive and study a reduced model describing discs with power law radial entropy profiles ("locally polytropic discs"), which are somewhat more realistic in that they possess physical (as opposed to numerical) surfaces. In all cases the fastest growing modes have very short wavelengths and are localised at the disc surfaces (if present), where the vertical shear is maximal. An additional class of modestly growing vertically global body modes is excited, corresponding to destabilised classical inertial waves ("r-modes"). We discuss the properties of both types of modes, and stress that those that grow fastest occur on the shortest available length scales (determined either by the numerical grid or the physical viscous length). This ill-posedness makes simulations of the instability difficult to interpret. We end with some brief speculation on the nonlinear saturation and resulting angular momentum transport.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The stress-pressure relationship in simulations of MRI-induced turbulence

We determine how MRI-turbulent stresses depend on gas pressure via a suite of unstratified shearing box simulations. Earlier numerical work reported only a very weak dependence at best, results that call into question the canonical alpha-disk model and the thermal stability results that follow from it. Our simulations, in contrast, exhibit a stronger relationship, and show that previous work was box-size limited: turbulent `eddies' were artificially restricted by the numerical domain rather than by the scale height. Zero-net-flux runs without physical diffusion coefficients yield a stress proportional to $P^{0.5}$, where P is pressure. The stresses are also proportional to the grid length and hence remain numerically unconverged. The same runs with physical diffusivities, however, give a result closer to an alpha-disk: the stress is proportional to $P^{0.9}$. Net-flux simulations without explicit diffusion exhibit stresses proportional to $P^{0.5}$, but stronger imposed fields weaken this correlation. In summary, compressibility is important for the saturation of the MRI, but the exact stress-pressure relationship is difficult to ascertain in local simulations because of numerical convergence issues and the influence of any imposed flux. As a consequence, the interpretation of thermal stability behaviour in local simulations is a problematic enterprise.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Thermodynamics of the dead-zone inner edge in protoplanetary disks

In protoplanetary disks, the inner boundary between the turbulent and laminar regions could be a promising site for planet formation, thanks to the trapping of solids at the boundary itself or in vortices generated by the Rossby wave instability. At the interface, the disk thermodynamics and the turbulent dynamics are entwined because of the importance of turbulent dissipation and thermal ionization. Numerical models of the boundary, however, have neglected the thermodynamics, and thus miss a part of the physics. The aim of this paper is to numerically investigate the interplay between thermodynamics and dynamics in the inner regions of protoplanetary disks by properly accounting for turbulent heating and the dependence of the resistivity on the local temperature. Using the Godunov code RAMSES, we performed a series of 3D global numerical simulations of protoplanetary disks in the cylindrical limit, including turbulent heating and a simple prescription for radiative cooling. We find that waves excited by the turbulence significantly heat the dead zone, and we subsequently provide a simple theoretical framework for estimating the wave heating and consequent temperature profile. In addition, our simulations reveal that the dead-zone inner edge can propagate outward into the dead zone, before staling at a critical radius that can be estimated from a mean-field model. The engine driving the propagation is in fact density wave heating close to the interface. A pressure maximum appears at the interface in all simulations, and we note the emergence of the Rossby wave instability in simulations with extended azimuth. Our simulations illustrate the complex interplay between thermodynamics and turbulent dynamics in the inner regions of protoplanetary disks. They also reveal how important activity at the dead-zone interface can be for the dead-zone thermodynamic structure.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hydrodynamic instability in warped astrophysical discs

Warped astrophysical discs are usually treated as laminar viscous flows, which have anomalous properties when the disc is nearly Keplerian and the viscosity is small: fast horizontal shearing motions and large torques are generated, which cause the warp to evolve rapidly, in some cases at a rate that is inversely proportional to the viscosity. However, these flows are often subject to a linear hydrodynamic instability, which may produce small-scale turbulence and modify the large-scale dynamics of the disc. We use a warped shearing sheet to compute the oscillatory laminar flows in a warped disc and to analyse their linear stability by the Floquet method. We find widespread hydrodynamic instability deriving from the parametric resonance of inertial waves. Even very small, unobservable warps in nearly Keplerian discs of low viscosity can be expected to generate hydrodynamic turbulence, or at least wave activity, by this mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Local and global dynamics of warped astrophysical discs

Astrophysical discs are warped whenever a misalignment is present in the system, or when a flat disc is made unstable by external forces. The evolution of the shape and mass distribution of a warped disc is driven not only by external influences but also by an internal torque, which transports angular momentum through the disc. This torque depends on internal flows driven by the oscillating pressure gradient associated with the warp, and on physical processes operating on smaller scales, which may include instability and turbulence. We introduce a local model for the detailed study of warped discs. Starting from the shearing sheet of Goldreich & Lynden-Bell, we impose the oscillating geometry of the orbital plane by means of a coordinate transformation. This warped shearing sheet (or box) is suitable for analytical and computational treatments of fluid dynamics, magnetohydrodynamics, etc., and it can be used to compute the internal torque that drives the large-scale evolution of the disc. The simplest hydrodynamic states in the local model are horizontally uniform laminar flows that oscillate at the orbital frequency. These correspond to the nonlinear solutions for warped discs found in previous work by Ogilvie, and we present an alternative derivation and generalization of that theory. In a companion paper we show that these laminar flows are often linearly unstable, especially if the disc is nearly Keplerian and of low viscosity. The local model can be used in future work to determine the nonlinear outcome of the hydrodynamic instability of warped discs, and its interaction with others such as the magnetorotational instability.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Local outflows from turbulent accretion disks

The aim of this paper is to investigate the properties of accretion disks threaded by a weak vertical magnetic field, with a particular focus on the interplay between MHD turbulence driven by the magnetorotational instability (MRI) and outflows that might be launched from the disk. For that purpose, we use a set of numerical simulations performed with the MHD code RAMSES in the framework of the shearing box model. We concentrate on the case of a rather weak vertical magnetic field such that the initial ratio beta0 of the thermal and magnetic pressures in the disk midplane equals 10^4. As reported recently, we find that MHD turbulence drives an efficient outflow out of the computational box. We demonstrate a strong sensitivity of that result to the box size: enlargements in the radial and vertical directions lead to a reduction of up to an order of magnitude in the mass-loss rate. Such a dependence prevents any realistic estimates of disk mass-loss rates being derived using shearing-box simulations. We find however that the flow morphology is robust and independent of the numerical details of the simulations. Its properties display some features and approximate invariants that are reminiscent of the Blandford & Payne launching mechanism, but differences exist. For the magnetic field strength considered in this paper, we also find that angular momentum transport is most likely dominated by MHD turbulence, the saturation of which scales with the magnetic Prandtl number, the ratio of viscosity and resistivity, in a way that is in good agreement with expectations based on unstratified simulations. This paper thus demonstrates for the first time that accretion disks can simultaneously exhibit MRI-driven MHD turbulence along with magneto-centrifugally accelerated outflows.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Hysteresis and thermal limit cycles in MRI simulations of accretion discs

The recurrent outbursts that characterise low-mass binary systems reflect thermal state changes in their associated accretion discs. The observed outbursts are connected to the strong variation in disc opacity as hydrogen ionises near 5000 K. This physics leads to accretion disc models that exhibit bistability and thermal limit cycles, whereby the disc jumps between a family of cool and low accreting states and a family of hot and efficiently accreting states. Previous models have parametrised the disc turbulence via an alpha (or `eddy') viscosity. In this paper we treat the turbulence more realistically via a suite of numerical simulations of the magnetorotational instability (MRI) in local geometry. Radiative cooling is included via a simple but physically motivated prescription. We show the existence of bistable equilibria and thus the prospect of thermal limit cycles, and in so doing demonstrate that MRI-induced turbulence is compatible with the classical theory. Our simulations also show that the turbulent stress and pressure perturbations are only weakly dependent on each other on orbital times; as a consequence, thermal instability connected to variations in turbulent heating (as opposed to radiative cooling) is unlikely to operate, in agreement with previous numerical results. Our work presents a first step towards unifying simulations of full MHD turbulence with the correct thermal and radiative physics of the outbursting discs associated with dwarf novae, low-mass X-ray binaries, and possibly young stellar objects.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The dynamics of inner dead-zone boundaries in protoplanetary disks

In protoplanetary disks, the inner radial boundary between the MRI turbulent (`active') and MRI quiescent (`dead') zones plays an important role in models of the disk evolution and in some planet formation scenarios. In reality, this boundary is not well-defined: thermal heating from the star in a passive disk yields a transition radius close to the star (<0.1 au), whereas if the disk is already MRI active, it can self-consistently maintain the requisite temperatures out to a transition radius of roughly 1 au. Moreover, the interface may not be static; it may be highly fluctuating or else unstable. In this paper, we study a reduced model of the dynamics of the active/dead zone interface that mimics several important aspects of a real disk system. We find that MRI-transition fronts propagate inward (a `dead front' suppressing the MRI) if they are initially at the larger transition radius, or propagate outward (an `active front' igniting the MRI) if starting from the smaller transition radius. In both cases, the front stalls at a well-defined intermediate radius, where it remains in a quasi-static equilibrium. We propose that it is this new, intermediate stalling radius that functions as the true boundary between the active and dead zones in protoplanetary disks. These dynamics are likely implicated in observations of variable accretion, such as FU Ori outbursts, as well as in those planet formation theories that require the accumulation of solid material at the dead/active interface.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The gravitational instability of a stream of co-orbital particles

We describe the dynamics of a stream of equally spaced macroscopic particles in orbit around a central body (e.g. a planet or star). A co-orbital configuration of small bodies may be subject to gravitational instability, which takes the system to a spreading, disordered and collisional state. We detail the linear instability's mathematical and physical features using the shearing sheet model and subsequently track its nonlinear evolution with local N-body simulations. This model provides a convenient tool with which to understand the gravitational and collisional dynamics of narrow belts, such as Saturn's F-ring and the streams of material wrenched from tidally disrupted bodies. In particular, we study the tendency of these systems to form long-lived particle aggregates. Finally, we uncover an unexpected connection between the linear dynamics of the gravitational instability and the magnetorotational instability.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The HBI in a quasi-global model of the intracluster medium

In this paper we investigate how convective instabilities influence heat conduction in the intracluster medium (ICM) of cool-core galaxy clusters. The ICM is a high-beta, weakly collisional plasma in which the transport of momentum and heat is aligned with the magnetic field. The anisotropy of heat conduction, in particular, gives rise to instabilities that can access energy stored in a temperature gradient of either sign. We focus on the heat-flux buoyancy-driven instability (HBI), which feeds on the outwardly increasing temperature profile of cluster cool cores. Our aim is to elucidate how the global structure of a cluster impacts on the growth and morphology of the linear HBI modes when in the presence of Braginskii viscosity, and ultimately on the ability of the HBI to thermally insulate cores. We employ an idealised quasi-global model, the plane-parallel atmosphere, which captures the essential physics -- e.g. the global radial profile of the cluster -- while letting the problem remain analytically tractable. Our main result is that the dominant HBI modes are localised to the the innermost (~<20%) regions of cool cores. It is then probable that, in the nonlinear regime, appreciable field-line insulation will be similarly localised. Thus, while radio-mode feedback appears necessary in the central few tens of kpc, heat conduction may be capable of offsetting radiative losses throughout most of a cool core over a significant fraction of the Hubble time. Finally, our linear solutions provide a convenient numerical test for the nonlinear codes that tackle the saturation of such convective instabilities in the presence of anisotropic transport.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Tidal disruption of satellites and formation of narrow rings

In this paper we investigate the formation of narrow planetary rings such as those found around Uranus and Saturn through the tidal disruption of a weak, gravitationally bound satellite that migrates within its Roche limit. Using $N$-body simulations, we study the behaviour of rubble piles placed on circular orbits at different distances from a central planet. We consider both homogeneous satellites and differentiated bodies containing a denser core. We show that the Roche limit for a rubble pile is closer to the planet than for a fluid body of the same mean density. The Roche limit for a differentiated body is also closer to the planet than for a homogeneous satellite of the same mean density. Within its Roche limit, a homogeneous satellite totally disrupts and forms a narrow ring. The initial stages of the disruption are similar to the evolution of a viscous fluid ellipsoid, which can be computed semi-analytically. On the other hand, when a differentiated satellite is just within the Roche limit only the mantle is disrupted. This process is similar to Roche-lobe overflow in interacting binary stars and produces two narrow rings on either side of a remnant satellite. We argue that the Uranian rings, and possibly their shepherd satellites, could have been formed through the tidal disruption of a number of protomoons that were formed inside the corotation radius of Uranus and migrated slowly inwards as a result of tidal interaction with the planet.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On linear dust-gas streaming instabilities in protoplanetary discs

We revisit, via a very simplified set of equations, a linear streaming instability (technically an overstability), which is present in, and potentially important for, dusty protoplanetary disks (Youdin & Goodman 2005). The goal is a better understanding of the physical origin of such instabilities, which are notoriously subtle. Rotational dynamics seem to be essential to this type of instability, which cannot be captured by one-dimensional Cartesian models. Dust `pileups' in moving pressure maxima are an important triggering mechanism of the instability, and drag feedback of dust upon the gas allows these maxima to be strengthened. Coriolis forces and the background drift counteract the effects of the pressure force.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Hydrodynamical simulations of viscous overstability in Saturn's rings

We perform axisymmetric hydrodynamical simulations that describe the nonlinear outcome of the viscous overstability in dense planetary rings. These simulations are particularly relevant for Cassini observations of fine-scale structure in Saturn's A and B-ring, which take the form of periodic microstructure on the 0.1 km scale, and irregular larger-scale structure on 1-10 km. Nonlinear wavetrains dominate all the simulations, and we associate them with the observed periodic microstructure. The waves can undergo small chaotic fluctuations in their phase and amplitude, and may be punctuated by more formidable `wave defects' distributed on longer scales. It is unclear, however, whether the defects are connected to the irregular larger-scale variations observed by Cassini. The long-term behaviour of the simulations is dominated by the imposed boundary conditions, and more generally by the limitations of the local model we use: the shearing box. When periodic boundary conditions are imposed, the system eventually settles on a uniform travelling wave of a predictable wavelength, while reflecting boundaries, and boundaries with buffer zones, maintain a disordered state. The simulations omit self-gravity, though we examine its influence in future work.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

MRI channel flows in vertically-stratified models of accretion disks

Simulations of the magnetorotational instability (MRI) in 'unstratified' shearing boxes exhibit powerful coherent flows, whereby the fluid vertically splits into countermoving planar jets or `channels'. Channel flows correspond to certain axisymmetric linear MRI modes, and their preponderance follows from the remarkable fact that they are approximate nonlinear solutions of the MHD equations in the limit of weak magnetic fields. We show in this paper, analytically and with one-dimensional numerical simulations, that this property is also shared by certain axisymmetric MRI modes in vertically-stratified shearing boxes. These channel flows rapidly capture significant amounts of magnetic and kinetic energy, and thus are vulnerable to secondary shear instabilities. We examine these parasites in the vertically stratified context, and estimate the maximum amplitudes that channels attain before they are destroyed. These estimates suggest that a dominant channel flow will usually drive the disk's magnetic field to thermal strengths. The prominence of these flows and their destruction place enormous demands on simulations, but channels in their initial stages also offer a useful check on numerical codes. These benchmarks are especially valuable given the increasing interest in the saturation of the stratified MRI. Lastly we speculate on the potential connection between 'run-away' channel flows and outburst behaviour in protostellar and dwarf nova disks.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Resistive double-diffusive instability in the dead-zones of protostellar disks

We outline a novel linear instability that may arise in the dead-zones of protostellar disks, and possibly the fluid interiors of planets and protoplanets. In essence it is an axisymmetric buoyancy instability, but one that would not be present in a purely hydrodynamical gas. The necessary ingredients for growth include a negative radial entropy gradient (of any magnitude), weak magnetic fields, and efficient resistive diffusion (in comparison with thermal diffusion). The character of the instability is local, axisymmetric, and double-diffusive, and it attacks lengths much shorter than the resistive scale. Like the axisymmetric convective instability, it draws its energy from the negative radial entropy gradient; but by utilising the diffusing magnetic field, it can negate the stabilising influence of rotation. Its nonlinear saturated state, while not transporting appreciable angular momentum, could drive radial and vertical mixing, which may influence the temperature structure of the disk, dust dynamics and, potentially, planet formation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The effect of the tachocline on differential rotation in the Sun

In this paper, we present a model for the effects of the tachocline on the differential rotation in the solar convection zone. The mathematical technique relies on the assumption that entropy is nearly constant ("well-mixed") in isorotation surfaces both outside and within the tachocline. The resulting solutions exhibit nontrivial features that strikingly resemble the true tachocline isorotation contours in unexpected detail. This strengthens the mathematical premises of the theory. The observed rotation pattern in the tachocline shows strong quadrupolar structure, an important feature that is explicitly used in constructing our solutions. The tachocline is treated locally as an interior boundary layer of small but finite thickness, and an explicit global solution is then constructed. A dynamical link can thus be established between the internal jump in the angular velocity at the tachocline and the spread of angular velocities observed near the solar surface. In general, our results suggest that the bulk of the solar convection zone is in thermal wind balance, and that simple quadrupolar stresses, local in radius, mediate the tachocline transition from differential rotation to uniform rotation in the radiative interior.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

MRI channel flows and their parasites

Local simulations of the magnetorotational instability (MRI) in accretion disks can exhibit recurrent coherent structures called channel flows. The formation and destruction of these structures may play a role in the development and saturation of MRI-induced turbulence, and consequently help us understand the time-dependent accretion behaviour of certain astrophysical objects. Previous investigations have revealed that channel solutions are attacked by various parasitic modes, foremost of which is an analogue of the Kelvin-Helmholtz instability. We revisit these instabilities and show how they relate to the classical instabilities of plasma physics, the kink and pinch modes. However, we argue that in most cases channels emerge from developed turbulence and are eventually destroyed by turbulent mixing, not by the parasites. The exceptions are the clean isolated channels which appear in systems near criticality or which emerge from low amplitude initial conditions. These structures inevitably achieve large amplitudes and are only then destroyed, giving rise to eruptive behaviour.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The viscous overstability, nonlinear wavetrains, and finescale structure in dense planetary rings

This paper addresses the fine-scale axisymmetric structure exhibited in Saturn's A and B-rings. We aim to explain both the periodic microstructure on 150-220m, revealed by the Cassini UVIS and RSS instruments, and the irregular variations in brightness on 1-10km, reported by the Cassini ISS. We propose that the former structures correspond to the peaks and troughs of the nonlinear wavetrains that form naturally in a viscously overstable disk. The latter variations on longer scales may correspond to modulations and defects in the wavetrains' amplitudes and wavelength. We explore these ideas using a simple hydrodynamical model which captures the correct qualitative behaviour of a disk of inelastically colliding particles, while also permitting us to make progress with analytic and semi-analytic techniques. Specifically, we calculate a family of travelling nonlinear density waves and determine their stability properties. Detailed numerical simulations that confirm our basic results will appear in a following paper.

0 citations0 reviewsNo code linkedOpen access