The dynamics of evolving fluid films in the viscous Stokes limit is relevant to various applications, such as the modelling of lipid bilayers in cells. While the governing equations were formulated by Scriven (1960), solving for the flow of a deformable viscous surface with arbitrary shape and topology has remained a challenge. In this study, we present a straightforward discrete model based on variational principles to address this long-standing problem. We replace the classical equations, which are expressed with tensor calculus in local coordinates, with a simple coordinate-free, differential-geometric formulation. The formulation provides a fundamental understanding of the underlying mechanics and translates directly to discretization. We construct a discrete analogue of the system using Onsager's variational principle, which, in a smooth context, governs the flow of a viscous medium. In the discrete setting, instead of term-wise discretizing the coordinate-based Stokes equations, we construct a discrete Rayleighian for the system and derive the discrete Stokes equations via the variational principle. This approach results in a stable, structure-preserving variational integrator that solves the system on general manifolds.

In this contribution, we develop a versatile formalism to derive unified two-phase models describing both the separated and disperse regimes as introduced by Loison et al. (Intl J. Multiphase Flow, vol. 177, 2024, 104857). It relies on the stationary action principle and interface geometric variables. This contribution provides a novel method to derive small-scale models for the dynamics of the interface geometry. They are introduced here on a simplified case where all the scales and phases have the same velocity and that does not take into account large-scale capillary forces. The derivation tools yield a proper mathematical framework through hyperbolicity and signed entropy evolution. The formalism encompasses a hierarchy of small-scale reduced-order models based on a statistical description at a mesoscopic kinetic level and is naturally able to include the description of a disperse phase with polydispersity in size. This hierarchy includes both a cloud of spherical droplets and non-spherical droplets experiencing a dynamical behaviour through incompressible oscillations. The associated small-scale variables are moments of a number density function resulting from the geometric method of moments (GeoMOM). This method selects moments as small-scale geometric variables compatible with the structure and dynamics of the interface; they are defined independently of the flow topology and, therefore, this model allows the coupling of the two-scale flow with an inter-scale transfer. It is shown, in particular, that the resulting dynamics provides partial closures for the interface area density equation obtained from the averaging approach.

Ciliated microorganisms near the base of the aquatic food chain either swim to encounter prey or attach at a substrate and generate feeding currents to capture passing particles. Here, we represent attached and swimming ciliates using a popular spherical model in viscous fluid with slip surface velocity that affords analytical expressions of ciliary flows. We solve an advection–diffusion equation for the concentration of dissolved nutrients, where the Péclet number ($Pe$) reflects the ratio of diffusive to advective time scales. For a fixed hydrodynamic power expenditure, we ask what ciliary surface velocities maximize nutrient flux at the microorganism's surface. We find that surface motions that optimize feeding depend on $Pe$. For freely swimming microorganisms at finite $Pe$, it is optimal to swim by employing a ‘treadmill’ surface motion, but in the limit of large $Pe$, there is no difference between this treadmill solution and a symmetric dipolar surface velocity that keeps the organism stationary. For attached microorganisms, the treadmill solution is optimal for feeding at $Pe$ below a critical value, but at larger $Pe$ values, the dipolar surface motion is optimal. We verified these results in open-loop numerical simulations and asymptotic analysis, and using an adjoint-based optimization method. Our findings challenge existing claims that optimal feeding is optimal swimming across all Péclet numbers, and provide new insights into the prevalence of both attached and swimming solutions in oceanic microorganisms.

Poiseuille flow is a fundamental flow in fluid mechanics and is driven by a pressure gradient in a channel. Although the rheology of active particle suspensions has been investigated extensively, knowledge of the Poiseuille flow of such suspensions is lacking. In this study, dynamic simulations of a suspension of active particles in Poiseuille flow, situated between two parallel walls, were conducted by Stokesian dynamics assuming negligible inertia. Active particles were modelled as spherical squirmers. In the case of inert spheres in Poiseuille flow, the distribution of spheres between the walls was layered. In the case of non-bottom-heavy squirmers, on the other hand, the layers collapsed and the distribution became more uniform. This led to a much larger pressure drop for the squirmers than for the inert spheres. The effects of volume fraction, swimming mode, swimming speed and the wall separation on the pressure drop were investigated. When the squirmers were bottom heavy, they accumulated at the channel centre in downflow, whereas they accumulated near the walls in upflow, as observed in former experiments. The difference in squirmer configuration alters the hydrodynamic force on the wall and hence the pressure drop and effective viscosity. In upflow, pusher squirmers induced a considerably larger pressure drop, while neutral and puller squirmers could even generate negative pressure drops, i.e. spontaneous flow could occur. While previous studies have reported negative viscosity of pusher suspensions, this study shows that the effective viscosity of bottom-heavy puller suspensions can be negative for Poiseuille upflow, which is a new finding. The knowledge obtained is important for understanding channel flow of active suspensions.