-
The solar dynamo begins near the surface
Authors:
Geoffrey M Vasil,
Daniel Lecoanet,
Kyle Augustson,
Keaton J Burns,
Jeffrey S Oishi,
Benjamin P Brown,
Nicholas Brummell,
Keith Julien
Abstract:
The Sun's magnetic dynamo cycle features a distinct pattern: a propagating region of sunspot emergence appears around 30 degrees latitude and vanishes near the equator every 11 years. Moreover, longitudinal flows called "torsional oscillations" closely shadow sunspot migration, undoubtedly sharing a common cause. Contrary to theories suggesting deep origins for these phenomena, helioseismology pin…
▽ More
The Sun's magnetic dynamo cycle features a distinct pattern: a propagating region of sunspot emergence appears around 30 degrees latitude and vanishes near the equator every 11 years. Moreover, longitudinal flows called "torsional oscillations" closely shadow sunspot migration, undoubtedly sharing a common cause. Contrary to theories suggesting deep origins for these phenomena, helioseismology pinpoints low-latitude torsional oscillations to the Sun's outer 5-10%, the "Near-Surface Shear Layer". Within this zone, inwardly increasing differential rotation coupled with a poloidal magnetic field strongly implicates the Magneto-Rotational Instability prominent in accretion-disk theory and observed in laboratory experiments. Together, these two facts prompt the general question: Is it possible that the solar dynamo is a near-surface instability? Here, we report strong affirmative evidence in stark contrast to traditional paradigms focusing on the deeper tachocline. Simple analytic estimates show that the near-surface magneto-rotational instability better explains the spatiotemporal scales of the torsional oscillations and inferred subsurface magnetic field amplitudes. State-of-the-art numerical simulations corroborate these estimates and, strikingly, reproduce hemispherical magnetic current helicity laws. The dynamo resulting from a well-understood near-surface phenomenon improves prospects for accurate predictions of full magnetic cycles and space weather, impacting Earth's electromagnetic infrastructure.
△ Less
Submitted 13 April, 2024; v1 submitted 11 April, 2024;
originally announced April 2024.
-
The photometric variability of massive stars due to gravity waves excited by core convection
Authors:
Evan H. Anders,
Daniel Lecoanet,
Matteo Cantiello,
Keaton J. Burns,
Benjamin A. Hyatt,
Emma Kaufman,
Richard H. D. Townsend,
Benjamin P. Brown,
Geoffrey M. Vasil,
Jeffrey S. Oishi,
Adam S. Jermyn
Abstract:
Massive stars die in catastrophic explosions, which seed the interstellar medium with heavy elements and produce neutron stars and black holes. Predictions of the explosion's character and the remnant mass depend on models of the star's evolutionary history. Models of massive star interiors can be empirically constrained by asteroseismic observations of gravity wave oscillations. Recent photometri…
▽ More
Massive stars die in catastrophic explosions, which seed the interstellar medium with heavy elements and produce neutron stars and black holes. Predictions of the explosion's character and the remnant mass depend on models of the star's evolutionary history. Models of massive star interiors can be empirically constrained by asteroseismic observations of gravity wave oscillations. Recent photometric observations reveal a ubiquitous red noise signal on massive main sequence stars; a hypothesized source of this noise is gravity waves driven by core convection. We present the first 3D simulations of massive star convection extending from the star's center to near its surface, with realistic stellar luminosities. Using these simulations, we make the first prediction of photometric variability due to convectively-driven gravity waves at the surfaces of massive stars, and find that gravity waves produce photometric variability of a lower amplitude and lower characteristic frequency than the observed red noise. We infer that the photometric signal of gravity waves excited by core convection is below the noise limit of current observations, so the red noise must be generated by an alternative process.
△ Less
Submitted 13 June, 2023;
originally announced June 2023.
-
The Puzzling Structure of Solar Convection: Window into the Dynamo
Authors:
Nicholas A. Featherstone,
Kyle C. Augustson,
Jonathan M. Aurnou,
Catherine Blume,
Benjamin P. Brown,
Nicholas Brummell,
Keaton J. Burns,
Michael A. Calkins,
Maria Camisassa,
Mausumi Dikpati,
Yuhong Fan,
J. R. Fuentes,
Gustavo Guerrero,
Bradley W. Hindman,
Keith Julien,
Irina N. Kitiashvili,
Lydia Korre,
Daniel Lecoanet,
Bhishek Manek,
Loren Matilsky,
Mark Miesch,
Nicholas J. Nelson,
Jeffrey S. Oishi,
Whitney T. Powers,
Matthias Rempel
, et al. (3 additional authors not shown)
Abstract:
The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective…
▽ More
The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective flows. Those flows are in-turn thought to sustain both the Sun's turbulent EMF and the large-scale flows of differential rotation and meridional circulation suspected of influencing the dynamo's organization and timing. Continued progress in this area is complicated by (i) inconsistencies between helioseismic measurements of convective and meridional flow made with different techniques and instruments, and (ii) a lack of high-latitude data for convection, differential rotation, and meridional flow. We suggest that the path forward to resolving these difficulties is twofold. First, the acquisition of long-term helioseismic and emissivity measurements obtained from a polar vantage point is vital to complete our picture of the Sun's outer convection zone. Second, sustained and expanded investment in theory-oriented and combined theory/observational research initiatives will be crucial to fully exploit these new observations and to resolve inconsistencies between existing measurements.
△ Less
Submitted 15 May, 2023;
originally announced May 2023.
-
The Stability of Prendergast Magnetic Fields
Authors:
Emma Kaufman,
Daniel Lecoanet,
Evan H. Anders,
Benjamin P. Brown,
Geoffrey M. Vasil,
Jeffrey S. Oishi,
Keaton J. Burns
Abstract:
Convection in massive main sequence stars generates large scale magnetic fields in their cores which persists as they evolve up the red giant branch. The remnants of these fields may take the form of the Prendergast magnetic field, a combination of poloidal and toroidal field components which are expected to stabilize each other. Previous analytic and numerical calculations did not find any eviden…
▽ More
Convection in massive main sequence stars generates large scale magnetic fields in their cores which persists as they evolve up the red giant branch. The remnants of these fields may take the form of the Prendergast magnetic field, a combination of poloidal and toroidal field components which are expected to stabilize each other. Previous analytic and numerical calculations did not find any evidence for instability of the Prendergast field over short timescales. In this paper, we present numerical simulations which show a long timescale, linear instability of this magnetic field. We find the instability to be robust to changes in boundary conditions and it is not stabilized by strong stable stratification. The instability is a resistive instability, and the growth rate has a power-law dependence on the resistivity, in which the growth rate decreases as the resistivity decreases. We estimate the growth rate of the instability in stars by extrapolating this power-law to stellar values of the resistivity. The instability is sufficiently rapid to destabilize the magnetic field on timescales shorter than the stellar evolution timescale, indicating that the Prendergast field is not a good model to use in studies of magnetic fields in stars.
△ Less
Submitted 28 September, 2022;
originally announced September 2022.
-
The Pencil Code, a modular MPI code for partial differential equations and particles: multipurpose and multiuser-maintained
Authors:
A. Brandenburg,
A. Johansen,
P. A. Bourdin,
W. Dobler,
W. Lyra,
M. Rheinhardt,
S. Bingert,
N. E. L. Haugen,
A. Mee,
F. Gent,
N. Babkovskaia,
C. -C. Yang,
T. Heinemann,
B. Dintrans,
D. Mitra,
S. Candelaresi,
J. Warnecke,
P. J. Käpylä,
A. Schreiber,
P. Chatterjee,
M. J. Käpylä,
X. -Y. Li,
J. Krüger,
J. R. Aarnes,
G. R. Sarson
, et al. (12 additional authors not shown)
Abstract:
The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the…
▽ More
The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the framework is general and can also be applied to situations not related to hydrodynamics or even PDEs, for example when just the message passing interface or input/output strategies of the code are to be used. The code can also evolve Lagrangian (inertial and noninertial) particles, their coagulation and condensation, as well as their interaction with the fluid.
△ Less
Submitted 17 September, 2020;
originally announced September 2020.
-
Single-hemisphere dynamos in M-dwarf stars
Authors:
Benjamin P. Brown,
Jeffrey S. Oishi,
Geoffrey M. Vasil,
Daniel Lecoanet,
Keaton J. Burns
Abstract:
M-dwarf stars below a certain mass are convective from their cores to their photospheres. These fully convective objects are extremely numerous, very magnetically active, and the likely hosts of many exoplanets. Here we study, for the first time, dynamo action in simulations of stratified, rotating fully convective M-dwarf stars. Importantly, we use new techniques to capture the correct full ball…
▽ More
M-dwarf stars below a certain mass are convective from their cores to their photospheres. These fully convective objects are extremely numerous, very magnetically active, and the likely hosts of many exoplanets. Here we study, for the first time, dynamo action in simulations of stratified, rotating fully convective M-dwarf stars. Importantly, we use new techniques to capture the correct full ball geometry down to the center of the star. We find surprising dynamo states in these systems, with the global-scale mean fields confined strongly to a single hemisphere, in contrast to prior stellar dynamo solutions. These hemispheric-dynamo stars are likely to have profoundly different interactions with their surroundings, with important implications for exoplanet habitability and stellar spindown.
△ Less
Submitted 5 August, 2020;
originally announced August 2020.
-
The Magnetorotational Instability Prefers Three Dimensions
Authors:
Jeffrey S. Oishi,
Geoffrey M. Vasil,
Morgan Baxter,
Andrew Swan,
Keaton J. Burns,
Daniel Lecoanet,
Benjamin P. Brown
Abstract:
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…
▽ More
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.
△ Less
Submitted 2 December, 2019;
originally announced December 2019.
-
Dedalus: A Flexible Framework for Numerical Simulations with Spectral Methods
Authors:
Keaton J. Burns,
Geoffrey M. Vasil,
Jeffrey S. Oishi,
Daniel Lecoanet,
Benjamin P. Brown
Abstract:
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 meth…
▽ More
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/.
△ Less
Submitted 2 December, 2019; v1 submitted 24 May, 2019;
originally announced May 2019.
-
The "Sphered Cube": A New Method for the Solution of Partial Differential Equations in Cubical Geometry
Authors:
Keaton J. Burns,
Daniel Lecoanet,
Geoffrey M. Vasil,
Jeffrey S. Oishi,
Benjamin P. Brown
Abstract:
A new gridding technique for the solution of partial differential equations in cubical geometry is presented. The method is based on volume penalization, allowing for the imposition of a cubical geometry inside of its circumscribing sphere. By choosing to embed the cube inside of the sphere, one obtains a discretization that is free of any sharp edges or corners. Taking full advantage of the simpl…
▽ More
A new gridding technique for the solution of partial differential equations in cubical geometry is presented. The method is based on volume penalization, allowing for the imposition of a cubical geometry inside of its circumscribing sphere. By choosing to embed the cube inside of the sphere, one obtains a discretization that is free of any sharp edges or corners. Taking full advantage of the simple geometry of the sphere, spectral bases based on spin-weighted spherical harmonics and Jacobi polynomials, which properly capture the regularity of scalar, vector and tensor components in spherical coordinates, can be applied to obtain moderately efficient and accurate numerical solutions of partial differential equations in the cube. This technique demonstrates the advantages of these bases over other methods for solving PDEs in spherical coordinates. We present results for a test case of incompressible hydrodynamics in cubical geometry: Rayleigh-Bénard convection with fully Dirichlet boundary conditions. Analysis of the simulations provides what is, to our knowledge, the first result on the scaling of the heat flux with the thermal forcing for this type of convection in a cube in a sphere.
△ Less
Submitted 29 March, 2019;
originally announced March 2019.
-
Breezing through the space environment of Barnard's Star b
Authors:
Julián D. Alvarado-Gómez,
Cecilia Garraffo,
Jeremy J. Drake,
Benjamin P. Brown,
Jeffrey S. Oishi,
Sofia P. Moschou,
Ofer Cohen
Abstract:
A physically realistic stellar wind model based on Alfvén wave dissipation has been used to simulate the wind from Barnard's Star and to estimate the conditions at the location of its recently discovered planetary companion. Such models require knowledge of the stellar surface magnetic field that is currently unknown for Barnard's Star. We circumvent this by considering the observed field distribu…
▽ More
A physically realistic stellar wind model based on Alfvén wave dissipation has been used to simulate the wind from Barnard's Star and to estimate the conditions at the location of its recently discovered planetary companion. Such models require knowledge of the stellar surface magnetic field that is currently unknown for Barnard's Star. We circumvent this by considering the observed field distributions of three different stars that constitute admissible magnetic proxies of this object. Under these considerations, Barnard's Star b experiences less intense wind pressure than the much more close-in planet Proxima~b and the planets of the TRAPPIST-1 system. The milder wind conditions are more a result of its much greater orbital distance rather than in differences in the surface magnetic field strengths of the host stars. The dynamic pressure experienced by the planet is comparable to present-day Earth values, but it can undergo variations by factors of several during current sheet crossings in each orbit. The magnetospause standoff distance would be $\sim$\,$20 - 40$\,\% smaller than that of the Earth for an equivalent planetary magnetic field strength.
△ Less
Submitted 29 March, 2019; v1 submitted 1 January, 2019;
originally announced January 2019.
-
Predicting the Rossby number in convective experiments
Authors:
Evan H. Anders,
Cathryn M. Manduca,
Benjamin P. Brown,
Jeffrey S. Oishi,
Geoffrey M. Vasil
Abstract:
The Rossby number is a crucial parameter describing the degree of rotational constraint on the convective dynamics in stars and planets. However, it is not an input to computational models of convection but must be measured ex post facto. Here, we report the discovery of a new quantity, the Predictive Rossby number, which is both tightly correlated with the Rossby number and specified in terms of…
▽ More
The Rossby number is a crucial parameter describing the degree of rotational constraint on the convective dynamics in stars and planets. However, it is not an input to computational models of convection but must be measured ex post facto. Here, we report the discovery of a new quantity, the Predictive Rossby number, which is both tightly correlated with the Rossby number and specified in terms of common inputs to numerical models. The Predictive Rossby number can be specified independent of Rayleigh number, allowing suites of numerical solutions to separate the degree of rotational constraint from the strength of the driving of convection. We examine the scaling of convective transport in terms of the Nusselt number and the degree of turbulence in terms of the Reynolds number of the flow, and we find scaling laws nearly identical to those in nonrotational convection at low Rossby number when the Predictive Rossby number is held constant. Finally, we describe the boundary layers as a function of increasing turbulence at constant Rossby number.
△ Less
Submitted 15 January, 2019; v1 submitted 11 December, 2018;
originally announced December 2018.
-
Accelerated evolution of convective simulations
Authors:
Evan H. Anders,
Benjamin P. Brown,
Jeffrey S. Oishi
Abstract:
High Peclet number, turbulent convection is a classic system with a large timescale separation between flow speeds and the thermal relaxation time. In this paper, we present a method of fast-forwarding through the long thermal relaxation of convective simulations, and we test the validity of this method. This accelerated evolution (AE) method involves measuring the dynamics of convection early in…
▽ More
High Peclet number, turbulent convection is a classic system with a large timescale separation between flow speeds and the thermal relaxation time. In this paper, we present a method of fast-forwarding through the long thermal relaxation of convective simulations, and we test the validity of this method. This accelerated evolution (AE) method involves measuring the dynamics of convection early in a simulation and using its characteristics to adjust the mean thermodynamic profile within the domain towards its evolved state. We study Rayleigh-Bénard convection as a test case for AE. Evolved flow properties of AE solutions are measured to be within a few percent of solutions which are reached through standard evolution (SE) over a full thermal timescale. At the highest values of the Rayleigh number at which we compare SE and AE, we find that AE solutions require roughly an order of magnitude fewer computing hours to evolve than SE solutions.
△ Less
Submitted 21 August, 2018; v1 submitted 17 July, 2018;
originally announced July 2018.
-
Tensor calculus in spherical coordinates using Jacobi polynomials, Part-II: Implementation and Examples
Authors:
Daniel Lecoanet,
Geoffrey M. Vasil,
Keaton J. Burns,
Benjamin P. Brown,
Jeffrey S. Oishi
Abstract:
We present a simulation code which can solve broad ranges of partial differential equations in a full sphere. The code expands tensorial variables in a spectral series of spin-weighted spherical harmonics in the angular directions and a scaled Jacobi polynomial basis in the radial direction, as described in Part-I. Nonlinear terms are calculated by transforming from the coefficients in the spectra…
▽ More
We present a simulation code which can solve broad ranges of partial differential equations in a full sphere. The code expands tensorial variables in a spectral series of spin-weighted spherical harmonics in the angular directions and a scaled Jacobi polynomial basis in the radial direction, as described in Part-I. Nonlinear terms are calculated by transforming from the coefficients in the spectral series to the value of each quantity on the physical grid, where it is easy to calculate products and perform other local operations. The expansion makes it straightforward to solve equations in tensor form (i.e., without decomposition into scalars). We propose and study several unit tests which demonstrate the code can accurately solve linear problems, implement boundary conditions, and transform between spectral and physical space. We then run a series of benchmark problems proposed in Marti et al (2014), implementing the hydrodynamic and magnetohydrodynamic equations. We are able to calculate more accurate solutions than reported in Marti et al 2014 by running at higher spatial resolution and using a higher-order timestepping scheme. We find the rotating convection and convective dynamo benchmark problems depend sensitively on details of timestepping and data analysis. We also demonstrate that in low resolution simulations of the dynamo problem, small changes in a numerical scheme can lead to large changes in the solution. To aid future comparison to these benchmarks, we include the source code used to generate the data, as well as the data and analysis scripts used to generate the figures.
△ Less
Submitted 24 April, 2018;
originally announced April 2018.
-
Convective dynamics and disequilibrium chemistry in the atmospheres of giant planets and brown dwarfs
Authors:
Baylee Bordwell,
Benjamin P. Brown,
Jeffrey S. Oishi
Abstract:
Disequilibrium chemical processes have a large effect upon the spectra of substellar objects. To study these effects, dynamical disequilibrium has been parameterized using the quench and eddy diffusion approximations, but little work has been done to explore how these approximations perform under realistic planetary conditions in different dynamical regimes. As a first step in addressing this prob…
▽ More
Disequilibrium chemical processes have a large effect upon the spectra of substellar objects. To study these effects, dynamical disequilibrium has been parameterized using the quench and eddy diffusion approximations, but little work has been done to explore how these approximations perform under realistic planetary conditions in different dynamical regimes. As a first step in addressing this problem, we study the localized, small scale convective dynamics of planetary atmospheres by direct numerical simulation of fully compressible hydrodynamics with reactive tracers using the Dedalus code. Using polytropically-stratified, plane parallel atmospheres in 2- and 3-D, we explore the quenching behavior of different abstract chemical species as a function of the dynamical conditions of the atmosphere as parameterized by the Rayleigh number. We find that in both 2- and 3-D, chemical species quench deeper than would be predicted based on simple mixing length arguments. Instead, it is necessary to employ length scales based on the chemical equilibrium profile of the reacting species in order to predict quench points and perform chemical kinetics modeling in 1-D. Based on the results of our simulations, we provide a new length scale, derived from the chemical scale height, which can be used to perform these calculations. This length scale is simple to calculate from known chemical data and makes reasonable predictions for our dynamical simulations.
△ Less
Submitted 8 February, 2018;
originally announced February 2018.
-
Perspectives on Reproducibility and Sustainability of Open-Source Scientific Software from Seven Years of the Dedalus Project
Authors:
Jeffrey S. Oishi,
Benjamin P. Brown,
Keaton J. Burns,
Daniel Lecoanet,
Geoffrey M. Vasil
Abstract:
As the Science Mission Directorate contemplates establishing an open code policy, we consider it timely to share our experiences as the developers of the open-source partial differential equation solver Dedalus. Dedalus is a flexible framework for solving partial differential equations. Its development team primarily uses it for studying stellar and planetary astrophysics. Dedalus was developed or…
▽ More
As the Science Mission Directorate contemplates establishing an open code policy, we consider it timely to share our experiences as the developers of the open-source partial differential equation solver Dedalus. Dedalus is a flexible framework for solving partial differential equations. Its development team primarily uses it for studying stellar and planetary astrophysics. Dedalus was developed originally for astrophysical fluid dynamics (AFD), though it has found a much broader user base, including applied mathematicians, plasma physicists, and oceanographers. Here, we will focus on issues related to open-source software from the perspective of AFD. We use the term AFD with the understanding that astrophysics simulations are inherently multi-physics: fluid dynamics coupled with some combination of gravitational dynamics, radiation transfer, relativity, and magnetic fields. In practice, a few well-known open-source simulation packages represent a large fraction of published work in the field. However, we will argue that an open-code policy should encompass not just these large simulation codes, but also the input files and analysis scripts. It is our interest that NASA adopt an open-code policy because without it, reproducibility in computational science is needlessly hampered.
△ Less
Submitted 24 January, 2018;
originally announced January 2018.
-
The weakly nonlinear magnetorotational instability in a local geometry
Authors:
S. E. Clark,
Jeffrey S. Oishi
Abstract:
The magnetorotational instability (MRI) is a fundamental process of accretion disk physics, but its saturation mechanism remains poorly understood despite considerable theoretical and computational effort. We present a multiple scales analysis of the non-ideal MRI in the weakly nonlinear regime -- that is, when the most unstable MRI mode has a growth rate asymptotically approaching zero from above…
▽ More
The magnetorotational instability (MRI) is a fundamental process of accretion disk physics, but its saturation mechanism remains poorly understood despite considerable theoretical and computational effort. We present a multiple scales analysis of the non-ideal MRI in the weakly nonlinear regime -- that is, when the most unstable MRI mode has a growth rate asymptotically approaching zero from above. Here, we develop our theory in a local, Cartesian channel. Our results confirm the finding by Umurhan et al. (2007) that the perturbation amplitude follows a Ginzburg-Landau equation. We further find that the Ginzburg-Landau equation will arise for the local MRI system with shear-periodic boundary conditions when the effects of ambipolar diffusion are considered. A detailed force balance for the saturated azimuthal velocity and vertical magnetic field demonstrates that even when diffusive effects are important, the bulk flow saturates via the combined processes of reducing the background shear and rearranging and strengthening the background vertical magnetic field. We directly simulate the Ginzburg-Landau amplitude evolution for our system and demonstrate the pattern formation our model predicts on long length and time scales. We compare the weakly nonlinear theory results to a direct numerical simulation of the MRI in a thin-gap Taylor Couette flow.
△ Less
Submitted 28 April, 2017; v1 submitted 5 October, 2016;
originally announced October 2016.
-
The weakly nonlinear magnetorotational instability in a global, cylindrical Taylor-Couette flow
Authors:
S. E. Clark,
Jeffrey S. Oishi
Abstract:
We conduct a global, weakly nonlinear analysis of the magnetorotational instability (MRI) in a Taylor-Couette flow. This is a multiscale perturbative treatment of the nonideal, axisymmetric MRI near threshold, subject to realistic radial boundary conditions and cylindrical geometry. We analyze both the standard MRI, initialized by a constant vertical background magnetic field, and the helical MRI,…
▽ More
We conduct a global, weakly nonlinear analysis of the magnetorotational instability (MRI) in a Taylor-Couette flow. This is a multiscale perturbative treatment of the nonideal, axisymmetric MRI near threshold, subject to realistic radial boundary conditions and cylindrical geometry. We analyze both the standard MRI, initialized by a constant vertical background magnetic field, and the helical MRI, with an azimuthal background field component. This is the first weakly nonlinear analysis of the MRI in a global Taylor-Couette geometry, as well as the first weakly nonlinear analysis of the helical MRI. We find that the evolution of the amplitude of the standard MRI is described by a real Ginzburg-Landau equation (GLE), while the amplitude of the helical MRI takes the form of a complex GLE. This suggests that the saturated state of the helical MRI may itself be unstable on long spatial and temporal scales.
△ Less
Submitted 28 April, 2017; v1 submitted 5 October, 2016;
originally announced October 2016.
-
Hybrid Adaptive Ray-Moment Method (HARM$^2$): A Highly Parallel Method for Radiation Hydrodynamics on Adaptive Grids
Authors:
Anna L. Rosen,
Mark R. Krumholz,
Jeffrey S. Oishi,
Aaron T. Lee,
Richard I. Klein
Abstract:
We present a highly-parallel multi-frequency hybrid radiation hydrodynamics algorithm that combines a spatially-adaptive long characteristics method for the radiation field from point sources with a moment method that handles the diffuse radiation field produced by a volume-filling fluid. Our Hybrid Adaptive Ray-Moment Method (HARM$^2$) operates on patch-based adaptive grids, is compatible with as…
▽ More
We present a highly-parallel multi-frequency hybrid radiation hydrodynamics algorithm that combines a spatially-adaptive long characteristics method for the radiation field from point sources with a moment method that handles the diffuse radiation field produced by a volume-filling fluid. Our Hybrid Adaptive Ray-Moment Method (HARM$^2$) operates on patch-based adaptive grids, is compatible with asynchronous time stepping, and works with any moment method. In comparison to previous long characteristics methods, we have greatly improved the parallel performance of the adaptive long-characteristics method by developing a new completely asynchronous and non-blocking communication algorithm. As a result of this improvement, our implementation achieves near-perfect scaling up to $\mathcal{O}(10^3)$ processors on distributed memory machines. We present a series of tests to demonstrate the accuracy and performance of the method.
△ Less
Submitted 6 July, 2016;
originally announced July 2016.
-
Turbulent Chemical Diffusion in Convectively Bounded Carbon Flames
Authors:
Daniel Lecoanet,
Josiah Schwab,
Eliot Quataert,
Lars Bildsten,
F. X. Timmes,
Keaton J. Burns,
Geoffrey M. Vasil,
Jeffrey S. Oishi,
Benjamin P. Brown
Abstract:
It has been proposed that mixing induced by convective overshoot can disrupt the inward propagation of carbon deflagrations in super-asymptotic giant branch stars. To test this theory, we study an idealized model of convectively bounded carbon flames with 3D hydrodynamic simulations of the Boussinesq equations using the pseudospectral code Dedalus. Because the flame propagation timescale is much l…
▽ More
It has been proposed that mixing induced by convective overshoot can disrupt the inward propagation of carbon deflagrations in super-asymptotic giant branch stars. To test this theory, we study an idealized model of convectively bounded carbon flames with 3D hydrodynamic simulations of the Boussinesq equations using the pseudospectral code Dedalus. Because the flame propagation timescale is much longer than the convection timescale, we approximate the flame as fixed in space, and only consider its effects on the buoyancy of the fluid. By evolving a passive scalar field, we derive a {\it turbulent} chemical diffusivity produced by the convection as a function of height, $D_{\rm t}(z)$. Convection can stall a flame if the chemical mixing timescale, set by the turbulent chemical diffusivity, $D_{\rm t}$, is shorter than the flame propagation timescale, set by the thermal diffusivity, $κ$, i.e., when $D_{\rm t}>κ$. However, we find $D_{\rm t}<κ$ for most of the flame because convective plumes are not dense enough to penetrate into the flame. Extrapolating to realistic stellar conditions, this implies that convective mixing cannot stall a carbon flame and that "hybrid carbon-oxygen-neon" white dwarfs are not a typical product of stellar evolution.
△ Less
Submitted 25 September, 2016; v1 submitted 29 March, 2016;
originally announced March 2016.
-
Tensor calculus in polar coordinates using Jacobi polynomials
Authors:
Geoffrey M. Vasil,
Keaton J. Burns,
Daniel Lecoanet,
Sheehan Olver,
Benjamin P. Brown,
Jeffrey S. Oishi
Abstract:
Spectral methods are an efficient way to solve partial differential equations on domains possessing certain symmetries. The utility of a method depends strongly on the choice of spectral basis. In this paper we describe a set of bases built out of Jacobi polynomials, and associated operators for solving scalar, vector, and tensor partial differential equations in polar coordinates on a unit disk.…
▽ More
Spectral methods are an efficient way to solve partial differential equations on domains possessing certain symmetries. The utility of a method depends strongly on the choice of spectral basis. In this paper we describe a set of bases built out of Jacobi polynomials, and associated operators for solving scalar, vector, and tensor partial differential equations in polar coordinates on a unit disk. By construction, the bases satisfy regularity conditions at r=0 for any tensorial field. The coordinate singularity in a disk is a prototypical case for many coordinate singularities. The work presented here extends to other geometries. The operators represent covariant derivatives, multiplication by azimuthally symmetric functions, and the tensorial relationship between fields. These arise naturally from relations between classical orthogonal polynomials, and form a Heisenberg algebra. Other past work uses more specific polynomial bases for solving equations in polar coordinates. The main innovation in this paper is to use a larger set of possible bases to achieve maximum bandedness of linear operations. We provide a series of applications of the methods, illustrating their ease-of-use and accuracy.
△ Less
Submitted 15 August, 2016; v1 submitted 25 September, 2015;
originally announced September 2015.
-
A Validated Nonlinear Kelvin-Helmholtz Benchmark for Numerical Hydrodynamics
Authors:
Daniel Lecoanet,
Michael McCourt,
Eliot Quataert,
Keaton J. Burns,
Geoffrey M. Vasil,
Jeffrey S. Oishi,
Benjamin P. Brown,
James M. Stone,
Ryan M. O'Leary
Abstract:
The nonlinear evolution of the Kelvin-Helmholtz instability is a popular test for code verification. To date, most Kelvin-Helmholtz problems discussed in the literature are ill-posed: they do not converge to any single solution with increasing resolution. This precludes comparisons among different codes and severely limits the utility of the Kelvin-Helmholtz instability as a test problem. The lack…
▽ More
The nonlinear evolution of the Kelvin-Helmholtz instability is a popular test for code verification. To date, most Kelvin-Helmholtz problems discussed in the literature are ill-posed: they do not converge to any single solution with increasing resolution. This precludes comparisons among different codes and severely limits the utility of the Kelvin-Helmholtz instability as a test problem. The lack of a reference solution has led various authors to assert the accuracy of their simulations based on ad-hoc proxies, e.g., the existence of small-scale structures. This paper proposes well-posed Kelvin-Helmholtz problems with smooth initial conditions and explicit diffusion. We show that in many cases numerical errors/noise can seed spurious small-scale structure in Kelvin-Helmholtz problems. We demonstrate convergence to a reference solution using both Athena, a Godunov code, and Dedalus, a pseudo-spectral code. Problems with constant initial density throughout the domain are relatively straightforward for both codes. However, problems with an initial density jump (which are the norm in astrophysical systems) exhibit rich behavior and are more computationally challenging. In the latter case, Athena simulations are prone to an instability of the inner rolled-up vortex; this instability is seeded by grid-scale errors introduced by the algorithm, and disappears as resolution increases. Both Athena and Dedalus exhibit late-time chaos. Inviscid simulations are riddled with extremely vigorous secondary instabilities which induce more mixing than simulations with explicit diffusion. Our results highlight the importance of running well-posed test problems with demonstrated convergence to a reference solution. To facilitate future comparisons, we include the resolved, converged solutions to the Kelvin-Helmholtz problems in this paper in machine-readable form.
△ Less
Submitted 27 December, 2021; v1 submitted 11 September, 2015;
originally announced September 2015.
-
Self-generated turbulence in magnetic reconnection
Authors:
Jeffrey S. Oishi,
Mordecai-Mark Mac Low,
David C. Collins,
Moeko Tamura
Abstract:
Classical Sweet-Parker models of reconnection predict that reconnection rates depend inversely on the resistivity, usually parameterized using the dimensionless Lundquist number ($\Lund$). We describe magnetohydrodynamic (MHD) simulations using a static, nested grid that show the development of a three-dimensional instability in the plane of a current sheet between reversing field lines without a…
▽ More
Classical Sweet-Parker models of reconnection predict that reconnection rates depend inversely on the resistivity, usually parameterized using the dimensionless Lundquist number ($\Lund$). We describe magnetohydrodynamic (MHD) simulations using a static, nested grid that show the development of a three-dimensional instability in the plane of a current sheet between reversing field lines without a guide field. The instability leads to rapid reconnection of magnetic field lines at a rate independent of $\Lund$ over at least the range $3.2\times 10^3 \lesssim \Lund \lesssim 3.2 \times 10^5$ resolved by the simulations. We find that this instability occurs even for cases with $\Lund \lesssim 10^4$ that in our models appear stable to the recently described, two-dimensional, plasmoid instability. Our results suggest that three-dimensional, MHD processes alone produce fast (resistivity independent) reconnection without recourse to kinetic effects or external turbulence. The unstable reconnection layers provide a self-consistent environment in which the extensively studied turbulent reconnection process can occur.
△ Less
Submitted 18 May, 2015;
originally announced May 2015.
-
Numerical Simulations of Internal Wave Generation by Convection in Water
Authors:
Daniel Lecoanet,
Michael Le Bars,
Keaton J. Burns,
Geoffrey M. Vasil,
Benjamin P. Brown,
Eliot Quataert,
Jeffrey S. Oishi
Abstract:
Water's density maximum at 4C makes it well suited to study internal gravity wave excitation by convection: an increasing temperature profile is unstable to convection below 4C, but stably stratified above 4C. We present numerical simulations of a water-like fluid near its density maximum in a two dimensional domain. We successfully model the damping of waves in the simulations using linear theory…
▽ More
Water's density maximum at 4C makes it well suited to study internal gravity wave excitation by convection: an increasing temperature profile is unstable to convection below 4C, but stably stratified above 4C. We present numerical simulations of a water-like fluid near its density maximum in a two dimensional domain. We successfully model the damping of waves in the simulations using linear theory, provided we do not take the weak damping limit typically used in the literature. In order to isolate the physical mechanism exciting internal waves, we use the novel spectral code Dedalus to run several simplified model simulations of our more detailed simulation. We use data from the full simulation as source terms in two simplified models of internal wave excitation by convection: bulk excitation by convective Reynolds stresses, and interface forcing via the mechanical oscillator effect. We find excellent agreement between the waves generated in the full simulation and the simplified simulation implementing the bulk excitation mechanism. The interface forcing simulations over excite high frequency waves because they assume the excitation is by the "impulsive" penetration of plumes, which spreads energy to high frequencies. However, we find the real excitation is instead by the "sweeping" motion of plumes parallel to the interface. Our results imply that the bulk excitation mechanism is a very accurate heuristic for internal wave generation by convection.
△ Less
Submitted 17 June, 2015; v1 submitted 9 December, 2014;
originally announced December 2014.
-
Conduction in low Mach number flows: Part I Linear & weakly nonlinear regimes
Authors:
Daniel Lecoanet,
Benjamin P. Brown,
Ellen G. Zweibel,
Keaton J. Burns,
Jeffrey S. Oishi,
Geoffrey M. Vasil
Abstract:
Thermal conduction is an important energy transfer and damping mechanism in astrophysical flows. Fourier's law - the heat flux is proportional to the negative temperature gradient, leading to temperature diffusion - is a well-known empirical model of thermal conduction. However, entropy diffusion has emerged as an alternative thermal conduction model, despite not ensuring the monotonicity of entro…
▽ More
Thermal conduction is an important energy transfer and damping mechanism in astrophysical flows. Fourier's law - the heat flux is proportional to the negative temperature gradient, leading to temperature diffusion - is a well-known empirical model of thermal conduction. However, entropy diffusion has emerged as an alternative thermal conduction model, despite not ensuring the monotonicity of entropy. This paper investigates the differences between temperature and entropy diffusion for both linear internal gravity waves and weakly nonlinear convection. In addition to simulating the two thermal conduction models with the fully compressible Navier-Stokes equations, we also study their effects in the reduced, "sound-proof" anelastic and pseudo-incompressible equations. We find that in the linear and weakly nonlinear regimes, temperature and entropy diffusion give quantitatively similar results, although there are some larger errors in the pseudo-incompressible equations with temperature diffusion due to inaccuracies in the equation of state. Extrapolating our weakly nonlinear results, we speculate that differences between temperature and entropy diffusion might become more important for strongly turbulent convection.
△ Less
Submitted 20 October, 2014;
originally announced October 2014.
-
Radial Stresses and Energy Transport in Accretion Disks
Authors:
Alexander Hubbard,
Colin P. McNally,
Jeffrey S. Oishi,
Wladimir Lyra,
Mordecai-Mark Mac Low
Abstract:
Early in the study of viscous accretion disks it was realized that energy transfers from distant sources must be important, not least because the flow at the disk midplane in the bulk of the disk is likely outwards, out of the gravitational potential well. If the source of the viscosity is powered by accretion, such as in the case of the magneto-rotational instability, such distant energy sources…
▽ More
Early in the study of viscous accretion disks it was realized that energy transfers from distant sources must be important, not least because the flow at the disk midplane in the bulk of the disk is likely outwards, out of the gravitational potential well. If the source of the viscosity is powered by accretion, such as in the case of the magneto-rotational instability, such distant energy sources must lie in the innermost regions of the disk, where accretion occurs even at the midplane. We argue here that modulations in this energy supply can alter the accretion rate on dynamical, rather than far longer viscous, time scales. This means that both the steady state value of and fluctuations in the inner disk's accretion rate, depending on the details of the inner boundary condition and occurring on the inner disk's rapid evolution time, can affect the outer disk. This is particularly interesting because observations have shown that disk accretion is not steady (e.g.~EX Lupi type objects). We also note that the power supplied to shearing boxes is set by the boxes themselves rather than the physical energy fluxes in a global disk. That is, their saturated magnetic field is not subject to the full set of energy constraints present in an actual disk. Our analysis suggests that large scale radial transport of energy has a critical impact on the evolution and variability of accretion disks.
△ Less
Submitted 12 December, 2014; v1 submitted 15 May, 2014;
originally announced May 2014.
-
Enzo: An Adaptive Mesh Refinement Code for Astrophysics
Authors:
The Enzo Collaboration,
Greg L. Bryan,
Michael L. Norman,
Brian W. O'Shea,
Tom Abel,
John H. Wise,
Matthew J. Turk,
Daniel R. Reynolds,
David C. Collins,
Peng Wang,
Samuel W. Skillman,
Britton Smith,
Robert P. Harkness,
James Bordner,
Ji-hoon Kim,
Michael Kuhlen,
Hao Xu,
Nathan Goldbaum,
Cameron Hummels,
Alexei G. Kritsuk,
Elizabeth Tasker,
Stephen Skory,
Christine M. Simpson,
Oliver Hahn,
Jeffrey S. Oishi
, et al. (4 additional authors not shown)
Abstract:
This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gra…
▽ More
This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.
△ Less
Submitted 8 July, 2013;
originally announced July 2013.
-
Magnetic Fields in Population III Star Formation
Authors:
Matthew J. Turk,
Jeffrey S. Oishi,
Tom Abel,
Greg Bryan
Abstract:
We study the buildup of magnetic fields during the formation of Population III star-forming regions, by conducting cosmological simulations from realistic initial conditions and varying the Jeans resolution. To investigate this in detail, we start simulations from identical initial conditions, mandating 16, 32 and 64 zones per Jeans length, and studied the variation in their magnetic field amplifi…
▽ More
We study the buildup of magnetic fields during the formation of Population III star-forming regions, by conducting cosmological simulations from realistic initial conditions and varying the Jeans resolution. To investigate this in detail, we start simulations from identical initial conditions, mandating 16, 32 and 64 zones per Jeans length, and studied the variation in their magnetic field amplification. We find that, while compression results in some amplification, turbulent velocity fluctuations driven by the collapse can further amplify an initially weak seed field via dynamo action, provided there is sufficient numerical resolution to capture vortical motions (we find this requirement to be 64 zones per Jeans length, slightly larger than, but consistent with previous work run with more idealized collapse scenarios). We explore saturation of amplification of the magnetic field, which could potentially become dynamically important in subsequent, fully-resolved calculations. We have also identified a relatively surprising phenomena that is purely hydrodynamic: the higher-resolved simulations possess substantially different characteristics, including higher infall-velocity, increased temperatures inside 1000 AU, and decreased molecular hydrogen content in the innermost region. Furthermore, we find that disk formation is suppressed in higher-resolution calculations, at least at the times that we can follow the calculation. We discuss the effect this may have on the buildup of disks over the accretion history of the first clump to form as well as the potential for gravitational instabilities to develop and induce fragmentation.
△ Less
Submitted 19 December, 2011;
originally announced December 2011.
-
Simulating the Common Envelope Phase of a Red Giant Using SPH and Uniform Grid Codes
Authors:
Jean-Claude Passy,
Orsola De Marco,
Chris L. Fryer,
Falk Herwig,
Steven Diehl,
Jeffrey S. Oishi,
Mordecai-Mark Mac Low,
Greg L. Bryan,
Gabriel Rockefeller
Abstract:
We use three-dimensional hydrodynamical simulations to study the rapid infall phase of the common envelope interaction of a red giant branch star of mass equal to 0.88 \msun and a companion star of mass ranging from 0.9 down to 0.1 \msun. We first compare the results obtained using two different numerical techniques with different resolutions, and find overall very good agreement. We then compare…
▽ More
We use three-dimensional hydrodynamical simulations to study the rapid infall phase of the common envelope interaction of a red giant branch star of mass equal to 0.88 \msun and a companion star of mass ranging from 0.9 down to 0.1 \msun. We first compare the results obtained using two different numerical techniques with different resolutions, and find overall very good agreement. We then compare the outcomes of those simulations with observed systems thought to have gone through a common envelope. The simulations fail to reproduce those systems in the sense that most of the envelope of the donor remains bound at the end of the simulations and the final orbital separations between the donor's remnant and the companion, ranging from 26.8 down to 5.9 \rsun, are larger than the ones observed. We suggest that this discrepancy vouches for recombination playing an essential role in the ejection of the envelope and/or significant shrinkage of the orbit happening in the subsequent phase.
△ Less
Submitted 22 September, 2011; v1 submitted 25 July, 2011;
originally announced July 2011.
-
Magnetorotational turbulence transports angular momentum in stratified disks with low magnetic Prandtl number but magnetic Reynolds number above a critical value
Authors:
Jeffrey S. Oishi,
Mordecai-Mark Mac Low
Abstract:
The magnetorotational instability (MRI) may dominate outward transport of angular momentum in accretion disks, allowing material to fall onto the central object. Previous work has established that the MRI can drive a mean-field dynamo, possibly leading to a self-sustaining accretion system. Recently, however, simulations of the scaling of the angular momentum transport parameter $\alphaSS$ with th…
▽ More
The magnetorotational instability (MRI) may dominate outward transport of angular momentum in accretion disks, allowing material to fall onto the central object. Previous work has established that the MRI can drive a mean-field dynamo, possibly leading to a self-sustaining accretion system. Recently, however, simulations of the scaling of the angular momentum transport parameter $\alphaSS$ with the magnetic Prandtl number $\Prandtl$ have cast doubt on the ability of the MRI to transport astrophysically relevant amounts of angular momentum in real disk systems. Here, we use simulations including explicit physical viscosity and resistivity to show that when vertical stratification is included, mean field dynamo action operates, driving the system to a configuration in which the magnetic field is not fully helical. This relaxes the constraints on the generated field provided by magnetic helicity conservation, allowing the generation of a mean field on timescales independent of the resistivity. Our models demonstrate the existence of a critical magnetic Reynolds number $\Rmagc$, below which transport becomes strongly $\Prandtl$-dependent and chaotic, but above which the transport is steady and $\Prandtl$-independent. Prior simulations showing $\Prandtl$-dependence had $\Rmag < \Rmagc$. We conjecture that this steady regime is possible because the mean field dynamo is not helicity-limited and thus does not depend on the details of the helicity ejection process. Scaling to realistic astrophysical parameters suggests that disks around both protostars and stellar mass black holes have $\Rmag >> \Rmagc$. Thus, we suggest that the strong $\Prandtl$ dependence seen in recent simulations does not occur in real systems.
△ Less
Submitted 12 September, 2011; v1 submitted 24 February, 2011;
originally announced February 2011.
-
A Multi-Code Analysis Toolkit for Astrophysical Simulation Data
Authors:
Matthew J. Turk,
Britton D. Smith,
Jeffrey S. Oishi,
Stephen Skory,
Samuel W. Skillman,
Tom Abel,
Michael L. Norman
Abstract:
The analysis of complex multiphysics astrophysical simulations presents a unique and rapidly growing set of challenges: reproducibility, parallelization, and vast increases in data size and complexity chief among them. In order to meet these challenges, and in order to open up new avenues for collaboration between users of multiple simulation platforms, we present yt (available at http://yt.enzoto…
▽ More
The analysis of complex multiphysics astrophysical simulations presents a unique and rapidly growing set of challenges: reproducibility, parallelization, and vast increases in data size and complexity chief among them. In order to meet these challenges, and in order to open up new avenues for collaboration between users of multiple simulation platforms, we present yt (available at http://yt.enzotools.org/), an open source, community-developed astrophysical analysis and visualization toolkit. Analysis and visualization with yt are oriented around physically relevant quantities rather than quantities native to astrophysical simulation codes. While originally designed for handling Enzo's structure adaptive mesh refinement (AMR) data, yt has been extended to work with several different simulation methods and simulation codes including Orion, RAMSES, and FLASH. We report on its methods for reading, handling, and visualizing data, including projections, multivariate volume rendering, multi-dimensional histograms, halo finding, light cone generation and topologically-connected isocontour identification. Furthermore, we discuss the underlying algorithms yt uses for processing and visualizing data, and its mechanisms for parallelization of analysis tasks.
△ Less
Submitted 15 November, 2010;
originally announced November 2010.
-
On the Stability of Dust-Laden Protoplanetary Vortices
Authors:
Philip Chang,
Jeffrey S. Oishi
Abstract:
The formation of planetesimals via gravitational instability of the dust layer in a protoplanetary disks demands that there be local patches where dust is concentrated by a factor of $\sim$ a few $\times 10^3$ over the background value. Vortices in protoplanetary disks may concentrate dust to these values allowing them to be the nurseries of planetesimals. The concentration of dust in the cores of…
▽ More
The formation of planetesimals via gravitational instability of the dust layer in a protoplanetary disks demands that there be local patches where dust is concentrated by a factor of $\sim$ a few $\times 10^3$ over the background value. Vortices in protoplanetary disks may concentrate dust to these values allowing them to be the nurseries of planetesimals. The concentration of dust in the cores of vortices increases the dust-gas ratio of the core compared to the background disk, creating a "heavy vortex." In this work, we show that these vortices are subject to an instability which we have called the heavy-core instability. Using Floquet theory, we show that this instability occurs in elliptical protoplanetary vortices when the gas-dust density of the core of the vortex is heavier than the ambient gas-dust density by a few tens of percent. The heavy-core instability grows very rapidly, with a growth timescale of a few vortex rotation periods. While the nonlinear evolution of this instability remains unknown, it will likely increase the velocity dispersion of the dust layer in the vortex because instability sets in well before sufficient dust can gather to form a protoplanetary seed. This instability may thus preclude vortices from being sites of planetesimal formation.
△ Less
Submitted 11 August, 2010; v1 submitted 14 July, 2010;
originally announced July 2010.
-
On Hydrodynamic Motions in Dead Zones
Authors:
Jeffrey S. Oishi,
Mordecai-Mark Mac Low
Abstract:
We investigate fluid motions near the midplane of vertically stratified accretion disks with highly resistive midplanes. In such disks, the magnetorotational instability drives turbulence in thin layers surrounding a resistive, stable dead zone. The turbulent layers in turn drive motions in the dead zone. We examine the properties of these motions using three-dimensional, stratified, local, shea…
▽ More
We investigate fluid motions near the midplane of vertically stratified accretion disks with highly resistive midplanes. In such disks, the magnetorotational instability drives turbulence in thin layers surrounding a resistive, stable dead zone. The turbulent layers in turn drive motions in the dead zone. We examine the properties of these motions using three-dimensional, stratified, local, shearing-box, non-ideal, magnetohydrodynamical simulations. Although the turbulence in the active zones provides a source of vorticity to the midplane, no evidence for coherent vortices is found in our simulations. It appears that this is because of strong vertical oscillations in the dead zone. By analyzing time series of azimuthally-averaged flow quantities, we identify an axisymmetric wave mode particular to models with dead zones. This mode is reduced in amplitude, but not suppressed entirely, by changing the equation of state from isothermal to ideal. These waves are too low-frequency to affect sedimentation of dust to the midplane, but may have significance for the gravitational stability of the resulting midplane dust layers.
△ Less
Submitted 2 September, 2009;
originally announced September 2009.
-
Supplementary Information for ``Rapid planetesimal formation in turbulent circumstellar discs''
Authors:
Anders Johansen,
Jeffrey S. Oishi,
Mordecai-Mark Mac Low,
Hubert Klahr,
Thomas Henning,
Andrew Youdin
Abstract:
This document contains refereed supplementary information for the paper ``Rapid planetesimal formation in turbulent circumstellar discs''. It contains 15 sections (\S1.1 -- \S1.15) that address a number of subjects related to the main paper. We describe in detail the Poisson solver used to find the self-potential of the solid particles, including a linear and a non-linear test problem (\S1.3). D…
▽ More
This document contains refereed supplementary information for the paper ``Rapid planetesimal formation in turbulent circumstellar discs''. It contains 15 sections (\S1.1 -- \S1.15) that address a number of subjects related to the main paper. We describe in detail the Poisson solver used to find the self-potential of the solid particles, including a linear and a non-linear test problem (\S1.3). Dissipative collisions remove energy from the motion of the particles by collisional cooling (\S1.4), an effect that allows gravitational collapse to occur in somewhat less massive discs (\S1.7). A resolution study of the gravitational collapse of the boulders is presented in \S1.6. We find that gravitational collapse can occur in progressively less massive discs as the grid resolution is increased, likely due to the decreased smoothing of the particle-mesh self-gravity solver with increasing resolution. In \S1.10 we show that it is in good agreement with the Goldreich & Ward (1973) stability analysis to form several-hundred-km-sized bodies, when the analysis is applied to 5 AU and to regions of increased boulder column density. \S11 is devoted to the measurement of random speeds and collision speeds between boulders. We find good agreement between our measurements and analytical theory for the random speeds, but the measured collision speeds are 3 times lower than expected from analytical theory. Higher resolution studies, and an improved analytical theory of collision speeds that takes into account epicyclic motion, will be needed to determine whether collision speeds have converged. In \S1.12 we present models with no magnetic fields. The boulder layer still exhibits strong clumping, due to the streaming instability, if the global solids-to-gas ratio is increased by a factor 3. Gravitational collapse occurs as readily as in magnetised discs.
△ Less
Submitted 29 August, 2007;
originally announced August 2007.
-
Rapid planetesimal formation in turbulent circumstellar discs
Authors:
Anders Johansen,
Jeffrey S. Oishi,
Mordecai-Mark Mac Low,
Hubert Klahr,
Thomas Henning,
Andrew Youdin
Abstract:
The initial stages of planet formation in circumstellar gas discs proceed via dust grains that collide and build up larger and larger bodies (Safronov 1969). How this process continues from metre-sized boulders to kilometre-scale planetesimals is a major unsolved problem (Dominik et al. 2007): boulders stick together poorly (Benz 2000), and spiral into the protostar in a few hundred orbits due t…
▽ More
The initial stages of planet formation in circumstellar gas discs proceed via dust grains that collide and build up larger and larger bodies (Safronov 1969). How this process continues from metre-sized boulders to kilometre-scale planetesimals is a major unsolved problem (Dominik et al. 2007): boulders stick together poorly (Benz 2000), and spiral into the protostar in a few hundred orbits due to a head wind from the slower rotating gas (Weidenschilling 1977). Gravitational collapse of the solid component has been suggested to overcome this barrier (Safronov 1969, Goldreich & Ward 1973, Youdin & Shu 2002). Even low levels of turbulence, however, inhibit sedimentation of solids to a sufficiently dense midplane layer (Weidenschilling & Cuzzi 1993, Dominik et al. 2007), but turbulence must be present to explain observed gas accretion in protostellar discs (Hartmann 1998). Here we report the discovery of efficient gravitational collapse of boulders in locally overdense regions in the midplane. The boulders concentrate initially in transient high pressures in the turbulent gas (Johansen, Klahr, & Henning 2006), and these concentrations are augmented a further order of magnitude by a streaming instability (Youdin & Goodman 2005, Johansen, Henning, & Klahr 2006, Johansen & Youdin 2007) driven by the relative flow of gas and solids. We find that gravitationally bound clusters form with masses comparable to dwarf planets and containing a distribution of boulder sizes. Gravitational collapse happens much faster than radial drift, offering a possible path to planetesimal formation in accreting circumstellar discs.
△ Less
Submitted 29 August, 2007;
originally announced August 2007.
-
Turbulent Torques on Protoplanets in a Dead Zone
Authors:
Jeffrey S. Oishi,
Mordecai-Mark Mac Low,
Kristen Menou
Abstract:
Migration of protoplanets in their gaseous host disks may be largely responsible for the observed orbital distribution of extrasolar planets. Recent simulations have shown that the magnetorotational turbulence thought to drive accretion in protoplanetary disks can affect migration by turning it into an orbital random walk. However, these simulations neglected the disk's ionization structure. Low…
▽ More
Migration of protoplanets in their gaseous host disks may be largely responsible for the observed orbital distribution of extrasolar planets. Recent simulations have shown that the magnetorotational turbulence thought to drive accretion in protoplanetary disks can affect migration by turning it into an orbital random walk. However, these simulations neglected the disk's ionization structure. Low ionization fraction near the midplane of the disk can decouple the magnetic field from the gas, forming a dead zone with reduced or no turbulence. Here, to understand the effect of dead zones on protoplanetary migration, we perform numerical simulations of a small region of a stratified disk with magnetorotational turbulence confined to thin active layers above and below the midplane. Turbulence in the active layers exerts decreased, but still measurable, gravitational torques on a protoplanet located at the disk midplane. We find a decrease of two orders of magnitude in the diffusion coefficient for dead zones with dead-to-active surface density ratios approaching realistic values in protoplanetary disks. This torque arises primarily from density fluctuations within a distance of one scale height of the protoplanet. Turbulent torques have correlation times of only $\sim 0.3$ orbital periods and apparently time-stationary distributions. These properties are encouraging signs that stochastic methods can be used to determine the orbital evolution of populations of protoplanets under turbulent migration. Our results indicate that dead zones may be dynamically distinct regions for protoplanetary migration.
△ Less
Submitted 20 July, 2007; v1 submitted 20 February, 2007;
originally announced February 2007.
-
Dynamical Expansion of H II Regions from Ultracompact to Compact Sizes in Turbulent, Self-Gravitating Molecular Clouds
Authors:
Mordecai-Mark Mac Low,
Jayashree Toraskar,
Jeffrey S. Oishi,
Tom Abel
Abstract:
The nature of ultracompact H II regions (UCHRs) remains poorly determined. In particular, they are about an order of magnitude more common than would be expected if they formed around young massive stars and lasted for one dynamical time, around 10^4 yr. We here perform three-dimensional numerical simulations of the expansion of an H II region into self-gravitating, radiatively cooled gas, both…
▽ More
The nature of ultracompact H II regions (UCHRs) remains poorly determined. In particular, they are about an order of magnitude more common than would be expected if they formed around young massive stars and lasted for one dynamical time, around 10^4 yr. We here perform three-dimensional numerical simulations of the expansion of an H II region into self-gravitating, radiatively cooled gas, both with and without supersonic turbulent flows. In the laminar case, we find that H II region expansion in a collapsing core produces nearly spherical shells, even if the ionizing source is off-center in the core. This agrees with analytic models of blast waves in power-law media. In the turbulent case, we find that the H II region does not disrupt the central collapsing region, but rather sweeps up a shell of gas in which further collapse occurs. Although this does not constitute triggering, as the swept-up gas would eventually have collapsed anyway, it does expose the collapsing regions to ionizing radiation. We suggest that these regions of secondary collapse, which will not all themselves form massive stars, may form the bulk of observed UCHRs. As the larger shell will take over 10^5 years to complete its evolution, this could solve the timescale problem. Our suggestion is supported by the ubiquitous observation of more diffuse emission surrounding UCHRs.
△ Less
Submitted 30 June, 2007; v1 submitted 19 May, 2006;
originally announced May 2006.
-
The Inability of Ambipolar Diffusion to set a Characteristic Mass Scale in Molecular Clouds
Authors:
Jeffrey S. Oishi,
Mordecai-Mark Mac Low
Abstract:
We investigate the question of whether ambipolar diffusion (ion-neutral drift) determines the smallest length and mass scale on which structure forms in a turbulent molecular cloud. We simulate magnetized turbulence in a mostly neutral, uniformly driven, turbulent medium, using a three-dimensional, two-fluid, magnetohydrodynamics (MHD) code modified from Zeus-MP. We find that substantial structu…
▽ More
We investigate the question of whether ambipolar diffusion (ion-neutral drift) determines the smallest length and mass scale on which structure forms in a turbulent molecular cloud. We simulate magnetized turbulence in a mostly neutral, uniformly driven, turbulent medium, using a three-dimensional, two-fluid, magnetohydrodynamics (MHD) code modified from Zeus-MP. We find that substantial structure persists below the ambipolar diffusion scale because of the propagation of compressive slow MHD waves at smaller scales. Contrary to simple scaling arguments, ambipolar diffusion thus does not suppress structure below its characteristic dissipation scale as would be expected for a classical diffusive process. We have found this to be true for the magnetic energy, velocity, and density. Correspondingly, ambipolar diffusion leaves the clump mass spectrum unchanged. Ambipolar diffusion appears unable to set a characteristic scale for gravitational collapse and star formation in turbulent molecular clouds.
△ Less
Submitted 12 October, 2005;
originally announced October 2005.