Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid characterization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41

XPCS Microrheology and Rheology of Sterically Stabilized Nanoparticle Dispersions in Aprotic Solvents

In this study, X-ray photon correlation spectroscopy (XPCS) microrheology and conventional bulk rheology were performed on silica nanoparticle dispersions associated with battery electrolyte applications to probe the properties of these specific complex materials and to explore the utility of XPCS microrheology in characterizing nanoparticle dispersions. Sterically stabilized shear-thickening electrolytes were synthesized by grafting poly(methyl methacrylate) chains onto silica nanoparticles. Coated silica dispersions containing 5–30 wt % nanoparticles dispersed in propylene carbonate were studied. In general, both XPCS microrheology and conventional rheology showed that coated silica dispersions were more viscous at higher concentrations, as expected. The complex viscosity of coated silica dispersions showed shear-thinning behavior over the frequency range probed by XPCS measurements. However, measurements using conventional mechanical rheometry yielded a shear viscosity with weak shear-thickening behavior for dispersions with the highest concentration of 30% particles. Our results indicate that there is a critical concentration needed for shear-thickening behavior, as well as appropriate particle size and surface polymer chain length, for this class of nanoparticle-based electrolytes. The results of this study can provide insights for comparing XPCS microrheology and bulk rheology for related complex fluids and whether XPCS microrheology can capture expected macroscopic rheological properties by probing small-scale particle dynamics.

36 MATERIALS SCIENCE↗

Magnetoelastic effects in the hyperhoneycomb Kitaev spin liquid

In recent years, several magnetic Mott insulators with strong spin-orbit coupling were suggested to be proximate to the Kitaev quantum spin liquid, whose one of the most exciting features is the fractionalization of spin excitations into itinerant Majorana fermions and static Z 2 fluxes. Unfortunately, the ground states of these systems cannot be easily captured by experiment, remaining featureless to conventional local probes. Here, we propose to study the signatures of fractionalized excitations by exploiting their coupling to the lattice vibrations, dubbed magnetoelastic coupling, which arises from the fact that the interaction between spins depends on the relative distance between them. We argue that the magnetoelastic coupling can lead to the distinct modification of the phonon dynamics, which can be observed by measuring renormalized phonon spectrum, the sound attenuation, and the phonon Hall viscosity. This makes the phonon dynamics a promising tool for the characterization and identification of quantum spin liquid phases. In this work, we focus on the magnetoelastic effects in the three-dimensional Kitaev model realized on the hyperhoneycomb lattice. Here, the hyperhoneycomb Kitaev spin liquid is particularly interesting since the strong Kitaev interaction was observed in the Kitaev magnet β-Li 2 IrO 3 , for which the spin-orbit entangled J eff = 1/2 moments of iridium ions form precisely the hyperhoneycomb lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Wave topology in Hall magnetohydrodynamics

Hall magnetohydrodynamics (HMHD) extends ideal MHD by incorporating the Hall effect via the induction equation, making it more accurate for describing plasma behavior at length scales below the ion skin depth. Despite its importance, a comprehensive description of the eigenmodes in HMHD has been lacking. In this work, we derive the complete spectrum and eigenvectors of HMHD waves and identify their underlying topological structure. We prove that the HMHD wave spectrum is homotopic to that of ideal MHD, consisting of three distinct branches: the slow magnetosonic-Hall waves, the shear Alfvén-Hall waves, and the fast magnetosonic-Hall waves, which continuously reduce to their ideal MHD counterparts in the limit of vanishing Hall parameter. Contrary to a recent claim [Mahajan, Sharma, and Lingam, Phys. Plasmas 31, 090701 (2024)], we find that HMHD does not admit any additional wave branches beyond those in ideal MHD. In conclusion, the key qualitative difference lies in the topological nature of the HMHD wave structure: it exhibits nontrivial topology characterized by a Weyl point—an isolated eigenmode degeneracy point—and associated nonzero Chern numbers of the eigenmode bundles over a 2-sphere in 𝐤-space surrounding the Weyl point.

Alfvén waves↗

Seaweed-derived fucoidans and rhamnan sulfates serve as potent anti-SARS-CoV-2 agents with potential for prophylaxis

Seaweeds represent a rich source of sulfated polysaccharides with similarity to heparan sulfate, a facilitator of myriad virus host cell attachment. For this reason, attention has been drawn to their antiviral activity, including the potential for anti-SARS-CoV-2 activity. We have identified and structurally characterized several fucoidan extracts, including those from different species of brown macroalga, and a rhamnan sulfate from a green macroalga species. A high molecular weight fucoidan extracted from Saccharina japonica (FSjRPI-27), and a rhamnan sulfate extracted from Monostroma nitidum (RSMn), showed potent competitive inhibition of spike glycoprotein receptor binding to a heparin-coated SPR chip. This inhibition was also observed in cell-based assays using hACE2 HEK-293 T cells infected by pseudotyped SARS-CoV-2 virus with IC50 values <1 μg/mL. Effectiveness was demonstrated in vivo using hACE2–transgenic mice. Intranasal administration of FSjRPI-27 showed protection when dosed 6 h prior to and at infection, and then every 2 days post-infection, with 100 % survival and no toxicity at 104 plaque-forming units per mouse vs. buffer control. At 5-fold higher virus dose, FSjRPI-27 reduced mortality and yielded reduced viral titers in bronchioalveolar fluid and lung homogenates vs. buffer control. Furthermore, these findings suggest the potential application of seaweed-based sulfated polysaccharides as promising anti-SARS-CoV-2 prophylactics.

60 APPLIED LIFE SCIENCES↗

Multifunctional Liquid Marble Compound Lenses

Mosquito compound eyes are elaborate multifunctional hierarchical structures. The presence of ordered curved features spanning length scales of nanometers to millimeters provides the mosquito eye with a wide field of view, an infinite depth of field, as well as antifogging properties. Developing bio-inspired compound lenses is challenging because of the need to mimic all characteristic curvatures along with their functionalities. Herein, we show how the curvature inherent to nanoparticles, emulsion droplets, and liquid marbles can be employed to mimic the hierarchical structure and functionality of mosquito compound eyes. At the nanometers to micron length scale we employ nanoparticle-stabilized emulsion droplets of photocurable oil to form microlenses with nanoscale surface features. After polymerization, the microlenses form a monolayer on an oil droplet to create an optically clear, millimeter scale, liquid marble that functions as a compound lens. Here, we characterize the optical and surface properties of the compound lenses and find that they reproduce the functionality of the mosquito eye. Furthermore, we exploit the mobility and reconfigurability of liquid marbles to create arrays (centimeter scale) of compound lenses and other types of functional lenses such as the Janus lens that magnifies the image acquired by the compound lens. Simple and scalable methods to create compound lenses could aid in the development of miniaturized advanced vision systems.

36 MATERIALS SCIENCE↗

X-ray imaging and radiation transport effects on cylindrical implosions

Magnetization of inertial confinement implosions is a promising means of improving their performance, owing to the potential reduction of energy losses within the target and mitigation of hydrodynamic instabilities. In particular, cylindrical implosions are useful for studying the influence of a magnetic field, thanks to their axial symmetry. In this work, we present experimental results from cylindrical implosions on the OMEGA-60 laser using a 40-beam, 14.5 kJ, 1.5 ns drive and an initial seed magnetic field of B 0 = 30 T along the axes of the targets, compared with reference results without an imposed B-field. Implosions were characterized using time-resolved x-ray imaging from two orthogonal lines of sight. We found that the data agree well with magnetohydrodynamic simulations, once radiation transport within the imploding plasma is considered. We show that for a correct interpretation of the data in these types of experiments, explicit radiation transport must be taken into account.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effective Permeability of a Nuclear Fuel Assembly

This report aids in the development of models to perform characterization studies of aerosol dispersal and deposition within a spent fuel cask system. Due to the complex geometry in a spent-fuel canister, direct simulation of buoyancy-driven flow through the fuel assemblies to model aerosol deposition within the fuel canister is computationally expensive. Identification of an effective permeability as given in this work for a nuclear fuel assembly greatly simplifies the requirements for thermal hydraulic computations. The results of computations performed using OpenFOAM® to solve the Navier-Stokes Equations for laminar flow are used to determine an effective permeability by applying Darcy's Law. The computations are validated against an analytical solution for the special case of an infinite array of pins for which the numerical and analytical solutions have excellent agreement. The effective permeability of a 1717 PWR nuclear fuel assembly in a basket without spacer grids is numerically determined to be 1.85010 -6 m 2 for the range of fluid viscosities and pressure drops expected in a spent fuel storage canister. However, the flow is not uniform on the scale of multiple pins. Instead, significantly higher velocities are attained in the space between the assembly and the basket walls compared to the flow between the fuel pins within the assembly. Comparison with an analytical solution for fully developed flow through an infinite array of pins shows that the larger spacing near the basket walls results in about a 20% larger permeability compared to the analytical solution which does not include the enhanced flow in the space between the assembly and basket wall, or entrance and exit effects. A preliminary assessment of turbulence effects shows that with a k-epsilon model, significantly higher flow velocities are attained between the fuel pins within the assembly compared to the flow velocity in the space between the assembly and the basket walls. This is the opposite of what is determined for laminar flow.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

FY2022 Progress on Imbibition Testing in Containment Science

Estimation of two-phase fluid flow properties is important to understand and predict water and gas movement through the vadose zone for agricultural, hydrogeological, and engineering applications, such as for vapor-phase contaminant transport and/or containment of noble gases in the subsurface. In this second progress report of FY22, we present two ongoing activities related to imbibition testing on volcanic rock samples. We present the development of a new analytical solution predicting the temperature response observed during imbibition into dry samples, as discussed in our previous first progress report for FY22. We also illustrate the use of a multi-modal capillary pressure distribution to simulate both early- and late-time imbibition data collected on tuff core that can exhibit multiple pore types. These FY22 imbibition tests were conducted for an extended period (i.e., far beyond the time required for the wetting front to reach the top of the sample), which is necessary for parameter estimation and characterization of two different pore types within the samples.

42 ENGINEERING↗

Direct observation of a magnetic-field-induced Wigner crystal

Wigner predicted that when the Coulomb interactions between electrons become much stronger than their kinetic energy, electrons crystallize into a closely packed lattice. A variety of two-dimensional systems have shown evidence for Wigner crystals (WCs). However, a spontaneously formed classical or quantum WC has never been directly visualized. Neither the identification of the WC symmetry nor direct investigation of its melting has been accomplished. Here we use high-resolution scanning tunnelling microscopy measurements to directly image a magnetic-field-induced electron WC in Bernal-stacked bilayer graphene and examine its structural properties as a function of electron density, magnetic field and temperature. At high fields and the lowest temperature, we observe a triangular lattice electron WC in the lowest Landau level. The WC possesses the expected lattice constant and is robust between filling factor ν ≈ 0.13 and ν ≈ 0.38 except near fillings where it competes with fractional quantum Hall states. Increasing the density or temperature results in the melting of the WC into a liquid phase that is isotropic but has a modulated structure characterized by the Bragg wavevector of the WC. At low magnetic fields, the WC unexpectedly transitions into an anisotropic stripe phase, which has been commonly anticipated to form in higher Landau levels. Furthermore, analysis of individual lattice sites shows signatures that may be related to the quantum zero-point motion of electrons in the WC lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

High-Speed and High-Quality Field Welding Repair Based on Advanced Non-Destructive Evaluation and Numerical Modeling

Creep strength-enhanced ferritic (CSEF) steels such as Grade 91 (9Cr-1Mo-V) and Grade 92 (Fe-9Cr-2W-0.5Mo) steels are widely used in the fossil-fuel-fired and nuclear power plants. The weld integrity of these steels is crucial for power plants' safe and reliable operations. Due to harsh service conditions, the steel weld can become susceptible to environmental degradation. Field welding repair is used to restore the degraded weld’s performance where a controlled temper-bead welding technique is commonly used to temper the freshly formed martensite during welding. However, knowledge of weld repairability is limited and experimental trial and error optimization to achieve desired microstructure and joint properties is expensive and time-consuming. Many existing computational models, e.g., finite element models, are limited to solving heat conduction equation and ignoring convective heat transfer due to molten metal flow. These models can result in over-prediction of peak temperatures of weld pool and heat-affected zone (HAZ), which in turn can affect the accuracy of tempering prediction. Moreover, these finite element models require an input of the deposit profiles in advance and thus limits the usability of these models. Here, a molten pool-based, multi-pass multi-layer model has been developed based on computational fluid dynamics (CFD) approach with the Volume of Fluid (VOF) method. The model calculates the bead formation, thereby eliminating the need for pre-determined bead profiles required by finite element models. For computational efficiency, a coordinate system attached to the moving heat source is utilized. A subroutine is developed to convert the temperature profiles in the reference frame stationary to the heat source to that stationary to the workpiece. The converted thermal cycles are then imported into a microstructure model to compute the tempering kinetics and resultant hardness using a Johnson-Mehl-Avrami-Kolmogorov (JMAK), and modified Grange-Baughman parameter. The modeling approach is first developed and validated on single- and multi-pass deposition of stainless steel filler metal onto a SA-533 high strength steel substrate. The models are then applied to a multi-pass V-groove repair weld of Grade 91 steel plate as well as directed energy deposition of Grade 92 steel. Non-destructive characterization of microstructures was performed on Grade 91 and 92 steel welds. Two welding processes, cold metal transfer (CMT) and flux-cored arc welding (FCAW), were investigated for the Grade 91 steel weld samples. For the Grade 92 weld samples, three different heat inputs (low, medium, and high) of gas tungsten arc welding (GTAW) were utilized to replicate traditional field welding processes. The non-destructive evaluation (NDE) method used for this research was immersion ultrasonic testing (UT) using a micro-resolution ultrasonic imaging methodology specifically designed to operate in the through-transmission configuration operating at 20 MHz of frequency. The system used a focused ultrasonic beam spot size diameter between 250-300 μm, and a 6 μm laser vibrometer spot size for detection, to produce highly defined images with longitudinal and mode-converted shear waves. From the micro-resolution ultrasonic C-scan images, three microstructural regions, i.e., weld metal (WM), HAZ, and base metal (BM), were clearly identifiable. Various levels of ultrasonic amplitudes distributed over the three regions were correlated with electron beam backscattered diffraction (EBSD) images using grain size, grain boundaries, and dislocation densities. The results showed that areas with relatively higher ultrasonic amplitude levels were associated with smaller grains and higher dislocation densities, while areas with lower amplitude levels were associated with larger grains and lower dislocation densities. In addition, ultrasonic velocity data obtained across the three different weld microstructural regions of Grade 91 test samples were correlated with optical metallographic images and hardness measurements. The results showed distinctive decreases in ultrasonic velocity and hardness over the HAZ region, where weld failures often occur during service.

36 MATERIALS SCIENCE↗

TEAMER: Experimental performance characterization of a shrouded axial-flow turbine

Sitkana has developed a shrouded hydrokinetic turbine with a modular, low-cost design that can be scaled to meet the needs of remote communities. With technical support from the University of Washington, Sitkana sought to experimentally characterize the mechanical power and structural loads of various 1:3.3 scale rotor geometries. In all, 11 different rotor geometries were characterized with variations in height-to-diameter ratio, blade number, and blade type (foiled versus flat). All tests were conducted in Reynolds-independent flow conditions in the Alice C. Tyler Flume at the University of Washington. Results allow Sitkana to (1) refine the optimal rotor geometry, (2) validate numerical models, and (3) predict power output for a full-scale system. This project is part of the TEAMER RFTS 8 (request for technical support) program.

16 TIDAL AND WAVE POWER↗

Concurrent Optimization of Capital Cost and Expected O&M

Concentrating solar power (CSP) technologies can utilize heat from concentrated sunlight from a field of tracking mirrors to generate electricity, reform fuel, provide process heat, or augment fossil plant heat sources. Electricity-generating power tower systems focus light from thousands of independent heliostats onto a thermal receiver, which uses the focused light to warm a heat transfer fluid (HTF), typically, a molten nitrate salt. The HTF is then sent to a power generation cycle or diverted into thermal energy storage (TES) for later use. Thermal storage is – in principle – a straightforward proposition. However, optimal utilization of a TES resource is complex and multi-faceted: thermal energy may be dispatched to produce electricity immediately upon first availability, or thermal energy may be reserved for next-day peak periods at risk of filling storage and dumping energy, or a portion of the thermal energy can be reserved to maintain equipment temperatures, reducing power cycle startup time, etc. Many possible dispatch permutations variously emphasize producing peak power, operating through transients, expediting daily startup, etc. The best operation strategy can change day-to-day throughout the year, depending on the weather and market pricing forecasts. The project we describe in this report develops a software package that allows users to explore design optimization, operations decisions, and performance characterization of concentrating solar power tower plants. Users interface with the tool through a scripting language, and results are reported in time series tables, plots, runtime logs, and design outputs. Users choose from a list of variables such as tower height, solar multiple, design-point irradiance, thermal storage size, etc., and specify information about the system using a list of parameters. The software can then optimize the specified variables to reduce the cost of energy produced by the system while meeting certain production requirements, accounting for uncertain weather and electricity price forecasts, and correcting for equipment failures or repair time. The software we develop is the first comprehensive design tool of its kind to incorporate all of these aspects while being deployed as open source.

14 SOLAR ENERGY↗

Instantaneous Three-Dimensional Temperature Measurements via Ultrafast Laser Spectroscopy with Structured Light

Detonations and flames are characterized by three-dimensional (3D) temperature fields, yet state-of- the-art temperature measurement techniques yield information at a point or along a line. The goal of the research documented here was to combine ultrafast laser spectroscopy and structured illumination to deliver an unprecedented measurement capability—three-dimensional, instantaneous temperature measurements in a gas-phase volume. To achieve this objective, different parts of the proposed technique were developed and tested independently. Structured illumination was used to image particulate matter (soot) in a turbulent flame at multiple planes using a single laser pulse and a single camera. Emission spectroscopy with structured detection was demonstrated for emission- based measurements of explosives with enhance dimensionality. Finally, an instrument for multi- planar laser-based temperature measurement technique was developed. Structured illumination techniques will continue to be developed for multi-dimensional and multi-parameter measurements. These new measurement capabilities will be important for heat transfer and fluid dynamic research areas.

42 ENGINEERING↗

Gaseous jet through an outward opening injector: Details of mixing characteristic and turbulence scales

Direct injection (DI) strategy of natural gas (NG) into internal combustion engines (ICE) has led to higher thermal efficiency and lower exhaust emissions. In order to thoroughly understand the most relevant phenomena affecting the performances of such engines, computational fluid dynamics (CFD) plays a key role as an accurate description of the jet evolution and interaction within the combustion chamber is required to that aim. Accurate description of high-pressure gaseous jets is rather challenging at high Mach numbers, as the injected gas is strongly under-expanded once in the ambient, giving room to shocks due to compressibility effects. Also the interaction between shock waves and mixing layers needs to be carefully represented with a multi-dimensional model, calling for substantial computational resources requirements. In this paper a numerical investigation of the behavior of a gaseous jet (Argon) through an outward opening injector has been carried out. A Large Eddy Simulation (LES) approach has been used in order to track the structures derived by the interaction of the injected fuel with the surrounding ambient. Although already good results were obtained using a Reynolds Averaged Navier-Stokes (RANS) approach, the adoption of LES is required to characterize more accurately the jet properties in terms of vortex structures and mixing effectiveness. Finally, the effect of the Nozzle Pressure Ratio (NPR) on the jet evolution has been highlighted in the paper, showing how a higher NPR would give a faster injection process, compromising however the homogeneity of the mixture.

33 ADVANCED PROPULSION SYSTEMS↗

Characterization of multi-regime reaction zones in a piloted inhomogeneous jet flame with local extinction

In this study, gradient free regime identification (GFRI) is applied to 1D Raman/Rayleigh/LIF measurements of temperature and major species from the intermediate velocity case of the Sydney piloted inhomogeneous jet flame series to better understand the structure of reaction zones and the downstream evolution of multi-regime characteristics. The GFRI approach allows local reaction zones to be detected and characterized as premixed, dominantly premixed, multi-regime, dominantly non-premixed, or non-premixed flame structures, based on flame markers (mixture fraction, chemical mode, and heat release rate) derived from the experimental data. The statistics of chemical mode zero-crossings, which mark premixed reaction zones, and the relative populations of flame structures are shown to be sensitive to the state of mixing in the near field of the flame and to the level of local extinction farther downstream. Multi-regime structures, where premixed and non-premixed reaction zones occur in close proximity and both contribute to overall heat release, account for nearly half the total population at streamwise locations within the first several jet diameters. There is a rapid transition within the near field whereby the relative population of non-premixed and dominantly non-premixed structures grows from 0.05 to nearly 0.5, and the population of premixed and dominantly premixed structures decreases correspondingly as fluid entering the reaction zone becomes progressively fuel-rich. Local extinction and re-ignition bring a resurgence in premixed-type structures, many of which occur at fuel-lean conditions. There are also modest populations of multi-regime structures, having chemical mode zero-crossings at lean conditions, which would not exist in a fully burning jet flame.

42 ENGINEERING↗

Machine Learning to Identify Geologic Factors Associated with Production in Geothermal Fields: A Case-Study Using 3D Geologic Data from Brady Geothermal Field and NMFk

In this paper, we present an analysis using unsupervised machine learning (ML) to identify the key geologic factors that contribute to the geothermal production in Brady geothermal field. Brady is a hydrothermal system in northwestern Nevada that supports both electricity production and direct use of hydrothermal fluids. Transmissive fuid-fow pathways are relatively rare in the subsurface, but are critical components of hydrothermal systems like Brady and many other types of fuid-fow systems in fractured rock. Here, we analyze geologic data with ML methods to unravel the local geologic controls on these pathways. The ML method, non-negative matrix factorization with k-means clustering (NMFk), is applied to a library of 14 3D geologic characteristics hypothesized to control hydrothermal circulation in the Brady geothermal field. Our results indicate that macro-scale faults and a local step-over in the fault system preferentially occur along production wells when compared to injection wells and non-productive wells. We infer that these are the key geologic characteristics that control the through-going hydrothermal transmission pathways at Brady. Our results demonstrate: (1) the specific geologic controls on the Brady hydrothermal system and (2) the efficacy of pairing ML techniques with 3D geologic characterization to enhance the understanding of subsurface processes. This submission includes the published journal article detailing this work, the published 3D geologic map of the Brady Geothermal Area used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity, 3D well data, along which geologic data were sampled for PCA analyses, and associated metadata file. This work was done using the GeoThermalCloud framework, which is part of SmartTensors (both are linked below).

15 GEOTHERMAL ENERGY↗

Influence of microstructure on replacement and porosity generation during experimental dolomitization of limestones

Replacement reactions commonly alter the multiscale pore structures of rocks during fluid-rock interactions. Analysis of these processes in various model fluid-rock systems during controlled laboratory experiments provides insights into the origins of microstructures found in natural materials. This study focused on understanding the effects of initial starting material permeability and resultant differences in transport pathways on porosity and mineralogical changes during limestone dolomitization.A series of replacement experiments (32–317 days in duration) have been conducted in which 1.59 cm (5/8 in.) diameter cores of two different limestones were reacted with saturated MgCl2 solutions at 200 °C. The Texas Cream (Austin Chalk) is a high-porosity, high-permeability limestone, whereas both the porosity and permeability of the Carthage Marble (Burlington Limestone) are relatively low. Altered limestones were imaged using scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM-EDX), Time-of-Flight Secondary Ion Mass Spectrometry (ToF–SIMS) and electron microprobe analysis (EMPA). A representative grain boundary of the low-porosity limestone was targeted for a focused ion beam (FIB) lift-out and characterized using transmission electron microscopy (TEM). These results were coupled with analyses of radial changes in the porosity distribution of the core derived from X-ray and neutron small- and ultra-small angle scattering ((U)SANS/(U)SAXS).The high-porosity/permeability limestone showed a four times faster bulk replacement rate than the lower-porosity/permeability material, and a different mechanism of porosity development. For the low-porosity limestone, a two-stage replacement occurred, with the reacted region of the core consisting of an inner rim in which the limestone was replaced by two calcite-dolomite solid solutions, and an outer rim in which the dolomite was replaced by magnesite. Elongated pores formed along grain boundaries at the initial limestone/dolomite reaction interface, and additional nanometer-scale porosity was formed at the secondary magnesite replacement rim. Grain boundaries were identified as preferential pathways for transport leading to dolomitization and a grain boundary diffusion rate was calculated based on microstructural characterization. In contrast, replacement in the high-porosity limestone was accompanied by porosity generation through replacement of individual grains by dolomite throughout the sample and, in longer runs, magnesite in outer parts of the core. These observations emphasize that both the mechanisms of the replacement reaction and the microstructure and chemistry of the replaced product are contingent on the initial structure of the starting material.

36 MATERIALS SCIENCE↗

Oxidation state and coordination environment of Pb in U-bearing minerals

This study examines the oxidation and coordination of lead (Pb) in the uranium (U) minerals uraninite (UO 2 ), coffinite (USiO 4 ), and brannerite (UTi 2 O 6 ), and minerals that have U as a minor constituent, zircon (ZrSiO 4 ) and titanite (CaTiSiO 5 ). The characterization of Pb was conducted through a combination of micro-synchrotron X-ray fluorescence mapping (µSXRF) and micro-X-ray absorption spectroscopy (µXAS), collected at the Pb L II and L III edge on characterized grains (electron microscopy and X-ray diffraction) from a range of geological settings and ages. For all minerals considered in this study, Pb L II and L III X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) data indicate that Pb exists as Pb 2+ . In the case of uraninite, DFT structure optimization calculations suggest that the substitution of Pb 2+ for U 4+ is facilitated by the stereochemically active electron lone pair indicative of Pb 2+ , which effectively takes the place of an oxygen vacancy (coupled substitution Pb 2+ + = U 4+ + O 2- ). In detail, however, the EXAFS and XANES data, supported by first principle calculations of the XANES associated with Pb 2+ replacing U 4+ in the crystal lattice of uraninite, indicate that Pb 2+ is present in a less distorted coordination than suggested by DFT, which may reflect a yet to be characterized link between interstitial oxygen defects and Pb 2+ incorporation, rather than a simple oxygen vacancy. For coffinite, the XANES data demonstrate that Pb exists mainly as Pb 2+ in the form of galena, indicating limited stability of the isomorphous substitution in coffinite. In contrast, Pb is homogeneously distributed in the metamict mineral brannerite, which displays a distinct XANES spectrum suggestive of Pb 2+ bonded to oxygens in a highly distorted coordination environment. Pb L II XANES indicate that radiogenic Pb exists predominantly as Pb 2+ within the crystal lattice of zircon and as distinct mineral PbO nano-scale inclusions or domains. Additionally, Pb L II XANES measurements also indicate that Pb 2+ exists within titanite although the coordination environment is unclear. Overall, these data suggest that the previous indications of the presence of Pb 4+ in radioactive minerals need to be reconsidered. The nature of radiogenic Pb 2+ incorporation in U-bearing minerals controls the mobilization of Pb in hydrothermal fluids and the utility of a mineral as a geochronometer. The Pb valence and coordination data presented in this study provide fundamental constraints imperative for our understanding of the controls governing the mobilization of radiogenic Pb during water-rock interaction and metamorphism.

58 GEOSCIENCES↗