Interfacial boundary conditions determined from empirical or ad hoc models remain the standard approach to model fluid flows over porous media, even in situations where the topology of the porous medium is known. We propose a non-empirical and accurate method to compute the effective boundary conditions at the interface between a porous surface and an overlying flow. Using a multiscale expansion (homogenization) approach, we derive a tensorial generalized version of the empirical condition suggested by Beavers & Joseph (J. Fluid Mech., vol. 30 (01), 1967, pp. 197-207). The components of the tensors determining the effective slip velocity at the interface are obtained by solving a set of Stokes equations in a small computational domain near the interface containing both free flow and porous medium. Using the lid-driven cavity flow with a porous bed, we demonstrate that the derived boundary condition is accurate and robust by comparing an effective model to direct numerical simulations. Finally, we provide an open source code that solves the microscale problems and computes the velocity boundary condition without free parameters over any porous bed.

The surface texture of materials plays a critical role in wettability, turbulence and transport phenomena. In order to design surfaces for these applications, it is desirable to characterise non-smooth and porous materials by their ability to exchange mass and momentum with flowing fluids. While the underlying physics of the tangential (slip) velocity at a fluid-solid interface is well understood, the importance and treatment of normal (transpiration) velocity and normal stress is unclear. We show that, when slip velocity varies at an interface above the texture, a non-zero transpiration velocity arises from mass conservation. The ability of a given surface texture to accommodate for a normal velocity of this kind is quantified by a transpiration length. We further demonstrate that normal momentum transfer gives rise to a pressure jump. For a porous material, the pressure jump can be characterised by so called resistance coefficients. By solving five Stokes problems, the introduced measures of slip, transpiration and resistance can be determined for any anisotropic non-smooth surface consisting of regularly repeating geometric patterns. The proposed conditions are a subset of effective boundary conditions derived from formal multi-scale expansion. We validate and demonstrate the physical significance of the effective conditions on two canonical problems -a lid-driven cavity and a turbulent channel flow, both with non-smooth bottom surfaces. † Email address for correspondence: ugis@mech.kth.se arXiv:1812.09401v2 [physics.flu-dyn]

Dispersion of low-density rigid particles with complex geometries is ubiquitous in both natural and industrial environments. We show that while explicit methods for coupling the incompressible Navier-Stokes equations and Newton's equations of motion are often sufficient to solve for the motion of cylindrical particles with low density ratios, for more complex particles -such as a body with a protrusion -they become unstable. We present an implicit formulation of the coupling between rigid body dynamics and fluid dynamics within the framework of the immersed boundary projection method. Similarly to previous work on this method, the resulting matrix equation in the present approach is solved using a block-LU decomposition. Each step of the block-LU decomposition is modified to incorporate the rigid body dynamics. We show that our method achieves second-order accuracy in space and first-order in time (third-order for practical settings), only with a small additional computational cost to the original method. Our implicit coupling yields stable solution for density ratios as low as 10 −4 . We also consider the influence of fictitious fluid located inside the rigid bodies on the accuracy and stability of our method.

Plants and animals use plumes, barbs, tails, feathers, hairs and fins to aid locomotion. Many of these appendages are not actively controlled, instead they have to interact passively with the surrounding fluid to generate motion. Here, we use theory, experiments and numerical simulations to show that an object with a protrusion in a separated flow drifts sideways by exploiting a symmetry-breaking instability similar to the instability of an inverted pendulum. Our model explains why the straight position of an appendage in a fluid flow is unstable and how it stabilizes either to the left or right of the incoming flow direction. It is plausible that organisms with appendages in a separated flow use this newly discovered mechanism for locomotion; examples include the drift of plumed seeds without wind and the passive reorientation of motile animals.

Despite the ubiquity of fluid flows interacting with porous and elastic materials, we lack a validated non-empirical macroscale method for characterizing the flow over and through a poroelastic medium. We propose a computational tool to describe such configurations by deriving and validating a continuum model for the poroelastic bed and its interface with the above free fluid. We show that, using stress continuity condition and slip velocity condition at the interface, the effective model captures the effects of small changes in the microstructure anisotropy correctly and predicts the overall behaviour in a physically consistent and controllable manner. Moreover, we show that the performance of the effective model is accurate by validating with fully microscopic resolved simulations. The proposed computational tool can be used in investigations in a wide range of fields, including mechanical engineering, bio-engineering and geophysics.