Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “APPROXIMATIONS”

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

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

At least 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A fourth-order phase-field fracture model: Formulation and numerical solution using a continuous/discontinuous Galerkin method

Modeling crack initiation and propagation in brittle materials is of great importance to be able to predict sudden loss of load-carrying capacity and prevent catastrophic failure under severe dynamic loading conditions. Second-order phase-field fracture models have gained wide adoption given their ability to capture the formation of complex fracture patterns, e.g. via crack merging and branching, and their suitability for implementation within the context of the conventional finite element method. Higher-order phase-field models have also been proposed to increase the regularity of the exact solution and thus increase the spatial convergence rate of its numerical approximation. However, they require special numerical techniques to enforce the necessary continuity of the phase field solution. In this paper, we derive a fourth-order phase-field model of fracture in two independent ways; namely, from Hamilton’s principle and from a higher-order micromechanics-based approach. The latter approach is novel, and provides a physical interpretation of the higher-order terms in the model. In addition, we propose a continuous/discontinuous Galerkin (C/DG) method for use in computing the approximate phase-field solution. This method employs Lagrange polynomial shape functions to guarantee -continuity of the solution at inter-element boundaries, and enforces the required regularity with the aid of additional variational and interior penalty terms in the weak form. Finally, the phase-field equation is coupled with the momentum balance equation to model dynamic fracture problems in hyper-elastic materials. Two benchmark problems are presented to compare the numerical behavior of the C/DG method with mixed finite element methods.

42 ENGINEERING↗
DOE OSTI2025

Osmotic control of the spacing of parallel shear cracks in shale growing subcritically in geologic past

The geological genesis of natural cracks in sedimentary rocks such as shale is a problem that needs to be understood to improve the technology of hydraulic fracturing as well as deep sequestration of harmful fluids. Why are the vertical natural cracks roughly parallel and equidistant, and why is the spacing roughly 10 cm rather than 1 cm or 100 cm? Fracture mechanics of critical cracks cannot answer this question. Neither can the material heterogeneity. The growth of critical parallel cracks is impossible because the relative crack face displacements would immediately localize into one crack, leading to an earthquake. The cracks must have formed, on the tectonic time scale, by a slow growth of subcritical shear cracks governed by the Charles-Evans law. The idea advanced here is that what controls the crack spacing is the balance between the reduction, due to shear dilatancy, of the concentration of ions such as Na + and Cl - in each fracture process zone (PFZ), which decelerates the cracks, and the restoration of ion concentration by diffusion of ions from the space between the cracks into the FPZ. This diffusion of water is driven mainly by the osmotic pressure gradient, which offsets the deceleration and depends strongly on the crack spacing. A simple analytical solution of the steady state is rendered possible by approximating the ion concentration profiles between adjacent cracks by parabolic arcs. Applying this theory to Woodford shale yields the approximate crack spacing of 10 cm, which is realistic. Furthermore, the stability of unlimited parallel mode II frictional crack growth is proven by examining the second variation of the free energy. Water concentration drop in the FPZ due to shear dilatancy and its restoration by water diffusion from the inter-crack space have similar effect, although probably much weaker.

42 ENGINEERING↗
DOE OSTI2024

Martensitic transformation induced strength-ductility synergy in additively manufactured maraging 250 steel by thermal history engineering

Maraging steels are known for their exceptional strength but suffer from limited work hardening and ductility. Here, in this study, we report an intermittent printing strategy to tailor the microstructure and mechanical properties of maraging 250 steel via tuning the thermal history during wire-arc directed energy deposition. By introducing a dwell time between adjacent layers, the maraging 250 steel is cooled below the martensite start temperature, triggering thermally-driven martensitic transformation during the printing process. Thermal cycling during subsequent layer deposition results in the formation of reverted austenite which shows a refined microstructure and induces elemental segregation between martensite and reverted austenite. The Ni enrichment in the austenite promotes stabilization of the reverted austenite upon cooling to room temperature. The reverted austenite is metastable during deformation, leading to strain-induced martensitic transformation under loading. Specifically, a 3 min interlayer dwell time produces a maraging 250 steel with approximately 8% reverted austenite, resulting in improved work hardening via martensitic transformation induced plasticity during deformation. Meanwhile, the higher cooling rate and refined prior austenite grains lead to substantially refined martensitic grains (by approximately fivefold) together with an increased dislocation density. With 3 min interlayer dwell time, the yield strength of the printed maraging 250 steel increases from 836 MPa to 990 MPa, and the uniform elongation is doubled from 3.2% to 6.5%. This intermittent deposition strategy demonstrates the potential to tune the microstructure of maraging steels for achieving strength-ductility synergy by engineering the thermal history during additive manufacturing.

Additive manufacturing↗
DOE OSTI2020

High spatial resolution thermal conductivity mapping of SiC/SiC composites

Herein, silicon carbide (SiC) ceramic matrix composites are being investigated as a new generation of fuel cladding materials due to its higher accident tolerance compared to Zircaloy. Characterization of the thermal conductivity of SiC constituents and their interfaces and interphases in SiC/SiC composites is needed as inputs to models of cladding performance. We used time-domain thermoreflectance (TDTR) to map the thermal conductivity with a spatial resolution of 2 µm. The SiC/SiC composite is comprised of Hi-Nicalon Type S fibers in a matrix made by chemical vapor infiltration (CVI matrix). The interphase material is a pyrolytic carbon/SiC multilayer. We report thermal conductivity maps of the SiC fiber and SiC matrix at temperatures of 25 °C, 90 °C, 164 °C, and 250 °C. The fiber has a uniform and isotropic thermal conductivity of 22 W m -1 K -1 at 25 °C; the thermal conductivity of fiber is approximately independent of temperature. The matrix thermal conductivity is not uniform and varies from 50 to 120 W m -1 K -1 at 25 °C; the thermal conductivity of the matrix scales approximately with the inverse of the square root of temperature, i.e., ∝ 1/T 1/2 . The thermal conductivity of the interphase at room temperature is 6 W m -1 K -1 . The calculated spatially averaged thermal conductivity of this SiC/SiC composite at 1000°C is 29 W m -1 K -1 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Evaluation of the anisotropic grain boundaries and surfaces of α -U via molecular dynamics

Alloys based on uranium-zirconium are gaining renewed interest as fuels for the Versatile Test Reactor and a number of microreactor designs. Implementing metallic fuel in reactors creates the need for robust descriptive and predictive fuel performance modeling. Here, the current state of metallic fuel performance modeling relies on empirical equations derived from historical experiments, which may be unreliable when applied outside of their temperature, power, and composition phase space. One area where such data is lacking is the irradiation behavior of alpha -U, specifically tearing and porosity formation at the early stages of irradiation. While grain boundaries likely play a key role in this fuel behavior, relatively little is known about grain boundaries in alpha-U. Thus, we evaluate the grain boundary, surface energy, and work of adhesion of alpha -U utilizing molecular dynamics. Symmetric tilt grain boundaries (STGBs) are analyzed with the tilt plane oriented along each major crystallographic axis, for a total of eighty unique grain boundaries. The effect of temperature, tilt plane, and misorientation angle on interfacial energies are analyzed. The interfacial energies typically increase with temperature and there is significant variance as a function of misorientation angle, irrespective of the tilt plane. At 500 K, the average surface energy (1.23 J/m 2 ) is approximately 1.5 times the grain boundary energy (0.79 J/m 2 ), and the work of adhesion is approximately twice the grain boundary energy (1.68 J/m 2 ). Orientations for the likely formation of twins and likely failure planes are identified.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Fission product distribution in irradiated safety-tested and as-irradiated AGR-2 TRISO particles

In this work, two tristructural isotropic (TRISO)-coated nuclear fuel particles were examined by electron probe microanalysis (EPMA) as part of the Advanced Gas Reactor program. The compacts’ average irradiation temperatures ranged from approximately 1260 to 1290 °C. One particle was examined in the as-irradiated condition, while the other was subject to 1600 °C post-irradiation safety testing. This study was undertaken to test a newly-developed EPMA technique to determine fission product masses in TRISO particles on a layer-by-layer basis, and to compare fission product distributions between an as-irradiated and safety-tested particle. Fission product concentration profiles were collected along two radii in each particle, with measured concentrations used to compute the fission product mass in each TRISO particle layer. These measured masses were then compared to those predicted from ORIGEN modeling calculations. Data collected from these measurements show that for these two particles, masses determined via EPMA were within ± 20% of the calculated masses for the rare-earth elements, Mo, Zr, Cs, I, and Pd. Elements that tend to be less homogeneously distributed include Sr, Te, Eu, Ag, and possibly Ba. Measured Ag masses differed by more than 40% from the calculated mass. Lanthanides other than Eu remain primarily within the fuel kernel in the as-irradiated particle but in the safety-tested particle these element masses were divided approximately equally between the kernel and kernel periphery. In both particles, the majority of Sr and Eu accumulated in the carbon-rich kernel periphery, although in the safety-tested particle, Sr and Eu accumulated farther from the fuel kernel than occurred with irradiation alone. A greater mass fraction of mobile elements, such as Cs and I accumulated in the buffer and IPyC in the safety-tested particle as compared to the as-irradiated particle. When fully developed and tested, this mass balance approach to TRISO particle analysis has the potential to provide insight into fuel behavior.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Comparison of unirradiated and irradiated AGR-2 TRISO fuel particle oxidation response

The silicon carbide (SiC) coating in a tristructural isotropic (TRISO) particle acts as a barrier to fission product release during reactor operation and accident scenarios. Oxidation and subsequent failure of the SiC layer during a rare air ingress event is a proposed mechanism for fission product release in a high-temperature gas-cooled reactor (HTGR). Although previous oxidation studies have analyzed unirradiated TRISO particle response, this study compared the oxidation behavior of irradiated and unirradiated TRISO particles from the second Advanced Gas Reactor Fuel Development and Qualification Program irradiation experiment (AGR-2). Particles with exposed SiC were subjected to six varying oxidizing tests in the Furnace for Irradiated TRISO Testing (FITT), examined for failure fraction with the Irradiated Microsphere Gamma Analyzer (IMGA) and characterized with focused ion beam and scanning/transmission electron microscopy techniques to analyze the oxide layer. Uncorrelated unirradiated particle failures throughout the series of exposures suggests that external factors inherent to the experiment increased particle failure sensitivity. However, irradiated particle observations indicated an increased failure response at 400 h 1400 °C in both 2% and 21% O 2 atmospheres above failure associated with external factors. Oxide thickness measurements after 400 h at 1400 °C revealed a greater oxidation rate than predicted by parabolic growth, which was attributed to the increased complexity of the oxide structure at longer exposure times. Altering the atmosphere from 21% to 2% O 2 reduced the average oxide thickness by approximately 12%–14% in both irradiated and unirradiated particles at 400 h 1400 °C. Altogether, the minor variations observed between irradiated and unirradiated particles in this study led to the conclusion that unirradiated TRISO particles can be used to approximate irradiated TRISO oxidation kinetics.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Modelica-based system modeling for studying control-related faults in chiller plants and boiler plants serving large office buildings

System modeling is critical when studying operation faults in chiller plants and boiler plants. However, current fault models have difficulties faithfully representing the operation of chiller plants and boiler plants under the effects of those faults, especially the control-related ones. In this study, we present a systematic method to develop high-fidelity models for approximating the behaviors of chiller plants and boiler plants under faulty conditions. Compared to existing ones, the resulting fault models have two advantages: first, they better characterize the dynamic patterns in the system operation. In those models, control architecture and control logic are faithfully implemented. Thus, they can be used to study control-related faults, such as incorrect staging control due to sensor bias and mistuned feedback control. Second, they are readily extensible and can support large-scale investigations to explore different faulty conditions/scenarios. Those models are established in a hierarchical structure while modules in each layer can be redeclared and parameterized at upper layers. In such a case, modifications to the models can be realized through model modifiers and the process can be easily streamlined with scripts. We applied the proposed models in a comprehensive fault impact evaluation of the 13 control-related faults of chiller and boiler plants. In this evaluation, the proposed plant model is coupled with the EnergyPlus thermal load model to study the impact of various faulty scenarios. Based on the evaluation results, we identified the faults that have the most significant impacts on the operation of the chiller and boiler plants, respectively. We also found that the relationship between the impacts of the studied faults and the severity level of the faults are highly non-linear

Huang, Sen↗
DOE OSTI2023

An ultrahigh-resolution E3SM land model simulation framework and its first application to the Seward Peninsula in Alaska

The availability of supercomputers and state-of-science datasets has made it possible to conduct large-scale land simulations at an ultrahigh-resolution. This study reported a computational framework for land surface simulation using the E3SM land model (ELM) at an unprecedented resolution (1 km x 1 km gridcell). The ultrahigh-resolution ELM (uELM) simulation framework includes three parts: (1) high-resolution atmospheric forcing and surface properties dataset generation, (2) massive gridcell-based simulation, and (3) large-scale simulation results analysis. Additionally, we implemented the uELM simulation framework and completed the first 1 km x 1 km terrestrial ecosystem simulation (from 1850 to 2014) over the Seward Peninsula in Alaska (78,000 km 2 ). The experiment contained two phases: a spin-up simulation and a transient simulation, and required five weeks of calculations using 320 cores in a 44-node Linux HPC computer. It created approximately 1.3 TB of data from the transient simulation alone (1850 - present). We selected sample results (monthly and daily simulation outputs) to illustrate the temporal and spatial variations of several variables in high-latitude Arctic ecosystems’ water, energy, and carbon cycles. At last, we summarized the lessons learned and proposed new developments for full-scale uELM simulations over the entire North American continent (approximately 22,000,000 km 2 ).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Application of Koopman operator for model-based control of fracture propagation and proppant transport in hydraulic fracturing operation

This work explores the application of the recently developed Koopman operator approach for model identification and feedback control of a hydraulic fracturing process. Controlling fracture propagation and proppant transport with precision is a challenge due in large part to the difficulty of constructing approximate models that accurately capture the characteristic moving boundary and highly-coupled dynamics exhibited by the process. Koopman operator theory is particularly attractive here as it offers a way to explicitly construct linear representations for even highly nonlinear dynamics. The method is data-driven and relies on lifting the states to an infinite-dimensional space of functions called observables where the dynamics are governed by a linear Koopman operator. Here this work considers two problems: (a) fracture geometry control, and (b) proppant concentration control. In both cases, an approximate linear model of the corresponding dynamics is constructed and used to design a model predictive controller (MPC). The manuscript shows that in the case of highly nonlinear dynamics, as observed in the proppant concentration, use of canonical functions in the observable basis fails. In such cases, a priori system knowledge can be leveraged to choose the required basis. The numerical experiments demonstrate that the Koopman linear model shows excellent agreement with the real system and successfully achieves the desired target values maximizing the oil and gas productivity. Additionally, due to its linear structure, the Koopman models allow convex MPC formulations that avoid any issues associated with nonlinear optimization.

42 ENGINEERING↗
DOE OSTI2024

O'Hare Airport roadway traffic prediction via data fusion and Gaussian process regression

This study proposes an approach of leveraging information gathered from multiple traffic data sources at different resolutions to obtain approximate inference on the traffic distribution of Chicago's O'Hare Airport area. Specifically, it proposes the ingestion of traffic datasets at different resolutions to build spatiotemporal models for predicting the distribution of traffic volume on the road network. Due to its good adaptability and flexibility for spatiotemporal data, the Gaussian process (GP) regression was employed to provide short-term forecasts using data collected by loop detectors (sensors) and supplemented by telematics data. The GP regression is used to make predictions of the distribution of the proportion of sensor data traffic volume represented by the telematics data for each location of the sensors. Consequently, the fitted GP model can be used to determine the approximate traffic distribution for a testing location outside of the training points. Policymakers in the transportation sector can find the results of this work helpful for making informed decisions relating to current and future transportation conditions in the area.

42 ENGINEERING↗
DOE OSTI2023

Analysis of backwash settings to maximize net water production in an engineering-scale ultrafiltration system for water reuse

Ultrafiltration (UF) has been widely utilized as water pretreatment for different applications especially in water reuse. The UF system operation is characterized by a filtration phase, where particles accumulate on the membrane surface resulting in an increase in the transmembrane pressure (TMP) and a cleaning phase, where foulants are removed through cleaning cycles including physical backwash and chemical-enhanced backwash (CEB). In this study, data from an engineering-scale UF system treating reclaimed wastewater were used to assess the impact of backwashing on the filtration process. TMP backwash trigger, backwash duration, and CEB frequency were purposely varied for a cycle-by-cycle investigation on the net water production, water recovery, initial operating TMP, and filtration cycle duration. As the TMP backwash trigger was varied between 62 and 145kPa, the maximum net water production (63 m 3 /d) was achieved at 103kPa and water recovery remained relatively constant at approximately 92%. Backwash durations of 45, 65, and 85s were performed where both net water production and water recovery yielded similar results (~63 m 3 /d and ~ 91%) compared with 103kPa TMP backwash trigger. The CEB frequency was also lowered from one every three backwashes (1/3) to 1/6 and 1/12 and resulted in decreased net water production and water recovery while the initial TMP increased. Interestingly, the total number of CEBs remained approximately constant regardless of their frequency. In conclusion, results suggest that CEB is an important fouling control process to maximize water production.

42 ENGINEERING↗
DOE OSTI2026

FIB-ToF-SIMS characterization of irradiated U-10Zr

Post-irradiation examination (PIE) is critical for the performance assessment and qualification of nuclear fuels. Secondary ion mass spectrometry (SIMS) is a powerful materials characterization technique that allows for elemental and isotopic mapping with a depth resolution greater than EDS and EPMA. However, it has not yet been applied to PIE of metallic nuclear fuel. Here, in this work, we characterize an fast neutron spectrum irradiated U-10Zr fuel sample using a time-of-flight SIMS (ToF-SIMS) system connected to a FIB/SEM system, which allows for flexible sample analysis compared to a dedicated ToF-SIMS instrument. Analysis of the resulting hyperspectral micrograph data was aided by the development of an unsupervised machine learning (ML) algorithm that iterates on existing methods to segment the 3D micrographic datasets based on the similarity of mass spectra. The results showed that the FIB-ToF-SIMS instrument was potentially capable of spatially resolving closed fission gas bubbles in 3D by continued ion sputtering of the analyzed volume. Additionally, the ML algorithm proved useful in revealing the chemical segregation of light fission products (those with an atomic mass between approximately 85–105 amu, such as ruthenium and rhodium) plus matrix zirconium, heavy fission products (those with an atomic mass between approximately 135–150 amu, such as the lanthanides) and uranium. Future studies are planned to conduct FIB-ToF-SIMS analysis on more irradiated U-Zr samples to study the constituent redistribution.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Electro-osmo-dialysis through nanoporous layers physically conjugated to micro-perforated ion-exchange membranes: Highly selective accumulation of trace coions

This study explores theoretically and numerically the phenomenon of counter-flow accumulation of trace coions during electro-osmo-dialysis through composite membranes containing nanoporous and perforated ion-exchange layers. The combination of charged nanopores and perforations allows electroosmosis, and ion concentration polarization appears near the ion-exchange layer. Accumulation of trace coions occurs because their advective and electromigration flow velocity components have opposite signs and are locally compensating in depletion regions where electric fields are high. This local compensation of flow components can give rise to high accumulation factors (>1000 in some cases) for trace coions at some distance from the interface between the nanoporous and ion-exchange layers. Moreover, the extent of this accumulation depends on the trace-coion diffusion coefficient. The resulting selectivity between different coions might serve as a basis for environmentally friendly separation of ionic species such as synthetic peptides or rare-earth elements. For example, a periodic process might include accumulation followed by species release into smaller volumes of solution. For both rare-earth ions and peptide-like species, the accumulation selectivity between ions with small differences in mobilities depends on system parameters such as nanoporous-layer thickness, distance between the perforations, and the zeta potential of the nanopore surface. Here, the use of optimal values of these membrane parameters may enable tuning of the electro-osmo-dialysis system for separation of specific ions. Modelling reveals that despite inhomogeneous patterns of fluid and ion flows close to the nanoporous/perforated ion-exchange interface, at relatively short distances from this interface (less than half the distance between the perforations) all the flows become practically one-dimensional (1D). Similarly, the concentration profiles of trace coions also become 1D. The flow one-dimensionality at short distances from the current-polarized interface enables a simple 1D approximation that is useful for an approximate analysis of trends in this multi-parameter system.

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

Numerical simulations of membrane distillation systems with actively heated membranes

Membrane distillation (MD) is a thermal desalination process that is gaining attention for treating hypersaline brines. One recent approach uses composite membranes with a thermally conductive layer that delivers heat to the membrane–feed interface, where it drives evaporation. Though this increases single-pass recovery, it comes with the challenge of promoting lateral heat conduction through the thin membrane. Herein, we develop a 3D, computational fluid dynamics (CFD) model that simulates conjugate heat, mass, and momentum transport in the feed channel and composite membrane. We then use the CFD to verify a simpler numerical model that approximates the feed velocity field analytically. We validate the numerical model experimentally, and use it to investigate the impacts of the conductive layer thickness, feed channel geometry, and operating conditions on temperature and concentration polarization, permeate production, and specific energy consumption. Overall, we find that lateral heating increases permeate production at the expense of increased concentration polarization. In extreme cases, the concentration increases more than four-fold along the membrane surface. Furthermore, we show however, that the feed flow rate and conductive layer thickness can be tailored to mitigate concentration polarization, for only a small reduction in permeate production.

42 ENGINEERING↗
DOE OSTI2024

Effect of stoichiometry and hydration level on water domain size and transport in poly(aryl piperidinium) alkaline anion-exchange membranes

Alkaline water electrolysis holds promise in decarbonizing the global economy by enabling renewable hydrogen production with non-precious group metal catalysts. Anion exchange membranes are an important component of alkaline water electrolyzers and would ideally be durable while allowing for high hydroxide conductivity. The poly(aryl piperidinium) (PAP) class of polymers has attracted recent interest due to their good mechanical robustness and high ionic conductivity. Here, in this work, we perform atomistic molecular dynamics (MD) simulations of several PAP polymers at experimentally relevant hydration levels and polymer ion exchange capacities (IECs) to gain nanoscale insight into their properties and to help elucidate the trade-offs that result from tuning the IECs through the polymer stoichiometry. Our MD-predicted macroscopic polymer properties were found to be in good agreement with experimentally available polymer swelling ratios, water-occupied volumes, X-ray scattering, and ionic conductivities. The models show that for hydration levels greater than 8H 2 O per cation a single water cluster will form that percolates through the system. The growth in water cluster size results in large polymer swelling, the creation of larger channels with widths of 7 Å or larger, and nanophase separation between the hydrophilic domains and the polymer with characteristic length scales of approximately 20–30 Å. The experimentally observed lack of a strong X-ray scattering peak at low wavevectors can be explained by a cancellation between the polymer-polymer/water-water and polymer-water correlations and not a loss in nanophase separation. The overlap in coordination environments of the hydroxide oxygen and polymer nitrogen atoms implies that vehicular diffusion between cationic groups could play a role in hydroxide transport. The polymers' hydroxide and water diffusion constants increase by approximately an order of magnitude between hydration levels of 8 and 20H 2 O per cation. However, there are diminishing returns in hydroxide diffusion constant once the IEC exceeds 2.4 meq/g.

08 HYDROGEN↗
DOE OSTI2025

Development of a high-throughput method for processing sponge-stick samples to detect viable Bacillus anthracis spores

Since the national validation of the sponge-stick based method for detection of Bacillus anthracis spores in environmental samples, there have not been focused efforts to address the low throughput nature of the method, which processes only one sample at one time. Sample processing remains a serious bottleneck for rapidly analyzing large numbers of samples expected from a biological warfare attack. Therefore, we developed a high-throughput method to simultaneously process multiple sponge-stick samples to be better prepared for rapid response and recovery after wide area anthrax incidents. In this method, sponges are placed in 50 mL tubes containing 25 mL extraction buffer and shaken to release spores, after which the suspension is recovered for analysis. Here, we determined that an additional extraction step, conducted in the same tubes with 10 mL buffer, further increased spore recovery from sponge-stick by approximately 10 %. We determined that orbital shaking and multi-tube vortexing were both more effective than reciprocating shaking for recovering spores. We conducted simultaneous processing of up to 12 sponge-stick samples and demonstrated comparable spore recovery efficiencies to the traditional low-throughput stomacher-based method (approximately 60 % recovery at 10 2 -spore level and 75 % recovery at 10 4 -spore level for both methods in three replicate experiments, P > 0.05 for two-tailed t-tests for each experiment and spore level). We also demonstrated that our high-throughput method could be integrated with Rapid Viability-Polymerase Chain Reaction (RV-PCR) analysis and could detect levels as low as 40 spores per sponge even when challenged by a PCR particulate contaminant.

Anthrax↗
DOE OSTI2026

Uncertainty guided online ensemble for non-stationary data streams in fusion science

Machine Learning (ML) is poised to play a pivotal role in the development and operation of next-generation fusion devices. Fusion data shows non-stationary behavior with distribution drifts, resulted by both experimental evolution and machine wear-and-tear. ML models assume stationary distribution and fail to maintain performance when encountered with such non-stationary data streams. Online learning techniques have been leveraged in other domains, however it has been largely unexplored for fusion applications. In this paper, we investigate online learning for continuous adaptation to drifting data streams in the prediction of Toroidal Field (TF) coils deflection at the DIII-D fusion facility. We further address the short-term performance degradation inherent to standard online learning, which arises because ground truth is unavailable at prediction time. To mitigate this issue, we propose an uncertainty-guided online ensemble framework. The method leverages the Deep Gaussian Process Approximation (DGPA) for calibrated uncertainty estimation and uses these uncertainty measures to guide a meta-algorithm that aggregates predictions from learners trained over different historical horizons. Our results show that online learning reduces prediction error by 80% compared to a static model. The online ensemble and the proposed uncertainty-guided ensemble further reduce error by approximately 6%, and 10% respectively, relative to standard single-model online learning, while also providing calibrated uncertainty estimates to support operational decision-making.

AI↗