Researcher profile

Sander G. Huisman

Sander G. Huisman contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
1topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

14 published item(s)

preprint2026arXiv

Settling of chiral particles in a turbulent flow

Chiral particles are experimentally investigated while settling inwater with various turbulence intensity levels. The locations and orientations of the particles are tracked over time, allowing the close investigation of the particles' settling dynamics. The generated turbulent flow is measured using laser Doppler anemometry (LDA), and the turbulence strength varies between experiments in the range $0 \leq Re_λ\leq 250$. Starting with quiescent particle settling, the chiral particle's orientation dynamics are studied, revealing a preferred alignment and a strong translation-rotation coupling. The particle chirality determines the preferred rotation direction, though the alignment and translation-rotation coupling gradually vanish with increasing turbulence. We identify multiple settling modes for the chiral particles, which are characterised by the evolution of the rotation angles. Finally, a theoretical model assuming a simplified chiral particle in Stokes flow clarifies the emergence of each settling mode.

preprint2026arXiv

Transition from classical to ultimate melting

Melting is omnipresent in nature and technology, with applications ranging from metallurgy, biology, food science, and latent thermal energy storage to oceanography, geophysics, and climate science, and occurring on all scales from sub-millimeter to global scales. The key objective is to understand the rate at which an object melts as a function of its size and of the ambient conditions. To achieve this it is important to be able to extrapolate from small scale experiments and observations to large or even global scales. This is done by scaling laws. However, these are only meaningful if there is no transition from one scaling relation to another one. Here we show, however, that for both fixed and freely-advected melting objects immersed in a turbulent flow a melting transition does exist, namely from slow melting at the small scales to fast melting at the large scales. We do so by controlled melting experiments and corresponding direct numerical simulations, covering four orders of magnitude in scale. The transition corresponds to the transition from a laminar-type boundary layer around the melting object to a turbulent-type boundary layer, i.e., from so-called classical turbulence to ultimate turbulence, with its enhanced transport properties. Our results thus provide a quantitative understanding of the flow physics of the melting process and thereby enable a better extrapolation and prediction of melt rates on large scales such as relevant in geophysics, oceanography, and climate science.

preprint2025arXiv

Collective effects of neighbouring melting ice objects

We present a study on the melting dynamics of neighbouring ice bodies by means of idealised simulations, focusing on collective effects, with the goal of obtaining fundamental insight into how collective interactions influence the melting of ice. Two neighbouring (vertically or horizontally aligned), square-shaped, and equally sized ice objects (size on the order of centimetres) are immersed in quiescent fresh water at a temperature of 20°C. By performing two-dimensional direct numerical simulations, and using the phase-field method to model the phase change, the collective melting of these objects is studied. When the objects are horizontally aligned, no significant influence of the neighbouring object on the melting time is observed. On the other hand, when vertically aligned, though the melting of the upper object is mostly unaffected, the melting time and the morphology of the lower ice body strongly depends on the initial inter-object distance. We report that the melting of the bottom object can be enhanced by more than 10%, or delayed more than 20%, displaying a non-monotonic dependence on the initial object size. We show that this behaviour results from a non-trivial competition between layering of cold fluid, which lowers the heat transfer, and convective flows, which favour mixing and heat transfer. For this melting in mixed convection, we were able to collapse our data onto a single curve.

preprint2022arXiv

Micro-droplet nucleation through solvent exchange in a turbulent buoyant jet

Solvent exchange is a process involving mixing between a good solvent with dissolved solute and a poor solvent. The process creates local oversaturation which causes the nucleation of minute solute droplets. Such ternary systems on a macro-scale have remained unexplored in the turbulent regime. We experimentally study the solvent exchange process by injecting mixtures of ethanol and trans-anethole into water, forming a turbulent buoyant jet in upward direction. Locally, turbulent mixing causes oversaturation of the trans-anethole following turbulent entrainment. We optically measure the concentration of the nucleated droplets using a light attenuation technique and find that the radial concentration profile is sub-Gaussian. In contrast to the entrainment-based models, the spatial evolution of the oversaturation reveals continuous droplet nucleation downstream and radially across the jet, which we attribute to the limited mixing capacity of the jet. Though we are far from a full quantitative understanding, this work extends the knowledge on solvent exchange into the turbulent regime, and brings in a novel type of flow, broadening the scope of multicomponent, multiphase turbulent jets with phase transition.

preprint2022arXiv

Recent developments of turbulent emulsions in Taylor-Couette flow

Emulsions are common in many natural and industrial settings. Recently, much attention has been put on understanding the dynamics of turbulent emulsions. This paper reviews some recent studies of emulsions in turbulent Taylor-Couette flow, mainly focusing on the statistics of the dispersed phase and the global momentum transport of the system. We first study the size distribution and the breakup mechanism of the dispersed droplets for turbulent emulsions with a low volume-fraction (dilute) of the dispersed phase. For systems with a high volume-fraction of the dispersed phase (dense), we address the detailed response of the global transport (effective viscosity) of the turbulent emulsion and its connection to the droplet statistics. Finally, we will discuss catastrophic phase inversions, which can happen when the volume fraction of the dispersed phase exceeds a critical value during dynamic emulsification. We end the manuscript with a summary and an outlook including some open questions for future research. This article is part of the theme issue `Taylor-Couette and Related Flows on the Centennial of Taylor's Seminal Philosophical Transactions Paper'.

preprint2022arXiv

Supergravitational turbulent thermal convection

High-Rayleigh number convective turbulence is ubiquitous in many natural phenomena and in industries, such as atmospheric circulations, oceanic flows, flows in the fluid core of planets, and energy generations. In this work, we present a novel approach to boost the Rayleigh number in thermal convection by exploiting centrifugal acceleration and rapidly rotating a cylindrical annulus to reach an effective gravity of 60 times Earth's gravity. We show that in the regime where the Coriolis effect is strong, the scaling exponent of Nusselt number versus Rayleigh number exceeds one-third once the Rayleigh number is large enough. The convective rolls revolve in prograde direction, signifying the emergence of zonal flow. The present findings open a new avenue on the exploration of high-Rayleigh number turbulent thermal convection and will improve the understanding of the flow dynamics and heat transfer processes in geophysical and astrophysical flows and other strongly rotating systems.

preprint2021arXiv

Broken mirror symmetry of tracer's trajectories in turbulence

Topological properties of physical systems play a crucial role in our understanding of nature, yet their experimental determination remains elusive. We show that the mean helicity, a dynamical invariant in ideal flows, quantitatively affects trajectories of fluid elements: the linking number of Lagrangian trajectories depends on the mean helicity. Thus, a global topological invariant and a topological number of fluid trajectories become related, and we provide an empirical expression linking them. The relation shows the existence of long-term memory in the trajectories: the links can be made of the trajectory up to a given time, with particles positions in the past. This property also allows experimental measurements of mean helicity.

preprint2020arXiv

Bubbly drag reduction using a hydrophobic inner cylinder in Taylor-Couette turbulence

In this study we experimentally investigate bubbly drag reduction in a highly turbulent flow of water with dispersed air at $5.0 \times 10^{5} \leq \text{Re} \leq 1.7 \times 10^{6}$ over a non-wetting surface containing micro-scale roughness. To do so, the Taylor-Couette geometry is used, allowing for both accurate global drag and local flow measurements. The inner cylinder - coated with a rough, hydrophobic material - is rotating, whereas the smooth outer cylinder is kept stationary. The crucial control parameter is the air volume fraction $α$ present in the working fluid. For small volume fractions ($α< {4}\,\%$), we observe that the surface roughness from the coating increases the drag. For large volume fractions of air ($α\geq 4\,\%$), the drag decreases compared to the case with both the inner and outer cylinders uncoated, i.e. smooth and hydrophilic, using the same volume fraction of air. This suggests that two competing mechanisms are at place: on the one hand the roughness invokes an extension of the log-layer - resulting in an increase in drag - and on the other hand there is a drag-reducing mechanism of the hydrophobic surface interacting with the bubbly liquid. The balance between these two effects determines whether there is overall drag reduction or drag enhancement. For further increased bubble concentration $α= {6}\,\%$ we find a saturation of the drag reduction effect. Our study gives guidelines for industrial applications of bubbly drag reduction in hydrophobic wall-bounded turbulent flows.

preprint2020arXiv

Catastrophic phase inversion in high-Reynolds number turbulent Taylor--Couette flow

Emulsions are omnipresent in the food industry, health care, and chemical synthesis. In this Letter the dynamics of meta-stable oil-water emulsions in highly turbulent ($10^{11}\leq\text{Ta}\leq 3\times 10^{13}$) Taylor--Couette flow, far from equilibrium, is investigated. By varying the oil-in-water void fraction, catastrophic phase inversion between oil-in-water and water-in-oil emulsions can be triggered, changing the morphology, including droplet sizes, and rheological properties of the mixture, dramatically. The manifestation of these different states is exemplified by combining global torque measurements and local in-situ laser induced fluorescence (LIF) microscopy imaging. Despite the turbulent state of the flow and the dynamic equilibrium of the oil-water mixture, the global torque response of the system is found to be as if the fluid were Newtonian, and the effective viscosity of the mixture was found to be several times bigger or smaller than either of its constituents.

preprint2020arXiv

Double maxima of angular momentum transport in $η=0.91$ TC turbulence

We use experiments and direct numerical simulations to probe the phase-space of low-curvature Taylor--Couette (TC) flow in the vicinity of the ultimate regime. The cylinder radius ratio is fixed at $η=r_i/r_o=0.91$. Non-dimensional shear drivings (Taylor numbers $\text{Ta}$) in the range $10^7\leq\text{Ta}\leq10^{11}$ are explored for both co- and counter-rotating configurations. In the $\text{Ta}$ range $10^8\leq\text{Ta}\leq10^{10}$, we observe two local maxima of the angular momentum transport as a function of the cylinder rotation ratio, which can be described as either as &#34;co-&#34; and &#34;counter-rotating&#34; due to their location or as &#34;broad&#34; and &#34;narrow&#34; due to their shape. We confirm that the broad peak is accompanied by the strengthening of the large-scale structures, and that the narrow peak appears once the driving (Ta) is strong enough. As first evidenced in numerical simulations by Brauckmann \emph{et al.}~(2016), the broad peak is produced by centrifugal instabilities and that the narrow peak is a consequence of shear instabilities. We describe how the peaks change with $\text{Ta}$ as the flow becomes more turbulent. Close to the transition to the ultimate regime when the boundary layers (BLs) become turbulent, the usual structure of counter-rotating Taylor vortex pairs breaks down and stable unpaired rolls appear locally. We attribute this state to changes in the underlying roll characteristics during the transition to the ultimate regime. Further changes in the flow structure around $\text{Ta}\approx10^{10}$ cause the broad peak to disappear completely and the narrow peak to move. This second transition is caused when the regions inside the BLs which are locally smooth regions disappear and the whole boundary layer becomes active.

preprint2020arXiv

Effect of axially varying sandpaper roughness on bubbly drag reduction in Taylor-Couette turbulence

We experimentally investigate the influence of alternating rough and smooth walls on bubbly drag reduction (DR). We apply rough sandpaper bands of width $s$ between $48.4\,mm$ and $148.5\,mm$, and roughness height $k = 695\,μm$, around the smooth inner cylinder (IC) of the Twente Turbulent Taylor-Couette facility. Between sandpaper bands, the IC is left uncovered over similar width $s$, resulting in alternating rough and smooth bands, a constant pattern in axial direction. We measure the DR in water that originates from introducing air bubbles to the fluid at (shear) Reynolds numbers $\textit{Re}_s$ ranging from $0.5 \times 10^6$ to $1.8 \times 10^6$. Results are compared to bubbly DR measurements with a completely smooth IC and an IC that is completely covered with sandpaper of the same roughness $k$. The outer cylinder is left smooth for all variations. Results are also compared to bubbly DR measurements where a smooth outer cylinder is rotating in opposite direction to the smooth IC. This counter rotation induces secondary flow structures that are very similar to those observed when the IC is composed of alternating rough and smooth bands. For the measurements with roughness, the bubbly DR is found to initially increase more strongly with $\textit{Re}_s$, before levelling off to reach a value that no longer depends on $\textit{Re}_s$. This is attributed to a more even axial distribution of the air bubbles, resulting from the increased turbulence intensity of the flow compared to flow over a completely smooth wall at the same $\textit{Re}_s$. The air bubbles are seen to accumulate at the rough wall sections in the flow. Here, locally, the drag is largest and so the drag reducing effect of the bubbles is felt strongest. Therefore, a larger maximum value of bubbly DR is found for the alternating rough and smooth walls compared to the completely rough wall.

preprint2020arXiv

Heat transfer enhancement in Rayleigh-Bénard convection using a single passive barrier

In this numerical study on Rayleigh-Bénard convection we seek to improve the heat transfer by passive means. To this end we introduce a single tilted conductive barrier centered in an aspect ratio one cell, breaking the symmetry of the geometry and to channel the ascending hot and descending cold plumes. We study the global and local heat transfer and the flow organization for Rayleigh numbers $10^5 \leq Ra \leq 10^9$ for a fixed Prandtl number of $Pr=4.3$. We find that the global heat transfer can be enhanced up to $18\%$, and locally around $800\%$. The averaged Reynolds number is always decreased when a barrier is introduced, even for those cases where the global heat transfer is increased. We map the entire parameter space spanned by the orientation and the size of a single barrier for $Ra=10^8$.

preprint2019arXiv

Controlling secondary flow in Taylor-Couette turbulence through spanwise-varying roughness

Highly turbulent Taylor-Couette flow with spanwise-varying roughness is investigated experimentally and numerically (direct numerical simulations (DNS) with an immersed boundary method (IBM)) to determine the effects of the spacing and axial width $s$ of the spanwise varying roughness on the total drag and {on} the flow structures. We apply sandgrain roughness, in the form of alternating {rough and smooth} bands to the inner cylinder. Numerically, the Taylor number is $\mathcal{O}(10^9)$ and the roughness width is varied between $0.47\leq \tilde{s}=s/d \leq 1.23$, where $d$ is the gap width. Experimentally, we explore $\text{Ta}=\mathcal{O}(10^{12})$ and $0.61\leq \tilde s \leq 3.74$. For both approaches the radius ratio is fixed at $η=r_i/r_o = 0.716$, with $r_i$ and $r_o$ the radius of the inner and outer cylinder respectively. We present how the global transport properties and the local flow structures depend on the boundary conditions set by the roughness spacing $\tilde{s}$. Both numerically and experimentally, we find a maximum in the angular momentum transport as function of $\tilde s$. This can be atributed to the re-arrangement of the large-scale structures triggered by the presence of the rough stripes, leading to correspondingly large-scale turbulent vortices.

preprint2019arXiv

Drag reduction in boiling Taylor-Couette turbulence

We create a highly controlled lab environment-accessible to both global and local monitoring-to analyse turbulent boiling flows and in particular their shear stress in a statistically stationary state. Namely, by precisely monitoring the drag of strongly turbulent Taylor-Couette flow (the flow in between two co-axially rotating cylinders, Reynolds number $\textrm{Re}\approx 10^6$) during its transition from non-boiling to boiling, we show that the intuitive expectation, namely that a few volume percent of vapor bubbles would correspondingly change the global drag by a few percent, is wrong. Rather, we find that for these conditions a dramatic global drag reduction of up to 45% occurs. We connect this global result to our local observations, showing that for major drag reduction the vapor bubble deformability is crucial, corresponding to Weber numbers larger than one. We compare our findings with those for turbulent flows with gas bubbles, which obey very different physics than vapor bubbles. Nonetheless, we find remarkable similarities and explain these.