Pulsatile poromechanics in layered soft media controls fluid flow and solute transport: from fundamentals to brain clearance

Technical Summary: Pulsatile Poromechanics in Layered Soft Media

Problem Statement
Soft porous media, particularly biological tissues like cartilage and the brain, frequently exhibit heterogeneous, layered structures with varying mechanical and fluid-flow properties. While the response of homogeneous media to static and periodic loading is well-characterized, the physical implications of macroscopic layering on nonlinear poromechanics and solute transport remain poorly understood. This gap is critical for understanding brain metabolic clearance, where pathological alterations (e.g., amyloid-β accumulation) may disrupt the layered architecture of grey and white matter. The authors hypothesize that the spatial arrangement of layers, rather than just their average properties, fundamentally reshapes stress, strain, fluid flow, and solute transport under large-deformation oscillations.

Methodology
The study employs a one-dimensional, nonlinear poromechanical model in Lagrangian coordinates to simulate a bilayer soft porous medium subjected to periodic, displacement-driven loading.

• Model Formulation: The model couples mass conservation for the fluid phase, Darcy's law with a deformation-dependent Kozeny-Carman permeability relation, and mechanical equilibrium using Hencky elasticity. Solute transport is modeled via conservation of mass, incorporating molecular diffusion and hydrodynamic dispersion.
• Experimental Design: To isolate the specific role of layering, the authors compare a homogeneous reference case against four distinct bilayer configurations. Crucially, they maintain a constant poroelastic timescale ($T_{PE}$) across all cases by adjusting material properties (porosity $\Phi_0$, permeability $k_0$, and p-wave modulus $M$) such that the product $k_0M$ remains constant within layers and the effective macroscopic $T_{PE}$ is preserved. This ensures that observed differences arise from spatial heterogeneity rather than changes in the global mechanical regime.
• Numerical Implementation: The system is solved numerically in MATLAB using compact finite differences in space and implicit Runge-Kutta integration in time.
• Physiological Application: The framework is applied to a realistic cerebral grey matter (GM)–white matter (WM) bilayer system. Parameters are derived from literature for healthy tissue, while pathological scenarios (P1 and P2) simulate amyloid-β induced reductions in stiffness, permeability, and porosity within the GM layer.

Key Contributions and Results
The study demonstrates that the spatial sequencing of layers relative to the loading boundary is a critical determinant of system response, independent of the poroelastic timescale.

• Poromechanical Localization: In bilayer media, the same $T_{PE}$ that governs homogeneous response leads to non-trivial localization or propagation patterns of strain and fluid flow.
  • Case (i.a): A proximal layer with high impedance (high $M$, low $k_0$) resists local deformation, effectively transmitting strain to the distal layer. This amplifies fluid and solute fluxes deep into the medium.
  • Case (i.b): A proximal layer with high deformability (low $M$, high $k_0$) localizes strains and fluid flow near the piston, significantly attenuating distal solute transport.
  • Porosity Variations (Series ii): Lower proximal porosity enhances localization and hinders transport, whereas higher proximal porosity mitigates localization, enhancing overall clearance.
• Solute Transport Dynamics: The spatial distribution of fluid volume changes dictates transport physics more than proximal deformability alone. Configurations that minimize localization and shift deformation toward the distal end (e.g., Case i.a and ii.b) maximize solute clearance. Conversely, configurations promoting localization (e.g., Case i.b and ii.a) reduce clearance efficiency.
• Pathological Implications: In the GM-WM application, progressive pathological stiffening and clogging of the GM layer (simulating amyloid-β effects) lead to a strong localization of fluid flux at the surface. This suppresses fluid flow into the white matter and significantly impairs solute transport, suggesting a mechanical mechanism for reduced waste clearance in disease.

Significance and Claims
The authors claim that their findings reveal a fundamental dependence of poromechanical response on the spatial ordering of tissue heterogeneities. They argue that layered architectures may provide functional benefits for cellular homeostasis by optimizing fluid flow and waste transport under oscillating boundary conditions (such as blood vessel pulsations).

• Biophysical Insight: The study suggests that the brain's layered structure is inherently sensitive to pathological shifts in material properties. A disruption in the mechanical balance between layers (e.g., GM stiffening) can alter the poromechanical equilibrium, potentially creating a "vicious circle" that further hinders waste accumulation.
• Engineering Application: The results offer a foundation for designing "smart" tissue engineering scaffolds with graded or layered architectures to control local stress, flow, and nutrient/waste transport.
• Limitations and Future Work: The authors note that the pathological parameter changes used in the GM-WM case are speculative and require experimental verification. They emphasize the need for precise in vivo measurements of permeability, stiffness, and porosity across disease stages. Future work aims to explore the upscaling of solute transport to derive effective transport estimations and to refine the cerebral model using dual-porosity poroelasticity.