Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “magnitude”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The impact of mineral reactive surface area variation on simulated mineral reactions and reaction rates

Reactive transport modeling is an essential tool to simulate complex geochemical reactions in porous media that can impact formation properties including porosity and permeability. However, simulating these reactions is challenging due to uncertainties in model parameters, particularly mineral surface areas. Imaging has emerged as a powerful means of estimating model parameters including porosity and mineral abundance, accessibility and accessible surface area. However, these parameters, particularly mineral accessible surface area, vary with image resolution. This work aims to enhance understanding of the impact of image resolution and other means of estimating mineral reactive surface area on simulated mineral reactions and reaction rates. Mineral surface areas calculated from images with resolutions of 0.34 μm and 5.71 μm were used to simulate mineral reactions in the context of geologic CO 2 sequestration in the Paluxy formation at the continuum scale. Additional simulations were carried out using BET surface areas collected from the literature and geometric surface areas. Simulations were run for 7300 days and mineral volume fractions and effluent ion concentrations tracked and compared. Variations in mineral surface areas measured from images are within 1 order of magnitude and yield similar simulation results, indicating the impact of image resolution on simulated reactions and reaction rates is minimum for the resolutions and sample considered. In comparison, surface areas obtained from BET and geometric approaches are 1–4 orders of magnitude higher than image-obtained surface areas and result in greater simulated reaction rates and extents. Minerals with high reaction rates (calcite and siderite) are most impacted by surface area values at short times where simulated mineral volume fractions at longer times agree relatively well, even for simulations with several orders of magnitude variation in surface area. Phases with lower reaction rates, such as K-feldspar and muscovite, are predominantly impacted over longer times where variations in surface areas impact reaction extents and porosity evolution. Overall, variations in surface areas due to image resolution are small and result in little variation in simulated reactions and reaction rates while there are significant variations in simulation results when other surface area estimates are used. These variations, however, depend on the reactivity of the mineral phase where surface areas of fast-reacting phases largely impact simulations at short (hours to days) timescales and surface areas of slower-reactive phases impact longer term simulations (years).

58 GEOSCIENCES↗
DOE OSTI2024

Health and air pollutant emission impacts of net zero CO2 by 2050 scenarios from the energy modeling forum 37 study

Carbon dioxide and non-greenhouse gas air pollutants are emitted from many of the same sources. Decarbonization actions thus typically yield air pollutant emission reductions, resulting in significant air quality benefits. Although several studies have highlighted this connection, including in the context of net zero carbon emission targets, substantial uncertainty remains regarding how alternative technological pathways to this goal will affect the spatial distribution and magnitude of air pollutants. Comprehensive multi-model and multi-scenario analyzes are needed to explore the relative impacts of alternative pathways. Here, our study begins to address this gap by leveraging the results from the recent Energy Modeling Forum 37 inter-model comparison exercise on U.S. decarbonization pathways. Comparing the results of the six teams who submitted air pollutant emissions suggests that strategies that target net zero U.S. carbon emissions would yield significant reductions in many air pollutants, and that this finding is generally robust across pathways. However, some energy sources, such as biomass and fossil fuels with carbon capture, will emit air pollutants and can potentially influence the magnitude, spatial distribution, and even sign of localized air pollutant emission changes. In the second part of this analysis, a simplified air quality and health impacts screening model is used to evaluate the air quality impacts in 2035 of sectoral emission changes from the three models that provided sectoral detail. Relative to a reference scenario, a net zero pathway is estimated to reduce fine particulate matter concentrations across the contiguous U.S., with health benefits from reduced mortality ranging from $\$$65 billion to $\$$250 billion in 2035 alone (2023$\$$s). These benefits would be expected to grow over time as the net zero trajectory becomes more stringent. Both the magnitude of potential benefits and the substantial variation of the projections across models underscore the need for an EMF-like inter-model comparison exercise focused on air quality.

Air pollutants↗
DOE OSTI2021

Variable frequency response testbed to validate protective relays up to 20 kHz☆

Microgrids operate in islanding or utility grid mode. In islanding mode, distributed energy sources with power electronics inverters can generate fault currents over a wide range of frequencies. To gain a better understanding of the nature of these events, it is crucial to measure currents at frequencies of up to 20 kHz at transient states with protective relays. Protective relays at these frequency responses should be tested prior to being installed in microgrids. However, testing at high-frequency currents cannot be done with Relay Test Equipment (RTE) that generate signals up to 50th Harmonic for 60 Hz (3 kHz). In this study, a frequency response method and testbed for high-frequency currents of up to 20.16 kHz was developed and tested. The testbed utilized a RTE and High-Frequency Power Source/ Amplifier. Furthermore, frequency relay experiments were run at frequencies of up to 1 and 20.16 kHz at a current magnitude of 0.1 and 0.5 A, respectively. The events were collected, and the current magnitude and frequency were plotted at time and frequency domain, respectively. The study of current magnitude and frequencies of up to 20.16 kHz for two overcurrent relays was performed and validated the variable frequency testbed and methodology.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

A data-driven sensor placement approach for detecting voltage violations in distribution systems

Stochastic fluctuations in power injections from distributed energy resources (DERs) combined with load variability can cause constraint violations (e.g., exceeded voltage limits) in electric distribution systems. To monitor grid operations, sensors are placed to measure important quantities such as the voltage magnitudes. Here, in this paper, we consider a sensor placement problem which seeks to identify locations for installing sensors that can capture all possible violations of voltage magnitude limits. We formulate a bilevel optimization problem that minimizes the number of sensors and avoids false sensor alarms in the upper level while ensuring detection of any voltage violations in the lower level. This problem is challenging due to the nonlinearity of the power flow equations and the presence of binary variables. Accordingly, we employ recently developed conservative linear approximations of the power flow equations that overestimate or underestimate the voltage magnitudes. By replacing the nonlinear power flow equations with conservative linear approximations, we can ensure that the resulting sensor locations and thresholds are sufficient to identify any constraint violations. Additionally, we apply various problem reformulations to significantly improve computational tractability while simultaneously ensuring an appropriate placement of sensors. Lastly, we improve the quality of the results via an approximate gradient descent method that adjusts the sensor thresholds. We demonstrate the effectiveness of our proposed method for several test cases, including a system with multiple switching configurations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

The solubility of ErPO 4 and Er speciation in hydrothermal fluids at varying pH and salinity between 350 and 450 °C

The rare earth elements (REE) are important for the green-energy transition and can be incorporated into the REE phosphates, such as xenotime-(Y), which also hosts heavy REE (Tb– Lu). Xenotime-(Y) is a common accessory mineral in metamorphic rocks and a range of mineral deposits where it controls the mobility of heavy REE, however, the impact of high temperature aqueous fluids on the behavior of heavy REE is largely unknown. Thermodynamic modeling can be utilized as a tool to predict the mobility of REE in hydrothermal aqueous fluids, but must be supported by accurate experimental data. Here, we measured the solubility of endmember synthetic xenotime-structured ErPO4 in NaCl-HCl-NaOH-bearing aqueous solutions at 350 °C and water vapor saturation pressure, at 400 and at 450 °C and 500 bar using batch-type Inconel reactors. Erbium speciation was investigated as a function of pH from 2.8 to 8, where Er chloride species are predominant at acidic conditions (pH <3) and Er hydroxyl complexes are predominant at near- neutral to alkaline conditions (pH >3). At pH 7–9, the measured ErPO 4 solubility (-9.8 to -7.5 log m Er ) is up to 2.5 orders of magnitude lower than thermodynamic predictions (-9.4 to -6.7 log m Er ) using existing thermodynamic databases. At pH 2–3, the predicted ErPO 4 solubility is ~0.5 orders of magnitude higher at 350 °C and ~1 order of magnitude lower at 450 °C compared to experimentally measured Er concentrations. The thermodynamic properties of aqueous Er species were therefore revised in this study. The partial molal Gibbs energy of formation (Δ f G 0 T,P ) for aqueous Er hydroxyl and chloride species are optimized using GEMSFITS and the logarithmic formation constants (logβ n ( Cl,OH) ) were derived at each experimental temperature and pressure. The updated thermodynamic properties for Er hydroxyl species (Er(OH) +2 , Er(OH) 2 + , and Er(OH) 3 0 ) show that their stability shifts to more acidic conditions at and below 400 °C. The Er chloride species (ErCl +2 and ErCl 2 + ) show increased stability compared to Er hydroxyl species at temperatures of 450 °C and 0.01 mol/kg NaCl. The updated thermodynamic properties are implemented into the GEM-Selektor modeling package to investigate the mobility of Er in saline hydrothermal fluids in equilibrium with alkaline rocks. Importantly, the updated properties for Er hydroxyl species result in low Er solubility at rock equilibrated pH conditions due to an expanded hydroxyl predominance zone, but lower aqueous complex stability overall, whereas previous models suggest greater stability for aqueous Er species. Furthermore, ErPO 4 solubility increases with decreasing temperature due to the deprotonation of HCl, which increases the acidity of hydrothermal fluids and the availability of Cl - to complex with the REE. These simulations highlight how fluid-rock reaction and temperature affect the mobility of REE in hydrothermal ore-forming systems.

58 GEOSCIENCES↗
DOE OSTI2021

Unraveling the implications of finite specimen size on the interpretation of dynamic experiments for polycrystalline aluminum through direct numerical simulations

Normal and Pressure-shear plate impact (NPI and PSPI) experiments are popular experimental techniques for studying the mean-field macroscopic behavior of polycrystalline metals under high-rate dynamic loading. However, since both configurations rely upon geometry for subjecting the specimen to high strain rates, these experiments often involve a limited specimen size. Moreover, because of the inherent heterogeneities present within polycrystalline metals, it is difficult to ascertain if the size of the specimen and/or regions where measurements are made are sufficiently large for making representative inferences about the mean-field macroscopic properties from single-point velocity measurements. In the present study, we quantify the expected measurement variability on observable point measurements in NPI and PSPI experiments by carrying out direct numerical simulations (DNS) of statistically representative polycrystalline microstructures subjected to dynamic compression and compression-shear loading. In particular, we consider the role of specific material heterogeneities (e.g. the grain-to-grain difference in size, crystallographic orientation) on dispersion in the normal and transverse particle velocity records and on local fluctuations in key state variables (e.g. velocity, accumulated plastic strain) by incorporating these effects directly into a representative synthetic microstructure geometry and crystalline description of pure polycrystalline aluminum. The form of the present study is a large parametric investigation, consisting of ten ensembles of one hundred simulations. Each of the thousand simulations reflects a randomly realized synthetic microstructure in one of five cases of decreasing average grain size for the two loading configurations. Our analysis of the DNS results demonstrates that for both of these experimental configurations, the grain size directly correlates with the coefficient of variation (CV) in simulated point measurements, showing a convergent decrease in CV to zero (i.e. particle velocity record approaches the mean-field value) with decreasing grain size. Remarkably, the magnitude of variations in the particle velocity record is shown to be largest where the deviatoric stresses are most significant. In the case of NPI, this occurs at the elastic and plastic wavefront, whereas, in the case of PSPI, the magnitude of fluctuations are approximately constant throughout the experimental window time. The reasoning for the scatter in particle velocity due to the heterogeneous microstructure is demonstrated to be dependent on the mechanisms for accommodating deformation and on the interaction of reflection waves generated at sites of heterogeneities occurring at the scale of grains. Lastly, we develop a power-law description for the magnitude of scattering versus characteristic length, which provides a statistical framework for assessing the required number of grains per characteristic specimen dimension for minimizing scatter within these two experimental configurations (NPI, PSPI).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Bioinspired nonwetting surfaces for corrosion inhibition over a range of temperature and corrosivity

Applications of superhydrophobic (SHS) and lubricant infused surfaces (LIS) involve exposure to corrosive environments from the acidic to the basic, at a range of temperatures, that are not fully characterized. Here, we present for the first time a multifactorial study of the effects of surface fabrication method, surface modification, surface functionalization time, temperature and pH of the immersion medium on the corrosion performance of nonwetting copper surfaces. Bioinspired SHS and LIS fabricated using facile methods of etching and electrodeposition are systematically assessed using potentiodynamic polarization measurements for their corrosion resistance in saline solution (pH ≈7) over a temperature range 23–85°C. SHS and LIS are shown to exhibit diminished corrosion rate, by up to two orders of magnitude, compared to bare copper surface. An Arrhenius model is developed for the first time, describing the temperature- dependent corrosion rate of SHS and LIS. Electrochemical impedance spectroscopy is used to show that corrosion resistance of LIS is larger by three orders of magnitude in extremely acidic (pH = 1) and by an order magnitude in extremely alkaline (pH = 14) media compared to bare copper surface. Etched LIS are generally more resistant to corrosion compared to SHS at all temperatures with excellent microstructural durability.

42 ENGINEERING↗
DOE OSTI2025

Implementation of a Mesh refinement algorithm into the quasi-static PIC code QuickPIC

Plasma-based acceleration (PBA) has emerged as a promising candidate for the accelerator technology used to build a future linear collider and/or an advanced light source. In PBA, a trailing or witness particle beam is accelerated in the plasma wave wakefield (WF) created by a laser or particle beam driver. The WF is often nonlinear and involves the crossing of plasma particle trajectories in real space and thus particle-in-cell methods are used. The distance over which the drive beam evolves is several orders of magnitude larger than the wake wavelength. This large disparity in length scales is amenable to the quasi-static approach. Three-dimensional (3D), quasi-static (QS), particle-in-cell (PIC) codes, e.g., QuickPIC, have been shown to provide high fidelity simulation capability with 2-4 orders of magnitude speedup over 3D fully explicit PIC codes. In PBA, the witness beam needs to be matched to the focusing forces of the WF to reduce the emittance growth. In some linear collider designs, the matched spot size of the witness beam can be 2 to 3 orders of magnitude smaller than the spot size (and wavelength) of the wakefield. Such an additional disparity in length scales is ideal for mesh refinement where the WF within the witness beam is described on a finer mesh than the rest of the WF. A mesh refinement scheme is described that has been implemented into the 3D QS PIC code, QuickPIC. Very fine (high) resolution is used in a small spatial region that includes the witness beam and progressively coarser resolutions in the rest of the simulation domain. A fast multigrid Poisson solver has been implemented for the field solve on the refined meshes and a Fast Fourier Transform (FFT) based Poisson solver is used for the coarse mesh. The code has been parallelized with both MPI and OpenMP, and the parallel scalability has also been improved by using pipelining. A preliminary adaptive mesh refinement technique is described to optimize the computational time for simulations with an evolving witness beam size. Several test problems are used to verify that the mesh refinement algorithm provides accurate results. Additionally, the results are benchmarked against highly resolved simulations exhibiting near-azimuthal symmetry, performed using QPAD—a novel hybrid QS PIC code that uses a PIC description in the coordinates (r, ct – z) and a gridless description in the azimuthal angle, Φ.

Linear collider↗
DOE OSTI2026

A hybrid Monte Carlo-deterministic second moment method with efficient variance reduction

In this work, we present a hybrid method that combines Monte Carlo with deterministic finite element methods to solve a linear Boltzmann transport equation. Our hybrid method runs orders of magnitude faster than Monte Carlo, without sacrificing accuracy, for a proxy problem from radiative transfer that contains both optically-thick and optically-thin material. We believe that this is the first demonstration of a hybrid Second Moment Method in more than one spatial dimension, the first to consider more than one material, and the first to use variance reduction. Our variance reduction approach arises from an asymptotic analysis in which we show that the magnitude of the scattering source grows without bound. We transform the problem to compute the deviation of the radiation intensity from isotropy. The magnitude of the source in the transformed problem is bounded, and the quality of the hybrid method solution is dramatically improved by a substantial reduction in the variance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Allometric-kinetic model predictions of concentration ratios for anthropogenic radionuclides across animal classes and food selection

Protection of the environment from radiation fundamentally relies on dose assessments for non-human biota. These dose assessments use measured or predicted concentrations of radionuclides in soil or water combined with Concentration Ratios (CRs) to estimate tissue concentrations in animals and plants, yet there is a paucity of CR data relative to the vast number of potential taxa and radioactive contaminants in the environment and their taxon-specific ecosystems. Because there are many taxa each having very different behaviors and biology and many possible bioavailable radionuclides, CRs have the potential to vary by orders-of-magnitude, as often seen in published data. Given the diversity of taxa, the International Commission on Radiological Protection (ICRP) has categorized non-human biota into 12 representative organisms (ROs), while the U.S. Department of Energy (DOE) uses the non-taxon specific categories of terrestrial, riparian, and aquatic animals. The question we examine here is: do these categorizations adequately “represent” the broad diversity of animals? To explore this question, here we analyze an Allometric-Kinetic (A-K) model to calculate radionuclide-specific CRs for common animal classes which are then further subcategorized into herbivores, omnivores, carnivores, and invertebrate detritivores. Comparisons in CRs among animal classes exhibited only small differences, but there was order of magnitude differences between herbivores, carnivores, and especially detritivores, for many radionuclides of interest. These findings suggest that the ICRP ROs and the DOE categories are reasonable, but their accuracy could be improved by including sub-categories related to animal dietary ecology. Finally, comparisons of A-K model predicted CR values to published CRs show order-of-magnitude variations, providing justification for additional studies of assimilation across radionuclides and animal classes.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2022

Soil frost controls streamflow generation processes in headwater catchments

The relationship between snowmelt and spring streamflow is changing under warming temperatures and diminishing snowpack. At the same time, the hydrologic connectivity across catchment landscape elements, such as snowpack and surface wetlands, can play a critical role in controlling the routing of snowmelt to streams. The role of hydrologic connectivity is important in headwater regions of the continental northern latitudes, where catchments have low topographic relief and seasonally frozen ground. Nevertheless, the effects of soil frost on the sequence, timing, and magnitudes of hydrologic events that drive the movement of water from a snowpack to a stream are not fully understood. Therefore, we examine two questions: First, what is the flowpath that snow melt and precipitation from spring rain events takes to generate spring streamflow, and second, what hydrologic, climatic, or landscape variables exert the most control on the magnitude of streamflow? In this work, we use long-term hydrological records from the two reference basins at the Marcell Experimental Forest in northern Minnesota to analyze the cascading effects across precipitation, snow, water table elevation, soil frost, and streamflow in peatland-dominated headwater catchments. We identify a sequence of fill-and-spill effects across the landscape that control the timing of spring streamflow generation. Then, we use stepwise regression to show that soil frost is a key supporting predictor for both the magnitude of streamflow in the spring as it adds significantly to the predictive power of precipitation and water table elevation. Our results highlight the importance of recognizing the role of soil frost, when present, on the partitioning of snowmelt between overland runoff and water table recharge during the critical snowmelt period, as well as the later partitioning between evapotranspiration and subsurface flows.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Shape polarization in the tin isotopes near N = 60 from precision g-factor measurements on short-lived 11/2 - isomers

The g factors of 11/2 - isomers in semimagic 109 Sn and 111 Sn (isomeric lifetimes $\tau$ = 2.9(3) ns and $\tau$ = 14.4(7) ns, respectively) were measured by an extension of the Time Differential Perturbed Angular Distribution technique, which uses LaBr3 detectors and the hyperfine fields of a gadolinium host to achieve precise measurements in a new regime of short-lived isomers. The results, g(11/2 - ; 109 Sn) = -0.186(8) and g(11/2 - ; 111 Sn) = -0.214(4), are significantly lower in magnitude than those of the 11/2 - isomers in the heavier isotopes and depart from the value expected for a near pure neutron h 11/2 configuration. Broken-symmetry density functional theory calculations applied to the sequence of 11/2 - states reproduce the magnitude and location of this deviation. The g(11/2 - ) values are affected by shape core polarization; the odd Oh 11/2 neutron couples to J π = 2 + , 4 + , 6 + ... configurations in the weakly-deformed effective core, causing a decrease in the g-factor magnitudes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Rewetting of soil: Revisiting the origin of soil CO 2 emissions

Rewetting dry soils is associated with a burst of microbial activity and mineralization, which manifests itself as a pulse in soil CO 2 emissions, long-known as the Birch effect. In arid and semi-arid systems, soil CO 2 losses upon rewetting at the end of extended dry periods can contribute a significant fraction to the overall carbon (C) budget. Microbial biomass is one of the sources of mineralized C, as was demonstrated over 30 years ago (Kieft et al., 1987). The present paper offers a perspective on how the field has progressed since the 1987 paper was published in Soil Biology & Biochemistry, what it means in terms of current concerns about global climate change, and the needs and potential emphases of future research. Many studies since 1987 have addressed the origin of this CO 2 pulse, finding multiple possible C sources involving both biotic and abiotic processes. We propose that the magnitude of the rewetting event (Δψ) determines the relative contribution from the array of substrates that contribute to the soil CO 2 pulse upon rewetting. The magnitude of the CO 2 pulse is likely related to soil physical characteristics and to the size of the available C pool, which is partly controlled by plants. Further, the relative contributions of the mechanisms generating soil CO 2 pulses upon rewetting are likely to be modified by climate change. To understand and predict the magnitude of soil CO 2 pulses upon rewetting, we advocate continued cross-disciplinary research involving soil microbial ecology, soil physics, soil chemistry, as well as increased integration and recognition of the importance of plant-soil interactions in controlling the soil C pools contributing to soil CO 2 pulses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A comparative analysis of different van der Waals treatments for molecular adsorption on the basal plane of 2H-MoS 2

The binding energy of hydrogen sulfide, ammonia, ethane, ethylene, butadiene, benzene, toluene, pyridine, pyrrole, and thiophene on the basal plane of the semi-conducting 2H-molybdenum sulfide (MoS 2 ) was calculated with the following flavors of Density Functional Theory (DFT): GGA-PW91, PBE-D3, vdW-DF, optPBE, optB86b, optB88, vdW-TS, and BEEF-vdW. The GGA-PW91 binding energies are negligible (<0.07 eV in magnitude) in all cases. The predictions with vdW-DF and PBE-D3 are the closest (error <0.05 eV) to the isosteric heats of adsorption calculated from reported temperature programmed desorption data for thiophene and butadiene. For all dispersion flavors examined here, the magnitude of the dispersion contribution to the binding energy increases linearly with the number of heavy atoms in the adsorbate, with each atom contributing 0.05 eV (BEEF-vdW) – 0.09 eV (optB88-vdW). Further, this implies that the calculated adsorption constants of molecules larger than acridine (i.e., comprising > 14 non-heavy atoms) can vary by more than four orders of magnitude at industrial conditions depending on the chosen method of dispersion correction. Further, dispersion effects fall off rapidly (>0.03 eV/ non-hydrogen atom/Å) as the adsorbate-surface distance increases.

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

Direct Measurement of the Selective Microwave-Induced Heating of Agglomerates of Dipolar Molecules: The Origin of and Parameters Controlling a Microwave Specific Superheating Effect

Agglomerates of polar molecules in nonpolar solvents are selectively heated by microwave radiation. The magnitude of the selective heating was directly measured by using the temperature dependence of the intensities of the Stokes and anti-Stokes bands in the Raman spectra of p-nitroanisole (pNA) and mesitylene. Under dynamic heating conditions, a large apparent temperature difference (ΔT) of over 100 °C was observed between the polar pNA solute and the nonpolar mesitylene solvent. This represents the first direct measurement of the selective microwave heating process. The magnitude of the selective microwave heating was affected by the properties of the agglomerated pNA. As the concentration of the pNA increases, the magnitude of the selective heating of the pNA was observed to decrease. This is explained by the tendency of the pNA dipoles to orient in an antiparallel fashion in the aggregates as measured by the Kirkwood g value, which decreased with increasing concentration. This effect reduces the net dipole moment of the agglomerates, which decreases the microwave absorption. After the radiation was terminated, the effective temperature of the dipolar molecules returned slowly to that of the medium. The slow heat transfer was modeled successfully by treating the solutions as a biphasic solvent/solute system. Based on modeling and the fact that the agglomerate can be heated above the boiling temperature of the solvent, an insulating layer of solvent vapor is suggested to form around the heated agglomerate, slowing convective heat transfer out of the agglomerate. Furthermore, this is an effect unique to microwave heating.

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

Characterizing Density and Spatial Distribution of Trap States in Ta 3 N 5 Thin Films for Rational Defect Passivation

Tantalum nitride (Ta 3 N 5 ) has gained significant attention as a potential photoanode material, yet it has been challenged by material quality issues. Defect-induced trap states are detrimental to the performance of any semiconductor material. Beyond influencing the performance of Ta 3 N 5 films, defects can also accelerate the degradation in water during desired electrochemical applications. Defect passivation has provided an enormous boost to the development of many semiconductor materials but is currently in its infancy for Ta 3 N 5 . This is in part due to a lack of experimental understanding regarding the spatial and energetic distribution of trap states throughout Ta 3 N 5 thin films. Here, we employ drive-level capacitance profiling (DLCP) to experimentally resolve the spatial and energetic distribution of trap states throughout Ta 3 N 5 thin films. The density of deeper energetic traps is found to reach ∼2.5 to 6 × 10 22 cm –3 at the interfaces of neat Ta 3 N 5 thin films, over an order of magnitude greater than the bulk. In addition to the spatial profile of deep trap states, we report neat Ta 3 N 5 thin films to be highly n-type in nature, owning a free carrier density of ∼9.74 × 10 17 cm –3 . This information, coupled with the present understanding of native oxide layers on Ta 3 N 5 , has facilitated the rational design of a targeted passivation strategy that simultaneously provides a means for catalyst immobilization. Loading catalyst via silatrane moieties suppresses the density of defects at the surface of Ta 3 N 5 thin films by two orders of magnitude, while also reducing the free carrier density of films by over one order of magnitude, effectively dedoping the films to ∼2.40 × 10 16 cm –3 . The surface passivation of Ta 3 N 5 films translates to suppressed defect-induced trapping and recombination of photoexcited carriers, as determined through absorption, photoluminescence, and transient photovoltage. Here, this illustrates how developing a deeper understanding of the distribution and influence of defects in Ta 3 N 5 thin films has the potential to guide future works and ultimately accelerate the integration and development of high-performance Ta 3 N 5 thin film devices.

defect passivation↗
DOE OSTI2020

Detection Limits and Near‐Field Ground Motions of Fast and Slow Earthquakes

Abstract We investigate theoretical limits to detection of fast and slow seismic events, and spatial variations of ground motion expected from M 6 earthquakes at short epicentral distances. The analyses are based on synthetic velocity seismograms calculated with the discrete wavenumber method assuming seismic velocities and attenuation properties of the crust in Southern California. The examined source properties include different magnitudes ( M −1.0 to M 6.0), static stress drops (0.1–10 MPa), and slow and fast ruptures (0.1–0.9 of shear wave velocity). For the M 6 events we also consider variations in rise times producing crack‐ and pulse‐type events and different rupture directivities. Slow events produce ground motion with considerably lower amplitude than corresponding regular fast earthquakes with the same magnitude, and hence are significantly more difficult to detect. The static stress drop and slip rise time also affect the maximum radiated seismic motion, and hence event detectability. Apart from geometrical factors, the saturation and depletion of seismic ground motion at short epicentral distances stem from radiation pattern, earthquake size (magnitude, stress drop), and rupture directivity. The rupture velocity, rise time, and directivity affect significantly the spatial pattern of the ground motions. The results can help optimizing detection of slow and fast small earthquakes and understand the spatial distribution of ground motion generated by large events.

58 GEOSCIENCES↗
DOE OSTI2020

Semiquantitative Estimates of Rainfall Variability During the 8.2 kyr Event in California Using Speleothem Calcium Isotope Ratios

A multiproxy record from a fast-growing stalagmite reveals variable hydroclimate on the California coast across the 8.2 kyr event and a precursor event likely caused by initial drainage of proglacial Lake Agassiz. Using speleothem δ 44 Ca, we develop the first semiquantitative estimates of paleorainfall variability for California through calibration with measurements of the modern climate and cave environment. We find that the magnitude of rainfall variability during the 8.2 kyr event approached the multiyear variability observable in the recent past (1950–2019) and the magnitude of variability during the precursor event likely exceeded this range. Additionally, we observe other instances of multidecadal variability comparable in magnitude to the precursor event during the record. Finally, our work suggests that speleothem calcium isotope ratios are a powerful semiquantitative means to reconstruct paleorainfall, although numerous factors must be assessed in each cave system before applying this approach.

58 GEOSCIENCES↗