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

Steven M. Tobias

Steven M. Tobias appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dynastro-ph.SRphysics.plasm-phastro-ph.EPphysics.geo-phcond-mat.mes-hallnlin.CDastro-ph.GAastro-ph.HEcond-mat.stat-mechmath-phmath.MPphysics.comp-phphysics.space-ph

Catalog footprint

What is connected

18works
14topics
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

18 published item(s)

preprint2025arXiv

On nonlinear transitions, minimal seeds and exact solutions for the geodynamo

Nearly fifty years ago, Roberts (1978) postulated that Earth's magnetic field, which is generated by turbulent motions of liquid metal in its outer core, likely results from a subcritical (finite-amplitude) dynamo instability characterised by a dominant balance between Coriolis, pressure and Lorentz forces. Here we numerically explore subcritical convective dynamo action in a spherical shell, using techniques from optimal control and dynamical systems theory to uncover the nonlinear dynamics of magnetic field generation. Through nonlinear optimisation, via direct-adjoint looping, we identify the minimal seed - the smallest magnetic field that attracts to a nonlinear dynamo solution. Additionally, using the Newton-hookstep algorithm, we converge stable and unstable travelling wave solutions to the governing equations. By combining these two techniques, complex nonlinear pathways between attracting states are revealed, providing insight into a potential subcritical origin of the geodynamo. This paper showcases these methods on the widely studied benchmark of Christensen et al. (2001), laying the foundations for future studies in more extreme and realistic parameter regimes. We show that the minimal seed reaches a nonlinear dynamo solution by first approaching an unstable travelling wave solution, which acts as an edge state separating a hydrodynamic solution from a magnetohydrodynamic one. Furthermore, by carefully examining the choice of cost functional, we establish a robust optimisation procedure that can systematically locate dynamo solutions on short time horizons with no prior knowledge of its structure.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Magnetoclinicity Instability

In strongly compressible magnetohydrodynamic turbulence, obliqueness between the large-scale density gradient and magnetic field gives an electromotive force mediated by density variance (intensity of density fluctuation). This effect is named ``magnetoclinicity'', and is expected to play an important role in large-scale magnetic-field generation in astrophysical compressible turbulent flows. Analysis of large-scale instability due to the magnetoclinicity effect shows that the mean magnetic-field perturbation is destabilised at large scales in the vicinity of strong mean density gradient in the presence of density variance.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Direct statistical simulation of the Lorenz63 system

We use direct statistical simulation (DSS) to find the low-order statistics of the well-known dynamical system, the Lorenz63 model. Instead of accumulating statistics from numerical simulation of the dynamical systems, we solve the equations of motion for the statistics themselves after closing them by making several different choices for the truncation. Fixed points of the statistics are obtained either by time evolving, or by iterative methods. Statistics so obtained are compared to those found by the traditional approach.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Optimizing the control of transition to turbulence using a Bayesian method

The nonlinear robustness of laminar plane Couette flow is considered under the action of in-phase spanwise wall oscillations by computing properties of the edge of chaos, i.e., the boundary of its basin of attraction. Three measures are used to quantify the chosen control strategy on laminar-to-turbulent transition: the kinetic energy of edge states (local attractors on the edge of chaos), the form of the minimal seed (least energetic perturbation on the edge of chaos), and the laminarization probability (the probability that a random perturbation from the laminar flow of given kinetic energy will laminarize). A novel Bayesian approach is introduced to enable the accurate computation of the laminarization probability at a fraction of the cost of previous methods. While the edge state and the minimal seed provide useful information about the dynamics of transition to turbulence, neither measure is particularly useful to judge the effectiveness of the control strategy since they are not representative of the global geometry of the edge. In contrast, the laminarization probability provides global information about the edge and can be used to evaluate the control effectiveness by computing a laminarization score (the expected laminarization probability) and the associated expected dissipation rate of the controlled flow. These two quantities allow for the determination of optimal control parameter values subject to desired constraints. The results discussed in the paper are expected to be applied to a wide range of transitional flows and control strategies aimed at suppressing or triggering transition to turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A probabilistic protocol for the assessment of transition and control

Transition to turbulence dramatically alters the properties of fluid flows. In most canonical shear flows, the laminar flow is linearly stable and a finite-amplitude perturbation is necessary to trigger transition. Controlling transition to turbulence is achieved via the broadening or narrowing of the basin of attraction of the laminar flow. In this paper, a novel methodology to assess the robustness of the laminar flow and the efficiency of control strategies is introduced. It relies on the statistical sampling of the phase space neighborhood around the laminar flow in order to assess the transition probability of perturbations as a function of their energy. This approach is applied to a canonical flow (plane Couette flow) and provides invaluable insight: in the presence of the chosen control, transition is significantly suppressed whereas plausible scalar indicators of the nonlinear stability of the flow, such as the edge state energy, do not provide conclusive predictions. The methodology presented here in the context of transition to turbulence is applicable to any nonlinear system displaying finite-amplitude instability.

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

Generation of shear flows and vortices in rotating anelastic convection

We consider the effect of stratification on systematic, large-scale flows generated in anelastic convection. We present results from three-dimensional numerical simulations of convection in a rotating plane layer in which the angle between the axis of rotation and gravity is allowed to vary. This model is representative of different latitudes of a spherical body. We consider two distinct parameter regimes: (i) weakly rotating and (ii) rapidly rotating. In each case, we examine the effect of stratification on the flow structure and heat transport properties focussing on the difference between Boussinesq and anelastic convection. Furthermore, we show that regimes (i) and (ii) generate very different large-scale flows and we investigate the role stratification has in modifying these flows. The stratified flows possess a net helicity not present in the Boussinesq cases which we suggest, when combined with the self-generated shear flows, could be important for dynamo action.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nontrivial topology in the continuous spectrum of a magnetized plasma

Classification of matter through topological phases and topological edge states between distinct materials has been a subject of great interest recently. While lattices have been the main setting for these studies, a relatively unexplored realm for this physics is that of continuum fluids. In the typical case of a fluid model with a point spectrum, nontrivial topology and associated edge modes have been observed previously. However, another possibility is that a continuous spectrum can coexist with the point spectrum. Here we demonstrate that a fluid plasma model can harbor nontrivial topology within its continuous spectrum, and that there are boundary modes at the interface between topologically distinct regions. We consider the ideal magnetohydrodynamics (MHD) model. In the presence of magnetic shear, we find nontrivial topology in the Alfvén continuum. For strong shear, the Chern number is $\pm 1$, depending on the sign of the shear. If the magnetic shear changes sign within the plasma, a topological phase transition occurs, and bulk-boundary correspondence then suggests a mode localized to the layer of zero magnetic shear. We confirm the existence of this mode numerically. Moreover, this reversed-shear Alfvén eigenmode (RSAE) is well known within magnetic fusion as it has been observed in several tokamaks. In examining how the MHD model might be regularized at small scales, we also consider the electron limit of Hall MHD. We show that the whistler band, which plays an important role in planetary ionospheres, has nontrivial topology. More broadly, this work raises the possibility that fusion devices could be carefully tailored to produce other topological states with potentially useful behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Topological Gaseous Plasmon Polariton in Realistic Plasma

Nontrivial topology in bulk matter has been linked with the existence of topologically protected interfacial states. We show that a gaseous plasmon polariton (GPP), an electromagnetic surface wave existing at the boundary of magnetized plasma and vacuum, has a topological origin that arises from the nontrivial topology of magnetized plasma. Because a gaseous plasma cannot sustain a sharp interface with discontinuous density, one must consider a gradual density falloff with scale length comparable or longer than the wavelength of the wave. We show that the GPP may be found within a gapped spectrum in present-day laboratory devices, suggesting that platforms are currently available for experimental investigation of topological wave physics in plasmas.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Parallel-in-time integration of Kinematic Dynamos

The precise mechanisms responsible for the natural dynamos in the Earth and Sun are still not fully understood. Numerical simulations of natural dynamos are extremely computationally intensive, and are carried out in parameter regimes many orders of magnitude away from real conditions. Parallelization in space is a common strategy to speed up simulations on high performance computers, but eventually hits a scaling limit. Additional directions of parallelization are desirable to utilise the high number of processor cores now available. Parallel-in-time methods can deliver speed up in addition to that offered by spatial partitioning but have not yet been applied to dynamo simulations. This paper investigates the feasibility of using the parallel-in-time algorithm Parareal to speed up initial value problem simulations of the kinematic dynamo, using the open source Dedalus spectral solver. Both the time independent Roberts and time dependent Galloway-Proctor 2.5D dynamos are investigated over a range of magnetic Reynolds numbers. Speed ups beyond those possible from spatial parallelization are found in both cases. Results for the Galloway-Proctor flow are promising, with Parareal efficiency found to be close to 0.3. Roberts flow results are less efficient, but Parareal still shows some speed up over spatial parallelization alone. Parallel in space and time speed ups of $\sim300$ were found for 1600 cores for the Galloway-Proctor flow, with total parallel efficiency of $\sim0.16$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Convection-driven kinematic dynamos at low Rossby and magnetic Prandtl numbers

Most large-scale planetary magnetic fields are thought to be driven by low Rossby number convection of a low magnetic Prandtl number fluid. Here kinematic dynamo action is investigated with an asymptotic, rapidly rotating dynamo model for the plane layer geometry that is intrinsically low magnetic Prandtl number. The thermal Prandtl number and Rayleigh number are varied to illustrate fundamental changes in flow regime, ranging from laminar cellular convection to geostrophic turbulence in which an inverse energy cascade is present. A decrease in the efficiency of the convection to generate a dynamo, as determined by an increase in the critical magnetic Reynolds number, is observed as the buoyancy forcing is increased. This decreased efficiency may result from both the loss of correlations associated with the increasingly disordered states of flow that are generated, and boundary layer behavior that enhances magnetic diffusion locally. We find that the spatial characteristics of $α$, and thus the large-scale magnetic field, is dependent only weakly on changes in flow behavior. However, our results are limited to the linear, kinematic dynamo regime, and future simulations including the Lorentz force are therefore necessary to assess the robustness of this result. In contrast to the large-scale magnetic field, the behavior of the small-scale magnetic field is directly dependent on, and therefore shows significant variations with, the small-scale convective flow field.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Convection-driven kinematic dynamos at low Rossby and magnetic Prandtl numbers: single mode solutions

The onset of dynamo action is investigated within the context of a newly developed low Rossby, low magnetic Prandtl number, convection-driven dynamo model. This multiscale model represents an asymptotically exact form of an $α^2$ mean field dynamo model in which the small-scale convection is represented explicitly by finite amplitude, single mode solutions. Both steady and oscillatory convection are considered for a variety of horizontal planforms. The kinetic helicity is observed to be a monotonically increasing function of the Rayleigh number. As a result, very small magnetic Prandtl number dynamos can be found for sufficiently large Rayleigh numbers. All dynamos are found to be oscillatory with an oscillation frequency that increases as the strength of the convection is increased and the magnetic Prandtl number is reduced. Kinematic dynamo action is strongly controlled by the profile of the helicity; single mode solutions which exhibit boundary layer behavior in the helicity show a decrease in the efficiency of dynamo action due to the enhancement of magnetic diffusion in the boundary layer regions. For a given value of the Rayleigh number, lower magnetic Prandtl number dynamos are excited for the case of oscillatory convection in comparison to steady convection. With regards to planetary dynamos, these results suggest that the low magnetic Prandtl number dynamos typical of liquid metals are more easily driven by thermal convection than by compositional convection.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Convective dynamo action in a spherical shell: symmetries and modulation

We consider dynamo action driven by three-dimensional rotating anelastic convection in a spherical shell. Motivated by the behaviour of the solar dynamo, we examine the interaction of hydromagnetic modes with different symmetries and demonstrate how complicated interactions between convection, differential rotation and magnetic fields may lead to modulation of the basic cycle. For some parameters, Type 1 modulation occurs by the transfer of energy between modes of different symmetries with little change in the overall amplitude, for other parameters, the modulation is of Type 2, where the amplitude is significantly affected (leading to grand minima in activity) without significant changes in symmetry. Most importantly, we identify the presence of "supermodulation" in the solutions, where the activity switches chaotically between Type 1 and Type 2 modulation, this is believed to be an important process in solar activity.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Mean flow generation in rotating anelastic two-dimensional convection

We investigate the processes that lead to the generation of mean flows in two-dimensional anelastic convection. The simple model consists of a plane layer that is rotating about an axis inclined to gravity. The results are two-fold: firstly we numerically investigate the onset of convection in three-dimensions, paying particular attention to the role of stratification and highlight a curious symmetry. Secondly, we investigate the mechanisms that drive both zonal and meridional flows in two dimensions. We find that, in general, non-trivial Reynolds stresses can lead to systematic flows and, using statistical measures, we quantify the role of stratification in modifying the coherence of these flows.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Multiscale Dynamo Model Driven by Quasi-geostrophic Convection

A convection-driven multiscale dynamo model is developed in the limit of low Rossby number for the plane layer geometry in which the gravity and rotation vectors are aligned. The small-scale fluctuating dynamics are described by a magnetically-modified quasi-geostrophic equation set, and the large-scale mean dynamics are governed by a diagnostic thermal wind balance. The model utilizes three timescales that respectively characterize the convective timescale, the large-scale magnetic evolution timescale, and the large-scale thermal evolution timescale. Distinct equations are derived for the cases of order one and low magnetic Prandtl number. It is shown that the low magnetic Prandtl number model is characterized by a magnetic to kinetic energy ratio that is asymptotically large, with ohmic dissipation dominating viscous dissipation on the large-scales. For the order one magnetic Prandtl number model the magnetic and kinetic energies are equipartitioned and both ohmic and viscous dissipation are weak on the large-scales; large-scale ohmic dissipation occurs in thin magnetic boundary layers adjacent to the solid boundaries. For both magnetic Prandtl number cases the Elsasser number is small since the Lorentz force does not enter the leading order force balance. The new models can be considered fully nonlinear, generalized versions of the dynamo model originally developed by Childress and Soward [Phys. Rev. Lett., \textbf{29}, p.837, 1972]. These models may be useful for understanding the dynamics of convection-driven dynamos in regimes that are only just becoming accessible to direct numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Energy Dissipation in Magnetohydrodynamic Turbulence: Coherent Structures or "Nanoflares"?

We investigate the intermittency of energy dissipation in magnetohydrodynamic (MHD) turbulence by identifying dissipative structures and measuring their characteristic scales. We find that the probability distribution of energy dissipation rates exhibits a power law tail with index very close to the critical value of -2.0, which indicates that structures of all intensities contribute equally to energy dissipation. We find that energy dissipation is uniformly spread among coherent structures with lengths and widths in the inertial range. At the same time, these structures have thicknesses deep within the dissipative regime. As the Reynolds number is increased, structures become thinner and more numerous, while the energy dissipation continues to occur mainly in large-scale coherent structures. This implies that in the limit of high Reynolds number, energy dissipation occurs in thin, tightly packed current sheets which nevertheless span a continuum of scales up to the system size, exhibiting features of both coherent structures and nanoflares previously conjectured as a coronal heating mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The dynamics and excitation of torsional waves in geodynamo simulations

The predominant force balance in rapidly rotating planetary cores is between Coriolis, pressure, buoyancy and Lorentz forces. This magnetostrophic balance leads to a Taylor state where the spatially averaged azimuthal Lorentz force is compelled to vanish on cylinders aligned with the rotation axis. Any deviation from this state leads to a torsional oscillation, signatures of which have been observed in the Earth's secular variation and are thought to influence length of day variations via angular momentum conservation. In order to investigate the dynamics of torsional oscillations, we perform several three-dimensional dynamo simulations in a spherical shell. We find torsional oscillations, identified by their propagation at the correct Alfvén speed, in many of our simulations. We find that the frequency, location and direction of propagation of the waves are influenced by the choice of parameters. Torsional waves are observed within the tangent cylinder and also have the ability to pass through it. Several of our simulations display waves with core travel times of 4 to 6 years. We calculate the driving terms for these waves and find that both the Reynolds force and ageostrophic convection acting through the Lorentz force are important in driving torsional oscillations.

0 citations0 reviewsNo code linkedOpen access