Source author record

Steven Tobias

Steven 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-dyn astro-ph.SR physics.comp-ph astro-ph.EP astro-ph.HE physics.ao-ph

Catalog footprint

What is connected

6works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2021arXiv

Efficiency gains of a multi-scale integration method applied to a scale-separated model for rapidly rotating dynamos

Numerical geodynamo simulations with parameters close to an Earth-like regime would be of great interest for understanding the dynamics of the Earth's liquid outer core and the associated geomagnetic field. Such simulations are far too computationally demanding owing to the large range in spatiotemporal scales. This paper explores the application of a multi-scale timestepping method to an asymptotic model for the generation of magnetic field in the fluid outer core of the Earth. The method is based on the heterogeneous multiscale modelling (HMM) strategy, which incorporates scale separation and utilizes several integrating models for the fast and slow fields. Quantitative comparisons between the multi-scale simulations and direct solution of the asymptotic model in the limit of rapid rotation and low Ekman number are performed. The multi-scale method accurately captures the varying temporal and spatial dynamics of the mean magnetic field at lower computational costs compared to the direct solution of the asymptotic model.

preprint2020arXiv

Performance of parallel-in-time integration for Rayleigh Bénard Convection

Rayleigh-Bénard convection (RBC) is a fundamental problem of fluid dynamics, with many applications to geophysical, astrophysical, and industrial flows. Understanding RBC at parameter regimes of interest requires complex physical or numerical experiments. Numerical simulations require large amounts of computational resources; in order to more efficiently use the large numbers of processors now available in large high performance computing clusters, novel parallelisation strategies are required. To this end, we investigate the performance of the parallel-in-time algorithm Parareal when used in numerical simulations of RBC. We present the first parallel-in-time speedups for RBC simulations at finite Prandtl number. We also investigate the problem of convergence of Parareal with respect to to statistical numerical quantities, such as the Nusselt number, and discuss the importance of reliable online stopping criteria in these cases.

preprint2016arXiv

Generalised quasi-linear approximation of the HMRI

Motivated by recent advances in Direct Statistical Simulation (DSS) of astrophysical phenomena such as out of equilibrium jets, we perform a Direct Numerical Simulation (DNS) of the helical magnetorotational instability (HMRI) under the generalised quasilinear approximation (GQL). This approximation generalises the quasilinear approximation (QL) to include the self-consistent interaction of large-scale modes, interpolating between fully nonlinear DNS and QL DNS whilst still remaining formally linear in the small scales. In this paper we address whether GQL can more accurately describe low-order statistics of axisymmetric HMRI when compared with QL by performing DNS under various degrees of GQL approximation. We utilise various diagnostics, such as energy spectra in addition to first and second cumulants, for calculations performed for a range of Reynolds and Hartmann numbers (describing rotation and imposed magnetic field strength respectively). We find that GQL performs significantly better than QL in describing the statistics of the HMRI even when relatively few large-scale modes are kept in the formalism. We conclude that DSS based on GQL (GCE2) will be significantly more accurate than that based on QL (CE2).

preprint2016arXiv

Three-Dimensional Rotating Couette Flow Via The Generalised Quasilinear Approximation

We examine the effectiveness of the Generalised Quasilinear (GQL) Approximation introduced by Marston et al (2016). This approximation splits the variables into large and small scales in directions where there is a translational symmetry and removes nonlinear interactions involving only small scales. We utilise as a paradigm problem three-dimensional, turbulent, rotating Couette flow. We compare the results obtained from Direct Numerical Solution of the equations with those from Quasilinear (QL) and GQL calculations. In this three-dimensional setting, there is a choice of cut-off wavenumber for the GQL approximation both in the streamwise and in the spanwise directions. We demonstrate that the GQL approximation significantly improves the accuracy of mean flows, spectra and two-point correlation functions over models that are quasilinear in any of the translationally invariant directions, even if only a few streamwise and spanwise modes are included. We argue that this provides significant support for a programme of Direct Statistical Simulation utilising the GQL approximation.

preprint2014arXiv

On Large-Scale Dynamo Action at High Magnetic Reynolds Number

We consider the generation of magnetic activity --- dynamo waves --- in the astrophysical limit of very large magnetic Reynolds number. We consider kinematic dynamo action for a system consisting of helical flow and large-scale shear. We demonstrate that large-scale dynamo waves persist at high $Rm$ if the helical flow is characterised by a narrow band of spatial scales and the shear is large enough. However for a wide band of scales the dynamo becomes small-scale with a further increase of $Rm$, with dynamo waves re-emerging only if the shear is then increased. We show that at high $Rm$ the key effect of the shear is to suppress small-scale dynamo action, allowing large-scale dynamo action to be observed. We conjecture that this supports a general "suppression principle" --- large-scale dynamo action can only be observed if there is a mechanism that suppresses the small-scale fluctuations.

preprint2013arXiv

Direct Statistical Simulation of Out-of-Equilibrium Jets

We present a Direct Statistical Simulation (DSS) of jet formation on a β-plane, solving for the statistics of a fluid flow via an expansion in cumulants. Here we compare an expansion truncated at second order (CE2) to statistics accumulated by direct numerical simulations (DNS). We show that, for jets near equilibrium, CE2 is capable of reproducing the jet structure (although some differences remain in the second cumulant). However as the degree of departure from equilibrium is increased (as measured by the zonostrophy parameter) the jets meander more and CE2 becomes less accurate. We discuss a possible remedy by inclusion of higher cumulants.