We experimentally and numerically investigate the temporal aspects of turbulent spots spreading in a plane Couette flow for transitional Reynolds numbers between 300 and 450. Spot growth rate, spot advection rate and large-scale flow intensity are measured as a function of time and Reynolds number. All these quantities show similar dynamics clarifying the role played by large-scale flows in the advection of the turbulent spot. The contributions of each possible growth mechanism, that is, growth induced by large-scale advection or growth by destabilization, are discussed for the different stages of the spot growth. A scenario that gathers all these elements is providing a better understanding of the growth dynamics of turbulent spots in plane Couette flow that should possibly apply to other extended shear flows.

Experiments and numerical simulations of cavitating flow inside a single-orifice nozzle are presented. The orifice is part of a closed flow circuit, with diesel fuel as the working fluid, designed to replicate the main flow pattern observed in high-pressure diesel injector nozzles. The focus of the present investigation is on cavitation structures appearing inside the orifice, their interaction with turbulence and the induced material erosion. Experimental investigations include high-speed shadowgraphy visualization, X-ray micro-computed tomography (micro-CT) of time-averaged volumetric cavitation distribution inside the orifice as well as pressure and flow rate measurements. The highly transient flow features that are taking place, such as cavity shedding, collapse and vortex cavitation (also known as ‘string cavitation’), have become evident from high-speed images. Additionally, micro-CT enabled the reconstruction of the orifice surface, which provided locations of cavitation erosion sites developed after sufficient operation time. The measurements are used to validate the presented numerical model, which is based on the numerical solution of the Navier–Stokes equation, taking into account compressibility of both the liquid and liquid–vapour mixture. Phase change is accounted for with a newly developed mass transfer rate model, capable of accurately predicting the collapse of vaporous structures. Turbulence is modelled using detached eddy simulation and unsteady features such as cavitating vortices and cavity shedding are observed and discussed. The numerical results show agreement within validation uncertainty with the obtained measurements.

The transport of passive scalars in turbulent flows over progressive water waves is studied using direct numerical simulation. A combined pseudo-spectral and finite-difference scheme on a wave-surface-fitted grid is used to simulate the flow and scalar fields above the wave surface. Three representative wave ages (i.e. wave-to-wind speed ratios) are considered, corresponding to slow, intermediate and fast wind-waves, respectively. For each wave condition, four Schmidt numbers are considered for the scalar transport. The presence of progressive surface waves is found to induce significant wave-phase-correlated variation to the scalar field, with the phase dependence varying with the wave age. The time- and plane-averaged profiles of the scalar over waves of various ages exhibit similar vertical structures as those found in turbulence over a flat wall, but with the von Kármán constant and effective wave surface roughness for the mean scalar profile exhibiting considerable variation with the wave age. The profiles of the root-mean-square scalar fluctuations and the horizontal scalar flux exhibit good scaling in the viscous sublayer that agrees with the scaling laws previously reported for flat-wall turbulence, but with noticeable wave-induced variation in the viscous wall region. The profiles of the vertical scalar flux in the viscous sublayer exhibit apparent discrepancies from the reported scaling law for flat-wall turbulence, due to a negative vertical flux region above the windward face of the wave crest. Direct observation and quadrant-based conditional averages indicate that the wave-dependent distributions of the scalar fluctuations and fluxes are highly correlated with the coherent vortical structures in the turbulence, which exhibit clear wave-dependent characteristics in terms of both shape and preferential location.

A solute is released from the surface of a sphere flowing freely in a cylindrical channel mimicking a modern drug delivery agent in a blood vessel. The solute then disperses by the combined action of advection and diffusion. We consider reflecting boundary conditions on the sphere and absorbing boundary conditions on the channel surface mimicking a biochemical reaction between the drug and endothelial cells on the vessel surface. The drug is released either instantaneously or continuously in time. The two key observables are the mean residence time in the flow before the drug is absorbed and the width over which it is spread on the vessel surface upon reaction. We numerically solve the Fokker–Planck equation for the time-dependent substance concentration combined with an analytical solution of the flow field. As expected, we find that the presence of the sphere leads to a substantial reduction in mean residence time and reaction width. Surprisingly, however, even in the limit of very large Péclet numbers (high velocities) the sphere-free case is not generally recovered. This observation can be attributed mainly to the small, but non-negligible radial flow component induced by the moving sphere. We further identify a strong influence of the release position which sharply separates two qualitatively different regimes. If the release position is between $\unicode[STIX]{x1D703}_{0}=0$ (front) and a critical $\unicode[STIX]{x1D703}_{c}$ the substance is quickly advected away from the sphere and its overall behaviour is similar to free diffusion in an empty channel. For release between $\unicode[STIX]{x1D703}_{c}$ and $\unicode[STIX]{x1D703}_{0}=\unicode[STIX]{x03C0}$ (tail), on the other hand, the substance is pushed towards the sphere leading to behaviour reminiscent of confined diffusion between two infinitely long cylinders. The critical position $\unicode[STIX]{x1D703}_{c}$ is generally smaller than $\unicode[STIX]{x03C0}/2$ which would correspond to an equatorial release position.

Motivated by the reopening mechanics of strongly collapsed airways, we study the steady propagation of an air finger through a collapsed oil-filled channel with a single compliant wall. In a previous study using fully compliant elastic tubes, a ‘pointed’ air finger was found to propagate at high speed and low pressure, which, if clinically accessible, offers the potential for rapid reopening of highly collapsed airways with minimal tissue damage (Heap & Juel Phys. Fluids, vol. 20 (8), 2008, 081702). The mechanism underlying the selection of that pointed finger, however, remained unexplained. In this paper, we identify the required selection mechanism by conducting an experimental study in a simpler geometry: a rigid rectangular Hele-Shaw channel with an elastic top boundary. The constitutive behaviour of this elasto-rigid channel is nonlinear and broadly similar to that of an elastic tube, but unlike the tube, the channel’s cross-section adopts self-similar shapes from the undeformed state to the point of first near wall contact. The ensuing simplification of the vessel geometry enables the systematic investigation of the reopening dynamics in terms of increasing initial collapse. We find that for low levels of initial collapse, a single centred symmetric finger propagates in the channel and its shape is set by the tip curvature. As the level of collapse increases, the channel cross-section develops a central region of near opposite wall contact, and the finger shape evolves smoothly towards a ‘flat-tipped’ finger whose geometry is set by the strong depth gradient near the channel walls. We show that the flat-tipped mode of reopening is analogous to the pointed finger observed in tubes. Its propagation is sustained by the vessel’s extreme cross-sectional profile at high collapse, while vessel compliance is necessary to stabilise it. A simple scaling argument based on the dissipated power reveals that this reopening mode is preferred at higher propagation speeds when it becomes favourable to displace the elastic channel wall rather than the viscous fluid.

A methodology is proposed here to estimate the stability characteristics of bluff-body wakes using local analysis under the assumption of weakly non-parallel flows. In this connection, a generalisation of the classic spatio-temporal stability analysis for fully three-dimensional flows is first described. Secondly, an additional higher-order correction term with respect to the common saddle-point global frequency estimation is included in the analysis. The proposed method is first validated for the case of the flow past a circular cylinder and then applied to two fully three-dimensional flows: the boundary layer flow over a wall-mounted hemispherical body and the wake flow past a fixed sphere. In all the cases considered, both the estimated unstable eigenvalue and the spatial shape of the associated eigenmode are determined by local stability analysis, and they are compared with the reference counterparts obtained at a definitely higher computational cost by a fully three-dimensional global stability analysis. It is shown that the results of local stability analysis, when the higher-order correction term is included, are in excellent agreement with those obtained by global stability analysis. It is also shown that the retained correction term is of crucial importance in this perspective, leading to a remarkable improvement in accuracy with respect to the classical saddle-point estimation.

The effect of large-scale motions (LSMs) on the turbulent/non-turbulent (T/NT) interface is examined in a turbulent boundary layer. Using flow fields from direct numerical simulation, the shape of the interface and near-interface statistics are evaluated conditional on the position of the LSM. The T/NT interface is identified using the vorticity magnitude and a streak detection algorithm is adopted to identify and track the LSMs. Two-point correlation and spectral analysis of variations in the interface height show that the spatial undulation of the interface is longer in streamwise wavelength than the boundary-layer thickness, and grows with the Reynolds number in a similar manner to the LSMs. The average variation in the interface height was evaluated conditional on the position of the LSMs. The result provides statistical evidence that the interface is locally modulated by the LSMs in both the streamwise and spanwise directions. The modulation is different when the coherent structure is high- versus low-speed motion: high-speed structures lead to a wedge-shaped deformation of the T/NT interface, which causes an anti-correlation between the angles of the interface and the internal shear layer. On the other hand, low-speed structures are correlated with crests in the interface. Finally, the sudden changes in turbulence statistics across the interface are in line with the changes in the population of low-speed structures, which consist of slower mean streamwise velocity and stronger turbulence than the high-speed counterparts.

The modulation of decaying isotropic turbulence by 45 000 spherical particles of Kolmogorov-length-scale size is studied using direct particle–fluid simulations, i.e. the flow field over each particle is fully resolved by direct numerical simulations of the conservation equations. A Cartesian cut-cell method is used by which the exchange of momentum and energy at the fluid–particle interfaces is strictly conserved. It is shown that the particles absorb energy from the large scales of the carrier flow while the small-scale turbulent motion is determined by the inertial particle dynamics. Whereas the viscous dissipation rate of the bulk flow is attenuated, the particles locally increase the level of dissipation due to the intense strain rate generated near the particle surfaces due to the crossing-trajectory effect. Analogously, the rotational motion of the particles decouples from the local fluid vorticity and strain-rate field at increasing particle inertia. The high level of dissipation is partially compensated by the transfer of momentum to the fluid via forces acting at the particle surfaces. The spectral analysis of the kinetic energy budget is supported by the average flow pattern about the particles showing a nearly universal strain-rate distribution. An analytical expression for the instantaneous rate of viscous dissipation induced by each particle is derived and subsequently verified numerically. Using this equation, the local balance of fluid kinetic energy around a particle of arbitrary shape can be precisely determined. It follows that two-way coupled point-particle models implicitly account for the particle-induced dissipation rate via the momentum-coupling terms; however, they disregard the actual length scales of the interaction. Finally, an analysis of the small-scale flow topology shows that the strength of vortex stretching in the bulk flow is mitigated due to the presence of the particles. This effect is associated with the energy conversion at small wavenumbers and the reduced level of dissipation at intermediate wavenumbers. Consequently, it damps the spectral flux of energy to the small scales.

A theoretical model is proposed in this work to study the scattering of sound waves from a serrated flat plat in the presence of a uniform flow, which is of both scientific significance and practical importance. The key contribution is the analytic and rigorous description of the scattering from the laterally periodic serrations by incorporating Fourier series expansions and the Wiener–Hopf method, which collectively give a closed-form analytical solution. To validate and demonstrate the model, a couple of test cases with some representative sinusoidal- and sawtooth-shaped serrations are studied by comparing with a commercial finite element solver. The comparisons show quite good agreement for various set-ups. The subsequent parametric studies further demonstrate the efficiency of the model and the effect of the serrations for noise control. Overall, the proposed theoretical model should be able to assist in studies of low-noise aerofoils and the silent flying capabilities of owls.

The problem of orthogonal penetration of a blade into the core of a vortex with non-zero axial flow was studied using a combination of scaling theory, a heuristic plug-flow model and full Navier–Stokes simulations. The particular focus of this paper was to understand the mechanics of the transient lift force that occurs during the initial penetration of the blade leading edge into the vortex core, and the relationship of this transient force to the steady-state lift force that develops due to the difference in vortex core radius over the blade surface. The three modelling approaches all lead to the conclusion that the maximum value of the lift coefficient for the transient blade penetration force is proportional to the impact parameter and inversely proportional to the axial flow parameter. This observation is used to develop a simple expression that collapses the predictions of the full Navier–Stokes simulations for lift coefficient over a wide range of parameter values.

Wave patterns in the Faraday instability have been studied for decades. Besides the rich wave dynamics observed at the interface, Faraday waves hide elusive flow patterns in the bulk – streaming patterns – which have not been studied experimentally. The streaming patterns are responsible for a net circulation in the flow, which is reminiscent of the circulation in convection cells. In this article, we analyse these streaming flows by conducting experiments in a Faraday-wave set-up using particle image velocimetry. To visualise the flows, we perform stroboscopic measurements to both generate trajectory maps and probe the streaming velocity field. We identify three types of patterns and experimentally show that identical Faraday waves can mask streaming patterns that are qualitatively very different. Next, we consider a three-dimensional model for streaming flows in quasi-inviscid fluids, whose key is the complex coupling occurring at all of the viscous boundary layers. This coupling yields modified boundary conditions in a three-dimensional Navier–Stokes formulation of the streaming flow. Numerical simulations based on this framework show reasonably good agreement, both qualitative and quantitative, with the velocity fields of our experiments. The model highlights the relevance of three-dimensional effects in the streaming patterns. Our simulations also reveal that the variety of streaming patterns is deeply linked to the boundary condition at the top interface, which may be strongly affected by the presence of contaminants.

A theoretical and numerical study of yield-stress fluid creeping flow about a particle is presented. Yield-stress fluids can hold rigid particles statically buoyant if the yield stress is large enough. In addressing sedimentation of rigid particles in viscoplastic fluids, we should know this critical ‘yield number’ beyond which there is no motion. As we get close to this limit, the role of viscosity becomes negligible in comparison to the plastic contribution in the leading order, since we are approaching the zero-shear-rate limit. Admissible stress fields in this limit can be found by using the characteristics of the governing equations of perfect plasticity (i.e. the sliplines). This approach yields a lower bound of the critical plastic drag force or equivalently the critical yield number. Admissible velocity fields also can be postulated to calculate the upper bound. This analysis methodology is examined for three families of particle shapes (ellipse, rectangle and diamond) over a wide range of aspect ratios. Numerical experiments of either resistance or mobility problems in a viscoplastic fluid validate the predictions of slipline theory and reveal interesting aspects of the flow in the yield limit (e.g. viscoplastic boundary layers). We also investigate in detail the cases of high and low aspect ratio of the particles.

The dynamics of a jet of an inviscid incompressible liquid spiralling out under the action of centrifugal forces is considered with both gravity and the surface tension taken into account. This problem is of direct relevance to a number of industrial applications, ranging from the spinning disc atomization process to nanofibre formation. The mathematical description of the flow by necessity requires the use of a local curvilinear non-orthogonal coordinate system centred around the jet’s baseline, and we present the general formulation of the problem without assuming that the jet is slender. To circumvent the inconvenience inherent in the non-orthogonality of the local coordinate system, the orthonormal Frenet basis is used in parallel with the local non-orthogonal basis, and the equation of motion, with the velocity considered with respect to the local coordinate system, is projected onto the Frenet basis. The variation of the latter along the baseline is then described by the Frenet equations which naturally brings the baseline’s curvature and torsion into the equations of motion. This technique allows one to handle different line-based non-orthogonal curvilinear coordinate systems in a straightforward and mathematically transparent way. An analysis of the slender-jet approximation that follows the general formulation shows how a set of ordinary differential equations describing the jet’s trajectory can be derived in two cases: $\mathit{We}=O(1)$ and $\unicode[STIX]{x1D716}\mathit{We}=O(1)$ as $\unicode[STIX]{x1D716}\rightarrow 0$, where $\unicode[STIX]{x1D716}$ is the ratio of characteristic length scales across and along the jet and $\mathit{We}$ is the Weber number. A one-dimensional model for the propagation of nonlinear peristaltic disturbances along the jet is derived in each of these cases. A critical review of the work published on this topic is presented showing where errors typically occur and how to identify and avoid them.

We develop a phase-field multiphase lattice Boltzmann model to systematically investigate the dynamic behaviour of a droplet passing through a microfluidic T-junction, especially focusing on the non-breakup of the droplet. Detailed information on the breakup and non-breakup is presented, together with the quantitative evolutions of driving and resistance forces as well as the droplet deformation characteristics involved. Through comparisons between cases of non-breakup and breakup, we find that the appearance of tunnels (the lubricating film between droplet and channel walls) provides a precondition for the final non-breakup of droplets, which slows down the droplet deformation rate and even induces non-breakup. The vortex flow formed inside droplets plays an important role in determining whether they break up or not. In particular, when the strength of vortex flow exceeds a critical value, a droplet can no longer break up. Additionally, more effort has been devoted to investigating the effects of viscosity ratio between disperse and continuous phases and width ratio between branch and main channels on droplet dynamic behaviours. It is found that a large droplet viscosity results in a small velocity gradient in a droplet, which restricts vortex generation and thus produces lower deformation resistance. Consequently, it is easier to break up a droplet with larger viscosity. Our work also reveals that a droplet in small branch channels tends to obstruct the channels and have small vortex flows, which induces easier breakup too. Eventually, several phase diagrams for droplet flow patterns are provided, and the corresponding power-law correlations ($l_{0}/w=\unicode[STIX]{x1D6FD}Ca^{b}$, where $l_{0}/w$ is dimensionless initial droplet length and $Ca$ is capillary number) are fitted to describe the boundaries between different flow patterns.

The generation of a temperature disturbance in a flow is accompanied by the production of acoustic waves (direct noise) and of an entropy perturbation. If this entropy perturbation is accelerated or decelerated (as is the case through a nozzle or flow restriction), additional acoustic waves are generated (indirect noise). Several studies have demonstrated this mechanism in controlled conditions by using entropy wave generators, in which entropy waves are generated and convected through a nozzle, leading to direct and indirect noise. An analytical analysis of the direct and indirect noise produced by the generation and acceleration of entropy waves in a reflective environment is presented. The effect of reverberation (repeated acoustic reflections) on low-frequency perturbations (characteristic of entropy wave generators) is determined analytically. These results are then implemented in a set of limit cases, showing the limit behaviours of such systems. The analytical model is applied to the case of the Cambridge entropy wave generator experiment, in which entropy waves are generated by an electric heater and accelerated through a subsonic orifice plate. Due to the clear time separation of direct and indirect noise in the experimental results, direct and indirect noise transfer functions can be extracted from the experimental data for the first time and compared directly with existing theoretical models. The backward-propagating indirect noise generated at an orifice plate is shown to be significantly higher than predicted by existing theoretical models for isentropic nozzles.

An experiment is performed using high-speed photography to elucidate the behaviours of jets formed by the interactions of two laser-induced tandem bubbles produced axisymmetrically for a range of dimensionless interaction parameters such as the bubble size ratio, $\unicode[STIX]{x1D709}$, the distance between the two cavitation bubbles, $l_{0}^{\ast }$, and the time difference in bubble generation, $\unicode[STIX]{x0394}\unicode[STIX]{x1D703}^{\ast }$. A strong interaction occurs for $l_{0}^{\ast }<1$. The first bubble produced (bubble A) deforms because of the rapid growth of the second bubble (bubble B) to create a pulsed conical jet, sometimes with spray formation at the tip, formed by the small amount of water confined between the two bubbles. This phenomenon is followed by bubble penetration, toroidal bubble collapse, and the subsequent fast contraction of bubble B accompanied by a fine jet. The formation mechanism of the conical jet is similar to that of a water spike developed in air from a deformed free surface of a single growing bubble; however, the pressures of the gases surrounding each type of jet differ. The jet behaviours can be controlled by manipulating the interaction parameters; the jet velocity is significantly affected by $\unicode[STIX]{x1D709}$ and $l_{0}^{\ast }$, but less so by $\unicode[STIX]{x0394}\unicode[STIX]{x1D703}^{\ast }$ for $\unicode[STIX]{x0394}\unicode[STIX]{x1D703}^{\ast }>\unicode[STIX]{x0394}\unicode[STIX]{x1D703}_{c}^{\ast }$ ($\unicode[STIX]{x0394}\unicode[STIX]{x1D703}_{c}^{\ast }$ being the critical birth-time difference). The optimum time of jet impact, at which bubble A reaches its maximum volume, depends on $\unicode[STIX]{x0394}\unicode[STIX]{x1D703}^{\ast }$. It is generally later for larger values of $\unicode[STIX]{x1D709}$. A pulsed jet could be used to create small pores in a cell membrane; therefore, the reported method may be useful for application in tandem-bubble sonoporation.

A new multiscale simulation methodology is introduced to facilitate computationally efficient simulations of high Reynolds number turbulence seen in wall-bounded flows. The scale splitting methodology uses traditional large eddy simulation (LES) with a wall model to simulate the larger scales which are subsequently enriched using a space–time compatible kinematic simulation. Computational feasibility and robustness of the methodology are investigated using two idealized problems that emulate turbulence within the planetary boundary layer (PBL), and a finite Reynolds number channel flow problem which serves to validate the methodology against direct numerical simulation. The space–time correlations and spectra generated using enriched LES show excellent agreement with LES conducted at high resolution for all three problems; thereby demonstrating the potential of this approach for high resolution PBL simulations with a drastic reduction in the computational costs when compared to the conventional approach.

We study numerically the inertial migration of a single rigid sphere and an oblate spheroid in straight square and rectangular ducts. A highly accurate interface-resolved numerical algorithm is employed to analyse the entire migration dynamics of the oblate particle and compare it with that of the sphere. Similarly to the inertial focusing of spheres, the oblate particle reaches one of the four face-centred equilibrium positions, however they are vertically aligned with the axis of symmetry in the spanwise direction. In addition, the lateral trajectories of spheres and oblates collapse into an equilibrium manifold before ending at the equilibrium positions, with the equilibrium manifold tangential to lines of constant background shear for both sphere and oblate particles. The differences between the migration of the oblate and sphere are also presented, in particular the oblate may focus on the diagonal symmetry line of the duct cross-section, close to one of the corners, if its diameter is larger than a certain threshold. Moreover, we show that the final orientation and rotation of the oblate exhibit chaotic behaviour for Reynolds numbers beyond a critical value. Finally, we document that the lateral motion of the oblate particle is less uniform than that of the spherical particle due to its evident tumbling motion throughout the migration. In a square duct, the strong tumbling motion of the oblate in the first stage of the migration results in a lower lateral velocity and consequently longer focusing length with respect to that of the spherical particle. The opposite is true in a rectangular duct where the higher lateral velocity of the oblate in the second stage of the migration, with negligible tumbling, gives rise to shorter focusing lengths. These results can help the design of microfluidic systems for bioapplications.

Whereas surface tension decreases linearly with temperature for most fluids, here we consider those that exhibit a well-defined minimum. Specifically, our study is motivated by dilute aqueous mixtures of long-chain alcohols, for which surface tension is typically assumed to be a quadratic function of temperature. Utilization of these so-called ‘self-rewetting fluids’ has grown significantly in the past decade, due to observations that heat transfer is enhanced in applications such as film boiling and pulsating heat pipes. With similar applications in mind, we investigate the dynamics of a thin film with quadratic surface tension which is subjected to a temperature modulation in the bounding gas. A model is developed within the framework of the long-wave approximation, and a time-averaged thermocapillary driving force for destabilization is uncovered that results from the nonlinear surface tension. Linear stability analysis of the nonlinear partial differential equation for the film thickness is used to determine the critical conditions at which this driving force destabilizes the film and numerical investigation of this evolution equation reveals that linearly unstable perturbations saturate to regular periodic solutions (when the modulational frequency is set properly). Properties of these flows such as bifurcation at critical points and long-domain flows, where multiple unstable linear modes interact, are also discussed.

We consider the influence of a smooth three-dimensional (3-D) indentation on the instability of an incompressible boundary layer by linear and nonlinear analyses. The numerical work was complemented by an experimental study to investigate indentations of approximately $11\unicode[STIX]{x1D6FF}_{99}$ and $22\unicode[STIX]{x1D6FF}_{99}$ width at depths of 45 %, 52 % and 60 % of $\unicode[STIX]{x1D6FF}_{99}$, where $\unicode[STIX]{x1D6FF}_{99}$ indicates 99% boundary layer thickness. For these indentations a separation bubble confined within the indentation arises. Upstream of the indentation, spanwise-uniform Tollmien–Schlichting (TS) waves are assumed to exist, with the objective to investigate how the 3-D surface indentation modifies the 2-D TS disturbance. Numerical corroboration against experimental data reveals good quantitative agreement. Comparing the structure of the 3-D separation bubble to that created by a purely 2-D indentation, there are a number of topological changes particularly in the case of the widest indentation; more rapid amplification and modification of the upstream TS waves along the symmetry plane of the indentation is observed. For the shortest indentations, beyond a certain depth there are then no distinct topological changes of the separation bubbles and hence on flow instability. The destabilising mechanism is found to be due to the confined separation bubble and is attributed to the inflectional instability of the separated shear layer. Finally for the widest width indentation investigated ($22\unicode[STIX]{x1D6FF}_{99}$), results of the linear analysis are compared with direct numerical simulations. A comparison with the traditional criteria of using $N$-factors to assess instability of properly 3-D disturbances reveals that a general indication of flow destabilisation and development of strongly nonlinear behaviour is indicated as $N=6$ values are attained. However $N$-factors, based on linear models, can only be used to provide indications and severity of the destabilisation, since the process of disturbance breakdown to turbulence is inherently nonlinear and dependent on the magnitude and scope of the initial forcing.