Max Planck University of Twente Center for Complex Fluid Dynamics

Wakata Y., Wang F., Sun C., Lohse D.

Volatile multicomponent liquid films show rich dynamics, due to the complex interplay of gradients in temperature and in solute concentrations. Here, we study the evaporation dynamics of a tricomponent liquid film, consisting of water, ethanol, and trans-anethole oil (known as “ouzo”). With the preferential evaporation of ethanol, cellular convective structures are observed both in the thermal patterns and in the nucleated oil droplet patterns. However, the feature sizes of these two patterns can differ, indicating dual instability mechanisms dominated by either temperature or solute concentration. Using numerical simulations, we quantitatively compare the contributions of temperature and solute concentration on the surface tension. Our results reveal that the thermal Marangoni effect predominates at the initial evaporation stage, resulting in cellular patterns in thermal images, while the solutal Marangoni effect gradually becomes dominant. By regulating the transition time of this thermal-solutal-induced bistability and the nucleation time of oil microdroplets in the ternary mixture, the oil droplet patterns can be well controlled. This capability not only enhances our understanding of the evaporation dynamics but also paves the way for precise manipulation of nucleation and deposition processes at larger scales.

Ren X., Xu D., Wang J., Chen S.

Diddens C., Rocha D.

We present a black-box method to numerically investigate the linear stability of arbitrary multi-physics problems. While the user just has to enter the system's residual in weak formulation, e.g. by a finite element method, all required discretized matrices are automatically assembled based on just-in-time generated and compiled C codes. Based on this method, entire phase diagrams in the parameter space can be obtained by bifurcation tracking and continuation at low computational costs. Particular focus is put on problems with moving domains, e.g. free surface problems in fluid dynamics, since a moving mesh introduces a plethora of complicated nonlinearities to the system. By symbolic differentiation before the code generation, however, these moving mesh problems are made accessible to bifurcation tracking methods. In a second step, our method is generalized to investigate symmetry-breaking instabilities of axisymmetric stationary solutions by effectively utilizing the symmetry of the base state. Each bifurcation type is validated on the basis of results reported in the literature on versatile fluid dynamics problems, for which we subsequently present novel results as well.

Yerragolam G.S., Howland C.J., Stevens R.J., Verzicco R., Shishkina O., Lohse D.

We provide scaling relations for the Nusselt number $Nu$ and the friction coefficient $C_{S}$ in sheared Rayleigh–Bénard convection, i.e. in Rayleigh–Bénard flow with Couette- or Poiseuille-type shear forcing, by extending the Grossmann & Lohse (J. Fluid Mech., vol. 407, 2000, pp. 27–56, Phys. Rev. Lett., vol. 86, 2001, pp. 3316–3319, Phys. Rev. E, vol. 66, 2002, 016305, Phys. Fluids, vol. 16, 2004, pp. 4462–4472) theory to sheared thermal convection. The control parameters for these systems are the Rayleigh number $Ra$ , the Prandtl number $Pr$ and the Reynolds number $Re_S$ that characterises the strength of the imposed shear. By direct numerical simulations and theoretical considerations, we show that, in turbulent Rayleigh–Bénard convection, the friction coefficients associated with the applied shear and the shear generated by the large-scale convection rolls are both well described by Prandtl's (Ergeb. Aerodyn. Vers. Gött., vol. 4, 1932, pp. 18–29) logarithmic friction law, suggesting some kind of universality between purely shear-driven flows and thermal convection. These scaling relations hold well for $10^6 \leq Ra \leq 10^8$ , $0.5 \leq Pr \leq 5.0$ , and $0 \leq Re_S \leq 10^4$ .

van Buuren D., Kant P., Meijer J.G., Diddens C., Lohse D.

A uniform solidification front undergoes nontrivial deformations when encountering an insoluble dispersed particle in a melt. For solid particles, the overall deformation characteristics are primarily dictated by heat transfer between the particle and the surrounding, remaining unaffected by the rate of approach of the solidification front. In this Letter we show that, conversely, when interacting with a droplet or a bubble, the deformation behavior exhibits entirely different and unexpected behavior. It arises from an interfacial dynamics which is specific to particles with free interfaces, namely, thermal Marangoni forces. Our Letter employs a combination of experiments, theory, and numerical simulations to investigate the interaction between the droplet and the freezing front and unveils its surprising behavior. In particular, we quantitatively understand the dependence of the front deformation Δ on the front propagation velocity, set by the strength of the applied thermal gradient, which, for larger front velocities (larger applied thermal gradients), can even revert from attraction (Δ<0) to repulsion (Δ>0). Published by the American Physical Society 2024

Yang R., Howland C.J., Liu H., Verzicco R., Lohse D.

Latent heat storage (LHS) has emerged as a promising solution for addressing the challenges of large-scale and long-term energy storage, offering a clean and reusable system. Being in the developmental stage, and with only limited theoretical predictions being available, there is a need to enhance the efficiency of LHS systems. In this study, we use numerical simulations to study the melting process of a phase-change material (PCM) in a rectangular domain. A significant enhancement in the melt rate of the PCM is found by a simple inclination of the domain, with a 60% increase in melt rate compared to a standard square unit without inclination through enhanced heat transfer. We establish a theoretical relation for the optimal PCM melt rate, based on the inclination angle and the domain aspect ratio, outlining the optimal conditions for the PCM melt rate, which is solely dependent on the domain geometry. We further propose a heat-transfer model that predicts the optimal aspect ratio for achieving the maximum melt rate as a function of the Rayleigh and Prandtl numbers. Published by the American Physical Society 2024

Engelen Y., Krysko D.V., Effimova I., Breckpot K., Versluis M., De Smedt S., Lajoinie G., Lentacker I.

Over the past decade, ultrasound (US) has gathered significant attention and research focus in the realm of medical treatments, particularly within the domain of anti-cancer therapies. This growing interest can be attributed to its non-invasive nature, precision in delivery, availability, and safety. While the conventional objective of US-based treatments to treat breast, prostate, and liver cancer is the ablation of target tissues, the introduction of the concept of immunogenic cell death (ICD) has made clear that inducing cell death can take different non-binary pathways through the activation of the patient's anti-tumor immunity. Here, we investigate high-intensity focused ultrasound (HIFU) to induce ICD by unraveling the underlying physical phenomena and resulting biological effects associated with HIFU therapy using an automated and fully controlled experimental setup. Our in-vitro approach enables the treatment of adherent cancer cells (B16F10 and CT26), analysis for ICD hallmarks and allows to monitor and characterize in real time the US-induced cavitation activity through passive cavitation detection (PCD). We demonstrate HIFU-induced cell death, CRT exposure, HMGB1 secretion and antigen release. This approach holds great promise in advancing our understanding of the therapeutic potential of HIFU for anti-cancer strategies.

Hack M.A., van der Linden M.N., Wijshoff H., Snoeijer J.H., Segers T.

Electrostatically stabilised colloidal particles destabilise when brought into contact with cations causing the particles to aggregate in clusters. When a drop with stabilised colloidal partices is deposited on a liquid film containing cations the delicate balance between the fluid-mechanical and physicochemical properties of the system governs the spreading dynamics and formation of colloidal particle clusters. High-speed imaging and digital holographic microscopy were used to characterise the spreading process. We reveal that a spreading colloidal drop evolves into a ring-shaped pattern after it is deposited on a thin saline water film. Clustered colloidal particles aggregate into larger trapezoidally-shaped 'supraclusters'. Using a simple model we show that the trapezoidal shape of the supraclusters is determined by the transition from inertial spreading dynamics to Marangoni flow. These results may be of interest to applications such as wet-on-wet inkjet printing, where particle destabilisation and hydrodynamic flow coexist.

Yang R., Chong K.L., Liu H., Verzicco R., Lohse D.

Melting and solidification in periodically time-modulated thermal convection are relevant for numerous natural and engineering systems, for example, glacial melting under periodic sun radiation and latent thermal energy storage under periodically pulsating heating. It is highly relevant for the estimation of melt rate and melt efficiency management. However, even the dynamics of a solid–liquid interface shape subjected to a simple sinusoidal heating has not yet been investigated in detail. In this paper, we offer a better understanding of the modulation frequency dependence of the melting and solidification front. We numerically investigate periodic melting and solidification in turbulent convective flow with the solid above and the melted liquid below, and sinusoidal heating at the bottom plate with the mean temperature equal to the melting temperature. We investigate how the periodic heating can prevent the full solidification, and the resulting flow structures and the quasi-equilibrium interface height. We further study the dependence on the heating modulation frequency. As the frequency decreases, we found two distinct regimes, which are ‘partially solid’ and ‘fully solid’. In the fully solid regime, the liquid freezes completely, and the effect of the modulation is limited. In the partially solid regime, the solid partially melts, and a steady or unsteady solid–liquid interface forms depending on the frequency. The interface height can be derived based on the energy balance through the interface. In the partially solid regime, the interface height oscillates periodically, following the frequency of modulation. Here, we propose a perturbation approach that can predict the dependency of the oscillation amplitude on the modulation frequency.

Piumini G., Assen M.P., Lohse D., Verzicco R.

We use three-dimensional direct numerical simulations of homogeneous isotropic turbulence in a cubic domain to investigate the dynamics of heavy, chiral, finite-size inertial particles and their effects on the flow. Using an immersed-boundary method and a complex collision model, four-way coupled simulations have been performed, and the effects of particle-to-fluid density ratio, turbulence strength and particle volume fraction have been analysed. We find that freely falling particles on the one hand add energy to the turbulent flow but, on the other hand, they also enhance the flow dissipation: depending on the combination of flow parameters, the former or the latter mechanism prevails, thus yielding enhanced or weakened turbulence. Furthermore, particle chirality entails a preferential angular velocity which induces a net vorticity in the fluid phase. As turbulence strengthens, the energy introduced by the falling particles becomes less relevant and stronger velocity fluctuations alter the solid phase dynamics, making the effect of chirality irrelevant for the large-scale features of the flow. Moreover, comparing the time history of collision events for chiral particles and spheres (at the same volume fraction) suggests that the former tend to entangle, in contrast to the latter which rebound impulsively.

Harte N.C., Obrist D., Caversaccio M., Lajoinie G.P., Wimmer W.

The cochlea, situated within the inner ear, is a spiral-shaped, liquid-filled organ responsible for hearing. The physiological significance of its shape remains uncertain. Previous research has scarcely addressed the occurrence of transverse flow within the cochlea, particularly in relation to its unique shape. This study aims to investigate the impact of the geometric features of the cochlea on fluid dynamics by characterizing transverse flow induced by harmonically oscillating axial flow in square ducts with curvature and torsion resembling human cochlear anatomy. We examined four geometries to investigate curvature and torsion effects on axial and transverse flow components. Twelve frequencies from 0.125 Hz to 256 Hz were studied, covering infrasound and low-frequency hearing, with mean inlet velocity amplitudes representing levels expected for normal conversation or louder situations. Our simulations show that torsion contributes significantly to transverse flow in unsteady conditions, and that its contribution increases with increasing oscillation frequency. Curvature alone has a small effect on transverse flow strength, which decreases rapidly with increasing frequency. Strikingly, the combined effect of curvature and torsion on transverse flow is greater than expected from a simple superposition of the two effects, especially when the relative contribution of curvature alone becomes negligible. These findings may be relevant to understanding physiological processes in the cochlea, including metabolite transport and wall shear stress. Further studies are needed to investigate possible implications for cochlear mechanics.

Nawijn C.L., Segers T., Lajoinie G., Berg S., Snipstad S., Davies C.D., Versluis M.

Ultrasound-triggered bubble-mediated local drug delivery has shown potential to increase therapeutic efficacy and reduce systemic side effects, by loading drugs into the microbubble shell and triggering delivery of the payload on demand using ultrasound. Understanding the behavior of the microbubbles in response to ultrasound is crucial for efficient and controlled release.

Hou X., Xu D., Wang J., Chen S.

Zhu X., Fu Y., De Paoli M.

We present a theory to describe the Nusselt number, $\operatorname {\mathit {Nu}}$ , corresponding to the heat or mass flux, as a function of the Rayleigh–Darcy number, $\operatorname {\mathit {Ra}}$ , the ratio of buoyant driving force over diffusive dissipation, in convective porous media flows. First, we derive exact relationships within the system for the kinetic energy and the thermal dissipation rate. Second, by segregating the thermal dissipation rate into contributions from the boundary layer and the bulk, which is inspired by the ideas of the Grossmann and Lohse theory (J. Fluid Mech., vol. 407, 2000; Phys. Rev. Lett., vol. 86, 2001), we derive the scaling relation for $\operatorname {\mathit {Nu}}$ as a function of $\operatorname {\mathit {Ra}}$ and provide a robust theoretical explanation for the empirical relations proposed in previous studies. Specifically, by incorporating the length scale of the flow structure into the theory, we demonstrate why heat or mass transport differs between two-dimensional and three-dimensional porous media convection. Our model is in excellent agreement with the data obtained from numerical simulations, affirming its validity and predictive capabilities.

Rudravaram S., Shukla R.P., Maheshwaram S.

Abstract In this work we report a vertically stacked nanosheet Field Effect Transistor (NSFET) in double gate configuration using transition metal dichalcogenide (TMD) based molybdenum disulphide (MoS2) as the conducting channel. The performance of the NSFET is analysed for number of channels, different channel thickness, different source/drain contacts. The performance of the device at different temperatures (T) also analysed. The proposed NSFET with three vertically stacked channels, exhibits a ON current (ION) of 30.6 μA μm−1, Subthreshold swing (SS) of 69 mV/dec and ON to OFF current ratio of more than 108 at Vds = 1V. Further the ION can be improved with multi-layer channel thickness. The performance of the vertically stacked MoS2 NSFET in junction less (JL) and inversion mode (IM) is compared, it is concluded from the simulations that JL vertically stacked MoS2 NSFET more immune to short channel effects such as threshold voltage (Vth) roll-off and drain induced barrier lowering (DIBL).