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

Markus Uhlmann

Markus Uhlmann 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.softmath.DSphysics.ao-phphysics.geo-ph

Catalog footprint

What is connected

19works
5topics
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

19 published item(s)

preprint2021arXiv

Direct numerical simulation of turbulent mass transfer at the surface of an open channel flow

We present direct numerical simulation results of turbulent open channel flow at bulk Reynolds numbers up to 12000, coupled with (passive) scalar transport at Schmidt numbers up to 200. Care is taken to capture the very large scale motions which appear already for relatively modest Reynolds numbers. The transfer velocity at the flat, free surface is found to scale with the Schmidt number to the power "-1/2", in accordance with previous studies and theoretical predictions for uncontaminated surfaces. The scaling of the transfer velocity with Reynolds number is found to vary, depending on the Reynolds number definition used. To compare the present results with those obtained in other systems, we define a turbulent Reynolds number at the edge of the surface-influenced layer. This allows us to probe the two-regime model of Theofanous [Turbulent mass transfer at free, gas-liquid interfaces, with applications to open-channel, bubble and jet flows. Int. J. Heat Mass Transfer 19, 613--624, 1976], which is found to correctly predict that small-scale vortices significantly affect the mass transfer for turbulent Reynolds numbers larger than 500. It is further established that the root-mean-square of the surface divergence is, on average, proportional to the mean transfer velocity. However, the spatial correlation between instantaneous surface divergence and transfer velocity tends to decrease with increasing Schmidt number and increase with increasing Reynolds number. The latter is shown to be caused by an enhancement of the correlation in high-speed regions, which in turn is linked to the spatial distribution of surface-parallel vortices.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Can preferential concentration of finite-size particles in plane Couette turbulence be reproduced with the aid of equilibrium solutions?

This work employs for the first time invariant solutions of the Navier-Stokes equations to study the interaction between finite-size particles and near-wall coherent structures. We consider horizontal plane Couette flow and focus on Nagata's upper-branch equilibrium solution (Nagata, 1990) at low Reynolds numbers where this solution is linearly stable. When adding a single heavy particle with a diameter equivalent to 2.5 wall units, we observe that the solution remains stable and is essentially unchanged away from the particle. This result demonstrates that it is technically feasible to utilize exact coherent structures in conjunction with particle-resolved DNS. While translating in the streamwise direction, the particle migrates laterally under the action of the quasi-streamwise vortices until it reaches the region occupied by the low-speed streak, where it attains a periodic state of motion. As a result of the ensuing preferential particle location, the time-average streamwise particle velocity differs from the plane-average fluid-phase velocity at the same wall-distance as the particle center, as previously observed in experiments and in numerical data for fully turbulent wall-bounded flows. Additional constrained simulations where the particle is maintained at a fixed spanwise position while freely translating in the other two directions reveal the existence of two equilibria located in the low-speed and in the high-speed streak, respectively, the former being an unstable point. A parametric study with different particle to fluid density ratios is conducted which shows how inertia affects the spanwise fluctuations of the periodic particle motion. Finally, we discuss a number of potential future investigations of solid particle dynamics which can be conducted with the aid of invariant solutions (exact coherent structures) of the Navier-Stokes equations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the ice-nucleating potential of warm hydrometeors in mixed-phase clouds

The question whether or not the presence of warm hydrometeors in clouds may play a significant role in the nucleation of new ice particles has been debated for several decades. While the early works of Fukuta and Lee (1986) and Baker (1991) indicated that it might be irrelevant, the more recent study of Prabhakaran et al. (2019) [arXiv:1906.06129] suggested otherwise. In this work, we are aiming to quantify the ice-nucleating potential using high-fidelity flow simulation techniques around a single hydrometeor and use favorable considerations to upscale the effects to a collective of ice particles in clouds. While we find that ice nucleation may be enhanced in the vicinity of a warm hydrometeor by several orders of magnitude and that the affected volume of air is much larger than previously estimated, it is very unlikely that this effect alone causes the rapid enhancement of ice nucleation observed in some types of clouds, mainly due to the low initial volumetric ice concentration. Nonetheless, it is suggested to implement this effect into existing cloud models in order to investigate second-order effects such as ice nucleus preactivation or enhancement after the onset of glaciation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the scaling of the instability of a flat sediment bed with respect to ripple-like patterns

We investigate the formation of subaqueous transverse bedforms in turbulent open channel flow by means of direct numerical simulations with fully-resolved particles. The main goal of the present analysis is to address the question whether the initial pattern wavelength scales with the particle diameter or with the mean fluid height. A previous study (Kidanemariam and Uhlmann, J. Fluid Mech., vol. 818, 2017, pp. 716-743) has observed a lower bound for the most unstable pattern wavelength in the range 75-100 times the particle diameter, which was equivalent to 3-4 times the mean fluid height. In the current paper, we vary the streamwise box length in terms of the particle diameter and of the mean fluid height independently in order to distinguish between the two possible scaling relations. For the chosen parameter range, the obtained results clearly exhibit a scaling of the initial pattern wavelength with the particle diameter, with a lower bound around a streamwise extent of approximately 80 particle diameters. In longer domains, on the other hand, patterns are observed at initial wavelengths in the range 150-180 times the particle diameter, which is in good agreement with experimental measurements. Variations of the mean fluid height, on the other hand, seem to have no significant influence on the most unstable initial pattern wavelength. Furthermore, for the cases with the largest relative submergence, we observe spanwise and streamwise sediment waves of similar amplitude to evolve and superimpose, leading to three-dimensional sediment patterns.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Voronoi tesselation analysis of sets of randomly placed finite-size spheres

The purpose of this note is to clarify the effect of the finite size of spherical particles upon the characteristics of their spatial distribution through a random Poisson process (RPP). This information is of special interest when using RPP data as a reference for the analysis of the spatial structure of a given (non-RPP) particulate system, in which case ignoring finite-size effects upon the former may yield misleading conclusions. We perform Monte Carlo simulations in triply-periodic spatial domains, and then analyze the particle-centered Voronoi tesselations for solid volume fractions ranging from 10^(-5) to 0.3. We show that the standard-deviation of these volumes decreases with the solid volume fraction, the deviation from the value of point sets being reasonably approximated by an exponential function. As can be expected, the domain size for which the random assemblies of finite-size particles are generated has a constraining effect if the number of particles per realization is chosen too small. This effect is quantified, and recommendations are given. We have also revisited the case of random point sets (i.e. the limit of vanishing particle diameter), for which we have confirmed the accuracy of the earlier data by Tanemura [Forma, 18(4):221-247, 2003].

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Interface-resolved direct numerical simulations of sediment transport in a turbulent oscillatory boundary layer

The flow within an oscillatory boundary layer, which approximates the flow generated by propagating sea waves of small amplitude close to the bottom, is simulated numerically by integrating Navier-Stokes and continuity equations. The bottom is made up of spherical particles, free to move, which mimic sediment grains. The approach allows to fully-resolve the flow around the particles and to evaluate the forces and torques that the fluid exerts on their surface. Then, the dynamics of sediments is explicitly computed by means of Newton-Euler equations. For the smallest value of the flow Reynolds number presently simulated, the flow regime turns out to fall in the intermittently turbulent regime such that turbulence appears when the free stream velocity is close to its largest values but the flow recovers a laminar like behaviour during the remaining phases of the cycle. For the largest value of the Reynolds number turbulence is significant almost during the whole flow cycle. The evaluation of the sediment transport rate allows to estimate the reliability of the empirical predictors commonly used to estimate the amount of sediments transported by the sea waves. For large values of the Shields parameter, the sediment flow rate during the accelerating phases does not differ from that observed during the decelerating phases. However, for relatively small values of the Shields parameter, the amount of moving particles depends not only on the bottom shear stress but also on flow acceleration. Moreover, the numerical results provide information on the role that turbulent eddies have on sediment dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Columnar structure formation of a dilute suspension of settling spherical particles in a quiescent fluid

The settling of heavy spherical particles in a column of quiescent fluid is investigated. The performed experiments cover a range of Galileo numbers ($110 \leq \text{Ga} \leq 310$) for a fixed density ratio of $Γ= ρ_p/ρ_f = 2.5$. In this regime the particles are known (M. Jenny, J. Dušek and G. Bouchet, Journal of Fluid Mechanics 508, 201 (2004).) to show a variety of motions. It is known that the wake undergoes several transitions for increasing $\text{Ga}$ resulting in particle motions that are successively: vertical, oblique, oblique oscillating, and finally chaotic. Not only does this change the trajectory of single, isolated, settling particles, but it also changes the dynamics of a swarm of particles as collective effects become important even for dilute suspensions, with volume fraction up to $Φ_V = \mathcal{O}\left(10^{-3}\right)$, which are investigated in this work. Multi-camera recordings of settling particles are recorded and tracked over time in 3 dimensions. A variety of analysis are performed and show a strong clustering behavior. The distribution of the cell areas of the Voronoï tessellation in the horizontal plane are compared to that of a random distribution of particles and shows clear clustering. Moreover, a negative correlation was found between the Voronoï area and the particle velocity; clustered particles fall faster. In addition, the angle between two adjacent particles and the vertical is calculated and compared to a homogeneous distribution of particles, clear evidence of vertical alignment of particles is found. The experimental findings are compared to simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Forcing homogeneous turbulence in DNS of particulate flow with interface resolution and gravity

We consider the case of finite-size spherical particles which are settling under gravity in a homogeneous turbulent background flow. Turbulence is forced with the aid of the random forcing method of Eswaran and Pope [Comput. Fluids, 16(3):257-278, 1988], while the solid particles are represented with an immersed-boundary method. The forcing scheme is used to generate isotropic turbulence in vertically elongated boxes in order to warrant better decorrelation of the Lagrangian signals in the direction of gravity. Since only a limited number of Fourier modes are forced, it is possible to evaluate the forcing field directly in physical space, thereby avoiding full-size transforms. The budget of box-averaged kinetic energy is derived from the forced momentum equations. Medium-sized simulations for dilute suspensions at low Taylor-scale Reynolds number $Re_λ=65$, small density ratio $ρ_p/ρ_f=1.5$ and for two Galileo numbers $Ga=0$ and 120 are carried out over long time intervals in order to exclude the possibility of slow divergence. It is shown that the results at zero gravity are fully consistent with previous experimental measurements and available numerical reference data. Specific features of the finite-gravity case are discussed with respect to a reduction of the average settling velocity, the acceleration statistics and the Lagrangian auto-correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Direct numerical simulation of pattern formation in subaqueous sediment

We present results of direct numerical simulation of incompressible fluid flow over a thick bed of mobile, spherically-shaped particles. The algorithm is based upon the immersed boundary technique for fluid-solid coupling and uses a soft-sphere model for the solid-solid contact. Two parameter points in the laminar flow regime are chosen, leading to the emergence of sediment patterns classified as `small dunes', while one case under turbulent flow conditions leads to `vortex dunes' with significant flow separation on the lee side. Wavelength, amplitude and propagation speed of the patterns extracted from the spanwise-averaged fluid-bed interface are found to be consistent with available experimental data. The particle transport rates are well represented by available empirical models for flow over a plane sediment bed in both the laminar and the turbulent regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Interface-resolved direct numerical simulation of the erosion of a sediment bed sheared by laminar channel flow

A numerical method based upon the immersed boundary technique for the fluid-solid coupling and on a soft-sphere approach for solid-solid contact is used to perform direct numerical simulation of the flow-induced motion of a thick bed of spherical particles in a horizontal plane channel. The collision model features a normal force component with a spring and a damper, as well as a damping tangential component, limited by a Coulomb friction law. The standard test case of a single particle colliding perpendicularly with a horizontal wall in a viscous fluid is simulated over a broad range of Stokes numbers, yielding values of the effective restitution coefficient in close agreement with experimental data. The case of bedload particle transport by laminar channel flow is simulated for 24 different parameter values covering a broad range of the Shields number. Comparison of the present results with reference data from the experiment of Aussillous et al. (J. Fluid Mech. 2013) yields excellent agreement. It is confirmed that the particle flow rate varies with the third power of the Shields number once the known threshold value is exceeded. The present data suggests that the thickness of the mobile particle layer (normalized with the height of the clear fluid region) increases with the square of the normalized fluid flow rate.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sedimentation of a dilute suspension of rigid spheres at intermediate Galileo numbers: the effect of clustering upon the particle motion

Direct numerical simulation of the gravity-induced settling of finite-size particles in triply-periodic domains has been performed under dilute conditions. For a single solid-to-fluid density ratio of 1.5 we have considered two values of the Galileo number corresponding to steady vertical motion (Ga=121) and to steady oblique motion (Ga=178) in the case of one isolated sphere. For the multi-particle system we observe strong particle clustering only in the latter case. The geometry and time scales related to clustering are determined from Voronoi tesselation and particle-conditioned averaging. As a consequence of clustering, the average particle settling velocity is increased by 12% as compared to the value of an isolated sphere; such a collective effect is not observed in the non-clustering case. By defining a local (instantaneous) fluid velocity average in the vicinity of the finite-size particles it is shown that the observed enhancement of the settling velocity is due to the fact that the downward fluid motion (with respect to the global average) which is induced in the cluster regions is preferentially sampled by the particles. It is further observed that the variance of the particle velocity is strongly enhanced in the clustering case. With the aid of a decomposition of the particle velocity it is shown that this increase is due to enhanced fluid velocity fluctuations (due to clustering) in the vicinity of the particles. Finally, we discuss a possible explanation for the observation of a critical Galileo number marking the onset of clustering under dilute conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

DNS of horizontal open channel flow with finite-size, heavy particles at low solid volume fraction

We have performed direct numerical simulation of turbulent open channel flow over a smooth horizontal wall in the presence of finite-size, heavy particles. The spherical particles have a diameter of approximately 7 wall units, a density of 1.7 times the fluid density and a solid volume fraction of 0.0005. The value of the Galileo number is set to 16.5, while the Shields parameter measures approximately 0.2. Under these conditions, the particles are predominantly located in the vicinity of the bottom wall, where they exhibit strong preferential concentration which we quantify by means of Voronoi analysis and by computing the particle-conditioned concentration field. As observed in previous studies with similar parameter values, the mean streamwise particle velocity is smaller than that of the fluid. We propose a new definition of the fluid velocity "seen" by finite-size particles based on an average over a spherical surface segment, from which we deduce in the present case that the particles are instantaneously lagging the fluid only by a small amount. The particle-conditioned fluid velocity field shows that the particles preferentially reside in the low-speed streaks, leading to the observed apparent lag. Finally, a vortex eduction study reveals that spanwise particle motion is significantly correlated with the presence of vortices with the corresponding sense of rotation which are located in the immediate vicinity of the near-wall particles.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spatial and temporal scales of force and torque acting on wall-mounted spherical particles in open channel flow

Data from direct numerical simulation of open channel flow over a geometrically rough wall at a bulk Reynolds number of 2900, generated by Chan-Braun et al. ["Force and torque acting on particles in a transitionally rough open-channel flow", J. Fluid Mech. 684, 441--474 (2011), 10.1017/jfm.2011.311] are further analysed with respect to the time and length scales of force and torque acting on the wall-mounted spheres. For the two sizes of spheres in a square arrangement (11 and 49 wall units in diameter, yielding hydraulically smooth and transitionally rough flow, respectively), the spatial structure of drag, lift and spanwise torque is investigated. The auto-correlation and spectra in time as well as the space-time correlation and convection velocities are presented and discussed. It is found that the statistics of spanwise particle torque are similar to those of shear stress at a smooth wall. Particle drag and lift are shown to differ from spanwise particle torque, exhibiting considerably smaller time and length scales; the convection velocities of drag and lift are somewhat larger than those of spanwise torque. Furthermore, correlations between the flow field and particle-related quantities are presented. The spatial structure of the correlation between streamwise velocity and drag/spanwise torque features elongated shapes reminiscent of buffer-layer streaks. The correlation between the pressure field and the particle drag exhibits two opposite-signed bulges on the upstream and downstream sides of a particle.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The motion of a single heavy sphere in ambient fluid: a benchmark for interface-resolved particulate flow simulations with significant relative velocities

Detailed data describing the motion of a rigid sphere settling in unperturbed fluid is generated by means of highly-accurate spectral/spectral-element simulations with the purpose of serving as a future benchmark case. A single solid-to-fluid density ratio of 1.5 is chosen, while the value of the Galileo number is varied from 144 to 250 such as to cover the four basic regimes of particle motion (steady vertical, steady oblique, oscillating oblique, chaotic). This corresponds to a range of the particle Reynolds number from 185 to 365. In addition to the particle velocity data, extracts of the fluid velocity field are provided, as well as the pressure distribution on the sphere's surface. Furthermore, the same solid-fluid system is simulated with a particular non-boundary-conforming approach, i.e. the immersed boundary method proposed by Uhlmann (JCP 2005), using various spatial resolutions. It is shown that the current benchmark case allows to adjust the resolution requirements for a given error tolerance in each flow regime.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

DNS of vertical plane channel flow with finite-size particles: Voronoi analysis, acceleration statistics and particle-conditioned averaging

We have performed a direct numerical simulation of dilute turbulent particulate flow in a vertical plane channel, fully resolving the phase interfaces. The flow conditions are the same as those in the main case of "Uhlmann, M., Phys. Fluids, vol. 20, 2008, 053305", with the exception of the computational domain length which has been doubled in the present study. The statistics of flow and particle motion are not significantly altered by the elongation of the domain. The large-scale columnar-like structures which had previously been identified do persist and they are still only marginally decorrelated in the prolonged domain. Voronoi analysis of the spatial particle distribution shows that the state of the dispersed phase can be characterized as slightly more ordered than random tending towards a homogeneous spatial distribution. It is also found that the p.d.f.'s of Lagrangian particle accelerations for wall-normal and spanwise directions follow a lognormal distribution as observed in previous experiments of homogeneous flows. The streamwise component deviates from this law presenting significant skewness. Finally, a statistical analysis of the flow in the near field around the particles reveals that particle wakes present two regions, a near wake where the velocity deficit decays as 1/x and a far wake with a decay of approximately 1/(x*x).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Force and torque acting on particles in a transitionally rough open channel flow

Direct numerical simulation of open channel flow over a geometrically rough wall has been performed at a bulk Reynolds number of approximately 2900. The wall consisted of a layer of spheres in a square arrangement. Two cases have been considered. In the first case the spheres are small (with diameter equivalent to 10.7 wall units) and the limit of the hydraulically smooth flow regime is approached. In the second case the spheres are more than three times larger (49.3 wall units) and the flow is in the transitionally rough flow regime. Special emphasis is given on the characterisation of the force and torque acting on a particle due to the turbulent flow. It is found that in both cases the mean drag, lift and spanwise torque are to a large extent produced at the top region of the particle surface. The intensity of the particle force fluctuations is significantly larger in the large-sphere case, while the trend differs for the fluctuations of the individual components of the torque. A simplified model is used to show that the torque fluctuations might be explained by the spheres acting as a filter with respect to the size of the flow scales which can effectively generate torque fluctuations. Fluctuations of both force and torque are found to exhibit strongly non-Gaussian probability density functions with particularly long tails, an effect which is more pronounced in the small-sphere case. Some implications of the present results for sediment erosion are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Interface-resolved DNS of vertical particulate channel flow in the turbulent regime

We have conducted a direct numerical simulation (DNS) study of dilute turbulent particulate flow in a vertical plane channel, considering thousands of finite-size rigid particles with resolved phase interfaces. The particle diameter corresponds to approximately 11 wall units and their terminal Reynolds number is set to 136. The fluid flow with bulk Reynolds number 2700 is directed upward, which maintains the particles suspended upon average. Two density ratios were simulated, differing by a factor of 4.5. The corresponding Stokes numbers of the two flow cases were O(10) in the near-wall region and O(1) in the outer flow. We have observed the formation of large-scale elongated streak-like structures with streamwise dimensions of the order of 8 channel half-widths and cross-stream dimensions of the order of one half-width. At the same time, we have found no evidence of significant formation of particle clusters, which suggests that the large structures are due to an intrinsic instability of the flow, triggered by the presence of the particles. It was found that the mean fluid velocity profile tends towards a concave shape, and the turbulence intensity as well as the normal stress anisotropy are strongly increased. The effect of varying the Stokes number while maintaining the buoyancy, particle size and volume fraction constant was relatively weak.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Significance of Simple Invariant Solutions in Turbulent Flows

Recent remarkable progress in computing power and numerical analysis is enabling us to fill a gap in the dynamical systems approach to turbulence. One of the significant advances in this respect has been the numerical discovery of simple invariant sets, such as nonlinear equilibria and periodic solutions, in well-resolved Navier--Stokes flows. This review describes some fundamental and practical aspects of dynamical systems theory for the investigation of turbulence, focusing on recently found invariant solutions and their significance for the dynamical and statistical characterization of low-Reynolds-number turbulent flows. It is shown that the near-wall regeneration cycle of coherent structures can be reproduced by such solutions. The typical similarity laws of turbulence, i.e. the Prandtl wall law and the Kolmogorov law for the viscous range, as well as the pattern and intensity of turbulence-driven secondary flow in a square duct can also be represented by these simple invariant solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Travelling-waves consistent with turbulence-driven secondary flow in a square duct

We present numerically determined travelling-wave solutions for pressure-driven flow through a straight duct with a square cross-section. This family of solutions represents typical coherent structures (a staggered array of counter-rotating streamwise vortices and an associated low-speed streak) on each wall. Their streamwise average flow in the cross-sectional plane corresponds to an eight vortex pattern much alike the secondary flow found in the turbulent regime.

0 citations0 reviewsNo code linkedOpen access