Researcher profile

Rich R. Kerswell

Rich R. Kerswell contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
2topics
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

11 published item(s)

preprint2021arXiv

A first coherent structure in elasto-inertial turbulence

Two dimensional channel flow simulations of FENE-P fluid in the elasto-inertial turbulence regime reveal distinct regimes ranging from chaos to a steady travelling wave which takes the form of an arrowhead structure. This coherent structure provides new insights in the polymer/flow interactions driving EIT, which are observed in a set of controlled numerical experiments and the study of transfer between elastic and turbulent kinetic energy.

preprint2021arXiv

Weakly nonlinear analysis of the viscoelastic instability in channel flow for finite and vanishing Reynolds numbers

The recently-discovered centre-mode instability of rectilinear viscoelastic shear flow (Garg et al. Phy. Rev. Lett. 121, 024502, 2018) has offered an explanation for the origin of elasto-inertial turbulence (EIT) which occurs at lower Weissenberg ($Wi$) numbers. In support of this, we show using weakly nonlinear analysis that the subcriticality found in Page et al. (Phys. Rev. Lett. 125, 154501, 2020) is generic across the neutral curve with the instability only becoming supercritical at low Reynolds ($Re$) numbers and high $Wi$. We demonstrate that the instability can be viewed as purely elastic in origin even for $Re=O(10^3)$, rather than `elasto-inertial', as the underlying shear does not energise the instability. It is also found that the introduction of a realistic maximum polymer extension length, $L_{max}$, in the FENE-P model moves the neutral curve closer to the inertialess $Re=0$ limit at a fixed ratio of solvent-to-solution viscosities, $β$. In the dilute limit ($β\rightarrow 1$) with $L_{max} =O(100)$, the linear instability can brought down to more physically-relevant $Wi\gtrsim 110$ at $β=0.98$, compared with the threshold $Wi=O(10^3)$ at $β=0.994$ reported recently by Khalid et al. (arXiv: 2103.06794) for an Oldroyd-B fluid. Again the instability is subcritical implying that inertialess rectilinear viscoelastic shear flow is nonlinearly unstable - i.e. unstable to finite amplitude disturbances - for even lower $Wi$.

preprint2020arXiv

Designing a minimal baffle to destabilise turbulence in pipe flows

Motivated by the results of recent experiments (Kühnen et al., Flow Turb. Combust., vol. 100, 2018, pp. 919-943), we consider the problem of designing a baffle (an obstacle to the flow) to relaminarise turbulence in pipe flows. Modelling the baffle as a spatial distribution of linear drag $\mathbf{f}(\mathbf{x},t)=-χ(\mathbf{x})\mathbf{u}_{tot}(\mathbf{x},t)$ within the flow ($\mathbf{u}_{tot}$ is the total velocity field and $χ\ge 0$ a scalar field), two different optimisation problems are considered to design $χ$ at a Reynolds number $Re=3000$. In the first, the smallest baffle defined in terms of a $L_1 $ norm of $χ$ is sought which minimises the viscous dissipation rate of the flow. In the second, a baffle which minimises the total energy consumption of the flow is treated. Both problems indicate that the baffle should be axisymmetric and radially localised near the pipe wall, but struggle to predict the optimal streamwise extent. A manual search finds an optimal baffle one radius long which is then used to study how the amplitude for relaminarisation varies with $Re$ up to $15\,000$. Large stress reduction is found at the pipe wall, but at the expense of an increased pressure drop across the baffle. Estimates are then made of the break even point downstream of the baffle where the stress reduction at the wall due to the relaminarised flow compensates for the extra drag produced by the baffle.

preprint2019arXiv

Searching turbulence for periodic orbits with dynamic mode decomposition

We present a new method for generating robust guesses for unstable periodic orbits (UPOs) by post-processing turbulent data using dynamic mode decomposition (DMD). The approach relies on the identification of near-neutral, repeated harmonics in the DMD eigenvalue spectrum from which both an estimate for the period of a nearby UPO and a guess for the velocity field can be constructed. In this way, the signature of a UPO can be identified in a short time series without the need for a near recurrence to occur, which is a considerable drawback to recurrent flow analysis, the current state-of-the-art. We first demonstrate the method by applying it to a known (simple) UPO and find that the period can be reliably extracted even for time windows of length one quarter of the full period. We then turn to a long turbulent trajectory, sliding an observation window through the time series and performing many DMD computations. Our approach yields many more converged periodic orbits (including multiple new solutions) than a standard recurrent flow analysis of the same data. Furthermore, it also yields converged UPOs at points where the recurrent flow analysis flagged a near recurrence but the Newton solver did not converge, suggesting that the new approach can be used alongside the old to generate improved initial guesses. Finally, we discuss some heuristics on what constitutes a "good" time window for the DMD to identify a UPO.

preprint2019arXiv

Stably stratified exact coherent structures in shear flow: the effect of Prandtl number

We examine how known unstable equilibria of the Navier-Stokes equations in plane Couette flow adapt to the presence of an imposed stable density difference between the two boundaries for varying values of the Prandtl number $Pr$, the ratio of viscosity to density diffusivity, and fixed moderate Reynolds number, $Re=400$. In the two asymptotic limits $Pr \to 0$ and $Pr \to \infty$, it is found that such solutions exist at arbitrarily high bulk stratification but for different physical reasons. In the $Pr \to 0$ limit, density variations away from a constant stable density gradient become vanishingly small as diffusion of density dominates over advection, allowing equilibria to exist for bulk Richardson number $Ri_b \lesssim O(Re^{-2}Pr^{-1})$. Alternatively, at high Prandtl numbers, density becomes homogenised in the interior by the dominant advection which creates strongly stable stratified boundary layers that recede into the wall as $Pr\to\infty$. In this scenario, the density stratification and the flow essentially decouple, thereby mitigating the effect of increasing $Ri_b$. An asymptotic analysis is presented in the passive scalar regime $Ri_b \lesssim O(Re^{-2})$, which reveals $O(Pr^{-1/3})$-thick stratified boundary layers with $O(Pr^{-2/9})$-wide eruptions, giving rise to density fingers of $O(Pr^{-1/9})$ length and $O(Pr^{-4/9})$ width that invade an otherwise homogeneous interior. Finally, increasing $Re$ to $10^5$ in this regime reveals that interior stably stratified density layers can form away from the boundaries, separating well-mixed regions.

preprint2019arXiv

The viscous Holmboe instability for smooth shear and density profiles

The Holmboe wave instability is one of the classic examples of a stratified shear instability, usually explained as the result of a resonance between a gravity wave and a vorticity wave. Historically, it has been studied by linear stability analyses at infinite Reynolds number, $Re$, and by direct numerical simulations at relatively low $Re$ in the regions known to be unstable from the inviscid linear stability results. In this paper, we perform linear stability analyses of the classical `Hazel model' of a stratified shear layer (where the background velocity and density distributions are assumed to take the functional form of hyperbolic tangents with different characteristic vertical scales) over a range of different parameters at finite $Re$, finding new unstable regions of parameter space. In particular, we find instability when the Richardson number is everywhere greater than $1/4$, where the flow would be stable at infinite $Re$ by the Miles-Howard theorem. We find unstable modes with no critical layer, and show that despite the necessity of viscosity for the new instability, the growth rate relative to diffusion of the background profile is maximised at large $Re$. We use these results to shed new light on the wave-resonance and over-reflection interpretations of stratified shear instability. We argue for a definition of Holmboe instability as being characterised by propagating vortices above or below the shear layer, as opposed to any reference to sharp density interfaces.

preprint2014arXiv

Fully localised nonlinear energy growth optimals in pipe flow

A new, fully-localised, energy growth optimal is found over large times and in long pipe domains at a given mass flow rate. This optimal emerges at a threshold disturbance energy below which a nonlinear version of the known (streamwise-independent) linear optimal (Schmid \& Henningson 1994) is selected, and appears to remain the optimal up until the critical energy at which transition is triggered. The form of this optimal is similar to that found in short pipes (Pringle et al.\ 2012) albeit now with full localisation in the streamwise direction. This fully-localised optimal perturbation represents the best approximation yet of the {\em minimal seed} (the smallest perturbation capable of triggering a turbulent episode) for `real' (laboratory) pipe flows.

preprint2014arXiv

The genesis of streamwise-localized solutions from globally periodic travelling waves in pipe flow

The aim in the dynamical systems approach to transitional turbulence is to construct a scaffold in phase space for the dynamics using simple invariant sets (exact solutions) and their stable and unstable manifolds. In large (realistic) domains where turbulence can co-exist with laminar flow, this requires identifying exact localized solutions. In wall-bounded shear flows the first of these has recently been found in pipe flow, but questions remain as to how they are connected to the many known streamwise-periodic solutions. Here we demonstrate the origin of the first localized solution in a modulational symmetry-breaking Hopf bifurcation from a known global travelling wave that has 2-fold rotational symmetry about the pipe axis. Similar behaviour is found for a global wave of 3-fold rotational symmetry, this time leading to two localized relative periodic orbits. The clear implication is that all global solutions should be expected to lead to more realistic localised counterparts through such bifurcations, which provides a constructive route for their generation.

preprint2012arXiv

Simple invariant solutions embedded in 2D Kolmogorov turbulence

We consider long simulations of 2D Kolmogorov turbulence body-forced by $\sin4y \ex$ on the torus $(x,y) \in [0,2π]^2$ with the purpose of extracting simple invariant sets or `exact recurrent flows' embedded in this turbulence. Each recurrent flow represents a sustained closed cycle of dynamical processes which underpins the turbulence. These are used to reconstruct the turbulence statistics in the spirit of Periodic Orbit Theory derived for certain types of low dimensional chaos. The approach is found to be reasonably successful at a low value of the forcing where the flow is close to but not fully in its asymptotic (strongly) turbulent regime. Here, a total of 50 recurrent flows are found with the majority buried in the part of phase space most populated by the turbulence giving rise to a good reproduction of the energy and dissipation probability density functions. However, at higher forcing amplitudes now in the asymptotic turbulent regime, the generated turbulence data set proves insufficiently long to yield enough recurrent flows to make viable predictions. Despite this, the general approach seems promising providing enough simulation data is available since it is open to extensive automation and naturally generates dynamically important exact solutions for the flow.

preprint2010arXiv

On Constructing the Minimal Seed of Turbulence: Nonlinear Transient Growth

Linear transient growth analysis is commonly used to suggest the structure of disturbances which are particularly efficient in triggering transition to turbulence in shear flows. We demonstrate that the addition of nonlinearity to the analysis can substantially change the prediction made in pipe flow from simple 2 dimensional streamwise rolls to a spanwise and cross-stream localised 3 dimensional state. This new nonlinear optimal is demonstrably more efficient in triggering turbulence than the linear optimal indicating that there are better ways to design perturbations to achieve transition.

preprint2010arXiv

Transient Growth in Shear Flows: Linearity vs Nonlinearity

Two approaches to the problem of transition to turbulence of shear flows are popular in the literature. The first is the linear one of transient growth which focuses on the likely form of the most 'dangerous' (lowest energy) turbulence-triggering disturbances. The second is the nonlinear calculation of the laminar-turbulent boundary which instead focuses on their typical amplitudes. We look to bridge the gap between these two perspectives by considering the fully nonlinear transient growth problem to estimate both the form and amplitude of the most dangerous disturbance. We thereby discover a new nonlinear optimal disturbance which outgrows the well-known linear optimal for the same initial energy and is crucially much more efficient in triggering turbulence. The conclusion is then that the most dangerous disturbance can differ markedly from what traditional linear transient growth analysis predicts.