Researcher profile

Keaton J. Burns

Keaton J. Burns contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
8topics
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

5 published item(s)

preprint2022arXiv

Inverse cascade suppression and shear layer formation in MHD turbulence subject to a guide field and misaligned rotation

Astrophysical plasmas are often subject to both rotation and large-scale background magnetic fields. Individually, each is known to two-dimensionalize the flow in the perpendicular plane. In realistic flows, both of these effects are simultaneously present and, importantly, need not be aligned. In this work, we numerically investigate three-dimensional forced magnetohydrodynamic (MHD) turbulence subject to the competing effects of global rotation and a perpendicular background magnetic field. We focus on the case of a strong background field and find that increasing the rotation rate from zero produces significant changes in the structure of the turbulent flow. Starting with a two-dimensional inverse cascade at zero rotation, the flow first transitions to a forward cascade of kinetic energy, then to a shear-layer dominated regime, and finally to a second shear-layer regime where the kinetic energy flux is strongly suppressed and the energy transfer is mediated by the induced magnetic field. We show that the first two transitions occur at distinct values of the Rossby number, and the third occurs at a distinct value of the Lehnert number. More generally, our results demonstrate that, when considering the simultaneous limits of strong rotation and a strong guide field, the order in which those limits are taken matters in the misaligned case.

preprint2020arXiv

Linearly forced fluid flow on a rotating sphere

We investigate generalized Navier-Stokes (GNS) equations that couple nonlinear advection with a generic linear instability. This analytically tractable minimal model for fluid flows driven by internal active stresses has recently been shown to permit exact solutions on a stationary two-dimensional sphere. Here, we extend the analysis to linearly driven flows on rotating spheres. We derive exact solutions of the GNS equations corresponding to time-independent zonal jets and superposed westward-propagating Rossby waves, qualitatively similar to those seen in planetary atmospheres. Direct numerical simulations with large rotation rates obtain statistically stationary states close to these exact solutions. The measured phase speeds of waves in the GNS simulations agree with analytical predictions for Rossby waves.

preprint2019arXiv

Dedalus: A Flexible Framework for Numerical Simulations with Spectral Methods

Numerical solutions of partial differential equations enable a broad range of scientific research. The Dedalus Project is a flexible, open-source, parallelized computational framework for solving general partial differential equations using spectral methods. Dedalus translates plain-text strings describing partial differential equations into efficient solvers. This paper details the numerical method that enables this translation, describes the design and implementation of the codebase, and illustrates its capabilities with a variety of example problems. The numerical method is a first-order generalized tau formulation that discretizes equations into banded matrices. This method is implemented with an object-oriented design. Classes for spectral bases and domains manage the discretization and automatic parallel distribution of variables. Discretized fields and mathematical operators are symbolically manipulated with a basic computer algebra system. Initial value, boundary value, and eigenvalue problems are efficiently solved using high-performance linear algebra, transform, and parallel communication libraries. Custom analysis outputs can also be specified in plain text and stored in self-describing portable formats. The performance of the code is evaluated with a parallel scaling benchmark and a comparison to a finite-volume code. The features and flexibility of the codebase are illustrated by solving several examples: the nonlinear Schrodinger equation on a graph, a supersonic magnetohydrodynamic vortex, quasigeostrophic flow, Stokes flow in a cylindrical annulus, normal modes of a radiative atmosphere, and diamagnetic levitation. The Dedalus code and the example problems are available online at http://dedalus-project.org/.

preprint2019arXiv

Low-frequency variability in massive stars: Core generation or surface phenomenon?

Bowman et al. (2019) reported low-frequency photometric variability in 164 O- and B-type stars observed with K2 and TESS. They interpret these motions as internal gravity waves, which could be excited stochastically by convection in the cores of these stars. The detection of internal gravity waves in massive stars would help distinguish between massive stars with convective or radiative cores, determine core size, and would provide important constraints on massive star structure and evolution. In this work, we study the observational signature of internal gravity waves generated by core convection. We calculate the \textit{wave transfer function}, which links the internal gravity wave amplitude at the base of the radiative zone to the surface luminosity variation. This transfer function varies by many orders of magnitude for frequencies $\lesssim 1 \, {\rm d}^{-1}$, and has regularly-spaced peaks near $1 \, {\rm d}^{-1}$ due to standing modes. This is inconsistent with the observed spectra which have smooth ``red noise'' profiles, without the predicted regularly-spaced peaks. The wave transfer function is only meaningful if the waves stay predominately linear. We next show that this is the case: low frequency traveling waves do not break unless their luminosity exceeds the radiative luminosity of the star; and, the observed luminosity fluctuations at high frequencies are so small that standing modes would be stable to nonlinear instability. These simple calculations suggest that the observed low-frequency photometric variability in massive stars is not due to internal gravity waves generated in the core of these stars. We finish with a discussion of (sub)surface convection, which produces low-frequency variability in low-mass stars, very similar to that observed in Bowman et al. (2019) in higher mass stars.

preprint2019arXiv

The Magnetorotational Instability Prefers Three Dimensions

The magnetorotational instability (MRI) occurs when a weak magnetic field destabilises a rotating, electrically conducting fluid with inwardly increasing angular velocity. The MRI is essential to astrophysical disk theory where the shear is typically Keplerian. Internal shear layers in stars may also be MRI unstable, and they take a wide range of profiles, including near-critical. We show that the fastest growing modes of an ideal magnetofluid are three-dimensional provided the shear rate, $S$, is near the two-dimensional onset value, $S_c$. For a Keplerian shear, three-dimensional modes are unstable above $S\approx0.10S_c$, and dominate the two-dimensional modes until $S\approx2.05S_{c}$. These three-dimensional modes dominate for shear profiles relevant to stars and at magnetic Prandtl numbers relevant to liquid-metal laboratory experiments. Significant numbers of rapidly growing three-dimensional modes remain well past $2.05S_{c}$. These finding are significant in three ways. First, weakly nonlinear theory suggests that the MRI saturates by pushing the shear rate to its critical value. This can happen for systems, like stars and laboratory experiments, that can rearrange their angular velocity profiles. Second, the non-normal character and large transient growth of MRI modes should be important whenever three-dimensionality exists. Finally, three-dimensional growth suggests direct dynamo action driven from the linear instability.