Researcher profile

Stephan Gekle

Stephan Gekle contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
8topics
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

5 published item(s)

preprint2022arXiv

On the accuracy and performance of the lattice Boltzmann method with 64-bit, 32-bit and novel 16-bit number formats

Fluid dynamics simulations with the lattice Boltzmann method (LBM) are very memory-intensive. Alongside reduction in memory footprint, significant performance benefits can be achieved by using FP32 (single) precision compared to FP64 (double) precision, especially on GPUs. Here, we evaluate the possibility to use even FP16 and Posit16 (half) precision for storing fluid populations, while still carrying arithmetic operations in FP32. For this, we first show that the commonly occurring number range in the LBM is a lot smaller than the FP16 number range. Based on this observation, we develop novel 16-bit formats - based on a modified IEEE-754 and on a modified Posit standard - that are specifically tailored to the needs of the LBM. We then carry out an in-depth characterization of LBM accuracy for six different test systems with increasing complexity: Poiseuille flow, Taylor-Green vortices, Karman vortex streets, lid-driven cavity, a microcapsule in shear flow (utilizing the immersed-boundary method) and finally the impact of a raindrop (based on a Volume-of-Fluid approach). We find that the difference in accuracy between FP64 and FP32 is negligible in almost all cases, and that for a large number of cases even 16-bit is sufficient. Finally, we provide a detailed performance analysis of all precision levels on a large number of hardware microarchitectures and show that significant speedup is achieved with mixed FP32/16-bit.

preprint2021arXiv

Analytic Solution to the Piecewise Linear Interface Construction Problem and its Application in Curvature Calculation for Volume-of-Fluid Simulation Codes

The plane-cube intersection problem has been around in literature since 1984 and iterative solutions to it have been used as part of piecewise linear interface construction (PLIC) in computational fluid dynamics simulation codes ever since. In many cases, PLIC is the bottleneck of these simulations regarding compute time, so a faster, analytic solution to the plane-cube intersection would greatly reduce compute time for such simulations. We derive an analytic solution for all intersection cases and compare it to the one previous solution from Scardovelli and Zaleski (Ruben Scardovelli and Stephane Zaleski. "Analytical relations connecting linear interfaces and volume fractions in rectangular grids". In: Journal of Computational Physics 164.1 (2000), pp. 228-237.), which we further improve to include edge cases and micro-optimize to reduce arithmetic operations and branching. We then extend our comparison regarding compute time and accuracy to include two different iterative solutions as well. We find that the best choice depends on the employed hardware platform: on the CPU, Newton-Raphson is fastest with vectorization while analytic solutions perform better without. The reason for this is that vectorization instruction sets do not include trigonometric functions as used in the analytic solutions. On the GPU, the fastest method is our optimized version of the analytic SZ solution. We finally provide details on one of the applications of PLIC: curvature calculation for the Volume-of-Fluid model used for free surface fluid simulations in combination with the lattice Boltzmann method.

preprint2021arXiv

Vortical flow structures induced by red blood cells in capillaries

Knowledge about the flow field of the plasma around the red blood cells in capillary flow is important for a physical understanding of blood flow and the transport of micro- and nanoparticles and molecules in the flowing plasma. We conduct an experimental study on the flow field around red blood cells in capillary flow that are complemented by simulations of vortical flow between red blood cells. Red blood cells were injected in a 10x12 micrometer rectangular microchannel at a low hematocrit and the flow field around a single or two cells have been characterized thanks to a highspeed camera and by tracking 250 nm nanoparticles in flow behaving as tracers. While the flow field around a steady croissant shape is found to be relatively similar to that of a rigid sphere, the flow field around a slipper shape exhibits a small vortex at the rear of the red blood cell. Even more pronounced are vortex-like structures observed in the central region between two neighboring croissants. Conclusions: The rotation frequency of the vortices is to a good approximation, inversely proportional to the distance between the cells. Our experimental data are confirmed and complemented by numerical simulations.

preprint2010arXiv

Supersonic Air Flow due to Solid-Liquid Impact

A solid object impacting on liquid creates a liquid jet due to the collapse of the impact cavity. Using visualization experiments with smoke particles and multiscale simulations we show that in addition a high-speed air-jet is pushed out of the cavity. Despite an impact velocity of only 1 m/s, this air-jet attains \emph{supersonic} speeds already when the cavity is slightly larger than 1 mm in diameter. The structure of the air flow resembles closely that of compressible flow through a nozzle -- with the key difference that here the "nozzle" is a \emph{liquid} cavity shrinking rapidly in time.

preprint2009arXiv

Nucleation threshold and deactivation mechanisms of nanoscopic cavitation nuclei

The acoustic nucleation threshold for bubbles trapped in cavities has theoretically been predicted within the crevice theory by Atchley & Prosperetti [J. Acoust. Soc. Am. 86, 1065-1084 (1989)]. Here, we determine this threshold experimentally, by applying a single pressure pulse to bubbles trapped in cylindrical nanoscopic pits ("artificial crevices") with radii down to 50 nm. By decreasing the minimum pressure stepwise, we observe the threshold for which the bubbles start to nucleate. The experimental results are quantitatively in excellent agreement with the theoretical predictions of Atchley & Prosperetti. In addition, we provide the mechanism which explains the deactivation of cavitation nuclei: gas diffusion together with an aspherical bubble collapse. Finally, we present superhydrophobic nuclei which cannot be deactivated, unless with a high-speed liquid jet directed into the pit.