This work focuses on diffusion and advection through cementitious media, the work arises from two research contracts undertaken at Loughborough University: Experiments to Demonstrate Chemical Containment funded by UK NDA and the SKIN project, funded by the European Atomic Energy Community's Seventh Framework Programme. Diffusion will be one of the most significant mechanisms controlling any radionuclide migration from a nuclear waste, deep geological disposal facility. Advection may also occur, particularly as the immediate post closure groundwater rebound and equilibration proceeds but is expected to be limited by effective siting and management during the operational phase of the facility. In this work advection is investigated at laboratory scale as a possible shorter timescale technique for providing insight into the much slower process of diffusion. Radial techniques for diffusion and advection have been developed and the developmental process is presented in some detail. Both techniques use a cylindrical sample geometry that allows the radionuclide of interest to be introduced into a core drilled through the centre of the test material. For diffusion the core is sealed and submerged in a container of receiving solution which is sampled and analysed as the radionuclide diffuses into it. For advection, a cell has been designed that allows inflow via the central core to pass through the sample in a radial manner and be collected as it exits from the outer surface. The radionuclide of interest can be injected directly into the central core without significant disturbance to the advective flow. Minor improvements continue to be made but both techniques have provided good quality, reproducible results. The majority of the work is concentrated on a potential cemetitious backfill known as NRVB (Nirex Reference Vault Backfill) this is a high porosity, high calcium carbonate content cementitious material. The radioisotopes used in this work are 3H (in tritiated water), 137Cs, 125I, 90Sr, 45Ca, 63Ni, 152Eu, 241Am along with U and Th salts. In addition the effect of cellulose degradation products (CDP) on radioisotope mobility was investigated by manufacturing solutions where paper tissues were degraded in water, at 80°C, in the absence of air and at high pH due to the presence of the components of NRVB. All diffusion experiments were carried out under a nitrogen atmosphere. All advection experiments were undertaken using an eluent reservoir pressurised with nitrogen where the system remained closed up to the point of final sample collection. Results for tritiated water and the monovalent ions of Cs and I were produced on a timescale of weeks to months for both diffusion and advection. The divalent ions of Sr, Ca and Ni produced results on a timescale of months to years. Variations of the experiments were undertaken using the CDP solutions. The effects of CDP were much more apparent at radiotracer concentration than the much higher radiotracer with non-active carrier, concentration. In the presence of CDP Cs, I and Ni were found to migrate more quickly; Sr and Ca were found to migrate more slowly. Additional Sr experiments were undertaken at elevated ionic strength to evaluate the effect of the higher dissolved solids content of the CDP solutions. Some of the results for HTO, Cs, I and Sr have been modelled using a simple numerical representation of the system in GoldSim to estimate effective diffusivity and partition coefficient. The diffusion model successfully produced outputs that were comparable to literature values. The advection model is not yet producing good matches with the observed data but it continues to be developed and more processes will be added as new results become available. Autoradiography has been used to visualise the radionuclide migration and several images are reproduced that show the fate of the radiotracers retained on the NRVB at the end of the experiments. As the experimental programme progressed it was clear that results could not be produced in a suitable timescale for Eu, Am U and Th. These experiments have been retained and will be monitored every six months until either diffusion is detected or the volume of receiving liquid is inadequate to ensure the NRVB is saturated.

Master of Science / Department of Mathematics / Nathan Albin / This work utilizes finite differences to approximate the first derivative of non-periodic smooth functions. Math literature indicates that stabilizing Partial Differential Equation solvers based on high order finite difference approximations of spatial derivatives of a non-periodic function becomes problematic near a boundary. Hagstrom and Hagstrom have discovered a method of introducing additional grid points near a boundary, which has proven to be effective in stabilizing Partial Differential Equation solvers. Hagstrom and Hagstrom demonstrated their method for the case of the one-dimensional advection equation using spatial derivative approximations of even orders up to twenty-second order. In this dissertation, we explore the efficacy of the Hagstrom and Hagstrom method for the same Partial Differential Equation with spatial derivative approximations of odd orders and orders higher than twenty-two and report the number and locations of additional grid points required for stability in each case.

Age-model derived sediment mass accumulation rates (MARs) are consistently higher than 230Th-normalized MARs in the Equatorial Pacific Ocean during the past 25 ka. The offset, being highest in the Panama Basin, suggests sediment redistribution in this region is prominent. I test the hypothesis that downslope transport of sediments from topographically highs that surround the Panama Basin is the cause of higher-than-expected xs230Th inventories in the deeper parts of the basin. There is little difference in xs230Th inventories between the highest and lowest reaches of the basin suggesting that the topographic highs did not serve as a source of xs230Th. A spatial analysis suggests that there may be an enhanced scavenging of xs230Th closest to the equator in productive waters. To examine whether lateral mixing of productive equatorial waters with adjacent waters delivers xs230Th to the Panama Basin, I measured dissolved 230Th in eight deep-water casts within the Guatemala, Panama, and Peru Basins along a meridional transect at ~86°W. Below 1000 m, the Panama Basin shows the highest deficit (~50%) of 230Th in deep waters assuming a reversible exchange of 230Th between dissolved and sinking particulate matter. Peru Basin waters have a larger range of dissolved 230Th concentrations (7.9-16.5 fg/kg) than that within Panama Basin waters (5.7-7.1 fg/kg). There is a progressive decrease, suggesting advection, in average dissolved deep-water (>1000 m) 230Th concentrations from the southernmost sites in the Peru Basin toward the Panama Basin. My calculations suggest that advected 230Th is between 15-30% of the total 230Th being produced within waters of the Panama Basin itself. In the Panama Basin, the averaged biogenic barium and opal MARs suggest that productivity was greater during the Holocene (0-13000 years) than that during the last glacial (13000-25000 years) suggesting higher productivity during the Holocene. Uauth, however, is higher in sediments deposited during the last glacial than in those deposited during the Holocene, suggesting that low bottom water oxygen contents rather than respiration of organic matter drives Uauth enrichment. This oxygen depletion during the last glacial suggests that bottom waters were enriched in respired carbon, which, in turn, could be a driver of lower glacial atmosphere pCO2 values.

Theoretical and numerical studies of chaotic mixing are performed to circumvent the difficulties of efficient mixing, which come from the lack of turbulence in microfluidic devices. In order to carry out efficient and accurate parametric studies and to identify a fully chaotic state, a spectral element algorithm for solution of the incompressible Navier-Stokes and species transport equations is developed. Using Taylor series expansions in time marching, the new algorithm employs an algebraic factorization scheme on multi-dimensional staggered spectral element grids, and extends classical conforming Galerkin formulations to nonconforming spectral elements. Lagrangian particle tracking methods are utilized to study particle dispersion in the mixing device using spectral element and fourth order Runge-Kutta discretizations in space and time, respectively. Comparative studies of five different techniques commonly employed to identify the chaotic strength and mixing efficiency in microfluidic systems are presented to demonstrate the competitive advantages and shortcomings of each method. These are the stirring index based on the box counting method, Poincare sections, finite time Lyapunov exponents, the probability density function of the stretching field, and mixing index inverse, based on the standard deviation of scalar species distribution. Series of numerical simulations are performed by varying the Peclet number (Pe) at fixed kinematic conditions. The mixing length (lm) is characterized as function of the Pe number, and lm ∝ ln(Pe) scaling is demonstrated for fully chaotic cases. Employing the aforementioned techniques, optimum kinematic conditions and the actuation frequency of the stirrer that result in the highest mixing/stirring efficiency are identified in a zeta potential patterned straight micro channel, where a continuous flow is generated by superposition of a steady pressure driven flow and time periodic electroosmotic flow induced by a stream-wise AC electric field. Finally, it is shown that the invariant manifold of hyperbolic periodic point determines the geometry of fast mixing zones in oscillatory flows in two-dimensional cavity.

Previous studies have suggested that advection from the Mexican Plateau (MP) may influence rainfall over Texas in spring and summer; generally air ascends over the cordillera and descends over the southern plains. The two mechanisms may link the northern Mexico drought to Texas drought. Observations and the Community Earth System Model are used in this study to describe the 2011 Texas-northern-Mexico drought and examine the role of the MP on the hydro-climate over the southern US, providing implications for the linkage between the MP and rainfall over Texas. A control run and three experimental runs were performed with prescribed sea surface temperatures and sea ice fractions. The results show that when the MP becomes dry, rainfall declines locally and downstream. During the spring, the dry air brought to Texas by prevailing westerly winds suppresses local convection; but dry air advection from the highlands has little influence on rainfall over Texas during the summer when Texas is no longer in the downstream areas. During the summer, a warmer MP draws moist air over the peripheral low elevation areas to the highlands; it bends the low-level jet towards the highlands and an anti-cyclonic flow anomaly forms over the southern US, which causes air to diverge and tends to reduce rainfall over the southern US. / text

Η παρούσα Μεταπτυχιακή εργασία πραγματεύεται την μεταφορά διατηρητικού ρύπου, ο οποίος εκρέει από γραμμική πηγή σταθερού μήκους, συνεχούς αλλά και στιγμιαίας έγχυσης, σε διαφορετικά βάθη ενός ανοικτού αγωγού με επίπεδο και υπό κλίση πυθμένα. Η επίλυση του προβλήματος της μεταφοράς ρύπου έγινε με τον εμπορικό κώδικα ANSYS FLUENT v.13 και στη συνέχεια συγκρίθηκε με τα αποτελέσματα που προκύπτουν από την αναλυτική λύση. / This master thesis is dealing with the mass transport problem of a conservative neutral pollutant inside an open channel flow. The open channel is flat bottommed with a mild slope. The pollutant outflows either continuously or instantaneously out of a linear source,with specified length, which every time is positioned on different depths. The mass transport problem was solved both numerically, by ANSYS FLUENT v.13 programme, and analytically in order to get those two types of solytions compared.

xi, 172 p. : ill. (some col.) / Accretion shock waves are present in many important astrophysical systems and have been a focus of research for decades. These investigations provide a large body of understanding as to the nature, characteristics, and evolutionary behaviors of accretion shock waves over a wide range of conditions. However, largely absent are investigations into the properties of accretion shock waves in the presence of strong magnetic fields. In such cases these strong magnetic fields can significantly alter the stability behaviors and evolution of the accretion shock wave through the production and propagation of magnetic waves as well as magnetically constrained advection. With strong magnetic fields likely found in a number of accretion shock systems, such as compact binary and protostellar systems, a better understanding of the behaviors of magnetic accretion shock waves is needed. A new magnetohydrodynamics simulation tool, IMOGEN, was developed to carry out an investigation of instabilities in strong, slow magnetic accretion shocks by modelling their long-term, nonlinear evolution. IMOGEN implements a relaxed, second-order, total variation diminishing, monotonic upwind scheme for conservation laws and incorporates a staggered-grid constrained transport scheme for magnetic advection. Through the simulated evolution of magnetic accretion shocks over a wide range of initial conditions, it has been shown, for sufficiently high magnetic field strengths, that magnetic accretion shocks are generally susceptible to corrugation instabilities, which arise in the presence of perturbations of the initial shock front. As these corrugation instabilities grow, they manifest as magnetic wave propagation in the upstream region of the accretion column, which propagate away from the accretion shock front, and as density columns, or fingers, that grow into the higher density downstream flow, defined and constrained by current loops created during the early evolution of the instability. / Committee in charge: Dr. James Schombert, Chair; Dr. James Imamura, Advisor; Dr. Alan Rempel, Member; Dr. John Toner, Member; Dr. Kent Stevens, Outside Member

Philosophiae Doctor - PhD / In Africa as well as in the world, South Africa plays an important role in the mining industry which dates back almost 120 years. Mining activities in South Africa mainly take place in Gauteng Province. Every year millions of tons of rocks are taken from underground, milled and processed to extract gold. Uranium is one of the minerals which is mined as a by-product of gold. The uranium bearing tailings are then placed on large dumps which are usually one or more kilometres in diameter in the environment. There are approximately 250 gold mine tailings dumps covering a total area of about 7000 ha. These gold mine tailings dumps contain considerable amounts of radium (226Ra) and have therefore been identified as large sources of radon (222Rn). Radon release from these tailings dumps pose health concerns for the surrounding communities. This study reports on the 222Rn soil-gas concentrations and exhalations measurements which were conducted at a non-operational mine dump (Kloof) which belongs to Carletonville Gold Field, Witwatersrand, South Africa. Radon soil-gas concentrations were measured using a continuous radon monitor (the Durridge RAD7). The radon soil gas concentration levels were measured at depths starting from 30 cm below ground/air interface up to 110 cm at intervals of 20 cm. These radon soil-gas measurements were performed at five different spots on the mine dump. The concentrations recorded ranged from to kBq.m-3. Furthermore, thirty four soil samples were taken from these spots for laboratory-based measurement. The soil samples were collected in the depth range 0-30 cm and measured using the low-background Hyper Pure Germanium (HPGe) -ray detector available at the Environmental Radioactivity Laboratory (ERL), iThemba LABS, Western Cape Province. The weighted average activity concentrations in the soils samples were Bq.kg-1, Bq.kg-1 and Bq.kg-1 for 238U, 40K and 232Th, respectively. 1 26 23 472 7 308 5 2551 18 The activity concentration of 40K and the decay series of 238U and 232Th for the top 30 cm of the accessible parts of the mine dump were measured using the MEDUSA (Multi-Element Detector for Underwater Sediment Activity) -ray detection system which was mounted on the front of a 44 vehicle, 0.6 m off the ground. The spectra acquired by the MEDUSA system were analysed and the activity concentrations were extracted using the Full Spectrum Analysis approach. The average activity concentrations obtained were 25975 Bq.kg-1 for 40K, 30940 Bq.kg-1 for 238U series and 185Bq.kg-1 for 232Th series for 2002 survey while 26075 Bq.kg-1 for 40K, 31040 Bq.kg-1 for 238U series and 20 5 Bq.kg-1 for 232Th series were recorded for 2010 survey. Moreover, a novel technique by means of the MEDUSA -ray detection system was utilised to map radon exhalation from the Kloof mine dump. In this technique the 214Bi concentration in the top 30 cm of the dump, which is used as a proxy for the 238U or 226Ra concentration, is measured even though there is not secular equilibrium between the radium and bismuth because radon has escaped. The actual activity concentrations of 40K and 226Ra are obtained from soil samples which were collected, sealed and measured in the laboratory settings after attaining secular equilibrium since radon gas does not escape. By comparing the ratios of the activity concentrations of the 214Bi obtained in the field using the MEDUSA -ray detector system to the activity concentrations of the 214Bi obtained in the laboratory using the HPGe -ray detector system yields a reasonable radon exhalation for a particular location in the mine dump. In this case it is assumed that the 40K and 232Th series activity concentrations are not or hardly affected by radon exhalation. By applying this new technique an average normalised radon flux of 0.120.02 and 0.120.02 Bq.m-2.s-1 were obtained for 2002 and 2010 surveys, respectively. The electrets and the IAEA standard formula were also utilised to calculated radon fluxes yielding 0.060.02 and 0.120.02 Bq.m-2.s-1 respectively.

Cette thèse présente un modèle bi-dimensionnel pour le transport atmosphérique à grande échelle, nommé Pangolin, conçu pour passer à l'échelle sur les achitectures parallèles. La version actuelle comporte une advection 2D ainsi qu'un schéma linéaire de chimie et servira de base pour un modèle de chimie-transport (MCT). Pour obtenir la conservation de la masse, un schéma en volume finis de type van Leer a été retenu pour l'advection et étendu au cas 2D en utilisant des opérateurs alternés. La conservation de la masse est assurée en corrigeant les vents en amont. Nous proposons une solution au problème "des pôles" de la grille régulière latitude-longitude grâce à une nouvelle grille préservant approximativement les aires des cellules et couvrant la sphère uniformément. La parallélisation du modèle se base sur l'advection et utilise un algorithme de décomposition de domaines spécialement adapté à la grille. Cela permet d'obtenir l'équilibrage de la charge de calcul avec MPI, une librairie d'échanges de messages. Pour que les performances soient à la hauteur sur les architectures parallèles actuelles et futures, les propriétés analytiques de la grille sont exploitées pour le schéma d'advection et la parallélisation en privilégiant le moindre coût des flops par rapport aux mouvement de données. Le modèle est validé sur des cas tests analytiques et comparé à des schémas de transport à l'aide d'un comparatif récemment publié. Pangolin est aussi comparé au MCT de Météo-France via un schéma linéaire d'ozone et l'utilisation de coordonnées isentropes. / We present in this thesis the development of a large-scale bi-dimensional atmospheric transport scheme designed for parallel architectures with scalability in mind. The current version, named Pangolin, contains a bi-dimensional advection and a simple linear chemistry scheme for stratospheric ozone and will serve as a basis for a future Chemistry Transport Model (CTM). For mass-preservation, a van Leer finite volume scheme was chosen for advection and extended to 2D with operator splitting. To ensure mass preservation, winds are corrected in a preprocessing step. We aim at addressing the "pole issue" of the traditional regular latitude-longitude by presenting a new quasi area-preserving grid mapping the sphere uniformly. The parallelization of the model is based on the advection operator and a custom domain-decomposition algorithm is presented here to attain load-balancing in a message-passing context. To run efficiently on current and future parallel architectures, algebraic features of the grid are exploited in the advection scheme and parallelization algorithm to favor the cheaper costs of flops versus data movement. The model is validated on algebraic test cases and compared to other state-of-theart schemes using a recent benchmark. Pangolin is also compared to the CTM of Météo-France, MOCAGE, using a linear ozone scheme and isentropic coordinates.

Numerical simulation and remote sensing have indicated that the northern half of the Red Sea has a significant role in the thermohaline circulation within the basin. However, very few studies with in situ observation have been performed in a region where the formation of Red Sea Outflow Water (RSOW) and occasionally of Red Sea Deep Water (RSDW) take place during the winter in the northern Red Sea (NRS). This study provides new insights into the seasonal variability and the mechanisms that drive the thermohaline circulation of the north half Red Sea using high-resolution glider observations combined with reanalysis and satellite datasets. The study describes the water masses characteristics, the mesoscale activity, and the forcing mechanisms. In addition, we examine the biogeochemical responses to the physical drivers in the northern half of the Red Sea and how these processes alter the marine ecosystem. During winter, the mesoscale eddy activity and heat fluxes create the necessary conditions for the formation of RSOW in the NRS. The cyclonic circulation elevates relatively denser water in the surface, which is exposed to the atmosphere exchange. Thus, it leads to subduction of the surface layer forming of RSOW. The subducted water has been characterized by high oxygen as it has recently been ventilated. In addition, chlorophyll fluorescence has subducted along the isopycnals, contributing to exporting material below the sunlit layer. After the formation of RSOW, a period of strong anticyclonic circulation was observed In late February, which stirred and mixed the advected waters from the south in the northern region. It is accompanied by heat flux transition, and at the periphery of the observed Anticyclonic Eddy, an uplifting of the densest water to the surface occurred. The presence of the anticyclonic circulation enables the water advection from the south and extends the time of the surface water for atmospheric exposure. In April, the warmer intrusion of fresher waters from the south dominated the eastern part of the NRS, reestablishing the cyclonic circulation. To the best of our knowledge, this is the first in situ observation in the NRS that captured the seasonal progression of the transition of heat flux in wintertime and water advection that terminates the formation of RSOW. A continuous supply of northward warmer, lower salinity near the coast from the south is observed throughout the summertime period. Strong stratification with surface mixed layers no deeper than 25-30 meters due to the advection of lower salinity surface water and local heating. Another change that occurred during the summer period is that the source of low salinity inflow into the region transitioned from Gulf of Aden Surface Water (GASW) to Gulf of Aden Intermediate Water (GAIW)—assuming that the inflow of GAIW began with the onset of the Southwest Monsoonal winds in the south. The summertime heating and along basin evaporation set up the system for the wintertime cooling and additional evaporation that contributes to the formation of RSOW and RSDW. The mixed layer Price-Weller-Pinkel (PWP) model (Price et al., 1986) is implemented to quantify the influence of local heat fluxes compared with horizontal advection of the Gulf of Aden Water on the upper layer. Simulation of the mixed layer showed that advection was the major contributor to the seasonally integrated heat content and mixed layer simulation in summer. In contrast to winter, the timing of the mesoscale eddy activity, significant cooling, and advection add complexity to the region. The difference in the heat content was significant, and the PWP predicted an increasing mixed layer depth, while the observed mixed layer depth remained relatively constant. The differences between the calculated and simulated heat content were minimum during the absence of the mesoscale eddy and advection from the south. Overall, the quantification suggests a complex relationship between atmospheric forcing and advection on the heat content and the mixed layer depth.