Over the past decade, the edge of chaos has proven to be a fruitful starting point for investigations of shear flows when the laminar base flow is linearly stable. Numerous computational studies of shear flows demonstrated the existence of states that separate laminar and turbulent regions of the state space. In addition, some studies determined invariant solutions that reside on this edge. In this paper, we study the unstable manifold of one such solution with the aid of continuous symmetry reduction, which we formulate here for the simultaneous quotiening of axial and azimuthal symmetries. Upon our investigation of the unstable manifold, we discover a previously unknown traveling-wave solution on the laminar-turbulent boundary with a relatively complex structure. By means of low-dimensional projections, we visualize different dynamical paths that connect these solutions to the turbulence. Our numerical experiments demonstrate that the laminar-turbulent boundary exhibits qualitatively different regions whose properties are influenced by the nearby invariant solutions.

The paper is devoted to the numerical analysis of wave patterns behind a shock wave propagating in a tube filled with a gaseous mixture. It is shown that the flow inside the boundary layer behind the shock wave is unstable, and the way the instability develops fully corresponds to the solution obtained for the boundary layer over a flat plate. Vortical perturbations inside the boundary layer determine the nonuniformity of the temperature field. In turn, exactly these nonuniformities define the way the ignition kernels arise in the combustible mixture after the reflected shock interaction with the boundary layer. In particular, the temperature nonuniformity determines the spatial limitations of probable ignition kernel position relative to the end wall and side walls of the tube. In the case of low-intensity incident shocks the ignition could start not farther than the point of first interaction between the reflected shock wave and roller vortices formed in the process of boundary layer development. Proposed physical mechanisms are formulated in general terms and can be used for interpretation of the experimental data in any systems with a delayed exothermal reaction start. It is also shown that contact surface thickening occurs due to its interaction with Tollmien-Schlichting waves. This conclusion is of importance for understanding the features of ignition in shock tubes operating in the over-tailored regime.

A regular hexagonal pattern of three-dimensional electroconvective flow induced by unipolar injection in dielectric liquids is numerically observed by solving the fully coupled governing equations using the lattice Boltzmann method. A small-amplitude perturbation in the form of a spatially periodic pattern of hexagonal cells is introduced initially. The transient development of convective cells that undergo a sequence of transitions agrees with the idea of flow seeking an optimal scale. Stable hexagonal convective cells and their subcritical bifurcation together with a hysteresis loop are clearly observed. In addition, the stability of the hexagonal flow pattern is analyzed in a wide range of relevant parameters, including the electric Rayleigh number $T$, nondimensional mobility $M$, and wave number $k$. It is found that centrally downflowing hexagonal cells, which are characterized by the central region being empty of charge, are preferred in the system.

Transition to turbulence in incompressible adverse-pressure-gradient (APG) boundary layers is investigated by direct numerical simulations. Purely two-dimensional instability waves develop on the inflectional base velocity profile. When the boundary layer is perturbed by isotropic turbulence from the free stream, streamwise elongated streaks form and may interact with the instability waves. Subsequent mechanisms that trigger transition depend on the intensity of the free-stream disturbances. All evidence from the present simulations suggest that the growth rate of instability waves is sufficiently high to couple with the streaks. Under very low levels of free-stream turbulence $\left(\sim 0.1\%\right)$, transition onset is highly sensitive to the inlet disturbance spectrum and is accelerated if the spectrum contains frequency–wave-number combinations that are commensurate with the instability waves. Transition onset and completion in this regime is characterized by formation and breakdown of $\mathrm{\Lambda }$ vortices, but they are more sporadic than in natural transition. Beneath free-stream turbulence with higher intensity (1–$2\%)$, bypass transition mechanisms are dominant, but instability waves are still the most dominant disturbances in wall-normal and spanwise perturbation spectra. Most of the breakdowns were by disturbances with critical layers close to the wall, corresponding to inner modes. On the other hand, the propensity of an outer mode to occur increases with the free-stream turbulence level. Higher intensity free-stream disturbances induce strong streaks that favorably distort the boundary layer and suppress the growth of instability waves. But the upward displacement of high amplitude streaks brings them to the outer edge of the boundary layer and exposes them to ambient turbulence. Consequently, high-amplitude streaks exhibit an outer-mode secondary instability.

Anisotropy in the turbulent stress tensor, which forms the basis of invariant analysis, is conducted using velocity time series measurements collected in the canopy sublayer (CSL) and the atmospheric surface layer (ASL). The goal is to assess how thermal stratification and surface roughness conditions simultaneously distort the scalewise relaxation towards isotropic state from large to small scales when referenced to homogeneous turbulence. To achieve this goal, conventional invariant analysis is extended to allow scalewise information about relaxation to isotropy in physical (instead of Fourier) space to be incorporated. The proposed analysis shows that the CSL is more isotropic than its ASL counterpart at large, intermediate, and small (or inertial) scales irrespective of the thermal stratification. Moreover, the small (or inertial) scale anisotropy is more prevalent in the ASL when compared to the CSL, a finding that cannot be fully explained by the intensity of the mean velocity gradient acting on all scales. Implications to the validity of scalewise Rotta and Lumley models for return to isotropy as well as advantages to using barycentric instead of anisotropy invariant maps for such scalewise analysis are discussed.

This paper addresses reconnection of vortex tubes, with particular focus on the topology of the vortex lines (field lines of the vorticity). This analysis of vortex line topology reveals key features of the reconnection process, such as the generation of many small flux rings, formed when reconnection occurs in multiple locations in the vortex sheet between the tubes. Consideration of three-dimensional reconnection principles leads to a robust measurement of the reconnection rate, even once instabilities break the symmetry. It also allows us to identify internal reconnection of vortex lines within the individual vortex tubes. Finally, the introduction of a third vortex tube is shown to render the vortex reconnection process fully three-dimensional, leading to a fundamental change in the topological structure of the process. An additional interesting feature is the generation of vorticity null points.

Vapor bubbles are formed in liquids by two mechanisms: evaporation (temperature above the boiling threshold) and cavitation (pressure below the vapor pressure). The liquid resists in these metastable (overheating and tensile, respectively) states for a long time since bubble nucleation is an activated process that needs to surmount the free energy barrier separating the liquid and the vapor states. The bubble nucleation rate is difficult to assess and, typically, only for extremely small systems treated at an atomistic level of detail. In this work a powerful approach, based on a continuum diffuse interface modeling of the two-phase fluid embedded with thermal fluctuations (fluctuating hydrodynamics), is exploited to study the nucleation process in homogeneous conditions, evaluating the bubble nucleation rates and following the long-term dynamics of the metastable system, up to the bubble coalescence and expansion stages. In comparison with more classical approaches, this methodology allows us on the one hand to deal with much larger systems observed for a much longer time than possible with even the most advanced atomistic models. On the other, it extends continuum formulations to thermally activated processes, impossible to deal with in a purely determinist setting.

The flow behind a model wind turbine under two different turbine operating regimes (region 2 for turbine operating at optimal condition with the maximum power coefficient and 1.4-deg pitch angle and region 3 for turbine operating at suboptimal condition with a lower power coefficient and 7-deg pitch angle) is investigated using wind tunnel experiments and numerical experiments using large-eddy simulation (LES) with actuator surface models for turbine blades and nacelle. Measurements from the model wind turbine experiment reveal that the power coefficient and turbine wake are affected by the operating regime. Simulations with and without a nacelle model are carried out for each operating condition to study the influence of the operating regime and nacelle on the formation of the hub vortex and wake meandering. Statistics and energy spectra of the simulated wakes are in good agreement with the measurements. For simulations with a nacelle model, the mean flow field is composed of an outer wake, caused by energy extraction by turbine blades, and an inner wake directly behind the nacelle, while for the simulations without a nacelle model, the central region of the wake is occupied by a jet. The simulations with the nacelle model reveal an unstable helical hub vortex expanding outward toward the outer wake, while the simulations without a nacelle model show a stable and columnar hub vortex. Because of the different interactions of the inner region of the wake with the outer region of the wake, a region with higher turbulence intensity is observed in the tip shear layer for the simulation with a nacelle model. The hub vortex for the turbine operating in region 3 remains in a tight helical spiral and intercepts the outer wake a few diameters further downstream than for the turbine operating in region 2. Wake meandering, a low-frequency large-scale motion of the wake, commences in the region of high turbulence intensity for all simulations with and without a nacelle model, indicating that neither a nacelle model nor an unstable hub vortex is a necessary requirement for the existence of wake meandering. However, further analysis of the wake meandering and instantaneous flow field using a filtering technique and dynamic mode decomposition show that the unstable hub vortex energizes the wake meandering. The turbine operating regime affects the shape and expansion of the hub vortex, altering the location of the onset of the wake meandering and wake meander oscillating intensity. Most important, the unstable hub vortex promotes a high-amplitude energetic meandering which cannot be predicted without a nacelle model.

Melting and, conversely, solidification processes in the presence of convection are key to many geophysical problems. An essential question related to these phenomena concerns the estimation of the (time-evolving) melting rate, which is tightly connected to the turbulent convective dynamics in the bulk of the melt fluid and the heat transfer at the liquid-solid interface. In this work, we consider a convective-melting model, constructed as a generalization of the Rayleigh-Bénard system, accounting for the basal melting of a solid. As the change of phase proceeds, a fluid layer grows at the heated bottom of the system and eventually reaches a turbulent convection state. By means of extensive lattice-Boltzmann numerical simulations employing an enthalpy formulation of the governing equations, we explore the model dynamics in two- and three-dimensional configurations. The focus of the analysis is on the scaling of global quantities like the heat flux and the kinetic energy with the Rayleigh number, as well as on the interface morphology and the effects of space dimensionality. Independently of dimensionality, we find that the convective-melting system behavior shares strong resemblances with that of the Rayleigh-Bénard one, and that the heat flux is only weakly enhanced with respect to that case. Such similarities are understood, at least to some extent, considering the resulting slow motion of the melting front (with respect to the turbulent fluid velocity fluctuations) and its generally little roughness (compared to the height of the fluid layer). Varying the Stefan number, accounting for the thermodynamical properties of the material, also seems to have only a mild effect, which implies the possibility of extrapolating results in numerically delicate low-Stefan setups from more convenient high-Stefan ones. Finally, we discuss the implications of our findings for the geophysically relevant problem of modeling Arctic ice melt ponds.

The growth of Tollmien-Schlichting (TS) waves is experimentally attenuated using a single-input and single-output (SISO) feedback system, where the TS wave packet is generated by a surface point source in a flat-plate boundary layer. The SISO system consists of a single wall-mounted hot wire as the sensor and a miniature speaker as the actuator. The actuation is achieved through a dual-slot geometry to minimize the cavity near-field effects on the sensor. The experimental setup to generate TS waves or wave packets is very similar to that used by Li and Gaster [J. Fluid Mech. 550, 185 (2006)]. The aim is to investigate the performance of the SISO control system in attenuating single-frequency, two-dimensional disturbances generated by these configurations. The necessary plant models are obtained using system identification, and the controllers are then designed based on the models and implemented in real-time to test their performance. Cancellation of the rms streamwise velocity fluctuation of TS waves is evident over a significant domain.

We report on experiments that consist of deforming a collection of monodisperse droplets produced by a microfluidic chip through a flow-focusing device. We show that a proper numerical modeling of the flow is necessary to access the stress applied by the latter on the droplet along its trajectory through the chip. This crucial step enables the full integration of the differential equation governing the dynamical deformation, and consequently the robust measurement of the interfacial tension by fitting the experiments with the calculated deformation. Our study thus demonstrates the feasibility of quantitative in situ rheology in microfluidic flows involving, e.g., droplets, capsules, or cells.

Concentration gradients in a fluid adjacent to a reactive surface due to contrast in surface reactivity generate convective flows. These flows result from contributions by electro- and diffusio-osmotic phenomena. In this study, we have analyzed reactive patterns that release and consume protons, analogous to bimetallic catalytic conversion of peroxide. Similar systems have typically been studied using either scaling analysis to predict trends or costly numerical simulation. Here, we present a simple analytical model, bridging the gap in quantitative understanding between scaling relations and simulations, to predict the induced potentials and consequent velocities in such systems without the use of any fitting parameters. Our model is tested against direct numerical solutions to the coupled Poisson, Nernst-Planck, and Stokes equations. Predicted slip velocities from the model and simulations agree to within a factor of $\approx 2$ over a multiple order-of-magnitude change in the input parameters. Our analysis can be used to predict enhancement of mass transport and the resulting impact on overall catalytic conversion, and is also applicable to predicting the speed of catalytic nanomotors.

In hydrodynamic topological transitions, one mass of fluid breaks into two or two merge into one. For example, in honey-drop formation when honey is dripping from a spoon, honey is extended to separate into two masses as the liquid neck bridging them thins down to the micron scale. At the moment when the topology changes due to the breakup, physical observables such as surface curvature locally diverge. Such singular dynamics has widely attracted physicists, revealing universality in self-similar dynamics, which shares much in common with critical phenomena in thermodynamics. Many experimental examples have been found, including an electric spout and vibration-induced jet eruption. However, only a few cases have been physically understood on the basis of equations that govern the singular dynamics and even in such a case the physical understanding is mathematically complicated, inevitably involving delicate numerical calculations. Here we study the breakup of air film entrained by a solid disk into viscous liquid in a confined space, which leads to formation, thinning, and breakup of the neck of air. As a result, we unexpectedly find that equations governing the neck dynamics can be solved analytically by virtue of two remarkable experimental features: Only a single length scale linearly dependent on time remains near the singularity and two universal scaling functions describing the singular neck shape and velocity field are both analytic. The present solvable case would be essential for a better understanding of the singular dynamics and will help reveal the physics of unresolved examples intimately related to daily-life phenomena and diverse practical applications.

Motivated by the emerging applications of liquid-infused surfaces (LIS), we study the drag reduction and robustness of transverse flows over two-dimensional microcavities partially filled with an oily lubricant. Using separate simulations at different scales, characteristic contact line velocities at the fluid-solid intersection are first extracted from nanoscale phase field simulations and then applied to micronscale two-phase flows, thus introducing a multiscale numerical framework to model the interface displacement and deformation within the cavities. As we explore the various effects of the lubricant-to-outer-fluid viscosity ratio ${\tilde{\mu }}_2/{\tilde{\mu }}_1$, the capillary number Ca, the static contact angle ${\theta }_s$, and the filling fraction of the cavity $\delta$, we find that the effective slip is most sensitive to the parameter $\delta$. The effects of ${\tilde{\mu }}_2/{\tilde{\mu }}_1$ and ${\theta }_s$ are generally intertwined but weakened if $\delta <1$. Moreover, for an initial filling fraction $\delta =0.94$, our results show that the effective slip is nearly independent of the capillary number when it is small. Further increasing Ca to about $0.01{\tilde{\mu }}_1/{\tilde{\mu }}_2$, we identify a possible failure mode, associated with lubricants draining from the LIS, for ${\tilde{\mu }}_2/{\tilde{\mu }}_1\lesssim 0.1$. Very viscous lubricants (e.g., ${\tilde{\mu }}_2/{\tilde{\mu }}_1>1$), however, are immune to such failure due to their generally larger contact line velocity.

Direct numerical simulations of turbulent duct flows with aspect ratios 1, 3, 5, 7, 10, and 14.4 at a center-plane friction Reynolds number ${\text{Re}}_{\tau ,c}\simeq 180$, and aspect ratios 1 and 3 at ${\text{Re}}_{\tau ,c}\simeq 360$, were carried out with the spectral-element code nek5000. The aim of these simulations is to gain insight into the kinematics and dynamics of Prandtl's secondary flow of the second kind and its impact on the flow physics of wall-bounded turbulence. The secondary flow is characterized in terms of the cross-plane component of the mean kinetic energy, and its variation in the spanwise direction of the flow. Our results show that averaging times of around 3000 convective time units (based on duct half-height $h$) are required to reach a converged state of the secondary flow, which extends up to a spanwise distance of around $\simeq 5h$ measured from the side walls. We also show that if the duct is not wide enough to accommodate the whole extent of the secondary flow, then its structure is modified as reflected through a different spanwise distribution of energy. Another confirmation of the extent of the secondary flow is the decay rate of kinetic energy of any remnant secondary motions for $z_c/h>5$ (where $z_c$ is the spanwise distance from the corner) in aspect ratios 7, 10, and 14.4, which exhibits a decreasing level of energy with increasing averaging time $t_a$, and in its rapid rate of decay given by $\sim t_a^{-1}$. This is the same rate of decay observed in a spanwise-periodic channel simulation, which suggests that at the core, the kinetic energy of the secondary flow integrated over the cross-sectional area, ${\langle K\rangle }_{yz}$, behaves as a random variable with zero mean, with rate of decay consistent with central limit theorem. Long-time averages of statistics in a region of rectangular ducts extending about the width of a well-designed channel simulation (i.e., extending about $\simeq 3h$ on each side of the center plane) indicate that ducts or experimental facilities with aspect ratios larger than 10 may, if properly designed, exhibit good agreement with results obtained from spanwise-periodic channel computations.

Nanomechanical resonators and sensors, operated in ambient conditions, often generate low-Mach-number oscillating rarefied gas flows. Cercignani [C. Cercignani, J. Stat. Phys. 1, 297 (1969)] proposed a variational principle for the linearized Boltzmann equation, which can be used to derive approximate analytical solutions of steady (time-independent) flows. Here we extend and generalize this principle to unsteady oscillatory rarefied flows and thus accommodate resonating nanomechanical devices. This includes a mathematical approach that facilitates its general use and allows for systematic improvements in accuracy. This formulation is demonstrated for two canonical flow problems: oscillatory Couette flow and Stokes' second problem. Approximate analytical formulas giving the bulk velocity and shear stress, valid for arbitrary oscillation frequency, are obtained for Couette flow. For Stokes' second problem, a simple system of ordinary differential equations is derived which may be solved to obtain the desired flow fields. Using this framework, a simple and accurate formula is provided for the shear stress at the oscillating boundary, again for arbitrary frequency, which may prove useful in application. These solutions are easily implemented on any symbolic or numerical package, such as Mathematica or matlab, facilitating the characterization of flows produced by nanomechanical devices and providing insight into the underlying flow physics.

Motivated by the recent discovery of subsurface oceans on planetary moons and the interest they have generated, we explore convective flows in shallow spherical shells of dimensionless gap width ${\varepsilon }^2\ll 1$ in the rapid rotation limit $\mathrm{E}\ll 1$, where $\mathrm{E}$ is the Ekman number. We employ direct numerical simulation (DNS) of the Boussinesq equations to compute the local heat flux $\mathrm{Nu}\left(\lambda \right)$ as a function of the latitude $\lambda$ and use the results to characterize the trapping of convection at low latitudes, around the equator. We show that these results are quantitatively reproduced by an asymptotically exact nonhydrostatic equatorial $\beta$-plane convection model at a much more modest computational cost than DNS. We identify the trapping parameter $\beta =\varepsilon {\mathrm{E}}^{-1}$ as the key parameter that controls the vigor and latitudinal extent of convection for moderate thermal forcing when $\mathrm{E}\sim \varepsilon$ and $\varepsilon \downarrow 0$. This model provides a theoretical paradigm for nonlinear investigations.

Eigenspectra and pseudospectra of the mean-linearized Navier-Stokes operator are used to characterize amplification mechanisms in laminar and turbulent flows in which linear mechanisms are important. Success of mean flow (linear) stability analysis for a particular frequency is shown to depend on whether two scalar measures of non-normality agree: (1) the product between the resolvent norm and the distance from the imaginary axis to the closest eigenvalue and (2) the inverse of the inner product between the most amplified resolvent forcing and response modes. If they agree, the resolvent operator can be rewritten in its dyadic representation to reveal that the adjoint and forward stability modes are proportional to the forcing and response resolvent modes at that frequency. Hence the real parts of the eigenvalues are important since they are responsible for resonant amplification and the resolvent operator is low rank when the eigenvalues are sufficiently separated in the spectrum. If the amplification is pseudoresonant, then resolvent analysis is more suitable to understand the origin of observed flow structures. Two test cases are studied: low Reynolds number cylinder flow and turbulent channel flow. The first deals mainly with resonant mechanisms, hence the success of both classical and mean stability analysis with respect to predicting the critical Reynolds number and global frequency of the saturated flow. Both scalar measures of non-normality agree for the base and mean flows, and the region where the forcing and response modes overlap scales with the length of the recirculation bubble. In the case of turbulent channel flow, structures result from both resonant and pseudoresonant mechanisms, suggesting that both are necessary elements to sustain turbulence. Mean shear is exploited most efficiently by stationary disturbances while bounds on the pseudospectra illustrate how pseudoresonance is responsible for the most amplified disturbances at spatial wavenumbers and temporal frequencies corresponding to well-known turbulent structures. Some implications for flow control are discussed.

The classical solution in the pressure-impulse theory for the inviscid, incompressible, and zero-surface-tension water impact of a flat plate at zero dead-rise angle exhibits both singular-in-time initial fluid acceleration, ${\left.\partial \mathbf{v}/\partial t\right|}_{t=0}\sim \delta \left(t\right)$, and a near-plate-edge spatial singularity in the velocity distribution, $\mathbf{v}\sim r^{-1/2}$, where $r$ is the distance from the plate edge. The latter velocity divergence also leads to the interface being stretched infinitely right after the impact, which is another nonphysical artifact. From the point of view of matched asymptotic analysis, this classical solution is a singular limit when three physical quantities achieve limiting values: sound speed $c_0\to \infty$, fluid kinematic viscosity $\nu \to 0$, and surface tension $\sigma \to 0$. This leaves open a question on how to resolve these singularities mathematically by including the neglected physical effects—compressibility, viscosity, and surface tension—first one by one and then culminating in the local compressible viscous solution valid for $t\to 0$ and $r\to 0$, demonstrating a nontrivial flow structure that changes with the degree of the bulk compressibility. In the course of this study, by starting with the general physically relevant formulation of compressible viscous flow, we clarify the parameter range(s) of validity of the key analytical solutions including classical ones (inviscid incompressible and compressible, etc.) and understand the solution structure, its intermediate asymptotics nature, characteristics influencing physical processes, and the role of potential and rotational flow components. In particular, it is pointed out that sufficiently close to the plate edge surface tension must be taken into account. Overall, the idea is to highlight the interesting physics behind the singularities in the pressure-impulse theory.