This dissertation consists of three research projects of kinetic models: a structure preserving scheme for Poisson-Nernst-Planck equations and two efficient spectral methods for multi-dimensional Boltzmann equation.<br><br>The Poisson-Nernst-Planck (PNP) equations is widely used to describe the dynamics of ion transport in ion channels. We introduce a structure-preserving semi-implicit finite difference scheme for the PNP equations in a bounded domain. A general boundary condition for the Poisson equation is considered. The fully discrete scheme is shown to satisfy the following properties: mass conservation, unconditional positivity, and energy dissipation (hence preserving the steady-state). <br><br>Numerical approximation of the Boltzmann equation presents a challenging problem due to its high-dimensional, nonlinear, and nonlocal collision operator. Among the deterministic methods, the Fourier-Galerkin spectral method stands out for its relative high accuracy and possibility of being accelerated by the fast Fourier transform. In this dissertation, we studied the state of the art in the fast Fourier method and discussed its limitation. Next, we proposed a new approach to implement the Fourier method, which can resolve those issues. <br><br>However, the Fourier method requires a domain truncation which is unphysical since the collision operator is defined in whole space R^d . In the last part of this dissertation, we introduce a Petrov-Galerkin spectral method for the Boltzmann equation in the unbounded domain. The basis functions (both test and trial functions) are carefully chosen mapped Chebyshev functions to obtain desired convergence and conservation properties. Furthermore, thanks to the close relationship of the Chebyshev functions and the Fourier cosine series, we can construct a fast algorithm with the help of the non-uniform fast Fourier transform (NUFFT).<br>

<p>Granular systems are ensembles of inelastic particles which dissipate energy during collisions. Granular systems serve as excellent models for a wide variety of materials such as sand, soils, corn, and powder. A rather remarkable property of granular systems is when excited, whether due to an interstitial fluid or via the boundaries, the granular particlesdisplay fluid-like behavior. As a result, there has been decades of granular research with the overarching goal of formulating a general granular hydrodynamic theory.</p> <p>However, the granular hydrodynamic theory is limited, and the underlying transport coefficients often require modifications which are based on empirical observations, and assuch, are system-specific. It is ideally better to devise a general theory which minimizes the information needed about the systema priori. The main thrust of the work undertaken shown here strives to develop such a model by using kinetic theory as the basis. More specifically, I investigate granular gases via the direct simulation Monte Carlo (DSMC) methodand modify the governing equations. In this thesis, two idealized cases of granular gases areconsidered: the homogeneous cooling state and a boundary-heated gas (or the pure conduc-tion case). In the former, the effects of polydispersity are probed. In the latter, the evolutionof the local hydrodynamics due to strong rarefaction effects are divulged. Additionally, amodified, more generalized constitutive relation for the heat flux is proposed and comparedwith DSMC results. Extensions of the DSMC method for dense granular gases and granulargases composed of non-spherical particles are also discussed.</p>

Nous étudions expérimentalement l'effondrement d'une colonne granulaire dans un liquide visqueux. Contrairement au cas sec, la morphologie des dépôts n'est principalement plus contrôlée par le rapport d'aspect initial du tas mais par la fraction volumique initiale de la masse granulaire. Deux régimes différents sont identifiés selon l'empilement initial. L'empilement lâche donne lieu à des dépôts minces et longs et la dynamique est rapide. Une surpression du liquide est mesurée sous de la colonne. Pour l'empilement dense, l'étalement final est deux fois moindre, le mouvement est lent et une dépression interstitielle est mesurée. Ces observations suggèrent que la dynamique de l'effondrement granulaire dans un fluide est fortement affectée par le comportement de la dilatance du milieu granulaire.Nous développons ensuite un modèle théorique basé sur des équations diphasiques moyennées dans l’épaisseur prenant en compte les mécanismes de dilatance. L’étude dimensionnelle de notre modèle permet de montrer que l’effondrement d’une colonne est contrôlé par trois paramètres sans dimension : le rapport d’aspect de la colonne, la fraction volumique initiale, et le nombre de grains dans l’épaisseur. On montre également que le temps caractéristique met en compétition le frottement visqueux et la gravité.De ce modèle, nous développons un algorithme de résolution lagrangien. Cette approche, grossière mais robuste, permet d’implanter s sans trop de difficulté. Le code est validé sur des configurations simples sur plan incliné avant de simuler l’effondrement de colonnes granulaires immergées dans la même gamme de paramètres que nos expériences. / Nous étudions expérimentalement l'effondrement d'une colonne granulaire dans un liquide visqueux. Contrairement au cas sec, la morphologie des dépôts n'est principalement plus contrôlée par le rapport d'aspect initial du tas mais par la fraction volumique initiale de la masse granulaire. Deux régimes différents sont identifiés selon l'empilement initial. L'empilement lâche donne lieu à des dépôts minces et longs et la dynamique est rapide. Une surpression du liquide est mesurée sous de la colonne. Pour l'empilement dense, l'étalement final est deux fois moindre, le mouvement est lent et une dépression interstitielle est mesurée. Ces observations suggèrent que la dynamique de l'effondrement granulaire dans un fluide est fortement affectée par le comportement de la dilatance du milieu granulaire.Nous développons ensuite un modèle théorique basé sur des équations diphasiques moyennées dans l’épaisseur prenant en compte les mécanismes de dilatance. L’étude dimensionnelle de notre modèle permet de montrer que l’effondrement d’une colonne est contrôlé par trois paramètres sans dimension : le rapport d’aspect de la colonne, la fraction volumique initiale, et le nombre de grains dans l’épaisseur. On montre également que le temps caractéristique met en compétition le frottement visqueux et la gravité.De ce modèle, nous développons un algorithme de résolution lagrangien. Cette approche, grossière mais robuste, permet d’implanter s sans trop de difficulté. Le code est validé sur des configurations simples sur plan incliné avant de simuler l’effondrement de colonnes granulaires immergées dans la même gamme de paramètres que nos expériences.

A granular material is a collection of discrete, solid particles of macroscopic size dispersed in an interstitial fluid, in which the fluid has an insignificant effect on the particle dynamics. Because they exhibit fascinating properties because of dissipative interactions, due to their importance in geophysical and industrial processes, flows of granular materials have been the focus of large amount of research involving physicists and engineers. A good understanding of the physics of granular materials is desired in order to design efficient processing and handling systems. Granular materials can be heaped like a solid, and can flow like a fluid. Though the two distinct regimes of granular flows are well described by kinetic theory (rapid flows) and plasticity theories (quasi-static), the intermediate dense flow regime, where collisional and frictional interactions are important, is not yet described successfully. In this thesis, we examine the applicability of kinetic theory for dense granular flows, the structure and dynamics in sheared inelastic hard disks systems and dynamics of sheared non-spherical particles. Two complementary simulation techniques, the discrete element (DE) technique for soft particles and the event driven (ED) simulation technique for hard particles, are used to examine the extent to which the dynamics of an unconfined dense granular flow can be well described by a hard particle model when the particle stiffness becomes large. First, we examine the average co-ordination number for the particles in the flow down an inclined plane using the DE technique using both linear and Hertzian contact models. The simulations show that the average co-ordination number decreases below 1 for values of the spring stiffness corresponding to real materials such as sand and glass, even when the angle of inclination is only 1olarger than the angle of repose. The results of the two simulation techniques for the Bagnold coefficients (ratio of stress and square of the strain rate) and the granular temperature (mean square of the fluctuating velocity) are found to be in quantitative agreement. In addition, we also conduct the comparison of the pre-collisional relative velocities of particles in contact. Since momentum is transported primarily by particle contacts in a dense flow, the relative velocity distribution is a sensitive comparison of the dynamics in the two simulation techniques. It is found that the relative velocity distribution in both simulation techniques are well approximated by an exponential distribution for small coefficients of restitution, indicating that the dynamics of a dense granular flow can be adequately described by a hard particle model. The structure and dynamics of the two-dimensional linear shear flow of inelastic disks at high area fractions are analysed. The event-driven simulation technique is used in the hard-particle limit, where the particles interact through instantaneous collisions. The structure (relative arrangement of particles) is analysed using the bond-orientational order parameter. It is found that the shear flow reduces the order in the system, and the order parameter in a shear flow is lower than that in a collection of elastic hard disks at equilibrium. The distribution of relative velocities between colliding particles is analysed. The relative velocity distribution undergoes a transition from a Gaussian distribution for nearly elastic particles, to an exponential distribution at low coefficients of restitution. However, the single-particle distribution function is close to a Gaussian in the dense limit, indicating that correlations between colliding particles have a strong influence on the relative velocity distribution. This results in a much lower dissipation rate than that predicted using the molecular chaos assumption, where the velocities of colliding particles are considered to be uncorrelated. The orientational ordering and dynamical properties of the shear flow of inelastic dumbbells in two dimensions are studied, as a first step towards examining the effect of shape on the properties of flowing granular materials. The dumbbells are smooth fused disks characterised by the ratio of the distance between centers (L) and the disk diameter (D), and the ratio (L/D)varies between 0 and 1 in our simulations. Area fractions studied are in the range 0.1 to 0.7, while coefficients of normal restitution from 0.99 to 0.6 are considered. The simulations are similar to the event driven simulations for circular disks, but the procedure for predicting collisions is much more complicated due to the non-circular shape of the particles and due to particle rotation. The average orientation is measured using an orientational order parameter S, which varies between 0 (for a perfectly disordered fluid) and 1 (for a fluid with the axis of all dumbbells in the same direction). It is found that there is a gradual increase in ordering as the area fraction is increased, as the aspect ratio is increased or as the coefficient of restitution is decreased, and the order parameter has a maximum value of about 0.5 for the highest area fraction and lowest coefficient of restitution considered here. However, there is no discontinuous nematic transition for all the parameters studied here. The axis of the dumbbells are preferentially oriented along the extensional axis (at an angle of 45ofrom the flow direction) at low area fraction, but the orientation is closer to the flow direction as the area fraction is increased. The orientation distribution is calculated, and it is found that the orientation distribution is well described by a function of the form P(θ) =(1/π)+ (2S/π)cos(2(θ−θp)), where θis the angle from the flow direction and θpis the principal orientation direction. The mean energy of the velocity fluctuations in the flow direction is found to be higher than that in the gradient direction and the rotational energy, though the difference decreases as the area fraction increases, due to the efficient collisional transfer of energy between the three directions. The distributions of the translational and rotational velocity are found to be Gaussian distributions to a very good approximation. The equation of state for the pressure is calculated, and it is found to be remarkably independent of the coefficient of restitution. The pressure and dissipation rate show relatively little variation when scaled by the collision frequency for all the area fractions studied here, indicating that the collision frequency determines the momentum transport and energy dissipation even at the lowest area fractions studied here. The mean angular velocity of the particles is examined in some detail. It is found that the mean angular velocity is equal to half the vorticity at low area fractions, but the magnitude of the mean angular velocity systematically decreases to less than half the vorticity as the area fraction is increased, even though the stress tensor is symmetric.

Dans le présent travail, nous étudions l'impact sur un obstacle rigide d'un écoulement transitoire à surface libre de fluide viscoplastique. Cette étude est conduite numériquement à l'aide de la méthode SPH (Smoothed Particle Hydrodynamics), en y intégrant le modèle rhéologique de Herschel-Bulkley. Le code employé est adapté à nos besoins et validé sur des cas test classiques. Les caractéristiques locales de l'écoulement à proximité de l'obstacle sont analysées et deux régimes d'impact sont mis en évidence en fonction de la pente d'écoulement. L'étude des pressions exercées sur l'obstacle, conduite spatialement et temporellement en fonction de ces régimes d'impact, nous permet de mettre en évidence les rôles respectifs des composantes gravitationnelle et cinétique de la pression. Nos résultats sont comparés systématiquement à des résultats expérimentaux issus de travaux précédents et sont cohérents avec ces derniers. Une étude comparative de nos écoulements de fluide viscoplastique avec des écoulements de matériau granulaires de propriétés similaires nous conduit à mettre en évidence des caractéristiques communes entre les deux matériaux. / In this work, we study the impact of a transient free-surface flow of viscoplastic fluid on a rigid obstacle. This study is conducted numerically using the SPH (Smoothed Particle Hydrodynamics) method, and the Herschel-Bulkley rheological model. The SPH code is adapted to our needs and validated on classic benchmarks. The local characteristics of the flow near the obstacle are analysed and two impact regimes are highlighted depending on the slope angle. By studying of the pressure exerted on the obstacle, both spatially and temporally, with regards to these impact regimes, we evidence the respective roles of the gravitational and kinetic components of the pressure. Our results are systematically compared with experimental data from a previous work and are shown to be consistent. A comparative study conducted on both our viscoplastic flows and flows of granular material of similar properties highlights common characteristics of the two materials.

Dans le présent travail, nous étudions l'impact sur un obstacle rigide d'un écoulement transitoire à surface libre de fluide viscoplastique. Cette étude est conduite numériquement à l'aide de la méthode SPH (Smoothed Particle Hydrodynamics), en y intégrant le modèle rhéologique de Herschel-Bulkley. Le code employé est adapté à nos besoins et validé sur des cas test classiques. Les caractéristiques locales de l'écoulement à proximité de l'obstacle sont analysées et deux régimes d'impact sont mis en évidence en fonction de la pente d'écoulement. L'étude des pressions exercées sur l'obstacle, conduite spatialement et temporellement en fonction de ces régimes d'impact, nous permet de mettre en évidence les rôles respectifs des composantes gravitationnelle et cinétique de la pression. Nos résultats sont comparés systématiquement à des résultats expérimentaux issus de travaux précédents et sont cohérents avec ces derniers. Une étude comparative de nos écoulements de fluide viscoplastique avec des écoulements de matériau granulaires de propriétés similaires nous conduit à mettre en évidence des caractéristiques communes entre les deux matériaux. / In this work, we study the impact of a transient free-surface flow of viscoplastic fluid on a rigid obstacle. This study is conducted numerically using the SPH (Smoothed Particle Hydrodynamics) method, and the Herschel-Bulkley rheological model. The SPH code is adapted to our needs and validated on classic benchmarks. The local characteristics of the flow near the obstacle are analysed and two impact regimes are highlighted depending on the slope angle. By studying of the pressure exerted on the obstacle, both spatially and temporally, with regards to these impact regimes, we evidence the respective roles of the gravitational and kinetic components of the pressure. Our results are systematically compared with experimental data from a previous work and are shown to be consistent. A comparative study conducted on both our viscoplastic flows and flows of granular material of similar properties highlights common characteristics of the two materials.

Ce travail est financé par le fond européen Interreg VA (projet "PowderReg"). L'optimisation du transport, du stockage et du mélange des dispersions granulaires passe par le contrôle de leur écoulement, par exemple en ajoutant des vibrations mécaniques. Ce travail permet d'apporter une meilleure compréhension de l'influence des vibrations sur la rhéologie apparente de dispersions granulaires modèles dans le régime quasi-statique. Des travaux expérimentaux et numériques sont réalisés afin de sonder les modifications de la dynamique locale ainsi que les hétérogénéités qui apparaissent lors de leur mise en écoulement. Deux géométries présentant des intérêts fondamentaux, géophysiques et industriels sont étudiées : la cellule de type Couette et le plan incliné. Dans un premier temps, des mesures iso-indice couplées à de la fluorescence induite par laser ont été réalisées sur des suspensions granulaires dans une cellule de type Couette sous vibrations. Les vibrations rendent la rhéologie locale en homogénéisant le système. Elles suppriment la contrainte seuil apparente et font apparaître un plateau Newtonien à bas gradient, intrinsèque à la dynamique locale qui est de nature diffusive. Il est apparu que le temps de réarrangement des particules, dépendant de l'intensité des vibrations, est relié au volume libre disponible autour de chaque particule. L'étude numérique de cette géométrie dans le cas d'une dispersion granulaire sèche a mis en évidence des résultats similaires. Dans un second temps, des écoulements de dispersions de grains secs sur un plan incliné vibrant ont été réalisés. Des travaux numériques sur cette même géométrie ont permis d'enrichir cette étude. Les résultats mettent en évidence l'existence de deux régimes sous vibrations. Le comportement dans le régime dominé par la gravité n'est pas influencé par les vibrations, et un profil de Bagnold est observé. Les vibrations permettent principalement de baisser la friction basale, influant alors la hauteur des dépôts. Dans le régime dominé par les vibrations, les écoulements sont déclenchés par les vibrations elles-mêmes. Il est apparu qu'elles induisent des fluctuations de vitesses, qui créées alors une température granulaire. Cette température permet d'activer les réorganisations à l'échelle de la particule, éliminant le seuil apparent responsable du blocage des écoulements et permettant leur contrôle au travers du taux de cisaillement / This work is funded by the European Founds Interreg VA ("PowderReg" project). The optimization of the transport, storage and mixture of granular dispersions involves their flow control, by adding mechanical vibrations, for example. The present work enables a better understanding of the influence of the vibrations on the apparent rheology of model granular dispersions in the quasi-static regime. Experimental and numerical studies were carried out in order to probe the modifications of the local dynamic and heterogeneities that appear during the flow. Two geometries, interesting for diverse applications in fundamental science, geophysics and industries, were studied: the Couette cell geometry and the inclined plane. In the first step, refractive-index matching technique, coupled to the planar laser induced fluorescence was used to make measurements in granular suspensions in a vibrated Couette cell geometry. Vibrations make the rheology local by homogenizing the system. They suppress the yield stress and result in the appearance of a Newtonian plateau at the low shear, intrinsic to the local dynamic, which is diffusive in nature. It appears that the rearrangement time of the particles, depending on the intensity of the vibrations, is linked to a free volume available around each particle. The numerical study of this geometry in the case of a dry granular dispersion shows similar results. In the second step, flows of dry granular dispersions down the inclined and vibrated plane were realized. This study was completed with a numerical work. The results demonstrated the existence of two distinct regimes under vibrations. The behavior in the gravity-driven regime is not affected by the vibrations and a Bagnold profile is observed. The vibrations mainly cause the decrease of the basal friction and thus influence the height of the deposits. In the vibration-driven regime, however, flows are triggered by the vibrations themselves. It appears that they induce velocity fluctuations that create a granular temperature. That temperature allows the activation of the reorganizations at the grain scale. This suppresses the apparent yield responsible for the flow jamming, and thus enables their control through the shear rate

Geophysical research has long been interdisciplinary, with many phenomena on the Earth's surface involving multiple, linked processes that are best understood using a combination of techniques. This is particularly true in the case of grain flows on sand dunes, in which the sedimentary stratification with which geologists are concerned arises from the granular processes investigated by physicists and engineers, and the water permeation that interests hydrologists and soil scientists determines the seismic velocities of concern to exploration geophysicists. In this dissertation, I describe four projects conducted for the degree of Doctor of Philosophy, using a combination of laboratory experimentation, fieldwork, numerical simulation, and mathematical modelling to link avalanching on dunes to its effects on stratification, on the permeation of water, and on seismic surveys. Firstly, I describe experiments on erodible, unbounded, grain piles in a channel, slowly supplied with additional grains, and I demonstrate that the behaviour of the consequent, discrete avalanches alternates between two regimes, typified by their size statistics. Reconciling the `self-organised criticality' that several authors have predicted for such a system with the hysteretic behaviour that others have observed, the system exhibits quasi-periodic, system-spanning avalanches in one regime, while in the other avalanches pass at irregular intervals and have a power-law size distribution. Secondly, I link this power-law size distribution to the strata emplaced by avalanches on bounded grain piles. A low inflow rate of grains into an experimental channel develops a pile, composed of strata in which blue-dyed, coarser grains overlie finer grains. Associating stopped avalanche fronts with the `trapped kinks' described by previous authors, I show that, in sufficiently large grain piles, mean stratum width increases linearly with distance downslope. This implies the possibility of interpreting paleodune height from the strata of aeolian sandstones, and makes predictions for the structure of avalanche-associated strata within active dunes. Thirdly, I discuss investigations of these strata within active, Qatari barchan dunes, using dye-infiltration to image strata in the field and extracting samples across individual strata with sub-centimetre resolution. Downslope increases in mean stratum width are evident, while measurements of particle size distributions demonstrate preferential permeation of water along substrata composed of finer particles, explaining the strata-associated, localised regions of high water content discovered by other work on the same dunes. Finally, I consider the effect of these within-dune variations in water content on seismic surveys for oil and gas. Having used high performance computing to simulate elastic wave propagation in the vicinity of an isolated, barchan sand dune, I demonstrate that such a dune acts as a resonator, absorbing energy from Rayleigh waves and reemitting it over an extensive period of time. I derive and validate a mathematical framework that uses bulk properties of the dune to predict quantitative properties of the emitted waves, and I demonstrate the importance of internal variations in seismic velocity, resulting from variations in water content.

Les études présentées dans ce mémoire portent sur la simulation numérique et l'analyse physique des milieux granulaires immergés dans un fluide. Des développements numériques ont été réalisés pour coupler la méthode Lattice Boltzmann pour la dynamique du fluide avec la méthode de Contact Dynamics en 2D et avec la méthode Molecular Dynamics en 3D pour la dynamique des grains. Ces outils numériques ont été utilisés pour étudier l'initiation des avalanches sur un plan incliné en fonction de la compacité initiale et de l'angle d'inclinaison en 3D. Les résultats sont en bon accord quantitatif avec les expériences et ont permis de mettre en évidence la stabilisation de la pente granulaire par une pression négative du fluide interstitielle induite par la dilatance, et l'évolution spatiotemporelle des grandeurs telles que la compacité et la déformation de cisaillement. Ces évolutions dans la phase de fluage qui précède la rupture de pente ont pu être mises à l'échelle par un modèle théorique incorporant la loi de Darcy et l'effet de la dilatance sur l'angle de frottement interne. L'analyse de la texture granulaire a révélé la distorsion du réseau des contacts pendant le fluage et la saturation de l'anisotropie comme un critère de rupture. La propagation des avalanches granulaires a été étudiée dans la configuration 2D pour deux géométries différentes : 1) l'effondrement et l'étalement d'une colonne sous son propre poids, 2) l'étalement d'une pente sous l'effet d'une énergie cinétique injectée. Nous avons en particulier montré que la distance et la durée d'étalement obéissent à des lois de puissances en fonction du rapport d'aspect initial ou de l'énergie injectée. Le fluide exerce deux effets contradictoires : réduire les temps de relaxation et lubrifier les contacts. Ces effets ont été analysés dans le régime visqueux en fonction des conditions initiales et la viscosité du fluide. / By means of numerical simulations, we investigate the behavior of granular materials immersed in a fluid. Numerical developments were carried out to interface the Lattice Boltzmann Method for fluid dynamics with the Contact Dynamics method in 2D and with the Molecular Dynamics method in 3D for the simulation of the grains. Extensive simulations were applied to study the initiation of avalanches in a granular bed inclined at an angle above its angle of repose as function of the initial packing fraction and slope angle in 3D. The results are in excellent quantitative agreement with the reported experimental data, and show the stabilization of the granular bed by a negative pore overpressure induced by the dilatancy of the bed and the spatiotemporal evolution of the packing fraction and the shear deformation. The time evolution of these variables during the creeping phase before slope failure is scaled by a theoretical model accounting for darcian drag forces and the effect of dilatancy on the internal friction coefficient. We also analyzed the granular microstructure, which shows a gradual distortion of the contact network during creep at nearly a constant connectivity and the saturation of the anisotropy at failure. The runout of granular avalanches were investigated in 2D for two different configurations : 1) the collapse of a granular column under its own weight and 2) the runout of a granular pile as a result of kinetic energy supplied directly to the grains. We find power-law dependence of the resulting runout lengths and times with respect to the initial geometry or energy of the system. The time scales are shown to be consequence of two competing effects of the fluid on the grains : reducing relaxation times by viscous friction and the lubricating the contacts between grains.

Bien que répandus dans l'environnement et dans l'industrie, les écoulements granulaires possèdent des caractéristiques particulières qui sont différentes de celles des écoulements de liquides ordinaires comme l'eau (fluides newtoniens). L'une de ces caractéristiques est de présenter un seuil d'écoulement, il est donc fréquent qu'une des frontières de l'écoulement soit constituée de grains au repos (frontière érodable). L'échange entre les deux états solide et liquide d'un écoulement granulaire est à l'origine du phénomène d'érosion et de déposition et constitue le cœur de cette thèse. Nous y effectuons à la fois des études expérimentales et des études théoriques en nous basant sur les modèles phénoménologiques concernant le taux d'échange proposés par Bouchaud, Cates, Ravi Prakash et Edwards et modifiés par Boutreux et de Gennes. A l'aide d'un dispositif expérimental permettant de lâcher sans vitesse initiale une colonne de grains sur une pente granulaire, nous avons quantifié l'importance relative des phénomènes d'érosion et de déposition notamment en déterminant l'angle neutre qui définit la pente critique pour laquelle ces deux processus s'équilibrent. Dans un deuxième temps nous avons appliqué le modèle proposé par Boutreux et de Gennes au cas de l'étalement d'une « marche » granulaire (marche du Sinaï). / Granular flows, which are common in nature and industry, have particular characteristics that are different from the ordinary flow of liquid (eg. water, oil, etc.) and are not yet well understood in the literature. This case underlines the lack of knowledge on the exchange rate between solid and liquid states of granular flows which is the origin of the erosion - deposition phenomenon and constitutes the core of this thesis. To address the issue of the exchange rate solid - liquid, this work is based mainly in two subjects. The first subject is to study the processes of erosion and deposition of a flow generated by the grains stocked in a reservoir. By releasing naturally with zero-velocity the grains to generate the granular flow, we quantified the relation between the erosion and deposition for determining the neutral angle which defines the slope where these two processes balance. The second subject is to analyze the fall of grains in Sinai step problem by studying the spread of the march. For both works, we conduct not only the experimental but also theoretical studies based on the phenomenological models of the exchange rate proposed by BCRE (Bouchaud, Cates, Ravi Prakash and Edwards) and BDG (Boutreux and de Gennes) with some modifications.