Search

Confined flows of particles can lead to clogging, and therefore failure, of various fluidic systems across many applications. As a result, design guidelines need to be developed to ensure that clogging is prevented or at least delayed. In this Letter, we investigate the influence of the angle of reduction in the cross section of the channel on the bridging of semidilute and dense non-Brownian suspensions of spherical particles. We observe a decrease of the clogging probability with the reduction of the constriction angle. This effect is more pronounced for dense suspensions close to the maximum packing fraction where particles are in contact in contrast to semidilute suspensions. We rationalize this difference in terms of arch selection. We describe the role of the constriction angle and the flow profile, providing insights into the distinct behavior of semidilute and dense suspensions. Published by the American Physical Society 2024

HypothesisThe dip coating of suspensions made of monodisperse non-Brownian spherical particles dispersed in a Newtonian fluid leads to different coating regimes depending on the ratio of the particle diameter to the thickness of the film entrained on the substrate. In particular, dilute particles dispersed in the liquid are entrained only above a threshold value of film thickness. In the case of anisotropic particles, in particular fibers, the smallest characteristic dimension will control the entrainment of the particle. Furthermore, it is possible to control the orientation of the anisotropic particles depending on the substrate geometry. In the thick film regime, the Landau-Levich-Derjaguin model remains valid if one account for the change in viscosity.ExperimentTo test the hypotheses, we performed dip-coating experiments with dilute suspensions of non-Brownian fibers with different length-to-diameter aspect ratios. We characterize the number of fibers entrained on the surface of the substrate as a function of the withdrawal velocity, allowing us to estimate a threshold capillary number below which all the particles remain in the liquid bath. Besides, we measure the angular distribution of the entrained fibers for two different substrate geometries: flat plates and cylindrical rods. We then measure the film thickness for more concentrated fiber suspensions.FindingsThe entrainment of the fibers on a flat plate and a cylindrical rod is primarily controlled by the smaller characteristic length of the fibers: their diameter. At first order, the entrainment threshold scales similarly to that of spherical particles. The length of the fibers only appears to have a minor influence on the entrainment threshold. No preferential alignment is observed for non-Brownian fibers on a flat plate, except for very thin films, whereas the fibers tend to align themselves along the axis of a cylindrical rod for a large enough ratio of the fiber length to the radius of the cylindrical rod. The Landau-Levich-Derjaguin law is recovered for more concentrated suspension by introducing an effective capillary number accounting for the change in viscosity.

Whether it’s pipes, highways, or arteries that are clogged, stopping the flow is always inconvenient and sometimes dangerous.

Debris flows are dense and fast-moving complex suspensions of soil and water that threaten lives and infrastructure. Assessing the hazard potential of debris flows requires predicting yield and flow behavior. Reported measurements of rheology for debris flow slurries are highly variable and sometimes contradictory due to heterogeneity in particle composition and volume fraction ([Formula: see text]) and also inconsistent measurement methods. Here we examine the composition and flow behavior of source materials that formed the postwildfire debris flows in Montecito, CA, in 2018, for a wide range of [Formula: see text] that encapsulates debris flow formation by overland flow. We find that shear viscosity and yield stress are controlled by the distance from jamming, [Formula: see text], where the jamming fraction [Formula: see text] is a material parameter that depends on grain size polydispersity and friction. By rescaling shear and viscous stresses to account for these effects, the data collapse onto a simple nondimensional flow curve indicative of a Bingham plastic (viscoplastic) fluid. Given the highly nonlinear dependence of rheology on [Formula: see text], our findings suggest that determining the jamming fraction for natural materials will significantly improve flow models for geophysical suspensions such as hyperconcentrated flows and debris flows.

Clogging can occur whenever a suspension of particles flows through a confined system. The formation of clogs is often correlated to a reduction in the cross-section of the channel. In this study, we consider the clogging by bridging, i.e., through the formation of a stable arch of particles at a constriction that hinders the transport of particles downstream of the clog. To characterize the role of the volume fraction of the suspension on the clogging dynamics, we study the flow of particulate suspensions through 3D-printed millifluidic devices. We systematically characterize the bridging of non-Brownian particles in a quasi-bidimensional system in which we directly visualize and track the particles as they flow and form arches at a constriction. We report the conditions for clogging by bridging when varying the constriction width to particle diameter ratio for different concentrations of the particles in suspension. We then discuss our results using a stochastic model to rationalize the influence of solid fraction on the probability of clogging. Understanding the mechanisms and conditions of clog formation is an important step for optimizing engineering design and developing more reliable dispensing systems.

In this study, we investigate the transition between the Newtonian and the viscoelastic regimes during the pinch-off of droplets of dilute polymer solutions and discuss its link to the coil-stretch transition. The detachment of a drop from a nozzle is associated with the formation of a liquid neck that causes the divergence of the local stress in a vanishingly small region. If the liquid is a polymer solution, this increasing stress progressively unwinds the polymer chains, up to a point where the resulting increase in the viscosity slows down drastically the thinning. This threshold to a viscoelastic behavior corresponds to a macroscopic strain rate. In the present study, we characterize the variations of with respect to the polymer concentration and molar weight, to the solvent viscosity, and to the nozzle size, i.e., the weight of the drop. We provide empirical scaling laws for these variations. We also analyze the thinning dynamics at the transition and show that it follows a self-similar dynamics controlled by the time scale c-1. This characteristic time is different and always shorter than the relaxation time of the polymer.

Dip coating consists of withdrawing a substrate from a bath to coat it with a thin liquid layer. This process is well understood for homogeneous fluids, but heterogeneities, such as particles dispersed in liquid, lead to more complex situations. Indeed, particles introduce a new length scale, their size, in addition to the thickness of the coating film. Recent studies have shown that, at first order, the thickness of the coating film for monodisperse particles can be captured by an effective capillary number based on the viscosity of the suspension, providing that the film is thicker than the particle diameter. However, suspensions involved in most practical applications are polydisperse, characterized by a wide range of particle sizes, introducing additional length scales. In this study, we investigate the dip coating of suspensions having a bimodal size distribution of particles. We show that the effective viscosity approach is still valid in the regime where the coating film is thicker than the diameter of the largest particles, although bidisperse suspensions are less viscous than monodisperse suspensions of the same solid fraction. We also characterize the intermediate regime that consists of a heterogeneous coating layer and where the composition of the film is different from the composition of the bath. A model to predict the probability of entraining the particles in the liquid film depending on their sizes is proposed and captures our measurements. In this regime, corresponding to a specific range of withdrawal velocities, capillarity filters the large particles out of the film.

SignificanceThe pinch-off of a liquid drop extruded from a nozzle is a canonical situation that involves a series of self-similar regimes ending in a finite-time singularity. This configuration allows for exploring capillary flows over a large range of scales. In the case of suspension drops, the presence of particles breaks the self-similarity by introducing a length scale that can be much larger than the particle diameter. This length scale is a signature of the heterogeneities and delimitates a regime, in which a continuum approach of a suspension can be used from a regime where the discrete nature of the particles is involved.

Clogging is a common obstacle encountered during the transport of suspensions and represents a significant energy and material cost across applications, including water purification, irrigation, biopharmaceutical processing, and aquifer recharge. Pulsatile pressure-driven flows can help mitigate clogging when compared to steady flows. Here, we study experimentally the influence of the amplitude of pulsation 0.25P0 ≤ δP ≤ 1.25P0, where P0 is the mean pressure, and of the frequency of pulsation 10-3 Hz ≤ f ≤ 10-1 Hz on clog mitigation in a microfluidic array of parallel channels using a dilute suspension of colloidal particles. The array geometry is representative of a classical filter, with parallel pores that clog over time, yielding a filter cake that continues to grow and can interact with other pores. We combine flow rate measurements with direct visualizations at the pore scale to correlate the observed clogging dynamics with the changes in flow rate. We observe that all pulsatile amplitudes at 0.1 Hz yield increased throughput compared to steady flows. The rearrangement of particles when subject to a dynamic shear environment can delay the clogging of a pore or even remove an existing clog. However, this benefit is drastically reduced at 10-2 Hz and disappears at 10-3 Hz as the pulsatile timescale becomes too large compared to the timescale associated with the clogging and the growth of the filter cakes in this system. The present study demonstrates that pulsatile flows are a promising method to delay clogging at both the pore and system scale.

Gravity-driven collapses involving large amounts of dense granular material, such as landslides, avalanches, or rock falls, in a geophysical context, represent significant natural hazards. Understanding their complex dynamics is hence a key concern for risk assessment. In the present work, we report experiments on the collapse of quasi-two-dimensional dry granular columns under the effect of gravity, where both the velocity at which the grains are released and the aspect ratio of the column are varied to investigate the dynamics of the falling grains. At high release velocity, classical power laws for the final deposit are recovered, meaning those are representative of a free-fall-like regime. For sufficiently high aspect ratios, the top of the column undergoes an overall free-fall-like motion. In addition, for all experiments, the falling grains also spread horizontally in a free-fall-like motion, and the characteristic time of spreading is related to the horizontal extension reached by the deposit at all altitudes. At low release velocity, a quasistatic state is observed, with scaling laws for the final geometry identical to those of the viscous regime of granular-fluid flow. The velocity at which the grains are released governs the collapse dynamics. Between these two asymptotic regimes, higher release velocity correlates with smaller impact on the collapse dynamics. The criterion V[over ¯]≥0.4sqrt[gH_{0}], where H_{0} is the initial height of the column, is found for the mean release velocity V[over ¯] not to influence the granular collapse.

Abstract:

When a droplet is generated, the ligament connecting the drop to the nozzle thins down and eventually pinches off. Adding solid particles to the liquid phase leads to a more complex dynamic, notably by increasing the shear viscosity. Moreover, it introduces an additional length scale to the system, the diameter of the particles, which eventually becomes comparable to the diameter of the ligament. In this paper, we experimentally investigate the thinning and pinch-off of drops of suspensions with two different sizes of particles. We characterize the thinning for different particle size ratios and different proportions of small particles. Long before the pinch-off, the thinning rate is that of an equivalent liquid whose viscosity is that of the suspension. Later, when the ligament thickness approaches the size of the large particles, the thinning accelerates and leads to an early pinch-off. We explain how the bidisperse particle size distribution lowers the viscosity by making the packing more efficient, which speeds up the thinning. This result can be used to predict the dynamics of droplet formation with bidisperse suspensions.

Abstract:

The formation of drops of a complex fluid, for instance including dissolvedpolymers and/or solid particles, has practical implications in severalindustrial and biophysical processes. In this Letter, we experimentallyinvestigate the generation of drops of a viscoelastic suspension, made ofnon-Brownian spherical particles dispersed in a dilute polymer solution. Usinghigh-speed imaging, we characterize the different stages of the detachment. Ourexperiments show that the particles primarily affect the initial Newtoniannecking by increasing the fluid viscosity. In the viscoelastic regime,particles do not affect the thinning until the onset of the blisteringinstability, which they accelerate. We find that the transition from one regimeto another, which corresponds to the coil-stretch transition of the polymerchains, strongly depends on the particle content. Considering that the presenceof rigid particles increase the deformation of the liquid phase, we propose anexpression for the local strain rate in the suspension, which rationalizes ourexperimental results. This method could enable the precise measurement of localstresses in particulate suspensions.

Abstract:

Modeling of tsunami waves generated by subaerial landslides is important to provide accurate hazard and risk assessments in coastal areas. We perform small-scale laboratory experiments where a tsunami-like wave is generated by the gravity-driven collapse of a subaerial granular column into water. We show that the maximal amplitude reached near-shore by the generated wave in our experiments is linked to the instantaneous immersed volume of grains and to the ultimate immersed deposit. Despite the differences in scale and geometry between our small-scale experiments and the larger-scale geophysical events, a rather good agreement is found between the experimental law and the field data. This approach offers an easy way to estimate the amplitude of paleo-tsunamis.

The withdrawal of a liquid or the translation of a liquid slug in a capillarytube leads to the deposition of a thin film on the inner wall. When particlesor contaminants are present in the liquid, they deposit and contaminate thetube if the liquid film is sufficiently thick. In this article, weexperimentally investigate the condition under which particles are depositedduring the air invasion in a capillary tube initially filled with a dilutesuspension. We show that the entrainment of particles in the film is controlledby the ratio of the particle and the tube radii and the capillary numberassociated with the front velocity. We also develop a model which suggestsoptimal operating conditions to avoid contamination during withdrawal of asuspension. This deposition mechanism can also be leveraged in coatingprocesses by controlling the deposition of particles on the inner walls ofchannels.

The fragmentation of liquid sheets produces a collection of droplets. Thesize distribution of the droplets has a considerable impact on the coatingefficiency of sprays and the transport of contaminants. Although many processescommonly used particulate suspensions, the influence of the particles on thespreading dynamics of the sheet and its subsequent fragmentation has so farbeen considered negligible. In this paper, we consider experimentally atransient suspension sheet that expands radially. We characterize the influenceof the particles on the dynamics of the liquid sheet and the fragmentationprocess. We highlight that the presence of particles modifies the thickness andreduces the stability of the liquid sheet. Our study suggests that particlescan significantly modify the dynamics of liquid films through capillaryeffects, even for volume fractions much smaller than the maximum packing.

We consider the rearrangements of a vertical twodimensional granular packinginduced by the withdrawal of an intruder. Here, we focus on the influence ofthe size of the intruder on the reorganization process. The long term evolutionof the granular packing is investigated as well as the avalanche dynamics thatcharacterize the short term rearrangements around the intruder. For smallenough intruder, we observe the formation of arches that periodicallydestabilize and influence the reorganization dynamics of the two-dimensionalpacking through larger rearrangement events.

An object withdrawn from a liquid bath is coated with a thin layer of liquid.Along with the liquid, impurities such as particles present in the bath can betransferred to the withdrawn substrate. Entrained particles locally modify thethickness of the film, hence altering the quality and properties of thecoating. In this study, we show that it is possible to entrain the liquid aloneand avoid contamination of the substrate, at sufficiently low withdrawalvelocity in diluted suspensions. Using a model system consisting of a plateexiting a liquid bath, we observe that particles can remain trapped in themeniscus which exerts a resistive capillary force to the entrainment. Wecharacterize different entrainment regimes as the withdrawal velocityincreases: from a pure liquid film, to a liquid film containing clusters ofparticles, and eventually individual particles. This capillary filtration is aneffective barrier against the contamination of substrates withdrawn from apolluted bath and finds application against biocontamination.

Withdrawing a plate from a suspension leads to the entrainment of a coating layer of fluid and particles on the solid surface. In this article, we study the Landau-Levich problem in the case of a suspension of non-Brownian particles at moderate volume fraction 10% < φ < 41%. We observe different regimes depending on the withdrawal velocity U, the volume fraction of the suspension φ, and the diameter of the particles 2a. Our results exhibit three coating regimes. (i) At small enough capillary number Ca, no particles are entrained, and only a liquid film coats the plate. (ii) At large capillary number, we observe that the thickness of the entrained film of suspension is captured by the Landau-Levich law using the effective viscosity of the suspension η(φ). (iii) At intermediate capillary numbers, the situation becomes more complicated with a heterogeneous coating on the substrate. We rationalize our experimental findings by providing the domain of existence of these three regimes as a function of the fluid and particles properties.

During the transport of colloidal suspensions in microchannels, thedeposition of particles can lead to the formation of clogs, typically atconstrictions. Once a clog is formed in a microchannel, advected particles forman aggregate upstream from the site of the blockage. This aggregate grows overtime, which leads to a dramatic reduction of the flow rate. In this paper, wepresent a model that predicts the growth of the aggregate formed upon cloggingof a microchannel. We develop an analytical description that captures the timeevolution of the volume of the aggregate, as confirmed by experiments performedusing a pressure-driven suspension flow in a microfluidic device. We show thatthe growth of the aggregate increases the hydraulic resistance in the channeland leads to a drop in the flow rate of the suspensions. We then derive a modelfor the growth of aggregates in multiple parallel microchannels where theclogging events are described using a stochastic approach. The aggregategrowths in the different channels are coupled. Our work illustrates thecritical influence of clogging events on the evolution of the flow rate inmicrochannels. The coupled dynamics of the aggregates described here forparallel channels is key to bridge clogging at the pore scale with macroscopicobservations of the flow rate evolution at the filter scale.

Wet granular aggregates are common precursors of construction materials, food, and health care products. The physical mechanisms involved in the mixing of dry grains with a wet substrate are not well understood and difficult to control. Here, we study experimentally the accretion of dry grains on a wet granular substrate by measuring the growth dynamics of the wet aggregate. We show that this aggregate is fully saturated and its cohesion is ensured by the capillary depression at the air-liquid interface. The growth dynamics is controlled by the liquid fraction at the surface of the aggregate and exhibits two regimes. In the viscous regime, the growth dynamics is limited by the capillary-driven flow of liquid through the granular packing to the surface of the aggregate. In the capture regime, the capture probability depends on the availability of the liquid at the saturated interface, which is controlled by the hydrostatic depression in the material. We propose a model that rationalizes our observations and captures both dynamics based on the evolution of the capture probability with the hydrostatic depression.

Various materials are made of long thin fibers that are randomly oriented to form a complex network in which drops of wetting liquid tend to accumulate at the nodes. The capillary force exerted by the liquid can bend flexible fibers, which in turn influences the morphology adopted by the liquid. In this paper, we investigate through a model situation the role of the fiber flexibility on the shape of a small volume of liquid on a pair of crossed flexible fibers. We characterize the liquid morphologies as we vary the volume of liquid, the angle between the fibers, and the length of the fibers. The drop morphologies previously reported for rigid crossed fibers, i.e., a drop, a column and a mixed morphology, are also observed on flexible crossed fibers with modified domains of existence. In addition, at small tilt angles between the fibers, a new behavior is observed: the fibers bend and collapse. Depending on the volume of liquid, a thin column with or without a drop is reported on the collapsed fibers. Our study suggests that the fiber flexibility adds a rich variety of behaviors that may be important for some applications.

Physical and chemical transformation processes in reactive granular media involve the reorganization of the structure. In this paper, we study experimentally the rearrangements of a two-dimensional (2D) granular packing undergoing a localized transformation. We track the position and evolution of all the disks that constitute the granular packing when either a large intruder shrinks in size or is pulled out of the granular structure. In the two situations the displacements at long time are similar to 2D quasistatic silo flows whereas the short-time dynamic is heterogeneous in both space and time. We observe an avalanchelike behavior with power-law distributed events uncorrelated in time. In addition, the instantaneous evolutions of the local solid fraction exhibit self-similar distributions. The averages and the standard deviations of the solid fraction variations can be rescaled, suggesting a single mechanism of rearrangement.

When droplets impact fibrous media, the liquid can be captured by the fibers or contact then break away. Previous studies have shown that the efficiency of drop capture by a rigid fiber depends on the impact velocity and a threshold velocity was defined below which the drop is captured. However, it is necessary to consider the coupling of elastic and capillary effects to achieve an improved understanding of the capture process for soft substrates. Here, we study experimentally the dynamics of a single drop impacting on a thin flexible fiber. Our results demonstrate that the threshold capture velocity depends on the flexibility of fibers in a non-monotonic way. We conclude that tuning the mechanical properties of fibers can optimize the efficiency of droplet capture.

Transport of liquid drops in fibrous media occurs in various engineering systems such as fog harvesting or cleaning of textiles. The ability to tune or to control liquid movement can increase the system efficiency and enable new engineering applications. In this letter, we experimentally investigate how partially wetting drops on a single fiber can be manipulated by vibrating the fiber. We show that a sliding motion along the fiber or a dripping of the drop can be triggered by standing waves. We identify the conditions on the drop volume, the fiber tilt angle, and the amplitude and frequency of oscillations to observe these different behaviors. Finally, we experimentally illustrate that vibrations can be used to control the transport and the collection of water drops along a fiber using a combination of the sliding and dripping transitions.

We characterize the different morphologies adopted by a drop of liquid placed on two randomly oriented fibers, which is a first step toward understanding the wetting of fibrous networks. The present work reviews previous modeling for parallel and touching crossed fibers and extends it to an arbitrary orientation of the fibers characterized by the tilting angle and the minimum spacing distance. Depending on the volume of liquid, the spacing distance between fibers and the angle between the fibers, we highlight that the liquid can adopt three different equilibrium morphologies: 1) a column morphology in which the liquid spreads between the fibers, 2) a mixed morphology where a drop grows at one end of the column or 3) a single drop located at the node. We capture the different morphologies observed using an analytical model that predicts the equilibrium configuration of the liquid based on the geometry of the fibers and the volume of liquid.

We investigate the equilibrium morphology of a finite volume of liquid placed on two parallel rigid fibers of different radii. As observed for identical radii fibers, the liquid is either in a column morphology or adopts a drop shape depending on the inter-fiber distance. However the cross-sectional area and the critical inter-fiber distance at which the transition occurs are both modified by the polydispersity of the fibers. Using energy considerations, we analytically predict the critical inter-fiber distance corresponding to the transition between the column and the drop morphologies. This distance depends both on the radii of the fibers and on the contact angle of the liquid. We perform experiments using a perfectly wetting liquid on two parallel nylon fibers: the results are in good agreement with our analytical model. The morphology of the capillary bridges between fibers of different radii is relevant to the modeling of large arrays of polydisperse fibers.

We report on a microfluidic method that allows measurement of a small concentration of large contaminants in suspensions of solid micrometer-scale particles. To perform the measurement, we flow the colloidal suspension through a series of constrictions, i.e., a microchannel of varying cross-section. We show and quantify the role of large contaminants in the formation of clogs at a constriction and the growth of the resulting filter cake. By measuring the time interval between two clogging events in an array of parallel microchannels, we are able to estimate the concentration of contaminants whose size is selected by the geometry of the microfluidic device. This technique for characterizing colloidal suspensions offers a versatile and rapid tool to explore the role of contaminants on the properties of the suspensions.

We investigate the wetting properties of the simplest element of an array ofrandom fibers: two rigid fibers crossing with an inclination angle and incontact with a droplet of a perfectly wetting liquid. We show experimentallythat the liquid adopts different morphologies when the inclination angle isincreased: a column shape, a mixed morphology state where a drop lies at theend of a column, or a drop centered at the node. An analytical model isprovided that predicts the wetting length as well as the presence of anon-symmetric state in the mixed morphology regime. The model also highlights asymmetry breaking at the transition between the column state and the mixedmorphology. The possibility to tune the morphology of the liquid could haveimportant implications for drying processes.

We report a new method to display the minute fluctuations induced by syringe pumps on microfluidic flows by using a liquid-liquid system with an ultralow interfacial tension. We demonstrate that the stepper motor inside the pump is a source of fluctuations in microfluidic flows by comparing the frequencies of the ripples observed at the interface to that of the pulsation of the stepper motor. We also quantify the fluctuations induced at different flow rates, using syringes of different diameters, and using different syringe pumps with different advancing distances per step. Our work provides a way to predict the frequency of the fluctuation that the driving syringe pump induces on a microfluidic system and suggests that syringe pumps can be a source of fluctuations in microfluidic flows, thus contributing to the polydispersity of the resulting droplets.

Immiscible aqueous phases, formed by dissolving incompatible solutes in water, have been used in green chemical synthesis, molecular extraction and mimicking of cellular cytoplasm. Recently, a microfluidic approach has been introduced to generate all-aqueous emulsions and jets based on these immiscible aqueous phases; due to their biocompatibility, these all-aqueous structures have shown great promises as templates for fabricating biomaterials. The physico-chemical nature of interfaces between two immiscible aqueous phases leads to unique interfacial properties, such as an ultra-low interfacial tension. Strategies to manipulate components and direct their assembly at these interfaces needs to be explored. In this paper, we review progress on the topic over the past few years, with a focus on the fabrication and stabilization of all-aqueous structures in a multiphase microfluidic platform. We also discuss future efforts needed from the perspectives of fluidic physics, materials engineering, and biology for fulfilling potential applications ranging from materials fabrication to biomedical engineering.

Copyright of Reflets de la Physique is the property of EDP Sciences and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

Tsunami-like waves can arise from large-scale geological events such as the collapse of a mountain flank or a volcano into a lake or an ocean. We here mimic this situation at the laboratory scale. The landslide is modeled by the sudden release of a rectangular granular column, which impacts a still water layer and generates a wave. Through experiments varying both the column dimensions and the initial water depth, three regimes of nonlinear waves of different shapes are reported. The amplitude and width of the generated waves are governed by a local Froude number based on the velocity of the granular front at the water surface.

The erosion of a cohesive soil by an impinging turbulent jet is observed, for instance, during the landing of a spacecraft or involved in the so-called jet erosion test. To provide a quantitative understanding of this situation for cohesive soils, we perform experiments using a model cohesion controlled granular material that allows us to finely tune the cohesion between particles while keeping the other properties constant. We investigate the response of this cohesive granular bed when subjected to an impinging normal turbulent jet. We characterize experimentally the effects of the cohesion on the erosion threshold and the development of the crater. We demonstrate that the results can be rationalized by introducing a cohesive Shields number that accounts for the inter-particles cohesion force. The results of our experiments highlight the crucial role of cohesion in erosion processes.