-
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.
-
Effective drag in rotating, poorly conducting plasma turbulence
Authors:
Santiago J. Benavides,
Keaton J. Burns,
Basile Gallet,
Glenn R. Flierl
Abstract:
Despite the increasing sophistication of numerical models of hot Jupiter atmospheres, the large time-scale separation required in simulating the wide range in electrical conductivity between the dayside and nightside has made it difficult to run fully consistent magnetohydrodynamic (MHD) models. This has led to many studies that resort to drag parametrizations of MHD. In this study, we revisit the…
▽ More
Despite the increasing sophistication of numerical models of hot Jupiter atmospheres, the large time-scale separation required in simulating the wide range in electrical conductivity between the dayside and nightside has made it difficult to run fully consistent magnetohydrodynamic (MHD) models. This has led to many studies that resort to drag parametrizations of MHD. In this study, we revisit the question of the Lorentz force as an effective drag by running a series of direct numerical simulations of a weakly rotating, poorly conducting flow in the presence of a misaligned, strong background magnetic field. We find that the drag parametrization fails once the time-scale associated with the Lorentz force becomes shorter than the dynamical time-scale in the system, beyond which the effective drag coefficient remains roughly constant, despite orders-of-magnitude variation in the Lorentz (magnetic) time-scale. We offer an improvement to the drag parametrization by considering the relevant asymptotic limit of low conductivity and strong background magnetic field, known as the quasi-static MHD approximation of the Lorentz force. This approximation removes the fast time-scale associated with magnetic diffusion, but retains a more complex version of the Lorentz force, which could be utilized in future numerical models of hot Jupiter atmospheric circulation.
△ Less
Submitted 20 October, 2022; v1 submitted 9 March, 2022;
originally announced March 2022.
-
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.
-
Linearly forced fluid flow on a rotating sphere
Authors:
Rohit Supekar,
Vili Heinonen,
Keaton J. Burns,
Jörn Dunkel
Abstract:
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 so…
▽ More
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.
△ Less
Submitted 9 April, 2020; v1 submitted 4 December, 2019;
originally announced December 2019.
-
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.
-
Low-frequency variability in massive stars: Core generation or surface phenomenon?
Authors:
Daniel Lecoanet,
Matteo Cantiello,
Eliot Quataert,
Louis-Alexandre Couston,
Keaton J. Burns,
Benjamin J. S. Pope,
Adam S. Jermyn,
Benjamin Favier,
Michael Le Bars
Abstract:
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, det…
▽ More
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.
△ Less
Submitted 12 November, 2019; v1 submitted 3 October, 2019;
originally announced October 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.
-
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.
-
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.
-
Conversion of Internal Gravity Waves into Magnetic Waves
Authors:
Daniel Lecoanet,
Geoffrey M Vasil,
Jim Fuller,
Matteo Cantiello,
Keaton J Burns
Abstract:
Asteroseismology probes the interiors of stars by studying oscillation modes at a star's surface. Although pulsation spectra are well understood for solar-like oscillators, a substantial fraction of red giant stars observed by Kepler exhibit abnormally low-amplitude dipole oscillation modes. Fuller et al. (2015) suggests this effect is produced by strong core magnetic fields that scatter dipole in…
▽ More
Asteroseismology probes the interiors of stars by studying oscillation modes at a star's surface. Although pulsation spectra are well understood for solar-like oscillators, a substantial fraction of red giant stars observed by Kepler exhibit abnormally low-amplitude dipole oscillation modes. Fuller et al. (2015) suggests this effect is produced by strong core magnetic fields that scatter dipole internal gravity waves (IGWs) into higher multipole IGWs or magnetic waves. In this paper, we study the interaction of IGWs with a magnetic field to test this mechanism. We consider two background stellar structures: one with a uniform magnetic field, and another with a magnetic field that varies both horizontally and vertically. We derive analytic solutions to the wave propagation problem and validate them with numerical simulations. In both cases, we find perfect conversion from IGWs into magnetic waves when the IGWs propagate into a region exceeding a critical magnetic field strength. Downward propagating IGWs cannot reflect into upward propagating IGWs because their vertical wavenumber never approaches zero. Instead, they are converted into upward propagating slow (Alfvenic) waves, and we show they will likely dissipate as they propagate back into weakly magnetized regions. Therefore, strong internal magnetic fields can produce dipole mode suppression in red giants, and gravity modes will likely be totally absent from the pulsation spectra of sufficiently magnetized stars.
△ Less
Submitted 26 December, 2016; v1 submitted 26 October, 2016;
originally announced October 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.
-
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.
-
FIRST, a fibered aperture masking instrument. I. First on-sky test results
Authors:
E. Huby,
G. Perrin,
F. Marchis,
S. Lacour,
T. Kotani,
G. Duchêne,
E. Choquet,
E. L. Gates,
J. M. Woillez,
O. Lai,
P. Fédou,
C. Collin,
F. Chapron,
V. Arslanyan,
K. J. Burns
Abstract:
In this paper we present the first on-sky results with the fibered aperture masking instrument FIRST. Its principle relies on the combination of spatial filtering and aperture masking using single-mode fibers, a novel technique that is aimed at high dynamic range imaging with high angular resolution. The prototype has been tested with the Shane 3-m telescope at Lick Observatory. The entrance pupil…
▽ More
In this paper we present the first on-sky results with the fibered aperture masking instrument FIRST. Its principle relies on the combination of spatial filtering and aperture masking using single-mode fibers, a novel technique that is aimed at high dynamic range imaging with high angular resolution. The prototype has been tested with the Shane 3-m telescope at Lick Observatory. The entrance pupil is divided into subpupils feeding single-mode fibers. The flux injection into the fibers is optimized by a segmented mirror. The beams are spectrally dispersed and recombined in a non-redundant exit configuration in order to retrieve all contrasts and phases independently. The instrument works at visible wavelengths between 600 nm and 760 nm and currently uses nine of the 30 43 cm subapertures constituting the full pupil. First fringes were obtained on Vega and Deneb. Stable closure phases were measured with standard deviations on the order of 1 degree. Closure phase precision can be further improved by addressing some of the remaining sources of systematic errors. While the number of fibers used in the experiment was too small to reliably estimate visibility amplitudes, we have measured closure amplitudes with a precision of 10 % in the best case. These first promising results obtained under real observing conditions validate the concept of the fibered aperture masking instrument and open the way for a new type of ground-based instrument working in the visible. The next steps of the development will be to improve the stability and the sensitivity of the instrument in order to achieve more accurate closure phase and visibility measurements, and to increase the number of sub-pupils to reach full pupil coverage.
△ Less
Submitted 22 March, 2012;
originally announced March 2012.