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 Latter

Henrik 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.SRastro-ph.GAastro-ph.COphysics.flu-dynphysics.plasm-phastro-ph.HEastro-ph.IM

Catalog footprint

What is connected

20works
8topics
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

20 published item(s)

preprint2023arXiv

Thermal hysteresis and front propagation in dense planetary rings

Saturn's rings are composed of icy grains, most in the mm to m size ranges, undergoing several collisions per orbit. Their collective behaviour generates a remarkable array of structure over many orders of magnitude, much of it not well understood. On the other hand, the collisional properties and parameters of individual ring particles are poorly constrained; usually N-body simulations and kinetic theory employ hard-sphere models with a coefficient of restitution $ε$ that is constant or a decreasing function of impact speed. Due to plastic deformation of surface regolith, however, it is likely that $ε$ will be more complicated, at the very least a non-monotonic function. We undertake N-body simulations with the REBOUND code with non-monotonic $ε$ laws to approximate surfaces that are friable but not sticking. Our simulations reveal that such ring models can support two thermally stable steady states for the same (dynamical) optical depth: a cold and a warm state. If the ring breaks up into radial bands of one or the other state, we find that warmer states tend to migrate into the colder states via a coherent travelling front. We also find stationary `viscous' fronts, which connect states of different optical depth, but the same angular momentum flux. We discuss these preliminary results and speculate on their implications for structure formation in Saturn's B and C-rings, especially with respect to structures that appear in Cassini images but not in occultations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gravitoturbulent dynamo in global simulations of gaseous disks

The turbulence driven by gravitational instabilities (GIs) can amplify magnetic fields in massive gaseous disks. This GI dynamo may appear in young circumstellar disks, whose weak ionization challenges other amplification routes, as well as in active galactic nuclei. Although regarded as a large-scale dynamo, only local simulations have so far described its kinematic regime. We study the GI dynamo in global magnetohydrodynamic (MHD) models of accretion disks, focusing on its kinematic phase. We perform resistive MHD simulations with the Pluto code for different radiative cooling times and electrical resistivities. A weak magnetic field seeds the dynamo, and we adopt mean-field and heuristic models to capture its essence. We recover the same induction process leading to magnetic field amplification as previously identified in local simulations. The dynamo is, however, global in nature, connecting distant annuli of the disk via a large-scale dynamo mode of a fixed growth rate. This large-scale amplification can be described by a mean-field model that does not rely on conventional alpha-Omega effects. When varying the disk parameters we find an optimal resistivity that facilitates magnetic amplification, whose magnetic Reynolds number, Rm < 10, is substantially smaller than in local simulations. Unlike local simulations, we find an optimal cooling rate and the existence of global oscillating dynamo modes. The nonlinear saturation of the dynamo puts the disk in a strongly magnetized turbulent state on the margins of the effective range of GI. In our simulations, the accretion power eventually exceeds the threshold required by local thermal balance against cooling, leaving the long-term nonlinear outcome of the GI dynamo uncertain.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The saturation of the VSI in protoplanetary disks via parametric instability

The vertical shear instability (VSI) is a robust and potentially important phenomenon in irradiated protoplanetary disks (PPDs), yet the mechanism by which it saturates remains poorly understood. Global simulations suggest that the non-linear evolution of the VSI is dominated by radially propagating inertial wavetrains (called `body modes'), but these are known to be susceptible to a parametric instability. In this paper, we propose that the global VSI saturates via this secondary instability, which initiates a redistribution of energy from the large scales to smaller-scale inertial waves, and finally into a turbulent cascade. We present an analytic theory of the instability in a simple idealised model that captures the main physical and mathematical details of the problem. In addition, we conduct numerical simulations with the SNOOPY code to consolidate the theory. Once the parametric instability prevails, the VSI is likely far more disordered and incoherent than current global simulations suggest. We also argue that it is challenging to capture parametric instability in global simulations unless the radial resolution is very fine, possibly $\sim 300$ grid cells per scale height in radius.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Magneto-Thermal Instability In Galaxy Clusters I: Theory and Two-Dimensional Simulations

Determining the origin of turbulence in galaxy clusters, and quantifying its transport of heat, is an outstanding problem, with implications for our understanding of their thermodynamic history and structure. As the dilute plasma of the intracluster medium (ICM) is magnetized, heat and momentum travel preferentially along magnetic field lines. This anisotropy triggers a class of buoyancy instabilities that destabilize the ICM, and whose turbulent motions can augment or impede heat transport. We focus on the magneto-thermal instability (MTI), which may be active in the periphery of galaxy clusters. We aim to take a fresh look at the problem and construct a general theory that explains the MTI saturation mechanism and provides scalings and estimates for the turbulent kinetic energy, magnetic energy, and heat flux. We simulate MTI turbulence with a Boussinesq code, SNOOPY, which, in contrast to previous work, allows us to perform an extensive sampling of the parameter space. In two dimensions the saturation mechanism involves an inverse cascade carrying kinetic energy from the short MTI injection scales to larger scales, where it is arrested by the stable entropy stratification; at a characteristic "buoyancy scale", the energy is dumped into large-scale g-modes, which subsequently dissipate. Consequently, the entropy stratification sets an upper limit on the size and strength of turbulent eddies. Meanwhile, the MTI conveys a substantial fraction of heat, despite the tangled geometry of the magnetic field. In a companion paper, these results are extended to three-dimensional flows, and compared to real cluster observations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Magneto-Thermal Instability In Galaxy Clusters II: Three-Dimensional Simulations

In the intracluster medium (ICM) of galaxies, exchanges of heat across magnetic field lines are strongly suppressed. This anisotropic heat conduction, in the presence of a large-scale temperature gradient, destabilizes the outskirts of galaxy clusters via the magneto-thermal instability (MTI), and thus supplies a source of observed ICM turbulence. In this paper we continue our investigation of the MTI with 3D simulations using the Boussinesq code SNOOPY. We redress two issues intrinsic to our previous 2D study: an inverse energy cascade and the impossibility of dynamo action. Contrary to 2D simulations, we find inconsequential transport of energy across scales (most energy is dissipated at the same scale as its injection), and that turbulent eddies are vertically elongated at or below the thermal conduction length, but relatively isotropic on larger scales. Similar to 2D, however, the saturated turbulent energy levels and the integral scale follow clear power-laws that depend on the thermal diffusivity, temperature gradient, and buoyancy frequency. We also show that the MTI amplifies magnetic fields, through a fluctuation dynamo, to equipartition strengths provided that the integral scale of MTI turbulence is larger than the viscous dissipation scale. Finally, we show that our scaling laws are consistent with extant observations of ICM turbulence if the thermal conductivity is reduced by a factor of $\sim 10$ from its Spitzer value, and that on global cluster scales the stable stratification significantly reduces the vertical elongation of MTI motions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Electric heating and angular momentum transport in laminar models of protoplanetary disks

The vertical temperature structure of a protoplanetary disk bears on several processes relevant to planet formation, such as gas and dust grain chemistry, ice lines and convection. The temperature profile is controlled by irradiation from the central star and by any internal source of heat such as might arise from gas accretion. We investigate the heat and angular momentum transport generated by the resistive dissipation of magnetic fields in laminar disks. We use local one-dimensional simulations to obtain vertical temperature profiles for typical conditions in the inner disk (0.5 to 4 au). Using simple assumptions for the gas ionization and opacity, the heating and cooling rates are computed self-consistently in the framework of radiative non-ideal magnetohydrodynamics. We characterize steady solutions that are symmetric about the midplane and which may be associated with saturated Hall-shear unstable modes. We also examine the dissipation of electric currents driven by global accretion-ejection structures. In both cases we obtain significant heating for a sufficiently high opacity. Strong magnetic fields can induce an order-unity temperature increase in the disk midplane, a convectively unstable entropy profile, and a surface emissivity equivalent to a viscous heating of $α\sim 10^{-2}$. These results show how magnetic fields may drive efficient accretion and heating in weakly ionized disks where turbulence might be inefficient, at least for a range of radii and ages of the disk.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Global simulations of self-gravitating magnetized protoplanetary disks

In the early stages of a protoplanetary disk, when its mass is a significant fraction of its star's, turbulence generated by gravitational instability (GI) should feature significantly in the disk's evolution. At the same time, the disk may be sufficiently ionised for magnetic fields to play some role in the dynamics. Though usually neglected, the impact of magnetism on the GI may be critical, with consequences for several processes: the efficiency of accretion, spiral structure formation, fragmentation, and the dynamics of solids. In this paper, we report on global three-dimensional magnetohydrodynamical simulations of a self-gravitating protoplanetary disk using the meshless finite mass (MFM) Lagrangian technique. We confirm that GI spiral waves trigger a dynamo that amplifies an initial magnetic field to nearly thermal amplitudes (plasma beta < 10), an order of magnitude greater than that generated by the magneto-rotational instability alone. We also determine the dynamo's nonlinear back reaction on the gravitoturbulent flow: the saturated state is substantially hotter, with an associated larger Toomre parameter and weaker, more 'flocculent' spirals. But perhaps of greater import is the dynamo's boosting of accretion via a significant Maxwell stress; mass accretion is enhanced by factors of several relative to either pure GI or pure MRI. Our simulations use ideal MHD, an admittedly poor approximation in protoplanetary disks, and thus future studies should explore the full gamut of non-ideal MHD. In preparation for that, we exhibit a small number of Ohmic runs that reveal that the dynamo, if anything, is stronger in a non-ideal environment. This work confirms that magnetic fields are a potentially critical ingredient in gravitoturbulent young disks, possibly controlling their evolution, especially via their enhancement of (potentially episodic) accretion.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On dust-gas gravitational instabilities in protoplanetary discs

In protoplanetary disks the aerodynamical friction between particles and gas induces a variety of instabilities that facilitate planet formation. Of these we examine the so-called `secular gravitational instability' (SGI) in the two-fluid approximation, deriving analytical expressions for its stability criteria and growth rates. Concurrently, we present a physical explanation of the instability that shows how it manifests upon an intermediate range of lengthscales exhibiting geostrophic balance in the gas component. In contrast to a single-fluid treatment, the SGI is quenched within a critical disk radius, as large as 10 AU and 30 AU for cm and mm sized particles respectively, although establishing robust estimates is hampered by uncertainties in the parameters (especially the strength of turbulence) and deficiencies in the razor-thin disk model we employ. It is unlikely, however, that the SGI is relevant for well-coupled dust. We conclude by applying these results to the question of planetesimal formation and the provenance of large-scale dust rings.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the survival of zombie vortices in protoplanetary discs

Recently it has been proposed that the zombie vortex instability (ZVI) could precipitate hydrodynamical activity and angular momentum transport in unmagnetised regions of protoplanetary discs, also known as "dead zones". In this letter we scrutinise, with high resolution 3D spectral simulations, the onset and survival of this instability in the presence of viscous and thermal physics. First, we find that the ZVI is strongly dependent on the nature of the viscous operator. Although the ZVI is easily obtained with hyper-diffusion, it is difficult to sustain with physical (second order) diffusion operators up to Reynolds numbers as high as 10^7. This sensitivity is probably due to the ZVI's reliance on critical layers, whose characteristic lengthscale, structure, and dynamics are controlled by viscous diffusion. Second, we observe that the ZVI is sensitive to radiative processes, and indeed only operates when the Peclet number is greater than a critical value ~10^4, or when the cooling time is longer than ~10 Omega^{-1}. As a consequence, the ZVI struggles to appear at R>0.3 AU in standard 0.01 M_sun T Tauri disc models, though younger more massive disks provide a more hospitable environment. Together these results question the prevalence of the ZVI in protoplanetary discs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The magnetorotational instability in debris-disc gas

Debris discs are commonly swathed in gas which can be observed in UV, in fine structure lines in FIR, and in resolved maps of CO emission. Carbon and oxygen are overabundant in such gas, but it is severely depleted in hydrogen. As a consequence, its ionisation fraction is remarkably high, suggesting magnetohydrodynamic (MHD) processes may be important. In particular, the gas may be subject to the magnetorotational instability (MRI), and indeed recent modelling of $β$ Pictoris requires an anomalous viscosity to explain the gas's observed radial structure. In this paper we explore the possibility that the MRI is active in debris-disc gas and responsible for the observed mass transport. We find that non-ideal MHD and dust-gas interactions play a subdominant role, and that linear instability is viable at certain radii. However, owing to low gas densities, the outer parts of the disc could be stabilised by a weak ambient magnetic field, though it is difficult to constrain such a field. Even if the MRI is stabilised by too strong a field, a magnetocentrifugal wind may be launched in its place and this could lead to equivalent (non-turbulent) transport. Numerical simulations of the vertically stratified MRI in conditions appropriate to the debris disc gas should be able to determine the nature of the characteristic behaviour at different radii, and decide on the importance of the MRI (and MHD more generally) on the evolution of these discs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the convective overstability in protoplanetary discs

This paper explores the driving of low-level hydrodynamical activity in protoplanetary-disc dead zones. A small adverse radial entropy gradient, ordinarily stabilised by rotation, excites oscillatory convection (`convective overstability') when thermal diffusion, or cooling, is neither too strong nor too weak. I revisit the linear theory of the instability, discuss its prevalence in protoplanetary discs, and show that unstable modes are exact nonlinear solutions in the local Boussinesq limit. Overstable modes cannot grow indefinitely, however, as they are subject to a secondary parametric instability that limits their amplitudes to relatively low levels. If parasites set the saturation level of the ensuing turbulence then the convective overstability is probably too weak to drive significant angular momentum transport or to generate vortices. But I also discuss an alternative, and far more vigorous, saturation route that generates radial `layers' or `zonal flows' (witnessed also in semiconvection). Numerical simulations are required to determine which outcome is favoured in realistic discs, and consequently how important the instability is for disc dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The ballistic transport instability in Saturn's rings II: nonlinear wave dynamics

The ejecta discharged by impacting meteorites can redistribute a planetary ring's mass and angular momentum. This `ballistic transport' of ring properties instigates a linear instability that could generate the 100--1000-km undulations observed in Saturn's inner B-ring and in its C-ring. We present semi-analytic results demonstrating how the instability sustains steadily travelling nonlinear wavetrains. At low optical depths, the instability produces approximately sinusoidal waves of low amplitude, which we identify with those observed between radii 77,000 and 86,000 km in the C-ring. On the other hand, optical depths of 1 or more exhibit hysteresis, whereby the ring falls into multiple stable states: the homogeneous background equilibrium or large-amplitude wave states. Possibly the `flat zones' and `wave zones' between radii 93,000 and 98,000 km in the B-ring correspond to the stable homogeneous and wave states, respectively. In addition, we test the linear stability of the wavetrains and show that only a small subset are stable. In particular, stable solutions all possess wavelengths greater than the lengthscale of fastest linear growth. We supplement our calculations with a weakly nonlinear analysis that suggests the C-ring reproduces some of the dynamics of the complex Ginzburg--Landau equation. In the third paper in the series, these results will be tested and extended with numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The ballistic transport instability in Saturn's rings III: numerical simulations

Saturn's inner B-ring and its C-ring support wavetrains of contrasting amplitudes but with similar length scales, 100-1000 km. In addition, the inner B-ring is punctuated by two intriguing `flat' regions between radii 93,000 km and 98,000 km in which the waves die out, whereas the C-ring waves coexist with a forest of plateaus, narrow ringlets, and gaps. In both regions the waves are probably generated by a large-scale linear instability whose origin lies in the meteoritic bombardment of the rings: the ballistic transport instability. In this paper, the third in a series, we numerically simulate the long-term nonlinear evolution of this instability in a convenient local model. Our C-ring simulations confirm that the unstable system forms low-amplitude wavetrains possessing a preferred band of wavelengths. B-ring simulations, on the other hand, exhibit localised nonlinear wave `packets' separated by linearly stable flat zones. Wave packets travel slowly while spreading in time, a result that suggests the observed flat regions in Saturn's B-ring are shrinking. Finally, we present exploratory runs of the inner B-ring edge which reproduce earlier numerical results: ballistic transport can maintain the sharpness of a spreading edge while building a `ramp' structure at its base. Moreover, the ballistic transport instability can afflict the ramp region, but only in low-viscosity runs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Vortex cycles at the inner edges of dead zones in protoplanetary disks

In protoplanetary disks, the inner boundary between the turbulent and laminar regions is a promising site for planet formation because solids may become trapped at the interface itself or in vortices generated by the Rossby wave instability. The disk thermodynamics and the turbulent dynamics at that location are entwined because of the importance of turbulent dissipation on thermal ionization and, conversely, of thermal ionisation on the turbulence. However, most previous work has neglected this dynamical coupling and have thus missed a key element of the physics in this region. In this paper, we aim to determine how the the interplay between ionization and turbulence impacts on the formation and evolution of vortices at the interface between the active and the dead zones. Using the Godunov code RAMSES, we have performed a 3D magnetohydrodynamic global numerical simulation of a cylindrical model of an MRI--turbulent protoplanetary disk, including thermodynamical effects as well as a temperature-dependant resistivity. The comparison with an analogous 2D viscous simulation has been extensively used to help identify the relevant physical processes and the disk's long-term evolution. We find that a vortex formed at the interface, due to Rossby wave instability, migrates inward and penetrates the active zone where it is destroyed by turbulent motions. Subsequently, a new vortex emerges a few tens of orbits later at the interface, and the new vortex migrates inward too. The sequence repeats itself, resulting in cycles of vortex formation, migration, and disruption. This behavior is successfully reproduced using two different codes. In this paper, we characterize this vortex life cycle and discuss its implications for planet formation at the dead/active interface. Our simulations highlight the importance of thermodynamical processes for the vortex evolution at the dead zone inner edge.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Large-scale N-body simulations of the viscous overstability in Saturn's rings

We present results from large-scale particle simulations of the viscous overstability in Saturn's rings. The overstability generates a variety of structure on scales covering a few hundred metres to several kilometres, including axisymmetric wavetrains and larger-scale modulations. Such patterns have been observed in Saturn's rings by the Cassini spacecraft. Our simulations model the collisional evolution of particles in a co-rotating patch of the disk. These are the largest N-body simulations of the viscous overstability yet performed. The radial box size is five orders of magnitude larger than a typical particle radius, and so describes a 20-50 km radial portion of the rings. Its evolution is tracked for more than 10,000 orbits. In agreement with hydrodynamics, our N-body simulations reveal that the viscous overstability exhibits a rich set of dynamics characterised by nonlinear travelling waves with wavelengths of a few hundred meters. In addition, wave defects, such as sources and shocks, punctuate this bed of waves and break them up into large-scale divisions of radial width ~5 km. We find that the wavelength of the travelling waves is positively correlated with the mean optical depth. In order to assess the role of the numerical boundary conditions and also background ring structure, we include simulations of broad spreading rings and simulations with a gradient in the background surface density. Overall, our numerical results and approach provide a tool with which to interpret Cassini occultation observations of microstructure in Saturn's rings. We present an example of such a synthetic occultation observation and discuss what features to expect. We also make the entire source code freely available.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The ballistic transport instability in Saturn's rings I: formalism and linear theory

Planetary rings sustain a continual bombardment of hypervelocity meteoroids that erode the surfaces of ring particles on time scales of 10^5 - 10^7 years. The debris ejected from such impacts re-accretes on to the ring, though often at a slightly different orbital radius from the point of emission. This `ballistic transport' leads to a rearrangement of the disk's mass and angular momentum, and gives rise to a linear instability that generates structure on relatively large scales. It is likely that the 100-km undulations in Saturn's inner B-ring and the plateaus and 1000-km waves in Saturn's C-ring are connected to the nonlinear saturation of the instability. In this paper the physical problem is reformulated so as to apply to a local patch of disk (the shearing sheet). This new streamlined model helps facilitate our physical understanding of the instability, and also makes more tractable the analysis of its nonlinear dynamics. We concentrate on the linear theory in this paper, showing that the instability is restricted to a preferred range of intermediate wavenumbers and optical depths. We subsequently apply these general results to the inner B-ring and the C-ring and find that in both regions the ballistic transport instability should be near marginality, a fact that may have important consequences for its prevalence and nonlinear development. Owing to damping via self-gravity wakes, the instability should not be present in the A-ring. A following paper will explore the instability's nonlinear saturation and how it connects to the observed large-scale structure.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Global model of differential rotation in the Sun

The isorotation contours of the solar convective zone (SCZ) show three distinct morphologies, corresponding to two boundary layers (inner and outer), and the bulk of the interior. Previous work has shown that the thermal wind equation together with informal arguments on the nature of convection in a rotating fluid could be used to deduce the shape of the isorotation surfaces in the bulk of the SCZ with great fidelity, and that the tachocline contours could also be described by relatively simple phenomenology. In this paper, we show that the form of these surfaces can be understood more broadly as a mathematical consequence of the thermal wind equation and a narrow convective shell. The analysis does not yield the angular velocity function directly, an additional surface boundary condition is required. But much can already be deduced without constructing the entire rotation profile. The mathematics may be combined with dynamical arguments put forth in previous works to the mutual benefit of each. An important element of our approach is to regard the constant angular velocity surfaces as an independent coordinate variable for what is termed the "residual entropy," a quantity that plays a key role in the equation of thermal wind balance. The difference between the dynamics of the bulk of the SCZ and the tachocline is due to a different functional form of the residual entropy in each region. We develop a unified theory for the rotational behavior of both the SCZ and the tachocline, using the solutions for the characteristics of the thermal wind equation. These characteristics are identical to the isorotation contours in the bulk of the SCZ, but the two deviate in the tachocline. The outer layer may be treated, at least descriptively, by similar mathematical techniques, but this region probably does not obey thermal wind balance.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Spherical single-roll dynamos at large magnetic Reynolds numbers

This paper concerns kinematic helical dynamos in a spherical fluid body surrounded by an insulator. In particular, we examine their behaviour in the regime of large magnetic Reynolds number $\Rm$, for which dynamo action is usually concentrated upon a simple resonant stream-surface. The dynamo eigensolutions are computed numerically for two representative single-roll flows using a compact spherical harmonic decomposition and fourth-order finite-differences in radius. These solutions are then compared with the growth rates and eigenfunctions of the Gilbert and Ponty (2000) large $\Rm$ asymptotic theory. We find good agreement between the growth rates when $\Rm>10^4$, and between the eigenfunctions when $\Rm>10^5$.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A comparison of local simulations and reduced models of MRI-induced turbulence

We run mean-field shearing-box numerical simulations with a temperature-dependent resistivity and compare them to a reduced dynamical model. Our simulations reveal the co-existence of two quasi-steady states, a `quiet' state and an `active' turbulent state, confirming the predictions of the reduced model. The initial conditions determine on which state the simulation ultimately settles. The active state is strongly influenced by the geometry of the computational box and the thermal properties of the gas. Cubic domains support permanent channel flows, bar-shaped domains exhibit eruptive behaviour, and horizontal slabs give rise to infrequent channels. Meanwhile, longer cooling time-scales lead to higher saturation amplitudes.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Inertial waves near corotation in 3D hydrodynamical disks

This paper concerns the interaction between non-axisymmetric inertial waves and their corotation resonances in a hydrodynamical disk. Inertial waves are of interest because they can localise in resonant cavities circumscribed by Lindblad radii, and as a consequence exhibit discrete oscillation frequencies that may be observed. It is often hypothesised that these trapped eigenmodes are affiliated with the poorly understood QPO phenomenon. We demonstrate that a large class of non-axisymmetric 3D inertial waves cannot manifest as trapped normal modes. This class includes any inertial wave whose resonant cavity contains a corotation singularity. Instead, these `singular' modes constitute a continuous spectrum and, as an ensemble, are convected with the flow, giving rise to shearing waves. Lastly, we present a simple demonstration of how the corotation singularity stabilizes three-dimensional perturbations in a slender torus.

0 citations0 reviewsNo code linkedOpen access