We consider the evolution of a spatially periodic, perturbed vortex sheet for small times after the formation of a curvature singularity at time $t=t_c$ as demonstrated by Moore (Proc. R. Soc. Lond. A, vol. 365, issue 1720, 1979, pp. 105–119). The Moore analysis is extended to provide the small-amplitude, full-sheet structure at $t=t_c$ for a general single-mode initial condition in terms of polylogarithmic functions, from which its asymptotic form near the singular point is determined. This defines an intermediate evolution problem for which the leading-order, and most singular, approximation is solved as a Taylor-series expansion in $\tau = t-t_c$, where coefficients are calculated by repeated differentiation of the defining Birkhoff–Rott (BR) equation. The first few terms are in good agreement with numerical calculation based on the full-sheet solution. The series is summed, providing an analytic continuation which shows sheet rupture at circulation $\varGamma =0^+$, $\tau >0^+$, but with non-physical features owing to the absence of end-tip sheet roll up. This is corrected by constructing an inner solution with $\varGamma < \tau$, as a perturbed similarity form with small parameter $\tau ^{1/2}$. Numerical solutions of both the inner, nonlinear zeroth-order and first-order linear BR equations are obtained whose outer limits match the intermediate solution. The composite solution shows sheet tearing at $\tau =0^+$ into two separate, rolled up algebraic spirals near the central singular point. Branch separation distance scales as $\tau$ with a non-local, $\tau ^{3/2}$ correction. Properties of the intermediate and inner solutions are discussed.

Motile bacteria play essential roles in biology that rely on their dynamic behaviours, including their ability to navigate, interact and self-organize. However, bacteria dynamics on fluid interfaces are not well understood. Swimmers adsorbed on fluid interfaces remain highly motile, and fluid interfaces are highly non-ideal domains that alter swimming behaviour. To understand these effects, we study flow fields generated by Pseudomonas aeruginosa PA01 in the pusher mode. Analysis of correlated displacements of tracers and bacteria reveals dipolar flow fields with unexpected asymmetries that differ significantly from their counterparts in bulk fluids. We decompose the flow field into fundamental hydrodynamic modes for swimmers in incompressible fluid interfaces. We find an expected force-doublet mode corresponding to propulsion and drag at the interface plane, and a second dipolar mode, associated with forces exerted by the flagellum on the cell body in the aqueous phase that are countered by Marangoni stresses in the interface. The balance of these modes depends on the bacteria's trapped interfacial configurations. Understanding these flows is broadly important in nature and in the design of biomimetic swimmers.

Motivated by wind blowing over water, we use asymptotic methods to study the evolution of short wavelength interfacial waves driven by the combined action of these flows. We solve the Rayleigh equation for the stability of the shear flow, and construct a uniformly valid approximation for the perturbed streamfunction, or eigenfunction. We then expand the real part of the eigenvalue, the phase speed, in a power series of the inverse wavenumber and show that the imaginary part is exponentially small. We give expressions for the growth rates of the Miles (J. Fluid Mech., vol. 3, 1957, pp. 185–204) and rippling (e.g. Young & Wolfe, J. Fluid Mech., vol. 739, 2014, pp. 276–307) instabilities that are valid for an arbitrary shear flow. The accuracy of the results is demonstrated by a comparison with the exact solution of the eigenvalue problem in the case when both the wind and the current have an exponential profile.

A novel boundary integral approach for the recovery of overhanging (or not) rotational water waves (with constant vorticity) from pressure measurements at the bottom is presented. The method is based on the Cauchy integral formula and on an Eulerian–Lagrangian formalism to accommodate overturning free surfaces. This approach eliminates the need to introduce a priori a special basis of functions, thus providing a general means of fitting the pressure data and, consequently, recovering the free surface. The effectiveness and accuracy of the method are demonstrated through numerical examples.

In continuation of our earlier work (Chen & Sreenivasan, J. Fluid Mech., vol. 908, 2021, R3; Chen & Sreenivasan, J. Fluid Mech., vol. 933, 2022a, A20 – together referred to as CS hereafter), we present a self-consistent Reynolds number asymptotics for wall-normal profiles of variances of streamwise and spanwise velocity fluctuations as well as root-mean-square pressure, across the entire flow region of channel and pipe flows and flat-plate boundary layers. It is first shown that, when normalized by peak values, the Reynolds number dependence and wall-normal variation of all three profiles can be decoupled, in excellent agreement with available data, sharing the common inner expansion of the type $\phi (y^+)=f_0(y^+)+f_1(y^+)/Re^{1/4}_\tau$, where $\phi$ is one of the quantities just mentioned, the functions $f_0$ and $f_1$ depend only on $y^+$, and $Re_\tau$ is the friction Reynolds number. Here, the superscript $+$ indicates normalization by wall variables. We show that this result is completely consistent with CS. Secondly, by matching the above inner expansion and the outer flow similarity form, a bounded variation $\phi (y^\ast )=\alpha _\phi -\beta _{\phi }y^{{\ast {1}/{4}}}$ is derived for the outer region where, for each $\phi$, the constants $\alpha _\phi$ and $\beta _{\phi }$ are independent of $Re_\tau$ and $y^\ast$ $\equiv y^+/Re_\tau$ – also in excellent agreement with simulations and experimental data. One of the predictions of the analysis is that, for asymptotically high Reynolds numbers, a finite plateau $\phi \approx \alpha _\phi$ appears in the outer region. This result sheds light on the intriguing issue of the outer shoulder of the variance of the streamwise velocity fluctuation, which should be bounded by the asymptotic plateau of approximately 10.

We numerically study the transverse flow-induced vibration (FIV) of elastically coupled tandem cylinders at Reynolds number $100$, using an in-house immersed boundary method-based solver in two-dimensional coordinates. While several previous studies considered tandem cylinders coupled through flow between them, a hitherto unexplored elastic coupling with fluid flow between them significantly influences FIV. We consider a wide range of gap ratio, reduced velocity, an equal mass ratio of both cylinders and zero damping. A systematic comparison between the classic elastically mounted tandem cylinders and elastically coupled cylinders is presented. The latter configuration exhibits two vibration modes, in-phase and out-of-phase, with corresponding natural frequencies approaching the Strouhal frequency of the system. We quantify variation of the following output variables with reduced velocity and gap ratios: cylinders’ displacement; fluid forces; amplitude spectral density of displacement and force signals; phase characteristics; energy harvesting potential; and discuss the wake characteristics using flow separation, pressure distribution, gap flow quantification, and dynamic mode decomposition characterization. The FIV response is classified into several regimes: initial desynchronization with and without gap vortices; final desynchronization; mixed mode; initial branch; lock-in; upper and lower branch; wake-induced vibration; galloping. We draw upon similarities of computed FIV characteristics with those of an isolated cylinder, in which the lower branch exhibits larger a amplitude than the upper branch. The elastically coupled cylinders show a galloping response similar to an isolated D-section cylinder. By invoking the elastic coupling, we demonstrate FIV suppression and augmentation for in-phase and out-of-phase systems. Our calculations show larger energy harvesting potential at reduced cost for elastically coupled cylinders.

Diapycnal mixing plays a key role in the thermohaline circulation of the deep ocean. Field observations have suggested that this mixing is intensified over the rough topography along the boundaries of the ocean. In this study, we experimentally explore the transport of salinity and tracer in a horizontally stirred, stratified fluid with a steady vertical buoyancy flux. The mechanical mixing occurs either uniformly across the tank or near a sidewall. To compare uniform and boundary mixing, first we explore the steady state dynamics and find that, in both cases, the surfaces of constant density are horizontal. With uniform mixing, vertical transport of buoyancy occurs uniformly throughout the tank while, with boundary mixing, transport is confined to the turbulent region and the interior space remains quiescent and does not mix or have any net movement. We then explore the transient evolution of the stratification when the source of buoyancy is removed from the base of the system. In the boundary mixing case, the resulting divergence of the turbulent diffusive flux in the boundary region leads to a reduction in the buoyancy of the fluid in the boundary region, and a net upflow develops in the boundary region. In turn, this drives a downwelling in the interior. Vertical gradients in the rate of downwelling lead to stretching of isopycnals and, together, these processes enable the interior stratification to evolve. The experiments highlight that, with boundary mixing, the main transport of buoyancy occurs near the boundaries even though the interior is stratified; independent measurements of tracer mixing in the interior show that this fluid may be quiescent and stratified even though there is a large flux being transported in the boundary region. This has important implications for the interpretation of mixing data in the ocean. In a companion paper, (Li and Woods Part 2), we explore the interaction of a net upwelling with such boundary mixing.

Turbine–wake and farm–atmosphere interactions influence wind farm power production. For large offshore farms, the farm–atmosphere interaction is usually the more significant effect. This study proposes an analytical model of the ‘momentum availability factor’ to predict the impact of farm–atmosphere interactions. It models the effects of net advection, pressure gradient forcing and turbulent entrainment, using steady quasi-one-dimensional flow assumptions. Turbulent entrainment is modelled by assuming self-similar vertical shear stress profiles. We used the model with the ‘two-scale momentum theory’ to predict the power of large finite-sized farms. The model compared well with existing results of large-eddy simulations of finite wind farms in conventionally neutral boundary layers. The model captured most of the effects of atmospheric boundary layer (ABL) height on farm performance by considering the undisturbed vertical shear stress profile of the ABL as an input. In particular, the model predicted the power of staggered wind farms with a typical error of 5 % or less. The developed model provides a novel way of predicting instantly the power of large wind farms, including the farm blockage effects. A further simplification of the model to predict analytically the ‘wind extractability factor’ is also presented. This study provides a novel framework for modelling farm–atmosphere interactions. Future studies can use the framework to better model large wind farms.

Evaporating sessile droplets are critical to many industrial applications and are also ubiquitous in nature. Two predominant evaporation models have emerged in the literature, one-sided and diffusion-limited, with differing assumptions on the evaporation process. Both models are used widely, and their predictions can differ greatly from each other, but the physical mechanisms responsible for these differences are not yet well understood. Here, we develop a lubrication-theory-based model of a thin evaporating sessile droplet, and compare predictions from both evaporation models to elucidate the origins of the differences in their predictions. For the one-sided model, we derive expressions for the droplet lifetime, show that in certain parameter regimes the total evaporation rate is proportional to the droplet surface area, and demonstrate that the contact line is always warmer than the bulk of the droplet. Furthermore, we show that differences in the structures of the evaporation models near the contact line lead to qualitatively different behaviour of the apparent contact angles and interface temperature profiles. The fundamental understanding gained from this work is expected to be helpful in determining which evaporation model is most appropriate for describing experimental observations.

In microfluidic systems, analyte/reagent mixing is essential to achieve rapid chemical/biochemical reactions. The low Reynolds number (Re) flow, however, makes passive mixing difficult and necessitates the use of some external stimulus to cause disturbance in the system. Here, we report how a periodic thermocapillary effect may interact with strategically patterned wall wettability and boost the mixing dynamics in a micro-confined ternary-liquid film system. Approximate, yet without compromise on the physics involved, analytical solutions to the energy and Navier–Stokes equations under creeping flow conditions are obtained to comprehend the thermal and hydrodynamic characteristics of the thermo-capillarity. In the binary fluid limit, our model is validated with the findings of Pendse & Esmaeeli (Intl J. Therm. Sci., vol. 49, 2010, pp. 1147–1155), and good agreement is found. The flow characteristics in the ternary-liquid system due to discrete wall temperatures is also demonstrated via finite-element-based numerical simulations, and agreement with the semi-analytical solution is noted. The qualitative study of the flow pattern reveals that enhanced mixing is obtained as a result of vortical motion created by the interplay of the periodic thermo-capillarity-driven interfacial flow and wall slip. The observation is further consolidated from the numerical simulation of the species transport equation. The species distribution so obtained is compared with fully developed laminar flow for the same inlet concentration. Our study further investigates the effects of the fluids’ relative thermal conductivity, wall slip length, relative film thickness and thermal (or slip) phase differences on the mixing efficiency of the proposed arrangement. While wall slip, relative film thickness and relative thermal conductivity regulate mixing in the top and bottom layers, phase difference determines how well the middle fluid layer mixes. Slip at the walls creates vortex distortions to get a strong churning effect. The lower thermal conductivity of the top and bottom fluids weakens the thermocapillary flow; however, the dominance of the patterned slip improves mixing in such cases.

In statistically stationary homogeneous incompressible turbulence, the average energy transfer rate balance which exists at diffusion/dissipation-dominated length scales does not reflect what actually happens locally in space and time. We use a highly resolved direct numerical simulation of forced periodic turbulence to shed some light on the actual fluctuating dynamics which occur at these very small scales and which are rubbed off by averaging. Even though the viscous diffusion in physical space averages to zero and fluctuates less intensely than all other terms (except the energy input rate) in the local (in space–time) two-point energy balance, it fundamentally cannot be neglected. The local unsteadiness and the interspace turbulence transport terms cannot be ignored either in the interscale energy dynamics in spite of the fact that they also average to zero.

The transition to turbulence in a plane Poiseuille flow of dilute polymer solutions is studied by direct numerical simulations of a finitely extensible nonlinear elastic fluid with the Peterlin closure. The range of Reynolds number ($Re$) $2000 \le Re \le 5000$ is studied but with the same level of elasticity in viscoelastic flows. The evolution of a finite-amplitude perturbation and its effects on the transition dynamics are investigated. A viscoelastic flow begins transition at an earlier time than its Newtonian counterparts, but the transition time appears to be insensitive to polymer concentration in the dilute and semi-dilute regimes studied. Increasing polymer concentration, however, decreases the maximum attainable energy growth during the transition process. The critical or minimum perturbation amplitude required to trigger transition is computed. Interestingly, both Newtonian and viscoelastic flows follow almost the same power-law scaling of $Re^\gamma$ with the critical exponent $\gamma \approx -1.25$, which is in close agreement with previous studies. However, a shift downward is observed for viscoelastic flow, suggesting that smaller perturbation amplitudes are required for the transition. A mechanism of the early transition is investigated by the evolution of wall-normal and spanwise velocity fluctuations and flow structure. The early growth of these fluctuations and the formation of quasi-streamwise vortices around low-speed streaks are promoted by polymers, hence causing an early transition. These vortical structures are found to support the critical exponent $\gamma \approx -1.25$. Once the transition process is completed, polymers play a role in dampening the wall-normal and spanwise velocity fluctuations and vortices to attain a drag-reduced state in viscoelastic turbulent flows.

In this work, the dynamics of two-dimensional rotating Janus drops in shear flow is studied numerically using a ternary-fluid diffuse interface method. The rotation of Janus drops is found to be closely related to their deformation. A new deformation parameter $D$ is proposed to assess the significance of the drop deformation. According to the maximum value of $D$ ($D_{max}$), the deformation of rotating Janus drops can be classified into linear deformation ($D_{max}\le 0.2$) and nonlinear deformation ($D_{max}> 0.2$). In particular, $D_{max}$ in the former depends linearly on the Reynolds and capillary numbers, which can be interpreted by a mass–spring model. Furthermore, the rotation period $t_R$ of a Janus drop is found to be more sensitive to the drop deformation than to the aspect ratio of the drop at equilibrium. By introducing a corrected shear rate and an aspect ratio of drop deformation, a rotation model for Janus drops is established based on Jeffery's theory for rigid particles, and it agrees well with our numerical results.

Advective dispersion of solutes in long thin axisymmetric channels is important to the analysis and design of a wide range of devices, including chemical separation systems and microfluidic chips. Despite extensive analysis of Taylor dispersion in various scenarios, most studies focus on long-term dispersion behaviour and cannot capture the transient evolution of the solute zone across the spatial variations in the channel. In the current study, we analyse the Taylor–Aris dispersion for arbitrarily shaped axisymmetric channels. We derive an expression for solute dynamics in terms of two coupled ordinary differential equations, which allow prediction of the time evolution of the mean location and axial (standard deviation) width of the solute zone as a function of the channel geometry. We compare and benchmark our predictions with Brownian dynamics simulations for a variety of cases, including linearly expanding/converging channels and periodic channels. We also present an analytical description of the physical regimes of transient positive versus negative axial growth of solute width. Finally, to further demonstrate the utility of the analysis, we demonstrate a method to engineer channel geometries to achieve desired solute width distributions over space and time. We apply the latter analysis to generate a geometry that results in a constant axial width and a second geometry that results in a sinusoidal axial variance in space.

A new resolvent-based method is developed to predict the space–time properties of the flow field. To overcome the deterioration of the prediction accuracy with increasing distance between the measurements and predictions in the resolvent-based estimation (RBE), the newly proposed method utilizes the RBE to estimate the relative energy distribution near the wall rather than the absolute energy directly estimated from the measurements. Using this extra information from RBE, the new method modifies the energy distribution of the spatially uniform and uncorrelated forcing that drives the flow system by minimizing the norm of the cross-spectral density tensor of the error matrix in the near-wall region in comparison with the RBE-estimated one, and therefore it is named as the resolvent-informed white-noise-based estimation (RWE) method. For validation, three time-resolved direct numerical simulation (DNS) datasets with the friction Reynolds numbers $Re_\tau = 180$, 550 and 950 are generated, with various locations of measurements ranging from the near-wall region ($y^+ = 40$) to the upper bound of the logarithmic region ($y/h \approx 0.2$, where h is the half-channel height) for the predictions. Besides the RWE, three existing methods, i.e. the RBE, the $\lambda$-model and the white-noise-based estimation (WBE), are also included for the validation. The performance of the RBE and scale-dependent model ($\lambda$-model) in predicting the energy spectra shows a strong dependence on the measurement locations. The newly proposed RWE shows a low sensitivity on $Re_{\tau }$ and the measurement locations, which may range from the near-wall region to the upper bound of the logarithmic region, and has a high accuracy in predicting the energy spectra. The RWE also performs well in predicting the space–time properties in terms of the correlation magnitude and the convection velocity. We further utilize the new method to reconstruct the instantaneous large-scale structures with measurements from the logarithmic region. Both the RWE and RBE perform well in estimating the instantaneous large-scale structure, and the RWE has smaller errors in the estimations near the wall. The structural inclination angles around $15^\circ$ are predicted by the RWE and WBE, which generally recover the DNS results.

Flow around curved tandem cylinders in the convex configuration has been studied by means of direct numerical simulations, for a Reynolds number of 500 and a nominal gap ratio of 3.0. Spanwise variation of flow regimes, as well as curvature-induced axial velocity, leads to an exceedingly complex vortex dynamics in the wake. Both parallel and oblique vortex shedding are observed. Oblique shedding is connected to repeated occurrences of dislocations. The dislocations are caused by two main mechanisms: frequency differences in the upper part of the curved geometry and shedding of gap vortices into the lower near wake. Both types of dislocations are closely associated with a mode switch in the gap. In parts of the gap, there is low-frequency quasi-periodic asymmetry of the gap vortices, where the flow is biased to one side of the gap for intervals of several wake vortex shedding periods. The switch from side to side is associated with a surge of the vertical velocity, and the frequency of the switch is similar to that of long-term variation of the recirculation length in the lower gap.