Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time factors”

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 1,009 records · Page 56

Artificial light at night: an underappreciated effect on phenology of deciduous woody plants

Artificial light at night (ALAN), an increasing anthropogenic driver, is widespread and shows rapid expansion with potential adverse impact on the terrestrial ecosystem. However, whether and to what extent does ALAN affect plant phenology, a critical factor influencing the timing of terrestrial ecosystem processes, remains unexplored due to limited ALAN observation. Here, we used the Black Marble ALAN product and phenology observations from USA National Phenology Network to investigate the impact of ALAN on deciduous woody plants phenology in the conterminous United States. We found that (1) ALAN significantly advanced the date of breaking leaf buds by 8.9 ± 6.9 days (mean ± SD) and delayed the coloring of leaves by 6.0 ± 11.9 days on average; (2) the magnitude of phenological changes was significantly correlated with the intensity of ALAN (P < 0.001); and (3) there was an interaction between ALAN and temperature on the coloring of leaves, but not on breaking leaf buds. We further showed that under future climate warming scenarios, ALAN will accelerate the advance in breaking leaf buds but exert a more complex effect on the coloring of leaves. This study suggests intensified ALAN may have far-reaching but underappreciated consequences in disrupting key ecosystem functions and services, which requires an interdisciplinary approach to investigate. Developing lighting strategies that minimize the impact of ALAN on ecosystems, especially those embedded and surrounding major cities, is challenging but must be pursued.

54 ENVIRONMENTAL SCIENCES↗

Frustrated magnetic interactions in FeSe

The structurally simplest high-temperature superconductor FeSe exhibits an intriguing superconducting nematic paramagnetic phase with unusual spin excitation spectra that are different from typical spin waves; thus, determining its effective magnetic exchange interactions is challenging. In this work, we report neutron scattering measurements of spin fluctuations of FeSe in the tetragonal paramagnetic phase. We show that the equal-time magnetic structure factor, S(Q), can be effectively modeled using the self-consistent Gaussian approximation calculation with highly frustrated nearest-neighbor (J 1 ) and next-nearest-neighbor (J 2 ) exchange couplings, and very weak further neighbor exchange interaction. Our results elucidate the frustrated magnetism in FeSe, which provides a natural explanation for the highly tunable superconductivity and nematicity in FeSe and related materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Fast and flexible analysis of direct dark matter search data with machine learning

We present the results from combining machine learning with the profile likelihood fit procedure, using data from the Large Underground Xenon (LUX) dark matter experiment. This approach demonstrates reduction in computation time by a factor of 30 when compared with the previous approach, without loss of performance on real data. We establish its flexibility to capture non-linear correlations between variables (such as smearing in light and charge signals due to position variation) by achieving equal performance using pulse areas with and without position-corrections applied. Its efficiency and scalability furthermore enables searching for dark matter using additional variables without significant computational burden. We demonstrate this by including a light signal pulse shape variable alongside more traditional inputs such as light and charge signal strengths. Furthermore, this technique can be exploited by future dark matter experiments to make use of additional information, reduce computational resources needed for signal searches and simulations, and make inclusion of physical nuisance parameters in fits tractable.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantum Enhanced Cavity QED Interferometer with Partially Delocalized Atoms in Lattices

Herein, we propose a quantum enhanced interferometric protocol for gravimetry and force sensing using cold atoms in an optical lattice supported by a standing-wave cavity. By loading the atoms in partially delocalized Wannier-Stark states, it is possible to cancel the undesirable inhomogeneities arising from the mismatch between the lattice and cavity fields and to generate spin squeezed states via a uniform one-axis twisting model. The quantum enhanced sensitivity of the states is combined with the subsequent application of a compound pulse sequence that allows us to separate atoms by several lattice sites. This, together with the capability to load small atomic clouds in the lattice at micrometric distances from a surface, make our setup ideal for sensing short-range forces. We show that for arrays of 10 4 atoms, our protocol can reduce the required averaging time by a factor of 10 compared to unentangled lattice-based interferometers after accounting for primary sources of decoherence.

74 ATOMIC AND MOLECULAR PHYSICS↗

Entanglement Purification and Protection in a Superconducting Quantum Network

High-fidelity quantum entanglement is a key resource for quantum communication and distributed quantum computing, enabling quantum state teleportation, dense coding, and quantum encryption. Any sources of decoherence in the communication channel, however, degrade entanglement fidelity, thereby increasing the error rates of entangled state protocols. Entanglement purification provides a method to alleviate these nonidealities by distilling impure states into higher-fidelity entangled states. In this work, we demonstrate the entanglement purification of Bell pairs shared between two remote superconducting quantum nodes connected by a moderately lossy, 1-meter long superconducting communication cable. We use a purification process to correct the dominant amplitude damping errors caused by transmission through the cable, with fractional increases in fidelity as large as 25%, achieved for higher damping errors. The best final fidelity the purification achieves is 94.09 ± 0.98%. In addition, we use both dynamical decoupling and Rabi driving to protect the entangled states from local noise, increasing the effective qubit dephasing time by a factor of 4, from 3 to 12 μs. These methods demonstrate the potential for the generation and preservation of very high-fidelity entanglement in a superconducting quantum communication network.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

LaRIS: Targeting Portability and Productivity for LAPACK Codes on Extreme Heterogeneous Systems by Using IRIS

In keeping with the trend of heterogeneity in high-performance computing, hardware manufacturers and vendors are developing new architectures and associated software stacks (e.g., libraries) to harness the best possible performance from commonly used kernels (e.g., linear algebra kernels). However, kernels tuned for one architecture are not portable to others. Moreover, the coexistence of different architectures in a single node makes orchestration difficult. To address these challenges, we introduce LaRIS, a portable framework for LAPACK functionalities. LaRIS ensures a separation between linear algebra algorithms and vendor-library kernels by using the IRIS run time and IRIS-BLAS library. Such abstraction at the algorithm level makes the implementation completely agnostic to the vendor library and architecture. LaRIS uses the IRIS run time to dynamically select the vendor-library kernel and suitable processor architecture at run time. Through LU factorization, we demonstrate that LaRIS can fully utilize different heterogeneous systems by launching and orchestrating different vendor-library kernels without any change in the source code.

Monil, M. A. H.↗

Laminar Burning Speed Measurements of Hydrogen/Natural Gas Mixtures

Abstract To address the shortage of fossil fuels and the environmental impact of exhaust emissions, research efforts in combustion seek to identify candidates for sustainable fuels. Hydrogen is an example, from which water is produced as the main combustion product. However, when used as neat hydrogen, it is very difficult to control hydrogen combustion due to its high flame speed, pre-ignition propensity, and flashback characteristics. One method to control hydrogen combustion is to blend it with a well-established fuel. Natural gas is a common fuel for use in engines or gas turbines. Natural gas is characterized by slow flame speeds and poor lean-burn ability; thus, engine power is decreased in the lean-burn region. Consequently, natural gas engines are typically operated at stoichiometric conditions to maximize output. Blending hydrogen to natural gas has several advantages: (i) Control of hydrogen’s ignition characteristics, (ii) extension of Natural gas’s lower flammability limit — making it easier to ignite, and (iii) reduced greenhouse gas emissions. In this study, the laminar burning speed (LBS) for neat natural gas and natural gas/hydrogen mixtures at 5 atm and 296 K were measured. Schlieren optical imaging was utilized to validate laminar flame conditions and further investigate the structure of natural gas/hydrogen flames. Results indicate that replacing 50% of natural gas with hydrogen increases the peak laminar burning velocity by a factor of ∼1.6 times. LBS simulations were conducted with NUI 1.1 and UCF NG/H2 mechanisms; Predictions by both mechanisms were found to be satisfactory for neat NG mixtures. However, for 50% NG/ 50% H2 mixture, the true LBS for fuel-lean conditions was under-predicted by both models, while predictions for fuel-rich conditions were satisfactory. Additionally, a literature review and similar experimental analysis will be performed for methane, ammonia, and hydrogen fuel blends to better understand the combustion phenomena of natural gas fuels.

Yovino, Louis↗

Electron beam injection from a hollow cathode plasma into a downstream reactive environment: Characterization of secondary plasma production and Si 3 N 4 and Si etching

A material etching system was developed by combining beam electron injection from a direct current (DC) hollow-cathode (HC) electron source with the downstream reactive environment of a remote CF 4 /O 2 low temperature plasma. The energy of the injected beam electrons is controlled using an acceleration electrode biased positively relative to the HC Argon discharge. For an acceleration voltage greater than the ionization potential of Ar, the extracted primary electrons can produce a secondary plasma in the process chamber. We characterized the properties of the secondary plasma by performing Langmuir probe measurements of the electron energy probability function (EEPF) 2.5 cm below the extraction ring. The data indicate the existence of two major groups of electrons, including electrons with a primary beam electron energy which varies as the acceleration voltage is varied, along with low energy electrons produced by ionization of the Ar gas atoms in the process chamber by the injected beam electrons. When combining the HC Ar beam electron with a remote CF 4 /O 2 electron cyclotron wave resonance (ECWR) plasma, the EEPF of both the low energy plasma electron and beam electron components decreases. Additionally, we studied surface etching of Si3N4 and poly-crystalline Si (poly-Si) thin films as a function of process parameters, including the acceleration voltage (0-70 V), discharge current of the HC discharge (1-2 A), pressure (2- 100 mTorr), source to substrate distance (2.5 to 5 cm), and feed gas composition (with or without CF 4 /O 2 ). The direction of the incident beam electrons was perpendicular to the surface. Si 3 N 4 and poly-crystalline silicon etching are seen and indicate an electronneutral synergy effect. No/little remote plasma spontaneous etching was observed for the conditions used in this study, and the etching is confined to the substrate area irradiated by the injected beam electrons. The electron etched Si 3 N 4 surface etching rate profile distribution is confined within a ~30 mm diameter circle, which is slightly broader than the area for which poly-Si etching is seen, and coincides closely with the spatial profile of beam electrons as determined by the Langmuir probe measurements. The magnitude of poly-Si etching rate is by a factor of two times smaller than the Si 3 N 4 etching rate. In this paper, we discuss possible explanations of the data and the role of surface charging.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling Isotope Separation in Electrochemical Lithium Deposition

Naturally occurring Li consists of two stable isotopes, 6 Li with an abundance of about 7.5%, and 7 Li making up the remainder with 92.5%. The development of a 6 Li enrichment technique, in terms of technical reliability and environmental safety to reach 6 Li future requirements, represents a key step in the roadmap for nuclear fusion energy supply worldwide. This paper uses finite element analysis-based models to simulate electrochemical Li isotope separation, which is an attractive method in terms of simplicity, safety, and scalability. In the model, we quantitatively analyze how different electrochemical factors including thermodynamics, charge-transfer kinetics, and diffusivities affect the separation process (separation factor), together with cell parameters, such as cell length and current density. The maximum separation factor of 1.128 could be obtained with the cell under the optimal thermodynamic, kinetic, and diffusive conditions, which is among the highest separation factors ever reported. Furthermore, these results will assist in designing the actual isotope separation setup with large separation factor and appropriate timing for sample collection.

07 ISOTOPE AND RADIATION SOURCES↗

Enhanced Carbon Dioxide Uptake in Drilled Hollow Core Fibers for Raman Spectroscopy

Fiber-based Raman spectroscopy enhances the Raman signal by maximizing the overlap of the optical field and the gas species. However, filling the hollow-core fiber (HCF) with gas requires time that is dependent on the fiber core diameter, fiber length, and pressure of the gas. At ambient pressure, the fiber gas uptake is driven by diffusion into the fiber ends, severely limiting the response time of the system. By laser drilling access holes to the core along the length of the fiber, the uptake time of the gas is reduced, improving the system response time. In this work, we study the carbon dioxide (CO 2 ) sensor dynamics based on Raman signal intensity generated in HCFs. The signal intensity versus gas concentration is characterized by controlling the CO 2 concentration in the surrounding environment of the fiber. Next, we characterize the gas uptake time in HCFs as a function of fiber length. Finally, we optimize the access hole configuration along the fiber, demonstrating reduced sensor uptake time by a factor of three.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparing Designed Training Sets to Optimize Multivariate Regression Models for Pr, Nd, and Nitric Acid Using Spectrophotometry

Chemometric regression models were developed for the quantification of praseodymium (Pr, 0–1000 µg/mL), neodymium (Nd, 0–1000 µg/mL), and nitric acid (HNO 3 , 0.1–5 M) using spectrophotometry. Designed calibration sets were composed of 20 samples each: 10 model points and 10 lack-of-fit (LOF) points. The D-optimal designs effectively minimized the number of samples required to build models, and each design resulted in similar prediction performance, suggesting that statistical design of experiments can provide a reliable framework for selecting training set samples in three-variable systems. Partial least squares regression (PLSR) models were validated against a one-factor-at-a-time validation set composed of 125 samples (three variables, five levels). The top PLS-1 models resulted in average percent root mean square error of prediction error values of 3.5%, 1.7%, and 1.2% for Pr(III), Nd(III), and HNO 3 , respectively. Power set augmentations of the model and LOF samples were investigated to optimize the number of training set samples. PLSR models built using just required model points (10) had similar predictive capabilities as models including the LOF points (20) but with fewer samples. The number of validation samples was also varied systematically to learn how many samples are needed to validate regression models. This work addresses long-standing questions in the field of chemometrics to help make this approach amenable to the near-real-time quantification of hazardous species in remote settings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Surrogate modeling of Cellular-Potts agent-based models as a segmentation task using the U-Net neural network architecture

The Cellular-Potts model is a powerful and ubiquitous framework for developing computational models for simulating complex multicellular biological systems. Cellular-Potts models (CPMs) are often computationally expensive due to the explicit modeling of interactions among large numbers of individual model agents and diffusive fields described by partial differential equations (PDEs). In this work, we develop a convolutional neural network (CNN) surrogate model using a U-Net architecture that accounts for periodic boundary conditions. We use this model to accelerate the evaluation of a mechanistic CPM previously used to investigate in vitro vasculogenesis. The surrogate model was trained to predict 100 computational steps ahead (Monte-Carlo steps, MCS), accelerating simulation evaluations by a factor of 562 times compared to single-core CPM code execution on CPU. Over short timescales of up to 3 recursive evaluations, or 300 MCS, our model captures the emergent behaviors demonstrated by the original Cellular-Potts model such as vessel sprouting, extension and anastomosis, and contraction of vascular lacunae. This approach demonstrates the potential for deep learning to serve as a step toward efficient surrogate models for CPM simulations, enabling faster evaluation of computationally expensive CPM simulations of biological processes.

97 MATHEMATICS AND COMPUTING↗

Advanced Solar and Load Forecasting Incorporating HD Sky Imaging (Phase III)

Due to rapidly changing sky conditions, the available solar irradiance for energy generation is subject to wide swings in amplitude. As the market penetration of solar energy continues to increase rapidly, these variations in solar power generation are beginning to have an impact on grid stability and increased wear on power switches. Solar power forecasting plays a critical role in operations for Independent System Operators (ISOs) and utilities. Accurate forecasts help maintain grid reliability, optimize generation from renewables, and reduce operating costs. Of particular interest to the ISOs and utilities are sudden changes in solar irradiance, termed “ramp events,” due to the movement of clouds. One significant impact of ramp events on the grid is additional ancillary service requirements necessary to manage such variability. Ramp events can also cause voltage fluctuations in the distribution grids and trigger actions of automated line equipment (e.g., tap changers), leading to additional maintenance costs. In high penetration solar regions, forecasts must be made for both transmission-connected and distribution-connected resources – either behind the meter or in front of it. Particularly for distributed solar energy resources, forecasting can be a challenge due to the lack of visibility of the energy resource. Brookhaven National Laboratory (BNL) has been working towards nowcasting Global Horizontal Irradiance (GHI) using low-cost technologies for several years now. The Solar NowCasting technology being developed by BNL is a 0 – 30 min solar “nowcasting” technology applicable to a scale covering both large generating facilities and residential, distributed solar installations that relies on a network of ground-based high-definition (HD) cameras and surface pyranometers. GHI forecasts are being produced for both regions using 8 ground-based HD cameras and at least 2 surface pyranometers. When stitched together, the camera images can be used to forecast the impact of clouds on available solar irradiance over a domain of ~50 km 2 . Phase I of the project was the engineering scale up conceptual design stage. Phase II was the initial field test and demonstration, in which the technology was scaled up by a factor of 20 times and successfully demonstrated in eastern Long Island. In Phase III, an additional forecasting network was added in upstate NY and both networks are currently being operated until at least a full year of data is gathered by the Albany network to allow for the collection of sufficient data to evaluate performance against the persistence and smart persistence models using the Solar Arbiter.

14 SOLAR ENERGY↗

304L Can Crush Validation Studies

Accurate prediction of ductile behavior of structural alloys up to and including failure is essential in component or system failure assessment, which is necessary for nuclear weapons alteration and life extensions programs of Sandia National Laboratories. Modeling such behavior requires computational capabilities to robustly capture strong nonlinearities (geometric and material), rate- dependent and temperature-dependent properties, and ductile failure mechanisms. This study's objective is to validate numerical simulations of a high-deformation crush of a stainless steel can. The process consists of identifying a suitable can geometry and loading conditions, conducting the laboratory testing, developing a high-quality Sierra/SM simulation, and then drawing comparisons between model and measurement to assess the fitness of the simulation in regards to material model (plasticity), finite element model construction, and failure model. Following previous material model calibration, a J 2 plasticity model with a microstructural BCJ failure model is employed to model the test specimen made of 304L stainless steel. Simulated results are verified and validated through mesh and mass-scaling convergence studies, parameter sensitivity studies, and a comparison to experimental data. The converged mesh and degree of mass-scaling are the mesh discretization with 140,372 elements, and a mass scaling with a target time increment of 1.0e-6 seconds and time step scale factor of 0.5, respectively. Results from the coupled thermal-mechanical explicit dynamic analysis are comparable to the experimental data. Simulated global force vs displacement (F/D) response predicts key points such as yield, ultimate, and kinks of the experimental F/D response. Furthermore, the final deformed shape of the can and field data predicted from the analysis are similar to that of the deformed can, as measured by 3D optical CMM scans and DIC data from the experiment.

36 MATERIALS SCIENCE↗

(U) A General-Purpose Code for Correlated Sampling Using Batch Statistics with MCNP6 for Fixed-Source Problems

Correlated sampling can be used to reduce the uncertainty of a difference of tallies by taking advantage of the negative covariance term in the sandwich formula. Booth first showed how correlated sampling can be applied with batch statistics using MCNP’s tally fluctuation chart (TFC) to reduce the uncertainty of a difference of tallies in fixed-source problems. Booth presented a problem in which a 1273% uncertainty in a difference was reduced to 8% by accounting for correlations. Researchers He and Su recently studied correlated sampling using the TFC in MCNP version 5. They determined that the code did not print enough digits in the TFC tally means for accurate batch statistics in some cases. After modifying the source code, they concluded that “correlated sampling can yield a standard deviation of about one magnitude smaller than that predicted by the direct, un-correlated simulation when the changes in system response are small (say about 1%), which is equivalent to saving in CPU time by a factor of 100. Such saving [sic] becomes less significant as the change in system response becomes larger.” He and Su provided the formulas needed to apply batch statistics to compute the correlated uncertainty of a difference of tallies. In this report, we follow up on their work by providing the formulas needed to apply batch statistics to compute the correlated uncertainty of a ratio of tallies and of a difference of two tallies divided by a third tally. We extend these formulas to differences and ratios of ratios. These formulas are applied to reduce the uncertainty associated with calculating a relative sensitivity. He and Su did not investigate the accuracy of their correlated sampling uncertainty estimates. We use their test problems and evaluate the accuracy of the uncertainty estimates by comparing with results obtained from random sampling, and, in simple cases, with theoretical values of the “exact” uncertainties. We find that the uncertainties obtained from batch statistics are accurate as long as at least 100 batches are used. We present a new computer code, COSUBS (COrrelated Sampling Using Batch Statistics), that reads MCNP6 TFCs and applies correlated sampling using batch statistics for the tally combinations that the user specifies. COSUBS is a very general tool that compares all TFCs for a base case and one or two perturbed cases. It computes uncertainties for ratios if given only a base case. This report is organized as follows. The equations to apply batch statistics to the difference of random tallies are reviewed in Sec. II. Section III presents the equations for applying batch statistics to a ratio of random tallies; this is useful for computing relative sensitivities using a one-sided finite difference and the relative sensitivity using the differential operator method. Section IV presents the equations for applying batch statistics to a difference of two random tallies divided by a third; this is useful for computing a relative sensitivities using a central difference. Section V presents the equations for applying batch statistics to a difference of two ratios with four random tallies. Section VI presents the equations for applying batch statistics to a one-sided finite difference estimate of the relative sensitivity of a ratio (this uses four random tallies). Section VII presents the equations for applying batch statistics to a central difference estimate of the relative sensitivity of a ratio (this uses six random tallies). Section VIII presents the equations for applying batch statistics to a sum of random tallies. Section IX discusses how to apply batch statistics using MCNP6. Section X presents COSUBS, describing its command-line options and logic. Sections XI through XVI present numerical results for various test problems. Section XVII is a summary and conclusions. Appendix A derives the theoretical Monte Carlo tally variance given certain assumptions; these variances are used to verify the batch statistics for some of the problems. Appendix B lists the MCNP6 input for the unperturbed example problem. Appendix C presents modifications made to MCNP6.3 to support this work.

97 MATHEMATICS AND COMPUTING↗

Relativistically Induced Transparency in Plasma (Final Research Report)

The major goal of this project was to study the process of relativistic transparency in plasmas. This process occurs when a very high intensity laser pulse interacts with a plasma. The electric field of the laser oscillates the electron speeds approaching the speed of light, increasing the effective mass of the electrons by the time-averaged Lorentz factor <γ>. The increase in the effective mass alters the critical density, the density at which the plasma becomes opaque to the electromagnetic wave of the laser. This project used experiments and particle-in-cell simulations to study this process in ultra-thin foil interactions. Experiments were performed through the LaserNetUS program on the Ohio State University laser system Scarlet. The transmitted and reflected light was studied as a function of the target thickness, the near-field profiles, the total energy and pulse duration were measured. Particle-in-cell simulations, using the OSIRIS code were used to model the interaction and better understand the dynamics. The results of both show current relativistic transparency theory is insufficient to predict the results, this is likely because the theories do not account for the laser energy absorption into the target and global plasma fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Exotic and Conventional Quarkonium Physics Prospects at Belle II

The Belle II experiment at the SuperKEKB energy-asymmetric $e^+ e^-$ collider is a substantial upgrade of the B factory facility at KEK in Tsukuba, Japan. It aims to record a factor of 50 times more data than its predecessor. The experiment completed a commissioning run in 2018, and began full operation in early 2019. Belle II is uniquely capable of studying the so-called "XYZ" particles: heavy exotic hadrons consisting of more than three quarks. First discovered by Belle, these now number in the dozens, and represent the emergence of a new category within quantum chromodynamics. This talk will present the prospects of Belle II to explore both exotic and conventional quarkonium physics.

Fulsom, Bryan G.↗

Oxidation of Dissolved Tetravalent Selenium by Birnessite: Se Isotope Fractionation and the Effects of pH and Birnessite Structure

Redox reactions control the mobility and bioavailability of selenium (Se) in biogeochemical systems, both modern and ancient. Se isotope ratio measurements (e.g., 82 Se/ 76 Se) have been developed to enhance understanding of biogeochemical transformations and transport of Se. Stable isotope ratios of many elements are known to be powerful indicators of redox reactions, and shifts in 82 Se/ 76 Se have been observed for Se reduction reactions. However, Se isotope shifts caused by naturally relevant oxidation reactions have not been published. Here, we report Se isotope fractionation factors for oxidation of Se(IV) by birnessite. Experiments were conducted at pH = 4.0 and 5.5, with two types of birnessite of contrasting composition at two concentrations of suspended birnessite. The results are consistent with a single 82 Se/ 76 Se fractionation factor, for all times during all experiments, of 0.99767 (±0.0035 2 s.d.). Expressed as ε, the fractionation is 2.33‰ (±0.08‰).

58 GEOSCIENCES↗