Max Planck University of Twente Center for Complex Fluid Dynamics

Inkjet printing is the most widespread technological application of microfluidics. It is characterized by its high drop productivity, small volumes, and extreme reproducibility. This review gives a synopsis of the fluid dynamics of inkjet printing and discusses the main challenges for present and future research. These lie both on the printhead side—namely, the detailed flow inside the printhead, entrained bubbles, the meniscus dynamics, wetting phenomena at the nozzle plate, and jet formation—and on the receiving substrate side—namely, droplet impact, merging, wetting of the substrate, droplet evaporation, and drying. In most cases the droplets are multicomponent, displaying rich physicochemical hydrodynamic phenomena. The challenges on the printhead side and on the receiving substrate side are interwoven, as optimizing the process and the materials with respect to either side alone is not enough: As the same ink (or other jetted liquid) is used and as droplet frequency and size matter on both sides, the process must be optimized as a whole.

The snapping shrimp (Alpheus heterochaelis) produces a loud snapping sound by an extremely rapid closure of its snapper claw. One of the effects of the snapping is to stun or kill prey animals. During the rapid snapper claw closure, a high-velocity water jet is emitted from the claw with a speed exceeding cavitation conditions. Hydrophone measurements in conjunction with time-controlled high-speed imaging of the claw closure demonstrate that the sound is emitted at the cavitation bubble collapse and not on claw closure. A model for the bubble dynamics based on a Rayleigh-Plesset-type equation quantitatively accounts for the time dependence of the bubble radius and for the emitted sound.

Taylor-Couette flow, the flow between two coaxial co-or counter-rotating cylinders, is one of the paradigmatic systems in the physics of fluids. The (dimensionless) control parameters are the Reynolds numbers of the inner and outer cylinders, the ratio of the cylinder radii, and the aspect ratio. One key response of the system is the torque required to retain constant angular velocities, which can be connected to the angular velocity transport through the gap. Whereas the low-Reynolds number regime was well explored in the 1980s and 1990s of the past century, in the fully turbulent regime major research activity developed only in the past decade. In this article, we review this recent progress in our understanding of fully developed Taylor-Couette turbulence from the experimental, numerical, and theoretical points of view. We focus on the parameter dependence of the global torque and on the local flow organization, including velocity profiles and boundary layers. Next, we discuss transitions between different (turbulent) flow states. We also elaborate on the relevance of this system for astrophysical disks (quasi-Keplerian flows). The review ends with a list of challenges for future research on turbulent Taylor-Couette flow.

Organs-on-chips are microfluidic systems with controlled, dynamic microenvironments in which cultured cells exhibit functions that emulate organ-level physiology. They can in principle be 'personalised' to reflect individual physiology, for example by including blood samples, primary human tissue, and cells derived from induced pluripotent stem cell-derived cells, as well as by tuning key physico-chemical parameters of the cell culture microenvironment based on personal health data. The personalised nature of such systems, combined with physiologically relevant read-outs, provides new opportunities for person-specific assessment of drug efficacy and safety, as well as personalised strategies for disease prevention and treatment; together, this is known as 'precision medicine'. There are multiple reports of how to personalise organs-on-chips, with examples including airway-on-a-chip systems containing primary patient alveolar epithelial cells, vessels-on-chips with shapes based on personal biomedical imaging data and lung-on-a-chip systems that can be exposed to various regimes of cigarette smoking. In addition, multi-organ chip systems even allow the systematic and dynamic integration of more complex combinations of personalised cell culture parameters. Current personalised organs-on-chips have not yet been used for precision medicine as such. The major challenges that affect the implementation of personalised organs-on-chips in precision medicine are related to obtaining access to personal samples and corresponding health data, as well as to obtaining data on patient outcomes that can confirm the predictive value of personalised organs-on-chips. We argue here that involving all biomedical stakeholders from clinicians and patients to pharmaceutical companies will be integral to transition personalised organs-on-chips to precision medicine.

Droplets abound in nature and technology. In general, they are multicomponent, and, when out of equilibrium, have gradients in concentration, implying flow and mass transport. Moreover, phase transitions can occur, in the form of evaporation, solidification, dissolution or nucleation of a new phase. The droplets and their surrounding liquid can be binary, ternary or contain even more components, with several in different phases. Since the early 2000s, rapid advances in experimental and numerical fluid dynamical techniques have enabled major progress in our understanding of the physicochemical hydrodynamics of such droplets, further narrowing the gap from fluid dynamics to chemical engineering and colloid and interfacial science, arriving at a quantitative understanding of multicomponent and multiphase droplet systems far from equilibrium, and aiming towards a one-to-one comparison between experiments and theory or numerics. This Perspective discusses examples of the physicochemical hydrodynamics of droplet systems far from equilibrium and the relevance of such systems for applications. Droplets in general are multicomponent and experience gradients in concentration, often leading to transport phenomena and phase transitions. This Perspective discusses recent progress on the physicochemical hydrodynamics of such droplet systems and their relevance for many important applications.

Very different types of scaling of the Nusselt number Nu with the Rayleigh number Ra have experimentally been found in the very large Ra regime beyond 1011. We understand and interpret these results by extending the unifying theory of thermal convection [Grossmann and Lohse, Phys. Rev. Lett. 86, 3316 (2001)] to the very large Ra regime where the kinetic boundary-layer is turbulent. The central idea is that the spatial extension of this turbulent boundary-layer with a logarithmic velocity profile is comparable to the size of the cell. Depending on whether the thermal transport is plume dominated, dominated by the background thermal fluctuations, or whether also the thermal boundary-layer is fully turbulent (leading to a logarithmic temperature profile), we obtain effective scaling laws of about Nu∝Ra0.14, Nu∝Ra0.22, and Nu∝Ra0.38, respectively. Depending on the initial conditions or random fluctuations, one or the other of these states may be realized. Since the theory is for both the heat flux Nu and the velocity amplitude Re, we can also give the scaling of the latter, namely, Re∝Ra0.42, Re∝Ra0.45, and Re∝Ra0.50 in the respective ranges.

Very fine sand is prepared in a well-defined and fully decompactified state by letting gas bubble through it. After turning off the gas stream, a steel ball is dropped on the sand. On impact of the ball, sand is blown away in all directions ("splash") and an impact crater forms. When this cavity collapses, a granular jet emerges and is driven straight into the air. A second jet goes downwards into the air bubble entrained during the process, thus pushing surface material deep into the ground. The air bubble rises slowly towards the surface, causing a granular eruption. In addition to the experiments and the discrete particle simulations we present a simple continuum theory to account for the void collapse leading to the formation of the upward and downward jets.

A paperclip can float on water. Drops of mercury do not spread on a surface. These capillary phenomena are macroscopic manifestations of molecular interactions and can be explained in terms of surface tension. We address several conceptual questions that are often encountered when teaching capillarity and provide a perspective that reconciles the macroscopic viewpoints from thermodynamics and fluid mechanics and the microscopic perspective from statistical physics.

To quantify the fate of respiratory droplets under different ambient relative humidities, direct numerical simulations of a typical respiratory event are performed. We found that, because small droplets (with initial diameter of 10 μm) are swept by turbulent eddies in the expelled humid puff, their lifetime gets extended by a factor of more than 30 times as compared to what is suggested by the classical picture by Wells, for 50% relative humidity. With increasing ambient relative humidity the extension of the lifetimes of the small droplets further increases and goes up to around 150 times for 90% relative humidity, implying more than 2 m advection range of the respiratory droplets within 1 sec. Employing Lagrangian statistics, we demonstrate that the turbulent humid respiratory puff engulfs the small droplets, leading to many orders of magnitude increase in their lifetimes, implying that they can be transported much further during the respiratory events than the large ones. Our findings provide the starting points for larger parameter studies and may be instructive for developing strategies on optimizing ventilation and indoor humidity control. Such strategies are key in mitigating the COVID-19 pandemic in the present autumn and upcoming winter.

This paper sketches my personal scientific bubble journey, starting with single-bubble sonoluminescence, continuing with sound emission and scattering of bubbles, cavitation, snapping shrimp, impact events, air entrainment, and surface micro- and nanobubbles, and finally arriving at effective force models for bubbles and dispersed bubbly two-phase flow.

Fluid turbulence is commonly associated with stronger drag, greater heat transfer, and more efficient mixing than in laminar flows. In many natural and industrial settings, turbulent liquid flows contain suspensions of dispersed bubbles and light particles. Recently, much attention has been devoted to understanding the behavior and underlying physics of such flows by use of both experiments and high-resolution direct numerical simulations. This review summarizes our present understanding of various phenomenological aspects of bubbly and buoyant particle–laden turbulent flows. We begin by discussing different dynamical regimes, including those of crossing trajectories and wake-induced oscillations of rising particles, and regimes in which bubbles and particles preferentially accumulate near walls or within vortical structures. We then address how certain paradigmatic turbulent flows, such as homogeneous isotropic turbulence, channel flow, Taylor–Couette turbulence, and thermally driven turbulence, are modified by the presence of these dispersed bubbles and buoyant particles. We end with a list of summary points and future research questions.

Here, we describe methods for combining impedance spectroscopy measurements with electrical simulation to reveal transepithelial barrier function and tissue structure of human intestinal epithelium cultured inside an organ-on-chip microfluidic culture device. When performing impedance spectroscopy measurements, electrical simulation enabled normalization of cell layer resistance of epithelium cultured statically in a gut-on-a-chip, which enabled determination of transepithelial electrical resistance (TEER) values that can be compared across device platforms. During culture under dynamic flow, the formation of intestinal villi was accompanied by characteristic changes in impedance spectra both measured experimentally and verified with simulation, and we demonstrate that changes in cell layer capacitance may serve as measures of villi differentiation. This method for combining impedance spectroscopy with simulation can be adapted to better monitor cell layer characteristics within any organ-on-chip in vitro and to enable direct quantitative TEER comparisons between organ-on-chip platforms which should help to advance research on organ function.

The evolution of electrogenerated gas bubbles during water electrolysis can significantly hamper the overall process efficiency. Promoting the departure of electrochemically generated bubbles during (water) electrolysis is therefore beneficial. For a single bubble, a departure from the electrode surface occurs when buoyancy wins over the downward-acting forces (e.g., contact, Marangoni, and electric forces). In this work, the dynamics of a pair of H2 bubbles produced during the hydrogen evolution reaction in 0.5 M H2SO4 using a dual platinum microelectrode system is systematically studied by varying the electrode distance and the cathodic potential. By combining high-speed imaging and electrochemical analysis, we demonstrate the importance of bubble–bubble interactions in the departure process. We show that bubble coalescence may lead to substantially earlier bubble departure as compared to buoyancy effects alone, resulting in considerably higher reaction rates at a constant potential. However, due to continued mass input and conservation of momentum, repeated coalescence events with bubbles close to the electrode may drive departed bubbles back to the surface beyond a critical current, which increases with the electrode spacing. The latter leads to the resumption of bubble growth near the electrode surface, followed by buoyancy-driven departure. While less favorable at small electrode spacing, this configuration proves to be very beneficial at larger separations, increasing the mean current up to 2.4 times compared to a single electrode under the conditions explored in this study.

Abstract We present results from large eddy simulations of extended wind‐farms for several turbine configurations with a range of different spanwise and streamwise spacing combinations. The results show that for wind‐farms arranged in a staggered configuration with spanwise spacings in the range ≈[3.5,8]D, where D is the turbine diameter, the power output in the fully developed regime depends primarily on the geometric mean of the spanwise and streamwise turbine spacings. In contrast, for the aligned configuration the power output in the fully developed regime strongly depends on the streamwise turbine spacing and shows weak dependence on the spanwise spacing. Of interest to the rate of wake recovery, we find that the power output is well correlated with the vertical kinetic energy flux, which is a measure of how much kinetic energy is transferred into the wind‐turbine region by the mean flow. A comparison between the aligned and staggered configurations reveals that the vertical kinetic energy flux is more localized along turbine columns for aligned wind‐farms than for staggered ones. This additional mixing leads to a relatively fast wake recovery for aligned wind‐farms. Copyright © 2015 John Wiley & Sons, Ltd.

A falling liquid drop, after impact on a rigid substrate, deforms and spreads, owing to the normal reaction force. Subsequently, if the substrate is nonwetting, the drop retracts and then jumps off. As we show here, not only is the impact itself associated with a distinct peak in the temporal evolution of the normal force, but also the jump-off, which was hitherto unknown. We characterize both peaks and elucidate how they relate to the different stages of the drop impact process. The time at which the second peak appears coincides with the formation of a Worthington jet, emerging through flow focusing. Even low-velocity impacts can lead to a surprisingly high second peak in the normal force, even larger than the first one, namely when the Worthington jet becomes singular due to the collapse of an air cavity in the drop.

The possible transition to the so-called ultimate regime, wherein both the bulk and the boundary layers are turbulent, has been an outstanding issue in thermal convection, since the seminal work by Kraichnan [Phys. Fluids 5, 1374 (1962)PFLDAS0031-917110.1063/1.1706533]. Yet, when this transition takes place and how the local flow induces it is not fully understood. Here, by performing two-dimensional simulations of Rayleigh-Bénard turbulence covering six decades in Rayleigh number Ra up to 10^{14} for Prandtl number Pr=1, for the first time in numerical simulations we find the transition to the ultimate regime, namely, at Ra^{*}=10^{13}. We reveal how the emission of thermal plumes enhances the global heat transport, leading to a steeper increase of the Nusselt number than the classical Malkus scaling Nu∼Ra^{1/3} [Proc. R. Soc. A 225, 196 (1954)PRLAAZ1364-502110.1098/rspa.1954.0197]. Beyond the transition, the mean velocity profiles are logarithmic throughout, indicating turbulent boundary layers. In contrast, the temperature profiles are only locally logarithmic, namely, within the regions where plumes are emitted, and where the local Nusselt number has an effective scaling Nu∼Ra^{0.38}, corresponding to the effective scaling in the ultimate regime.

Photonic crystals (PhCs) display photonic stop bands (PSBs) and at the edges of these PSBs transport light with reduced velocity, enabling the PhCs to confine and manipulate incident light with enhanced light-matter interaction. Intense research has been devoted to leveraging the optical properties of PhCs for the development of optical sensors for bioassays, diagnosis, and environmental monitoring. These applications have furthermore benefited from the inherently large surface area of PhCs, giving rise to high analyte adsorption and the wide range of options for structural variations of the PhCs leading to enhanced light-matter interaction. Here, we focus on bottom-up assembled PhCs and review the significant advances that have been made in their use as label-free sensors. We describe their potential for point-of-care devices and in the review include their structural design, constituent materials, fabrication strategy, and sensing working principles. We thereby classify them according to five sensing principles: sensing of refractive index variations, sensing by lattice spacing variations, enhanced fluorescence spectroscopy, surface-enhanced Raman spectroscopy, and configuration transitions.

Direct numerical simulations of highly turbulent Rayleigh-B\'enard convection up to Rayleigh numbers Ra = ${10}^{9}$ in horizontally periodic domains with aspect ratios up to \ensuremath{\Gamma} = 128, find superstructures of very large-scale and long living coherent structures.

The development of needle-free drug injection systems is of great importance to global healthcare. However, in spite of its great potential and research history over many decades, these systems are not commonly used. One of the main problems is that existing methods use diffusive jets, which result in scattered penetration and severe deceleration of the jets, causing frequent pain and insufficient penetration. Another long-standing challenge is the development of accurate small volume injections. In this paper we employ a novel method of needle-free drug injection, using highly-focused high speed microjets, which aims to solve these challenges. We experimentally demonstrate that these unique jets are able to penetrate human skin: the focused nature of these microjets creates an injection spot smaller than a mosquito's proboscis and guarantees a high percentage of the liquid being injected. The liquid substances can be delivered to a much larger depth than conventional methods, and create a well-controlled dispersion pattern. Thanks to the excellent controllability of the microjet, small volume injections become feasible. Furthermore, the penetration dynamics is studied through experiments performed on gelatin mixtures (human soft tissue equivalent) and human skin, agreeing well with a viscous stress model which we develop. This model predicts the depth of the penetration into both human skin and soft tissue. The results presented here take needle-free injections a step closer to widespread use.