Simulating acoustically-actuated flows in complex microchannels using the volume penalization technique

TL;DR

This paper introduces a volume penalization method for simulating acoustically-driven flows in complex microchannels, effectively handling nonlinear interactions and complex geometries with scalable solvers and novel force evaluation techniques.

Contribution

The paper presents a novel volume penalization approach combined with a perturbation method and efficient solvers for simulating acoustically-actuated microfluidic flows in complex geometries, validated against body-fitted methods.

Findings

High accuracy in complex geometries demonstrated

Efficient scalable solvers for first- and second-order problems

Novel contour integration for acoustic radiation force evaluation

Abstract

We present a volume penalization technique for simulating acoustically-actuated flows in geometrically complex microchannels. Using a perturbation approach, the nonlinear response of an acoustically-actuated compressible Newtonian fluid moving over obstacles or flowing in a geometrically complex domain is segregated into two sub-problems: a harmonic first-order problem and a time-averaged second-order problem, where the latter utilizes forcing terms and boundary conditions arising from the first-order solution. This segregation results in two distinct volume penalized systems of equations. The no-slip boundary condition at the fluid-solid interface is enforced by prescribing a zero structure velocity for the first-order problem, while spatially varying Stokes drift -- which depends on the gradient of the first-order solution -- is prescribed as the structure velocity for the second-order problem. The harmonic first-order system is solved via MUMPS direct solver, whereas the steady state second-order system is solved iteratively using a novel projection method-based preconditioner. The preconditioned iterative solver for the second-order system is demonstrated to be highly effective and scalable with respect to increasing penalty force and grid resolution, respectively. A novel contour integration technique to evaluate the acoustic radiation force on an immersed object is also proposed. Through test cases featuring representative microfluidic geometries, we demonstrate excellent agreement between the volume penalized and body-fitted grid. We also identify suitable penalty factors and interfacial smearing widths to accurately capture the first- and second-order solutions. These results provide first-of-its-kind empirical evidence of the efficacy of the volume penalization method for simulating acoustic streaming problems that have so far been analyzed using body-fitted methods in literature.