Acoustic Streaming in A Microchannel Cross Section
Acoustic Streaming in A Microchannel Cross Section
Acoustic Streaming in A Microchannel Cross Section
Recent advances in the fabrication of microfluidic systems require handling of live cells and other
micro particles as well as mixing. All this can, for example, be achieved using acoustic radiation
forces and the viscous drag from the streaming flow.
Streaming: Due to the nonlinear terms in the Navier-Stokes equations, harmonic perturbation of
the flow will lead to a net time-averaged flow called acoustic streaming. Acoustic streaming is a
second order (nonlinear) acoustic effect. The effect can be simulated in two ways: either by direct
simulation solving the nonlinear Navier-Stokes equations, or as shown here by separation of time
scales (2nd order perturbation).
Radiation Force: Due to nonlinear terms in the governing equations momentum can be transferred
from an acoustic field to particles. This results in a net force acting on the particles – the acoustic
radiation force.
The trajectory of particles in devices will be governed by the balance between the viscous drag
force (from the streaming flow) and the acoustic radiation force. This model shows how to include
both using COMSOL Multiphysics.
Background and Motivation
This model first solves the 1st order acoustic field using the Thermoviscous Acoustics, Frequency
Domain interface. This solves the full thermoviscous acoustic problem and properly resolves the
acoustic boundary layer. The streaming flow is driven in this layer.
Secondly, the 2nd order time-averaged net flow is solved. This is achieved by adding the necessary
source terms (they stem from the first order fields) to the Laminar Flow interface.
Finally, the path of polystyrene micro particles in the microchannel cross section is modeled using
the Particle Tracing for Fluid Flow interface. The forces on the particles are added using the built-in
Acoustophoretic Radiation Force and the Drag Force features.
The model presented here is based on:
• P. B. Muller, R. Barnkob, M. J. Herring Jensen, and H. Bruus, “A numerical study of microparticle acoustophoresis driven by acoustic
radiation forces and streaming-induced drag forces,” Lab. Chip., 12, 4617–4627 (2012).
The Acoustophoretic Radiation Force feature including both thermal and viscous effects as presented
in this paper:
• J. T. Karlsen and H. Bruus, “Forces acting on a small particle in an acoustical field in a thermoviscous fluid,” Phys. Rev. E 92, 043010
(2015).
Theory
The 1st order acoustic fields (p1,u1,T1) are solved using the Thermoviscous Acoustics, Frequency
Domain interface of the Acoustics Module.
The streaming flow (time averaged) (<p2>,<u2>) is solved using the Laminar Flow physics interface of
the CFD Module by adding the appropriate time-averaged 1st order sources: a mass source and a
volume force.
The sources are not actual physical sources but part of solving the full non-linear equations. They
appear due to the solution technique.
System
H = 100 µm
W = 150 µm
Acoustic Fields
Pressure: p1
f0 = 5 MHz
Density: 1
Acoustic Fields
Temperature: T1
The streaming flow (time averaged) (<p2>,<u2>) in a microchannel cross section is solved using the
Thermoviscous Acoustics and the Laminar Flow physics interface (of the Acoustics and the CFD
modules).
Note that the bulk viscosity term has been disregarded in the streaming flow calculation. The effect
of this term is very small in this application.
The solution is achieved using second order perturbation theory and time-averaging the second
order equations manually.
The movement of micro particles in the channel due to the streaming flow (drag force) and
radiation forces is modeled using the Particle Tracing for Fluid Flow interface.
The nature and the particle movement is governed by the balance between drag and radiation
forces. This depends on the particle size a and the material parameters of the fluid and particles.