We present an improved method for computing incompressible viscous flow around suspended rigid particles using a fixed and uniform computational grid. The main idea is to incorporate PeskinÕs regularized delta function approach [Acta Numerica 11 (2002) 1] into a direct formulation of the fluid-solid interaction force in order to allow for a smooth transfer between Eulerian and Lagrangian representations while at the same time avoiding strong restrictions of the time step. This technique was implemented in a finite-difference and fractional-step context. A variety of two-and three-dimensional simulations are presented, ranging from the flow around a single cylinder to the sedimentation of 1000 spherical particles. The accuracy and efficiency of the current method are clearly demonstrated.

Recent remarkable progress in computing power and numerical analysis is enabling us to fill a gap in the dynamical systems approach to turbulence. One of the significant advances in this respect has been the numerical discovery of simple invariant sets, such as nonlinear equilibria and periodic solutions, in well-resolved Navier--Stokes flows. This review describes some fundamental and practical aspects of dynamical systems theory for the investigation of turbulence, focusing on recently found invariant solutions and their significance for the dynamical and statistical characterization of low-Reynolds-number turbulent flows. It is shown that the near-wall regeneration cycle of coherent structures can be reproduced by such solutions. The typical similarity laws of turbulence, i.e. the Prandtl wall law and the Kolmogorov law for the viscous range, as well as the pattern and intensity of turbulence-driven secondary flow in a square duct can also be represented by these simple invariant solutions.Comment: To appear in Annual Review of Fluid Mechanics, Vol. 44, 201

The behaviour of turbulent shear flow over a mass-neutral permeable wall is studied numerically. The transpiration is assumed to be proportional to the local pressure fluctuations. It is first shown that the friction coefficient increases by up to 40% over passively porous walls, even for relatively small porosities. This is associated with the presence of large spanwise rollers, originating from a linear instability which is related both to the Kelvin–Helmholtz instability of shear layers, and to the neutral inviscid shear waves of the mean turbulent profile. It is shown that the rollers can be forced by patterned active transpiration through the wall, also leading to a large increase in friction when the phase velocity of the forcing resonates with the linear eigenfunctions mentioned above. Phase-lock averaging of the forced solutions is used to further clarify the flow mechanism. This study is motivated by the control of separation in boundary layers.

We have conducted a direct numerical simulation (DNS) study of dilute turbulent particulate flow in a vertical plane channel, considering thousands of finite-size rigid particles with resolved phase interfaces. The particle diameter corresponds to approximately 11 wall units and their terminal Reynolds number is set to 136. The fluid flow with bulk Reynolds number 2700 is directed upward, which maintains the particles suspended upon average. Two density ratios were simulated, differing by a factor of 4.5. The corresponding Stokes numbers of the two flow cases were O(10) in the near-wall region and O(1) in the outer flow. We have observed the formation of large-scale elongated streak-like structures with streamwise dimensions of the order of 8 channel half-widths and cross-stream dimensions of the order of one half-width. At the same time, we have found no evidence of significant formation of particle clusters, which suggests that the large structures are due to an intrinsic instability of the flow, triggered by the presence of the particles. It was found that the mean fluid velocity profile tends towards a concave shape, and the turbulence intensity as well as the normal stress anisotropy are strongly increased. The effect of varying the Stokes number while maintaining the buoyancy, particle size and volume fraction constant was relatively weak.