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.

Abstract: The thesis work is part of a project that aims to develop a reliable CFD model to investigate the fluid-dynamics of a fluidized bubbling bed during gasification of refuse derived fuel (RDF) from sorted municipal solid waste (MSW). Gasification is a thermochemical process that converts carbon-containing materials into syngas. In this specific context scaling up is challenging because it implies dealing with a complex chemistry combined to heat and mass transfer phenomena in a multi-phase fluid environment. CFD modeling could represent a potential tool to predict the impact of the reactor configuration and operating conditions on gas yield, composition and potential contaminants. Validation of CFD simulations for such systems has been so far possible using different sophisticated experimental tools, allowing to link the model with experimental data. However, such high tech equipment may not always be available, especially at industrial scale. Hence, this work focuses on investigating the accuracy and numerical sensitivity of two different CFD models employed in the characterization of dense solid-particle flows in bubbling fluidized beds. The key parameter adopted to describe and quantify the dynamic behavior of this multiphase system is the power spectral density (PSD) distribution of pressure fluctuations. This PSD function was used to assess the accuracy of CFD models using one set of operating condition. The same type of analysis, extended to a wider range of operating conditions, may lead to a robust validation of the numerical models presented in this work. In spite of his measurement simplicity, pressure drop data present a strong connection with the bed fluid-dynamics and its interpretation could help to improve the fluidized bed technologies very fast, pushing CFD models closer to applications. / Résumé : Le but de ce projet est de développer un modèle CFD fiable pour étudier la dynamique des fluides d'un lit fluidisé en régime bullant pendant la gazéification de combustibles solides de récupération (CSR) triés à partir de déchets solides municipaux (DSM). La gazéification est un processus thermochimique qui convertit les matériaux contenant du carbone en gaz de synthèse. La mise à l'échelle est difficile dans ce cas car elle implique une chimie complexe combinée aux phénomènes de transfert de chaleur et de masse dans un environnement fluide multiphasique. La modélisation CFD représente un outil potentiel pour prédire l'impact de la configuration du réacteur et des conditions de fonctionnement sur le rendement, la composition et les contaminants potentiels du gaz. La validation des simulations CFD pour de tels systèmes a été jusqu'à présent possible grâce à l’utilisation de différents outils expérimentaux sophistiqués, permettant de lier le modèle aux données expérimentales. Toutefois, un tel équipement de pointe n’est pas toujours disponible, en particulier à l'échelle industrielle. Par conséquent, ce travail se concentre sur l'étude de la précision et de la sensibilité numérique de deux modèles CFD différents, utilisés dans la caractérisation des flux de particules solides denses dans les lits fluidisés bouillonnants. Le paramètre clé adopté pour décrire et quantifier le comportement dynamique de ce système multiphase est la distribution de la densité spectrale de puissance (DSP) des fluctuations de pression. La fonction DSP a été utilisée pour évaluer la précision des modèles CFD en utilisant un ensemble de conditions de fonctionnement. Le même type d'analyse, étendu à une plus large gamme de conditions de fonctionnement, peut conduire à une validation robuste des modèles numériques présentés dans ce travail. En dépit de sa simplicité de mesure, les données de chute de pression présentent une importante corrélation avec les lits fluidisés, de plus, leur interprétation pourrait aider à améliorer ces technologies très rapidement, poussant les modèles CFD plus près des applications.

<p>  </p> <p>Direct reduction ironmaking has gained popularity as a low carbon alternative to the typical blast furnace ironmaking route. A popular method of producing direct reduced iron is through the reduction of iron ore pellets in a reduction shaft furnace. Critical to this process is the use of a reliable continuous pellet feed system to provide a steady flow of pellets to the furnace. Therefore, any disruption in pellet flow can have a significant negative impact on the production rate of iron. </p> <p><br></p> <p>An iron ore pellet feed system for a direct reduction ironmaking furnace is jamming during winter operation. The pellets are jamming in a hopper at the top of the feed system above the furnace, and a hot gas, that seals off the furnace flue gas, flows counter to the pellets. A computational model of the feed system is built utilizing the discrete element method and computational fluid dynamics, using Siemen’s commercial multiphysics software Star-CCM+, to study the conditions that cause the jam to occur. The study is divided into six parts: pellet bulk flow calibration, computational cost reduction, modeling of the baseline operation, modeling the effect of moisture, development of a thermal model, and investigation of the minimal amount of icy and wet material to jam the system. The findings show that the location of jamming during operation matches the area in the simulation where it is most likely to occur, and that moisture alone is unlikely to result in jamming. Results indicate that the system will jam when charged with a minimum of 15% icy pellets, and when charged with 10% icy together with 5% wet pellets. Experimental work is recommended to validate the findings and to calibrate the simulations accordingly.</p>

Roterugnar är cylindriska kärl som används för att höja materials temperaturer i en kontinuerlig process som kallas för kalcinering. Roterugnar kan tillämpas i olika processer såsom reduktion av oxidmalm samt återvinning av farligt avfall. Fördelen med roterugnar ligger i dess förmåga att hantera råmaterial som sträcker sig från slam till granulära material med en mängd olika partikelstorlekar, och därigenom upprätthålla distinkta miljöer såsom en bädd av fasta partiklar som samexisterar med ett oxiderande fribord. Sex olika bäddbeteende har dokumenterats med avseende på fyllningsgrad samt Froude nummer. Syftet med denna studie var att utveckla en tvådimensionell suspensions modell med CFD genom att använda den kommersiella mjukvaran COMSOL 5.5 för att simulera de två faser, gas och fast, som en blandad fas efter verk av Philips et. al., Physics of Fluids A: Fluid Dynamics 4.1 (1992) 30-40 och Acrivos &amp; Zhang., International Journal Multiphase Flow 20.3 (1994) 579-591. Denna modell undersöktes genom att jämföra den med de dokumenterade flödesregimerna samt genom parameter som partikelstorlek, partikeldensitet och viskositeten hos gas i flödesregimen känd som rullande läge. Dessutom undersöktes temperaturprofilen för den roterande ugnen genom att utforska hur blandningsvariationer av den fasta bädden i den roterande ugnen påverkas av värmeöverföringen när värme tillförs från väggen under rullande läge. Resultaten av den tvådimensionella suspension modellen visade att det var bara möjligt att simulera glidläge korrekt; andra lägen kunde inte beskrivas som dokumenterat i litteraturen. Det indikeras att vilovinkeln och viskösa krafter i den roterande ugnen var låga vilket resulterade i att suspensions modellen inte kunde avbilda exakt de återstående flödesregimerna som dokumenterat. Till exempel avbildades rullningsläget mer likt forsandeläge då partiklarna fall fritt efter höjning av bädden. Partikelstorlek och partikeldensitet har visat sig ha en betydande påverkan på suspensions modellen eftersom de viskösa krafterna blir låga för en partikelstorlek och partikeldensitet under 0,4 mm respektive 1500 kg/m3. Angående gasens viskositet visades det sig att ju närmare värdet 2.055e-3 (Pa*s) den blev desto större blev sedimentationsflödet vilket resulterade i att bäddpartiklarna dras ner och förblir där. Suspensions modellen kunde således simulera en fast och flytande fas och inte en gasfas som avsett. Slutligen visade temperaturanalysen att påverkan av den termiska konduktiviteten var mer signifikant än den specifika värmekapaciteten i intervallet 1 - 50 (W/(m*K)) respektive 300 - 800 (J/(kg*K)) på grund av den tid det tog att nå en homogen temperaturprofil. / Rotary kilns are cylindrical vessels used to raise materials temperature in a continuous process known as calcination. Rotary kilns find application in various processes such as reduction of oxide ore and hazardous waste reclamation. The advantage of the rotary kiln lies in its ability to handle feedstock ranging from slurries to granular materials with a variety of particle size, thereby maintaining distinct environments such as a bed of solid particles coexisting with an oxidising freeboard. Six different bed behaviours within the kiln have been documented with respect to the filling degree and Froude number. The aim of this study was to develop a two-dimensional suspension model with CFD by using the commercial software COMSOL 5.5 to simulate the two phases, gas and solid, as a mixed phase, following the works of Philips et. al., Physics of Fluids A:  Fluid Dynamics 4.1 (1992) 30-40 and Acrivos &amp; Zhang., International Journal Multiphase Flow 20.3 (1994) 579-591. This model was investigated by comparing it against the documented flow regimes as well as through parameters such as particle size, particle density and viscosity of gas in the flow regime known as rolling mode. In addition, the temperature profile of the rotary kiln was investigated by exploring how the mixture variation of the solid bed within the rotary kiln affects the heat transfer when heat is supplied from the wall during a rolling mode. The results of the two-dimensional suspension model showed that it was only possible to simulate the slipping mode accurately; others mode could not be described as documented in literature. It is indicated that the angle of repose and viscous forces within the rotary kiln were low resulting in the suspension model not being able to accurately depict the remaining flow regimes as documented. For instance, the rolling mode was depicted more as a cataracting mode due to the free fall of particles after elevation of the bed. The particle size and the particle density were found to have a significant impact on the suspension model as the viscous forces became low for a particle size and particle density below 0.4 mm and 1500 kg/m3 respectively. As for the viscosity of gas it was found that the closer it got to the value 2.055e-3 (Pa*s) the sedimentation flux became too large resulting in the bed particles being pulled down and remaining there. Thus, the suspension model could simulate a solid and liquid phase and not a gas phase as intended. Lastly, the temperature analysis revealed that the impact of the thermal conductivity was more significant than the specific heat capacity in the range of 1 - 50 (W/(m*K)) and 300 - 800 (J/(kg*K)) respectively, due to the time it took to reach a homogeneous temperature profile.