Source author record

Henrik Bruus

Henrik Bruus appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2020arXiv

Particle-size-dependent acoustophoretic motion and depletion of micro- and nanoparticles at long time scales

We present three-dimensional measurements of size-dependent acoustophoretic motion of microparticles with diameters from 4.8 um down to 0.5 um suspended in either homogeneous or inhomogeneous fluids inside a glass-silicon microchannel and exposed to a standing ultrasound wave. To study the cross-over from radiation force dominated to streaming dominated motion as the particle size is decreased, we extend previous studies to long time scales, where the particles smaller than the cross-over size move over distances comparable to the channel width. We observe a particle-size-dependent particle depletion at late times for the particles smaller than the cross-over size. The mechanisms behind this depletion in homogeneous fluids are rationalized by numerical simulations which take the Brownian motion into account. Experimentally, the particle trajectories in inhomogeneous fluids show focusing in the bulk of the microchannel at early times, even for the particles below the critical size, which clearly demonstrates the potential to manipulate submicrometer particles.

preprint2020arXiv

Suppression of acoustic streaming in shape-optimized channels

Acoustic streaming is an ubiquitous phenomenon resulting from time-averaged nonlinear dynamics in oscillating fluids. In this theoretical study, we show that acoustic streaming can be suppressed by two orders of magnitude in major regions of a fluid by optimizing the shape of its confining walls. Remarkably, the acoustic pressure is not suppressed in this shape-optimized cavity, and neither is the acoustic radiation force on suspended particles. This basic insight may lead to applications, such as acoustophoretic handling of nm-sized particles, which is otherwise impaired by acoustic~streaming.

preprint2020arXiv

Theory of acoustic trapping of microparticles in capillary tubes

We present a semi-analytical theory for the acoustic fields and particle-trapping forces in a viscous fluid inside a capillary tube with arbitrary cross section and ultrasound actuation at the walls. We find that the acoustic fields vary axially on a length scale proportional to the square root of the quality factor of the two-dimensional (2D) cross-section resonance mode. This axial variation is determined analytically based on the numerical solution to the eigenvalue problem in the 2D cross section. The analysis is developed in two steps: First, we generalize a recently published expression for the 2D standing-wave resonance modes in a rectangular cross section to arbitrary shapes, including the viscous boundary layer. Second, based on these 2D modes, we derive analytical expressions in three dimensions for the acoustic pressure, the acoustic radiation and trapping force, as well as the acoustic energy flux density. We validate the theory by comparison to three-dimensional numerical simulations.

preprint2019arXiv

Bulk-driven acoustic streaming at resonance in closed microcavities

Bulk-driven acoustic (Eckart) streaming is the steady flow resulting from the time-averaged acoustic energy flux density in the bulk of a viscous fluid. In simple cases, like the one-dimensional single standing-wave resonance, this energy flux is negligible, and therefore the bulk-driven streaming is often ignored relative to the boundary-driven (Rayleigh) streaming in the analysis of resonating acoustofluidic devices with length scales comparable to the acoustic wavelength. However, in closed acoustic microcavities with viscous dissipation, two overlapping resonances may be excited at the same frequency as a double mode. In contrast to single modes, the double modes can support a steady rotating acoustic energy flux density and thus a corresponding rotating bulk-driven acoustic streaming. We derive analytical solutions for the double modes in a rectangular box-shaped cavity including the viscous boundary layers, and use them to map out possible rotating patterns of bulk-driven acoustic streaming. Remarkably, the rotating bulk-driven streaming may be excited by a non-rotating actuation, and we determine the optimal geometry that maximizes this excitation. In the optimal geometry, we finally simulate a horizontal 2-by-2, 4-by-4, and 6-by-6 streaming-roll pattern in a shallow square cavity. We find that the high-frequency 6-by-6 streaming-roll pattern is dominated by the bulk-driven streaming as opposed to the low-frequency 2-by-2 streaming pattern, which is dominated by the boundary-driven streaming.

preprint2019arXiv

Three-dimensional numerical modeling of surface acoustic wave devices: Acoustophoresis of micro- and nanoparticles including streaming

Surface acoustic wave (SAW) devices form an important class of acoustofluidic devices, in which the acoustic waves are generated and propagate along the surface of a piezoelectric substrate. Despite their wide-spread use, only a few fully three-dimensional (3D) numerical simulations have been presented in the literature. In this paper, we present a 3D numerical simulation taking into account the electromechanical fields of the piezoelectric SAW device, the acoustic displacement field in the attached elastic material, in which the liquid-filled microchannel is embedded, the acoustic fields inside the microchannel, as well as the resulting acoustic radiation force and streaming-induced drag force acting on micro- and nanoparticles suspended in the microchannel. A specific device design is presented, for which the numerical predictions of the acoustic resonances and the acoustophoretic repsonse of suspended microparticles in 3D are successfully compared with experimental observations. The simulation provides a physical explanation of the the observed qualitative difference between devices with an acoustically soft and hard lid in terms of traveling and standing waves, respectively. The simulations also correctly predict the existence and position of the observed in-plane streaming flow rolls. The presented simulation model may be useful in the development of SAW devices optimized for various acoustofluidic tasks.

preprint2016arXiv

Modeling of microdevices for SAW-based acoustophoresis --- a study of boundary conditions

We present a finite-element method modeling of acoustophoretic devices consisting of a single, long, straight, water-filled microchannel surrounded by an elastic wall of either borosilicate glass (pyrex) or the elastomer polydimethylsiloxane (PDMS) and placed on top of a piezoelectric transducer that actuates the device by surface acoustic waves (SAW). We compare the resulting acoustic fields in these full solid-fluid models with those obtained in reduced fluid models comprising of only a water domain with simplified, approximate boundary conditions representing the surrounding solids. The reduced models are found to only approximate the acoustically hard pyrex systems to a limited degree for large wall thicknesses and not at all for the acoustically soft PDMS systems.

preprint2016arXiv

The acoustic force density acting on inhomogeneous fluids in acoustic fields

We present a theory for the acoustic force density acting on inhomogeneous fluids in acoustic fields on time scales that are slow compared to the acoustic oscillation period. The acoustic force density depends on gradients in the density and compressibility of the fluid. For microfluidic systems, the theory predicts a relocation of the inhomogeneities into stable field-dependent configurations, which are qualitatively different from the horizontally layered configurations due to gravity. Experimental validation is obtained by confocal imaging of aqueous solutions in a glass-silicon microchip.

preprint2015arXiv

A sharp-interface model of electrodeposition and ramified growth

We present a sharp-interface model of two-dimensional ramified growth during quasi-steady electrodeposition. Our model differs from previous modeling methods in that it includes the important effects of extended space-charge regions and nonlinear electrode reactions. The model is validated by comparing its behavior in the initial stage with the predictions of a linear stability analysis.

preprint2015arXiv

A theoretical study of time-dependent, ultrasound-induced acoustic streaming in microchannels

Based on first- and second-order perturbation theory, we present a numerical study of the temporal build-up and decay of unsteady acoustic fields and acoustic streaming flows actuated by vibrating walls in the transverse cross-sectional plane of a long straight microchannel under adiabatic conditions and assuming temperature-independent material parameters. The unsteady streaming flow is obtained by averaging the time-dependent velocity field over one oscillation period, and as time increases, it is shown to converge towards the well-known steady time-averaged solution calculated in the frequency domain. Scaling analysis reveals that the acoustic resonance builds up much faster than the acoustic streaming, implying that the radiation force may dominate over the drag force from streaming even for small particles. However, our numerical time-dependent analysis indicates that pulsed actuation does not reduce streaming significantly due to its slow decay. Our analysis also shows that for an acoustic resonance with a quality factor Q, the amplitude of the oscillating second-order velocity component is Q times larger than the usual second-order steady time-averaged velocity component. Consequently, the well-known criterion v << c for the validity of the perturbation expansion is replaced by the more restrictive criterion v << c/Q. Our numerical model is available in the supplemental material in the form of Comsol model files and Matlab scripts.

preprint2015arXiv

Continuum Nanofluidics

This paper introduces the fundamental continuum theory governing momentum transport in isotropic nanofluidic flows. The theory is an extension to the classical Navier-Stokes equation, which includes coupling between translational and rotational degrees of freedom, as well as non-local response functions that incorporates spatial correlations. The continuum theory is compared with molecular dynamics simulation data for both relaxation processes and fluid flows showing excellent agreement on the nanometer length scale. We also present practical tools to estimate when the extended theory should be used.

preprint2015arXiv

Forces acting on a small particle in an acoustical field in a thermoviscous fluid

We present a theoretical analysis of the acoustic radiation force on a single small particle, either a thermoviscous fluid droplet or a thermoelastic solid particle, suspended in a viscous and heat-conducting fluid medium. Our analysis places no restrictions on the length scales of the viscous and thermal boundary layer thicknesses $δ_\mathrm{s}$ and $δ_\mathrm{t}$ relative to the particle radius $a$, but it assumes the particle to be small in comparison to the acoustic wavelength $λ$. This is the limit relevant to scattering of sound and ultrasound waves from micrometer-sized particles. For particles of size comparable to or smaller than the boundary layers, the thermoviscous theory leads to profound consequences for the acoustic radiation force. Not only do we predict forces orders of magnitude larger than expected from ideal-fluid theory, but for certain relevant choices of materials, we also find a sign change in the acoustic radiation force on different-sized but otherwise identical particles. This phenomenon may possibly be exploited in handling of submicrometer-sized particles such as bacteria and vira in lab-on-a-chip systems.

preprint2015arXiv

Morphological instability during steady electrodeposition at overlimiting currents

We present a linear stability analysis of a planar metal electrode during steady electrodeposition. We extend the previous work of Sundstrom and Bark by accounting for the extended space-charge density, which develops at the cathode once the applied voltage exceeds a few thermal voltages. In accordance with Chazalviel's conjecture, the extended space-charge region is found to greatly affect the morphological stability of the electrode. To supplement the numerical solution of the stability problem, we have derived analytical expressions valid in limit of low and high voltage, respectively.

preprint2014arXiv

A numerical study of thermoviscous effects in ultrasound-induced acoustic streaming in microchannels

We present a numerical study of thermoviscous effects on the acoustic streaming flow generated by an ultrasound standing-wave resonance in a long straight microfluidic channel containing a Newtonian fluid. These effects enter primarily through the temperature and density dependence of the fluid viscosity. The resulting magnitude of the streaming flow is calculated and characterized numerically, and we find that even for thin acoustic boundary layers, the channel height affects the magnitude of the streaming flow. For the special case of a sufficiently large channel height we have successfully validated our numerics with analytical results from 2011 by Rednikov and Sadhal for a single planar wall. We analyze the time-averaged energy transport in the system and the time-averaged second-order temperature perturbation of the fluid. Finally, we have made three main changes in our previously published numerical scheme to improve the numerical performance: (i) The time-averaged products of first-order variables in the time-averaged second-order equations have been recast as flux densities instead of as body forces. (ii) The order of the finite element basis functions has been increased in an optimal manner. (iii) Based on the International Association for the Properties of Water and Steam (IAPWS 1995, 2008, and 2011), we provide accurate polynomial fits in temperature for all relevant thermodynamic and transport parameters of water in the temperature range from 10 C to 50 C.

preprint2014arXiv

Acoustic interaction forces between small particles in an ideal fluid

We present a theoretical expression for the acoustic interaction force between small spherical particles suspended in an ideal fluid exposed to an external acoustic wave. The acoustic interaction force is the part of the acoustic radiation force on one given particle involving the scattered waves from the other particles. The particles, either compressible liquid droplets or elastic microspheres, are considered to be much smaller than the acoustic wavelength. In this so-called Rayleigh limit, the acoustic interaction forces between the particles are well approximated by gradients of pair-interaction potentials with no restriction on the inter-particle distance. The theory is applied to studies of the acoustic interaction force on a particle suspension in either standing or traveling plane waves. The results show aggregation regions along the wave propagation direction, while particles may attract or repel each other in the transverse direction. In addition, a mean-field approximation is developed to describe the acoustic interaction force in an emulsion of oil droplets in water.

preprint2014arXiv

Concentration polarization, surface currents, and bulk advection in a microchannel

We present a comprehensive analysis of salt transport and overlimiting currents in a microchannel during concentration polarization. We have carried out full numerical simulations of the coupled Poisson-Nernst-Planck-Stokes problem governing the transport and rationalized the behaviour of the system. A remarkable outcome of the investigations is the discovery of strong couplings between bulk advection and the surface current; without a surface current, bulk advection is strongly suppressed. The numerical simulations are supplemented by analytical models valid in the long channel limit as well as in the limit of negligible surface charge. By including the effects of diffusion and advection in the diffuse part of the electric double layers, we extend a recently published analytical model of overlimiting current due to surface conduction.

preprint2013arXiv

Transport-limited water splitting at ion-selective interfaces during concentration polarization

We present an analytical model of salt- and water-ion transport across an ion-selective interface based on an assumption of local equilibrium of the water-dissociation reaction. The model yields current-voltage characteristics and curves of water-ion current versus salt-ion current, which are in qualitative agreement with experimental results published in the literature. The analytical results are furthermore in agreement with direct numerical simulations. As part of the analysis, we find approximate solutions to the classical problem of pure salt transport across an ion-selective interface. These solutions provide closed-form expressions for the current-voltage characteristics, which include the overlimiting current due to the development of an extended space charge region. Finally, we discuss how the addition of an acid or a base affects the transport properties of the system and thus provide predictions accessible to further experimental tests of the model.

preprint2013arXiv

Ultrasound-induced acoustophoretic motion of microparticles in three dimensions

We derive analytical expressions for the three-dimensional (3D) acoustophoretic motion of spherical microparticles in rectangular microchannels. The motion is generated by the acoustic radiation force and the acoustic streaming-induced drag force. In contrast to the classical theory of Rayleigh streaming in shallow, infinite, parallel-plate channels, our theory does include the effect of the microchannel side walls. The resulting predictions agree well with numerics and experimental measurements of the acoustophoretic motion of polystyrene spheres with nominal diameters of 0.537 um and 5.33 um. The 3D particle motion was recorded using astigmatism particle tracking velocimetry under controlled thermal and acoustic conditions in a long, straight, rectangular microchannel actuated in one of its transverse standing ultrasound-wave resonance modes with one or two half-wavelengths. The acoustic energy density is calibrated in situ based on measurements of the radiation dominated motion of large 5-um-diam particles, allowing for quantitative comparison between theoretical predictions and measurements of the streaming induced motion of small 0.5-um-diam particles.

preprint2012arXiv

Acoustic radiation- and streaming-induced microparticle velocities determined by micro-PIV in an ultrasound symmetry plane

We present micro-PIV measurements of suspended microparticles of diameters from 0.6 um to 10 um undergoing acoustophoresis in an ultrasound symmetry plane in a microchannel. The motion of the smallest particles are dominated by the Stokes drag from the induced acoustic streaming flow, while the motion of the largest particles are dominated by the acoustic radiation force. For all particle sizes we predict theoretically how much of the particle velocity is due to radiation and streaming, respectively. These predictions include corrections for particle-wall interactions and ultrasonic thermoviscous effects, and they match our measurements within the experimental uncertainty. Finally, we predict theoretically and confirm experimentally that the ratio between the acoustic radiation- and streaming-induced particle velocities is proportional to the square of the particle size, the actuation frequency and the acoustic contrast factor, while it is inversely proportional to the kinematic viscosity.

preprint2012arXiv

Diving with microparticles in acoustic fields

Sound can move particles. A good example of this phenomenon is the Chladni plate, in which an acoustic wave is induced in a metallic plate and particles migrate to the nodes of the acoustic wave. For several years, acoustophoresis has been used to manipulate microparticles in microscopic scales. In this fluid dynamics video, submitted to the 30th Annual Gallery of Fluid Motion, we show the basic mechanism of the technique and a simple way of visualize it. Since acoustophoretic phenomena is essentially a three-dimensional effect, we employ a simple technique to visualize the particles in 3D. The technique is called Astigmatism Particle Tracking Velocimetry and it consists in the use of cylindrical lenses to induce a deformation in the particle shape, which will be then correlated with its distance from the observer. With this method we are able to dive with the particles and observe in detail particle motion that would otherwise be missed. The technique not only permits visualization but also precise quantitative measurements that can be compared with theory and simulations.

preprint2011arXiv

On the forces acting on a small particle in an acoustical field in a viscous fluid

We calculate the acoustic radiation force from an ultrasound wave on a compressible, spherical particle suspended in a viscous fluid. Using Prandtl--Schlichting boundary-layer theory, we include the kinematic viscosity of the solvent and derive an analytical expression for the resulting radiation force, which is valid for any particle radius and boundary-layer thickness provided that both of these length scales are much smaller than the wavelength of the ultrasound wave (mm in water at MHz frequencies). The acoustophoretic response of suspended microparticles is predicted and analyzed using parameter values typically employed in microchannel acoustophoresis.

preprint2009arXiv

A numerical analysis of finite Debye-length effects in induced-charge electro-osmosis

For a microchamber filled with a binary electrolyte and containing a flat un-biased center electrode at one wall, we employ three numerical models to study the strength of the resulting induced-charge electro-osmotic (ICEO) flow rolls: (i) a full nonlinear continuum model resolving the double layer, (ii) a linear slip-velocity model not resolving the double layer and without tangential charge transport inside this layer, and (iii) a nonlinear slip-velocity model extending the linear model by including the tangential charge transport inside the double layer. We show that compared to the full model, the slip-velocity models significantly overestimate the ICEO flow. This provides a partial explanation of the quantitative discrepancy between observed and calculated ICEO velocities reported in the literature. The discrepancy increases significantly for increasing Debye length relative to the electrode size, i.e. for nanofluidic systems. However, even for electrode dimensions in the micrometer range, the discrepancies in velocity due to the finite Debye length can be more than 10% for an electrode of zero height and more than 100% for electrode heights comparable to the Debye length.

preprint2009arXiv

Kilohertz microfluidics as an analytical tool for determining dynamic characteristics of microfluidic systems

The advent in recent years of highly parallelized microfluidic chemical reaction systems necessitates an understanding of all fluid dynamic time scales including the often neglected millisecond time scale of the inertia of the liquid. We propose the use of harmonically oscillating microfluidics in the low kilohertz range as an analytical tool for the deduction of these time scales. Furthermore, we suggest the use of systems-level equivalent circuit theory as an adequate theory of the behavior of the system. A novel pressure source capable of operation in the desired frequency range is presented for this generic analysis. As a proof of concept, we study the fairly complex system of water-filled interconnected elastic microfluidic tubes containing a large, trapped air-bubble and driven by a pulsatile pressure difference. We demonstrate good agreement between the systems-level model and the experimental results, allowing us to determine the dynamic time scales of the system.

preprint2009arXiv

Strongly nonlinear dynamics of electrolytes in large ac voltages

We study the response of a model micro-electrochemical cell to a large ac voltage of frequency comparable to the inverse cell relaxation time. To bring out the basic physics, we consider the simplest possible model of a symmetric binary electrolyte confined between parallel-plate blocking electrodes, ignoring any transverse instability or fluid flow. We analyze the resulting one-dimensional problem by matched asymptotic expansions in the limit of thin double layers and extend previous work into the strongly nonlinear regime, which is characterized by two novel features - significant salt depletion in the electrolyte near the electrodes and, at very large voltage, the breakdown of the quasi-equilibrium structure of the double layers. The former leads to the prediction of "ac capacitive desalination", since there is a time-averaged transfer of salt from the bulk to the double layers, via oscillating diffusion layers. The latter is associated with transient diffusion limitation, which drives the formation and collapse of space-charge layers, even in the absence of any net Faradaic current through the cell. We also predict that steric effects of finite ion sizes (going beyond dilute solution theory) act to suppress the strongly nonlinear regime in the limit of concentrated electrolytes, ionic liquids and molten salts. Beyond the model problem, our reduced equations for thin double layers, based on uniformly valid matched asymptotic expansions, provide a useful mathematical framework to describe additional nonlinear responses to large ac voltages, such as Faradaic reactions, electro-osmotic instabilities, and induced-charge electrokinetic phenomena.

preprint2009arXiv

Topology and shape optimization of induced-charge electro-osmotic micropumps

For a dielectric solid surrounded by an electrolyte and positioned inside an externally biased parallel-plate capacitor, we study numerically how the resulting induced-charge electro-osmotic (ICEO) flow depends on the topology and shape of the dielectric solid. In particular, we extend existing conventional electrokinetic models with an artificial design field to describe the transition from the liquid electrolyte to the solid dielectric. Using this design field, we have succeeded in applying the method of topology optimization to find system geometries with non-trivial topologies that maximize the net induced electro-osmotic flow rate through the electrolytic capacitor in the direction parallel to the capacitor plates. Once found, the performance of the topology optimized geometries has been validated by transferring them to conventional electrokinetic models not relying on the artificial design field. Our results show the importance of the topology and shape of the dielectric solid in ICEO systems and point to new designs of ICEO micropumps with significantly improved performance.

preprint2008arXiv

Multi-directional sorting modes in deterministic lateral displacement devices

Deterministic lateral displacement (DLD) devices separate micrometer-scale particles in solution based on their size using a laminar microfluidic flow in an array of obstacles. We investigate array geometries with rational row-shift fractions in DLD devices by use of a simple model including both advection and diffusion. Our model predicts novel multi-directional sorting modes that could be experimentally tested in high-throughput DLD devices containing obstacles that are much smaller than the separation between obstacles.

preprint2007arXiv

Acoustic resonances in microfluidic chips: full-image micro-PIV experiments and numerical simulations

We show that full-image micro-PIV analysis in combination with images of transient particle motion is a powerful tool for experimental studies of acoustic radiation forces and acoustic streaming in microfluidic chambers under piezo-actuation in the MHz range. The measured steady-state motion of both large 5 um and small 1 um particles can be understood in terms of the acoustic eigenmodes or standing ultra-sound waves in the given experimental microsystems. This interpretation is supported by numerical solutions of the corresponding acoustic wave equation.

preprint1996arXiv

Energy level statistics of the two-dimensional Hubbard model at low filling

The energy level statistics of the Hubbard model for $L \times L$ square lattices (L=3,4,5,6) at low filling (four electrons) is studied numerically for a wide range of the coupling strength. All known symmetries of the model (space, spin and pseudospin symmetry) have been taken into account explicitly from the beginning of the calculation by projecting into symmetry invariant subspaces. The details of this group theoretical treatment are presented with special attention to the nongeneric case of L=4, where a particular complicated space group appears. For all the lattices studied, a significant amount of levels within each symmetry invariant subspaces remains degenerated, but except for L=4 the ground state is nondegenerate. We explain the remaining degeneracies, which occur only for very specific interaction independent states, and we disregard these states in the statistical spectral analysis. The intricate structure of the Hubbard spectra necessitates a careful unfolding procedure, which is thoroughly discussed. Finally, we present our results for the level spacing distribution, the number variance $Σ^2$, and the spectral rigidity $Δ_3$, which essentially all are close to the corresponding statistics for random matrices of the Gaussian ensemble independent of the lattice size and the coupling strength. Even very small coupling strengths approaching the integrable zero coupling limit lead to the Gaussian ensemble statistics stressing the nonperturbative nature of the Hubbard model.

preprint1996arXiv

Parametric correlation of Coulomb blockade conductance peaks in chaotic quantum dots

We investigate the autocorrelator of conductance peak heights for quantum dots in the Coulomb blockade regime. Analytical and numerical results based on Random Matrix Theory are presented and compared to exact numerical calculations based on a simple dynamical model. We consider the case of preserved time-reversal symmetry, which is realized experimentally by varying the shape of the quantum dot in the absence of magnetic fields. Upon a proper rescaling, the correlator becomes independent of the details of the system and its form is solely determined by symmetry properties and the number of channels in the leads. The magnitude of the scaling parameter is estimated by a semiclassical approach.