Magnetic fields are observed in massive Ap/Bp stars and are presumably present in the radiative zone of solar-like stars. To date, there is no clear understanding of the dynamics of the magnetic field in stably stratified layers. A purely toroidal magnetic field configuration is known to be unstable, developing mainly non-axisymmetric modes. Rotation and a poloidal field component may lead to stabilization. Here we perform global MHD simulations with the EULAG-MHD code to explore the evolution of a toroidal magnetic field located in a layer whose Brunt-Väisälä frequency resembles the lower solar tachocline. Our numerical experiments allow us to explore the initial unstable phase as well as the long-term evolution of such field. During the first Alfven cycles, we observe the development of the Tayler instability with the prominent longitudinal wavenumber, m = 1. Rotation decreases the growth rate of the instability and eventually suppresses it. However, after a stable phase, energy surges lead to the development of higher-order modes even for fast rotation. These modes extract energy from the initial toroidal field. Nevertheless, our results show that sufficiently fast rotation leads to a lower saturation energy of the unstable modes, resulting in a magnetic topology with only a small fraction of poloidal field, which remains steady for several hundreds of Alfven traveltimes. The system then becomes turbulent and the field is prone to turbulent diffusion. The final toroidal–poloidal configuration of the magnetic field may represent an important aspect of the field generation and evolution in stably stratified layers.
We study the properties of oscillatory double-diffusive convection (ODDC) in the presence of a uniform vertical background magnetic field. ODDC takes place in stellar regions that are unstable according to the Schwarzschild criterion and stable according to the Ledoux criterion (sometimes called semiconvective regions), which are often predicted to reside just outside the core of intermediate-mass main-sequence stars. Previous hydrodynamic studies of ODDC have shown that the basic instability saturates into a state of weak wave-like convection, but that a secondary instability can sometimes transform it into a state of layered convection, where layers then rapidly merge and grow until the entire region is fully convective. We find that magnetized ODDC has very similar properties overall, with some important quantitative differences. A linear stability analysis reveals that the fastest-growing modes are unaffected by the field, but that other modes are. Numerically, the magnetic field is seen to influence the saturation of the basic instability, overall reducing the turbulent fluxes of temperature and composition. This in turn affects layer formation, usually delaying it, and occasionally suppressing it entirely for sufficiently strong fields. Further work will be needed, however, to determine the field strength above which layer formation is actually suppressed in stars. Potential observational implications are briefly discussed.
more » « less- Award ID(s):
- 1908338
- NSF-PAR ID:
- 10369396
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 935
- Issue:
- 1
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 33
- Size(s):
- ["Article No. 33"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
more » « less
Abstract Understanding Jupiter's present‐day interior structure and dynamics is key to constraining planetary accretion models. In particular, the extent of stable stratification (i.e., non‐convective regions) in the planet strongly influences long‐term cooling processes, and may record primordial heavy element gradients from early in a planet's formation. Because the Galileo entry probe measured a subsolar helium abundance, Jupiter interior models often invoke an outer stably stratified region due to helium rain. Additionally, Juno gravity data suggest a deeper, potentially stratified dilute core extending halfway through the planet. However, fits to Jupiter's gravitational data are non‐unique, and outstanding uncertainty over the equations of state for hydrogen and helium remain. Here, we use high‐resolution numerical magnetohydrodynamic simulations of Jupiter's magnetic field to place constraints on the extent of stable stratification within the planet. We find that compared to traditional interior models, an upper stably stratified layer between 0.9 and 0.95 Jupiter radii (
R J ) helps to explain both Jupiter's dipolar magnetic field and zonal winds. In contrast, an extended dilute core that is entirely stably stratified (no convective layers) yields significantly worse fits to both. However, our models with extended deep stratification still generate dipolar magnetic fields if an upper stratified region is also present. Overall, we find that a planet with a dilute core i.e., strongly stably stratified is increasingly challenging to reconcile with Jupiter's magnetic field and winds. Thus if a dilute core is present, alternative modalities such as a fully convective dilute core, a complex multilayered interior structure, or double diffusive convection may be required. -
null (Ed.)A theoretical model is developed which illustrates the dynamics of layering instability, frequently realized in ocean regions with active fingering convection. Thermohaline layering is driven by the interplay between large-scale stratification and primary double-diffusive instabilities operating at the microscale – temporal and spatial scales set by molecular dissipation. This interaction is described by a combination of direct numerical simulations and an asymptotic multiscale model. The multiscale theory is used to formulate explicit and dynamically consistent flux laws, which can be readily implemented in large-scale analytical and numerical models. Most previous theoretical investigations of thermohaline layering were based on the flux-gradient model, which assumes that the vertical transport of density components is uniquely determined by their local background gradients. The key deficiency of this approach is that layering instabilities predicted by the flux-gradient model have unbounded growth rates at high wavenumbers. The resulting ultraviolet catastrophe precludes the analysis of such basic properties of layering instability as its preferred wavelength or the maximal growth rate. The multiscale model, on the other hand, incorporates hyperdiffusion terms that stabilize short layering modes. Overall, the presented theory carries the triple advantage of (i) offering an explicit description of the interaction between microstructure and layering modes, (ii) taking into account the influence of non-uniform stratification on microstructure-driven mixing, and (iii) avoiding unphysical behaviour of the flux-gradient laws at small scales. While the multiscale approach to the parametrization of time-dependent small-scale processes is illustrated here on the example of fingering convection, we expect the proposed technique to be readily adaptable to a wide range of applications.more » « less
-
more » « less
Abstract Surface waves on Earth's magnetopause have a controlling effect upon global magnetospheric dynamics. Since spacecraft provide sparse in situ observation points, remote sensing these modes using ground‐based instruments in the polar regions is desirable. However, many open conceptual questions on the expected signatures remain. Therefore, we provide predictions of key qualitative features expected in auroral, ionospheric, and ground magnetic observations through both magnetohydrodynamic theory and a global coupled magnetosphere‐ionosphere simulation of a magnetopause surface eigenmode. These show monochromatic oscillatory field‐aligned currents (FACs), due to both the surface mode and its non‐resonant Alfvén coupling, are present throughout the magnetosphere. The currents peak in amplitude at the equatorward edge of the magnetopause boundary layer, not the open‐closed boundary as previously thought. They also exhibit slow poleward phase motion rather than being purely evanescent. We suggest the upward FAC perturbations may result in periodic auroral brightenings. In the ionosphere, convection vortices circulate the poleward moving FAC structures. Finally, surface mode signals are predicted in the ground magnetic field, with ionospheric Hall currents rotating perturbations by approximately (but not exactly) 90° compared to the magnetosphere. Thus typical dayside magnetopause surface modes should be strongest in the East‐West ground magnetic field component. Overall, all ground‐based signatures of the magnetopause surface mode are predicted to have the same frequency across
L ‐shells, amplitudes that maximize near the magnetopause's equatorward edge, and larger latitudinal scales than for field line resonance. Implications in terms of ionospheric Joule heating and geomagnetically induced currents are discussed. -
more » « less
Abstract Both the core collapse of rotating massive stars, and the coalescence of neutron star (NS) binaries result in the formation of a hot, differentially rotating NS remnant. The timescales over which differential rotation is removed by internal angular-momentum transport processes (
viscosity ) have key implications for the remnant’s long-term stability and the NS equation of state (EOS). Guided by a nonrotating model of a cooling proto-NS, we estimate the dominant sources of viscosity using an externally imposed angular-velocity profile Ω(r ). Although the magneto-rotational instability provides the dominant source of effective viscosity at large radii, convection and/or the Tayler–Spruit dynamo dominate in the core of merger remnants whered Ω/dr ≥ 0. Furthermore, the viscous timescale in the remnant core is sufficiently short that solid-body rotation will be enforced faster than matter is accreted from rotationally supported outer layers. Guided by these results, we develop a toy model for how the merger remnant core grows in mass and angular momentum due to accretion. We find that merger remnants with sufficiently massive and slowly rotating initial cores may collapse to black holes via envelope accretion, even when the total remnant mass is less than the usually considered threshold ≈1.2M TOVfor forming a stable solid-body rotating NS remnant (whereM TOVis the maximum nonrotating NS mass supported by the EOS). This qualitatively new picture of the post-merger remnant evolution and stability criterion has important implications for the expected electromagnetic counterparts from binary NS mergers and for multimessenger constraints on the NS EOS.
