Physics of Fluids publishes in traditional areas of fluid dynamics as well as in novel and emerging areas of the field including dynamics of gases, liquids, and complex or multiphase liquids. Physics of Fluids is published with the cooperation of the APS Division of Fluid Dynamics.
Editor's Picks

The RichtmyerMeshkov instability (RMI) is investigated using the Direct Simulation Monte Carlo (DSMC) method of molecular gas dynamics. Due to the inherent statistical noise and the significant computational requirements, DSMC is hardly ever applied to hydrodynamic flows. Here, DSMC RMI simulations are performed to quantify the shockdriven growth of a singlemode perturbation on the interface between two atmosphericpressure monatomic gases prior to reshocking as a function of the Atwood and Mach numbers. The DSMC results qualitatively reproduce all features of the RMI and are in reasonable quantitative agreement with existing theoretical and empirical models. Consistent with previous work in this field, the DSMC simulations indicate that RMI growth follows a universal behavior.

Lubrication analysis is used to determine analytical expressions for the elements of the resistance matrix describing the interaction of two rigid cylindrical particles in twodimensional shear flow in a symmetrically confined channel geometry. The developed model is valid for nonBrownian particles in a lowReynoldsnumber flow between two sliding plates with thin gaps between the two particles and also between the particles and the walls. Using this analytical model, a comprehensive overview of the dynamics of interacting cylindrical particles in shear flow is presented. With only hydrodynamic interactions, rigid particles undergo a reversible interaction with no crossstreamline migration, irrespective of the confinement value. However, the interaction time of the particle pair substantially increases with confinement, and at the same time, the minimum distance between the particle surfaces during the interaction substantially decreases with confinement. By combining our purely hydrodynamic model with a simple on/off nonhydrodynamic attractive particle interaction force, the effects of confinement on particle aggregation are qualitatively mapped out in an aggregation diagram. The latter shows that the range of initial relative particle positions for which aggregation occurs is increased substantially due to geometrical confinement. The interacting particle pair exhibits tangential and normal lubrication forces on the sliding plates, which will contribute to the rheology of confined suspensions in shear flow. Due to the combined effects of the confining walls and the particle interaction, the particle velocities and resulting forces both tangential and perpendicular to the walls exhibit a nonmonotonic evolution as a function of the orientation angle of the particle pair. However, by incorporating appropriate scalings of the forces, velocities, and doublet orientation angle with the minimum free fraction of the gap height and the plate speed, master curves for the forces versus orientation angle can be constructed.

Bulk acoustic wave devices are typically operated in a resonant state to achieve enhanced acoustic amplitudes and high acoustofluidic forces for the manipulation of microparticles. Among other loss mechanisms related to the structural parts of acoustofluidic devices, damping in the fluidic cavity is a crucial factor that limits the attainable acoustic amplitudes. In the analytical part of this study, we quantify all relevant loss mechanisms related to the fluid inside acoustofluidic microdevices. Subsequently, a numerical analysis of the timeharmonic viscoacoustic and thermoviscoacoustic equations is carried out to verify the analytical results for 2D and 3D examples. The damping results are fitted into the framework of classical linear acoustics to set up a numerically efficient device model. For this purpose, all damping effects are combined into an acoustofluidic loss factor. Since some components of the acoustofluidic loss factor depend on the acoustic mode shape in the fluid cavity, we propose a twostep simulation procedure. In the first step, the loss factors are deduced from the simulated mode shape. Subsequently, a second simulation is invoked, taking all losses into account. Owing to its computational efficiency, the presented numerical device model is of great relevance for the simulation of acoustofluidic particle manipulation by means of acoustic radiation forces or acoustic streaming. For the first time, accurate 3D simulations of realistic microdevices for the quantitative prediction of pressure amplitudes and the related acoustofluidic forces become feasible.

We propose a novel strategy for designing chaotic micromixers using curved channels confined between two flat planes. The location of the separatrix between the Dean vortices, induced by centrifugal forces, is dependent on the location of the maxima of axial velocity. An asymmetry in the axial velocity profile can change the location of the separatrix. This is achieved physically by introducing slip alternatingly at the top and bottom walls. This leads to streamline crossing and Lagrangian chaos. An approximate analytical solution of the velocity field is obtained using perturbation theory. This is used to find the Lagrangian trajectories of fluid particles. Poincare sections taken at periodic locations in the axial direction are used to study the extent of chaos. We study two microchannel designs, called circlet and serpentine, in which the Dean vortices in adjacent half cells are corotating and counterrotating, respectively. The extent of mixing, at low Re and low slip length, is shown to be greater in the serpentine case. Wide channels are observed to have much better mixing than tall channels; an important observation not made for separatrix flows till now. Eulerian indicators are used to gauge the extent of mixing, with varying slip length, and it is shown that an optimum slip length exists which maximizes the mixing in a particular geometry. Once the parameter space of relatively high mixing is identified, detailed variance computations are carried out to identify the detailed features.

Townsend’s attached eddy hypothesis forms the basis of an established model of the logarithmic layer in wallbounded turbulent flows in which this inertially dominated region is characterised by a hierarchy of geometrically selfsimilar eddying motions that scale with their distance to the wall. The hypothesis has gained considerable support from high Reynolds number measurements of the secondorder moments of the fluctuating velocities. Recently, Meneveau and Marusic [“Generalized logarithmic law for highorder moments in turbulent boundary layers,” J. Fluid Mech. 719, R1 (2013)] presented experimental evidence that all evenordered moments of the streamwise velocity will exhibit a logarithmic dependence on the distance from the wall. They demonstrated that this was consistent with the attached eddy hypothesis, so long as the velocity distribution is assumed to be Gaussian (which allows the use of the central limit theorem). In this paper, we derive this result from the attached eddy model without assuming a Gaussian velocity distribution, and find that such logarithmic behaviours are valid in the large Reynolds number limit. We also revisit the physical and mathematical basis of the attached eddy hypothesis, in order to increase rigour and minimise the assumptions required to apply the hypothesis. To this end, we have extended the proof of Campbell’s theorem to apply to the velocity field corresponding to a forest of variously sized eddies that are randomly placed on the wall. This enables us to derive all moments of the velocity in the logarithmic region, including crosscorrelations between different components of the velocity. By contrast, previous studies of the attached eddy hypothesis have considered only the mean velocity and its second order moments. From this, we obtain qualitatively correct skewnesses and flatnesses for the spanwise and wallnormal fluctuations. The issue of the Reynolds number dependence of von Kármán’s constant is also addressed.