Abstract
Efficient mixing is a persistent bottleneck in agricultural and agrochemical processing, where rapid and uniform mixing must be achieved under laminar flow with low energy input and gentle shear. Inspired by peristaltic transport in biological systems, this study investigates a bio-inspired flexible-wall squeezing mixer and establishes a two-dimensional computational framework to quantify how periodic wall deformation governs scalar homogenization in a flexible conduit. An Arbitrary Lagrangian-Eulerian dynamic mesh approach is implemented to resolve moving boundaries and to prescribe actuation, enabling the systematic evaluation of the separate and coupled effects of peak wall-normal velocity amplitude A and actuation frequency f on mixing performance. Mixing effectiveness is quantified using a variance-based mixing index MI and a sustained-threshold mixing time ts, and response surface methodology is employed to map the A-f design space and interpret the roles of time-dependent shear, interfacial stretching and folding, and vortex intensification. Relative to a non-actuated baseline, a peak wall-normal velocity amplitude of 3 × 10(-3) m s(-1) at 2 Hz reduces t(s) by 21.3%. At fixed f = 3 Hz, increasing A from 1 × 10(-3) to 4 × 10(-3) m s(-1) shortens t(s) by 10.2%, while at fixed A = 3 × 10(-3) m s(-1), raising f from 1 to 5 Hz further decreases t(s) by 6.6% with diminishing gains at the lowest frequencies. The response surface identifies an operating optimum at A = 4 × 10(-3) m s(-1) and f = 5 Hz, achieving a peak MI of 0.9557 and a minimum t(s) of 7.81 s. A periodically squeezed physical mixing loop was further examined using fluorescence imaging to assess outlet homogeneity trends. The stabilized outlet coefficient of variation (CV) decreased from about 0.65 without squeezing to 0.60 at 1 Hz and 10 mm s(-1), 0.58 at 2 Hz and 10 mm s(-1), and 0.54 at 2 Hz and 30 mm s(-1), indicating that stronger and faster actuation improves outlet uniformity. The numerical and experimental results are therefore interpreted jointly as mechanistic and trend-level evidence, while a rigorous quantitative prediction for the cylindrical compliant device will require future three-dimensional, compliance-resolved simulations and broader experimental benchmarking.