Uma climatologia de 21 anos dos Distúrbios Ondulatórios de Leste (DOLs) sobre a região NEB foi realizada com o intuito de obter um melhor entendimento dos processos dinâmicos e sinóticos do ciclo de vida destes sistemas, incluindo gênese, crescimento/decaimento, trajetória e dissipação. Adicionalmente, foi avaliada a eficiência do modelo de mesoescala WRF em simular este tipo de sistema. A identificação dos DOLs foi obtida de forma subjetiva através de imagens de satélite no canal infravermelho e campos de linhas de corrente e vorticidade relativa nos níveis de 1000, 850, 700, 500 e 200 hPa da reanálise do ERA-Interim. Neste período foram identificados 518 eventos de DOLs, onde 97% (3%) dessas ondas atingiram (não atingiram) a região do NEB, 64% (36%) foram convectivas (não convectivas) e 14% (86%) atingiram a região da Amazônia. Os principais sistemas que deram origem aos DOLs foram: Zona de Convergência Intertropical (ZCIT), Vórtice Ciclônico de Altos Níveis (VCAN), Frentes Frias (FF) e conglomerados convectivos provenientes da costa oeste da África (AF). Além destes, a interação entre os sistemas AF/FF, AF/ZCIT, AF/VCAN e ZCIT/FF estiveram associadas a sua formação. Em média, foram observadas aproximadamente 25 ondas por ano, com máxima (mínima) frequência compreendida entre os meses de Março a Agosto (Setembro a Fevereiro) e, com pronunciada variabilidade interanual. O ciclo de vida dos DOLs foi avaliado de forma objetiva a partir de um método automático de identificação e rastreio (TracKH), onde dos 518 eventos identificados na análise subjetiva, conseguiu capturar 342 ondas ( ~ 66%). A partir desta detecção, uma análise foi feita no nível 850 hPa para determinar as características típicas do seu ciclo de vida, que mostra um pico proeminente entre as longitudes de 35°W e 15°W e latitudes de 20°S e 5°N associado a densidade de gênese. A trajetória e dissipação se concentraram sobre a costa leste do NEB, entre os estados de Alagoas e Rio Grande do Norte, porém a dissipação diminuía a medida que adentrava ao continente. Os padrões sinóticos associados aos DOLs foram analisados através da anomalia de composição durante o período de máxima (úmido) e mínima (seco) frequência desde 3 dias antes até um dia após os DOLs atingirem a costa do NEB. Durante o período úmido, a circulação associada aos DOLs apresentou anomalia ciclônica e confluente, vorticidade negativa e convergência até médios níveis, enquanto que em 200 hPa apresentou apenas a característica do cavado. Por outro lado, para o período seco, estas características foram observadas somente em baixos níveis. Anomalias negativas de movimentos verticais e temperatura e positivas de umidade foram observadas associadas aos DOLs nas duas estações, porém atingindo maiores altitudes durante o período seco. A composição de precipitação indicou que os DOLs são responsáveis por um acréscimo de 16% (4%) durante a estação úmida (seca). As características típicas dessas ondas foram: período médio de 8 (73) dias, comprimento de onda de aproximadamente 4500 (5500) km e velocidade de fase da ordem de 6,5 (0,9) m.s-1, para o período úmido (seco). O modelo WRF simulou os padrões sinóticos, a precipitação e as características típicas associadas aos DOLs de forma coerente, sendo comparado aos resultados do Climate Forecast System Reanalysis (CFSR), porém com maior intensidade, para os eventos de 11 e 17 de junho de 2006. / A 21-year climatology of Easterly Wave Disturbances (EWDs) over NEB region was constructed in order to obtain a better understanding of dynamic and synoptic processes life cycle of these systems, including genesis, growth / decay, trajectory and dissipation. Moreover, we evaluated the efficiency of WRF mesoscale model to simulate this type of system. The identification of EWDs was obtained subjectively through satellite images in infrared channel and fields of streamlines and relative vorticity at the levels 1000, 850, 700, 500 and 200 hPa from ERA-Interim reanalysis. During this period, 518 EWDs were identified, where 97% (3%) of these waves hit (not hit) the NEB region, 64% (36%) were convective (non-convective) and 14% (86%) reached the Amazon region. The main systems that gave rise to DOLs were: Intertropical Convergence Zone (ITCZ), Upper-Tropospheric Cyclonic Vortices (UTCV), Cold Fronts (FF) and convective clusters from the west coast of Africa (AF). In addition, the interaction between systems AF/FF, AF/ITCZ, AF/UTCV and ITCZ/FF were associated with their formation. On average, we observed approximately 25 waves per year, with maximum (minimum) frequency between the months of March to August (September to February) and, with pronounced interannual variability. The life cycle of EWDs was evaluated objectively using an automatic method for the identification and tracking (TracKH). From the 518 events identified in the subjective analysis, TracKH was able to capture 342 waves (~ 66%). From this detection, an analysis was made at 850 hPa level to determine the characteristics of their life cycle, which shows a prominent peak around of 35°W-15°W and 20°S-5°N associated with density genesis. The trajectory and dissipation have concentrated over east coast of the NEB, between the states of Alagoas and Rio Grande do Norte, but the dissipation decreased once the systems enter the continent. The synoptic patterns associated with EWDs were analyzed by composing anomaly during the period of maximum (wet) and minimum (dry) frequency from 3 days before until one day after the EWDs reaching the NEB coast. During the wet period, the circulation presented cyclonic and confluent anomaly, negative vorticity and convergence at all levels except at 200 hPa which only showed a trough characteristic while for the dry season, this feature was only observed at low levels. Negative anomalies of vertical movements and temperature and positive humidity associated with EWDs were observed in both seasons, but reaching higher elevations during the dry period. The precipitation composition indicated that the EWDs are responsible for an increase of 16% (4%) during the wet season (dry). Typical characteristics of these waves were: mean period of 8 (73) days, wavelength of about 4500 (5500) km and phase velocity of about 6.5 (0.9) m.s-1 for the wet period (dry). The WRF model simulated the synoptic patterns, precipitation and the typical features associated with EWDs in a coherently manner when compared to the results of the Climate Forecast System Reanalysis (CFSR), but with greater intensity to the events of 11 and 17 June 2006.

L'objectif de cette thèse est de bâtir un modèle prédictif d'une unité industrielle continue de dissolution de minerai de nickel/ et précipitation de fer simultanées. La méthodologie adoptée consiste dans un premier temps à réaliser des expériences en réacteur fermé ou semi-fermé pour identifier et modéliser séparément les différents phénomènes en jeu, puis dans un second temps à construire un modèle d'unité continu les rassemblant tous. Ce travail a permis certaines avancées tant sur le système retenu que sur les méthodes et modèles adoptés. Tout d'abord la thermodynamique des solutions Na-Fe(III)-H2SO4 aux alentours de 100°C a été clarifiée et la constante de solubilité de la natrojarosite Na0,84H0,16Fe2,90(SO4)2(OH)5,7 a été déduite. Ensuite les paramètres influant sur les précipitations des jarosite de sodium et potassium ont été mis en évidence et leur cinétiques de croissance ont pu être déterminées grâce à l'utilisation originale de la méthode des caractéristiques. Puis les cinétiques de dissolution des minerais ont été obtenues, en prenant en compte différentes phases du minerai et leurs granulométries. Après avoir déduit tous les paramètres nécessaires, un modèle original, permettant de simuler une cascade de réacteur de dissolution/précipitation avec recyclage a été construit puis validé. Cet outil de conception, couplé à une étude technico-économique peut permettre d'optimiser le procédé. / The goal of this thesis is to build a predictive model for a continuous industrial unit combining simultaneously nickel ore dissolution and iron precipitation. A two steps method was used. First experiments in batch or semi-batch reactors were done in order to understand and model separately the different phenomenon that take place. Then all the equations and associated constants were used to build a model. This work allowed some advances both on the system studied and on the methods used. First thermodynamics of Na-Fe(III)-H2SO4 solutions in the 70-100°C temperature range was clarified and natrojarosite solubility constant Na0,84H0,16Fe2,90(SO4)2(OH)5,7 was deduced. Secondly parameters acting on sodium and potassium jarosite precipitation were highlighted and their growth kinetics were deduced through an original use of caracterisctics method. Then ore dissolution kinetics were found, taking into account different ores phases and their granulometry. After deducing all the necessary parameters, an original model allowing to simulate a cascade of dissolution/precipitation reactors with recyling loop was build and validated. This conception tool, coupled with a technico-economic study allows the optimisation of such a process.

Cette thèse traite de la morphodynamique pluri-décennale des environnements sableux soumis à l’action combinée des courants de marée et des vagues. Dans ces travaux, une approche de modélisation numérique basée sur les processus hydro-sédimentaires est utilisée, celle-ci étant combinée à des données collectées au sein du laboratoire LIENSs au cours de ces dernières années sur différents sites ateliers (i.e. données bathymétriques, hydrodynamiques, et sismiques). Afin de réaliser des simulations morphodynamiques de ce type d’environnement côtier, deux principaux développements ont été apportés au modèle : (1) la prise en compte de l’hétérogénéité granulométrique du sédiment, et (2) le calcul de l’évolution du fond suivant un schéma numérique basé sur le formalisme WENO et adapté à un maillage de type non-structuré. La première application du système de modélisation a concerné l'évolution rétrospective du banc de sable estuarien de la Longe de Boyard, situé en baie de Marennes-Oléron, sur la période 1960 à 2000. Les résultats de ce travail ont tout d'abord montré qu'en dépit de la morphologie globale du banc indiquant une dominance de la marée, la contribution des vagues était essentielle pour expliquer son évolution pluri-décennale. La rythmicité de l’accrétion sédimentaire observée et simulée dans la partie sud du banc a ensuite pu être analysée d'un point de vue stratigraphique, grâce à la prise en compte de l'hétérogénéité granulométrique du sédiment dans le modèle. Le lien entre la rythmicité des variations granulométriques du sédiment dans cette zone et la saisonnalité du climat de vagues a ainsi été mis en évidence. La seconde application du système de modélisation a concerné l’évolution de la flèche sableuse d'Arçay, située sur le littoral vendéen. Les résultats de simulation morphodynamique suggèrent que l’évolution de cette flèche, principalement contrôlée par les vagues, présente un comportement combinant allo-cyclicité et auto-cyclicité. / This work considers the pluri-decadal morphodynamic modeling of sandy coastal environments subjected to tides and waves. A process-based numerical model is used, together with a collection of field data from the LIENSs laboratory (i.e. bathymetric, hydrodynamic, and seismic data). To simulate the morphodynamic of these environments, two main numerical developments have been added to the model : (1) the sediment heterogeneity, and (2) the bed evolution computation following a WENO-based scheme adapted to unstructured grids. The first model application concerns the 40-year hindcast (period 1960 to 2000) of a wide estuary mouth sandbank located in the Marennes-Oléron bay : the Longe de Boyard sandbank. Numerical results suggest that this sandbank long-term evolution is strongly controlled by waves, in spite of its global tide-dominated morphology. Rhythmic nature of sediment accretion, which is observed and modeled in the south part of the bank, was then analyzed from a stratigraphic point of view thanks to the integration of sediment heterogeneity into the model. Grain size rhythmic variations in this area appeared to be related to local wave climate seasonality. The second model application considers the Arçay sandspit evolution (Vendée coast). Morphodynamic results suggest that this sandspit evolution, mainly controlled by waves, is characterized by a combination of both autogenic and external-influenced behaviors.

This focuses on hydrological transport in shallow catchments with topography-driven flow paths. The thesis gives new insight to kinematic pathway models for estimation of tracer discharge at the catchment outlet. A semi-analytical methodology is presented for transient travel time and age distributions referred to as "kinematic pathway approach“(KPA) that accounts for dispersion at two levels of morphological and macro-dispersion. Macro-dispersion and morphological dispersion components are reflected in KPA by assuming an effective Péclet number and topographically driven pathway length distributions, respectively. The kinematic measure of the transport, defined as a characteristic velocity of water flow through the catchment is obtained from the overall water balance in the catchment. To include transformation process in its simplest form of linear decay/degradation a framework is presented that solves one-dimensional reactive transport with numerically simulated travel times as the independent variable. The proposed KPA and coupled transport framework for quantifying tracer discharge at the shallow catchment outlet are applied to two selected catchments in Sweden. KPA is applied to modeling of a 23-year long chloride data series for the Kringlan catchment whereas the implantation of the framework for quantifying natural attenuation is illustrated for the Forsmark catchment. Numerical simulations of Forsmark catchment advective travel times are obtained by means of particle tracking using the fully-integrated flow model MIKE SHE. The KPA is found to provide reasonable estimates of tracer discharge distribution when considering the transport controlled by hillslope processes associated with short topographically driven flow paths to adjacent discharge zones, e.g. rivers and lakes. Simulated natural attenuation for Forsmark is also estimated well provided that the pathway length distribution is skewed toward short pathway lengths. This fact is indicative of the controlling impact of topography on flow path length and travel time distributions in shallow catchments. Our work has shown that the pathway (Lagrangian) methodologies are promising as predictive tools for hydrological transport. / <p>QC 20170928</p>

Abstract In the design and development of unit operations in chemical engineering, experimental testing is often very expensive or even impossible to perform. In these kinds of situations, numerical simulations offer a good approach to study process characteristics. Typically in chemical engineering, data-based modeling is applied to study the process. This requires many experiments for tuning the model parameters and validating the model. In a phenomenology-based approach, the evolution of the system is dictated by fluid and particle transport equations. These equations are independent of the process, and can therefore be applied in various systems. However, depending on the system, there are several aspects that have to be taken into account in order to choose the correct model for the problem in hand. In this work, computational fluid dynamics (CFD) and discrete element method (DEM) modeling have been applied in different unit operations in the field of chemical engineering. CFD was applied in preventing sedimentation in a tube heat exchanger, estimating the cooling efficiency of a vertical water jet onto a hot metal plate, and studying the formation of the slag free open-eye area on the steel ladle. For comparison, DEM was applied in the continuous high-shear granulation of pharmaceutical powder. The different models used in this work are reviewed, and the results are presented from the point of view of model and process development. The grid aspects in CFD simulations and the termination criteria for DEM and CFD simulations are also studied. Based on the results of this work, phenomenological modeling can be considered to be an efficient tool for unit operation design. Together with experimental work, different modeling strategies offer a powerful tool for the design and development of unit operations. / Tiivistelmä Kemiantekniikan yksikköprosessien suunnittelussa kokeellinen tutkimus on usein erittäin kallista ja joskus jopa mahdotonta toteuttaa. Tällöin mallinnus tarjoaa hyvän lähestymistavan yksikköprosessin ominaisuuksien tutkimiselle. Tyypillisesti kemianteollisuudessa, datapohjaista mallinnusmenetelmiä on käytetty systeemin ominaisuuksien tutkimiseksi. Tämä menetelmä vaatii kuitenkin paljon koetoistoja mallin parametrien virittämiseksi ja mallin validoimiseksi. Ilmiöpohjaisessa mallinnuksessa, systeemin aikakehitys määräytyy fluidi- ja partikkelivirtauksia määräävien kuljetusyhtälöiden perusteella. Nämä yhtälöt ovat prosessista riippumattomia, jolloin niitä voidaan soveltaa yleisesti mihin tahansa systeemiin. Riippuen kuitenkin tutkittavasta yksikköprosessista, eri asioita on otettava huomioon, jotta oikea malli voidaan valita kuvaamaan systeemiä. Tässä työssä virtauslaskentaa (computational fluid dynamics, CFD) ja partikkelimenetelmää (discrete element method, DEM) on käytetty erilaisten kemiantekniikan yksikköprosessien tutkimuksessa. CFD:n avulla on tutkittu putkilämmönvaihtimen sedimentaation ehkäisyä, laminaarisen suorasuihkujäähdytyksen tehokkuutta teräslevyn jäähdytyksessä, sekä senkkaprosessissa teräksen pinnalle ilmestyvän kuonasilmäkkeen muodostumista. DEM mallinnusta käytettiin lääkejauheen jatkuvatoimisen rakeistuksen tutkimuksessa. Mallinnuksessa käytetyt mallit esitellään ja niiden tulokset esitellään malli- ja prosessikehityksen näkökulmasta. Työssä on otettu myös esille mallinnustekniset asiat CFD:n vaatiman laskentahilan ja laskennan lopetuskriteerien näkökulmasta. Työssä esitettyjen tulosten perusteella voidaan todeta, että ilmiöpohjainen mallinnus on tehokas työkalu yksikköprosessien suunnittelussa. Yhdessä kokeellisen tutkimuksen kanssa, eri mallinnusmenetelmät tehostavat yksikköprosessien suunnittelua ja kehitystä.

Dans les chaînes de montagnes, les séismes de magnitudes intermédiaires à fortes (Mw>6) déclenchent systématiquement un grand nombre de glissements de terrain responsables de l'introduction de volumes massifs de sédiments dans le réseau fluviatile. L'évacuation progressive de ces sédiments hors de la zone épicentrale affecte la dynamique des rivières et provoque des aléas hydro-sédimentaires dans les plaines alluviales (avulsion des rivières, crues...). La quantification des transferts sédimentaires est essentielle pour mieux comprendre l'évolution des paysages à court et moyen terme (de l'heure au siècle) et permettre une gestion efficace des risques dans les zones d'accumulation. Cependant, les flux de sédiments grossiers étant difficiles à mesurer, les facteurs contrôlant l'évacuation des glissements de terrain restent à ce jour mal compris. Cette thèse a donc porté sur l'étude, via la modélisation, des paramètres influençant la mobilisation des glissements de terrain, la préservation de la capacité de transport la transition entre gorge et plaine alluviale et la dynamique court terme des cônes alluviaux soumis à de forts apports sédimentaires. Les approches développées sont appliquées au contexte de la côte Ouest de la Nouvelle Zélande où la probabilité d'occurrence d'un séisme de magnitude 8 est de 50% dans les 50 ans à venir. Cette problématique à été abordée analytiquement et via une approche numérique avec le modèle 2D d'évolution des paysages et des rivières, Eros. Avec l'approche analytique, nous démontrons que la conservation de la capacité de transport long terme à la transition entre gorges et plaines alluviales est généralement réalisée par le passage à un système en tresse. Nous identifions aussi la variabilité des débits comme facteur dominant de la capacité de transport long terme comparé à l'effet de la végétation riparienne. Avec l'approche numérique, nous utilisons Eros qui est composé 1. d'un modèle hydrodynamique 2D, 2. d'un modèle de transport/dépôt de sédiments et 3. de modèles gérant les flux latéraux d'érosion et de dépôt. La combinaison de ces éléments permet l'émergence de diverses géométries de rivières alluviales (droites/sinueuses ou en tresses) en fonction des forçages externes qu'elles subissent (débit d'eau, flux sédimentaires). L'application d'Eros à des cas naturels a nécessité la validation et la calibration de ses paramètres principaux à l'aide: 1. de solutions analytiques et 2. de la reproduction morphodynamique de systèmes naturels, tel que l'évolution de la rivière Poerua en Nouvelle Zélande suite au glissement de terrain du Mont Adams. Dans la partie aval du bassin, les simulations numériques démontrent les capacités du modèle 1) à prédire efficacement l'évolution de plaines alluviales soumises à plusieurs scénario d'apports sédimentaires massifs et 2) à générer des cartes de risques probabilistes. Dans la partie amont du bassin, les résultats mettent en évidence le rôle clef de la réduction dynamique de largeur des rivières par rapport à la largeur de la gorge fluviatile, sur l'accélération de l'évacuation des sédiments issus des glissements de terrain. Une loi unique caractérisant les temps d'export d'une distribution de glissements de terrain peut être définie en fonction du rapport entre volume de sédiment et capacité de transport initiale de la rivière, permettant ainsi d'estimer leur temps de résidence moyen à 5-30 ans pour un scénario de séisme de Mw=8 beaucoup plus faibles que ceux estimés précédemment (~100 ans). L'approche numérique développée dans ce travail suggère que l'étude de la réponse des chaînes de montagnes à un forçage sismique fort ne peut être effectuée efficacement qu'avec un modèle 2D capable de prendre en compte les non-linéarités entre écoulements des rivières, leurs géométries et le transport sédimentaire. Les résultats obtenus permettent une meilleure caractérisation de la dynamique des paysages à l'échelle du cycle sismique et des aléas à court terme. / In mountainous areas, intermediate to large earthquakes (Mw > 6) systematically trigger a large number of landslides supplying the fluvial network with massive volumes of sediment. The progressive evacuation of the sediment out of the epicentral area alters river dynamics and may cause hydro-sedimentary hazards in alluvial plains (river avulsion, inundations, bank erosion, ...). The quantification of sediment transfers is critical to better understand landscape evolution on short timescales (i.e. hours to centuries) and improve hazard management in deposition areas. However, the factors controlling the coarse sediment transfers are still poorly known due to a lack of field measurements and adequate numerical models. The aim of this work is thus to study, via numerical modeling, the parameters influencing landslides evacuation, the transport capacity variations at the gorge/alluvial plain transition and the short-term dynamics and hazards of alluvial fans. This work is set up in the context of the West Coast of New Zealand (NZ) which presents a 50% probability to experience a magnitude 8 earthquake in the next 50 years. This problematic has been addressed analytically and via a numerical approach. Using the analytical approach, we demonstrate that the conservation of long-term transport capacity at the bedrock gorge and alluvial plain transition usually implies the channel narrowing in the alluvial part that is generally realized by a transition to a braided system. We identify discharge variability as the dominant factor of alluvial river long term transport capacity compared to riparian vegetation. To explore the role of channel self-organization on coarse sediment transport, we use Eros, a 2D morphodynamic model able to simulate landscape evolution improved by a new 2D hydrodynamic model. Combined with a sediment transport/deposition model and lateral fluxes modeling (bank erosion and transverse deposition), Eros allows for the emergence of diverse alluvial river regimes and geometries (e.g. straight/sinuous and braided channels) as a function of the external forcing experienced by the river (water and sediment fluxes). The application of Eros on natural cases has required the validation and calibration of its principal parameters using analytical solutions and the morphodynamic reproduction of natural systems such as the evolution of the Poerua river in New Zealand following the Mount Adams landslide. In the downstream part of the catchment, the ensemble numerical simulations demonstrate Eros abilities to 1) efficiently predict the morphodynamic evolution of alluvial fans submitted to different scenarios of large sediment supplies and 2) generate probabilistic risk maps. In the upstream part, the results highlight the dominant role of dynamic river narrowing reducing export times of landslide-derived sediments. We define a new law characterizing export times as a function of landslide volume and pre-landslide transport capacity that predicts mean residence times for a M8 earthquake in a mountain range of 5-30 yr, much lower than previous estimations of ~ 100 yr. The numerical approach developed in this work suggests that the study of mountain ranges response to severe landslide disruption can only be addressed with a 2D model able to account for the non-linearities between river flow, channel geometry and sediment transport. The results allow for a better characterization of landscape dynamics at the scale of a seismic cycle and hydro-sedimentary hazards in the short term.

Dans cette thèse, nous avons développé un code de simulation numérique directe pour l’étude des écoulements particulaires. Le schéma numérique est basé sur une technique de domaines fictifs. Le code est validé sur de nombreux cas test puis nous l’avons utilisé pour étudier la sédimentation de particules bidimensionnelles en milieu confiné. Trois cas ont été analysés : sédimentation d’une particule unique, d’un doublet de particules et d’un grand nombre de particules. Dans le premier cas nous retrouvons le phénomène de survitesse qui apparaît pour une particule excentrée à bas nombre de Reynolds. Nous montrons que cette survitesse est très sensible à l’inertie du fluide : elle diminue lorsqu’on augmente le nombre de Reynolds. Cet effet est retardé par le confinement. Dans le cas d’un doublet de particules, nous retrouvons les comportements complexes observés dans la littérature (hystérésis, cascade sous-harmonique et chaos). Nous montrons qu’une nouvelle série de bifurcations et un nouvel attracteur apparaissent pour des particules plus pesantes. Il s’agit là d’une transition vers le chaos par la voie de la quasi-périodicité. Nous donnons le diagramme de bifurcation étendu. La nouvelle branche correspond à une structure horizontale qui conduit à une sédimentation lente. Dans le cas d’un grand nombre de particules, nous montrons que la vitesse de chute de l’interface fluide-particules suit une loi de type Richardson-Zaki, mais avec un exposant d’environ 4. Comme pour des sphères, la valeur de cet exposant dépend du confinement. Enfin, nous observons un phénomène de blocage, inattendu pour des particules non-cohésives, dû au caractère bidimensionnel de la suspension / In the present thesis, a fully-resolved numerical code has been developed for the analysis of particle-laden flows. A fictitious domain method is used. First, this numerical tool has been validated by using classical benchmarks. It has then been used to simulate the complex sedimentation of particles in three generic two-dimensional configurations: a single particle, a particle pair and a large number of particles in a confined domain. In the first case, the peak-velocity of an off-centred inclusion is recovered at low-Reynolds number. It is shown that this peak-velocity is very sensitive to fluid inertia: the peak-velocity decreases when the Reynolds number increases. This effect is delayed by the confinement. The very complex dynamics of a pair of particles sedimenting in a confined domain, observed in the litterature, is recovered (hysteresis, period-doubling cascade and chaos). It is shown that a new series of bifurcations, leading to a new attractor, emerges when the non-dimensional particle weight is increased. This new transition corresponds to a quasi-periodic route. The extended bifurcation diagram is given. The new branch discovered in this work corresponds to a nearly horizontal particle doublet, with a slow settling velocity. In the case of the settling of large number of particles, a RZ-like law is recovered for the sedimentation velocity of the fluid-particle interface. The exponent is close to 4, in contrast with the case of spheres. Finally, the sedimentation velocity at the end of the settling process is observed to be significantly reduced, like for cohesive sediments. This unexpected behaviour is related to the two-dimensionality of the suspension

The geological repository for final storage of spent nuclear fuel, envisioned by the Swedish Nuclear Fuel and Management Company (SKB), relies on several barriers: copper canisters deposited in holes in the floor of underground tunnels in deep bedrock, embedded in a buffer of compacted bentonite. The initially unsaturated buffer would take up water from the surrounding rock mass and swell to seal any potential gap. This initial two-phase (gas and liquid) regime with two components (air and water) may impact the final density, swelling pressure and biogeochemical conditions in the buffer. A main objective of this work is to identify factors and mechanisms that govern deposition hole inflow and bentonite wetting under the prevailing two-phase flow conditions in sparsely fractured bedrock. For this purpose, we use the numerical code TOUGH2 to perform two-phase flow simulations, conditioned by a companion field experiment (the Bentonite Rock Interaction Experiment or BRIE) performed in a 417 m deep tunnel of the Äspö Hard Rock Laboratory in southeastern Sweden. The models predict a significant de-saturation of the rock wall, which was confirmed by field data. To predict the early buffer wetting rates and patterns, the position of local flowing fractures and estimates of local rock matrix permeability appear more important than the total open hole groundwater inflow. A global sensitivity analysis showed that the buffer wetting time and the persistence of unsaturated conditions over extended periods of time in the rock depend primarily on the local fracture positions, rock matrix permeability, ventilation conditions in the tunnel and pressure far in the rock. Dismantling photographs from BRIE were used to reconstruct a fine-scale snapshot of saturation at the bentonite/rock interface, showing tremendous spatial variability. The high level of heterogeneity in the rock generates complex two-phase flow phenomena (air trapping, dissolution), which need to be accounted for in buffer design and rock suitability criteria. In particular, results suggest that uncertainties regarding two-phase flow behavior are relatively high close to residual air saturation, which may also have important implications for other applications involving two-phase flows, such as geological storage of carbon dioxide. / <p>At the time of the doctoral defense, the following papers were unpublished and had a status as follows: Paper 1: Manuscript. Paper 4: Manuscript.</p>

Shippagan Gully is a tidal inlet located near Shippagan, New Brunswick on the Gulf of Saint Lawrence. It is a particularly complex tidal inlet due to the fact that its tidal lagoon transects the Acadian peninsula and is open to the Bay des Chaleurs at its opposite end. As such, two open boundaries with phase lagged tidal cycles drive flow through the inlet, alternating direction with each tide and reaching velocities in excess of 2m/s. Hydrodynamic and morphological processes at the site are further complicated by the presence of a highly variable wave climate. Presently, shipping practices through the inlet are limited due to continual sedimentation within and immediately offshore from Shippagan Gully. As such, an extensive field study, desktop analysis and numerical and morphological modeling of Shippagan Gully have been conducted in order to provide guidance for future works. Modeling was conducted using the CMS-Wave and CMS-Flow numerical modeling system.

This doctorate thesis focuses on investigating the parameter space involved in numerically modeling the bandgap engineering of a GaInP/InGaAs/Ge lattice matched multi-junction solar cell (MJSC) using InAs/InGaAs quantum dots (QDs) in the middle sub-cell. The simulation environment – TCAD Sentaurus – solves the semiconductor equations using finite element and finite difference methods throughout well-defined meshes in the device to simulate the optoelectronic behavior first for single junction solar cells and subsequently for MJSCs with and without quantum dots under concentrated illumination of up to 1000 suns’ equivalent intensity. The MJSC device models include appropriate quantum tunneling effects arising in the tunnel junctions which serve as transparent sub-cell interconnects. These tunneling models are calibrated to measurements of AlGaAs/GaAs and AlGaAs/AlGaAs tunnel junctions reaching tunneling peak current densities above 1000 A/cm^2. Self-assembled InAs/GaAs quantum dots (QDs) are treated as an effective medium through a description of appropriate generation and recombination processes. The former includes analytical expressions for the absorption coefficient that amalgamates the contributions from the quantum dot, the InAs wetting layer (WL) and the bulk states. The latter includes radiative and non-radiative lifetimes with carrier capture and escape considerations from the confinement potentials of the QDs. The simulated external quantum efficiency was calibrated to a commercial device from Cyrium Technologies Inc., and required 130 layers of the QD effective medium to match the contribution from the QD ground state. The current – voltage simulations under standard testing conditions (1 kW/cm^2, T=298 K) demonstrated an efficiency of 29.1%, an absolute drop of 1.5% over a control structure. Although a 5% relative increase in photocurrent was observed, a 5% relative drop in open circuit voltage and an absolute drop of 3.4% in fill factor resulted from integrating lower bandgap nanostructures with shorter minority carrier lifetimes. However, these results are considered a worst case scenario since maximum capture and minimum escape rates are assumed for the effective medium model. Decreasing the band offsets demonstrated an absolute boost in efficiency of 0.5% over a control structure, thus outlining the potential benefits of using nanostructures in bandgap engineering MJSCs.