Nutrient uptake by a squirmer: The critical role of closed streamlines
Phys. Rev. Fluids 10, 074502 – Published 14 July, 2025
DOI: https://doi.org/10.1103/j4cb-s83l
Abstract
We consider the problem of determining the rate of nutrient uptake by a self-propelled spherical squirmer, a model for a microorganism that generates motion along its surface. The tangential velocity induced by the squirmer varies along the surface in such a manner that a region of fluid recirculation is formed on the rear part of the squirmer. This problem was formulated by Magar et al. [Q. J. Mech. Appl. Math. 56, 65 (2003)] who investigated it in detail numerically covering a large range of Peclet numbers, being the radius of the squirmer, is its swimming speed, and is the diffusivity of the nutrient in the surrounding liquid. Their numerical results for large suggested that Sherwood number , the nondimensional mass transfer rate, increases in proportion to but they were unable to theoretically determine the constant of proportionality. The present study reexamines this problem in the limit of large and small Reynolds numbers. The concentration in the bulk of the recirculating region approaches a constant whose value is determined by requiring that the mass transfer of nutrient across the streamline dividing the open and closed streamline regions must equal the mass transfer to the squirmer from the closed streamlines. The analysis of the nutrient concentration in the vicinity of the dividing streamline is complicated by the fact that this streamline originates from the stagnation point at the squirmer surface. The boundary layers formed at the nutrient surface on the either side of this stagnation point continue along the dividing streamline resulting in a discontinuity in the nutrient flux. As a consequence, new boundary layers form within these two larger boundary layers. When the squirmer motion induces the recirculation region ahead of itself, the distribution near the dividing streamline is simpler to determine but that near the squirmer surface involves formation of new boundary layers inside the boundary layers returning from the dividing streamline. Analytical expressions for large are derived in both cases and shown to be in agreement with the numerical results of Magar et al.