Showing 1–37 of 37 results for author: Versluis, M

Role of surfactants on droplet formation in piezoacoustic inkjet printing across microsecond-to-second timescales

Authors: Maaike Rump, Christian Diddens, Uddalok Sen, Michel Versluis, Detlef Lohse, Tim Segers

Abstract: In piezo acoustic drop-on-demand inkjet printing a single droplet is produced for each piezo driving pulse. This droplet is typically multicomponent, including surfactants to control the spreading and drying of the droplet on the substrate. However, the role of these surfactants on the droplet formation process remains rather elusive. Surfactant concentration gradients may manifest across microsec… ▽ More In piezo acoustic drop-on-demand inkjet printing a single droplet is produced for each piezo driving pulse. This droplet is typically multicomponent, including surfactants to control the spreading and drying of the droplet on the substrate. However, the role of these surfactants on the droplet formation process remains rather elusive. Surfactant concentration gradients may manifest across microsecond-to-second timescales, spanning both the rapid ejection of ink from the nozzle exit and the comparatively slower idling timescale governing the firing of successive droplets. In the present work, we study the influence of surfactants on droplet formation across 6 orders of magnitude in time. To this end, we visualize the microsecond droplet formation process using stroboscopic laser-induced fluorescence microscopy while we vary the nozzle idle time. Our results show that increasing the idle time up to O(1) s affects only the break-up dynamics of the inkjet but not its velocity. By contrast, for idle times $>$ O(1) s, both the break-up dynamics are altered and the velocity of the inkjet increases. We show that the increased velocity results from a decreased surface tension of the ejected droplet, which we extracted from the observed shape oscillations of the jetted droplets in flight. The measured decrease in surface tension is surprising as the $μ$s timescale of droplet formation is much faster than the typical ms-to-s timescale of surfactant adsorption. By varying the bulk surfactant concentration, we show that the fast decrease in surface tension results from a local surfactant concentration increase to more than 200 times the CMC. Our results suggest that a local high concentration of surfactant allows for surfactant adsorption to the interface of an inkjet at the $μ$s-to-ms timescale, which is much faster than the typical ms-to-s timescale of surfactant adsorption. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: 13 pages, 8 figures

Vaporization dynamics of a super-heated water-in-oil droplet: modeling and numerical solution

Authors: Muhammad Saeed Saleem, Michel Versluis, Guillaume Lajoinie

Abstract: The study of vapor bubble growth following droplet vaporization in a superheated liquid involves research areas such as hydrodynamics, heat transfer, mass transfer, and thermodynamics. The interplay between these multiscale aspects is strongly dependent on the geometry, the thermodynamic response, and the local physical properties of the system. To understand the role of each aspect of this comple… ▽ More The study of vapor bubble growth following droplet vaporization in a superheated liquid involves research areas such as hydrodynamics, heat transfer, mass transfer, and thermodynamics. The interplay between these multiscale aspects is strongly dependent on the geometry, the thermodynamic response, and the local physical properties of the system. To understand the role of each aspect of this complex mechanism we model super-heated droplet vaporization by coupling the equation of motion for bubble growth with the thermodynamics of phase change and heat transfer through the convection-diffusion equation. The semi-analytical model is validated with the analytical description for vapor bubble growth dominated either by inertia (Rayleigh) or by thermal diffusion (Plesset-Zwick), depending on droplet radius and degree of superheat. The effect of a mismatch of the thermal properties between the host liquid and the droplet is shown to be relevant only for low superheating, above which an increase in thermal diffusivity leads to a reduction in the rate of vaporization. At medium to high superheating, the droplet vaporizes completely without relying on thermal diffusion. At the point of complete vaporization, the potential energy within the system drives the bubble overshoots, which vary based on the droplet size and degree of superheat. △ Less

Submitted 26 April, 2024; originally announced April 2024.

A unifying Rayleigh-Plesset-type equation for bubbles in viscoelastic media

Authors: Alexandros T. Oratis, Kay Dijs, Guillaume Lajoinie, Michel Versluis, Jacco H. Snoeijer

Abstract: Understanding the ultrasound pressure-driven dynamics of microbubbles confined in viscoelastic materials is relevant for multiple biomedical applications, ranging from contrast-enhanced ultrasound imaging to ultrasound-assisted drug delivery. The volumetric oscillations of spherical bubbles is analyzed using the Rayleigh-Plesset equation, which describes the conservation of mass and momentum in th… ▽ More Understanding the ultrasound pressure-driven dynamics of microbubbles confined in viscoelastic materials is relevant for multiple biomedical applications, ranging from contrast-enhanced ultrasound imaging to ultrasound-assisted drug delivery. The volumetric oscillations of spherical bubbles is analyzed using the Rayleigh-Plesset equation, which describes the conservation of mass and momentum in the surrounding medium. Several studies have considered an extension of the Rayleigh-Plesset equation for bubbles embedded into viscoelastic media, but these are restricted to a particular choice of constitutive model and/or to small deformations. Here, we derive a unifying equation applicable to bubbles in viscoelastic media with arbitrary complex moduli and that can account for large bubble deformations. To derive this equation, we borrow concepts from finite-strain theory. We validate our approach by comparing the result of our model to previously published results and extend it to show how microbubbles behave in arbitrary viscoelastic materials. In particular, we use our viscoelastic Rayleigh-Plesset model to compute the bubble dynamics in benchmarked viscoelastic liquids and solids. △ Less

Submitted 1 November, 2023; originally announced November 2023.

Comments: 13 pages, 6 figures

An anthropomorphic thyroid phantom for ultrasound-guided radiofrequency ablation of nodules

Authors: Tim Boers, Wyger Brink, Leonardo Bianchi, Paola Saccomandi, Johan van Hespen, Germen Wennemars, Sicco Braak, Michel Versluis, Srirang Manohar

Abstract: Background: Needle-based procedures such as fine needle aspiration (FNA) and thermal ablation, are often applied for thyroid nodule diagnosis and therapeutic purposes, respectively. With blood vessels and nerves nearby, these procedures can pose risks in damaging surrounding critical structures. Purpose: The development and validation of innovative strategies to manage these risks require a test o… ▽ More Background: Needle-based procedures such as fine needle aspiration (FNA) and thermal ablation, are often applied for thyroid nodule diagnosis and therapeutic purposes, respectively. With blood vessels and nerves nearby, these procedures can pose risks in damaging surrounding critical structures. Purpose: The development and validation of innovative strategies to manage these risks require a test object with well-characterized physical properties. For this work, we focus on the application of ultrasound-guided thermal radio-frequency ablation (RFA). Methods: We have developed an anthropomorphic phantom mimicking the thyroid and surrounding anatomical and physiological structures that are relevant to ultrasound-guided thermal ablation. The phantom was composed of a mixture of polyacrylamide, water, and egg white extract and was cast using molds in multiple steps. The thermal, acoustical, and electrical characteristics were experimentally validated. The ablation zones were analyzed via non-destructive T2-weighted MRI scans utilizing the relaxometry changes of coagulated egg albumen, and the temperature distribution was monitored using an array of fiber Bragg sensors. Results: The physical properties of the phantom were verified both on ultrasound as well as its response to thermal ablation. The final temperature achieved (92°C), the median percentage of the nodule ablated (82.1%), the median volume ablated outside the nodule (0.8 mL), and the median number of critical structures affected (0) were quantified. Conclusion: An anthropomorphic phantom that can provide a realistic model for development and training in ultrasound-guided needle-based thermal interventions for thyroid nodules has been presented. In the future, this model can also be extended to novel needle-based diagnostic procedures. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: 19 pages, 10 figures, 3 tables

Journal ref: Med Phys. 2023

A Review of Machine Learning Applications for the Proton Magnetic Resonance Spectroscopy Workflow

Authors: Dennis M. J. van de Sande, Julian P. Merkofer, Sina Amirrajab, Mitko Veta, Ruud J. G. van Sloun, Maarten J. Versluis, Jacobus F. A. Jansen, Johan S. van den Brink, Marcel Breeuwer

Abstract: This literature review presents a comprehensive overview of machine learning (ML) applications in proton magnetic resonance spectroscopy (MRS). As the use of ML techniques in MRS continues to grow, this review aims to provide the MRS community with a structured overview of the state-of-the-art methods. Specifically, we examine and summarize studies published between 2017 and 2023 from major journa… ▽ More This literature review presents a comprehensive overview of machine learning (ML) applications in proton magnetic resonance spectroscopy (MRS). As the use of ML techniques in MRS continues to grow, this review aims to provide the MRS community with a structured overview of the state-of-the-art methods. Specifically, we examine and summarize studies published between 2017 and 2023 from major journals in the magnetic resonance field. We categorize these studies based on a typical MRS workflow, including data acquisition, processing, analysis, and artificial data generation. Our review reveals that ML in MRS is still in its early stages, with a primary focus on processing and analysis techniques, and less attention given to data acquisition. We also found that many studies use similar model architectures, with little comparison to alternative architectures. Additionally, the generation of artificial data is a crucial topic, with no consistent method for its generation. Furthermore, many studies demonstrate that artificial data suffers from generalization issues when tested on in-vivo data. We also conclude that risks related to ML models should be addressed, particularly for clinical applications. Therefore, output uncertainty measures and model biases are critical to investigate. Nonetheless, the rapid development of ML in MRS and the promising results from the reviewed studies justify further research in this field. △ Less

Submitted 16 May, 2023; originally announced May 2023.

Resonance behavior of a compliant piezo-driven inkjet channel with an entrained microbubble

Authors: Hans Reinten, Yogesh Jethani, Arjan Fraters, Roger Jeurissen, Detlef Lohse, Michel Versluis, Tim Segers

Abstract: Microbubbles entrained in a piezo-driven drop-on-demand (DOD) printhead disturb the acoustics of the microfluidic ink channel and thereby the jetting behavior. Here, the resonance behavior of an ink channel as a function of the microbubble size and the number of bubbles is studied through theoretical modeling and experiments. The system is modeled as a set of two coupled harmonic oscillators: one… ▽ More Microbubbles entrained in a piezo-driven drop-on-demand (DOD) printhead disturb the acoustics of the microfluidic ink channel and thereby the jetting behavior. Here, the resonance behavior of an ink channel as a function of the microbubble size and the number of bubbles is studied through theoretical modeling and experiments. The system is modeled as a set of two coupled harmonic oscillators: one corresponding to the compliant ink channel and one to the microbubble. The predicted and measured eigenfrequencies are in excellent agreement. It was found that the resonance frequency is independent of the bubble size as long as the compliance of the bubble dominates over that of the piezo actuator. An accurate description of the eigenfrequency of the coupled system requires the inclusion of the increased inertance of the entrained microbubble due to confinement. We show that the inertance of a confined bubble can be accurately obtained by using a simple potential flow approach. The model is further validated by the excellent agreement between the modeled and measured microbubble resonance curves. The present work therefore provides physical insight in the coupled dynamics of a compliant ink channel with an entrained microbubble. △ Less

Submitted 24 February, 2023; originally announced February 2023.

Selective evaporation at the nozzle exit in piezoacoustic inkjet printing

Authors: Maaike Rump, Uddalok Sen, Roger Jeurissen, Hans Reinten, Michel Versluis, Detlef Lohse, Christian Diddens, Tim Segers

Abstract: In practical applications of inkjet printing the nozzles in a printhead have intermittent idle periods, during which ink can evaporate from the nozzle exit. Inks are usually multicomponent where each component has its own characteristic evaporation rate resulting in concentration gradients within the ink. These gradients may directly and indirectly (via Marangoni flows) alter the jetting process a… ▽ More In practical applications of inkjet printing the nozzles in a printhead have intermittent idle periods, during which ink can evaporate from the nozzle exit. Inks are usually multicomponent where each component has its own characteristic evaporation rate resulting in concentration gradients within the ink. These gradients may directly and indirectly (via Marangoni flows) alter the jetting process and thereby its reproducibility and the resulting print quality. In the present work, we study selective evaporation from an inkjet nozzle for water-glycerol mixtures. Through experiments, analytical modeling, and numerical simulations, we investigate changes in mixture composition with drying time. By monitoring the acoustics within the printhead, and subsequently modeling the system as a mass-spring-damper system, the composition of the mixture can be obtained as a function of drying time. The results from the analytical model are validated using numerical simulations of the full fluid mechanical equations governing the printhead flows and pressure fields. Furthermore, the numerical simulations reveal that the time independent concentration gradient we observe in the experiments is due to the steady state of water flux through the printhead. Finally, we measure the number of drop formation events required in this system before the mixture concentration within the nozzle attains the initial (pre-drying) value, and find a stronger than exponential trend in the number of drop formations required. These results shed light on the complex physiochemical hydrodynamics associated with the drying of ink at a printhead nozzle, and help in increasing the stability and reproducibility of inkjet printing. △ Less

Submitted 27 February, 2023; v1 submitted 13 May, 2022; originally announced May 2022.

Comments: For supplementary movie, see https://www.youtube.com/watch?v=a9PF8gOvABw

Vorticity-induced flow-focusing leads to bubble entrainment in an inkjet printhead: synchrotron X-ray and volume-of-fluid visualizations

Authors: Maaike Rump, Youssef Saade, Uddalok Sen, Kamel Fezzaa, Michel Versluis, Detlef Lohse, Tim Segers

Abstract: The oscillatory flows present in an inkjet printhead can lead to strong deformations of the air-liquid interface at the nozzle exit. Such deformations may lead to an inward directed air jet with bubble pinch-off and the subsequent entrainment of an air bubble, which is highly detrimental to the stability of inkjet printing. Understanding the mechanisms of bubble entrainment is therefore crucial in… ▽ More The oscillatory flows present in an inkjet printhead can lead to strong deformations of the air-liquid interface at the nozzle exit. Such deformations may lead to an inward directed air jet with bubble pinch-off and the subsequent entrainment of an air bubble, which is highly detrimental to the stability of inkjet printing. Understanding the mechanisms of bubble entrainment is therefore crucial in improving print stability. In the present work, we use ultrafast X-ray phase-contrast imaging and direct numerical simulations based on the Volume-of-Fluid method to study the mechanisms underlying the bubble entrainment in a piezo-acoustic printhead. We first demonstrate good agreement between experiments and numerics. We then show the different classes of bubble pinch-off obtained in experiments, and that those were also captured numerically. The numerical results are then used to show that the baroclinic torque, which is generated at the gas-liquid interface due to the misalignment of density and pressure gradients, results in a flow-focusing effect that drives the formation of the air jet from which a bubble can pinch-off. △ Less

Submitted 8 August, 2022; v1 submitted 13 May, 2022; originally announced May 2022.

Comments: For supplementary movies, visit: https://www.youtube.com/playlist?list=PLmU3YKXYzq7oEm9qKD9bL_eWKSNP3uGtw

Super-Resolved Microbubble Localization in Single-Channel Ultrasound RF Signals Using Deep Learning

Authors: Nathan Blanken, Jelmer M. Wolterink, Hervé Delingette, Christoph Brune, Michel Versluis, Guillaume Lajoinie

Abstract: Recently, super-resolution ultrasound imaging with ultrasound localization microscopy (ULM) has received much attention. However, ULM relies on low concentrations of microbubbles in the blood vessels, ultimately resulting in long acquisition times. Here, we present an alternative super-resolution approach, based on direct deconvolution of single-channel ultrasound radio-frequency (RF) signals with… ▽ More Recently, super-resolution ultrasound imaging with ultrasound localization microscopy (ULM) has received much attention. However, ULM relies on low concentrations of microbubbles in the blood vessels, ultimately resulting in long acquisition times. Here, we present an alternative super-resolution approach, based on direct deconvolution of single-channel ultrasound radio-frequency (RF) signals with a one-dimensional dilated convolutional neural network (CNN). This work focuses on low-frequency ultrasound (1.7 MHz) for deep imaging (10 cm) of a dense cloud of monodisperse microbubbles (up to 1000 microbubbles in the measurement volume, corresponding to an average echo overlap of 94%). Data are generated with a simulator that uses a large range of acoustic pressures (5-250 kPa) and captures the full, nonlinear response of resonant, lipid-coated microbubbles. The network is trained with a novel dual-loss function, which features elements of both a classification loss and a regression loss and improves the detection-localization characteristics of the output. Whereas imposing a localization tolerance of 0 yields poor detection metrics, imposing a localization tolerance corresponding to 4% of the wavelength yields a precision and recall of both 0.90. Furthermore, the detection improves with increasing acoustic pressure and deteriorates with increasing microbubble density. The potential of the presented approach to super-resolution ultrasound imaging is demonstrated with a delay-and-sum reconstruction with deconvolved element data. The resulting image shows an order-of-magnitude gain in axial resolution compared to a delay-and-sum reconstruction with unprocessed element data. △ Less

Submitted 9 April, 2022; originally announced April 2022.

Comments: Accepted for publication in IEEE Transactions on Medical Imaging

US Velocimetry in Participants with Aortoiliac Occlusive Disease

Authors: Stefan Engelhard, Majorie van Helvert, Jason Voorneveld, Johan G Bosch, Guillaume PR Lajoinie, Michel Versluis, Erik Groot Jebbink, Michel MPJ Reijnen

Abstract: The accurate quantification of blood flow in aortoiliac arteries is challenging but clinically relevant because local flow patterns can influence atherosclerotic disease. To investigate the feasibility and clinical application of two-dimensional blood flow quantification using high-frame-rate contrast-enhanced US (HFR-CEUS) and particle image velocimetry (PIV), or US velocimetry, in participants w… ▽ More The accurate quantification of blood flow in aortoiliac arteries is challenging but clinically relevant because local flow patterns can influence atherosclerotic disease. To investigate the feasibility and clinical application of two-dimensional blood flow quantification using high-frame-rate contrast-enhanced US (HFR-CEUS) and particle image velocimetry (PIV), or US velocimetry, in participants with aortoiliac stenosis. In this prospective study, participants with a recently diagnosed aortoiliac stenosis underwent HFR-CEUS measurements of the pre- and poststenotic vessel segments. Two-dimensional quantification of blood flow was achieved by performing PIV analysis, which was based on pairwise cross-correlation of the HFR-CEUS images. Visual inspection of the entire data set was performed by five observers to evaluate the ability of the technique to enable adequate visualization of blood flow. The contrast-to-background ratio and average vector correlation were calculated. In two participants who showed flow disturbances, the flow complexity and vorticity were calculated. Results: 35 participants were included. Visual scoring showed that flow quantification was achieved in 41 of 42 locations. In 25 locations, one or multiple issues occurred that limited optimal flow quantification, including loss of correlation during systole, shadow regions, a short vessel segment in the image plane, and loss of contrast during diastole. In the remaining 16 locations, optimal quantification was achieved. The contrast-to-background ratio was higher during systole than during diastole, whereas the vector correlation was lower. Flow complexity and vorticity were high in regions with disturbed flow. Blood flow quantification with US velocimetry is feasible in patients with an aortoiliac stenosis, but several challenges must be overcome before implementation into clinical practice. △ Less

Submitted 7 February, 2022; originally announced February 2022.

Journal ref: Radiology 301(2), 332-338 (2021)

Time-resolved velocity and pressure field quantification in a flow-focusing device for ultrafast microbubble production

Authors: Sarah Cleve, Christian Diddens, Tim Segers, Guillaume Lajoinie, Michel Versluis

Abstract: Flow-focusing devices have gained great interest in the past decade, due to their capability to produce monodisperse microbubbles for diagnostic and therapeutic medical ultrasound applications. However, up-scaling production to industrial scale requires a paradigm shift from single chip operation to highly parallelized systems. Parallelization gives rise to fluidic interactions between nozzles tha… ▽ More Flow-focusing devices have gained great interest in the past decade, due to their capability to produce monodisperse microbubbles for diagnostic and therapeutic medical ultrasound applications. However, up-scaling production to industrial scale requires a paradigm shift from single chip operation to highly parallelized systems. Parallelization gives rise to fluidic interactions between nozzles that, in turn, may lead to a decreased monodispersity. Here, we study the velocity and pressure field fluctuations in a single flow-focusing nozzle during bubble production. We experimentally quantify the velocity field inside the nozzle at 100 ns time resolution, and a numerical model provides insight into both the oscillatory velocity and pressure fields. Our results demonstrate that, at the length scale of the flow focusing channel, the velocity oscillations propagate at fluid dynamical time scale (order of microseconds) whereas the dominant pressure oscillations are linked to the bubble pinch-off and propagate at a much faster time scale (order of nanoseconds). △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 30 pages, 7 figures

Higher-order meniscus oscillations driven by flow-focusing leading to bubble pinch off and entrainment in a piezo acoustic inkjet nozzle

Authors: Arjan Frates, Maaike Rump, Roger Jeurissen, Marc van den Berg, Youri de Loore, Hans Reinten, Herman Wijshoff, Devaraj van der Meer, Detlef Lohse, Michel Versluis, Tim Segers

Abstract: The stability of high-end piezo-acoustic drop-on-demand (DOD) inkjet printing is sometimes compromised by the entrainment of an air bubble inside the ink channel. Here, bubble pinch-off from an acoustically driven meniscus is studied in an optically transparent DOD printhead as a function of the driving waveform. We show that bubble pinch-off follows from low-amplitude higher-order meniscus oscill… ▽ More The stability of high-end piezo-acoustic drop-on-demand (DOD) inkjet printing is sometimes compromised by the entrainment of an air bubble inside the ink channel. Here, bubble pinch-off from an acoustically driven meniscus is studied in an optically transparent DOD printhead as a function of the driving waveform. We show that bubble pinch-off follows from low-amplitude higher-order meniscus oscillations on top of the global high-amplitude meniscus motion that drives droplet formation. In a certain window of control parameters, phase inversion between the low and high frequency components leads to the enclosure of an air cavity and bubble pinch-off. Although phenomenologically similar, bubble pinch-off is not a result of capillary wave interaction such as observed in drop impact on a liquid pool. Instead, we reveal geometrical flow focusing as the mechanism through which at first, an outward jet is formed on the retracted concave meniscus. When the subsequent high-frequency pressure wave hits the now toroidal-shaped meniscus, it accelerates the toroidal ring outward resulting in the formation of an air cavity that can pinch off. The critical control parameters for pinch off are the pulse timing and amplitude. To cure the bubble entrainment problem, the threshold for bubble pinch-off can be increased by suppressing the high frequency acoustic waves through appropriate waveform design. The present work therefore aids the improvement of the stability of inkjet printers through a physical understanding of meniscus instabilities. △ Less

Submitted 7 May, 2021; originally announced May 2021.

Feedback-controlled microbubble generator producing one million monodisperse bubbles per second

Authors: Benjamin van Elburg, Gonzalo Collado Lara, Gert-Wim Bruggert, Tim Segers, Michel Versluis, Guillaume Lajoinie

Abstract: Monodisperse lipid-coated microbubbles are a promising route to unlock the full potential of ultrasound contrast agents for medical diagnosis and therapy. Here, we present a stand-alone lab-on-a-chip instrument that allows microbubbles to be formed with high monodispersity at high production rates. Key to maintaining a long-term stable, controlled, and safe operation of the microfluidic device wit… ▽ More Monodisperse lipid-coated microbubbles are a promising route to unlock the full potential of ultrasound contrast agents for medical diagnosis and therapy. Here, we present a stand-alone lab-on-a-chip instrument that allows microbubbles to be formed with high monodispersity at high production rates. Key to maintaining a long-term stable, controlled, and safe operation of the microfluidic device with full control over the output size distribution is an optical transmission based measurement technique that provides real-time information on production rate and bubble size. We feed the data into a feedback loop and demonstrate that this system can control the on-chip bubble radius (2.5 to 20 $μ$m) and the production rate up to $10^6$ bubbles/s. The freshly formed phospholipid-coated bubbles stabilize after their formation to a size approximately two times smaller than their initial on-chip bubble size without loss of monodispersity. The feedback control technique allows for full control over the size distribution of the agent and can aid the development of microfluidic platforms operated by non-specialist end users. △ Less

Submitted 23 February, 2021; originally announced April 2021.

Comments: We've have added the supplementary information with the submission in the same format

Journal ref: Review of Scientific Instruments 92, 035110 (2021)

High-frequency acoustic droplet vaporization is initiated by resonance

Authors: Guilllaume Lajoinie, Tim Segers, Michel Versluis

Abstract: Vaporization of low-boiling point droplets has numerous applications in combustion, process engineering and in recent years, in clinical medicine. However, the physical mechanisms governing the phase conversion are only partly explained. Here, we show that an acoustic resonance can arise from the large speed of sound mismatch between a perfluorocarbon microdroplet and its surroundings. The fundame… ▽ More Vaporization of low-boiling point droplets has numerous applications in combustion, process engineering and in recent years, in clinical medicine. However, the physical mechanisms governing the phase conversion are only partly explained. Here, we show that an acoustic resonance can arise from the large speed of sound mismatch between a perfluorocarbon microdroplet and its surroundings. The fundamental resonance mode obeys a unique relationship $kR$~$\sim$ 0.65 between droplet size and driving frequency that leads to a 3-fold pressure amplification inside the droplet. Classical nucleation theory shows that this pressure amplification increases the nucleation rate by several orders of magnitude. These findings are confirmed by high-speed imaging performed at a timescale of ten nanoseconds. The optical recordings demonstrate that droplets exposed to intense acoustic waves generated by inter-digital-transducers nucleate only if they match the theoretical resonance size. △ Less

Submitted 23 February, 2021; originally announced February 2021.

Comments: the package also contains the supplementary meterial, in the same format as the paper

Journal ref: Phys. Rev. Lett. 126, 034501, Published 19 January 2021

The retraction of jetted slender viscoelastic liquid filaments

Authors: Uddalok Sen, Charu Datt, Tim Segers, Herman Wijshoff, Jacco H. Snoeijer, Michel Versluis, Detlef Lohse

Abstract: Long and slender liquid filaments are produced during inkjet printing, which can subsequently either retract to form a single droplet, or break up to form a primary droplet and one or more satellite droplets. These satellite droplets are undesirable since they degrade the quality and reproducibility of the print. Existing strategies for the suppression of satellite droplet formation include, among… ▽ More Long and slender liquid filaments are produced during inkjet printing, which can subsequently either retract to form a single droplet, or break up to form a primary droplet and one or more satellite droplets. These satellite droplets are undesirable since they degrade the quality and reproducibility of the print. Existing strategies for the suppression of satellite droplet formation include, among others, adding viscoelasticity to the ink. In the present work, we aim to improve the understanding of the role of viscoelasticity in suppressing satellite droplets in inkjet printing. We demonstrate that very dilute viscoelastic aqueous solutions (concentrations $\sim$ 0.003\% wt. polyethylene oxide (PEO), corresponding to nozzle Deborah number $De_{n}$ $\sim$ 3) can suppress satellite droplet formation. Furthermore, we show that, for a given driving condition, upper and lower bounds of polymer concentration exist, within which satellite droplets are suppressed. Satellite droplets are formed at concentrations below the lower bound, while jetting ceases for concentrations above the upper bound (for fixed driving conditions). Moreover, we observe that, with concentrations in between the two bounds, the filaments retract at velocities larger than the corresponding Taylor-Culick velocity for the Newtonian case. We show that this enhanced retraction velocity can be attributed to the elastic tension due to polymer stretching, which builds up during the initial jetting phase. These results shed some light on the complex interplay between inertia, capillarity, and viscoelasticity for retracting liquid filaments, which is important for the stability and quality of inkjet printing of polymer solutions. △ Less

Submitted 9 October, 2021; v1 submitted 17 February, 2021; originally announced February 2021.

Comments: 22 pages, 7 figures

Journal ref: J. Fluid Mech. 929 (2021) A25

Evaporation-Induced Crystallization of Surfactants in Sessile Multicomponent Droplets

Authors: Yaxing Li, Valentin Salvator, Herman Wijshoff, Michel Versluis, Detlef Lohse

Abstract: Surfactants have been widely studied and used in controlling droplet evaporation. In this work, we observe and study the crystallization of sodium dodecyl sulfate (SDS) within an evaporating glycerol-water mixture droplet. The crystallization is induced by the preferential evaporation of water, which decreases the solubility of SDS in the mixture. As a consequence, the crystals shield the droplet… ▽ More Surfactants have been widely studied and used in controlling droplet evaporation. In this work, we observe and study the crystallization of sodium dodecyl sulfate (SDS) within an evaporating glycerol-water mixture droplet. The crystallization is induced by the preferential evaporation of water, which decreases the solubility of SDS in the mixture. As a consequence, the crystals shield the droplet surface and cease the evaporation. The universality of the evaporation characteristics for a range of droplet sizes is revealed by applying a diffusion model, extended by Raoult's law. To describe the nucleation and growth of the crystals, we employ the 2-dimensional crystallization model of Weinberg [J. Non-Cryst. Solids 1991, 134, 116]. The results of this model compare favorably to our experimental results. Our findings may inspire the community to reconsider the role of high concentration of surfactants in a multicomponent evaporation system. △ Less

Submitted 27 July, 2020; originally announced July 2020.

Rayleigh-Taylor instability by segregation in an evaporating multi-component microdroplet

Authors: Yaxing Li, Christian Diddens, Tim Segers, Herman Wijshoff, Michel Versluis, Detlef Lohse

Abstract: The evaporation of multi-component droplets is relevant to various applications but challenging to study due to the complex physicochemical dynamics. Recently, Li (2018) reported evaporation-triggered segregation in 1,2-hexanediol-water binary droplets. In this present work, we added 0.5 wt% silicone oil into the 1,2-hexanediol-water binary solution. This minute silicone oil concentration dramatic… ▽ More The evaporation of multi-component droplets is relevant to various applications but challenging to study due to the complex physicochemical dynamics. Recently, Li (2018) reported evaporation-triggered segregation in 1,2-hexanediol-water binary droplets. In this present work, we added 0.5 wt% silicone oil into the 1,2-hexanediol-water binary solution. This minute silicone oil concentration dramatically modifies the evaporation process as it triggers an early extraction of the 1,2-hexanediol from the mixture. Surprisingly, we observe that the segregation of 1,2-hexanediol forms plumes, rising up from the rim of the sessile droplet towards the apex during the droplet evaporation. By orientating the droplet upside down, i.e., by studying a pendant droplet, the absence of the plumes indicates that the flow structure is induced by buoyancy, which drives a Rayleigh-Taylor instability (i.e., driven by density differences & gravitational acceleration). From micro-PIV measurement, we further prove that the segregation of the non-volatile component (1,2-hexanediol) hinders the evaporation near the contact line, which leads to a suppression of the Marangoni flow in this region. Hence, on long time scales, gravitational effects play the dominant role in the flow structure, rather than Marangoni flows. We compare the measurement of the evaporation rate with the diffusion model of Popov (2005), coupled with Raoult's law and the activity coefficient. This comparison indeed confirms that the silicone-oil-triggered segregation of the non-volatile 1,2-hexanediol significantly delays the evaporation. With an extended diffusion model, in which the influence of the segregation has been implemented, the evaporation can be well described. △ Less

Submitted 28 May, 2020; originally announced May 2020.

Microuidics control the ballistic energy of thermocavitation liquid jets for needle-free injections

Authors: Loreto Oyarte Galvez, Arjan Fraters, Herman L. Offerhaus, Michel Versluis, Ian W. Hunter, David Fernandez Rivas

Abstract: Illuminating a water solution with a focused continuous wave laser produces a strong local heating of the liquid that leads to the nucleation of bubbles, also known as thermocavitation. During the growth of the bubble, the surrounding liquid is expelled from the constraining microfluidic channel through a nozzle, creating a jet. The characteristics of the resulting liquid jet was imaged using ultr… ▽ More Illuminating a water solution with a focused continuous wave laser produces a strong local heating of the liquid that leads to the nucleation of bubbles, also known as thermocavitation. During the growth of the bubble, the surrounding liquid is expelled from the constraining microfluidic channel through a nozzle, creating a jet. The characteristics of the resulting liquid jet was imaged using ultra-fast imaging techniques. Here, we provide a phenomenological description of the jet shapes and velocities, and compare them with a Boundary Integral numerical model. We define the parameter regime, varying jet speed, taper geometry and liquid volume, for optimal printing, injection and spray applications. These results are important for the design of energy-efficient needle-free jet injectors based on microfluidic thermocavitation. △ Less

Submitted 2 March, 2020; originally announced March 2020.

Giant and explosive plasmonic bubbles by delayed nucleation

Authors: Yuliang Wang, Mikhail E. Zaytsev, Guillaume Lajoinie, Hai Le The, Jan C. T. Eijkel, Albert van den Berg, Michel Versluis, Bert M. Weckhuysen, Xuehua Zhang, Harold J. W. Zandvliet, Detlef Lohse

Abstract: When illuminated by a laser, nano particles immersed in water can very quickly and strongly heat up, leading to the nucleation of so called vapor bubbles, which have huge application potential in e.g. solar light-harvesting, catalysis, and for medical applications.Whilst the long-time behavior of such bubbles has been well-studied, here, by employing ultra-high-speed imaging, we reveal the nucleat… ▽ More When illuminated by a laser, nano particles immersed in water can very quickly and strongly heat up, leading to the nucleation of so called vapor bubbles, which have huge application potential in e.g. solar light-harvesting, catalysis, and for medical applications.Whilst the long-time behavior of such bubbles has been well-studied, here, by employing ultra-high-speed imaging, we reveal the nucleation and early life phase of these bubbles. After some delay time after beginning of the illumination, a giant bubble explosively grows, up to a maximal radius of 80 um, and collapses again within 200 us (bubble life phase 1). The maximal bubble volume remarkably increases with decreasing laser power P. To explain this behavior, we measure the delay time from the beginning of the illumination up to nucleation,which drastically increases with decreasing laser power, leading to less total dumped energy . This dumped energy E shows a universal linear scaling relation. This finding supports that the initial giant bubble is a pure vapor bubble. In contrast, the delay time does depend on the gas concentration of the water, as gas pockets in the water facilitate an earlier vapor bubble nucleation, which leads to smaller delay times and lower bubble nucleation temperatures. After the collapse of the initial giant bubbles, first much smaller oscillating bubbles form out of the remaining gas nuclei (bubble life phase 2, up to typically 10 ms). Subsequently the known vaporization dominated growth phase takes over and the bubble stabilizes (life phase 3). In the final life phase 4 the bubble slowly grows by gas expelling due to heating of the surrounding. Our findings on the explosive growth and collapse during the early life phase of a vapor bubble have strong bearings on possible applications of such bubbles, affecting their risk assessment. △ Less

Submitted 12 March, 2019; v1 submitted 4 March, 2019; originally announced March 2019.

Gravitational effect in evaporating binary microdroplets

Authors: Yaxing Li, Christian Diddens, Pengyu Lv, Herman Wijshoff, Michel Versluis, Detlef Lohse

Abstract: The flow in an evaporating glycerol-water binary sub-millimeter droplet with Bond number Bo $\ll$ 1 is studied both experimentally and numerically. First, we measure the flow fields near the substrate by micro-PIV for both sessile and pendant droplets during the evaporation process, which surprisingly show opposite radial flow directions -- inward and outward, respectively. This observation clearl… ▽ More The flow in an evaporating glycerol-water binary sub-millimeter droplet with Bond number Bo $\ll$ 1 is studied both experimentally and numerically. First, we measure the flow fields near the substrate by micro-PIV for both sessile and pendant droplets during the evaporation process, which surprisingly show opposite radial flow directions -- inward and outward, respectively. This observation clearly reveals that in spite of the small droplet size, gravitational effects play a crucial role in controlling the flow fields in the evaporating droplets. We theoretically analyze that this gravity-driven effect is triggered by the lower volatility of glycerol which leads to a preferential evaporation of water then the local concentration difference of the two components leads to a density gradient that drives the convective flow. We show that the Archimedes number Ar is the nondimensional control parameter for the occurrence of the gravitational effects. We confirm our hypothesis by experimentally comparing two evaporating microdroplet systems, namely a glycerol-water droplet and a 1,2-propanediol-water droplet. We obtain different Ar, larger or smaller than a unit by varying a series of droplet heights, which corresponds to cases with or without gravitational effects, respectively. Finally, we simulate the process numerically, finding good agreement with the experimental results and again confirming our interpretation. △ Less

Submitted 5 February, 2019; originally announced February 2019.

Partial renal coverage in EVAR causes unfavorable renal flow patterns in an infrarenal aneurysm model

Authors: Lennart van de Velde, Esmé J. Donselaar, Erik Groot Jebbink, Jorrit T. Boersen, Guillaume P. J. Lajoinie, Jean-Paul M. de Vries, Clark J. Zeebregts, Michel Versluis, Michel M. P. J. Reijnen

Abstract: Objective: To achieve an optimal sealing zone during EVAR, the intended positioning of the proximal end of the endograft fabric should be as close as possible to the most caudal edge of the renal arteries. Some endografts exhibit a small offset between the radiopaque markers and the proximal fabric edge. Unintended partial renal artery coverage may thus occur. This study investigates the consequen… ▽ More Objective: To achieve an optimal sealing zone during EVAR, the intended positioning of the proximal end of the endograft fabric should be as close as possible to the most caudal edge of the renal arteries. Some endografts exhibit a small offset between the radiopaque markers and the proximal fabric edge. Unintended partial renal artery coverage may thus occur. This study investigates the consequences of partial coverage on renal flow patterns and wall shear stress. Methods: In-vitro models of an abdominal aortic aneurysm were used to visualize pulsatile flow using 2D particle image velocimetry under physiologic resting conditions. One model served as control and two models were stented with an Endurant endograft, one without and one with partial renal artery coverage with 1.3 mm of stent fabric extending beyond the marker (16\% area coverage). The magnitude and oscillation of wall shear stress, relative residence time and backflow in the renal artery were analyzed. Results: In both stented models, a region along the caudal renal artery wall presented with low and oscillating wall shear stress, not present in the control model. A region with very low wall shear stress (<0.1 Pa) was present in the model with partial coverage over a length of 7 mm, compared to a length of 2 mm in the model without renal coverage. Average renal backflow area percentage in the renal artery incrementally increased from control (0.9%) to the stented model without (6.4%) and with renal coverage (18.8%). Conclusion: In this flow model partial renal coverage after EVAR causes low and marked oscillations in wall shear stress, potentially promoting atherosclerosis and subsequent renal artery stenosis. Awareness of the device-dependent offset between the fabric edge and the radiopaque markers is therefore important in endovascular practice. △ Less

Submitted 3 January, 2019; originally announced January 2019.

Comments: 8 pages, 6 figures, supplementary video at publisher

Journal ref: J.Vasc.Surg 67 (2018) 1585-1594

Microdroplets nucleation by dissolution of a multicomponent drop in a host liquid

Authors: Huanshu Tan, Christian Diddens, Ali Akash Mohammed, Junyi Li, Michel Versluis, Xuehua Zhang, Detlef Lohse

Abstract: Multicomponent liquid drops in a host liquid are very relevant in various technological applications. Their dissolution or growth dynamics is complex. Differences in solubility between the drop components combined with the solutal Marangoni effect and natural convections contribute to this complexity, which can be even further increased in combination with the ouzo effect, i.e., the spontaneous nu… ▽ More Multicomponent liquid drops in a host liquid are very relevant in various technological applications. Their dissolution or growth dynamics is complex. Differences in solubility between the drop components combined with the solutal Marangoni effect and natural convections contribute to this complexity, which can be even further increased in combination with the ouzo effect, i.e., the spontaneous nucleation of microdroplets due to composition-dependent miscibilities in a ternary system. The quantitive understanding of this combined process is important for applications in industry, particularly for modern liquid-liquid microextraction processes. In this work, as a model system, we experimentally and theoretically explore water/ethanol drops dissolving in anethole oil. During the dissolution, we observed two types of microdroplets nucleation, namely water microdroplet nucleation in the surrounding oil at drop midheight and oil microdroplet nucleation in the aqueous drop, again at midheight. The nucleated oil microdroplets are driven by Marangoni flows inside the aqueous drop and evolve into microdroplets rings. A one-dimensional multiphase and multicomponent diffusion model in combination with thermodynamical equilibrium theory is proposed to predict the behavior of spontaneous emulsification, i.e. the microdroplet nucleation, that is triggered by diffusion. A scale analysis together with experimental investigations of the fluid dynamics of the system reveals that both the solutal Marangoni flow inside the drop and the buoyancy-driven flow in the host liquid influence the diffusion-triggered emulsification process. Our work provides a physical understanding of the microdroplet nucleation by dissolution of a multicomponent drop in a host liquid. △ Less

Submitted 8 November, 2018; originally announced November 2018.

Optical verification and in-vitro characterization of two commercially available acoustic bubble counters for cardiopulmonary bypass systems

Authors: Tim Segers, Marco C. Stehouwer, Filip M. J. J. de Somer, Bastian A. de Mol, Michel Versluis

Abstract: Gaseous microemboli (GME) introduced during cardiac surgery are considered as a potential source of morbidity, which has driven the development of the first bubble counters. Two new generation bubble counters claim correct sizing and counting of GME. This in-vitro study aims to validate the accuracy of two bubble counters using monodisperse bubbles in a highly controlled setting at low GME concent… ▽ More Gaseous microemboli (GME) introduced during cardiac surgery are considered as a potential source of morbidity, which has driven the development of the first bubble counters. Two new generation bubble counters claim correct sizing and counting of GME. This in-vitro study aims to validate the accuracy of two bubble counters using monodisperse bubbles in a highly controlled setting at low GME concentrations. Monodisperse GME with a radius of 43 μm were produced in a microfluidic chip. Directly after their formation, they were injected one-by-one into the BCC200 and the EDAC sensors. GME size and count were optically verified using high-speed imaging. During best-case scenarios or low GME concentrations of GME with a size of 43 μm in radius in an in-vitro setup, the BCC200 overestimates GME size by a factor of 2 to 3 while the EDAC underestimates the average GME size by at least a factor of two. The BCC200 overestimates the GME concentration by approximately 20% while the EDAC overestimates the concentration by nearly one order of magnitude. Nevertheless, the calculated total GME volume is only over-predicted by a factor 2 since the EDAC underestimates the actual GME size. For the BCC200, the total GME volume was over-predicted by 25 times due to the over-estimation of GME size. The measured errors in the absolute sizing/counting of GME do not imply that all results obtained using the bubble counters are insignificant or invalid. A relative change in bubble size or bubble concentration can accurately be measured. However, care must be taken in the interpretation of the results and their absolute values. Moreover, the devices cannot be used interchangeably when reporting GME activity. Nevertheless, both devices can be used to study the relative air removal characteristics of CPB components or for the quantitative monitoring of GME production during CPB interventions. △ Less

Submitted 24 July, 2018; originally announced July 2018.

Journal ref: Perfusion, 2018, 33(1), 16-24

Monodisperse versus polydisperse ultrasound contrast agents: nonlinear response, sensitivity, and deep tissue imaging potential

Authors: Tim Segers, Pieter Kruizinga, Maarten Kok, Guillaume Lajoinie, Nico de Jong, Michel Versluis

Abstract: Monodisperse microbubble ultrasound contrast agents have been proposed to further increase the signal-to-noise-ratio of contrast enhanced ultrasound imaging. Here, the sensitivity of a polydisperse preclinical agent was compared experimentally to that of its size- and acoustically-sorted derivatives by using narrowband pressure- and frequency-dependent scattering and attenuation measurements. The… ▽ More Monodisperse microbubble ultrasound contrast agents have been proposed to further increase the signal-to-noise-ratio of contrast enhanced ultrasound imaging. Here, the sensitivity of a polydisperse preclinical agent was compared experimentally to that of its size- and acoustically-sorted derivatives by using narrowband pressure- and frequency-dependent scattering and attenuation measurements. The sorted monodisperse agents showed up to a two orders of magnitude increase in sensitivity, i.e. in the average scattering cross-section per bubble. Moreover, we demonstrate here, for the first time, that the highly nonlinear response of acoustically sorted microbubbles can be exploited to confine scattering and attenuation to the focal region of ultrasound fields used in clinical imaging. This property is a result of minimal prefocal scattering and attenuation and can be used to minimize shadowing effects in deep tissue imaging. Moreover, it potentially allows for more localized therapy using microbubbles through the spatial control of resonant microbubble oscillations. △ Less

Submitted 12 June, 2018; originally announced June 2018.

Journal ref: Ultrasound in Medicine & Biology Volume 44, Issue 7, July 2018, Pages 1482-1492 Ultrasound in Medicine & Biology Volume 44, Issue 7, July 2018, Pages 1482-1492 Ultrasound in Medicine and Biology, 44(7), 2018, 1482-1492

Three year outcome of the Covered Endovascular Reconstruction of the Aortic Bifurcation (CERAB) technique for aortoiliac occlusive disease

Authors: Kim Taeymans, Erik Groot Jebbink, Suzanne Holewijn, Jasper Martens, Michel Versluis, Peter Goverde, Michel Reijnen

Abstract: Objective:The objective of the current study was to demonstrate the three year outcome of the CERAB technique for the treatment of extensive aortoiliac occlusive disease (AIOD). Methods:Between February 2009 and July 2016, all patients treated with the CERAB technique for AIOD were identified in the local databases of two centers and analyzed. Demographics and lesion characteristics were scored. F… ▽ More Objective:The objective of the current study was to demonstrate the three year outcome of the CERAB technique for the treatment of extensive aortoiliac occlusive disease (AIOD). Methods:Between February 2009 and July 2016, all patients treated with the CERAB technique for AIOD were identified in the local databases of two centers and analyzed. Demographics and lesion characteristics were scored. Follow-up (FU) consisted of clinical assessment, duplex ultrasound and ankle brachial indices (ABI). Patency rates and clinically driven target lesion revascularization (CD-TLR) were calculated by Kaplan-Meier analysis. Results: 130 patients (69 male and 61 female) were treated of which 68% patients were diagnosed with intermittent claudication and 32% suffered from critical limb ischemia. The vast majority (89%) were TASC-II D lesions and the remaining were TASC-II B and C (both 5%). Median follow-up was 24 months (range 0-67 months). The technical success rate was 97% and 67% of cases were performed completely percutaneously. The ABI improved significantly from 0.65 plusminus 0.22 preoperatively to 0.88 plusminus 0.15 after the procedure. The 30-day minor and major complication rate was 33% and 7%. The median hospital stay was 2 days (range 1-76 days). At 1 and 3-years FU 94% and 96% of the patients clinically improved at least 1 Rutherford category (2% and 0% unchanged, 4% and 4% worsened). Limb salvage rate at 1-year was 98% and 97% at three year follow-up. Primary, primary-assisted and secondary patency was 86%/91%/97% at 1-year, 84%/89%/97% at 2-year and 82%/87%/97% at 3-year FU. Freedom from CD-TLR was 87% at 1-year and 86% at both 2 and 3-year FU. Conclusions:The CERAB technique is a safe and feasible technique for the treatment of extensive aortoiliac occlusive disease with good three year results regarding patency and clinical improvement. △ Less

Submitted 11 June, 2018; originally announced June 2018.

Evaporation-triggered segregation of sessile binary droplets

Authors: Yaxing Li, Pengyu Lv, Christian Diddens, Huanshu Tan, Herman Wijshoff, Michel Versluis, Detlef Lohse

Abstract: Droplet evaporation of multicomponent droplets is essential for various physiochemical applications, e.g. in inkjet printing, spray cooling and microfabrication. In this work, we observe and study phase segregation of an evaporating sessile binary droplet, consisting of a mixture of water and a surfactant-like liquid (1,2-hexanediol). The phase segregation (i.e., demixing) leads to a reduced water… ▽ More Droplet evaporation of multicomponent droplets is essential for various physiochemical applications, e.g. in inkjet printing, spray cooling and microfabrication. In this work, we observe and study phase segregation of an evaporating sessile binary droplet, consisting of a mixture of water and a surfactant-like liquid (1,2-hexanediol). The phase segregation (i.e., demixing) leads to a reduced water evaporation rate of the droplet and eventually the evaporation process ceases due to shielding of the water by the non-volatile 1,2-hexanediol. Visualizations of the flow field by particle image velocimetry and numerical simulations reveal that the timescale of water evaporation at the droplet rim is faster than that of the Marangoni flow, which originates from the surface tension difference between water and 1,2-hexanediol, eventually leading to segregation. △ Less

Submitted 30 January, 2018; originally announced January 2018.

Journal ref: Phys. Rev. Lett. 120, 224501 (2018)

Apparatus to control and visualize the impact of a high-energy laser pulse on a liquid target

Authors: Alexander L. Klein, Detlef Lohse, Michel Versluis, Hanneke Gelderblom

Abstract: We present an experimental apparatus to control and visualize the response of a liquid target to a laser-induced vaporization. We use a millimeter-sized drop as target and present two liquid-dye solutions that allow a variation of the absorption coefficient of the laser light in the drop by seven orders of magnitude. The excitation source is a Q-switched Nd:YAG laser at its frequency-doubled wavel… ▽ More We present an experimental apparatus to control and visualize the response of a liquid target to a laser-induced vaporization. We use a millimeter-sized drop as target and present two liquid-dye solutions that allow a variation of the absorption coefficient of the laser light in the drop by seven orders of magnitude. The excitation source is a Q-switched Nd:YAG laser at its frequency-doubled wavelength emitting nanosecond pulses with energy densities above the local vaporization threshold. The absorption of the laser energy leads to a large-scale liquid motion at timescales that are separated by several orders of magnitude, which we spatiotemporally resolve by a combination of ultra-high-speed and stroboscopic high-resolution imaging in two orthogonal views. Surprisingly, the large-scale liquid motion at upon laser impact is completely controlled by the spatial energy distribution obtained by a precise beam-shaping technique. The apparatus demonstrates the potential for accurate and quantitative studies of laser-matter interactions. △ Less

Submitted 7 September, 2017; originally announced September 2017.

Comments: Submitted to Review of Scientific Instruments

Evaporating pure, binary and ternary droplets: thermal effects and axial symmetry breaking

Authors: Christian Diddens, Huanshu Tan, Pengyu Lv, Michel Versluis, J. G. M. Kuerten, Xuehua Zhang, Detlef Lohse

Abstract: The Greek aperitif Ouzo is not only famous for its specific anise-flavored taste, but also for its ability to turn from a transparent miscible liquid to a milky-white colored emulsion when water is added. Recently, it has been shown that this so-called Ouzo effect, i.e. the spontaneous emulsification of oil microdroplets, can also be triggered by the preferential evaporation of ethanol in an evapo… ▽ More The Greek aperitif Ouzo is not only famous for its specific anise-flavored taste, but also for its ability to turn from a transparent miscible liquid to a milky-white colored emulsion when water is added. Recently, it has been shown that this so-called Ouzo effect, i.e. the spontaneous emulsification of oil microdroplets, can also be triggered by the preferential evaporation of ethanol in an evaporating sessile Ouzo drop, leading to an amazingly rich drying process with multiple phase transitions [H. Tan et al., Proc. Natl. Acad. Sci. USA 113(31) (2016) 8642]. Due to the enhanced evaporation near the contact line, the nucleation of oil droplets starts at the rim which results in an oil ring encircling the drop. Furthermore, the oil droplets are advected through the Ouzo drop by a fast solutal Marangoni flow. In this article, we investigate the evaporation of mixture droplets in more detail, by successively increasing the mixture complexity from pure water over a binary water-ethanol mixture to the ternary Ouzo mixture (water, ethanol and anise oil). In particular, axisymmetric and full three-dimensional finite element method simulations have been performed on these droplets to discuss thermal effects and the complicated flow in the droplet driven by an interplay of preferential evaporation, evaporative cooling and solutal and thermal Marangoni flow. By using image analysis techniques and micro-PIV measurements, we are able to compare the numerically predicted volume evolutions and velocity fields with experimental data. The Ouzo droplet is furthermore investigated by confocal microscopy. It is shown that the oil ring predominantly emerges due to coalescence. △ Less

Submitted 21 June, 2017; originally announced June 2017.

Journal ref: J. Fluid Mech. 823, 470-497, (2017)

Self-wrapping of an ouzo drop induced by evaporation on a superamphiphobic surface

Authors: Huanshu Tan, Christian Diddens, Michel Versluis, Hans-Jürgen Butt, Detlef Lohse, Xuehua Zhang

Abstract: Evaporation of multi-component drops is crucial to various technologies and has numerous potential applications because of its ubiquity in nature. Superamphiphobic surfaces, which are both superhydrophobic and superoleophobic, can give a low wettability not only for water drops but also for oil drops. In this paper, we experimentally, numerically and theoretically investigate the evaporation proce… ▽ More Evaporation of multi-component drops is crucial to various technologies and has numerous potential applications because of its ubiquity in nature. Superamphiphobic surfaces, which are both superhydrophobic and superoleophobic, can give a low wettability not only for water drops but also for oil drops. In this paper, we experimentally, numerically and theoretically investigate the evaporation process of millimetric sessile ouzo drops (a transparent mixture of water, ethanol, and trans-anethole) with low wettability on a superamphiphobic surface. The evaporation-triggered ouzo effect, i.e. the spontaneous emulsification of oil microdroplets below a specific ethanol concentration, preferentially occurs at the apex of the drop due to the evaporation flux distribution and volatility difference between water and ethanol. This observation is also reproduced by numerical simulations. The volume decrease of the ouzo drop is characterized by two distinct slopes. The initial steep slope is dominantly caused by the evaporation of ethanol, followed by the slower evaporation of water. At later stages, thanks to Marangoni forces the oil wraps around the drop and an oil shell forms. We propose an approximate diffusion model for the drying characteristics, which predicts the evaporation of the drops in agreement with experiment and numerical simulation results. This work provides an advanced understanding of the evaporation process of ouzo (multi-component) drops. △ Less

Submitted 5 June, 2017; originally announced June 2017.

Comments: 41 pages, 8 figures

Oblique drop impact onto a deep liquid pool

Authors: Marise V. Gielen, Pascal Sleutel, Jos Benschop, Michel Riepen, Victoria Voronina, Claas Willem Visser, Detlef Lohse, Jacco H. Snoeijer, Michel Versluis, Hanneke Gelderblom

Abstract: Oblique impact of drops onto a solid or liquid surface is frequently observed in nature. Most studies on drop impact and splashing, however, focus on perpendicular impact. Here, we study oblique impact of 100-micrometer drops onto a deep liquid pool, where we quantify the splashing threshold, maximum cavity dimensions and cavity collapse by high-speed imaging above and below the water surface. Gra… ▽ More Oblique impact of drops onto a solid or liquid surface is frequently observed in nature. Most studies on drop impact and splashing, however, focus on perpendicular impact. Here, we study oblique impact of 100-micrometer drops onto a deep liquid pool, where we quantify the splashing threshold, maximum cavity dimensions and cavity collapse by high-speed imaging above and below the water surface. Gravity can be neglected in these experiments. Three different impact regimes are identified: smooth deposition onto the pool, splashing in the direction of impact only, and splashing in all directions. We provide scaling arguments that delineate these regimes by accounting for the drop impact angle and Weber number. The angle of the axis of the cavity created below the water surface follows the impact angle of the drop irrespectively of the Weber number, while the cavity depth and its displacement with respect to the impact position do depend on the Weber number. Weber number dependency of both the cavity depth and displacement is modeled using an energy argument. △ Less

Submitted 25 July, 2017; v1 submitted 17 November, 2016; originally announced November 2016.

Journal ref: Phys. Rev. Fluids 2, 083602 (2017)

Taming Acoustic Cavitation

Authors: David Fernandez Rivas, Bram Verhaagen, Oscar R. Enriquez, Michel Versluis, Andrea Prosperetti, Han Gardeniers, Detlef Lohse

Abstract: In this fluid dynamics video we show acoustic cavitation occurring from pits etched on a silicon surface. By immersing the surface in a liquid, gas pockets are entrapped in the pits which upon ultrasonic insonation, are observed to shed cavitation bubbles. Modulating the driving pressure it is possible to induce different behaviours based on the force balance that determines the interaction among… ▽ More In this fluid dynamics video we show acoustic cavitation occurring from pits etched on a silicon surface. By immersing the surface in a liquid, gas pockets are entrapped in the pits which upon ultrasonic insonation, are observed to shed cavitation bubbles. Modulating the driving pressure it is possible to induce different behaviours based on the force balance that determines the interaction among bubbles and the silicon surface. This system can be used for several applications like sonochemical water treatment, cleaning of surfaces with deposited materials such as biofilms. △ Less

Submitted 12 October, 2012; originally announced October 2012.

Bubble size prediction in co-flowing streams

Authors: Wim van Hoeve, Benjamin Dollet, José M. Gordillo, Michel Versluis, Detlef Lohse

Abstract: In this paper, the size of bubbles formed through the breakup of a gaseous jet in a co-axial microfluidic device is derived. The gaseous jet surrounded by a co-flowing liquid stream breaks up into monodisperse microbubbles and the size of the bubbles is determined by the radius of the inner gas jet and the bubble formation frequency. We obtain the radius of the gas jet by solving the Navier-Stokes… ▽ More In this paper, the size of bubbles formed through the breakup of a gaseous jet in a co-axial microfluidic device is derived. The gaseous jet surrounded by a co-flowing liquid stream breaks up into monodisperse microbubbles and the size of the bubbles is determined by the radius of the inner gas jet and the bubble formation frequency. We obtain the radius of the gas jet by solving the Navier-Stokes equations for low Reynolds number flows and by minimization of the dissipation energy. The prediction of the bubble size is based on the system's control parameters only, i.e. the inner gas flow rate $Q_i$, the outer liquid flow rate $Q_o$, and the tube radius $R$. For a very low gas-to-liquid flow rate ratio ($Q_i / Q_o \rightarrow 0$) the bubble radius scales as $r_b / R \propto \sqrt{Q_i / Q_o}$, independently of the inner to outer viscosity ratio $η_i/η_o$ and of the type of the velocity profile in the gas, which can be either flat or parabolic, depending on whether high-molecular-weight surfactants cover the gas-liquid interface or not. However, in the case in which the gas velocity profiles are parabolic and the viscosity ratio is sufficiently low, i.e. $η_i/η_o \ll 1$, the bubble diameter scales as $r_b \propto (Q_i/Q_o)^β$, with $β$ smaller than 1/2. △ Less

Submitted 1 March, 2011; originally announced March 2011.

Journal ref: Europhys. Lett. 94, 64001 (2011)

Microbubble formation and pinch-off scaling exponent in flow-focusing devices

Authors: Wim van Hoeve, Benjamin Dollet, Michel Versluis, Detlef Lohse

Abstract: We investigate the gas jet breakup and the resulting microbubble formation in a microfluidic flow-focusing device using ultra high-speed imaging at 1 million frames/s. In recent experiments [Dollet et al., Phys. Rev. Lett. 100, 034504 (2008)] it was found that in the final stage of the collapse the radius of the neck scales with time with a 1/3 power-law exponent, which suggested that gas inertia… ▽ More We investigate the gas jet breakup and the resulting microbubble formation in a microfluidic flow-focusing device using ultra high-speed imaging at 1 million frames/s. In recent experiments [Dollet et al., Phys. Rev. Lett. 100, 034504 (2008)] it was found that in the final stage of the collapse the radius of the neck scales with time with a 1/3 power-law exponent, which suggested that gas inertia and the Bernoulli suction effect become important. Here, ultra high-speed imaging was used to capture the complete bubble contour and quantify the gas flow through the neck. It revealed that the resulting decrease in pressure, due to Bernoulli suction, is too low to account for an accelerated pinch-off. The high temporal resolution images enable us to approach the final moment of pinch-off to within 1 μs. We observe that the final moment of bubble pinch-off is characterized by a scaling exponent of 0.41 +/- 0.01. This exponent is approximately 2/5, which can be derived, based on the observation that during the collapse the neck becomes less slender, due to the exclusive driving through liquid inertia. △ Less

Submitted 28 February, 2011; originally announced February 2011.

Breakup of diminutive Rayleigh jets

Authors: Wim van Hoeve, Stephan Gekle, Jacco H. Snoeijer, Michel Versluis, Michael P. Brenner, Detlef Lohse

Abstract: Discharging a liquid from a nozzle at sufficient large velocity leads to a continuous jet that due to capillary forces breaks up into droplets. Here we investigate the formation of microdroplets from the breakup of micron-sized jets with ultra high-speed imaging. The diminutive size of the jet implies a fast breakup time scale $τ_\mathrm{c} = \sqrt{ρr^3 / γ}$ of the order of 100\,ns{}, and require… ▽ More Discharging a liquid from a nozzle at sufficient large velocity leads to a continuous jet that due to capillary forces breaks up into droplets. Here we investigate the formation of microdroplets from the breakup of micron-sized jets with ultra high-speed imaging. The diminutive size of the jet implies a fast breakup time scale $τ_\mathrm{c} = \sqrt{ρr^3 / γ}$ of the order of 100\,ns{}, and requires imaging at 14 million frames per second. We directly compare these experiments with a numerical lubrication approximation model that incorporates inertia, surface tension, and viscosity [Eggers and Dupont, J. Fluid Mech. 262, 205 (1994); Shi, Brenner, and Nagel, Science 265, 219 (1994)]. The lubrication model allows to efficiently explore the parameter space to investigate the effect of jet velocity and liquid viscosity on the formation of satellite droplets. In the phase diagram we identify regions where the formation of satellite droplets is suppressed. We compare the shape of the droplet at pinch-off between the lubrication approximation model and a boundary integral (BI) calculation, showing deviations at the final moment of the pinch-off. Inspite of this discrepancy, the results on pinch-off times and droplet and satellite droplet velocity obtained from the lubrication approximation agree with the high-speed imaging results. △ Less

Submitted 1 November, 2010; originally announced November 2010.

Journal ref: Phys. Fluids 22, 122003 (2010)

Collapse of Non-Axisymmetric Cavities

Authors: Oscar R. Enriquez, Ivo R. Peters, Stephan Gekle, Laura Schmidt, Michel Versluis, Devaraj van der Meer, Detlef Lohse

Abstract: A round disk with a harmonic disturbance impacts on a water surface and creates a non-axisymmetric cavity which collapses under the influence of hydrostatic pressure. We use disks deformed with mode m=2 to m=6. For all mode numbers we find clear evidence for a phase inversion of the cavity wall during the collapse. We present a fluid dynamics video showing high speed imaging of different modes,… ▽ More A round disk with a harmonic disturbance impacts on a water surface and creates a non-axisymmetric cavity which collapses under the influence of hydrostatic pressure. We use disks deformed with mode m=2 to m=6. For all mode numbers we find clear evidence for a phase inversion of the cavity wall during the collapse. We present a fluid dynamics video showing high speed imaging of different modes, pointing out the characteristic features during collapse. △ Less

Submitted 14 October, 2009; originally announced October 2009.

On the sound of snapping shrimp

Authors: Michel Versluis, Barbara Schmitz, Anna von der Heydt, Detlef Lohse

Abstract: Fluid dynamics video: Snapping shrimp produce a snapping sound by an extremely rapid closure of their snapper claw. Our high speed imaging of the claw closure has revealed that the sound is generated by the collapse of a cavitation bubble formed in a fast flowing water jet forced out from the claws during claw closure. The produced sound originates from the cavitation collapse of the bubble. At… ▽ More Fluid dynamics video: Snapping shrimp produce a snapping sound by an extremely rapid closure of their snapper claw. Our high speed imaging of the claw closure has revealed that the sound is generated by the collapse of a cavitation bubble formed in a fast flowing water jet forced out from the claws during claw closure. The produced sound originates from the cavitation collapse of the bubble. At collapse a short flash of light is emitted, just as in single bubble sonoluminescence. A model based on the Rayleigh-Plesset equation can quantitatively account for the visual and acoustical observations. △ Less

Submitted 13 December, 2007; originally announced December 2007.

Comments: Fluid dynamics video

Leaping shampoo and the stable Kaye effect

Authors: Michel Versluis, Cor Blom, Devaraj van der Meer, Ko van der Weele, Detlef Lohse

Abstract: Shear-thinning fluids exhibit surprisingly rich behaviour. One example is the Kaye effect which occurs when a thin stream of a solution of polyisobutylene in Decalin is poured into a dish of the fluid. As pouring proceeds, a small stream of liquid occasionally leaps upward from the heap. This surprising effect, which lasts only a second or so, is named after its first observer A. Kaye, who could… ▽ More Shear-thinning fluids exhibit surprisingly rich behaviour. One example is the Kaye effect which occurs when a thin stream of a solution of polyisobutylene in Decalin is poured into a dish of the fluid. As pouring proceeds, a small stream of liquid occasionally leaps upward from the heap. This surprising effect, which lasts only a second or so, is named after its first observer A. Kaye, who could offer no explanation for this behaviour. Later, Collyer and Fischer suggested from 250 frames per second cine recordings that the fluid must be highly shear thinning as well as elastic and 'pituitous'. In addition, they concluded that a rigid surface is required to back the reflected liquid stream. While the words bouncing and reflection are associated with non-continuous and elastic effects, we will show here that the Kaye effect is in fact a continuous flow phenomenon. We show that the Kaye effect works for many common fluids, including shampoos and liquid soaps. We reveal its physical mechanism (formation, stability and disruption) through high-speed imaging. The measurements are interpreted with a simple theoretical model including only the shear thinning behaviour of the liquid; elastic properties of the liquid play no role. We show that the Kaye effect can be stable and that it can be directed. We even demonstrate a stable Kaye effect on a thin soap film excluding the necessity of a rigid backing surface. △ Less

Submitted 22 March, 2006; originally announced March 2006.

Comments: 24 pages, 7 figures

Journal ref: J. Stat. Mech. (2006) P07007