On the rectification of oscillatory flows by flexible leaflets in a confined geometry

Problem Statement
Fluid-structure interactions (FSI) in compliant channels and confined geometries are fundamental to biological transport and microfluidic design. While local nonlinearities induced by single flexible structures are well-characterized, predicting the collective behavior of multiple interacting structures under time-varying flows remains a significant challenge. Existing literature predominantly focuses on steady-state regimes or single structures, often neglecting the rich temporal dynamics that emerge from the interaction of multiple compliant elements. Specifically, the mechanisms governing how arrays of flexible leaflets rectify oscillatory flows in low-Reynolds-number environments—breaking the inherent time-reversal symmetry of the Stokes regime to generate net unidirectional transport—are not fully understood.

Methodology
The authors investigate this problem using a minimal two-dimensional model of a channel containing a bed of asymmetric leaflets anchored to the bottom wall. The system is driven by the vertical oscillation of the top wall, simulating a squeeze flow. The study employs a fully coupled numerical approach combining:

1. Lattice Boltzmann Method (LBM): Used to solve the incompressible Navier-Stokes equations in the low-Mach-number limit. To handle the time-dependent nature of the problem without prohibitive time-step restrictions, a dual-time-stepping (DTS) scheme accelerated by a multigrid approach is utilized.
2. Immersed Boundary Method (IBM): Used to integrate the flexible leaflet geometry into the fluid domain. Leaflets are modeled as rigid plates rotating around a pivot, attached to a rotational spring, responding to hydrodynamic torques.
3. Analytical Framework: A continuum-limit model is developed to isolate the influence of leaflet elasticity and density. This model utilizes local torque balance and global angular momentum conservation to predict flow rectification in the steady limit, which is subsequently extended via a quasi-steady approximation to the oscillatory regime.

Key Contributions and Results
The study characterizes the generation of net fluid transport through the interplay of leaflet elasticity, density, and channel geometry. Key findings include:

• Role of the Elasto-Viscous Number ($\eta$): The system's behavior is governed by the dimensionless elasto-viscous number, defined as the ratio of viscous hydrodynamic forces to the restorative elastic forces of the leaflets. The net flow exhibits a non-monotonic dependence on $\eta$:
  • Stiff leaflets ($\eta \ll 1$): Deformations are negligible, resulting in minimal flow rectification.
  • Optimal flexibility: An intermediate value of $\eta$ maximizes transport. Here, leaflets deform sufficiently to create geometric asymmetry between contraction and relaxation phases without losing their structural integrity.
  • Excessively soft leaflets ($\eta \gg 1$): Leaflets deform excessively (often past 90 degrees), reducing the geometric asymmetry and diminishing net transport.
• Collective Interactions and Density: High leaflet densities maximize collective interactions and net transport. As density increases, the discrete leaflet array approaches a continuum limit where the effective boundary profile becomes smooth. The net flow rate increases with density until saturating, a behavior captured by a phenomenological scaling function $Q(\eta, \phi) \sim \phi/(\phi + K(\eta))$.
• Steady vs. Dynamic Regimes: In the steady limit, the net flow is independent of the Reynolds number and oscillation frequency, depending solely on $\eta$ and density. However, in the fully oscillatory regime, the system exhibits a phase lag between the wall motion and the leaflet response. This dynamic delay is controlled by a dimensionless timescale ratio $T_r$, representing the ratio of the leaflet's viscous relaxation time to the wall oscillation period. When $T_r$ is significant, transient configurations arise that diminish the net rectified transport, a phenomenon not captured by purely quasi-static models.
• Geometric Dependence: While the quasi-steady model accurately predicts flow in the high-density limit, it fails to capture the dependence of net flow on channel length ($L_x$) observed in simulations, highlighting the importance of dynamic effects in shorter channels.

Significance
The paper establishes a foundational framework for analyzing how collective slender structures interact dynamically within viscous environments. By developing an analytical continuum model validated against fully coupled FSI simulations, the authors provide a predictive tool for understanding flow rectification driven by flexible arrays. The work demonstrates that collective interactions can break the symmetries of the Stokes regime to generate unidirectional transport, offering insights relevant to biological fluid transport systems and the design of passive microfluidic rectifiers. The identification of an optimal elasto-viscous number and the characterization of dynamic phase lags provide specific design parameters for controlling flow in compliant geometries.