BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Stephan Gekle

Stephan Gekle appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dyncond-mat.softphysics.comp-phBiological PhysicsCell BehaviorComputational Geometrycond-mat.stat-mechDistributed, Parallel, and Cluster Computingnlin.CDphysics.ao-phphysics.atom-ph

Catalog footprint

What is connected

16works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

16 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hydrodynamic interaction between particles near elastic interfaces

We present an analytical calculation of the hydrodynamic interaction between two spherical particles near an elastic interface such as a cell membrane. The theory predicts the frequency dependent self- and pair-mobilities accounting for the finite particle size up to the 5th order in the ratio between particle diameter and wall distance as well as between diameter and interparticle distance. We find that particle motion towards a membrane with pure bending resistance always leads to mutual repulsion similar as in the well-known case of a hard-wall. In the vicinity of a membrane with shearing resistance, however, we observe an attractive interaction in a certain parameter range which is in contrast to the behavior near a hard wall. This attraction might facilitate surface chemical reactions. Furthermore, we show that there exists a frequency range in which the pair-mobility for perpendicular motion exceeds its bulk value, leading to short-lived superdiffusive behavior. Using the analytical particle mobilities we compute collective and relative diffusion coefficients. The appropriateness of the approximations in our analytical results is demonstrated by corresponding boundary integral simulations which are in excellent agreement with the theoretical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Long-lived anomalous thermal diffusion induced by elastic cell membranes on nearby particles

The physical approach of a small particle (virus, medical drug) to the cell membrane represents the crucial first step before active internalization and is governed by thermal diffusion. Using a fully analytical theory we show that the stretching and bending of the elastic membrane by the approaching particle induces a memory in the system which leads to anomalous diffusion, even though the particle is immersed in a purely Newtonian liquid. For typical cell membranes the transient subdiffusive regime extends beyond 10 ms and can enhance residence times and possibly binding rates up to 50\%. Our analytical predictions are validated by numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Mobility of an axisymmetric particle near an elastic interface

Using a fully analytical theory, we compute the leading order corrections to the translational, rotational and translation-rotation coupling mobilities of an arbitrary axisymmetric particle immersed in a Newtonian fluid moving near an elastic cell membrane that exhibits resistance towards stretching and bending. The frequency-dependent mobility corrections are expressed as general relations involving separately the particle's shape-dependent bulk mobility and the shape-independent parameters such as the membrane-particle distance, the particle orientation and the characteristic frequencies associated with shearing and bending of the membrane. This makes the equations applicable to an arbitrary-shaped axisymmetric particle provided that its bulk mobilities are known, either analytically or numerically. For a spheroidal particle, these general relations reduce to simple expressions in terms of the particle's eccentricity. We find that the corrections to the translation-rotation coupling mobility are primarily determined by bending, whereas shearing manifests itself in a more pronounced way in the rotational mobility. We demonstrate the validity of the analytical approximations by a detailed comparison with boundary integral simulations of a truly extended spheroidal particle. They are found to be in a good agreement over the whole range of applied frequencies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Particle mobility between two planar elastic membranes: Brownian motion and membrane deformation

We study the motion of a solid particle immersed in a Newtonian fluid and confined between two parallel elastic membranes possessing shear and bending rigidity. The hydrodynamic mobility depends on the frequency of the particle motion due to the elastic energy stored in the membrane. Unlike the single-membrane case, a coupling between shearing and bending exists. The commonly used approximation of superposing two single-membrane contributions is found to give reasonable results only for motions in the parallel, but not in the perpendicular direction. We also compute analytically the membrane deformation resulting from the motion of the particle, showing that the presence of the second membrane reduces deformation. Using the fluctuation-dissipation theorem we compute the Brownian motion of the particle, finding a long-lasting subdiffusive regime at intermediate time scales. We finally assess the accuracy of the employed point-particle approximation via boundary-integral simulations for a truly extended particle. They are found to be in excellent agreement with the analytical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Air flow in a collapsing cavity

We experimentally study the airflow in a collapsing cavity created by the impact of a circular disk on a water surface. We measure the air velocity in the collapsing neck in two ways: Directly, by means of employing particle image velocimetry of smoke injected into the cavity and indirectly, by determining the time rate of change of the volume of the cavity at pinch-off and deducing the air flow in the neck under the assumption that the air is incompressible. We compare our experiments to boundary integral simulations and show that close to the moment of pinch-off, compressibility of the air starts to play a crucial role in the behavior of the cavity. Finally, we measure how the air flow rate at pinch-off depends on the Froude number and explain the observed dependence using a theoretical model of the cavity collapse.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Collapse of a non-axisymmetric, impact-created air cavity in water

The axisymmetric collapse of a cylindrical air cavity in water follows a universal power law with logarithmic corrections. Nonetheless, it has been suggested that the introduction of a small azimuthal disturbance induces a long term memory effect, reflecting in oscillations which are no longer universal but remember the initial condition. In this work, we create non-axisymmetric air cavities by driving a metal disc through an initially-quiescent water surface and observe their subsequent gravity-induced collapse. The cavities are characterized by azimuthal harmonic disturbances with a single mode number $m$ and amplitude $a_m$. For small initial distortion amplitude (1 or 2% of the mean disc radius), the cavity walls oscillate linearly during collapse, with nearly constant amplitude and increasing frequency. As the amplitude is increased, higher harmonics are triggered in the oscillations and we observe more complex pinch-off modes. For small amplitude disturbances we compare our experimental results with the model for the amplitude of the oscillations by Schmidt et al. (2009) and the model for the collapse of an axisymmetric impact-created cavity previously proposed by Bergmann et al. (2009b). By combining these two models we can reconstruct the three-dimensional shape of the cavity at any time before pinch-off.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Breakup of diminutive Rayleigh jets

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Collapse of Non-Axisymmetric Cavities

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Generation and Breakup of Worthington Jets After Cavity Collapse

Helped by the careful analysis of their experimental data, Worthington (1897) described roughly the mechanism underlying the formation of high-speed jets ejected after the impact of an axisymmetric solid on a liquid-air interface. In this work we combine detailed boundary-integral simulations with analytical modeling to describe the formation and break-up of such Worthington jets in two common physical systems: the impact of a circular disc on a liquid surface and the release of air bubbles from an underwater nozzle. We first show that the jet base dynamics can be predicted for both systems using our earlier model in Gekle, Gordillo, van der Meer and Lohse. Phys. Rev. Lett. 102 (2009). Nevertheless, our main point here is to present a model which allows us to accurately predict the shape of the entire jet. Good agreement with numerics and some experimental data is found. Moreover, we find that, contrarily to the capillary breakup of liquid cylinders in vacuum studied by Rayleigh, the breakup of stretched liquid jets at high values of both Weber and Reynolds numbers is not triggered by the growth of perturbations coming from an external source of noise. Instead, the jet breaks up due to the capillary deceleration of the liquid at the tip which produces a corrugation to the jet shape. This perturbation, which is self-induced by the flow, will grow in time promoted by a capillary mechanism. We are able to predict the exact shape evolution of Worthington jets ejected after the impact of a solid object - including the size of small droplets ejected from the tip due to a surface-tension driven instability - using as the single input parameters the minimum radius of the cavity and the flow field before the jet emerges.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Controlled impact of a disk on a water surface: Cavity dynamics

In this paper we study the transient surface cavity which is created by the controlled impact of a disk of radius h0 on a water surface at Froude numbers below 200. The dynamics of the transient free surface is recorded by high speed imaging and compared to boundary integral simulations. An excellent agreement is found between both. The flow surrounding the cavity is measured with high speed particle image velocimetry and is found to also agree perfectly with the flow field obtained from the simulations. We present a simple model for the radial dynamics of the cavity based on the collapse of an infinite cylinder. This model accounts for the observed asymmetry of the radial dynamics between the expansion and contraction phase of the cavity. It reproduces the scaling of the closure depth and total depth of the cavity which are both found to scale roughly proportional to Fr^{1/2} with a weakly Froude number dependent prefactor. In addition, the model accurately captures the dynamics of the minimal radius of the cavity, the scaling of the volume Vbubble of air entrained by the process, namely Vbubble/h0^3 proportional (1 + 0.26Fr^{1/2})Fr^{1/2}, and gives insight into the axial asymmetry of the pinch-off process.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Hydrogen atom in crossed electric and magnetic fields: Phase space topology and torus quantization via periodic orbits

A hierarchical ordering is demonstrated for the periodic orbits in a strongly coupled multidimensional Hamiltonian system, namely the hydrogen atom in crossed electric and magnetic fields. It mirrors the hierarchy of broken resonant tori and thereby allows one to characterize the periodic orbits by a set of winding numbers. With this knowledge, we construct the action variables as functions of the frequency ratios and carry out a semiclassical torus quantization. The semiclassical energy levels thus obtained agree well with exact quantum calculations.

0 citations0 reviewsNo code linkedOpen access