Digital elevation models (DEMs) and their derivatives such as slope, flow direction and flow accumulation maps, are used frequently as inputs to hydrologic and nonpoint source modeling. The depressions which are frequently present in DEMs may represent the actual topography, but are often the result of errors. Creating a depression-free surface is commonly required prior to deriving flow direction, flow accumulation, flow network, and watershed boundary maps. The objectives of this study were: 1) characterize the occurrence of depressions in 30m USGS DEMs and assess correlations to watershed topographic characteristics, and 2) compare the performance of two algorithms used to remove depressions and delineate flow networks from DEMs. Sixty-six watersheds were selected to represent a range of topographic conditions characteristic of the Piedmont and Mountain and Valley regions of Virginia. Analysis was based on USGS 30m DEMs with elevations in integer meters. With few exceptions watersheds fell on single 7.5minute USGS quadrangle sheets, ranged in size from 450 to 3000 hectares, and had average slopes ranging from 3 to 20 percent. ArcView (3.1) with the Spatial Analyst (1.1) extension was used to summarize characteristics of each watershed including slope, elevation range, elevation standard deviation, curvature, channel slope, and drainage density. TOPAZ (ver 1.2) and ArcView were each used to generate a depression-free surface, flow network and watershed area. Characteristics of the areas 'cut' and 'filled' by the algorithms were compared to topographic characteristics of the watersheds. Blue line streams were digitized from scanned USGS 7.5minute topographic maps (DRGs) then rasterized at 30 m for analysis of distance from the derived flow networks. The removal of depressions resulted in changes in elevation values in 0 - 11% of the cells in the watersheds. The percentage of area changed was higher in flatter watersheds. Changed elevation cells resulted in changes in two to three times as many cells in derivative flow direction, flow accumulation and slope grids. Mean fill depth by watershed ranged from 0 to 10 m, with maximum fill depths up to 40 m. In comparison with ArcView, TOPAZ, on average affected 30% fewer cells with less change in elevation. The significance of the difference between ArcView and TOPAZ decreased as watershed slope increased. A spatial assessment of the modified elevation and slope cells showed that depressions in the DEMs occur predominantly on or along the flow network. Flow networks derived by ArcView and TOPAZ were not significantly different from blue line streams digitized from the USGS quadrangles as indicated by a paired t test. Watershed area delineated by ArcView and TOPAZ was different for almost all watersheds, but was generally within 1%. Conclusions from this study are: 1) The depressions in 30 m DEMs can make up a significant portion of the area especially for flatter watersheds; 2) The TOPAZ algorithm performed better than ArcView in minimizing the area modified in the process of creating a depressionless surface, particularly in flatter topography; 3) Areas affected by removing depressions are predominantly adjacent to the stream network; 4) For every elevation cell changed, slopes are changed for two to three cells, on average; and 5) ArcView and TOPAZ derived flow networks closely matched the blue line streams. / Master of Science

Landslides are one of the most widespread, yet difficult to predict, geohazards in the world. Theyrepresent a major threat to society, which is why it is so important to understand and study them.Landslide detection is one of the approaches that help scientists to gain a better knowledge of landslides,how they form and the characteristics they hold. In this work, a new landslide detection technique istested in the Klarälven Region, Sweden. The method involves using multi-temporal DEMs generatedfrom LiDAR point clouds from 2010 and 2020. Once two respective DEMs are generated from theoriginal point cloud data sets, a DEM of Difference (DoD) is obtained by subtracting the 2010 DEMfrom the 2020 DEM. The DoD highlights areas of vertical change over the ten-year period. In order todetect potential landslides, a Level of detection (LoD) threshold is applied, as well as a second upperthreshold filtering out elevation difference resulting from anthropogenic activities. Lastly, a spatialthreshold is applied to eliminate any remaining potential errors clusters that are too small to belandslides. Once this stage is completed, the study area is surveyed, focusing on certain areacharacterised by particular slope angles and soil types most likely to exhibit mass movement. Thishelped in the analysis and interpretation of the data. In order to discuss the results as fully as possible,four types of region were studied in more detail using four representative examples. Since a few doubtfulcases were found, potential sources of error are also being explored. Despite minor inaccuracies, themethod appears to be a relevant technique for detecting differences in altitude, even if it also considersunnatural mass movements.

Die überarbeitete PhD-Dissertation (2017) des irischen Religionswissenschaftlers Brendan McNamara von der University of Cork befasst sich mit der religiösen Aufbruchstimmung im Großbritannien des Jahrzehnts vor dem Ausbruch des Ersten Weltkrieges.

This thesis presents a computational framework that models the coupled behavior between sand deposits and their surrounding fluid. Particle deposits that form in gas turbine engines and industrial burners, can change flow dynamics and heat transfer, leading to performance degradation and impacting durability. The proposed coupled framework allows insight into the coupled behavior of sand deposits at high temperatures with the flow, which has not been available previously. The coupling is done by using a CFD-DEM framework in which a physics based collision model is used to predict the post-collision state-of-the-sand-particle. The collision model is sensitive to temperature dependent material properties of sand. Particle deposition is determined by the particle's softening temperature and the calculated coefficient of restitution of the collision. The multiphase treatment facilitates conduction through the porous deposit and the coupling between the deposit and the fluid field. The coupled framework was first used to model the behavior of softened sand particles in a laminar impinging jet flow field. The temperature of the jet and the impact surface were varied(T^* = 1000 – 1600 K), to observe particle behavior under different temperature conditions. The Reynolds number(Rejet = 20, 75, 100) and particle Stokes numbers (Stp = 0.53, 0.85, 2.66, 3.19) were also varied to observe any effects the particles' responsiveness had on deposition and the flow field. The coupled framework was found to increase or decrease capture efficiency, when compared to an uncoupled simulation, by as much as 10% depending on the temperature field. Deposits that formed on the impact surface, using the coupled framework, altered the velocity field by as much as 130% but had a limited effect on the temperature field. Simulations were also done that looked at the formation of an equilibrium deposit when a cold jet impinged on a relatively hotter surface, under continuous particle injection. An equilibrium deposit was found to form as deposited particles created a heat barrier on the high temperature surface, limiting more particle deposition. However, due to the transient nature of the system, the deposit temperature increased once deposition was halted. Further particle injection was not performed, but it can be predicted that the formed deposit would begin to grow again. Additionally, a Large-Eddy Simulation (LES) simulation, with the inclusion of the Smagorinsky subgrid model, was performed to observe particle deposition in a turbulent flow field. Deposition of sand particles was observed as a turbulent jet (Re jet=23000,T_jet^*= 1200 K) impinged on a hotter surface(T_surf^*= 1600 K). Differences between the simulated flow field and relevant experiments were attributed to differing jet exit conditions and impact surface thermal conditions. The deposit was not substantive enough to have a significant effect on the flow field. With no difference in the flow field, no difference was found in the capture efficiency between the coupled and decoupled frameworks. / Master of Science / Particle deposits can form in a wide range of environments leading to altered performance. In applications, such as jet engines, particles are heated to critically high temperatures. At these high temperatures, the particles can soften, and begin to exhibit characteristics of both a liquid and a solid. Overtime as these softened particles aggregate on a wall, a deposit will begin to form. These deposits alter the geometry resulting in changes in fluid temperature and velocity. This change in fluid behavior will affect the rate of particle deposition that happens in the future. There has been limited work that has looked at the coupled behavior between a deposit and its surrounding fluid, experimentally or computationally. The purpose of this research was to develop a framework that models the deposition of softened particles, and the coupled behavior between deposits and the fluid. This research was able to show that the presence of a deposit could change its surrounding fluid’s velocity and temperature significantly. Differences in the rate of particle deposition also occurred when a deposit had formed on a surface. These results show the importance of capturing the relationship between deposits and the surrounding fluid. With further development, this proposed framework can provide insight into altered gas turbine performance and can lead to improved maintenance plans.

Recently observed cumulation of unexpected collapses of slope-distant waste dumps in lignite mining areas of eastern germany re-initiated research of soil liquefaction. Especially it turned the question of internal initials that correspond to water rise. Parallel to laboritory tests and field experiments a micromechanical model should be developed, which can reproduce processes in the soil during saturation. In first approximation a partly saturated soil consists of two phases: the soil particles and the pore fluid. For micromechanical modeling a coupling of discontinuum particles) and continuum (fluid) is required. The soil particles can be simulated with the Discrete-Element-Method (DEM). For the pore fluid, which is assumed to be a mixture of liquid and gaseous fractions, Pore scale model with Finite Volumes (PFV) is used. At low water content liquid bridges (meniscii) arise between the particles that cause an apparent cohesion. The effect of the meniscii is considered by a correspondingly contact law in the DEM model. During the saturation of a soil both, cohesive effect and fluid bulk modulus, are reduced. In addition buoyancy acts on the particles during the process. The micromechanical modeling approach has the advantage, that just a few model parameters are needed. The numerical model shows pore fluid pressures during saturation process, that leads to a reduction of effective stress. It is investigated how much the reduction is regarding porosity, degree of saturation, stress conditions and grain shape. Furthermore the influence of model parameters as well as hydromechanics is investigated. The investigations are completed with another series of experiments under special conditions like integration of macropores, horizontal fixed model boundaries and abrupt saturation.

Regional landslide risk, as it is most commonly defined, is a product of the following: hazard, vulnerability and exposed population. The first objective of this research project is to estimate the regional landslide hazard level by calculating its probability of slope failure based on maximum slope angles, as estimated using data provided by digital elevation models (DEM). Furthermore, it addresses the impact of DEM resolution on perceived slope angles, using local averaging theory, by comparing the results predicted from DEM datasets of differing resolutions. Although the likelihood that a landslide will occur can be predicted with a hazard assessment model, the extent of the damage inflicted upon a region is a function of vulnerability. This introduces the second objective of this research project: vulnerability assessment. The third and final objective concerns the impact of urbanization and population growth on landslide risk levels.

Les analyses structurale et mécanique des pentes rocheuses constituent des éléments clés pour l'évaluation de leur stabilité. L'utilisation complémentaire de la photogrammétrie et des modèles numériques qui couplent les réseaux discrets de discontinuités (DFN selon son sigle en anglais) avec la méthode des éléments discrets (DEM selon son sigle en anglais), présente une méthodologie qui peut être utilisée pour évaluer le comportement mécanique des configurations tridimensionnelles de terrain pour lesquelles l'existence de discontinuités non persistantes peut être supposée. La stabilité des masses rocheuses est généralement supposée être contrôlée par la résistance au cisaillement le long des plans de discontinuité. Si les discontinuités sont non persistantes, avec leur continuité interrompue par la présence de ponts rocheux (portions de roche intacte reliant la masse rocheuse au massif), leur résistance apparente augmente considérablement. Dans ce cas, la contribution des ponts rocheux localisés entre ces discontinuités doit être prise en compte dans l'analyse de stabilité. La déstabilisation progressive des massifs rocheux dans lesquels des discontinuités non persistantes sont présentes, peut être étudiée par des simulations numériques réalisées à l'aide de l'approche DEM. La roche intacte est représentée comme un assemblage de particules (ou éléments discrets) liées entre elles par des contacts dont les lois de comportement spécifiques peuvent être calibrées pour représenter correctement le comportement de la roche. L'intérêt de la méthode est qu'elle permet de simuler l'initiation de la rupture et sa propagation à l'intérieur de la matrice rocheuse du fait de la rupture des contacts cohésifs entre les particules. De plus, les discontinuités préexistantes peuvent être prises en compte explicitement dans le modèle en utilisant une loi de contact ad hoc qui assure un comportement mécanique représentatif des plans de discontinuité. Des analyses de stabilité ont été effectuées et ont mis en évidence le rôle des ponts rocheux dans la génération de nouvelles surfaces de rupture qui peuvent se développer à travers des mécanismes de rupture mixte en traction et en cisaillement. On peut considérer la formulation de Jennings comme l'une des premières méthodes d'analyse de la stabilité des pentes rocheuses qui évaluent la résistance au glissement comme une combinaison pondérée des résistances mécaniques des ponts rocheux et des plans de discontinuité. Sa validité a été discutée et systématiquement comparée aux résultats obtenus à partir de simulations numériques. Il a pu être montré que la formulation de Jennings perd sa validité dès que la rupture des ponts rocheux intervient majoritairement par des mécanismes de traction. Une formulation complémentaire a alors été proposée. En ce qui concerne l'étude de la stabilité des massifs rocheux sur site, il a été montré que l'association entre les données issues de la photogrammétrie en haute résolution et l'approche DFN-DEM peut être utilisée pour identifier des scénarios de rupture. L'analyse en retour de cas réels a montré que les surfaces de rupture peuvent être simulées comme le résultat de mécanismes combinant la fracturation des ponts rocheux et le glissement le long des discontinuités préexistantes. La rupture d'un dièdre qui a eu lieu dans une mine de charbon australienne, a été utilisée pour valider cette méthodologie. Des simulations numériques ont été réalisées pour déterminer les scénarios pour lesquels les surfaces de rupture simulées et celles repérées sur le terrain, peuvent être utilisés pour calibrer les paramètres de résistance du modèle numérique. Le travail présenté ici répond à un besoin plus général visant à améliorer la gestion des risques naturels et miniers liés aux masses rocheuses instables. La méthodologie proposée constitue une alternative robuste dédiée à renforcer la fiabilité des analyses de stabilité pour les pentes rocheuses fracturées à structure complexe. / Structural and mechanical analyses of rock mass are key components for rock slope stability assessment. The complementary use of photogrammetric techniques and numerical models coupling discrete fracture networks (DFN) with the discrete element method (DEM) provides a methodology that can be applied to assess the mechanical behaviour of realistic three-dimensional (3D) configurations for which fracture persistence cannot be assumed. The stability of the rock mass is generally assumed to be controlled by the shear strength along discontinuity planes present within the slope. If the discontinuities are non–persistent with their continuity being interrupted by the presence of intact rock bridges, their apparent strength increases considerably. In this case, the contribution of the rock bridges located in-between these discontinuities have to be accounted for in the stability analysis. The progressive failure of rock slope involving non–persistent discontinuities can be numerically investigated based upon simulations performed using a DEM approach. The intact material is represented as an assembly of bonded particles interacting through dedicated contact laws that can be calibrated to properly represent the behaviour of the rock material. The advantage of the method is that it enables to simulate fracture initiation and propagation inside the rock matrix as a result of inter-particle bond breakage. In addition, pre–existing discontinuities can be explicitly included in the model by using a modified contact logic that ensures an explicit and constitutive mechanical behaviour of the discontinuity planes. Stability analyses were carried out with emphasis on the contribution of rock bridges failure through a mixed shear-tensile failure process, leading to the generation of new failure surfaces. Jennings' formulation being considered to be one of the first rock slope stability analysis that evaluates the resistance to sliding as a weighted combination of both, intact rock bridges and discontinuity planes strengths, its validity was discussed and systematically compared to results obtained from numerical simulations. We demonstrate that the validity of Jennings' formulation is limited as soon as tensile failure becomes predominant and an alternative formulation is proposed to assess the resulting equivalent strength. Regarding field slope stability, we show that the combination of high resolution photogrammetric data and DFN-DEM modelling can be used to identify valid model scenarios of unstable wedges and blocks daylighting at the surface of both natural and engineered rock slopes. Back analysis of a real case study confirmed that failure surfaces can be simulated as a result of both fracture propagation across rock bridges and sliding along pre-existing discontinuities. An identified wedge failure that occurred in an Australian coal mine was used to validate the methodology. Numerical simulations were undertaken to determine in what scenarios the measured and predicted failure surfaces can be used to calibrate strength parameters in the model. The work presented here is part of a more global need to improve natural and mining hazards management related to unstable rock masses. We believe that the proposed methodology can strengthen the basis for a more comprehensive stability analysis of complex fractured rock slopes.

La modélisation multi-échelle multi-physique de la solidification d'alliages métalliques demande de combiner des phénomènes à l'échelle macroscopique du produit et microscopiques à l'échelle des structures de solidification. Dans cette thèse, l'empilement aléatoire des grains équiaxes avec des morphologies typiques de solidification est étudié. Nous mettons tout d'abord en évidence les paramètres hydrodynamiques adimensionnels qui régissent l'empilement de grains équiaxes : le nombre de Stokes, St, le nombre d'Archimède, Ar, et le rapport entre le temps caractéristique de la croissance et le temps caractéristique du mouvement, Γ. Un dispositif expérimental a été conçu par similitude hydrodynamique avec le phénomène réel de l'empilement de la solidification afin d'étudier l'influence de la géométrie des grains équiaxes et l'influence des conditions hydrodynamiques sur la fraction d'empilement. En outre, un outil numérique basé sur le méthode des éléments discrets a été développé pour compléter le travail expérimental de détermination de : la fraction d'empilement locale, le nombre de particules voisines en contact et l'orientation des particules. Des fractions d'empilement entre environ 0,53 et 0,67 ont été mesurées et calculées pour les grains sphériques non-cohésifs, alors que des valeurs allant jusqu'à environ 0,30 sont trouvées pour les grains dendritiques non-cohésifs. Enfin, nous étudions la dynamique de l'empilement, qui est la transition d'un régime de sédimentation à l'équilibre mécanique. L'évolution des variables comme la fraction locale de solide, le nombre de particules voisines en contact et l'orientation du grain en fonction du temps est présentée / Solidification multiphase multiscale modeling of metal alloys is based on the combination of the phenomena at the macroscopic scale of the product and at the microscopic scale of the solidification structures. In this thesis, the random packing of the typical equiaxed grain morphologies in metal alloy solidification is investigated. Firstly, we highlight the hydrodynamic dimensionless parameters governing the grain packing in the melt: the Stokes number, St, the Archimedes number, Ar, and the growth-to-motion ratio, Γ. Subsequently, an experimental setup is designed by hydrodynamic similarity with the actual solidification packing phenomenon in order to investigate the influence of the equiaxed grain geometry and the hydrodynamic conditions on the average solid packing fraction. Additionally, a numerical Discrete Element Method tool is developed to complement the experimental work by accessing to those granular variables which result difficult to be experimentally obtained such as the local packing fraction, the contacting neighbors and the particle orientation. Packing fractions between approximately 0.53 and 0.67 are measured and computed for the spherical noncohesive grains, for different hydrodynamic, frictional and polydispersity conditions, whereas values down to approximately 0.30 are found for noncohesive dendrite envelopes. Finally, we investigate the packing dynamics, which is the transition from a sedimentation regime to the mechanical equilibrium (packing). The evolution of the local solid fraction, contacting neighbors, mechanical contacts and grain orientation are given

Dans la pratique actuelle du génie géotechnique, les essais de pénétration tels que les CPT, SPT, Panda sont largement utilisés pour caractériser mécaniquement les sols, au travers notamment d’une caractéristique de rupture appelée résistance de pointe. Par ailleurs, les dernières évolutions technologiques apportées aux essais de pénétration dynamique (Panda 3) permettent d’obtenir pendant chaque impact une courbe charge–enfoncement donnant la charge en pointe en fonction de l’enfoncement à partir de la mesure et du découplage des ondes générées durant l’essai. L’exploitation de cette nouvelle courbe fournit des informations non seulement sur la résistance de pointe dynamique mais également sur des paramètres mécaniques complémentaires mis en jeu pendant l’enfoncement de la pointe. L’objectif de cette thèse est de développer un modèle numérique en 2D capable de reproduire les signaux pénétrométriques obtenus expérimentalement par essais de type statique ou dynamique. Ce modèle est basé sur la méthode des éléments discrets avec une loi de contact linéaire simple. Une fois le modèle validé, une étude paramétrique a été réalisée en jouant essentiellement sur les modes d’application de la sollicitation (vitesse d’impact ou de pénétration), la granulométrie du matériau ainsi que l’arrangement granulaire (variation de la densité). Outre l’influence de ces paramètres sur les signaux pénétrométriques et la résistance de pointe mesurée, une attention particulière est portée sur l’analyse micromécanique : dissipation d’énergie dans le milieu, évolution des chaines de force, orientations des contacts. Cette analyse nécessite de développer des outils numériques spécifiques afin de mieux comprendre le mécanisme de l’enfoncement et tenter d’expliquer la réponse mécanique macroscopique obtenue. L’effet de la vitesse n’influence significativement que sur les essais de pénétration statiques et dynamiques en régime d’écoulement dense. A vitesse d’enfoncement comparable, il n’y a aucune différence significative au niveau microscopique entre les deux modes de sollicitation statique et dynamique. En ce qui concerne l’influence des caractéristiques du matériau, les résultats obtenus par le modèle numérique conforment aux celui réel lors que le frottement entre particules ou la compacité du milieu varie. Concernant la granulométrie, la variation de la courbe charge-enfoncement et la force de pointe dynamique augmente lorsque le diamètre moyen augmente. / In the field of in situ mechanical characterization of soils, penetration tests are commonly used. Penetration tests measure the properties of soils in the domain of large deformations. The tip resistances, deduced from pile driving theory, can be measured either in dynamic conditions (q d ) either in static conditions (q c ). Recently, the measurement technique in dynamic conditions has been improved and it is now possible to record the whole response of the soil during one impact in terms of tip force and penetration distance. The exploitation of this new curve provides information not only on dynamic tip resistance but also on additional mechanical parameters involved during the driving of the tip. The objective of this work is to develop a numerical model in 2D able to reproduce the penetrometric record obtained experimentally by static or dynamic penetration tests. This model is based on the discrete element method with a simple linear contact model. After the validation of the model, a parametric study was performed essentially on the loading type (static or dynamic), the penetration rate, the particle size of the granular material and the arrangement (density variation). Besides the influence of these parameters on the penetrometer signals and the tip resistance, a particular attention was focused on micromechanical analysis: energy dissipation in the medium, force chain evolution, contact orientation. This analysis requires the development of specific numerical tools to better understand the penetration mechanism and try to explain the macroscopic mechanical response obtained. The penetration rate influences significantly only in the dense flow regime on the static and dynamic penetration tests. There is no significant microscopic difference between static and dynamic penetration tests with similar penetration rates. Regarding the influence of the characteristics of the material, the numerical results obtained conform to the real results when the particle friction or the compactness of the medium varies. Concerning the particle size, the dynamic signal variation and the dynamic tip force increases when the average particle diameter increases.

Le travail réalisé en collaboration entre l’entreprise CNIM et le laboratoire IUSTI a permis la miseen place d’un outil de simulation numérique afin d’étudier à l’échelle locale différents phénomènescouplés qui se produisent dans un réacteur de gazéification en lit fluidisé dense. L’approche choisie,DEM (Discrete Element Method), est basée sur le suivi de paquets de particules. Les résultats desimulation ont été comparés à des mesures expérimentales réalisées dans les laboratoires IUSTIet LERMAB : mesure de pression dans un lit fluidisé peu profond et caractérisation de la réactionde pyrolyse par suivis temporels de la masse volumique par une méthode innovante et de latempérature de la particule de biomasse ainsi que la composition des gaz produits. À plus grandeéchelle, les prédictions numériques du code ont été comparées à des mesures expérimentales deprofil de vitesse des particules et de mélange et ségrégation issues de la littérature. Les résultatsnumériques du lit fluidisé bisolide sans et avec réactions chimiques, obtenus sont en bon accordavec les mesures expérimentales. Ils ont permis de mettre en évidence le fort couplage entre lesphénomènes hydrodynamiques et thermochimiques.Ces travaux ouvrent de nouvelles perspectives tant sur le plan de l’expérimentation que de la simulation numérique où l’approche DEM retenue a montré un fort potentiel. L’extrapolation dumodèle DEM pour la simulation de la gazéification à des réacteurs industriels reste encore un défiau regard des ressources informatiques. / Energy production from green and renewable resources, such as biomass, are currently experiencinga significant growth. Thermochemical conversion of this biomass by gasification is a process usedfor over a century but still requires significant developments in terms of rentability optimizationand quality improvement of products gases.The work carried out in collaboration between the CNIM company and the IUSTI laboratoryallowed the establishment of a numerical simulation tool to study locally different coupled phenomenaoccurring in a dense fluidized bed gasification reactor. The chosen approach, DEM (DiscreteElement Method), is based on the monitoring of particle packets. The simulation results werecompared to experimental measurements realised in IUSTI and LERMAB laboratories : measuringpressure in a shallow fluidized bed and characterization of pyrolysis reaction with the timetracking of particle density by an innovative method, temperature and the composition of theproducts gases. On a larger scale, the numerical code predictions were compared with velocity,particles mixing and segregation profiles from experimental measurements of the literature. Theobtained numerical results of bi-solid fluidized bed with and without chemical reactions are ingood agreement with the experimental measurements. It helped to highlight the strong couplingbetween hydrodynamic and thermochemical phenomena.This work opens up new perspectives on the experimental plan and numerical simulation whichDEM approach has shown great potential. The extrapolation of the DEM model for the simulationof gasification industrial reactors remains a challenge in terms of computer resources.