Diffusion-driven transport in clayrock formations

Clay mineral-rich sedimentary formations are currently under investigation to evaluate their potential use as host formations for installation of deep underground disposal facilities for radioactive waste (e.g. Boom Clay (BE), Opalinus Clay (CH), Callovo-Oxfordian argillite (FR)). The ultimate safety of the corresponding repository concepts depends largely on the capacity of the host formation to limit the flux towards the biosphere of radionuclides (RN) contained in the waste to acceptably low levels. Data for diffusion-driven transfer in these formations shows extreme differences in the measured or modelled behaviour for various radionuclides, e. g. between halogen RN (Cl-36, I-129) and actinides (U-238,U-235, Np-237, Th-232, etc.), which result from major differences between RN of the effects on transport of two phenomena: diffusion and sorption. This paper describes recent research aimed at improving understanding of these two phenomena, focusing on the results of studies carried out during the EC Funmig IP on clayrocks from the above three formations and from the Boda formation (HU). Project results regarding phenomena governing water, cation and anion distribution and mobility in the pore volumes influenced by the negatively-charged surfaces of clay minerals show a convergence of the modelling results for behaviour at the molecular scale and descriptions based on electrical double layer models. Transport models exist which couple ion distribution relative to the clay-solution interface and differentiated diffusive characteristics. These codes are able to reproduce the main trends in behaviour observed experimentally, e.g. D-e(anion) < D-e(HTO) < D-e(cation) and D-e(anion) variations as a function of ionic strength and material density. These trends are also well-explained by models of transport through ideal porous matrices made up of a charged surface material. Experimental validation of these models is good as regards monovalent alkaline cations, in progress for divalent electrostatically-interacting cations (e.g. Sr2+) and still relatively poor for 'strongly sorbing', high K-d cations. Funmig results have clarified understanding of how clayrock mineral composition, and the corresponding organisation of mineral grain assemblages and their associated porosity, can affect mobile solute (anions, HTO) diffusion at different scales (mm to geological formation). In particular, advances made in the capacity to map clayrock mineral grain-porosity organisation at high resolution provide additional elements for understanding diffusion anisotropy and for relating diffusion characteristics measured at different scales. On the other hand, the results of studies focusing on evaluating the potential effects of heterogeneity on mobile species diffusion at the formation scale tend to show that there is a minimal effect when compared to a homogeneous property model. Finally, the results of a natural tracer-based study carried out on the Opalinus Clay formation increase confidence in the use of diffusion parameters measured on laboratory scale samples for predicting diffusion over geological time-space scales. Much effort was placed on improving understanding of coupled sorption-diffusion phenomena for sorbing cations in clayrocks. Results regarding sorption equilibrium in dispersed and compacted materials for weakly to moderately sorbing cations (Sr2+, Cs+, Co2+) tend to show that the same sorption model probably holds in both systems. It was not possible to demonstrate this for highly sorbing elements such as Eu(III) because of the extremely long times needed to reach equilibrium conditions, but there does not seem to be any clear reason why such elements should not have similar behaviour. Diffusion experiments carried out with Sr2+, Cs+ and Eu(III) on all of the clayrocks gave mixed results and tend to show that coupled diffusion-sorption migration is much more complex than expected, leading generally to greater mobility than that predicted by coupling a batch-determined K-d and Ficks law based on the diffusion behaviour of HTO. If the K-d measured on equivalent dispersed systems holds as was shown to be the case for Sr, Cs (and probably Co) for Opalinus Clay, these results indicate that these cations have a D-e value higher than HTO (up to a factor of 10 for Cs+). Results are as yet very limited for very moderate to strongly sorbing species (e.g. Co(II), Eu(III), Cu(II)) because of their very slow transfer characteristics. (C) 2011 Elsevier Ltd. All rights reserved.

On the seasonal variations of salinity of the tropical Atlantic mixed layer

The physical processes controlling the mixed layer salinity (MLS) seasonal budget in the tropical Atlantic Ocean are investigated using a regional configuration of an ocean general circulation model. The analysis reveals that the MLS cycle is generally weak in comparison of individual physical processes entering in the budget because of strong compensation. In evaporative regions, around the surface salinity maxima, the ocean acts to freshen the mixed layer against the action of evaporation. Poleward of the southern SSS maxima, the freshening is ensured by geostrophic advection, the vertical salinity diffusion and, during winter, a dominant contribution of the convective entrainment. On the equatorward flanks of the SSS maxima, Ekman transport mainly contributes to supply freshwater from ITCZ regions while vertical salinity diffusion adds on the effect of evaporation. All these terms are phase locked through the effect of the wind. Under the seasonal march of the ITCZ and in coastal areas affected by river (7°S:15°N), the upper ocean freshening by precipitations and/or runoff is attenuated by vertical salinity diffusion. In the eastern equatorial regions, seasonal cycle of wind forced surface currents advect freshwaters, which are mixed with subsurface saline water because of the strong vertical turbulent diffusion. In all these regions, the vertical diffusion presents an important contribution to the MLS budget by providing, in general, an upwelling flux of salinity. It is generally due to vertical salinity gradient and mixing due to winds. Furthermore, in the equator where the vertical shear, associated to surface horizontal currents, is developed, the diffusion depends also on the sheared flow stability.

Estimation and calibration of the water isotope differential diffusion length in ice core records

Palaeoclimatic information can be retrieved from the diffusion of the stable water isotope signal during firnification of snow. The diffusion length, a measure for the amount of diffusion a layer has experienced, depends on the firn temperature and the accumulation rate. We show that the estimation of the diffusion length using power spectral densities (PSDs) of the record of a single isotope species can be biased by uncertainties in spectral properties of the isotope signal prior to diffusion. By using a second water isotope and calculating the difference in diffusion lengths between the two isotopes, this problem is circumvented. We study the PSD method applied to two isotopes in detail and additionally present a new forward diffusion method for retrieving the differential diffusion length based on the Pearson correlation between the two isotope signals. The two methods are discussed and extensively tested on synthetic data which are generated in a Monte Carlo manner. We show that calibration of the PSD method with this synthetic data is necessary to be able to objectively determine the differential diffusion length. The correlation-based method proves to be a good alternative for the PSD method as it yields precision equal to or somewhat higher than the PSD method. The use of synthetic data also allows us to estimate the accuracy and precision of the two methods and to choose the best sampling strategy to obtain past temperatures with the required precision. In addition to application to synthetic data the two methods are tested on stable-isotope records from the EPICA (European Project for Ice Coring in Antarctica) ice core drilled in Dronning Maud Land, Antarctica, showing that reliable firn temperatures can be reconstructed with a typical uncertainty of 1.5 and 2 °C for the Holocene period and 2 and 2.5 °C for the last glacial period for the correlation and PSD method, respectively.

Diffusion boundary layers ameliorate the negative effects of ocean acidification on the temperate coralline macroalga Arthrocardia corymbosa

Anthropogenically-modulated reductions in pH, termed ocean acidification, could pose a major threat to the physiological performance, stocks, and biodiversity of calcifiers and may devalue their ecosystem services. Recent debate has focussed on the need to develop approaches to arrest the potential negative impacts of ocean acidification on ecosystems dominated by calcareous organisms. In this study, we demonstrate the role of a discrete (i.e. diffusion) boundary layer (DBL), formed at the surface of some calcifying species under slow flows, in buffering them from the corrosive effects of low pH seawater. The coralline macroalga Arthrocardia corymbosa was grown in a multifactorial experiment with two mean pH levels (8.05 'ambient' and 7.65 a worst case 'ocean acidification' scenario projected for 2100), each with two levels of seawater flow (fast and slow, i.e. DBL thin or thick). Coralline algae grown under slow flows with thick DBLs (i.e., unstirred with regular replenishment of seawater to their surface) maintained net growth and calcification at pH 7.65 whereas those in higher flows with thin DBLs had net dissolution. Growth under ambient seawater pH (8.05) was not significantly different in thin and thick DBL treatments. No other measured diagnostic (recruit sizes and numbers, photosynthetic metrics, %C, %N, %MgCO3) responded to the effects of reduced seawater pH. Thus, flow conditions that promote the formation of thick DBLs, may enhance the subsistence of calcifiers by creating localised hydrodynamic conditions where metabolic activity ameliorates the negative impacts of ocean acidification.

Temperature Evolution and Light Species Diffusion in Armor and Structural Material for Inertial Fusion Reactor Chambers: a Case for HiPER 4a

One of the most advance designs for HiPER fusion reactor is a spherical chamber 10 m in diameter based on dry wall concept. In this system, the first wall will have to withstand short energy pulses of 5 to 20 MJ at a repetition rate of 0.5-10 Hz mostly in form of X-rays and charged particles. To avoid melting of the inner surface, the first wall consists on a thin armor attached to the structural material. Thickness (th) and material of each layer have to be chosen to assure the proper functioning of the facility during its planned lifetime.

Optimizing phosphorus diffusion for photovoltaic applications: peak doping, inactive phosphorus,gettering, and contact formation

The phosphosilicate glass (PSG), fabricated by tube furnace diffusion using a POCl3 source, is widely used as a dopant source in the manufacturing of crystalline silicon solar cells. Although it has been a widely addressed research topic for a long time, there is still lack of a comprehensive understanding of aspects such as the growth, the chemical composition, possible phosphorus depletion, the resulting in-diffused phosphorus profiles, the gettering behavior in silicon, and finally the metal-contact formation. This paper addresses these different aspects simultaneously to further optimize process conditions for photovoltaic applications. To do so, a wide range of experimental data is used and combined with device and process simulations, leading to a more comprehensive interpretation. The results show that slight changes in the PSG process conditions can produce high-quality emitters. It is predicted that PSG processes at 860 °C for 60 min in combination with an etch-back and laser doping from PSG layer results in high-quality emitters with a peak dopant density Npeak = 8.0 × 1018 cm−3 and a junction depth dj = 0.4 μm, resulting in a sheet resistivityρsh = 380 Ω/sq and a saturation current-density J0 below 10 fA/cm2. With these properties, the POCl3 process can compete with ion implantation or doped oxide approaches.

Procedure for calculating the chloride diffusion coefficient and surface concentration from a profile having a maximum beyond the concrete surface

Chlorides induce local corrosion in the steel reinforcements when reaching the bar surface. The measurement of the rate of ingress of these ions, is made by mathematically fitting the so called “error function equation” into the chloride concentration profile, obtaining so the diffusion coefficient and the chloride concentration at the concrete surface. However, the chloride profiles do not always follow Fick’s law by having the maximum concentration at the concrete surface, but often the profile shows a maximum concentration more in the interior, which indicates a different composition and performance of the most external concrete layer with respect to the internal zones. The paper presents a procedure prepared during the time of the RILEM TC 178-TMC: “Testing and modeling chloride penetration in concrete”, which suggests neglecting the external layer where the chloride concentration increases and using the maximum as an “apparent” surface concentration, called C max and to fit the error function equation into the decreasing concentration profile towards the interior. The prediction of evolution should be made also from the maximum.

(Table 1) Conditions of vertical near-bottom profiling and the main parameters of the benthic boundary layer in several regions of the Black Sea shelf

Measurements of 14 vertical profiles of currents and hydrological parameters in the near-bottom layer with depth resolution of 0.1 m were carried out in several regions of the Black Sea shelf, at five points over the continental slope, and in three deep water regions. The upper boundary of the benthic boundary layer (BBL) was reliably determined at a point at distance from 5-7 to 35-40 m from the bottom where the gradients of density and current velocity changed. Experimental data obtained were used to determine the coefficient of bottom friction, friction velocity, coefficients of vertical diffusion of momentum and density, and vertical fluxes of temperature and salinity in the BBL.

Comparison of Lagrangian Coherent Structures and Relative Dispersion in a Mixing Layer

Thesis (Master's)--University of Washington, 2016-06

Ion beam modification and analysis of metal/polymer bi-layer thin films

A set of varying-thickness Au-films were thermally evaporated onto poly(styrene-co-acrylonitrile) thin film surfaces. The Au/PSA bi-layer targets were then implanted with 50 keV N+ ions to a fluence of 1 × 1016 ions/cm2 to promote metal-to-polymer adhesion and to enhance their mechanical and electrical performance. Electrical conductivity measurements of the implanted Au/PSA thin films showed a sharp percolation behavior versus the pre-implant Au-film thickness with a percolation threshold near the nominal thickness of 44 Å. The electrical conductivity results are discussed along with the film microstructure and the elemental diffusion/mixing within the Au/PSA interface obtained by scanning electron microscopy (SEM) and ion beam analysis techniques (RBS and ERD).

The influence of excipients on the diffusion of ibuprofen and paracetamol in gastric mucus

The aim of this study was to examine the diffusion of commonly administered analgesics, ibuprofen and paracetamol, through gastric mucus. As ibuprofen and paracetamol are often formulated with alkalising excipients, or are commonly co-administered with antacids that have been demonstrated to alter their absorption, diffusion was also studied in the presence of a range of soluble and insoluble antacids or buffering agents. The effect of pH, which has been demonstrated to modify the properties of mucus, was also studied. Mucus was a significant barrier to diffusion for both drugs, compared to an unstirred aqueous layer with diffusion rates significantly lower in the presence of a mucus barrier for both drugs; ibuprofen diffusion also demonstrated a significant increase in the lag time. Paracetamol diffusion was not significantly affected by addition of any antacid, whereas ibuprofen rates were affected and the diffusion lag time for ibuprofen was significantly reduced in all cases. Isolated increases in pH increased the rate and reduced the lag time for ibuprofen diffusion. It was shown that mucus acts as a passive barrier in the case of paracetamol diffusion, and an interactive barrier to ibuprofen diffusion. Changes in mucus viscosity at different pH values may be responsible for the observed changes in ibuprofen diffusion rate. © 2004 Elsevier B.V. All rights reserved.

Ionic and molecular diffusion in cementitious materials

The work described in this thesis is an attempt to provide improved understanding of the effects of several factors affecting diffusion in hydrated cement pastes and to aid the prediction of ionic diffusion processes in cement-based materials. Effect of pore structure on diffusion was examined by means of comparative diffusion studies of quaternary ammonium ions with different ionic radii. Diffusivities of these ions in hydrated pastes of ordinary portland cement with or without addition of fly ash were determined by a quasi-steady state technique. The restriction of the pore geometry on diffusion was evaluated from the change of diffusivity in response to the change of ionic radius. The pastes were prepared at three water-cement ratios, 0.35, 0.50 and 0.65. Attempts were made to study the effect of surface charge or the electrochemical double layer at the pore/solution interface on ionic diffusion. An approach was to evaluate the zeta potentials of hydrated cement pastes through streaming potential measurements. Another approach was the comparative studies of the diffusion kinetics of chloride and dissolved oxygen in hydrated pastes of ordinary portland cement with addition of 0 and 20% fly ash. An electrochemical technique for the determination of oxygen diffusivity was also developed. Non-steady state diffusion of sodium potassium, chloride and hydroxyl ions in hydrated ordinary portland cement paste of water-cement ratio 0.5 was studied with the aid of computer-modelling. The kinetics of both diffusion and ionic binding were considered for the characterization of the concentration profiles by Fick's first and second laws. The effect of the electrostatic interactions between ions on the overall diffusion rates was also considered. A general model concerning the prediction of ionic diffusion processes in cement-based materials has been proposed.

Kinetics and coverage dependent reaction mechanisms of the copper atomic layer deposition from copper dimethylamino-2-propoxide and diethylzinc

Atomic layer deposition (ALD) has been recognized as a promising method to deposit conformal and uniform thin film of copper for future electronic devices. However, many aspects of the reaction mechanism and the surface chemistry of copper ALD remain unclear. In this paper, we employ plane wave density functional theory (DFT) to study the transmetalation ALD reaction of copper dimethylamino-2-propoxide [Cu(dmap)2] and diethylzinc [Et2Zn] that was realized experimentally by Lee et al. [ Angew. Chem., Int. Ed. 2009, 48, 4536−4539]. We find that the Cu(dmap)2 molecule adsorbs and dissociates through the scission of one or two Cu–O bonds into surface-bound dmap and Cu(dmap) fragments during the copper pulse. As Et2Zn adsorbs on the surface covered with Cu(dmap) and dmap fragments, butane formation and desorption was found to be facilitated by the surrounding ligands, which leads to one reaction mechanism, while the migration of ethyl groups to the surface leads to another reaction mechanism. During both reaction mechanisms, ligand diffusion and reordering are generally endothermic processes, which may result in residual ligands blocking the surface sites at the end of the Et2Zn pulse, and in residual Zn being reduced and incorporated as an impurity. We also find that the nearby ligands play a cooperative role in lowering the activation energy for formation and desorption of byproducts, which explains the advantage of using organometallic precursors and reducing agents in Cu ALD. The ALD growth rate estimated for the mechanism is consistent with the experimental value of 0.2 Å/cycle. The proposed reaction mechanisms provide insight into ALD processes for copper and other transition metals.

Nano-Engineering of Densely Packed Electrochemical Energy Storage Architectures by Atomic Layer Deposition

Nanostructures are highly attractive for future electrical energy storage devices because they enable large surface area and short ion transport time through thin electrode layers for high power devices. Significant enhancement in power density of batteries has been achieved by nano-engineered structures, particularly anode and cathode nanostructures spatially separated far apart by a porous membrane and/or a defined electrolyte region. A self-aligned nanostructured battery fully confined within a single nanopore presents a powerful platform to determine the rate performance and cyclability limits of nanostructured storage devices. Atomic layer deposition (ALD) has enabled us to create and evaluate such structures, comprised of nanotubular electrodes and electrolyte confined within anodic aluminum oxide (AAO) nanopores. The V2O5- V2O5 symmetric nanopore battery displays exceptional power-energy performance and cyclability when tested as a massively parallel device (~2billion/cm2), each with ~1m3 volume (~1fL). Cycled between 0.2V and 1.8V, this full cell has capacity retention of 95% at 5C rate and 46% at 150C, with more than 1000 charge/discharge cycles. These results demonstrate the promise of ultrasmall, self-aligned/regular, densely packed nanobattery structures as a testbed to study ionics and electrodics at the nanoscale with various geometrical modifications and as a building block for high performance energy storage systems[1, 2]. Further increase of full cell output potential is also demonstrated in asymmetric full cell configurations with various low voltage anode materials. The asymmetric full cell nanopore batteries, comprised of V2O5 as cathode and prelithiated SnO2 or anatase phase TiO2 as anode, with integrated nanotubular metal current collectors underneath each nanotubular storage electrode, also enabled by ALD. By controlling the amount of lithium ion prelithiated into SnO2 anode, we can tune full cell output voltage in the range of 0.3V and 3V. This asymmetric nanopore battery array displays exceptional rate performance and cyclability. When cycled between 1V and 3V, it has capacity retention of approximately 73% at 200C rate compared to 1C, with only 2% capacity loss after more than 500 charge/discharge cycles. With increased full cell output potential, the asymmetric V2O5-SnO2 nanopore battery shows significantly improved energy and power density. This configuration presents a more realistic test - through its asymmetric (vs symmetric) configuration – of performance and cyclability in nanoconfined environment. This dissertation covers (1) Ultra small electrochemical storage platform design and fabrication, (2) Electron and ion transport in nanostructured electrodes inside a half cell configuration, (3) Ion transport between anode and cathode in confined nanochannels in symmetric full cells, (4) Scale up energy and power density with geometry optimization and low voltage anode materials in asymmetric full cell configurations. As a supplement, selective growth of ALD to improve graphene conductance will also be discussed[3]. References: 1. Liu, C., et al., (Invited) A Rational Design for Batteries at Nanoscale by Atomic Layer Deposition. ECS Transactions, 2015. 69(7): p. 23-30. 2. Liu, C.Y., et al., An all-in-one nanopore battery array. Nature Nanotechnology, 2014. 9(12): p. 1031-1039. 3. Liu, C., et al., Improving Graphene Conductivity through Selective Atomic Layer Deposition. ECS Transactions, 2015. 69(7): p. 133-138.

Detailed analysis of a conservative difference approximation for the time fractional diffusion equation

Diffusion equations that use time fractional derivatives are attractive because they describe a wealth of problems involving non-Markovian Random walks. The time fractional diffusion equation (TFDE) is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α ∈ (0, 1). Developing numerical methods for solving fractional partial differential equations is a new research field and the theoretical analysis of the numerical methods associated with them is not fully developed. In this paper an explicit conservative difference approximation (ECDA) for TFDE is proposed. We give a detailed analysis for this ECDA and generate discrete models of random walk suitable for simulating random variables whose spatial probability density evolves in time according to this fractional diffusion equation. The stability and convergence of the ECDA for TFDE in a bounded domain are discussed. Finally, some numerical examples are presented to show the application of the present technique.