Researcher profile

Steven A. Balbus

Steven A. Balbus contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

18 published item(s)

preprint2021arXiv

A Poynting theorem formulation for the gravitational wave stress pseudo tensor

A very simple and physical derivation of the conservation equation for the propagation of gravitational radiation is presented. The formulation is exact. The result takes the readily recognisable and intuitive form of a Poynting-style equation, in which the outward propagation of stress-energy is directly related to the volumetric equivalent of a radiation reaction force acting back upon the sources, including the purely gravitational contribution to the sources. Upon averaging, the emergent pseudo tensor for the gravitational radiation is in exact agreement with that found by much more labour-intensive methods.

preprint2020arXiv

The spectral evolution of disc dominated tidal disruption events

We perform a detailed numerical and analytical study of the properties of observed light curves from relativistic thin discs, focussing on observational bands most appropriate for comparison with tidal disruption events (TDEs). We make use of asymptotic expansion techniques applied to the spectral emission integral, using time dependent disc temperature profiles appropriate for solutions of the relativistic thin disc equation. Rather than a power law associated with bolometric disc emission $L \sim t^{-n}$, the observed X-ray flux from disc-dominated TDEs will typically have the form of a power law multiplied by an exponential (see eq. 91). While precise details are somewhat dependent on the nature of the ISCO stress and disc-observer orientational angle, the general form of the time-dependent flux is robust and insensitive to the exact disc temperature profile. We present numerical fits to the UV and X-ray light curves of ASASSN-14li, a particularly well observed TDE. This modelling incorporates strong gravity optics. The full 900 days of ASASSN-14li X-ray observations are very well fit by a simple relativistic disc model, significantly improving upon previous work. The same underlying model also fits the final 1000 days of ASASSN-14li observations in three different UV bandpasses. Finally, we demonstrate that the analytic formulae reproduce the properties of full numerical modelling at both UV and X-ray wavelengths with great fidelity.

preprint2016arXiv

Surprises in astrophysical gasdynamics

Much of astrophysics consists of the study of ionised gas under the influence of gravitational and magnetic fields. Thus, it is not possible to understand the astrophysical universe without a detailed knowledge of the dynamics of magnetised fluids. Fluid dynamics is, however, a notoriously tricky subject, in which it is all too easy for one's a priori intuition to go astray. In this review, we seek to guide the reader through a series of illuminating yet deceptive problems, all with an enlightening twist. We cover a broad range of topics including the instabilities acting in accretion discs, the hydrodynamics governing the convective zone of the Sun, the magnetic shielding of a cooling galaxy cluster, and the behaviour of thermal instabilities and evaporating clouds. The aim of this review is to surprise and intrigue even veteran astrophysical theorists with an idiosynchratic choice of problems and counterintuitive results. At the same time, we endeavour to bring forth the fundamental ideas, to set out important assumptions, and to describe carefully whatever novel techniques may be appropriate to the problem at hand. By beginning at the beginning, and analysing a wide variety of astrophysical settings, we seek not only to make this review suitable for fluid dynamic veterans, but to engage novice recruits as well with what we hope will be an unusual and instructive introduction to the subject.

preprint2016arXiv

The Goldreich-Schubert-Fricke instability in stellar radiative zones

The Goldreich-Schubert-Fricke (GSF) instability is a rotational instability that is thought to contribute to the transfer of angular momentum in differentially rotating stars. It has been included in recent codes of stellar evolution in a diffusion-like approximation, under the assumption that the kinematic viscosity $ν$ is unimportant for the development of the instability. As noted previously by other authors, for most stellar applications this may not be a valid approximation. We discuss this issue in detail, solving the dispersion relation of the perturbed modes for realistic values of $ν$ in the bulk of the radiative zone of the Sun and of three red giant stars at various evolutionary stages. We find that the GSF instability is triggered only in layers of extremely strong shear. In a simple case study, we also investigate the effect of a small deviation from axisymmetry or a small background magnetic field. We find that, like the viscosity, these have a stabilising effect. We conclude that this instability is probably far less efficient in transporting angular momentum than is often assumed, and may not even be present.

preprint2016arXiv

The radiative zone of the Sun and the tachocline: stability of baroclinic patterns of differential rotation

Barotropic rotation and radiative equilibrium are mutually incompatible in stars. The issue is often addressed by allowing for a meridional circulation, but this is not devoid of theoretical complications. Models of rotation in the Sun which maintain strict radiative equilibrium, making use of the observation that the Sun is not in a state of barotropic rotation, have recently been suggested. To investigate the dynamical behaviour of these solutions, we study the local stability of stratified, weakly magnetized, differentially rotating fluids to non-axisymmetric perturbations. Finite heat conductivity, kinematic viscosity, and resistivity are present. The evolution of local embedded perturbations is governed by a set of coupled, ordinary differential equations with time-dependent coefficients. Two baroclinic models of rotation for the upper radiative zone and tachocline are studied: (i) an interpolation based on helioseismology data, (ii) a theoretical solution directly compatible with radiative equilibrium. The growth of the local Goldreich-Schubert-Fricke instability appears to be suppressed, largely because of the viscosity. An extensive exploration of wavenumber space is carried out, with and without a magnetic field. Although we easily find classical local instabilities when they ought formally to be present, for the Sun the analysis reveals neither unstable solutions, nor even solutions featuring a large transient growth. We have not ruled out larger scale or nonlinear instabilities, nor have we rigorously proven local stability. But rotational configurations in close agreement with observations, generally thought to be vulnerable to the classic local Goldreich-Schubert-Fricke instability, do appear to be locally stable under rather general circumstances.

preprint2015arXiv

Inefficient Driving of Bulk Turbulence by Active Galactic Nuclei in a Hydrodynamic Model of the Intracluster Medium

Central jetted active galactic nuclei (AGN) appear to heat the core regions of the intracluster medium (ICM) in cooling-core galaxy clusters and groups, thereby preventing a cooling catastrophe. However, the physical mechanism(s) by which the directed flow of kinetic energy is thermalized throughout the ICM core remains unclear. We examine one widely discussed mechanism whereby the AGN induces subsonic turbulence in the ambient medium, the dissipation of which provides the ICM heat source. Through controlled inviscid 3-d hydrodynamic simulations, we verify that explosive AGN-like events can launch gravity waves (g-modes) into the ambient ICM which in turn decay to volume-filling turbulence. In our model, however, this process is found to be inefficient, with less than 1% of the energy injected by the AGN activity actually ending up in the turbulence of the ambient ICM. This efficiency is an order of magnitude or more too small to explain the observations of AGN-feedback in galaxy clusters and groups with short central cooling times. Atmospheres in which the g-modes are strongly trapped/confined have an even lower efficiency since, in these models, excitation of turbulence relies on the g-modes' ability to escape from the center of the cluster into the bulk ICM. Our results suggest that, if AGN-induced turbulence is indeed the mechanism by which the AGN heats the ICM core, its driving may rely on physics beyond that captured in our ideal hydrodynamic model.

preprint2014arXiv

Dynamical, biological, and anthropic consequences of equal lunar and solar angular radii

The nearly equal lunar and solar angular sizes as subtended at the Earth is generally regarded as a coincidence. This is, however, an incidental consequence of the tidal forces from these bodies being comparable. Comparable magnitudes implies strong temporal modulation, as the forcing frequencies are nearly but not precisely equal. We suggest that on the basis of paleogeographic reconstructions, in the Devonian period, when the first tetrapods appeared on land, a large tidal range would accompany these modulated tides. This would have been conducive to the formation of a network of isolated tidal pools, lending support to A.S. Romer's classic idea that the evaporation of shallow pools was an evolutionary impetus for the development of chiridian limbs in aquatic tetrapodomorphs. Romer saw this as the reason for the existence of limbs, but strong selection pressure for terrestrial navigation would have been present even if the limbs were aquatic in origin. Since even a modest difference in the Moon's angular size relative to the Sun's would lead to a qualitatively different tidal modulation, the fact that we live on a planet with a Sun and Moon of close apparent size is not entirely coincidental: it may have an anthropic basis.

preprint2014arXiv

On the high frequency spectrum of a classical accretion disc

We derive simple and explicit expressions for the high frequency spectrum of a classical accretion disc. Both stress-free and finite stress inner boundaries are considered. A classical accretion disc spectrum with a stress-free inner boundary departs from a Wien spectrum at large $ν$, scaling as $ν^{2.5}$ (as opposed to $ν^3$) times the usual exponential cut-off. If there is finite stress at the inner disc boundary, the maximum disc temperature generally occurs at this edge, even at relatively modest values of the stress. In this case, the high frequency spectrum is proportional to $ν^2$ times the exponential cut-off. If the temperature maximum is a local hot spot, instead of an axisymmetric ring, then an interior maximum produces a $ν^2$ prefactor while an edge maximum yields $ν^{1.5}$. Because of beaming effects, these latter findings should pertain to a classical relativistic disc. The asymptotics are in general robust and independent of the detailed temperature profile, provided only that the liberated free energy of differential rotation is dissipated locally, and may prove useful beyond the strict domain of classical disc theory. As observations continue to improve with time, our findings suggest the possibility of using the high energy spectral component of black hole candidates as a signature prediction of classical theory, as well as an diagnostic of the stress at the inner regions of an accretion disc.

preprint2012arXiv

The stability of stratified, rotating systems and the generation of vorticity in the Sun

We examine the linear behavior of three-dimensional Lagrangian displacements in a stratified, shearing background. The isentropic and iso-rotation surfaces of the equilibrium flow are assumed to be axisymmetric, but otherwise fully two-dimensional. Three-dimensional magnetic fields are included in the perturbation equations; however the equilibrium is assumed to be well-described by purely hydrodynamic forces. The model, in principle very general, is used to study the behavior of fluid displacements in an environment resembling the solar convection zone. Some very suggestive results emerge. All but high-latitude displacements align themselves with the observed surfaces of constant angular velocity. The tendency for the angular velocity to remain constant with depth in the bulk of the convective zone, together with other critical features of the rotation profile, emerge from little more than a visual inspection of the governing equation. In the absence of a background axial angular velocity gradient, displacements exhibit no poleward bias, suggesting that solar convection "plays-off" of prexisting shear rather than creates it. We argue that baroclinic vorticity of precisely the right order is generated at the radiative/convective zone boundary due to centrifugal distortion of equipotential surfaces that is not precisely followed by isothermal surfaces. If so, many features of the Sun's internal rotation become more clear, including: i) the general appearance of the tachocline; ii) the extension of differential rotation well into the radiative zone; iii) the abrupt change of morphology of convective zone isorotation surfaces; and iv) the inability of current numerical simulations to reproduce the solar rotation profile without imposed entropy boundary conditions.

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.

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.

preprint2010arXiv

Differential rotation in fully convective stars

Under the assumption of thermal wind balance and effective entropy mixing in constant rotation surfaces, the isorotational contours of the solar convective zone may be reproduced with great fidelity. Even at this early stage of development, this helioseismology fit may be used to put a lower bound on the midlatitude {\em radial} solar entropy gradient, which in good accord with standard mixing length theory. In this paper, we generalize this solar calculation to fully convective stars (and potentially planets), retaining the assumptions of thermal wind balance and effective entropy mixing in isorotational surfaces. It is found that each isorotation contour is of the form $R^2 = A+BΦ(r)$, where $R$ is the radius from the rotation axis, $Φ(r)$ is the (assumed spherical) gravitational potential, and $A$ and $B$ are constant along the contour. This result is applied to simple models of fully convective stars. Both solar-like surface rotation profiles (angular velocity decreasing toward the poles) as well as "antisolar" profiles (angular velocity increasing toward the poles) are modeled; the latter bear some suggestive resemblance to numerical simulations. We also perform exploratory studies of zonal surface flows similar to those seen in Jupiter and Saturn. In addition to providing a practical framework for understanding the results of large scale numerical simulations, our findings may also prove useful in dynamical calculations for which a simple but viable model for the background rotation profile in a convecting fluid is needed. Finally, our work bears directly on an important goal of the CoRoT program: to elucidate the internal structure of rotating, convecting stars.

preprint2010arXiv

Radiative and dynamic stability of a dilute plasma

We analyze the linear stability of a dilute, hot plasma, taking into account the effects of stratification and anisotropic thermal conduction. The work is motivated by attempts to understand the dynamics of the intracluster medium in galaxy clusters. We show that magnetic field configurations that nominally stabilize either the heat-flux driven buoyancy instability (associated with a positive thermal gradient) or the magnetothermal instability (negative thermal gradient) can lead to previously unrecognized g-mode overstabilities. The driving source of the overstability is either radiative cooling (positive temperature gradient) or the heat flux itself (negative temperature gradient). While the implications of these overstabilities have yet to be explored, we speculate that the cold fronts observed in many relaxed galaxy clusters may be related to their non-linear evolution.

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.

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.

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.

preprint2009arXiv

Magnetostrophic MRI in the Earth's Outer Core

We show that a simple, modified version of the Magnetorotational Instability (MRI) can develop in the outer liquid core of the Earth, in the presence of a background shear. It requires either thermal wind, or a primary instability, such as convection, to drive a weak differential rotation within the core. The force balance in the Earth's core is very unlike classical astrophysical applications of the MRI (such as gaseous disks around stars). Here, the weak differential rotation in the Earth core yields an instability by its constructive interaction with the planet's much larger rotation rate. The resulting destabilising mechanism is just strong enough to counteract stabilizing resistive effects, and produce growth on geophysically interesting timescales. We give a simple physical explanation of the instability, and show that it relies on a force balance appropriate to the Earth's core, known as magnetostrophic balance.

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.