Source author record

Thorsten Poeschel

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

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

24 published item(s)

preprint2026arXiv

Dilatancy-induced surface deformation in dense cohesive granular media

When granular materials with interstitial liquid bridges are sheared in a split-bottom cell, a localized shear band develops, accompanied by a surface elevation. Cohesion, governed by the surface tension of the interstitial liquid, enhances dilatancy in dense cohesive packings, leading to expansion within the shear band and the emergence of a surface elevation. Surface deformation is observed not only in cohesive systems with high particle density and large liquid surface tension, but also in those with lower values of these parameters. The equivalent Bond number arises as a key control parameter for the surface deformation, shaping both the evolution of the surface profile and the packing density. At higher shear rates, inertial effects dominate dilatancy, resulting in less pronounced surface deformation.

preprint2016arXiv

Stable algorithm for event detection in event-driven particle dynamics: logical states

Following the recent development of a stable event-detection algorithm for hard-sphere systems, the implications of more complex interaction models are examined. The relative location of particles leads to ambiguity when it is used to determine the interaction state of a particle in stepped potentials, such as the square-well model. To correctly predict the next event in these systems, the concept of an additional state that is tracked separately from the particle position is introduced and integrated into the stable algorithm for event detection.

preprint2014arXiv

Stable algorithm for event detection in event-driven particle dynamics

Event-Driven Particle Dynamics is a fast and precise method to simulate particulate systems of all scales. In this work it is demonstrated that, despite the high accuracy of the method, the finite machine precision leads to simulations entering invalid states where the dynamics are undefined. A general event-detection algorithm is proposed which handles these situations in a stable and efficient manner. This requires a definition of the dynamics of invalid states and leads to improved algorithms for event-detection in hard-sphere systems.

preprint2013arXiv

A numerical study of the Navier-Stokes transport coefficients for 2D granular hydrodynamics

A numerical study is presented to analyze the thermal mechanisms of unsteady, supersonic granular flow, by means of hydrodynamic simulations of the Navier-Stokes granular equations. For this purpose a paradigmatic problem in granular dynamics such as the Faraday instability is selected. Two different approaches for the Navier-Stokes transport coefficients for granular materials are considered, namely the traditional Jenkins-Richman theory for moderately dense quasi-elastic grains, and the improved Garzó-Dufty-Lutsko theory for arbitrary inelasticity, which we also present here. Both solutions are compared with event-driven simulations of the same system under the same conditions, by analyzing the density, the temperature and the velocity field. Important differences are found between the two approaches leading to interesting implications. In particular, the heat transfer mechanism coupled to the density gradient which is a distinctive feature of inelastic granular gases, is responsible for a major discrepancy in the temperature field and hence in the diffusion mechanisms.

preprint2013arXiv

Energy Dissipation in Driven Granular Matter in the Absence of Gravity

We experimentally investigate the energy dissipation rate in sinusoidally driven boxes which are partly filled by granular material under conditions of weightlessness. We identify two different modes of granular dynamics, depending on the amplitude of driving, $A$. For intense forcing, A>A_0, the material is found in the collect-and-collide regime where the center of mass of the granulate moves synchronously with the driven container while for weak forcing, A<A_0, the granular material exhibits gas-like behavior. Both regimes correspond to different dissipation mechanisms, leading to different scaling with amplitude and frequency of the excitation and with the mass of the granulate. For the collect-and-collide regime, we explain the dependence on frequency and amplitude of the excitation by means of an effective one-particle model. For both regimes, the results may be collapsed to a single curve characterizing the physics of granular dampers.

preprint2013arXiv

Event-driven Molecular Dynamics of Soft Particles

The dynamics of dissipative soft-sphere gases obeys Newton's equation of motion which are commonly solved numerically by (force-based) Molecular Dynamics schemes. With the assumption of instantaneous, pairwise collisions, the simulation can be accelerated considerably using event-driven Molecular Dynamics, where the coefficient of restitution is derived from the interaction force between particles. Recently it was shown, however, that this approach may fail dramatically, that is, the obtained trajectories deviate significantly from the ones predicted by Newton's equations. In this paper, we generalize the concept of the coefficient of restitution and derive a numerical scheme which, in the case of dilute systems and frictionless interaction, allows us to perform highly efficient event-driven Molecular Dynamics simulations even for non-instantaneous collisions. We show that the particle trajectories predicted by the new scheme agree perfectly with the corresponding (force-based) Molecular Dynamics, except for a short transient period whose duration corresponds to the duration of the contact. Thus, the new algorithm solves Newton's equations of motion like force-based MD while preserving the advantages of event-driven simulations.

preprint2013arXiv

Numerical modeling of the wind flow over a transverse dune

Transverse dunes, which form under unidirectional winds and have fixed profile in the direction perpendicular to the wind, occur on all celestial objects of our solar system where dunes have been detected. Here we perform a numerical study of the average turbulent wind flow over a transverse dune by means of computational fluid dynamics simulations. We find that the length of the zone of recirculating flow at the dune lee --- the {\em{separation bubble}} --- displays a surprisingly strong dependence on the wind shear velocity, $u_{\ast}$: it is nearly independent of $u_{\ast}$ for shear velocities within the range between $0.2\,$m$$s and $0.8\,$m$$s but increases linearly with $u_{\ast}$ for larger shear velocities. Our calculations show that transport in the direction opposite to dune migration within the separation bubble can be sustained if $u_{\ast}$ is larger than approximately $0.39\,$m$$s, whereas a larger value of $u_{\ast}$ (about $0.49\,$m$$s) is required to initiate this reverse transport.

preprint2013arXiv

Origins of barchan dune asymmetry: insights from numerical simulations

Barchan dunes --- crescent-shaped dunes that form in areas of unidirectional winds and low sand availability --- commonly display an asymmetric shape, with one limb extended downwind. Several factors have been identified as potential causes for barchan dune asymmetry on Earth and Mars: asymmetric bimodal wind regime, topography, influx asymmetry and dune collision. However, the dynamics and potential range of barchan morphologies emerging under each specific scenario that leads to dune asymmetry are far from being understood. In the present work, we use dune modeling in order to investigate the formation and evolution of asymmetric barchans. We find that a bimodal wind regime causes limb extension when the divergence angle between primary and secondary winds is larger than $90^{\circ}$, whereas the extended limb evolves into a seif dune if the ratio between secondary and primary transport rates is larger than 25%. Calculations of dune formation on an inclined surface under constant wind direction also lead to barchan asymmetry, however no seif dune is obtained from surface tilting alone. Asymmetric barchans migrating along a tilted surface move laterally, with transverse migration velocity proportional to the slope of the terrain. Limb elongation induced by topography can occur when a barchan crosses a topographic rise. Furthermore, transient asymmetric barchan shapes with extended limb also emerge during collisions between dunes or due to an asymmetric influx. Our findings can be useful for making quantitative inference on local wind regimes or spatial heterogeneities in transport conditions of planetary dune fields hosting asymmetric barchans.

preprint2011arXiv

Coefficient of Restitution as a Fluctuating Quantity

The coefficient of restitution of a spherical particle in contact with a flat plate is investigated as a function of the impact velocity. As an experimental observation we notice non-trivial (non-Gaussian) fluctuations of the measured values. For a fixed impact velocity, the probability density of the coefficient of restitution, $p(ε)$, is formed by two exponential functions (one increasing, one decreasing) of different slope. This behavior may be explained by a certain roughness of the particle which leads to energy transfer between the linear and rotational degrees of freedom.

preprint2011arXiv

Collision of Viscoelastic Spheres: Compact Expressions for the Coefficient of Normal Restitution

The coefficient of restitution of colliding viscoelastic spheres is analytically known as a complete series expansion in terms of the impact velocity where all (infinitely many) coefficients are known. While beeing analytically exact, this result is not suitable for applications in efficient event-driven Molecular Dynamics (eMD) or Monte Carlo (MC) simulations. Based on the analytic result, here we derive expressions for the coefficient of restitution which allow for an application in efficient eMD and MC simulations of granular Systems.

preprint2011arXiv

Oblique Impact of Frictionless Spheres: On the Limitations of Hard Sphere Models for Granular Dynamics

When granular systems are modeled by frictionless hard spheres, particle-particle collisions are considered as instantaneous events. This implies that while the velocities change according to the collision rule, the positions of the particles are the same before and after such an event. We show that depending on the material and system parameters, this assumption may fail. For the case of viscoelastic particles we present a universal condition which allows to assess whether the hard-sphere modeling and, thus, event-driven Molecular Dynamics simulations are justified.

preprint2010arXiv

Movers and shakers: Granular damping in microgravity

The response of an oscillating granular damper to an initial perturbation is studied using experiments performed in microgravity and granular dynamics mulations. High-speed video and image processing techniques are used to extract experimental data. An inelastic hard sphere model is developed to perform simulations and the results are in excellent agreement with the experiments. The granular damper behaves like a frictional damper and a linear decay of the amplitude is bserved. This is true even for the simulation model, where friction forces are absent. A simple expression is developed which predicts the optimal damping conditions for a given amplitude and is independent of the oscillation frequency and particle inelasticities.

preprint2010arXiv

Three Balls Problem Revisited - On the Limitations of Event-Driven Modeling

If a tennis ball is held above a basket ball with their centers vertically aligned, and the balls are released to collide with the floor, the tennis ball may rebound at a surprisingly high speed. We show in this article that the simple textbook explanation of this effect is an oversimplification, even for the limit of perfectly elastic particles. Instead, there may occur a rather complex scenario including multiple collisions which may lead to a very different final velocity as compared with the velocity resulting from the oversimplified model.

preprint2006arXiv

Translations and Rotations are correlated in Granular Gases

In a granular gas of rough particles the axis of rotation is shown to be correlated with the translational velocity of the particles. The average relative orientation of angular and linear velocities depends on the parameters which characterise the dissipative nature of the collision. We derive a simple theory for these correlations and validate it with numerical simulations for a wide range of coefficients of normal and tangential restitution. The limit of smooth spheres is shown to be singular: even an arbitrarily small roughness of the particles gives rise to orientational correlations.

preprint2003arXiv

Giant fluctuations at a granular phase separation threshold

We investigate a phase separation instability that occurs in a system of nearly elastically colliding hard spheres driven by a thermal wall. If the aspect ratio of the confining box exceeds a threshold value, granular hydrostatics predict phase separation: the formation of a high-density region coexisting with a low-density region along the wall that is opposite to the thermal wall. Event-driven molecular dynamic simulations confirm this prediction. The theoretical bifurcation curve agrees with the simulations quantitatively well below and well above the threshold. However, in a wide region of aspect ratios around the threshold, the system is dominated by fluctuations, and the hydrostatic theory breaks down. Two possible scenarios of the origin of the giant fluctuations are discussed.

preprint2003arXiv

Kinetic Integrals in the Kinetic Theory of dissipative gases

The kinetic theory of gases, including Granular Gases, is based on the Boltzmann equation. Many properties of the gas, from the characteristics of the velocity distribution function to the transport coefficients may be expressed in terms of functions of the collision integral which we call kinetic integrals. Although the evaluation of these functions is conceptually straightforward, technically it is frequently rather cumbersome. We report here a method for the analytical evaluation of kinetic integrals using computer algebra. We apply this method for the computation of some properties of Granular Gases, ranging from the moments of the velocity distribution function to the transport coefficients. For their technical complexity most of these quantities cannot be computed manually.

preprint2002arXiv

A Statistical Approach to Vehicular Traffic

A two-dimensional cellular automaton is introduced to model the flow and jamming of vehicular traffic in cities. Each site of the automaton represents a crossing where a finite number of cars can wait approaching the crossing from each of the four directions. The flow of cars obeys realistic traffic rules. We investigate the dependence of the average velocity of cars on the global traffic density. At a critical threshold for the density the average velocity reduces drastically caused by jamming. For the low density regime we provide analytical results which agree with the numerical results.

preprint2002arXiv

Entropy and Long range correlations in literary English

We investigated long range correlations in two literary texts, Moby Dick by H. Melville and Grimm's tales. The analysis is based on the calculation of entropy like quantities as the mutual information for pairs of letters and the entropy, the mean uncertainty, per letter. We further estimate the number of different subwords of a given length n. Filtering out the contributions due to the effects of the finite length of the texts, we find correlations ranging to a few hundred letters. Scaling laws for the mutual information (decay with a power law), for the entropy per letter (decay with the inverse square root of n) and for the word numbers (stretched exponential growth with n and with a power law of the text length) were found.

preprint2002arXiv

Guessing probability distributions from small samples

We propose a new method for the calculation of the statistical properties, as e.g. the entropy, of unknown generators of symbolic sequences. The probability distribution p(k) of the elements k of a population can be approximated by the frequencies f(k) of a sample provided the sample is long enough so that each element k occurs many times. Our method yields an approximation if this precondition does not hold. For a given f(k) we recalculate the Zipf-ordered probability distribution by optimization of the parameters of a guessed distribution. We demonstrate that our method yields reliable results.

preprint2002arXiv

Hydrodynamics of dense granular systems

The properties of dense granular systems are analyzed from a hydrodynamical point of view, based on conservation laws for the particle number density and linear momentum. We discuss averaging problems associated with the nature of such systems and the peculiarities of the sources of noise. We perform a quantitative study by combining analytical methods and numerical results obtained by ensemble-averaging of data on creep during compaction and molecular dynamics simulations of convective flow. We show that numerical integration of the hydrodynamic equations gives the expected evolution for the time-dependent fields.

preprint2002arXiv

Size segregation and convection

The size segregation of granular materials in a vibrating container is investigated using Molecular Dynamics. We find that the rising of larger particles is accompanied by the existence of convection cells even in the case of the lowest possible frequencies. The convection can, however, also be triggered by the larger particle itself. The possibility of rising through this mechanism strongly depends on the depth of the larger particle.