The effects of changing the trailing edge shape on the wake and propulsive performance of a pitching rigid panel are examined experimentally. The panel aspect ratio is $\text{AR}=1$, and the trailing edges are symmetric chevron shapes with convex and concave orientations of varying degree. Concave trailing edges delay the natural vortex bending and compression of the wake, and the mean streamwise velocity field contains a single jet. Conversely, convex trailing edges promote wake compression and produce a quadfurcated wake with four jets. As the trailing edge shape changes from the most concave to the most convex, the thrust and efficiency increase significantly.

The modal and nonmodal linear stability of the flow in a microchannel with either one or both walls coated with a superhydrophobic material is studied. The topography of the bounding wall(s) has the shape of elongated microridges with arbitrary alignment with respect to the direction of the mean pressure gradient. The superhydrophobic walls are modeled using the Navier slip condition through a slip tensor, and the results depend parametrically on the slip length and orientation angle of the ridges. The stability analysis is carried out in the temporal framework; the modal analysis is performed by solving a generalized eigenvalue problem, and the nonmodal, optimal perturbation analysis is done with an adjoint optimization approach. We show theoretically and verify numerically that Squire's theorem does not apply in the present settings, despite the fact that Squire modes are found to be always damped. The most notable result is the appearance of a streamwise wall-vortex mode at very low Reynolds numbers when the ridges are sufficiently inclined with respect to the mean pressure gradient, in the case of a single superhydrophobic wall. When two walls are rendered water repellent, the exponential growth of the instability results from either a two-dimensional or a three-dimensional Orr-Sommerfeld mode, depending on the ridges' orientation and amplitude. Nonmodal results for either one or two superhydrophobic wall(s) display but a mild modification of the no-slip case.

We quantify the diffusiophoresis and electrophoresis of a uniformly charged, spherical colloid in a binary electrolyte using modified Poisson-Nernst-Planck equations that account for steric repulsion between finite sized ions. Specifically, we utilize the Bikerman (Bik) lattice gas model and the Carnahan-Starling (CS) and Boublik-Mansoori-Carnahan-Starling-Leland (BMCSL) equations of state for monodisperse and polydisperse, respectively, hard spheres. We compute the phoretic mobility for weak applied fields using an asymptotic approach for thin diffuse layers, where ion steric effects are expected to be most prevalent. The thin diffuse layer limit requires ${\lambda }_D/R\to 0$, where ${\lambda }_D$ is the Debye screening length and $R$ is the particle radius; this limit is readily attained for micron-sized colloids in concentrated electrolytic solutions. It is well known that the classic Poisson-Boltzmann (PB) model for pointlike, noninteracting ions leads to a prediction of a maximum in both the diffusiophoretic and electrophoretic mobilities with increasing particle zeta potential (at fixed ${\lambda }_D/R$). In contrast, we find that ion sterics essentially eliminate this maximum (for reasonably attainable zeta potentials) and increase the mobility relative to PB. Next, we consider the more experimentally relevant case of a particle with a constant surface charge density and vary the electrolyte concentration, neglecting charge regulation on surface active sites. Rather surprisingly, there is little difference between the predictions of the four models (PB, Bik, CS, and BMCSL) for electrophoretic mobility in concentrated solutions, at reasonable surface charge densities ($\sim 1–10\mu \mathrm{C}/{\mathrm{cm}}^2$). This is because as the concentration increases, the zeta potential is reduced (to below the thermal voltage for concentrations above about 1 M) and therefore the diffuse layer structure is largely unaffected by ion sterics. For gradients of symmetric electrolytes (equal diffusivities, charge, and size) diffusiophoresis is also essentially unaffected by ion sterics, with a mobility that approaches zero with increasing concentration, just as in electrophoresis. For gradients of asymmetric electrolytes, the difference in diffusivities of the cation and anions leads to an induced electric field that acts on the charged particle. Importantly, we show that ion sterics leads to an excess contribution to the induced electric field, which increases rapidly with concentration. This increase overwhelms the accompanying decrease in zeta potential. The result is the diffusiophoretic mobility increases with concentration, rather than approaching zero. Therefore, diffusiophoresis could be an appealing alternative transport mechanism to electrophoresis in concentrated electrolyte solutions.

We investigate the large-scale transport of inertial particles. We derive explicit analytic expressions for the eddy diffusivities for generic Stokes times. These latter expressions are exact for any shear flow while they correspond to the leading contribution either in the deviation from the shear flow geometry or in the Péclet number of general random Gaussian velocity fields. Our explicit expressions allow us to investigate the role of inertia for such a class of flows and to make exact links with the analogous transport problem for tracer particles.

Topological features of vortices in terms of local flow geometry in an isotropic homogeneous turbulence and relationships between the topology and development of a vortex are investigated. Swirlity and sourcity represent the unidirectionality and intensities of the azimuthal and radial flows, respectively. Skewness reflects the symmetry quantity of the vortical flow. Together, the three quantities characterize the details of the invariant flow geometry. The analysis shows that flow symmetry is correlated with swirlity, and associated with its development or decay, especially for weak vortices. Stronger vortices exhibit lower correlations between the two and lack sufficient flow symmetry for inflow in all directions. This clarifies why highly intense vortices cannot attain effective vortex stretching. Vortex stretching is formulated in terms of the radial flow in the swirl plane. It shows that stretching is not entirely characterized by the eigenvalues of the rate-of-strain tensor, and depends on the eigenvalues of the radial flow classified by sourcity. Sourcity shows that numerous vortices, which have been classified as inflow (convergent) vortices using the complex eigenvalues of the velocity gradient tensor, have a mixture of inflow and outflow. These vortices break the orthogonality of the vortical axis to the swirl plane during vortex stretching. Stretching in vortices with complete inflow increases vorticity normal to the swirl plane effectively and improves axis orthogonality. Sourcity classifies these characteristics, and flow symmetry and sourcity are important quantities for effective stretching. The present topological analysis provides important details of vortical flow features in turbulent flows or realistic vortex models, and sourcity extracts the specific vortical region supported by the flow geometry itself.

Extensional rheometry and Brownian dynamics simulations of flexible polymer solutions confirm predictions based on blob concepts that coil-stretch hysteresis in extensional flows increases with concentration, reaching a maximum at the critical overlap concentration $c^{*}$ before progressively vanishing in the semidilute regime. These observations demonstrate that chain stretching strengthens intermolecular hydrodynamic screening in dilute solutions, but weakens it in semidilute solutions. Flow can thus strongly modify the concentration dependence of viscoelastic properties of polymer solutions.

We investigate the dynamics of two superposed layers with density contrast flowing countercurrent inside a channel, when the lower layer is much thinner than the wavelength of interfacial waves. We apply a low-dimensional film model to the bottom (heavier) layer and introduce a fast and efficient method to predict the onset of flow reversal in this phase. We study three vertical scenarios with different applied pressure gradients and compare the temporal growth rates of linear and weakly nonlinear waves to the Orr-Sommerfeld problem and to the weakly nonlinear theory, respectively. At the loading point, i.e., when a large wave hump stands at the interface, our spatiotemporal analysis shows that the system is absolutely unstable. We then present profiles of nonlinear saturated waves, pressure field, and streamline distribution in agreement with direct numerical simulation. The reduced model presented here allows us to explore the effect of the upper-layer speed on the wave pattern, showing that the wave profile is very sensitive when the mean film thickness, rather than the liquid flow rate, is maintained constant in the simulation. In addition, we show the strong effect of surface tension on both the maximum wave hump and the crest steepness before the loading point. Finally, we reveal how the nonlinear wave speed affects the vortex distribution within the lower layer by analyzing the stream function under different scenarios.

Proper orthogonal decomposition (POD) is applied to distinct data sets in order to characterize the propagation of error arising from basis truncation in the description of turbulence. Experimental data from stereo particle image velocimetry measurements in a wind turbine array and direct numerical simulation data from a fully developed channel flow are used to illustrate dependence of the anisotropy tensor invariants as a function of POD modes used in low-order descriptions. In all cases, ensembles of snapshots illuminate a variety of anisotropic states of turbulence. In the near wake of a model wind turbine, the turbulence field reflects the periodic interaction between the incoming flow and rotor blade. The far wake of the wind turbine is more homogenous, confirmed by the increased magnitude of the anisotropy factor. By contrast, the channel flow exhibits many anisotropic states of turbulence. In the inner layer of the wall-bounded region, one observes one-component turbulence at the wall; immediately above, the turbulence is dominated by two components, with the outer layer showing fully three-dimensional turbulence, conforming to theory for wall-bounded turbulence. The complexity of flow descriptions resulting from truncated POD bases can be greatly mitigated by severe basis truncations. However, the current work demonstrates that such simplification necessarily exaggerates the anisotropy of the modeled flow and, in extreme cases, can lead to the loss of three-dimensionality. Application of simple corrections to the low-order descriptions of the Reynolds stress tensor significantly reduces the residual root-mean-square error. Similar error reduction is seen in the anisotropy tensor invariants. Corrections of this form reintroduce three-dimensionality to severe truncations of POD bases. A threshold for truncating the POD basis based on the equivalent anisotropy factor for each measurement set required many more modes than a threshold based on energy. The mode requirement to reach the anisotropy threshold after correction is reduced by a full order of magnitude for all example data sets, ensuring that economical low-dimensional models account for the isotropic quality of the turbulence field.

We experimentally investigate the propulsive characteristics of a pitching NACA 0012 airfoil section, with emphasis on thrust and propulsive efficiency, at a Reynolds number of $1.7\times {10}^4$. For the sake of mechanical simplicity, we consider an airfoil restricted to a single actuator in the pitching direction. We examine the effect of changing the airfoil's axis of rotation, finding that contrary to Garrick's linear theory, there exists a pitching axis near the airfoil that maximizes propulsive efficiency. Next, we examine the effect of placing passive springs on the airfoil in the heave (transverse) direction using our Cyber-Physical Fluid Dynamics technique. This elastic heaving motion allows the airfoil to combine pitching and heaving modes while being actuated only in the pitching direction. Two sets of dynamics are considered: one case where the airfoil is weighted unevenly and pitched about its center of mass (so that the resulting heaving motion is independent of inertial forces), and another case where the airfoil's center of mass is fixed at its centroid. For pitching at an amplitude of $8^{\circ }$ and a reduced frequency $k$ of two, we find that elastic heave produces a maximum propulsive efficiency of 35%, compared to 25% without any heave motion. Further, while operating at the same efficiency as the static-pivot case, we find that passive heaving greatly increases the magnitude of the airfoil's thrust. The airfoil configurations with highest propulsive efficiency generally involve pitching near or ahead of the airfoil's leading edge.

The linear response of the wall-shear stress of a turbulent flow to wall waviness is analyzed in the context of a comparison between existing experiments, direct numerical simulations, and analytical approximations. The spectral region where the response is largest is found to be amenable to a simplified quasilaminar analysis. The end result is a parameterless description of this phenomenon that completely captures its physics in a single analytical formula, a Padé approximation of the response function.

We consider the stability of a binary liquid film flowing down a heated incline. A theoretical model is implemented which captures the Soret effect and the dependence of surface tension on both temperature and solutal concentration. The model also allows for variation in the density of the liquid mixture with thermal and solutal differences. A linear stability analysis is performed with asymptotic and numerical results being obtained. The coupling of the effect of a variable density with the thermosolutal-Marangoni instability and the Soret effect is investigated. Good agreement with previous results for the constant density case is found.

A weighted residual approximation is used to derive a reduced model to analyze the stability of a binary liquid film flowing down a heated incline discussed in Part I. Nonlinear simulations are carried out and reveal a permanent wave structure arising in unstable flow configurations. In addition, a weakly nonlinear analysis is performed and yields a Benney-type equation which has been solved both analytically and numerically. Good agreement between the linear stability analysis, weighted residual model, and weakly nonlinear analysis is found.

Direct numerical simulations are carried out to investigate the effect of the wall temperature on the behavior of oblique shock wave turbulent boundary layer interactions at free-stream Mach number 2.28 and shock angle of the wedge generator $\phi =8^{\circ }$. Five values of the wall-to-recovery-temperature ratio $\left(T_w/T_r\right)$ are considered, corresponding to cold, adiabatic, and hot wall thermal conditions. We show that the main effect of cooling is to decrease the characteristic scales of the interaction in terms of upstream influence and extent of the separation bubble. The opposite behavior is observed in the case of heating, which produces a marked dilatation of the interaction region. The distribution of the Stanton number shows that a strong amplification of the heat transfer occurs across the interaction, with the maximum thermal and dynamic loads found for the case of the cold wall. The analysis reveals that the fluctuating heat flux exhibits a strong intermittent behavior, characterized by scattered spots with extremely high values compared to the mean. Furthermore, the analogy between momentum and heat transfer, typical of compressible, wall-bounded, equilibrium turbulent flows, does not apply for most of the interaction domain. The premultiplied spectra of the wall heat flux do not show any evidence of the influence of the low-frequency shock motion, and the primary mechanism for the generation of peak heating is found to be linked with the turbulence amplification in the interaction region.

The present study attempts to change the regime transitions and heat transfer properties in rotating Rayleigh-Bénard convection by injecting $\sim 100-\mu \mathrm{m}$-diam particles in the flow. The particles start settling out of the fluid immediately after injection and separate entirely from the fluid over a period of several hours. The particles deposit on the top and bottom surfaces, forming porous layers with nonideal thermal properties, and, as expected, decrease the heat flux. The reduction in heat transfer is a result of the inability of the layers to respond rapidly enough to fluid temperature fluctuations. However, in the rotation-dominated geostrophic regime, the heat transfer normalized by its nonrotating value is higher in the presence of the particle layers than without them. Direct numerical simulations with ideal heat transfer walls indicate that the temperature fluctuations in the bulk become slower under the damping effect of rotation, in contrast with those in the boundary layers, which do not show any damping until the flow transitions to the geostrophic regime. In this regime, the dominant time scale of the near-wall fluid temperature fluctuations increases substantially as the rotation rate is increased. It is thus likely that the time response of the particle layers in relation to that of the nearby fluid improves only in the geostrophic regime, which is reflected in a relatively larger heat transfer.

The boundary layer structure of the velocity and temperature fields in turbulent Rayleigh-Bénard flows in closed cylindrical cells of unit aspect ratio is revisited from a transitional and turbulent viscous boundary layer perspective. When the Rayleigh number is large enough, the dynamics at the bottom and top plates can be separated into an impact region of downwelling plumes, an ejection region of upwelling plumes, and an interior region away from the side walls. The latter is dominated by the shear of the large-scale circulation (LSC) roll, which fills the whole cell and continuously varies its orientation. The working fluid is liquid mercury or gallium at a Prandtl number $\text{Pr}=0.021$ for Rayleigh numbers $3\times {10}^5\le \text{Ra}\le 4\times {10}^8$. The generated turbulent momentum transfer corresponds to macroscopic flow Reynolds numbers with $1.8\times {10}^3\le \text{Re}\le 4.6\times {10}^4$. In highly resolved spectral element direct numerical simulations, we present the mean profiles of velocity, Reynolds stress, and temperature in inner viscous units and compare our findings with convection experiments and channel flow data. The complex three-dimensional and time-dependent structure of the LSC in the cell is compensated by a plane-by-plane symmetry transformation which aligns the horizontal velocity components and all its derivatives with the instantaneous orientation of the LSC. As a consequence, the torsion of the LSC is removed, and a streamwise direction in the shear flow can be defined. It is shown that the viscous boundary layers for the largest Rayleigh numbers are highly transitional and obey properties that are directly comparable to transitional channel flows at friction Reynolds numbers ${\text{Re}}_{\tau }\lesssim {10}^2$. The transitional character of the viscous boundary layer is also underlined by the strong enhancement of the fluctuations of the wall stress components with increasing Rayleigh number. An extrapolation of our analysis data suggests that the friction Reynolds number ${\text{Re}}_{\tau }$ in the velocity boundary layer can reach values of 200 for $\text{Ra}\gtrsim {10}^{11}$. Thus the viscous boundary layer in a liquid metal flow would become turbulent at a much lower Rayleigh number than for turbulent convection in gases and gas mixtures.

We report details of the fluid motion formed within and above grooves when a laminar continuous phase fluid flows over a second immiscible fluid enclosed in a grooved microchannel. Vortical structures within the transverse grooves were caused by a slip velocity at the fluid-fluid interface and act as “fluid bearings” on the boundary to lubricate the flow of the continuous phase. We investigated the drag reduction in the laminar flow in the microchannel by measuring slip at the boundaries and calculating an effective slip length, taking into account the influence of the effect of the viscosity ratio of the two fluids on the effective slip length. The “fluid bearing” effect can be used to transport high viscosity fluids using low viscosity fluids trapped in cavities to reduce drag.

Singular solutions to the Euler equations could provide essential insight into the formation of very small scales in highly turbulent flows. Previous attempts to find singular flow structures have proven inconclusive. We reconsider the problem of interacting vortex tubes, for which it has long been observed that the flattening of the vortices inhibits sustained self-amplification of velocity gradients. Here we consider an iterative mechanism, based on the transformation of vortex filaments into sheets and their subsequent instability back into filaments. Elementary fluid mechanical arguments are provided to support the formation of a singular structure via this iterated mechanism, which we analyze based on a simplified model of filament interactions.

The present work investigates the effectiveness of the control strategy introduced by Schoppa and Hussain [Phys. Fluids 10, 1049 (1998)] as a function of Reynolds number ($\text{Re}$). The skin-friction drag reduction method proposed by these authors, consisting of streamwise-invariant, counter-rotating vortices, was analyzed by Canton et al. [Flow, Turbul. Combust. 97, 811 (2016)] in turbulent channel flows for friction Reynolds numbers (${\text{Re}}_{\tau }$) corresponding to the value of the original study (i.e., 104) and 180. For these $\text{Re}$, a slightly modified version of the method proved to be successful and was capable of providing a drag reduction of up to 18%. The present study analyzes the Reynolds number dependence of this drag-reducing strategy by performing two sets of direct numerical simulations (DNS) for ${\text{Re}}_{\tau }=360$ and 550. A detailed analysis of the method as a function of the control parameters (amplitude and wavelength) and $\text{Re}$ confirms, on the one hand, the effectiveness of the large-scale vortices at low $\text{Re}$ and, on the other hand, the decreasing and finally vanishing effectiveness of this method for higher $\text{Re}$. In particular, no drag reduction can be achieved for ${\text{Re}}_{\tau }=550$ for any combination of the parameters controlling the vortices. For low Reynolds numbers, the large-scale vortices are able to affect the near-wall cycle and alter the wall-shear-stress distribution to cause an overall drag reduction effect, in accordance with most control strategies. For higher $\text{Re}$, instead, the present method fails to penetrate the near-wall region and cannot induce the spanwise velocity variation observed in other more established control strategies, which focus on the near-wall cycle. Despite the negative outcome, the present results demonstrate the shortcomings of the control strategy and show that future focus should be on methods that directly target the near-wall region or other suitable alternatives.

In this paper, we present results from an experimental study into turbulent Rayleigh-Bénard convection forced externally by periodically modulated unidirectional rotation rates. We find that the azimuthal rotation velocity $\dot{\theta }\left(t\right)$ and thermal amplitude $\delta \left(t\right)$ of the large-scale circulation (LSC) are modulated by the forcing, exhibiting a variety of dynamics including increasing phase delays and a resonant peak in the amplitude of $\dot{\theta }\left(t\right)$. We also focus on the influence of modulated rotation rates on the frequency of occurrence $\eta$ of stochastic cessation or reorientation events, and on the interplay between such events and the periodically modulated response of $\dot{\theta }\left(t\right)$. Here we identify a mechanism by which $\eta$ can be amplified by the modulated response, and these normally stochastic events can occur with high regularity. We provide a modeling framework that explains the observed amplitude and phase responses, and we extend this approach to make predictions for the occurrence of cessation events and the probability distributions of $\dot{\theta }\left(t\right)$ and $\delta \left(t\right)$ during different phases of a modulation cycle, based on an adiabatic approach that treats each phase separately. Last, we show that such periodic forcing has consequences beyond influencing LSC dynamics, by investigating how it can modify the heat transport even under conditions where the Ekman pumping effect is predominant and strong enhancement of heat transport occurs. We identify phase and amplitude responses of the heat transport, and we show how increased modulations influence the average Nusselt number.