BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Steven M. Tobias

Steven M. Tobias contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dynastro-ph.SRphysics.plasm-phastro-ph.EPcond-mat.mes-hallnlin.CDcond-mat.stat-mechphysics.comp-phphysics.geo-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 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