Search NASA⌕ Search

SEARCH · Search NASA

Results for “Element removal”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The continuum and leading twist limits of parton distribution functions in lattice QCD

In this study, we present continuum limit results for the unpolarized parton distribution function of the nucleon computed in lattice QCD. This study is the first continuum limit using the pseudo-PDF approach with Short Distance Factorization for factorizing lattice QCD calculable matrix elements. Our findings are also compared with the pertinent phenomenological determinations. Inter alia, we are employing the summation Generalized Eigenvalue Problem (sGEVP) technique in order to optimize our control over the excited state contamination which can be one of the most serious systematic errors in this type of calculations. A crucial novel ingredient of our analysis is the parameterization of systematic errors using Jacobi polynomials to characterize and remove both lattice spacing and higher twist contaminations, as well as the leading twist distribution. This method can be expanded in further studies to remove all other systematic errors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Insulation Activation and Contamination in Low-Power Experiments

Among the diverse advanced nuclear reactor concepts being developed are designs that expect to operate at much higher temperatures than conventional light-water reactors, and therefore require alternative materials for components capable of withstanding these extreme temperatures. Insulation is one such component where alternative high-temperature industrial composites are being considered. Low-power experiments can provide an opportunity for advancing our understanding of material behavior under irradiation; however, these experiments are not without occupational hazards. The focus of this paper is the potential for activation of insulating composites when exposed to neutron fluence during reactor experiments, and to discuss the subsequent contamination of the testing area. This paper will also discuss practical measures for preventing the creation and inhalation of activated particulate matter (i.e., dust), such as engineering and administrative controls and personal protective equipment. This paper evaluates the suitability of several insulation materials for use in reactor experiments based on their activation from neutron irradiation. The materials considered are Pyrogel XTE by Aspen Aerogel, Cerablanket and Kaowool by Morgan Advanced Materials, Thermo-12 Gold by Johns Manville, and Maxsil CF6-2000 by McAllister Mills. Although these materials are not yet widely used in the nuclear industry, they are common in other industries due to their thermal insulating properties. These materials produce dust while handling that can pose health and safety hazards, which could be further exacerbated by activation. The analysis presented in this paper utilizes a deterministic method developed to predict the activation source term of the provided insulation materials. The method calculates the atomic density of each element within the composites and simulates their exposure to a specified neutron flux density representative of low power experimental conditions. During the irradiation phase and throughout the extensive cooling period, neutron absorptions and nuclear decays are tracked for relevant isotopes. The residual activity within each composite is then assessed to identify the material with the lowest potential radiological hazard to workers. Dust will be created during installation and removal of insulating materials and may also be generated during the experiment. Contamination of surrounding surfaces due to dispersion is highly likely unless preventative measures are taken. The results of this analysis are used to estimate an internal dose from inhalation of activated dust, which is significantly lower than the allowable whole-body dose. However, the risk of internal exposure may increase if the dust becomes a transport vector for nuclear material. The goals of this paper are to enhance the understanding of the potential activation and contamination hazards from the use of these insulating materials in low-power nuclear experiments and to foster safer work environments through informed decision-making and strategic planning. This work ultimately contributes to a broader understanding of how new non-nuclear materials for advanced reactors impact the radiological safety of co-located workers, helping to ensure the responsible management of radioactive materials during installation and decommissioning activities.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2021

Degradation Diagnostics from the Subsurface of Lithium-Ion Battery Electrodes

Despite the long-established rocking-chair theory of lithium-ion batteries (LIBs), developing novel characterization methodology with higher spatiotemporal resolution facilitates a better understanding of the solid electrolyte interphase studies to shape the reaction mechanisms. In this work, we develop a Xenon ion plasma focused ion beam (Xe + PFIB)-based characterization technique to probe the cross-sectional interface of both ternary cathode and graphite anode electrodes, with the focus on revealing the chemical composition and distribution underneath the electrode surface by in-depth analysis of secondary ions. Particularly, the lithium fluoride is detected in the pristine cathode prior to contact with the electrolyte, reflecting that the electrode degradation is in the form of the loss of lithium inventory during electrode preparation. This degradation is related to the hydrolysis of the cathode material and the decomposition of the PVDF binder. Here, through the quantitative analysis of the transition-metal degradation products, manganese is found to be the dominant element in the newly formed inactive fluoride deposition on the cathode, while no transition metal signal can be found inside the anode electrode. These insights at high resolution implemented via a PFIB-based characterization technique not only enrich the understanding of the degradation mechanism in the LIBs but also identify and enable a high-sensitivity methodology to obtain the chemical survey at the subsurface, which will help remove the capacity-fade observed in most LIBs.

25 ENERGY STORAGE↗
DOE OSTI2023

Tracking Active Phase Behavior on Boron Nitride during the Oxidative Dehydrogenation of Propane Using Operando X-ray Raman Spectroscopy

Hexagonal boron nitride (hBN) is a highly selective catalyst for the oxidative dehydrogenation of propane (ODHP) to propylene. Using a variety of ex situ characterization techniques, the activity of the catalyst has been attributed to the formation of an amorphous boron oxyhydroxide surface layer. The ODHP reaction mechanism proceeds via a combination of surface mediated and gas phase propagated radical reactions with the relative importance of both depending on the surface-to-void-volume ratio. Here we demonstrate the unique capability of operando X-ray Raman spectroscopy (XRS) to investigate the oxyfunctionalization of the catalyst under reaction conditions (1 mm outer diameter reactor, 500 to 550 °C, P = 30 kPa C 3 H 8 , 15 kPa O 2 , 56 kPa He). We probe the effect of a water cofeed on the surface of the activated catalyst and find that water removes boron oxyhydroxide from the surface, resulting in a lower reaction rate when the surface reaction dominates and an enhanced reaction rate when the gas phase contribution dominates. Computational description of the surface transformations at an atomic-level combined with high precision XRS spectra simulations with the OCEAN code rationalize the experimental observations. Finally, this work establishes XRS as a powerful technique for the investigation of light element-containing catalysts under working conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Binding of noble gas atoms by superhalogens

Because of their closed shells, noble gas (Ng) atoms (Ng = Ne, Ar, Kr, and Xe) seldom take part in chemical reactions, yet finding such mechanisms not only is of scientific interest but also has practical significance. Following a recent work by Mayer et al. [Proc. Natl. Acad. Sci. U. S. A. 116, 8167–8172 (2019)] on the room temperature binding of Ar to a superelectrophilic boron site embedded in a negative ion complex, B 12 (CN) 11 - , we have systematically studied the effect of cluster size and terminal ligands on the interaction of Ng by focusing on B 12 X 11 (Ng) (X = H, CN, and BO) and B 12 X 10 (Ng) 2 (X = CN and BO) whose stabilities are governed by the Wade–Mingos rule and on C 5 BX 5 (Ng) (X = H, F, and CN) and C 4 B 2 (CN) 4 (Ng) 2 whose stabilities are governed by the Huckel’s aromaticity rule. Our conclusion, based on density functional theory, is that both the cluster size and the terminal ligands matter—the interaction between the cluster and the Ng atoms becomes stronger with increasing cluster size and the electron affinity of the terminal ligands. Our studies also led to a counter-intuitive finding—removing multiple terminal ligands can enable electrophilic centers to bind multiple Ng atoms simultaneously without compromising their binding strength.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

LDRD conclusion poster - Synthesizing Heterometallic Uranium Single Crystals to Understand the Influence of the Secondary Metals on Uranyl Axial Bond Strength

Understanding how transition metals influenced the chemistry of lanthanide and actinide (f-element) materials is critical for advancing separation technologies, materials design, and coordination chemistry. This project examined how incorporating first-row transition metals affected the structural and spectroscopic properties of f-element coordination polymers. In uranium(VI)-based systems synthesized with 2,6-pyridinedicarboxylic acid (PDC) ligands, single-crystal X-ray diffraction and Raman spectroscopy revealed that the presence of transition metals shortened the uranyl axial bond and induced a blue shift in its symmetric stretching vibration—evidence of increased bond strength. Electronic structure analysis, including Density of States (DOS) calculations using density functional theory (DFT), revealed altered orbital overlaps and highlighted the role of transition metal d-orbitals in modulating bonding. Raman modes were modeled using truncated structural fragments in collaboration with the University of Notre Dame, and although the predicted frequencies were lower than experimental values, they remained within expected ranges. In parallel, similar experiments with cerium (Ce) in the presence of cobalt (Co) and PDC demonstrated multi-step single-crystal-to-single-crystal transformations—behavior not observed in the uranium systems. Initial products included light yellow, orange, and polycrystalline materials. Single-crystal X-ray diffraction studies, conducted in collaboration with the Colorado School of Mines, identified the yellow phase as monometallic Ce(PDC)2(H2O)2·4H2O and the orange phase as heterometallic Ce2Co(PDC)4(H2O)6. After standing in solution for one week, both phases fully transformed into a dark yellow crystalline phase, [Ce3(PDC)5(H2O)8].6(H2O). Remarkably, this transformation was reversible—disturbing the equilibrium by removing some crystals caused reversion to the initial Ce(PDC)2(H2O)2·4H2O phase, highlighting dynamic behavior. All three structures were previously unreported. Solid-state UV-visible and Raman spectroscopy further distinguished these phases, revealing ligand-to-metal charge transfer involving Ce and characteristic d–d transitions from Co(II). The precise mechanism driving these transformations remained unclear; however, pH-dependent experiments confirmed that the transformation did not occur when the pH decreased. Overall, the project demonstrated that transition metals could be employed to tune bonding interactions, structural dimensionality, and optical properties in f-element materials, establishing new pathways for designing functional heterometallic systems. The work resulted in several novel structural discoveries and fostered productive collaborations with the University of Notre Dame and the Colorado School of Mines.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2025

NewLife Nuclear - An Environmentally and Economically Minded Solution for Fusion Energy Waste Handling

Energy demand is rising as a result of innovative and increasingly more energy intensive processes coming to fruition, particularly through the recent interest in the development of AI data centers as well as manufacturing with the push towards increasing domestic manufacturing interest. Fusion energy can provide virtually limitless energy to support this increase in energy demand. Fusion energy concepts, largely classified as magnetic fusion energy (MFE) and inertial fusion energy (IFE) are being pursued, each having unique challenges to overcome before the successful deployment of electricity to the grid. Achieving fusion ignition on the National Ignition Facility, first in December 2022, and eight times since, has demonstrated the scientific viability of the IFE approach. Meanwhile, MFE test stands continue to improve confinement times, making meaningful strides in progressing towards experimental scientific viability. In each of these approaches, an emphasis is placed on generating more power out of the system than what is required to power the system. An under-researched area applicable to both IFE and MFE is handling activated waste coming out of fusion energy systems, both in the course of normal daily operations, as well as in intermittent periods as structural materials may need to be replaced. In the context of an IFE plant system, commonly discussed plant designs suggest targets are ignited within a chamber at a rate of up to one million targets per day. Between each shot, the chamber housing the ignition event will clear a portion of the chamber – resulting in a mixture of vaporized target gas, target debris, and other materials being expelled from the chamber [source]. Additionally, IFE system concepts typically discuss the modularization of plant designs, which are expected to be replaced periodically as the components degrade over time. This would result in the irradiated chamber structure materials, likely metals and alloys, needing to be removed and safely stored. In MFE plant systems, while targets are not ignited at a repetition rate with the frequent chamber clearing as is expected in IFE plant systems, it is anticipated that portions of the confinement area interfacing with the hot plasma will need to be replaced periodically. In each system, without additional investment and research into alternative processing and recycling methods, the result is storing irradiated materials, and other elements in a safe containment area until they are no longer activated. – resulting in significant waste both economic and environmental.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Coal-derived electrically conductive asphalt pavements for snow/ice melting: From laboratory to field

Timely removal of ice and snow from roads is critical to safe, fast, and uninterrupted transportation networks in cold regions. Constructing electrically conductive asphalt pavements to melt the ice and snow on the roads through resistive eating is an emerging alternative technology to traditional snow/ice removal approaches such as utilizing snowplow machines and deicing chemicals. Carbon-based fibers and fillers including carbon fiber and graphite have been widely reported to make electrically conductive hot mix asphalt mixtures for pavement snow/ice-melting applications. This study aimed to develop and demonstrate a novel type of electrically conductive asphalt pavements for snow/ice-melting, which utilizes electrically conductive cold mix asphalt (CMA) mixtures incorporating coal-derived carbon-based coke aggregate as resistive heating elements. Both laboratory experiments and field tests were conducted to investigate the electrical, mechanical, and thermal properties of such electrically conductive asphalt mixtures and pavements. The laboratory experiment results indicated that the electrically conductive CMA mixtures incorporating coke aggregate had sufficiently high electrical conductivity and satisfactory mechanical performance and the pavement prototype slab utilizing a thin layer of such CMA mixtures could successfully raise the pavement surface temperatures to 8.3–11.7 °C rom a low temperature of –5 °C with an input power density of 473 W/m 2 . The field test results showed that the full-scale coal-derived electrically conductive asphalt pavements were able to increase the pavement surface temperatures when electricity was applied, but the magnitude of temperature increase was highly dependent on the power density. Furthermore, it is promising to use coke aggregate to construct coal-derived electrically conductive asphalt pavements for snow/ice melting.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A reliable workflow for improving nanoscale X-ray fluorescence tomographic analysis on nanoparticle-treated HeLa cells

Abstract Scanning X-ray fluorescence (XRF) tomography provides powerful characterization capabilities in evaluating elemental distribution and differentiating their inter- and intra-cellular interactions in a three-dimensional (3D) space. Scanning XRF tomography encounters practical challenges from the sample itself, where the range of rotation angles is limited by geometric constraints, involving sample substrates or nearby features either blocking or converging into the field of view. This study aims to develop a reliable and efficient workflow that can (1) expand the experimental window for nanoscale tomographic analysis of local areas of interest within a laterally extended specimen, and (2) bridge 3D analysis at micrometer and nanoscales on the same specimen. We demonstrate the workflow using a specimen of HeLa cells exposed to iron oxide core and titanium dioxide shell (Fe3O4/TiO2) nanocomposites. The workflow utilizes iterative and multiscale XRF data collection with intermediate sample processing by focused ion beam (FIB) sample preparation between measurements at different length scales. Initial assessment combined with precise sample manipulation via FIB allows direct removal of sample regions that are obstacles to both incident X-ray beam and outgoing XRF signals, which considerably improves the subsequent nanoscale tomography analysis. This multiscale analysis workflow has advanced bio-nanotechnology studies by providing deep insights into the interaction between nanocomposites and single cells at a subcellular level as well as statistical assessments from measuring a population of cells.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2025

Geometric Interpretation of the Cluster Location Problem Part II: Application to the Pahala, Hawaii, Earthquake Sequence

In the companion “Theory” article, we presented a new framing of the seismic location problem in terms of differential geometry (Harris et al., 2025). From that viewpoint, we developed a “project and correct” approach for estimating the relative locations of earthquakes. Here, in this study, we use project and correct to estimate high-precision relative locations of events from an earthquake sequence beneath the town of Pahala, Hawaii, using high-precision correlation-derived picks. The sequence was active from 2020 through 2022 and produced many highly correlated signals at Hawaii Volcano Observatory (HVO) stations on the island of Hawaii. The data we inverted consisted of 2882 events with observations at 5 HVO stations. For comparison with the travel-time image, we also produced conventional hypocenter solutions using both the Bayesloc program (Myers et al., 2007, 2009) and a purpose-built double-difference code. There were obvious structural elements in the resulting image, the resolution of which we used to test the performance of the project and the correct algorithm. For the projection step, we first produced a 3D local basis using an singular value decomposition (SVD) of the 2882 groups of times. Projection of the travel-time vectors into this basis resulted in an image with structures similar to those produced by our conventional locators, but with distortion as predicted by theory. Removing the distortion requires an inverse operator generated from the metric tensor at the geometric centroid of the events. We compared two approaches to obtaining such an inverse operator. The first uses an estimate of the geographic centroid of the event cloud from the centroid of the travel-time data. The second approach uses the centroid of the conventionally produced locations. The first approach produces a corrected image very similar to the conventional results, but with a rotation. The corrected image produced using the conventionally derived centroid is a near-exact match to the conventional locations.

Dodge, Douglas A. [Lawrence Livermore National Lab↗
DOE OSTI2020

Fourier Analysis-Based Evolutionary Multi-Objective Multiphysics Optimization of Liquid-Cooled Heat Sinks

Optimal heat dissipation in power modules can significantly increase their power density. Removing the generated heat is critical for capturing the benefits of advanced semiconductor materials and improving the reliability of the device operation. This article proposes a design optimization method for liquid-cooled heat sinks that use a Fourier analysis-based tool and an evolutionary optimization algorithm to optimize the heat sink geometry for specified objectives. The optimized heat sink geometry is then compared with state-of-the-art solutions in literature based on finite element analysis of different designs. The proposed methodology can develop complex geometries that outperform the conventional heat sink geometries.

Sahu, Raj↗
DOE OSTI2026

Understanding_the_deactivation_mechanisms_of_ethanol_conversion_over_Cu-Y_Beta_catalyst

Direct conversion of bioethanol to C₃⁺olefins is a promising pathway for sustainable aviation fuel (SAF) production, but catalyst deactivation limits long-term operation. The stability and deactivation mechanisms of multifunctional Cu–Y/Beta zeolite catalysts were investigated for ethanol-to-olefins conversion over 300 h time-on-stream in the presence of H2. Catalytic testing reveals progressive losses in ethanol conversion and C₃⁺ olefin selectivity accompanied by increased acetaldehyde formation. The catalyst testing studies correlate with a suite of characterizations of fresh, spent, and regenerated catalysts to identify the deactivation factors. The loss of Y Lewis acid sites is the primary deactivation element. Reversible acid site deactivation is caused by coke deposition, which blocks Y-derived Lewis acid sites responsible for aldol condensation, MPV reduction, and alcohol dehydration. Minor irreversible deactivation is observed and possibly results from hydrothermal dehydroxylation of Y–silanol interactions, resulting in permanent loss of Lewis acidity without zeolite framework degradation or Y aggregation. Cu sites undergo limited agglomeration into small nanoparticles but contribute insignificantly to catalyst deactivation, under the investigated time frame. Oxidative regeneration removes coke and redistributes Cu sites, leading to full recovery of the initial catalytic performance though the Y Lewis acid sites are unable to fully recover. These findings establish Lewis acid site degradation as the primary deactivation mechanism impacting long-term catalyst stability

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Actinium purification from elemental components of stainless steel and aluminum

Radioisotopes of actinium are valuable because of their potential use in the medical industry. Actinium-225 shows promise for treating disease—like cancer—via targeted alpha therapy (TAT), and the longer-lived actinium-227 is the parent of two radionuclides with potential application in TAT radiopharmaceuticals, thorium-227 and radium-223. Continuing progress in the development of these medical applications requires robust and diverse methods for extracting and purifying actinium from a wide range of matrices. Further, to define the strengths and limitations of separation methods commonly employed in actinium processing, we characterized the performance of ion exchange and extraction chromatographic steps for removing contaminants found in stainless steel (chromium, iron, nickel, and silicon) and aluminum from actinium-227. While AG-MP1 anion exchange resin with HCl (aq) successfully removed iron, cation exchanger AG 50W-X8 removed most other contaminants. The most persistent contaminant was aluminum, which was removed using a DGA Normal extraction resin. These results are presented within the context of applying that methodology toward actinium purification strategies.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Validation and Parametric Investigations of an Internal Permanent Magnet Motor Using a Lumped Parameter Thermal Model

One of the key challenges for the electric vehicle industry is to develop high-power-density electric motors. Achieving higher power density requires efficient heat removal from inside the motor. In order to improve thermal management, a multiphysics modeling framework that is able to accurately predict the behavior of the motor, while being computationally efficient, is essential. This paper first presents a detailed validation of a lumped parameter thermal network (LPTN) model of an Internal Permanent Magnet synchronous motor within the commercially available MOTOR-CAD modeling environment. The validation is based on temperature comparison with experimental data and with more detailed finite element analysis (FEA). All critical input parameters of the LPTN are considered in detail for each layer of the stator, especially the contact resistances between the impregnation, liner, laminations, and housing. Finally, a sensitivity analysis for each of the critical input parameters is provided. A maximum difference of 4% - for the highest temperature in the slot-winding and the end-winding - was found between the LPTN and the experimental data. Comparing the results from the LPTN and the FEA model, the maximum difference was 2% for the highest temperature in the slot-winding and end-winding. As for the LPTN sensitivity analysis, the thermal parameter with the highest sensitivity was found to be the liner-to-lamination contact resistance.

42 ENGINEERING↗
DOE OSTI2022

Functional Data Analysis for Extracting the Intrinsic Dimensionality of Spectra: Application to Chemical Homogeneity in the Open Cluster M67

High-resolution spectroscopic surveys of the Milky Way have entered the Big Data regime and have opened avenues for solving outstanding questions in Galactic archeology. However, exploiting their full potential is limited by complex systematics, whose characterization has not received much attention in modern spectroscopic analyses. In this work, we present a novel method to disentangle the component of spectral data space intrinsic to the stars from that due to systematics. Using functional principal component analysis on a sample of 18,933 giant spectra from APOGEE, we find that the intrinsic structure above the level of observational uncertainties requires ≈10 functional principal components (FPCs). Our FPCs can reduce the dimensionality of spectra, remove systematics, and impute masked wavelengths, thereby enabling accurate studies of stellar populations. To demonstrate the applicability of our FPCs, we use them to infer stellar parameters and abundances of 28 giants in the open cluster M67. We employ Sequential Neural Likelihood, a simulation-based Bayesian inference method that learns likelihood functions using neural density estimators, to incorporate non-Gaussian effects in spectral likelihoods. By hierarchically combining the inferred abundances, we limit the spread of the following elements in M67: Fe ≲ 0.02 dex; C ≲ 0.03 dex; O, Mg, Si, Ni ≲ 0.04 dex; Ca ≲ 0.05 dex; N, Al ≲ 0.07 dex (at 68% confidence). Our constraints suggest a lack of self-pollution by core-collapse supernovae in M67, which has promising implications for the future of chemical tagging to understand the star formation history and dynamical evolution of the Milky Way.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Eu speciation in apatite at 1 bar: An experimental study of valence-state partitioning by XANES, lattice strain, and Eu/Eu* in basaltic systems

Abstract Partition coefficients for rare earth elements (REEs) between apatite and basaltic melt were determined as a function of oxygen fugacity (fO2; iron-wüstite to hematite-magnetite buffers) at 1 bar and between 1110 and 1175 °C. Apatite-melt partitioning data for REE3+ (La, Sm, Gd, Lu) show near constant values at all experimental conditions, while bulk Eu becomes more incompatible (with an increasing negative anomaly) with decreasing fO2. Experiments define three apatite calibrations that can theoretically be used as redox sensors. The first, a XANES calibration that directly measures Eu valence in apatite, requires saturation at similar temperature-composition conditions to experiments and is defined by: ( E u 3 + ∑ E u ) Apatite = 1 1 + 10 - 0.10 ± 0.01 × l o g ⁡ ( f o 2 ) - 1.63 ± 0.16 . The second technique involves analysis of Sm, Eu, and Gd in both apatite and coexisting basaltic melt (glass), and is defined by: ( Eu E u * ) D Sm × Gd = 1 1 + 10 - 0.15 ± 0.03 × l o g ⁡ ( f o 2 ) - 2.46 ± 0.41 . The third technique is based on the lattice strain model and also requires analysis of REE in both apatite and basalt. This calibration is defined by ( Eu E u * ) D lattice strain = 1 1 + 10 - 0.20 ± 0.03 × l o g ⁡ ( f o 2 ) - 3.03 ± 0.42 . The Eu valence-state partitioning techniques based on (Sm×Gd) and lattice strain are virtually indistinguishable, such that either methodology is valid. Application of any of these calibrations is best carried out in systems where both apatite and coexisting glass are present and in direct contact with one another. In holocrystalline rocks, whole rock analyses can be used as a guide to melt composition, but considerations and corrections must be made to either the lattice strain or Sm×Gd techniques to ensure that the effect of plagioclase crystallization either prior to or during apatite growth can be removed. Similarly, if the melt source has an inherited either a positive or negative Eu anomaly, appropriate corrections must also be made to lattice strain or Sm×Gd techniques that are based on whole rock analyses. This being the case, if apatite is primary and saturates from the parent melt early during the crystallization sequence, these corrections may be minimal. The partition coefficients for the REE between apatite and melt range from a maximum DEu3+ = 1.67 ± 0.25 (as determined by lattice strain) to DLu3+ = 0.69 ± 0.10. The REE partition coefficient pattern, as observed in the Onuma diagram, is in a fortuitous situation where the most compatible REE (Eu3+) is also the polyvalent element used to monitor fO2. These experiments provide a quantitative means of assessing Eu anomalies in apatite and how they be used to constrain the oxygen fugacity of silicate melts.

Geochemistry & Geophysics↗
DOE OSTI2020

Fourier Analyses of High-Order Continuous and Discontinuous Galerkin Methods

In this paper, we present a Fourier analysis of wave propagation problems subject to a class of continuous and discontinuous discretizations using high-degree Lagrange polynomials. This allows us to obtain explicit analytical formulas for the dispersion relation and group velocity and, for the first time to our knowledge, characterize analytically the emergence of gaps in the dispersion relation at specific wavenumbers, when they exist, and compute their specific locations. Wave packets with energy at these wavenumbers will fail to propagate correctly, leading to significant numerical dispersion. We also show that the Fourier analysis generates mathematical artifacts, and we explain how to remove them through a branch selection procedure conducted by analysis of eigenvectors and associated reconstructed solutions. The higher frequency eigenmodes, named erratic in this study, are also investigated analytically and numerically.

97 MATHEMATICS AND COMPUTING↗