BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Adrian J. Barker

Adrian J. Barker appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.EPphysics.flu-dynastro-ph.HEphysics.ao-ph

Catalog footprint

What is connected

19works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

19 published item(s)

preprint2022arXiv

Tidal dissipation due to inertial waves can explain the circularization periods of solar-type binaries

Tidal dissipation is responsible for circularizing the orbits and synchronizing the spins of solar-type close binary stars, but the mechanisms responsible are not fully understood. Previous work has indicated that significant enhancements to the theoretically-predicted tidal dissipation rates are required to explain the observed circularization periods ($P_\mathrm{circ}$) in various stellar populations, and their evolution with age. This was based partly on the common belief that the dominant mechanism of tidal dissipation in solar-type stars is turbulent viscosity acting on equilibrium tides in convective envelopes. In this paper we study tidal dissipation in both convection and radiation zones of rotating solar-type stars following their evolution. We study equilibrium tide dissipation, incorporating a frequency-dependent effective viscosity motivated by the latest hydrodynamical simulations, and inertial wave (dynamical tide) dissipation, adopting a frequency-averaged formalism that accounts for the realistic structure of the star. We demonstrate that the observed binary circularization periods can be explained by inertial wave (dynamical tide) dissipation in convective envelopes. This mechanism is particularly efficient during pre-main sequence phases, but it also operates on the main sequence if the spin is close to synchronism. The predicted $P_\mathrm{circ}$ due to this mechanism increases with main-sequence age in accord with observations. We also demonstrate that both equilibrium tide and internal gravity wave dissipation are unlikely to explain the observed $P_\mathrm{circ}$, even during the pre-main sequence, based on our best current understanding of these mechanisms. Finally, we advocate more realistic dynamical studies of stellar populations that employ tidal dissipation due to inertial waves.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Angular momentum transport, layering, and zonal jet formation by the GSF instability: nonlinear simulations at a general latitude

We continue our investigation into the nonlinear evolution of the Goldreich-Schubert-Fricke (GSF) instability in differentially rotating radiation zones. This instability may be a key player in transporting angular momentum in stars and giant planets, but its nonlinear evolution remains mostly unexplored. In a previous paper we considered the equatorial instability, whereas here we simulate the instability at a general latitude for the first time. We adopt a local Cartesian Boussinesq model in a modified shearing box for most of our simulations, but we also perform some simulations with stress-free, impenetrable, radial boundaries. We first revisit the linear instability and derive some new results, before studying its nonlinear evolution. The instability is found to behave very differently compared with its behaviour at the equator. In particular, here we observe the development of strong zonal jets ("layering" in the angular momentum), which can considerably enhance angular momentum transport, particularly in axisymmetric simulations. The jets are, in general, tilted with respect to the local gravity by an angle that corresponds initially with that of the linear modes, but which evolves with time and depends on the strength of the flow. The instability transports angular momentum much more efficiently (by several orders of magnitude) than it does at the equator, and we estimate that the GSF instability could contribute to the missing angular momentum transport required in both red giant and subgiant stars. It could also play a role in the long-term evolution of the solar tachocline and the atmospheric dynamics of hot Jupiters.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Convection with Misaligned Gravity and Rotation: Simulations and Rotating Mixing Length Theory

We present numerical simulations, using two complementary setups, of rotating Boussinesq thermal convection in a three-dimensional Cartesian geometry with misaligned gravity and rotation vectors. This model represents a small region at a non-polar latitude in the convection zone of a star or planet. We investigate the effects of rotation on the bulk properties of convection at different latitudes, focusing on determining the relation between the heat flux and temperature gradient. We show that our results may be interpreted using rotating mixing length theory (RMLT). The simplest version of RMLT (due to Stevenson) considers the single mode that transports the most heat. This works reasonably well in explaining our results, but there is a systematic departure from these predictions (up to approximately $30\%$ in the temperature gradient) at mid-latitudes. We develop a more detailed treatment of RMLT that includes the transport afforded by multiple modes, and we show that this accounts for most of the systematic differences. We also show that convectively-generated zonal flows and meridional circulations are produced in our simulations, and that their properties depend strongly on the dimensions of the box. These flows also affect the heat transport, contributing to departures from RMLT at some latitudes. However, we find the theoretical predictions of the multi-mode theory for the mid-layer temperature gradient, the root-mean-square (RMS) vertical velocity, the RMS temperature fluctuation, and the spatial spectrum of the heat transport at different latitudes, are all in reasonably good agreement with our numerical results when zonal flows are small.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Convective turbulent viscosity acting on equilibrium tidal flows: new frequency scaling of the effective viscosity

Turbulent convection is thought to act as an effective viscosity ($ν_E$) in damping tidal flows in stars and giant planets. However, the efficiency of this mechanism has long been debated, particularly in the regime of fast tides, when the tidal frequency ($ω$) exceeds the turnover frequency of the dominant convective eddies ($ω_c$). We present the results of hydrodynamical simulations to study the interaction between tidal flows and convection in a small patch of a convection zone. These simulations build upon our prior work by simulating more turbulent convection in larger horizontal boxes, and here we explore a wider range of parameters. We obtain several new results: 1) $ν_E$ is frequency-dependent, scaling as $ω^{-0.5}$ when $ω/ω_c \lesssim 1$, and appears to attain its maximum constant value only for very small frequencies ($ω/ω_c \lesssim 10^{-2}$). This frequency-reduction for low frequency tidal forcing has never been observed previously. 2) The frequency-dependence of $ν_E$ appears to follow the same scaling as the frequency spectrum of the energy (or Reynolds stress) for low and intermediate frequencies. 3) For high frequencies ($ω/ω_c\gtrsim 1-5$), $ν_E\propto ω^{-2}$. 4) The energetically-dominant convective modes always appear to contribute the most to $ν_E$, rather than the resonant eddies in a Kolmogorov cascade. These results have important implications for tidal dissipation in convection zones of stars and planets, and indicate that the classical tidal theory of the equilibrium tide in stars and giant planets should be revisited. We briefly touch upon the implications for planetary orbital decay around evolving stars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficiency of tidal dissipation in slowly rotating fully convective stars or planets

Turbulent convection is thought to act as an effective viscosity in damping equilibrium tidal flows, driving spin and orbital evolution in close convective binary systems. Compared to mixing-length predictions, this viscosity ought to be reduced when the tidal frequency $|ω_t|$ exceeds the turnover frequency $ω_{cν}$ of the dominant convective eddies, but the efficiency of this reduction has been disputed. We reexamine this long-standing controversy using direct numerical simulations of an idealized global model. We simulate thermal convection in a full sphere, and externally forced by the equilibrium tidal flow, to measure the effective viscosity $ν_E$ acting on the tidal flow when $|ω_t|/ω_{cν} \gtrsim 1$. We demonstrate that the frequency reduction of $ν_E$ is correlated with the frequency spectrum of the (unperturbed) convection. For intermediate frequencies below those in the turbulent cascade ($|ω_t|/ω_{cν} \sim 1-5$), the frequency spectrum displays an anomalous $1/ω^α$ power law that is responsible for the frequency-reduction $ν_E \propto 1/|ω_t|^α$, where $α< 1$ depends on the model parameters. We then get $|ν_E| \propto 1/|ω_t|^δ$ with $δ> 1$ for higher frequencies, and $δ=2$ is obtained for a Kolmogorov turbulent cascade. A generic $|ν_E| \propto 1/|ω_t|^{2}$ suppression is next found for higher frequencies within the dissipation range of the convection (but with negative values). Our results indicate that a better knowledge of the frequency spectrum of convection is necessary to accurately predict the efficiency of tidal dissipation in stars and planets resulting from this mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tidal dissipation in evolving low-mass and solar-type stars with predictions for planetary orbital decay

We study tidal dissipation in stars with masses in the range $0.1-1.6 M_\odot$ throughout their evolution, including turbulent effective viscosity acting on equilibrium tides and inertial waves in convection zones, and internal gravity waves in radiation zones. We consider a range of stellar evolutionary models and incorporate the frequency-dependent effective viscosity acting on equilibrium tides based on the latest simulations. We compare the tidal flow and dissipation obtained with the conventional equilibrium tide, which is strictly invalid in convection zones, finding that the latter typically over-predicts the dissipation by a factor of 2-3. Dissipation of inertial waves is computed using a frequency-averaged formalism accounting for realistic stellar structure for the first time, and is the dominant mechanism for binary circularization and synchronization on the main sequence. Dissipation of gravity waves in the radiation zone assumes these waves to be fully damped (e.g.~by wave breaking), and is the dominant mechanism for planetary orbital decay. We calculate the critical planetary mass required for wave breaking as a function of stellar mass and age, and show that this mechanism predicts destruction of many hot Jupiters but probably not Earth-mass planets on the main sequence. We apply our results to compute tidal quality factors following stellar evolution, and tidal evolutionary timescales, for the orbital decay of hot Jupiters, and the spin synchronization and circularization of binary stars. We also provide predictions for shifts in transit arrival times due to tidally-driven orbital decay of hot Jupiters that may be detected with NGTS, TESS or PLATO.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Turbulent viscosity acting on the equilibrium tidal flow in convective stars

Convection is thought to act as a turbulent viscosity in damping tidal flows and in driving spin and orbital evolution in close convective binary systems. This turbulent viscosity should be reduced, compared to mixing-length predictions, when the forcing (tidal) frequency $|ω_t|$ exceeds the turnover frequency $ω_{cv}$ of the dominant convective eddies. However, two contradictory scaling laws have been proposed and this issue remains highly disputed. To revisit this controversy, we conduct the first direct numerical simulations (DNS) of convection interacting with the equilibrium tidal flow in an idealized global model of a low-mass star. We present direct computations of the turbulent effective viscosity, $ν_E$, acting on the equilibrium tidal flow. We unexpectedly report the coexistence of the two disputed scaling laws, which reconciles previous theoretical (and numerical) findings. We recover the universal quadratic scaling $ν_E \propto (|ω_t|/ω_{cv})^{-2}$ in the high-frequency regime $|ω_t|/ω_{cv} \gg 1$. Our results also support the linear scaling $ν_E \propto (|ω_t|/ω_{cv})^{-1}$ in an intermediate regime with $1 \leq |ω_t|/ω_{cv} = \mathcal{O}(10)$. Both regimes may be relevant to explain the observed properties of close binaries, including spin synchronization of solar-type stars and the circularization of low-mass stars. The robustness of these two regimes of tidal dissipation, and the transition between them, should be explored further in more realistic models. A better understanding of the interaction between convection and tidal flows is indeed essential to correctly interpret observations of close binary stars and short-period planetary orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Wave propagation in semi-convective regions of giant planets

Recent observations of Jupiter and Saturn suggest that heavy elements may be diluted in the gaseous envelope, providing a compositional gradient that could stabilise ordinary convection and produce a stably-stratified layer near the core of these planets. This region could consist of semi-convective layers with a staircase-like density profile, which have multiple convective zones separated by thin stably-stratified interfaces, as a result of double-diffusive convection. These layers could have important effects on wave propagation and tidal dissipation that have not been fully explored. We analyse the effects of these layers on the propagation and transmission of internal waves within giant planets, extending prior work in a local Cartesian model. We adopt a simplified global Boussinesq planetary model in which we explore the internal waves in a non-rotating spherical body. We begin by studying the free modes of a region containing semi-convective layers. We then analyse the transmission of internal waves through such a region. The free modes depend strongly on the staircase properties, and consist of modes with both internal and interfacial gravity wave-like behaviour. We determine the frequency shifts of these waves as a function of the number of steps to explore their potential to probe planetary internal structures. We also find that wave transmission is strongly affected by the presence of a staircase. Very large-wavelength waves are transmitted efficiently, but small-scale waves are only transmitted if they are resonant with one of the free modes. The effective size of the core is therefore larger for non-resonant modes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nonlinear hydrodynamical evolution of eccentric Keplerian discs in two dimensions: validation of secular theory

We perform global two-dimensional hydrodynamical simulations of Keplerian discs with free eccentricity over thousands of orbital periods. Our aim is to determine the validity of secular theory in describing the evolution of eccentric discs, and to explore their nonlinear evolution for moderate eccentricities. Linear secular theory is found to correctly predict the structure and precession rates of discs with small eccentricities. However, discs with larger eccentricities (and eccentricity gradients) are observed to precess faster (retrograde relative to the orbital motion), at a rate that depends on their eccentricities (and eccentricity gradients). We derive analytically a nonlinear secular theory for eccentric gas discs, which explains this result as a modification of the pressure forces whenever eccentric orbits in a disc nearly intersect. This effect could be particularly important for highly eccentric discs produced in tidal disruption events, or for narrow gaseous rings; it might also play a role in causing some of the variability in superhump binary systems. In two dimensions, the eccentricity of a moderately eccentric disc is long-lived and persists throughout the duration of our simulations. Eccentric modes are however weakly damped by their interaction with non-axisymmetric spiral density waves (driven by the Papaloizou-Pringle instability, which occurs in our idealised setup with solid walls), as well as numerical diffusion.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nonlinear tides in a homogeneous rotating planet or star: global modes and elliptical instability

We revisit the global modes and instabilities of homogeneous rotating ellipsoidal fluid masses, which are the simplest global models of rotationally and tidally deformed gaseous planets or stars. The tidal flow in a short-period planet may be unstable to the elliptical instability, a hydrodynamic instability that can drive tidal evolution. We perform a global (and local WKB) analysis to study this instability using the elegant formalism of Lebovitz & Lifschitz. We survey the parameter space of global instabilities with harmonic orders $\ell\leq 5$, for planets with spins that are purely aligned (prograde) or anti-aligned (retrograde) with their orbits. In general, the instability has a much larger growth rate if the planetary spin and orbit are anti-aligned rather than aligned. We have identified a violent instability for anti-aligned spins outside of the usual frequency range for the elliptical instability (when $\frac{n}Ω\lesssim -1$, where $n$ and $Ω$ are the orbital and spin angular frequencies, respectively) if the tidal amplitude is sufficiently large. We also explore the instability in a rigid ellipsoidal container, which is found to be quantitatively similar to that with a realistic free surface. Finally, we study the effect of rotation and tidal deformation on mode frequencies. We find that larger rotation rates and larger tidal deformations both decrease the frequencies of the prograde sectoral surface gravity modes. This increases the prospect of their tidal excitation, potentially enhancing the tidal response over expectations from linear theory. In a companion paper, we use our results to interpret global simulations of the elliptical instability.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nonlinear tides in a homogeneous rotating planet or star: global simulations of the elliptical instability

I present results from the first global hydrodynamical simulations of the elliptical instability in a tidally deformed gaseous planet (or star) with a free surface. The elliptical instability is potentially important for tidal evolution of the shortest-period hot Jupiters. I model the planet as a spin-orbit aligned or anti-aligned, and non-synchronously rotating, tidally deformed, homogeneous fluid body. A companion paper presented an analysis of the global modes and instabilities of such a planet. Here I focus on the nonlinear evolution of the elliptical instability. This is observed to produce bursts of turbulence that drive the planet towards synchronism with its orbit in an erratic manner. If the planetary spin is initially anti-aligned, the elliptical instability also drives spin-orbit alignment on a similar timescale as the spin synchronisation. The instability generates differential rotation inside the planet in the form of zonal flows, which play an important role in the saturation of the instability, and in producing the observed burstiness. These results are broadly consistent with the picture obtained using a local Cartesian model (where columnar vortices played the role of zonal flows). I also simulate the instability in a container that is rigid (but stress-free) rather than free, finding broad quantitative agreement. The dissipation resulting from the elliptical instability could explain why the shortest-period hot Jupiters tend to have circular orbits inside about 2-3 days, and predicts spin-synchronisation (and spin-orbit alignment) out to about 10-15 days. However, other mechanisms must be invoked to explain tidal circularisation for longer orbital periods.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On turbulence driven by axial precession and tidal evolution of the spin-orbit angle of close-in giant planets

The spin axis of a rotationally deformed planet is forced to precess about its orbital angular momentum vector, due to the tidal gravity of its host star, if these directions are misaligned. This induces internal fluid motions inside the planet that are subject to a hydrodynamic instability. We study the turbulent damping of precessional fluid motions, as a result of this instability, in the simplest local computational model of a giant planet (or star), with and without a weak internal magnetic field. Our aim is to determine the outcome of this instability, and its importance in driving tidal evolution of the spin-orbit angle in precessing planets (and stars). We find that this instability produces turbulent dissipation that is sufficiently strong that it could drive significant tidal evolution of the spin-orbit angle for hot Jupiters with orbital periods shorter than about 10-18 days. If this mechanism acts in isolation, this evolution would be towards alignment or anti-alignment, depending on the initial angle, but the ultimate evolution (if other tidal mechanisms also contribute) is expected to be towards alignment. The turbulent dissipation is proportional to the cube of the precession frequency, so it leads to much slower damping of stellar spin-orbit angles, implying that this instability is unlikely to drive evolution of the spin-orbit angle in stars (either in planetary or close binary systems). We also find that the instability-driven flow can act as a system-scale dynamo, which may play a role in producing magnetic fields in short-period planets.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the vertical-shear instability in astrophysical discs

We explore the linear stability of astrophysical discs exhibiting vertical shear, which arises when there is a radial variation in the temperature or entropy. Such discs are subject to a "vertical-shear instability", which recent nonlinear simulations have shown to drive hydrodynamic activity in the MRI-stable regions of protoplanetary discs. We first revisit locally isothermal discs using the quasi-global reduced model derived by Nelson et al. (2013). This analysis is then extended to global axisymmetric perturbations in a cylindrical domain. We also derive and study a reduced model describing discs with power law radial entropy profiles ("locally polytropic discs"), which are somewhat more realistic in that they possess physical (as opposed to numerical) surfaces. In all cases the fastest growing modes have very short wavelengths and are localised at the disc surfaces (if present), where the vertical shear is maximal. An additional class of modestly growing vertically global body modes is excited, corresponding to destabilised classical inertial waves ("r-modes"). We discuss the properties of both types of modes, and stress that those that grow fastest occur on the shortest available length scales (determined either by the numerical grid or the physical viscous length). This ill-posedness makes simulations of the instability difficult to interpret. We end with some brief speculation on the nonlinear saturation and resulting angular momentum transport.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hydrodynamic instability in eccentric astrophysical discs

Eccentric Keplerian discs are believed to be unstable to three-dimensional hydrodynamical instabilities driven by the time-dependence of fluid properties around an orbit. These instabilities could lead to small-scale turbulence, and ultimately modify the global disc properties. We use a local model of an eccentric disc, derived in a companion paper, to compute the nonlinear vertical ("breathing mode") oscillations of the disc. We then analyse their linear stability to locally axisymmetric disturbances for any disc eccentricity and eccentricity gradient using a numerical Floquet method. In the limit of small departures from a circular reference orbit, the instability of an isothermal disc is explained analytically. We also study analytically the small-scale instability of an eccentric neutrally stratified polytropic disc with any polytropic index using a WKB approximation. We find that eccentric discs are generically unstable to the parametric excitation of small-scale inertial waves. The nonlinear evolution of these instabilities should be studied in numerical simulations, where we expect them to lead to a decay of the disc eccentricity and eccentricity gradient as well as to induce additional transport and mixing. Our results highlight that it is essential to consider the three-dimensional structure of eccentric discs, and their resulting vertical oscillatory flows, in order to correctly capture their evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Local and global dynamics of eccentric astrophysical discs

We formulate a local dynamical model of an eccentric disc in which the dominant motion consists of elliptical Keplerian orbits. The model is a generalization of the well known shearing sheet, and is suitable for both analytical and computational studies of the local dynamics of eccentric discs. It is spatially homogeneous in the horizontal dimensions but has a time-dependent geometry that oscillates at the orbital frequency. We show how certain averages of the stress tensor in the local model determine the large-scale evolution of the shape and mass distribution of the disc. The simplest solutions of the local model are laminar flows consisting of a (generally nonlinear) vertical oscillation of the disc. Eccentric discs lack vertical hydrostatic equilibrium because of the variation of the vertical gravitational acceleration around the eccentric orbit, and in some cases because of the divergence of the orbital velocity field associated with an eccentricity gradient. We discuss the properties of the laminar solutions, showing that they can exhibit extreme compressional behaviour for eccentricities greater than about $0.5$, especially in discs that behave isothermally. We also derive the linear evolutionary equations for an eccentric disc that follow from the laminar flows in the absence of a shear viscosity. In a companion paper we show that these solutions are linearly unstable and we determine the associated growth rates and unstable modes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Theory and simulations of rotating convection

We study thermal convection in a rotating fluid in order to better understand the properties of convection zones in rotating stars and planets. We first derive mixing-length theory for rapidly-rotating convection, arriving at the results of Stevenson (1979) via simple physical arguments. The theory predicts the properties of convection as a function of the imposed heat flux and rotation rate, independent of microscopic diffusivities. In particular, it predicts the mean temperature gradient; the rms velocity and temperature fluctuations; and the size of the eddies that dominate heat transport. We test all of these predictions with high resolution three-dimensional hydrodynamical simulations of Boussinesq convection in a Cartesian box. The results agree remarkably well with the theory across more than two orders of magnitude in rotation rate. For example, the temperature gradient is predicted to scale as the rotation rate to the 4/5th power at fixed flux, and the simulations yield $0.75\pm 0.06$. We conclude that the mixing length theory is a solid foundation for understanding the properties of convection zones in rotating stars and planets.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonlinear evolution of the elliptical instability in the presence of weak magnetic fields

We investigate whether the elliptical instability is important for tidal dissipation in gaseous planets and stars. In a companion paper, we found that the conventional elliptical instability results in insufficient dissipation because it produces long-lived vortices that then quench further instability. Here, we study whether the addition of a magnetic field prevents those vortices from forming, and hence leads to enhanced dissipation. We present results from magnetohydrodynamical simulations that evolve the elliptical instability in a local patch of a rotating planet or star, in the presence of a weak magnetic field. We find that magnetic fields do indeed prevent vortices from forming, and hence greatly enhance the steady state dissipation rate. In addition, the resulting turbulence acts as a small-scale dynamo, amplifying the initially weak field. The inferred tidal dissipation is potentially important at short orbital periods. For example, it can circularise hot Jupiters with orbital periods shorter than 2.5 days, and synchronise their spins with their orbits out to 6 days. However, it appears unable to account for the hot Jupiters that appear to have been circularised out to six to ten day orbital periods. It also cannot account for the inferred circularisation of many close binary stars.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonlinear evolution of the tidal elliptical instability in gaseous planets and stars

Tidally distorted rotating stars and gaseous planets are subject to a well-known linear fluid instability -- the elliptical instability. It has been proposed that this instability might drive enough energy dissipation to solve the long-standing problem of the origin of tidal dissipation in stars and planets. But the nonlinear outcome of the elliptical instability has yet to be investigated in the parameter regime of interest, and the resulting turbulent energy dissipation has not yet been quantified. We do so by performing three dimensional hydrodynamical simulations of a small patch of a tidally deformed fluid planet or star subject to the elliptical instability. We show that when the tidal deformation is weak, the nonlinear outcome of the instability leads to the formation of long-lived columnar vortices aligned with the axis of rotation. These vortices shut off the elliptical instability, and the net result is insufficient energy dissipation to account for tidal dissipation. However, further work is required to account for effects neglected here, including magnetic fields, turbulent convection, and realistic boundary conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Magnetic buoyancy instabilities in the presence of magnetic flux pumping at the base of the solar convection zone

We perform idealised numerical simulations of magnetic buoyancy instabilities in a model of the solar tachocline. We introduce a simplified model of magnetic flux pumping in an upper layer (the convection zone), and study the effects of its inclusion on the evolution of buoyancy instabilities in a lower layer (the radiative interior). We study its effects on the instability of both a preconceived magnetic slab and of a shear-generated magnetic layer. In the former, we find that in the regime in which the downward pumping velocity is comparable with the Alfven speed of the magnetic layer, flux pumping is able to hold back the bulk of the magnetic field, with only small pockets of strong field able to rise into the upper layer. In simulations in which the magnetic layer is generated by shear, we find that the shear velocity is not necessarily required to exceed that of the pumping (therefore the kinetic energy of the shear is not required to exceed that of the overlying convection), for strong localised pockets of magnetic field to be produced which can rise into the upper layer. This is because magnetic flux pumping acts to store the field below the interface, allowing it to be amplified both by the shear, and by vortical fluid motions, until pockets of field can achieve sufficient strength to rise into the upper layer. In addition, we find that the interface between the two layers is a natural location for the production of strong vertical gradients in the magnetic field. If these gradients are sufficiently strong to allow the development of magnetic buoyancy instabilities, strong shear is not necessarily required to drive them (c.f. previous work by Vasil & Brummell). We find that the addition of magnetic flux pumping appears to be able to assist shear-driven magnetic buoyancy in producing strong flux concentrations that can rise into the convection zone from the radiative interior.

0 citations0 reviewsNo code linkedOpen access