Longitudinal to Transverse Metachronal Wave Transitions in an In Vitro Model of Ciliated Bronchial Epithelium.

Myriads of cilia beat on ciliated epithelia, which are ubiquitous in life. When ciliary beats are synchronized, metachronal waves emerge, whose direction of propagation depends on the living system in an unexplained way. We show on a reconstructed human bronchial epithelium in vitro that the direction of propagation is determined by the ability of mucus to be transported at the epithelial surface. Numerical simulations show that longitudinal waves maximize the transport of mucus while transverse waves, observed when the mucus is rigid and still, minimize the energy dissipated by the cilia.

Self-Induced Rayleigh-Taylor Instability in Segregating Dry Granular Flows.

Dry-granular material flowing on rough inclines can experience a self-induced Rayleigh-Taylor (RT) instability followed by the spontaneous emergence of convection cells. For this to happen, particles are different in size and density; the larger particles are denser but still segregate toward the surface. When the flow is initially made of two layers of particles (dense particles above), a RT instability develops during the flow. When the flow is initially made of one homogeneous layer mixture, the granular segregation leads to the formation of an unstable layer of large, dense particles at the surface, that subsequently destabilizes in a RT plume pattern. The unstable density gradient has been only induced by the motion of the granular matter. This self-induced Rayleigh-Taylor instability and the two-layer RT instability are studied using two different methods: experiments and simulations. At last, contrary to the usual fluid behavior where the RT instability relaxes into two superimposed stable layers of fluid, the granular flow evolves to a pattern of alternated bands corresponding to recirculation cells analogous to Rayleigh-Bénard convection cells where segregation sustains the convective motion.

Improved partially saturated method for the lattice Boltzmann pseudopotential multicomponent flows.

This paper extends the partially saturated method (PSM), used for curved or complex walls, to the lattice Boltzmann (LB) pseudopotential multicomponent model and adapts the wetting boundary condition to model the contact angle. The pseudopotential model is widely used for various complex flow simulations due to its simplicity. To simulate the wetting phenomenon within this model, the mesoscopic interaction force between the boundary fluid and solid nodes is used to mimic the microscopic adhesive force between the fluid and the solid wall, and the bounce-back (BB) method is normally adopted to achieve the no-slip boundary condition. In this paper, the pseudopotential interaction forces are computed with eighth-order isotropy since fourth-order isotropy leads to the condensation of the dissolved component on curved walls. Due to the staircase approximation of curved walls in the BB method, the contact angle is sensitive to the shape of corners on curved walls. Furthermore, the staircase approximation makes the movement of the wetting droplet on curved walls not smooth. To solve this problem, the curved boundary method may be used, but due to the interpolation or extrapolation process, most curved boundary conditions suffer from massive mass leakage when applied to the LB pseudopotential model. Through three test cases, it is found that the improved PSM scheme is mass conservative, that nearly identical static contact angles are observed on flat and curved walls under the same wetting condition, and that the movement of a wetting droplet on curved and inclined walls is smoother compared to the usual BB method. The present method is expected to be a promising tool for modeling flows in porous media and in microfluidic channels.

Mechanisms for recirculation cells in granular flows in rotating cylindrical rough tumblers.

Friction at the endwalls of partially filled horizontal rotating tumblers induces curvature and axial drift of particle trajectories in the surface flowing layer. Here we describe the results of a detailed discrete element method study of the dry granular flow of monodisperse particles in three-dimensional cylindrical tumblers with endwalls and cylindrical wall that can be either smooth or rough. Endwall roughness induces more curved particle trajectories, while a smooth cylindrical wall enhances drift near the endwall. This drift induces recirculation cells near the endwall. The use of mixed roughness (cylindrical wall and endwalls having different roughness) shows the influence of each wall on the drift and curvature of particle trajectories as well as the modification of the free surface topography. The effects act in opposite directions and have variable magnitude along the length of the tumbler such that their sum determines both direction of net drift and the recirculation cells. Near the endwalls, the dominant effect is always the endwall effect, and the axial drift for surface particles is toward the endwalls. For long enough tumblers, a counter-rotating cell occurs adjacent to each of the endwall cells having a surface drift toward the center because the cylindrical wall effect is dominant there. These cells are not dynamically coupled with the two endwall cells. The competition between the drifts induced by the endwalls and the cylindrical wall determines the width and drift amplitude for both types of cells.

Multigrid dual-time-stepping lattice Boltzmann method.

The lattice Boltzmann method often involves small numerical time steps due to the acoustic scaling (i.e., scaling between time step and grid size) inherent to the method. In this work, a second-order dual-time-stepping lattice Boltzmann method is proposed in order to avoid any time-step restriction. The implementation of the dual time stepping is based on an external source in the lattice Boltzmann equation, related to the time derivatives of the macroscopic flow quantities. Each time step is treated as a pseudosteady problem. The convergence rate of the steady lattice Boltzmann solver is improved by implementing a multigrid method. The developed solver is based on a two-relaxation time model coupled to an immersed-boundary method. The reliability of the method is demonstrated for steady and unsteady laminar flows past a circular cylinder, either fixed or towed in the computational domain. In the steady-flow case, the multigrid method drastically increases the convergence rate of the lattice Boltzmann method. The dual-time-stepping method is able to accurately reproduce the unsteady flows. The physical time step can be freely adjusted; its effect on the simulation cost is linear, while its impact on the accuracy follows a second-order trend. Two major advantages arise from this feature. (i) Simulation speed-up can be achieved by increasing the time step while conserving a reasonable accuracy. A speed-up of 4 is achieved for the unsteady flow past a fixed cylinder, and higher speed-ups are expected for configurations involving slower flow variations. Significant additional speed-up can also be achieved by accelerating transients. (ii) The choice of the time step allows us to alter the range of simulated timescales. In particular, increasing the time step results in the filtering of undesired pressure waves induced by sharp geometries or rapid temporal variations, without altering the main flow dynamics. These features may be critical to improve the efficiency and range of applicability of the lattice Boltzmann method.

Hydrodynamic model of directional ciliary-beat organization in human airways.

In the lung, the airway surface is protected by mucus, whose transport and evacuation is ensured through active ciliary beating. The mechanisms governing the long-range directional organization of ciliary beats, required for effective mucus transport, are much debated. Here, we experimentally show on human bronchial epithelium reconstituted in-vitro that the dynamics of ciliary-beat orientation is closely connected to hydrodynamic effects. To examine the fundamental mechanisms of this self-organization process, we build a two-dimensional model in which the hydrodynamic coupling between cilia is provided by a streamwise-alignment rule governing the local orientation of the ciliary forcing. The model reproduces the emergence of the mucus swirls observed in the experiments. The predicted swirl sizes, which scale with the ciliary density and mucus viscosity, are in agreement with in-vitro measurements. A transition from the swirly regime to a long-range unidirectional mucus flow allowing effective clearance occurs at high ciliary density and high mucus viscosity. In the latter case, the mucus flow tends to spontaneously align with the bronchus axis due to hydrodynamic effects.

Explicit and viscosity-independent immersed-boundary scheme for the lattice Boltzmann method.

Viscosity independence of lattice-Boltzmann methods is a crucial issue to ensure the physical relevancy of the predicted macroscopic flows over large ranges of physical parameters. The immersed-boundary (IB) method, a powerful tool that allows one to immerse arbitrary-shaped, moving, and deformable bodies in the flow, suffers from a boundary-slip error that increases as a function of the fluid viscosity, substantially limiting its range of application. In addition, low fluid viscosities may result in spurious oscillations of the macroscopic quantities in the vicinity of the immersed boundary. In this work, it is shown mathematically that the standard IB method is indeed not able to reproduce the scaling properties of the macroscopic solution, leading to a viscosity-related error on the computed IB force. The analysis allows us to propose a simple correction of the IB scheme that is local, straightforward and does not involve additional computational time. The derived method is implemented in a two-relaxation-time D2Q9 lattice-Boltzmann solver, applied to several physical configurations, namely, the Poiseuille flow, the flow around a cylinder towed in still fluid, and the flow around a cylinder oscillating in still fluid, and compared to a noncorrected immersed-boundary method. The proposed correction leads to a major improvement of the viscosity independence of the solver over a wide range of relaxation times (from 0.5001 to 50), including the correction of the boundary-slip error and the suppression of the spurious oscillations. This improvement may considerably extend the range of application of the IB lattice-Boltzmann method, in particular providing a robust tool for the numerical analysis of physical problems involving fluids of varying viscosity interacting with solid geometries.

Evidence of reverse and intermediate size segregation in dry granular flows down a rough incline.

In a dry granular flow, size segregation had been shown to behave differently for a mixture containing a few large particles with a size ratio above 5 [N. Thomas, Phys. Rev. E 62, 961 (2000)1063-651X10.1103/PhysRevE.62.961]. For moderately large size ratios, large particles migrate to an intermediate depth in the bed: this is called "intermediate segregation." For the largest size ratios, large particles migrate down to the bottom of the flow: this is called "reverse segregation," in contrast with surface segregation. As the reversal and intermediate depth values depend on the fraction of particles, this numerical study mainly uses one single large tracer. Small fractions of large beads are also computed showing the link between single tracer behavior and collective segregation process. For each device (half-filled rotating tumbler and rough plane), two (2D) and three (3D) dimensional cases are distinguished. In the tumbler, the trajectories of a large tracer show that it reaches a constant depth during the flowing phase. For large size ratios, this depth is intermediate. A progressive sinking of the depth is obtained when the size ratio is increased. The largest size ratios correspond to tracers being at the bottom of the flowing layer. All 3D simulation results are in quantitative agreement with the experimental surface, intermediate, and reverse-segregation results. In the flow down a rough incline, a large tracer reaches an equilibrium depth during flow. For large size ratios, the depth is inside the bed, at an intermediate position, and for the largest size ratios, this depth is reverse, located near the bottom. Results are slightly different for a thin or a thick flow. For 3D thick flows, the reversal between surface and bottom positions occurs within a short range of size ratios: no tracer stabilizes near half-height and two reachable intermediate depth layers exist, below the surface and above the bottom reverse layer. For 3D thin flows, all intermediate depths are reachable by a tracer, depending on the size ratio. The numerical study of larger fractions of tracers (5% or 10%) shows the three segregation patterns (surface, intermediate, reverse) corresponding to the three types of equilibrium depth. The reversal is smoother than for a single tracer, and happens around the size ratio 4.5, in good agreement with experiments.

Transport and Mixing Induced by Beating Cilia in Human Airways.

The fluid transport and mixing induced by beating cilia, present in the bronchial airways, are studied using a coupled lattice Boltzmann-Immersed Boundary solver. This solver allows the simulation of both single and multi-component fluid flows around moving solid boundaries. The cilia are modeled by a set of Lagrangian points, and Immersed Boundary forces are computed onto these points in order to ensure the no-slip velocity conditions between the cilia and the fluids. The cilia are immersed in a two-layer environment: the periciliary layer (PCL) and the mucus above it. The motion of the cilia is prescribed, as well as the phase lag between two cilia in order to obtain a typical collective motion of cilia, known as metachronal waves. The results obtained from a parametric study show that antiplectic metachronal waves are the most efficient regarding the fluid transport. A specific value of phase lag, which generates the larger mucus transport, is identified. The mixing is studied using several populations of tracers initially seeded into the pericilary liquid, in the mucus just above the PCL-mucus interface, and in the mucus far away from the interface. We observe that each zone exhibits different chaotic mixing properties. The larger mixing is obtained in the PCL layer where only a few beating cycles of the cilia are required to obtain a full mixing, while above the interface, the mixing is weaker and takes more time. Almost no mixing is observed within the mucus, and almost all the tracers do not penetrate the PCL layer. Lyapunov exponents are also computed for specific locations to assess how the mixing is performed locally. Two time scales are introduced to allow a comparison between mixing induced by fluid advection and by molecular diffusion. These results are relevant in the context of respiratory flows to investigate the transport of drugs for patients suffering from chronic respiratory diseases.

Recirculation cells for granular flow in cylindrical rotating tumblers.

To better understand the velocity field and flowing layer structure, we have performed a detailed discrete element method study of the flow of monodisperse particles in a partially filled three-dimensional cylindrical rotating tumblers. Similar to what occurs near the poles in spherical and conical tumblers, recirculation cells (secondary flows) develop near the flat endwalls of a cylindrical tumbler in which particles near the surface drift axially toward the endwall, while particles deeper in the flowing layer drift axially toward the midlength of the tumbler. Another recirculation cell with the opposite sense develops next to each endwall recirculation cell, extending to the midlength of the tumbler. For a long enough tumbler, each endwall cell is about one quarter of the tumbler diameter in length. Endwall cells are insensitive to tumbler length and relatively insensitive to rotation speed (so long as the flowing layer remains flat and continuously flowing) or fill level (from 25% to 50% full). However, for shorter tumblers (0.5 to 1.0 length/diameter aspect ratio) the endwall cell size does not change much, while center cells reduce their size and eventually disappear for the shortest tumblers. For longer tumblers (length/diameter aspect ratio larger than 2), a stagnation zone appears in between the central cells. These results provide insight into the mixing of monodisperse particles in rotating cylindrical tumblers as well as the frictional effects of the tumbler endwalls.

Influence of rough and smooth walls on macroscale granular segregation patterns.

Size bidisperse granular materials in a spherical tumbler segregate into two different patterns of three bands with either small particles at the equator and large particles at the poles or vice versa, depending upon the fill level in the tumbler. Here we use discrete element method simulations with supporting qualitative experiments to explore the effect of the tumbler wall roughness on the segregation pattern, modeling the tumbler walls as either a closely packed monolayer of fixed particles resulting in a rough wall or a frictional geometrically smooth wall. Even though the tumbler wall is in contact with the flowing layer only at its periphery, the impact of wall roughness is profound. Smooth walls tend toward a small-large-small (SLS) band pattern at the pole-equator-pole at all but the highest fill fractions; rough walls tend toward a large-small-large (LSL) band pattern at all but the lowest fill fractions. This comes about because smooth walls induce poleward axial drift of small particles and an equator-directed drift for large particles, resulting in an SLS band pattern. On the other hand, rough walls result in both sizes of particles moving poleward at the surface of the flow. Due to radial segregation, small particles percolate lower in the flowing layer and when arriving near the pole are caught in the return current drift that carries them back toward the equator incrementally with each passage through the flowing layer, while large particles remain at the surface near the pole, resulting in an LSL band pattern. The tendency toward either of the two segregation patterns depends on the fill level in the tumbler and the roughness of the tumbler's bounding wall.

Influence of rough and smooth walls on macroscale flows in tumblers.

Walls in discrete element method simulations of granular flows are sometimes modeled as a closely packed monolayer of fixed particles, resulting in a rough wall rather than a geometrically smooth wall. An implicit assumption is that the resulting rough wall differs from a smooth wall only locally at the particle scale. Here we test this assumption by considering the impact of the wall roughness at the periphery of the flowing layer on the flow of monodisperse particles in a rotating spherical tumbler. We find that varying the wall roughness significantly alters average particle trajectories even far from the wall. Rough walls induce greater poleward axial drift of particles near the flowing layer surface but decrease the curvature of the trajectories. Increasing the volume fill level in the tumbler has little effect on the axial drift for rough walls but increases the drift while reducing curvature of the particle trajectories for smooth walls. The mechanism for these effects is related to the degree of local slip at the bounding wall, which alters the flowing layer thickness near the walls, affecting the particle trajectories even far from the walls near the equator of the tumbler. Thus, the proper choice of wall conditions is important in the accurate simulation of granular flows, even far from the bounding wall.

Slow axial drift in three-dimensional granular tumbler flow.

Models of monodisperse particle flow in partially filled three-dimensional tumblers often assume that flow along the axis of rotation is negligible. We test this assumption, for spherical and double cone tumblers, using experiments and discrete element method simulations. Cross sections through the particle bed of a spherical tumbler show that, after a few rotations, a colored band of particles initially perpendicular to the axis of rotation deforms: particles near the surface drift toward the pole, while particles deeper in the flowing layer drift toward the equator. Tracking of mm-sized surface particles in tumblers with diameters of 8-14 cm shows particle axial displacements of one to two particle diameters, corresponding to axial drift that is 1-3% of the tumbler diameter, per pass through the flowing layer. The surface axial drift in both double cone and spherical tumblers is zero at the equator, increases moving away from the equator, and then decreases near the poles. Comparing results for the two tumbler geometries shows that wall slope causes axial drift, while drift speed increases with equatorial diameter. The dependence of axial drift on axial position for each tumbler geometry is similar when both are normalized by their respective maximum values.