The possibility of subcritical behaviour in the geostrophic turbulence regime of rapidly rotating thermally driven convection is explored. In this regime a non-local inverse energy transfer may compete with the more traditional and local direct cascade. We show that, even for control parameters for which no inverse cascade has previously been observed, a subcritical transition towards a large-scale vortex state can occur when the system is initialized with a vortex dipole of finite amplitude. This new example of bistability in a turbulent flow, which may not be specific to rotating convection, opens up new avenues for studying energy transfer in strongly anisotropic three-dimensional flows such as atmospheric or oceanic circulations.

Viscoelastic flow past a cylinder is a classic benchmark problem that is not completely understood. Using novel three-dimensional (3D) holographic particle velocimetry, we report three main discoveries of the elastic instability upstream of a single cylinder in viscoelastic channel flow. First, we observe that upstream vortices initiate at the corner between the cylinder and the wall, and grow with increasing flow rate. Second, beyond a critical Weissenberg number, the flow upstream becomes unsteady and switches between two bistable configurations, leading to symmetry breaking in the cylinder axis direction that is highly 3D in nature. Lastly, we find that the disturbance of the elastic instability propagates relatively far upstream via an elastic wave, and is weakly correlated with that in the cylinder wake. The wave speed and the extent of the instability increase with Weissenberg number, indicating an absolute instability in viscoelastic fluids.

Horizontal density layers are commonly observed in stratified turbulence. Recent work (e.g. Taylor & Zhou, J. Fluid Mech., vol. 823, 2017, R5) has reinvigorated interest in the Phillips instability (PI), by which density layers form via negative diffusion if the turbulent buoyancy flux weakens as stratification increases. Theoretical understanding of PI is incomplete, in part because it remains unclear whether and by what mechanism the flux-gradient relationship for a given example of turbulence has the required negative-diffusion property. Furthermore, the difficulty of analysing the flux-gradient relation in evolving turbulence obscures the operating mechanism when layering is observed. These considerations motivate the search for an example of PI that can be analysed clearly. Here PI is shown to occur in two-dimensional Boussinesq sheared stratified turbulence maintained by stochastic excitation. PI is analysed using the second-order S3T closure of statistical state dynamics, in which the dynamics is written directly for statistical variables of the turbulence. The predictions of S3T are verified using nonlinear simulations. This analysis provides theoretical underpinning of PI based on the fundamental equations of motion that complements previous analyses based on phenomenological models of turbulence.

Fluid mechanics and rheology involve many unsolved challenges related to the transport mechanisms of mass, momentum and energy – especially when it comes to realistic, industrially relevant materials. Very interesting are suspensions or granular fluids with solid, particulate ingredients that feature contact mechanics on the micro-scale, which affect the transport properties on the continuum- or macro-scale. Their unique ability to behave as either fluid, or solid or both, can be quantified by non-Newtonian rheological rules, and results in interesting mechanisms such as super-diffusion, shear thickening, fluid–solid transitions (jamming) or relaxation/creep. Focusing on the steady state flow of a granular fluid, one can attempt to answer a long-standing question: how do realistic material properties such as dissipation, stiffness, friction or cohesion influence the rheology of a granular fluid? In a recent paper Macaulay & Rognon (J. Fluid Mech., vol. 858, 2019, R2) shed new light on the effect cohesion can have on mass transport in sheared, sticky granular fluids. On top of the usual diffusive, stochastic modes of transport, cohesion can create and stabilise clusters of particles into bigger agglomerates that carry particles over large distances – either ballistically in the dilute regime, or by their rotation in the dense regime. Importantly, these clusters must not only be larger than the particles (defining the intermediate, meso-scale), but they must also have a finite lifetime, in order to be able to exchange mass with each other, which can seriously enhance transport in sticky granular fluids by rotection, i.e. a combination of rotation and convection.

The unsteady axisymmetric flow through a circular aperture in a thin plate subjected to harmonic forcing (for instance under the effect of an incident acoustic wave) is a classical problem first considered by Howe (Proc. R. Soc. Lond. A, vol. 366, 1979, pp. 205–223), using an inviscid model. The purpose of this work is to reconsider this problem through a numerical resolution of the incompressible linearized Navier–Stokes equations (LNSE) in the laminar regime, corresponding to $Re=[500,5000]$. We first compute a steady base flow which allows us to describe the vena contracta phenomenon in agreement with experiments. We then solve a linear problem allowing us to characterize both the spatial amplification of the perturbations and the impedance (or equivalently the Rayleigh conductivity), which is a key quantity to investigate the response of the jet to acoustic forcing. Since the linear perturbation is characterized by a strong spatial amplification, the numerical resolution requires the use of a complex mapping of the axial coordinate in order to enlarge the range of Reynolds number investigated. The results show that the impedances computed with $Re\gtrsim 1500$ collapse onto a single curve, indicating that a large Reynolds number asymptotic regime is effectively reached. However, expressing the results in terms of conductivity leads to substantial deviation with respect to Howe’s model. Finally, we investigate the case of finite-amplitude perturbations through direct numerical simulations (DNS). We show that the impedance predicted by the linear approach remains valid for amplitudes up to order $10^{-1}$, despite the fact that the spatial evolution of the perturbations in the jet is strongly nonlinear.

The stability properties of selected flow configurations, usually denoted as base flows, can be significantly altered by small modifications of the flow, which can be caused, for instance, by a non-intrusive passive control. This aspect is amply demonstrated in the literature by ad hoc sensitivity studies which, however, focus on configurations characterised by a steady base flow. Nevertheless, several flow configurations of interest are characterised by a time-periodic base flow. To this purpose, we propose here an original theoretical framework suitable to quantify the effects of base-flow variations in the stability properties of saturated time-periodic limit cycles. In particular, starting from a Floquet analysis of the linearised Navier–Stokes equations and using adjoint methods, it is possible to estimate the variation of a selected Floquet exponent caused by a generic structural perturbation of the base-flow equations. This link is expressed concisely using the adjoint operators coming from the analysis, and the final result, when applied to spatially localised disturbances, is used to build spatial sensitivity and control maps. These maps identify the regions of the flow where the placement of a infinitesimal small object produces the largest effect on the Floquet exponent and may also provide a quantification of this effect. Such analysis brings useful insights both for passive control strategies and for further characterising the investigated instability. As an example of application, the proposed analysis is applied here to the three-dimensional flow instabilities in the wake past a circular cylinder. This is a classical problem which has been widely studied in the literature. Nevertheless, by applying the proposed analysis we derive original results comprising a further characterisation of the instability and related control maps. We finally show that the control maps obtained here are in very good agreement with control experiments documented in the literature.

We perform an experimental and theoretical study to investigate the interaction between an external harmonic excitation and a self-excited oscillatory mode ($f_{n0}$) of a prototypical thermoacoustic system, a horizontal Rijke tube. Such an interaction can lead to forced synchronization through the routes of phase locking or suppression. We characterize the transition in the synchronization behaviour of the forcing and the response signals of the acoustic pressure while the forcing parameters, i.e. amplitude ($A_{f}$) and frequency ($f_{f}$) of forcing are independently varied. Further, suppression is categorized into synchronous quenching and asynchronous quenching depending upon the value of frequency detuning ($|\,f_{n0}-f_{f}|$). When the applied forcing frequency is close to the natural frequency of the system, the suppression in the amplitude of the self-excited oscillation is known as synchronous quenching. However, this suppression is associated with resonant amplification of the forcing signal, leading to an overall increase in the response amplitude of oscillations. On the other hand, an almost 80 % reduction in the root mean square value of the response oscillation is observed when the system is forced for a sufficiently large value of the frequency detuning (only for $f_{f}<f_{n0}$). Such a reduction in amplitude occurs due to asynchronous quenching where resonant amplification of the forcing signal does not occur, as the frequency detuning is significantly high. Further, the results from a reduced-order model developed for a horizontal Rijke tube show a qualitative agreement with the dynamics observed in experiments. The relative phase between the acoustic pressure ($p^{\prime }$) and the heat release rate ($\dot{q}^{\prime }$) oscillations in the model explains the occurrence of maximum reduction in the pressure amplitude due to asynchronous quenching. Such a reduction occurs when the positive coupling between $p^{\prime }$ and $\dot{q}^{\prime }$ is disrupted and their interaction results in overall acoustic damping, although both of them oscillate at the forcing frequency. Our study on the phenomenon of asynchronous quenching thus presents new possibilities to suppress self-sustained oscillations in fluid systems in general.

Currently, the most important preparation process for porous polymer membranes is the phase inversion process. While applied for several decades in industry, the mechanism that leads to diverse morphology is not fully understood today. In this work, we present time resolved experiments using light microscopy that indicate viscous fingering during the early stage of pore formation in porous polymer membranes. Numerical simulations using the smoothed particle hydrodynamics method are also performed based on Cahn–Hilliard and Navier–Stokes equations to investigate the formation of viscous fingers in miscible and immiscible systems. The comparison of pore formation characteristics in the experiment and simulation shows that immiscible viscous fingering is present; however, it is only relevant in specific preparation set-ups similar to Hele-Shaw cells. In experiments, we also observe the formation of Liesegang rings. Enabling diffusive mass transport across the immiscible interface leads to Liesegang rings in the simulation. We conclude that further investigations of Liesegang pattern as a relevant mechanism in the formation of morphology in porous polymer membranes are necessary.

Rectangular air jets of aspect ratio $2$ are studied at $Re=UD_{h}/\unicode[STIX]{x1D708}=17\,750$ using particle image velocimetry and hot-wire anemometry as they develop naturally or under acoustic forcing. The velocity spectra and the spatial theory of linear stability characterize the fundamental ($f_{n}$) and subharmonic ($f_{n}/2$) modes corresponding to the Kelvin–Helmholtz roll-up and vortex pairing, respectively. The rectangular cross-section of the jet deforms into elliptic/circular shapes downstream due to axis switching. Despite the apparent rotation of the vortex rings or the jet cross-section, the axis-switching phenomenon occurs due to reshaping into rounder geometries. By enhancing the vortex pairing, excitation at $f_{n}/2$ shortens the potential core, increases the jet spread rate and eliminates the overshoot typically observed in the centreline velocity fluctuations. Unlike circular jets, the skewness and kurtosis of the rectangular jets demonstrate elevated anisotropy/intermittency levels before the end of the potential core. The axis-switching location is found to be variable by the acoustic control of the relative expansion/contraction rates of the shear layers in the top (longer edge), side (shorter edge) and diagonal views. The self-induced vortex deformations are demonstrated by the spatio-temporal evolution of the phase-locked three-dimensional ring structures. The curvature-induced velocities are found to reshape the vortex ring by imposing nonlinear azimuthal perturbations occurring at shorter wavelengths with time/space evolution. Eventually, the multiple high-curvature/high-velocity regions merge into a single mode distribution. In the plane of the top view, the self-induced velocity distribution evolves symmetrically while the tilted ring results in the asymmetry of the azimuthal perturbations in the side view as the side closer to the acoustic source rolls up in more upstream locations.

The near wakes of flows past single- and multi-scale arrays of bars are studied by means of planar laser induced fluorescence (PLIF) and particle image velocimetry (PIV). The aim of this research is to better understand dispersion of passive scalar downstream of the multi-scale turbulence generator. In particular, the focus is on plausible manifestations of the space-scale unfolding (SSU) mechanism, which is often considered in the literature as the reason for the enhancement of the turbulent scalar flux in flows past fractal grids (i.e. specific multi-scale turbulence generators). The analysis of qualitative and quantitative PLIF results, as well as the simultaneously acquired PIV results, confirms the appearance of a physical scenario resembling the SSU mechanism. Unlike the anticipation of the literature, however, this scenario applies to some extent also to the flow past the single-scale obstacle. Application of a triple decomposition technique (which splits the acquired fields into their means, a number of coherent fluctuations and their stochastic parts) and a conditional-averaging technique reveals that the SSU mechanism is active in the vicinity of an intersection point between two adjacent wakes and is driven almost exclusively by coherent fluctuations associated with the larger of the intersecting wakes. This suggests that the SSU mechanism is related to the coherent fluctuations embedded in the flow rather than to the fine-scale turbulence and its underlying integral length scale, as proposed in previous works.

A restricted nonlinear (RNL) model, obtained by partitioning the state variables into streamwise-averaged quantities and superimposed perturbations, is used in order to track the exact coherent state in plane channel flow investigated by Toh & Itano (J. Fluid Mech., vol. 481, 2003, pp. 67–76). When restricting nonlinearities to quadratic interaction of the fluctuating part into the streamwise-averaged component, it is shown that the coherent structure and its dynamics closely match results from direct numerical simulation (DNS), even if only a single streamwise Fourier mode is retained. In particular, both solutions exhibit long quiescent phases, spanwise shifts and bursting events. It is also shown that the dynamical trajectory passes close to equilibria that exhibit either low- or high-drag states. When statistics are collected at times where the friction velocity peaks, the mean flow and root-mean-square profiles show the essential features of wall turbulence obtained by DNS for the same friction Reynolds number. For low-drag events, the mean flow profiles are related to a universal asymptotic state called maximum drag reduction (Xi & Graham, Phys. Rev. Lett., vol. 108, 2012, 028301). Hence, the intermittent nature of self-sustaining processes in the buffer layer is contained in the dynamics of the RNL model, organized in two exact coherent states plus an asymptotic turbulent-like attractor. We also address how closely turbulent dynamics approaches these equilibria by exploiting a DNS database associated with a larger domain.

The problem of homogeneous isotropic turbulence (HIT) is revisited within the analytical framework of the Navier–Stokes equations, with a view to assessing rigorously the consequences of the scale invariance (an exact property of the Navier–Stokes equations) for any Reynolds number. The analytical development, which is independent of the 1941 (K41) and 1962 (K62) theories of Kolmogorov for HIT for infinitely large Reynolds number, is applied to the transport equations for the second- and third-order moments of the longitudinal velocity increment, $(\unicode[STIX]{x1D6FF}u)$. Once the normalised equations and the constraints required for complying with the scale-invariance property of the equations are presented, results derived from these equations and constraints are discussed and compared with measurements. It is found that the fluid viscosity, $\unicode[STIX]{x1D708}$, and the mean kinetic energy dissipation rate, $\overline{\unicode[STIX]{x1D716}}$ (the overbar denotes spatial and/or temporal averaging), are the only scaling parameters that make the equations scale-invariant. The analysis further leads to expressions for the distributions of the skewness and the flatness factor of $(\unicode[STIX]{x1D6FF}u)$ and shows that these distributions must exhibit plateaus (of different magnitudes) in the dissipative and inertial ranges, as the Taylor microscale Reynolds number $Re_{\unicode[STIX]{x1D706}}$ increases indefinitely. Also, the skewness and flatness factor of the longitudinal velocity derivative become constant as $Re_{\unicode[STIX]{x1D706}}$ increases; this is supported by experimental data. Further, the analysis, backed up by experimental evidence, shows that, beyond the dissipative range, the behaviour of $\overline{(\unicode[STIX]{x1D6FF}u)^{n}}$ with $n=2$, 3 and 4 cannot be represented by a power law of the form $r^{\unicode[STIX]{x1D701}_{n}}$ when the Reynolds number is finite. It is shown that only when $Re_{\unicode[STIX]{x1D706}}\rightarrow \infty$ can an $n$-thirds law (i.e. $\overline{(\unicode[STIX]{x1D6FF}u)^{n}}\sim r^{\unicode[STIX]{x1D701}_{n}}$, with $\unicode[STIX]{x1D701}_{n}=n/3$) emerge, which is consistent with the onset of a scaling range.

The unsteady dynamics of wall-mounted flexible plates under inclined flows was fundamentally described using theoretical arguments and experiments under various Cauchy numbers $Ca=\unicode[STIX]{x1D70C}_{f}bL^{3}U_{0}^{2}/(EI)\in [7,81]$ (where $\unicode[STIX]{x1D70C}_{f}$ is the fluid density, $b$ and $L$ are the plate width and length, $U_{0}$ is the incoming velocity, $E$ is Young’s modulus and $I$ is the second moment of the area) and inclination angles $\unicode[STIX]{x1D6FC}$. Three-dimensional particle tracking velocimetry and a high-resolution force sensor were used to characterize the evolution of the plate dynamics and aerodynamic force. We show the existence of three distinctive, dominant modes of tip oscillations, which are modulated by the structure dynamic and flow instability. The first mode is characterized by small-amplitude, planar fluttering-like motions occurring under a critical Cauchy number, $Ca=Ca_{c}$. Past this condition, the motions are dominated by the second mode consisting of unsteady twisting superimposed onto the fluttering patterns. The onset of this mode is characterized by a sharp increase of the force fluctuation intensity. At sufficiently high $Ca$ and $\unicode[STIX]{x1D6FC}$, the plate may undergo a third mode given by large-scale tip orbits about the mean bending. Using the equation of motion and first-order approximations, we propose a formulation to estimate $Ca_{c}$ as a function of $\unicode[STIX]{x1D6FC}$; it exhibits solid agreement with experiments.

We study the trajectory and the maximum diving depth of floating axisymmetric streamlined bodies impacting water with a vertical velocity. Three different types of underwater trajectory can be observed. For a centre of mass of the projectile located close to its leading edge, the trajectory is either straight at low velocity or y-shaped at high velocity. When the centre of mass is far from the leading edge, the trajectory has a U-shape, independent of the initial velocity. We first characterize experimentally the aerodynamic properties of the projectile and then solve the equations of motion to recover the three types of trajectories. We finally discuss the transitions between the different regimes.

We study the generation of resonantly growing mean flow by weakly nonlinear internal wave beams. With a perturbational expansion, we construct analytic solutions for three-dimensional internal wave beams, exact up to first-order accuracy in the viscosity parameter. We specifically focus on the subtleties of wave beam generation by oscillating boundaries, such as wave makers in laboratory set-ups. The exact solutions to the linearized equations allow us to derive an analytic expression for the mean vertical vorticity production term, which induces a horizontal mean flow. Whereas mean flow generation associated with viscous beam attenuation – known as streaming – has been described before, we are the first to also include a peculiar inviscid mean flow generation in the vicinity of the oscillating wall, resulting from line vortices at the lateral edges of the oscillating boundary. Our theoretical expression for the mean vertical vorticity production is in good agreement with earlier laboratory experiments, for which the previously unrecognized inviscid mean flow generation mechanism turns out to be significant.

The drag reduction characteristics and sustainability bounds of superhydrophobic (SH) surfaces in high Reynolds number turbulent flows are investigated using results from direct numerical simulation (DNS) and scaling-law analysis. The DNS studies were performed, using lattice Boltzmann methods, in turbulent channel flows at bulk Reynolds numbers of $Re_{b}=3600$ ($Re_{\unicode[STIX]{x1D70F}_{0}}\approx 222$) and $Re_{b}=7860$ ($Re_{\unicode[STIX]{x1D70F}_{0}}\approx 442$) with SH longitudinal microgrooves or SH aligned microposts on the walls. Surface microtexture geometrical parameters corresponding to microgroove widths or micropost spacings of $4\lesssim g^{+0}\lesssim 128$ in base flow wall units and solid fractions of $1/64\leqslant \unicode[STIX]{x1D719}_{s}\leqslant 1/2$ were investigated at interface protrusion angles of $\unicode[STIX]{x1D703}_{p}=0^{\circ }$ and $\unicode[STIX]{x1D703}_{p}=-30^{\circ }$. Analysis of the governing equations and DNS results shows that the magnitude of drag reduction is not only a function of the geometry and size of the surface microtexture in wall units, but also the Reynolds number of the base flow. A Reynolds number independent measure of drag reduction can be constructed by parameterizing the magnitude of drag reduction in terms of the friction coefficient of the base flow and the shift, $(B-B_{0})$, in the intercept of a logarithmic law representation of the mean velocity profile in the flow with SH walls compared to the base flow, where $(B-B_{0})$ is Reynolds number independent. The scaling laws for $(B-B_{0})$, in terms of the geometrical parameters of the surface microtexture in wall units, are presented for SH longitudinal microgrooves and aligned microposts. The same scaling laws are found to also apply to liquid-infused (LI) surfaces as long as the viscosity ratios are large, $N\equiv \unicode[STIX]{x1D707}_{o}/\unicode[STIX]{x1D707}_{i}\gtrsim 10$. These scaling laws, in conjunction with the parametrization of drag reduction in terms of $(B-B_{0})$, allow for a priori prediction of the magnitude of drag reduction with SH or LI surfaces in turbulent flow at any Reynolds number. For the most stable of these SH surface microtextures, namely, longitudinal microgrooves, the pressure stability bounds of the SH surface under the pressure loads of turbulent flow are investigated. It is shown that the pressure stability bounds of SH surfaces are also significantly curtailed with increasing Reynolds number of the flow. Using these scaling laws, the narrow range of SH surface geometrical parameters which can yield large drag reduction as well as sustainability in high Reynolds number turbulent flows is identified.

We examine the flows induced by horizontally modulated, vertically confined (or guided), internal wavepackets in a stratified, Boussinesq fluid. The wavepacket induces both an Eulerian flow and a Stokes drift, which together determine the Lagrangian transport of passive tracers. We derive equations describing the wave-induced flows in arbitrary stable stratification and consider four special cases: a two-layer fluid, symmetric and asymmetric piecewise constant (‘top-hat’) stratification and, more representative of the ocean, exponential stratification. In a two-layer fluid, the Stokes drift is positive everywhere with the peak value at the interface, whereas the Eulerian flow is negative and uniform with depth for long groups. Combined, the net depth-integrated Lagrangian transport is zero. If one layer is shallower than the other, the wave-averaged interface displaces into that layer making the Eulerian flow in that layer more negative and the Eulerian flow in the opposite layer more positive so that the depth-integrated Eulerian transports are offset by the same amount in each layer. By contrast, in continuous stratification the depth-integrated transport due to the Stokes drift and Eulerian flow are each zero, but the Eulerian flow is singular if the horizontal phase speed of the induced flow equals the group velocity of the wavepacket, giving rise to a single resonance in uniform stratification (McIntyre, J. Fluid Mech., vol. 60, 1973, pp. 801–811). In top-hat stratification, this single resonance disappears, being replaced by multiple resonances occurring when the horizontal group velocity of the wavepacket matches the horizontal phase speed of higher-order modes. Furthermore, if the stratification is not vertically symmetric, then the Eulerian induced flow varies as the inverse squared horizontal wavenumber for shallow waves, the same as for the asymmetric two-layer case. This ‘infrared catastrophe’ also occurs in the case of exponential stratification suggesting significant backward near-surface transport by the Eulerian induced flow for modulated oceanic internal modes. Numerical simulations are performed confirming these theoretical predictions.

The mixing induced by breaking internal gravity waves is an important contributor to the ocean’s energy budget, shaping, inter alia, nutrient supply, water mass transformation and the large-scale overturning circulation. Much of the energy input into the internal wave field is supplied by the conversion of barotropic tides at rough bottom topography, which hence needs to be described realistically in internal gravity wave models and mixing parametrisations based thereon. A new semi-analytical method to describe this internal wave forcing, calculating not only the total conversion but also the direction of this energy flux, is presented. It is based on linear theory for variable stratification and finite depth, that is, it computes the energy flux into the different vertical modes for two-dimensional, subcritical, small-amplitude topography and small tidal excursion. A practical advantage over earlier semi-analytical approaches is that the new one gives a positive definite conversion field. Sensitivity studies using both idealised and realistic topography allow the identification of suitable numerical parameter settings and corroborate the accuracy of the method. This motivates the application to the global ocean in order to better account for the geographical distribution of diapycnal mixing induced by low-mode internal gravity waves, which can propagate over large distances before breaking. The first results highlight the significant differences of energy flux magnitudes with direction, confirming the relevance of this more detailed approach for energetically consistent mixing parametrisations in ocean models. The method used here should be applicable to any physical system that is described by the standard wave equation with a very wide field of sources.

An exhaustive description of the dynamics under shear flow of a large number of red blood cells in a dilute regime is proposed, which highlights and takes into account the dispersion in cell properties within a given blood sample. Physiological suspending fluid viscosity is considered, a configuration surprisingly seldom considered in experimental studies, as well as a more viscous fluid that is a reference in the literature. Stable and unstable flipping motions well described by Jeffery orbits or modified Jeffery orbits are identified, as well as transitions to and from tank-treading motion in the more viscous suspending fluid case. Hysteresis loops upon shear rate increase or decrease are highlighted for the transitions between unstable and stable orbits as well as for the transition between flipping and tank-treading. We identify which of the characteristic parameters of motion and of the transition thresholds depend on flow stress only or also on suspending fluid viscosity.

The steady laminar annular spread of a thin liquid film generated by a circular jet which impinges perpendicularly in direction of gravity on the centre of a rotating disc is examined both analytically and numerically. Matched asymptotic expansions of the flow quantities provide the proper means for studying the individual flow regimes arising due to the largeness of the Reynolds number formed with the radius of the jet, its slenderness and the relative magnitude of the centrifugal body force. This is measured by a suitably defined Rossby number, $Ro$. The careful analysis of jet impingement predicts a marked influence of gravity and surface tension on the film flow, considered in the spirit of a shallow-water approach, only through the vorticity imposed by the jet flow. Accordingly, associated downstream conditions are disregarded as the local Froude and Weber numbers are taken to be sufficiently large. Hence, the parabolic problem shaped from the governing equations in a rigorous manner describes the strongly supercritical spread of a developed viscous film past an infinite disc, essentially controlled by $Ro$. Its numerical solutions are discussed for a wide range of values of $Ro$. The different flow regimes reflecting varying effects of viscous shear and centrifugal force are elucidated systematically to clarify the surprising richness of flow phenomena. Special attention is paid to the cases $Ro\gg 1$ and $Ro\ll 1$. The latter, referring to relatively high disc spin, implies a delicate breakdown of the asymptotic flow structure, thus requiring a specific analytical and numerical treatment. Finally, the impact of gravity and capillarity and thus of the disc edge on the film flow is envisaged in brief.