Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Control Samples”

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 217 records · Page 12

Action Effectiveness Monitoring and Research of Dredged Material Placement at Woodland Islands (Final Report)

Benthic invertebrates play critical ecosystem roles including the breakdown of organic matter, sediment mixing, and nutrient cycling. In the lower Columbia River and estuary, benthic invertebrates provide a foraging resource for threatened and endangered juvenile salmon. Driven by the goal to create low velocity, shallow water, and riparian shrub habitats to benefit juvenile salmon, the USACE placed 237,000 CY of dredged material, resulting in the deposition of 13.5 acres of sand on the off-channel margins at Woodland Islands. The Pacific Northwest National Laboratory (PNNL) designed and implemented an action effectiveness research study to understand how dredged material placement at Woodland Islands affected sediment conditions and benthic invertebrate assemblages. A Before-After-Control-Impact, or BACI, study design was used to evaluate the response of environmental conditions and benthic invertebrates to dredged material placement. The design included one impact site and two control sites. Sampling occurred for two years prior to dredged material placement and two years after. Our study found that spatial variation was a significant factor for both environmental and biological response variables which suggests local conditions are important considerations for mechanisms affecting benthic assemblages. We found that off-channel habitats across all locations sampled were producing benthic invertebrates, many of which are common prey items for juvenile salmon and steelhead—e.g., insects, chironomids, crustaceans, and corophium. There was a significant BACI effect (i.e., an effect at the dredged material placement site, relative to conditions at the control sites) for concentration of carbon and ammonium in sediment, but not for phosphorous concentrations or for percent sand. The estimated abundance for the three invertebrate response variables—total abundance, total chironomid abundance, and total corophium abundance—was significantly lower at the impact site after dredged material placement, compared to the control sites. At Woodland Island, the estimated mean abundance for all invertebrates combined decreased 28% after dredged material placement. Estimated mean abundances of chironomid and corophium decreased by 8% and 88%, respectively. While invertebrate abundances were lower after dredged placement, the composition of benthic invertebrates was similar before and after placement suggesting that as the new habitat feature evolves recolonization will likely follow. These findings provide a foundation for understanding potential benefits and consequences of repurposing dredged material for habitat creation in the LCRE.

54 ENVIRONMENTAL SCIENCES↗

Operator Relaxation and the Optimal Depth of Classical Shadows

Classical shadows are a powerful method for learning many properties of quantum states in a sample-efficient manner, by making use of randomized measurements. Here we study the sample complexity of learning the expectation value of Pauli operators via “shallow shadows,” a recently proposed version of classical shadows in which the randomization step is effected by a local unitary circuit of variable depth t. Here we show that the shadow norm (the quantity controlling the sample complexity) is expressed in terms of properties of the Heisenberg time evolution of operators under the randomizing (“twirling”) circuit—namely the evolution of the weight distribution characterizing the number of sites on which an operator acts nontrivially. For spatially contiguous Pauli operators of weight k, this entails a competition between two processes: operator spreading (whereby the support of an operator grows over time, increasing its weight) and operator relaxation (whereby the bulk of the operator develops an equilibrium density of identity operators, decreasing its weight). From this simple picture we derive (i) an upper bound on the shadow norm which, for depth t~log⁡(k), guarantees an exponential gain in sample complexity over the t=0 protocol in any spatial dimension, and (ii) quantitative results in one dimension within a mean-field approximation, including a universal subleading correction to the optimal depth, found to be in excellent agreement with infinite matrix product state numerical simulations. Our Letter connects fundamental ideas in quantum many-body dynamics to applications in quantum information science, and paves the way to highly optimized protocols for learning different properties of quantum states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Vegetation Warming Experiment: Leaf Mass Area, Leaf Carbon and Nitrogen Content, Utqiagvik (Barrow), Alaska, 2018

Leaf mass per area (LMA), leaf carbon and nitrogen content of vegetation (Arctagrostis latifolia, Petasites frigidus) within warming chambers and paired control plots. Samples include gas exchange measurement leaves and an end of season bulk harvest. The files included in this data package are in .csv format, and include 2 data files and 4 metadata files. See related datasets for plant physiology, phenology and environmental conditions of these samples. These data were collected in 2018 as part of a series of single-season warming experiments on tundra vegetation on the Barrow Environmental Observatory (BEO), Utqiagvik, Alaska. A different plant species was targeted each year, over four experimental years from 2017–2021. Each year, five warming chambers and paired ambient control plots were deployed from around the time of snowmelt in mid-June through to mid-September. Average seasonal warming of 3–4°C was achieved using Zero Power Warming (ZPW) chambers (Lewin et al, 2017).The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Cooling charge determination for a local thermal management system

A system and method for estimating the state of charge, i.e., the amount of cooling remaining, in a latent heat thermal storage module of a local environmental control system. The local environmental control system stores cooling energy at night by freezing a latent heat phase change material (PCM), and then releasing the stored cooling energy by removing heat from room air flowing through the PCM during the day. Differential temperature sensors are used measure the air temperature drop across the latent thermals storage unit and the fan voltage, which is directly related to fan air flow rate, is sampled. A controller is programmed to determine the charging state from these measurements.

Khalifa, H. Ezzat↗

Controlling radioisotope proportions when randomly sampling from Dirichlet distributions in PyRIID

As machine learning models for radioisotope quantification become more powerful, likewise the need for high-quality synthetic training data grows as well. For problem spaces that involve estimating the relative isotopic proportions of various sources in gamma spectra it is necessary to generate training data that accurately represents the variance of proportions encountered. In this report, we aim to provide guidance on how to target a desired variance of proportions which are randomly when using the PyRIID Seed Mixer, which samples from a Dirichlet distribution. We provide a method for properly parameterizing the Dirichlet distribution in order to maintain a constant variance across an arbitrary number of dimensions, where each dimension represents a distinct source template being mixed. We demonstrate that our method successfully parameterizes the Dirichlet distribution to target a specific variance of proportions, provided that several conditions are met. This allows us to follow a principled technique for controlling how random mixture proportions are generated which are then used downstream in the synthesis process to produce the final, noisy gamma spectra.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Measurement of the e + e − → π + π − π 0 cross section in the energy range 0.62–3.50 GeV at Belle II

We report a measurement of the e + e − → π + π − π 0 cross section in the energy range from 0.62 to 3.50 GeV using an initial-state radiation technique. We use an e + e − data sample corresponding to 191 fb − 1 of integrated luminosity, collected at a center-of-mass energy at or near the ϒ ( 4 S ) resonance with the Belle II detector at the SuperKEKB collider. Signal yields are extracted by fitting the two-photon mass distribution in e + e − → π + π − π 0 γ events, which involve a π 0 → γ γ decay and an energetic photon radiated from the initial state. Signal efficiency corrections with an accuracy of 1.6% are obtained from several control data samples. The uncertainty on the cross section at the ω and ϕ resonances is dominated by the systematic uncertainty of 2.2%. The resulting cross sections in the 0.62–1.80 GeV energy range yield a μ 3 π = [ 48.91 ± 0.23 ( stat ) ± 1.07 ( syst ) ] × 10 − 10 for the leading-order hadronic vacuum polarization contribution to the muon anomalous magnetic moment. This result differs by 2.5 standard deviations from the most precise current determination. Published by the American Physical Society 2024

Adachi, I. (ORCID:0000000322870173)↗

High Throughput In-Line Coating Metrology Development for Solid Oxide Fuel Cell Manufacturing

Coatings play key roles in solid oxide fuel cell (SOFC) stack durability. For example, diffusion barrier coatings on Cr-containing interconnect and balance of plant (BOP) components protect electrodes from Cr poisoning over the long operational lifetimes (>10,000 hours) of the fuel cell stack. Common defects in coatings, such as cracks, pinholes, and porosity, result in a failure to protect the electrodes, resulting in shorter operational lifetime and thus higher cost. It is very unlikely, even in the best coating process, that all these defects can be mitigated, hence identifying critical defects in parts, and removing defective parts from production before they can damage the stack, becomes paramount. Furthermore, these quality control techniques must be operational in the production/assembly line (in-line), i.e., high throughput and non-destructive, and cost effective. Redox Power Systems, LLC (Redox) together with the National Renewable Energy Laboratory (NREL) developed much needed high throughput, in-line metrology techniques for protective coatings. The overall goal of the project is to lower cost while increasing robustness, reliability, and endurance of SOFC stacks. To accomplish this, we had several objectives, including: to identify key coating and substrate defects that lead to coating failure through the use of detailed characterization methods (e.g., microscopy, XRD, EDS, electrochemistry); to assess capabilities of in-line metrology techniques, e.g., optical profilometry (Redox) and thermography (NREL), to probe these defects, or evidence thereof; demonstrate long-term performance of “defect-free” protective coatings, as identified by in-line metrology, in solid oxide fuel cell (SOFC) stack operation. In the first part of this project, the ability to identify key defects expected to lead to coating and SOFC degradation using in-line metrology tools were evaluated. Coated interconnect samples with controlled defect types and populations were tested under conditions similar to SOFC operation, followed by detailed post-test analysis to reveal the defects responsible for observed degradation. In the second part of the project, the optimal in-line metrology techniques and methodologies were used to map the defect distribution in full-size interconnects with critical defects intentionally allowed to exist in some cases. These interconnects underwent SOFC testing for extended periods (up to ~3,000 hours) followed by post-test analysis to evaluate the effectiveness of in-line metrology techniques in mitigating MCO coating related degradation. Key accomplishments in this project included the following: Demonstrated ASR of < 0.05 ohm-cm 2 at 650 °C for 1,000 hours with low defect (determined by in-line metrology) interconnect samples (average ASR=37 milliohms-cm 2 after over 1,000 hours). Demonstrated that low defect coatings on interconnects (as screened using in-line metrology) have low volatilization of chromium at ~650 °C for 1,000 hours as detected using Cr-getter material (< 5 at% increase above baseline); 1022 hour duration tests under humidified, elevated temperature (750 °C rather than 650 °C) compared a base case against different coating thicknesses. Demonstrated capability to identify initial key defects of interest with in-line metrology techniques using up to 8 cm by 10 cm having coatings with and without intentional defects of interest using thermal imaging and optical profilometry. Correlated key defects identified using metrology techniques with observed coating performance (e.g., ASR and Cr volatility). Conducted several 4 cm by 4 cm cell tests using MCO-interconnects that were pre-screened using some of the metrology techniques developed in the project (e.g., optical profilometry). An analysis of ASR measurements were able to show that defect-free coatings resulted in the anticipated performance in the cell tests.

01 COAL, LIGNITE, AND PEAT↗

Controlled Solidification of Liquids Within the SFC Primary Phase Field of the “Fe{sub 2}O{sub 3}”-CaO-SiO{sub 2} System in Air

To provide fundamental information on the phases and microstructures formed during sintering, a liquid with a bulk composition within the silico ferrite of calcium (SFC) primary phase field in the ternary “Fe{sub 2}O{sub 3}”-CaO-SiO{sub 2} system in air was solidified at a controlled rate. Samples of a bulk composition with a CaO/SiO{sub 2} ratio of 4.00 and 69.24 wt pct Fe{sub 2}O{sub 3}, were cooled from 1623 K (1350 °C) at 2 K/s, with samples quenched at temperatures between 1513 K (1240 °C) to 1453 K (1180 °C). The silico ferrite of calcium and aluminium I (SFCA-I) and Ca{sub 7.2}Fe{sub 0.8}{sup 2+}Fe{sub 30}{sup 3+}O{sub 53} phases were observed to form an intergrowth (‘SFC-I’) rather than the anticipated SFC phase. Solidification was found to occur in three stages, Liquid + ‘SFC-I’, Liquid + ‘SFC-I’ + C{sub 2}S + CF{sub 2}, and C{sub 2}S + CF{sub 2} + CF, where C{sub 2}S denotes dicalcium silicate, CF denotes calcium ferrite and CF{sub 2} denotes calcium diferrite. The phases formed and the solidification sequence differ from those predicted under equilibrium and Scheil–Gulliver Cooling. Although not directly applicable to industrial operations, this research clearly shows that the formation of both the SFCA and SFCA-I phase in iron ore sinters is controlled by kinetic processes rather than equilibrium conditions.

36 MATERIALS SCIENCE↗

Competing Magnetic Interactions and Field-Induced Metamagnetic Transition in Highly Crystalline Phase-Tunable Iron Oxide Nanorods

The inherent existence of multi phases in iron oxide nanostructures highlights the significance of them being investigated deliberately to understand and possibly control the phases. Here, the effects of annealing at 250 °C with a variable duration on the bulk magnetic and structural properties of high aspect ratio biphase iron oxide nanorods with ferrimagnetic Fe3O4 and antiferromagnetic α-Fe2O3 are explored. Increasing annealing time under a free flow of oxygen enhanced the α-Fe2O3 volume fraction and improved the crystallinity of the Fe3O4 phase, identified in changes in the magnetization as a function of annealing time. A critical annealing time of approximately 3 h maximized the presence of both phases, as observed via an enhancement in the magnetization and an interfacial pinning effect. This is attributed to disordered spins separating the magnetically distinct phases which tend to align with the application of a magnetic field at high temperatures. The increased antiferromagnetic phase can be distinguished due to the field-induced metamagnetic transitions observed in structures annealed for more than 3 h and was especially prominent in the 9 h annealed sample. Our controlled study in determining the changes in volume fractions with annealing time will enable precise control over phase tunability in iron oxide nanorods, allowing custom-made phase volume fractions in different applications ranging from spintronics to biomedical applications.

36 MATERIALS SCIENCE↗

Shade signals alter the expression of circadian clock genes in newly-formed bioenergy sorghum internodes

Stem internodes of bioenergy sorghum inbred R.07020 are longer at high plant density (shade) than at low plant density (control). Initially, the youngest newly-formed subapical stem internodes of shade-treated and control plants are comparable in length. However, full-length internodes of shade-treated plants are three times longer than the internodes of the control plants. To identify the early molecular events associated with internode elongation in response to shade, we analyzed the transcriptome of the newly-formed internodes of shade-treated and control plants sampled between 4 and 6 hr after the start of the light period (14 hr light/10 hr dark). Sorghum genes homologous to the Arabidopsis shade marker genes ATHB2 and PIL1 were not differentially expressed. The results indicate that shade signals promote internode elongation indirectly because sorghum internodes are not illuminated and grow while enclosed with leaf sheaths. Sorghum genes homologous to the Arabidopsis morning-phased circadian clock genes LHY, RVE, and LNK were downregulated and evening-phased genes such as TOC1, PRR5, and GI were upregulated in young internodes in response to shade. We hypothesize that a change in the function or patterns of expression of the circadian clock genes is the earliest molecular event associated with internode elongation in response to shade in bioenergy sorghum. Increased expression of CycD1, which promotes cell division, and decreased expression of cell wall-loosening and MBF1-like genes, which promote cell expansion, suggest that shade signals promote internode elongation in bioenergy sorghum in part through increasing cell number by delaying transition from cell division to cell expansion.

54 ENVIRONMENTAL SCIENCES↗

Towards high-precision parton distributions from lattice QCD via distillation

We apply the Distillation spatial smearing program to the extraction of the unpolarized isovector valence PDF of the nucleon. The improved volume sampling and control of excited-states afforded by distillation leads to a dramatically improved determination of the requisite Ioffe-time Pseudo-distribution (pITD). The impact of higher-twist effects is subsequently explored by extending the Wilson line length present in our non-local operators to one half the spatial extent of the lattice ensemble considered. The valence PDF is extracted by analyzing both the matched Ioffe-time Distribution (ITD), as well as a direct matching of the pITD to the PDF. Through development of a novel prescription to obtain the PDF from the pITD, we establish a concerning deviation of the pITD from the expected DGLAP evolution of the pseudo-PDF. The presence of DGLAP evolution is observed once more following introduction of a discretization term into the PDF extractions. Observance and correction of this discrepancy further highlights the utility of distillation in such structure studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The physics of desiccation cracks 2: Modeling and prediction of the crack patterns

Here, this paper extends the experimental results of our companion paper by modeling and predicting the onset and pattern formation of desiccation cracks in geomaterials (Ruoyu et al., 2023). Thin-layer silt samples in controlled atmospheric conditions were tested to obtain the surface strain maps with the digital image correlation (DIC) method during dehydration. Support experiments, including consolidation and displacement-controlled triaxial tests, were conducted for the properties of geomaterials. These experimental results were used to validate a viscoplastic theoretical model by comparing cracks and singularities locations that distribute following the Cnoidal wave pattern. A critical value exists in the viscoplastic model that determines the number of singularities, which accurately predicts the number of cracks in the experiment. The critical values contain two crucial parameters: rate sensitivity and λ. Rate sensitivity describes the rate-dependent stress–strain relation, while λ shows the ability of pore pressure redistribution under the external mechanical loading rate. These results provide a new view to analyze the desiccation cracks considering the rate-dependent viscoplasticity.

42 ENGINEERING↗

Development of a versatile, high-temperature, high-throughput ion irradiation system

We report ion irradiation has long been used as a surrogate for neutron irradiation experiments which are comparatively slower, more expensive, and often can leave materials activated and difficult to handle post-irradiation. However, with the growing use of combinatorial synthesis methods and high-throughput, automated characterization techniques, ion irradiation facilities will need to upgrade their capabilities to be able to accommodate these new material platforms. To address this challenge, a high-throughput ion irradiation system has been developed at the University of Wisconsin-Madison Ion Beam Laboratory (IBL). The system is capable of accommodating large sample arrays of arbitrary geometries and heat samples individually using an infrared laser to enable high-temperature, high-throughput ion irradiation while mitigating the risk of annealing the irradiation damage in nearby samples. To control stage motion, laser power, and beam current measurement, among other tasks, the program Chronos was developed which enables automated high-throughput irradiation experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

CO 2 –Brine Contact Angle Measurements on Navajo, Nugget, Bentheimer, Bandera Brown, Berea, and Mt. Simon Sandstones

In this study, contact angles were measured for CO 2 bubbles on six different sandstones (Navajo, Nugget, Bentheimer, Bandera Brown, Berea, and Mt. Simon) that could potentially represent properties of CO 2 storage depositional environments. The impacts of pressure and temperature were studied by focusing on the CO 2 phase behavior in three different scenarios: gaseous, liquid, and supercritical conditions. Despite controlling the sample preparation and cleanliness, CO 2 -brine equilibration conditions, and pressure and temperature, there were inconsistencies in contact angle trends that could largely be attributed to natural sample heterogeneity resulting from localized variations in topography, surface roughness, and mineral composition across the surface. Despite these variations, the analysis of 298 measurements from this study showed that 81% of the contact angles were below 40° representing strongly water-wet to (moderately) water-wet behaviors. Also, 17.3% of the measurements were between 40° and 60° (weakly water-wet) and primarily belonged to small CO 2 bubbles (< 500 μm) that were heavily impacted by localized surface heterogeneity in natural sandstone samples. Additionally, only 1.7% of the measurements had angles greater than 60° and could be classified as CO 2 -wet behavior, but these measurements occurred on extremely small bubbles (100-200 μm) that were very dependent on localized surface heterogeneity. While some of the detailed physics of the CO 2 -brine-sandstone interface is still poorly understood, from an application standpoint, the sandstones of this study were best characterized as strongly water-wet to (moderately) water-wet.

03 NATURAL GAS↗

Injection Parameters That Promote Branching of Hydraulic Cracks

Fluid injection into rock formations can either produce complex branched hydraulic fractures, create simple planar fractures, or be dominated by porous diffusion. Currently, the optimum injection parameters to create branched fractures are unknown. We conducted repeatable hydraulic fracturing experiments using analog-rock samples with controlled heterogeneity to quantify the fluid parameters that promote fracture branching. A large range of injection rates and fluid viscosities were used to investigate their effects on induced fracture patterns. Paired with a simple analytical model, our results identify the threshold at which fracture transitions from an isolated planar crack to branched cracks when closed natural fractures exist. These results demonstrate that this transition can be controlled by injection rate and fluid viscosity. In relation to the field practices, the present model predicts slickwater and lower viscosity fluid injections promote fracture branching, with the Marcellus shale used as an example.

58 GEOSCIENCES↗

Development of an ERT‐Based Framework for Bentonite Buffers Monitoring From Laboratory Tests: 1. Characterizing Thermal–Hydrological–Mechanical Processes

Abstract Bentonite clay is widely used in engineered barrier systems for the permanent disposal of high‐level radioactive waste due to its low permeability, high swelling capacity, and thermal stability. However, the complex thermal‐hydrological‐mechanical (THM) processes induced by heating from decaying radioactive waste and hydration from surrounding rock can lead to heterogeneous changes that are difficult to measure and predict. This study develops an Electrical Resistivity Tomography (ERT)‐based framework for monitoring THM processes, progressing from sample‐scale to bench‐scale tests, to inform field‐scale applications. Sample‐scale tests analyzed small bentonite samples under controlled variations in water content, temperature, and porosity to establish fundamental resistivity relationships. Bench‐scale tests involved larger bentonite columns subjected to heating (up to 200°C) and hydration under controlled pressure, simulating repository conditions. ERT measurements, complemented by X‐ray CT imaging, temperature monitoring, and tracing sensors, revealed coupled THM processes, such as hydration‐induced compression, swelling, and thermal gradients, leading to complex resistivity patterns. The results demonstrate the potential of ERT for capturing THM‐induced resistivity changes, though challenges remain in upscaling and quantitative analysis. This study evaluates laboratory test capabilities and proposes future improvements for understanding THM‐induced resistivity responses. A conceptual framework for ERT implementation in field‐scale monitoring is presented, synthesizing findings from both scales and exploring how ERT data can inform long‐term modeling and reduce prediction uncertainties. Overall, this ERT‐based framework offers a robust method for monitoring bentonite buffers, aiding in early issue detection and supporting the safe long‐term disposal of radioactive waste in geological repositories, while highlighting the need for future development. Plain Language Summary Bentonite clay is crucial in engineered barrier systems (EBS) for containing high‐level radioactive waste due to its ability to absorb water, swell, seal and remain stable under high temperatures. When bentonite absorbs water and heats up from radioactive decay, it experiences complex changes in its physical and mechanical properties. Understanding these changes is important for ensuring the long‐term safety and effectiveness of EBS. This study used Electrical Resistivity Tomography (ERT), a non‐invasive method that measures electrical conductivity to monitor these changes during laboratory experiments. The ERT data revealed significant variations in resistivity corresponding to changes in water content, temperature, and density, providing detailed spatial and temporal insights into the behavior of bentonite. These findings enhance our ability to predict the long‐term performance of bentonite barriers, ensuring the safe containment of radioactive waste. By improving our understanding of bentonite's behavior, this research supports the development of more reliable and effective barrier systems for radioactive waste disposal, protecting the environment and public health. Key Points ERT monitoring was employed to capture resistivity changes in bentonite during controlled heating and hydration experiments, providing insights into THM processes ERT data reveal significant resistivity changes correlated with water content, temperature, and mechanical effects, enhancing the understanding of THM dynamics in bentonite This study explores the potential of the framework for application in field‐scale EBS monitoring, emphasizing the need for integrating additional geophysical methods for comprehensive subsurface imaging

Chen, Hang↗

Automated calculation and convergence of defect transport tensors

Defect diffusion is a key process in materials science and catalysis, but as migration mechanisms are often too complex to enumerate a priori, calculation of transport tensors typically have no measure of convergence and require significant end-user intervention. These two bottlenecks prevent high-throughput implementations essential to propagate model-form uncertainty from interatomic interactions to predictive simulations. In order to address these issues, we extend a massively parallel accelerated sampling scheme, autonomously controlled by Bayesian estimators of statewide sampling completeness, to build atomistic kinetic Monte Carlo models on a state-space irreducible under exchange and space group symmetries. Focusing on isolated defects, we derive analytic expressions for drift and diffusion coefficients, providing a convergence metric by calculating the Kullback–Leibler divergence across the ensemble of diffusion processes consistent with the sampling uncertainty. The autonomy and efficacy of the method is demonstrated on surface trimers in tungsten and Hexa-interstitials in magnesium oxide, both of which exhibit complex, correlated migration mechanisms.

36 MATERIALS SCIENCE↗

Shock compression of cemented tungsten carbides to 100 GPa: Structure of shock waves, Hugoniot relations, and strength

Plate impact experiments are conducted on cemented tungsten carbides (WC) with a 3.7 and 6.0 wt. % cobalt binder to better understand their dynamic, high-pressure response to 100 GPa. The measured wave profiles show propagation of steady structured waves. Standard impedance matching procedures are used to determine the Hugoniot relations in the shock velocity–particle velocity (U s –v p ) and Hugoniot stress–specific volume (P–V/V o ) planes. The Hugoniot elastic limit of the samples is controlled by ductility of the Co binder and is determined to be 4.45 ± 0.29 GPa for cemented WC with 3.7 wt. % cobalt and 3.72 ± 0.24 GPa for a 6.0 wt. % cobalt binder. Both grades show a non-linear U s –v p relationship depending on whether the particle velocity is in the strength dominated or hydrodynamic regime. In the strength dominated regime, a non-linear decrease in U s is observed as v p increases from ambient to the material’s hydrodynamic limit. In the hydrodynamic regime, the U s –v p Hugoniot is linear and is determined to be U s = 4.97(±0.006)+1.446(±0.018)v p km/s for WC with 3.7 wt. % Co and U s = 4.93(±0.006)+1.463(±0.017)v p km/s for 6 wt. % Co. Both WC grades indicate shear-stress hardening with mean stress immediately after yield, followed by pressure softening, and then a sharp fall in stress carrying capacity as the mean stress is increased to ≈70 GPa (hydrodynamic limit) and beyond. This behavior is in contrast to pure WC ceramics, which show continued shear-stress hardening with mean stress to ≈80 GPa.

36 MATERIALS SCIENCE↗