SEARCH · Engineering Papers

Results for “mean current vector”

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.

SourceResearch connectionFrom yearThrough year

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

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Data-Driven Cyber-Attack Detection for PV Farms via Time-Frequency Domain Features

The internetworking of grid-connected power electronics converters (PECs) in photovoltaic (PV) farms has inevitably expanded the cyber-attack surfaces. Here this paper presents a comprehensive study on cyber-attack detection and diagnosis for PEC-enabled PV farms via single waveform sensor to distinguish between normal conditions, open-circuit faults, short-circuit faults, and cyber-attacks. To our knowledge, this has not been attempted before. Firstly, we propose frequency-domain magnitude-based residuals to identify short-circuit faults and a time-domain mean current vector-based feature to distinguish open-circuit faults from other threats. These features can fully reflect the specific physical characteristics of PV farms during threat duration. Secondly, unlike micro phasor measurement units (µPMU) and raw electric waveform-based methods, the proposed innovative features can address novel cyber-attacks that are excluded from the training process. Thirdly, an online hardware-in-the-loop (HIL) testbed using the OPAL-RT real-time digital simulator has verified the effectiveness. The monitoring system runs in real-time while using HIL as an operational solar farm and a National Instruments (NI) data acquisition card as the electric waveform sensor at the point of coupling.

42 ENGINEERING↗

DOE OSTI2020

Reconstruction for Liquid Argon TPC Neutrino Detectors Using Parallel Architectures.

Neutrinos are particles that interact rarely, so identifying them re-quires large detectors which produce lots of data. Processing this data with thecomputing power available is becoming more dicult as the detectors increasein size to reach their physics goals. In liquid argon time projection chambers(TPCs) the charged particles from neutrino interactions produce ionization elec-trons which drift in an electric field towards a series of collection wires, and thesignal on the wires is used to reconstruct the interaction. The MicroBooNEdetector currently collecting data at Fermilab has 8000 wires, and planned fu-ture experiments like DUNE will have 100 times more, which means that thetime required to reconstruct an event will scale accordingly. Modernization ofliquid argon TPC reconstruction code, including vectorization, parallelizationand code portability to GPUs, will help to mitigate these challenges. The liquidargon TPC hit finding algorithm within theLArSoftframework used acrossmultiple experiments has been vectorized and parallelized. This increases thespeed of the algorithm on the order of ten times within a standalone version onIntel architectures. This new version has been incorporated back intoLArSoftso that it can be generally used. These methods will also be applied to otherlow-level reconstruction algorithms of the wire signals such as the deconvolu-tion. The applications and performance of this modernized liquid argon TPCwire reconstruction will be presented.

Berkman, Sophie↗

DOE OSTI2023

A review on recent machine learning applications for imaging mass spectrometry studies

Imaging mass spectrometry (IMS) is a powerful analytical technique widely used in biology, chemistry, and materials science fields that continue to expand. IMS provides a qualitative compositional analysis and spatial mapping with high chemical specificity. The spatial mapping information can be 2D or 3D depending on the analysis technique employed. Due to the combination of complex mass spectra coupled with spatial information, large high-dimensional datasets (hyperspectral) are often produced. Therefore, the use of automated computational methods for an exploratory analysis is highly beneficial. The fast-paced development of artificial intelligence (AI) and machine learning (ML) tools has received significant attention in recent years. These tools, in principle, can enable the unification of data collection and analysis into a single pipeline to make sampling and analysis decisions on the go. There are various ML approaches that have been applied to IMS data over the last decade. Here, in this review, we discuss recent examples of the common unsupervised (principal component analysis, non-negative matrix factorization, k-means clustering, uniform manifold approximation and projection), supervised (random forest, logistic regression, XGboost, support vector machine), and other methods applied to various IMS datasets in the past five years. The information from this review will be useful for specialists from both IMS and ML fields since it summarizes current and representative studies of computational ML-based exploratory methods for IMS.

47 OTHER INSTRUMENTATION↗

DOE OSTI2026

Transplatformer: translating toxicogenomic profiles between generations of platforms

Background Transcriptomic profiling technologies have advanced the analysis of biological and toxicological responses. However, substantial differences in probe design, dynamic range, gene coverage, and preprocessing pipelines across platforms introduce artifacts that limit cross-study integration and hinder the reuse of historical datasets. We aim to develop computational methods for accurate cross-platform translation to maximize the value of legacy resources. Results We present TransPlatformer a deep learning framework for translating gene expression profiles across heterogeneous toxicogenomics platforms. TransPlatformer employs a novel attention-based architecture to map high-dimensional fold-change vectors from legacy microarray technologies to current platforms. Models are trained and evaluated using DrugMatrix, spanning three technological generations. We investigate mixed-tissue, single-tissue, and cross-tissue training paradigms and benchmark performance against multilayer perceptron and matrix-completion baselines. In mixed-tissue training, TransPlatformer achieves a greater than 50% reduction in mean absolute error (0.043 vs. 0.09) and nearly doubles Pearson correlation ( ≈ 0.71 vs. 0.37) relative to baseline methods. Importantly, TransPlatformer preserves rare but biologically meaningful over- and under-expressed signals, with mean absolute error below 0.22. Single-tissue models yield further improvements for well-represented organs, such as a 10% reduction in liver mean absolute error, while underscoring the need for data augmentation strategies in low-sample tissues.ra Conclusions TransPlatformer provides an effective and scalable computational solution for cross-platform transcriptomic translation. By enabling biologically faithful harmonization of gene expression data, the proposed approach facilitates the reuse of legacy toxicogenomics datasets, enhances downstream biomarker discovery, and supports more reproducible predictive modeling in toxicology.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Beyond the two-point correlation: Constraining primordial non-Gaussianity with density-perturbation moments

Constraining primordial non-Gaussianity (PNG) on the large-scale cosmic structure (LSS) is an important step in understanding properties of the early Universe, specifically in distinguishing between different inflationary models. Measuring PNG relies on evaluating the scale-dependent correlations in the density field. New summary statistics beyond the two- and three-point correlation functions in configuration space and their Fourier-space counterparts, the power- and bispectrum may provide increased sensitivity. We introduce a new method for extracting the PNG signal imprinted on the LSS by using the first three Gaussian moments of the normalized correlation in density perturbations, evaluated on varying distance scales. We aim to assess this method’s sensitivity to local PNG, parameterized by f NL . We performed spherical convolutions on a range of scales on dark-matter-halo simulations to measure the scale-dependent correlations in the density field. From these, we computed the first three moments and compared them to a model expectation vector, parameterized to the second power in f NL . Our method provides about 21% improvement in sensitivity to f NL with respect to using the two-point correlation function alone. Notably, we find that the second moment alone carries nearly as much constraining power as the mean, highlighting the potential of higher order statistics. Given its simplicity and efficiency, this framework is well suited for application to current and upcoming large-scale surveys such as the Dark Energy Spectroscopic Instrument (DESI).

early universe↗

DOE OSTI2026

A Machine Learning based Approach of Estimating Equivalent Circuit Model Parameters at Different SoCs of Li-ion Batteries from Voltage Relaxation

Abstract: In this study, an approach of estimating the equivalent circuit model (ECM) parameters for Li-ion batteries (LIBs) is proposed based on the voltage value at different intervals while relaxing the LIB after discharge. The typical approach for estimating ECM parameters of a LIB is to conduct electrochemical impedance spectroscopy (EIS) measurements at different frequencies and fit them to a predefined circuit model, which requires additional measuring arrangements and specialized devices. The proposed methodology utilizes four different voltages at 0s, 60s, 360s, and 1800s alongside the specific state of charge (SoC) value for a specific constant discharge current value of ~1C until the relaxation stage to train and evaluate three regression-based machine learning models— Support Vector Regression (SVR), Extreme Gradient Boosting (XGBoost), and Gaussian Process Regression (GPR)—for estimating the ECM parameters of the selected model. Bayesian optimization is employed for hyperparameter tuning to achieve optimal performance for all the regressor models, among which, the GPR provided the best performance with the root-mean-squared error (RMSE) of less than 4x10-4 on average for the resistive components and less than 0.27 for capacitive components with excellent R2 scores. The simplicity of the approach enables it to eliminate the need for sophisticated measuring equipment and computation power.

Sagar, Md. Samiul [The University of Alabama (UA)]↗

DOE OSTI2021

Application of Dislocation Theory to Minimize Defects in Artificial Solids Built with Nanocrystal Building Blocks

Conspectus. Oriented atomic attachment of colloidal inorganic nanocrystals represents a powerful synthetic method for preparing complex inorganic superstructures. Examples include fusion of nanocrystals into dimer and superlattice structures. If the attachment were perfect throughout, then the resulting materials would have single crystal-like alignment of the individual nanocrystals' atomic lattices. While individual colloidal nanocrystals typically are free of many defects, there are a multitude of pathways that can generate defects upon nanocrystal attachment. These attachment generated defects are typically undesirable, and thus developing strategies to favor defect-free attachment or heal defective interfaces are essential. There may also be some cases where attachment-derived defects are desirable. In this Account, we summarize our current understanding of how these defects arise, in order to offer guidance to those who are designing nanocrystal derived solids. The small size of inorganic nanocrystals means short diffusion lengths to the surface, which favor the formation of nanocrystal building blocks with pristine atomic structures. Upon attachment, however, there are numerous pathways that can lead to atomic scale defects, and bulk crystal dislocation theory provides an invaluable guide to understanding these phenomena. As an example, an atomic step edge can be incorporated into the interface leading to an extra half-plane of atoms, known as an edge dislocation. These dislocations can be well described by the Burgers vector description of dislocations, which geometrically identifies planes in which a dislocation can move. Our in situ measurements have verified that bulk dislocation theory predictions for 1D defects hold true at few-nanometer length scales in PbTe and CdSe nanocrystal interfaces. Ultimately, the applicability of dislocation theory to nanocrystal attachment enables the predictive design of attachment to prevent or facilitate healing of defects upon nanocrystal attachment. We applied similar logic to understand formation of planar (2D) defects such as stacking faults upon nanocrystal attachment. Again concepts from bulk crystal defect crystallography can identify attachment pathways that can prevent or deterministically form planar defects upon nanocrystal attachment. The concepts we discuss work well for identifying favorable attachment geometries for nanocrystal pairs ; however it is currently unclear how to translate these ideas to near-simultaneous multiparticle attachment. Geometric frustration, which prevents nanocrystal rotation, and yet-to-be considered defect generation pathways unique to multiparticle attachment complicate defect-free superlattice attachment. New imaging methods now allow for the direct observation of local attachment trajectories and may enable improved understanding of such multiparticle phenomena. With further refinement, a unified framework for understanding and ultimately eliminating structural defects in fused nanocrystal superstructures may well be achievable in coming years.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2021

Powering the Blue Economy: A Survey of Station-Keeping Methods for Mooringless Platforms

The term “ocean platform” is used to reference everything from stationary, typically moored, buoys to mobile water vehicles, whether they operate on the ocean’s surface or underwater. For certain applications for which relatively stationary station-keeping conditions are desired, the use of mooring systems is not always a viable alternative either for economic, environmental, regulatory, or otherwise practical reasons, or a combination thereof, (e.g., short deployments, sensitive ecosystems, very deep project sites). Maintaining a platform at a single waypoint or reference location without being moored would require additional control systems and a power source to counteract the drift forces that would naturally displace it. Mobile platforms, which are usually untethered except for remotely operated vehicles, typically require energy input to power their station-keeping capabilities so that they hold or control their location in the ocean. Currently, most of these platforms use combustion engines or batteries for this purpose, which, depending on the specific systems, may be costly, pollute the environment, or create limitations on the length of the deployment. However, powering this kind of platforms with surrounding renewable resources (waves, currents, winds, or sun) has been identified as a promising solution to expand their application. The intent of this report is to investigate station-keeping methods for various ocean platforms that are not moored or otherwise anchored to the ocean floor, or another platform or vessel, paying particular interest to technologies that use marine renewable resources to power their operation, because that is of particular interest to the U.S. Department of Energy’s Powering the Blue Economy (PBE) initiative. As a first step, 72 articles and technical reports related to mooringless station-keeping methods were collected for review. The preliminary literature review provided a broad overview of common themes across the literature from which a descriptive methodology for analyzing various platforms was developed. That is, station-keeping methods were categorized based on their predominant energy source and consumption (renewable, nonrenewable, or hybrid if the platform uses renewable and nonrenewable resources equally), and their localization strategy (drift reduction, “path-planning or “waypoint-holding”). In addition, platform types were segregated into the following groups: buoys, surface drifters, and unoccupied surface vehicles (USVs); offshore renewable energy systems; and unoccupied underwater vehicles (UUVs). The main types of station-keeping methods encountered in this report achieve their intended localization strategy by means of drift mitigation, steering, and/or propulsion. Drift mitigation is commonly accomplished via drogues and sea anchors. Stand-along steering subsystems use control surfaces (e.g., ship rudder, wing sail, etc.) that react to ocean currents, waves, or winds to provide varying-degrees of course adjustments. Combined steering and propulsion subsystems include differential thrusters, directional thrusters separate from a primary thruster that cause the platform to pitch up/down or yaw clockwise/counterclockwise, or vectored thrusters that direct the propulsion in a range of directions relative to the platform’s local coordinate system. Propulsion is often achieved by running a motor and applying active control strategies but can also involve buoyancy shifts and using sails to generate lifting forces that propel a platform in a desired direction. Future research is primarily expected to take place in the form of a technoeconomic analysis that would aim to determine the technological viability, cost, and added value of mooringless station-keeping use cases identified through this research, including docking for UUV recharging or for georeferencing drifter buoys, deep-sea floating wind farms, U.S. Navy sonar arrays, and a Pacific Ocean wave buoy network.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

Exploring the parameter space of MagLIF implosions using similarity scaling. I. Theoretical framework

Magneto-inertial fusion concepts, such as the magnetized liner inertial fusion (MagLIF) platform, constitute an alternative path for achieving ignition and significant fusion yields in the laboratory. The space of experimental input parameters defining a MagLIF load is highly multi-dimensional, and the implosion itself is a complex event involving many physical processes. In the first paper of this series, we develop a simplified analytical model that identifies the main physical processes at play during a MagLIF implosion. Using non-dimensional analysis, we determine the most important dimensionless parameters characterizing MagLIF implosions and provide estimates of such parameters using typical fielded or experimentally observed quantities for MagLIF. Here, we then show that MagLIF loads can be “incompletely” similarity scaled, meaning that the experimental input parameters of MagLIF can be varied such that many (but not all) of the dimensionless quantities are conserved. Based on similarity-scaling arguments, we can explore the parameter space of MagLIF loads and estimate the performance of the scaled loads. Then, in the follow-up papers of this series, we test the similarity-scaling theory for MagLIF loads against simulations for two different scaling “vectors,” which include current scaling and rise-time scaling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Isovector Axial Charge and Form Factors of Nucleons from Lattice QCD

A survey of the calculations of the isovector axial vector form factor of the nucleon using lattice QCD is presented. Attention is paid to statistical and systematic uncertainties, in particular those due to excited state contributions. Based on a comparison of results from various collaborations, a case is made that lattice results are consistent within 10%. A similar level of uncertainty is in the axial charge gAu−d, the mean squared axial charge radius ⟨rA2⟩, the induced pseudoscalar charge gP∗, and the pion–nucleon coupling gπNN. Even with the current methodology, a significant reduction in errors is expected over the next few years with higher statistics data on more ensembles closer to the physical point. Lattice QCD results for the form factor GA(Q2) are compatible with those obtained from the recent MINERνA experiment but lie 2–3σ higher than the phenomenological extraction from the old ν–deuterium bubble chamber scattering data for Q2>0.3 GeV2. Current data show that the dipole ansatz does not have enough parameters to fit the form factor over the range 0≤Q2≤1 GeV2, whereas even a z2 truncation of the z expansion or a low order Padé are sufficient. Looking ahead, lattice QCD calculations will provide increasingly precise results over the range 0≤Q2≤1 GeV2, and MINERνA-like experiments will extend the range to Q2∼2 GeV2 or higher. Nevertheless, improvements in lattice methods to (i) further control excited state contributions and (ii) extend the range of Q2 are needed.

Astronomy & Astrophysics↗

DOE OSTI2024

A three-dimensional model of terrain-induced updrafts for movement ecology studies

Abstract Background Spatially explicit simulation models of animal movements through the atmosphere necessarily require a representation of the spatial and temporal variation of atmospheric conditions. In particular, for movements of soaring birds that rely extensively on vertical updrafts to avoid flapping flight, accurate and reliable estimation of the vertical component of wind is critical. The interaction between wind and complex terrain shapes both the horizontal and vertical wind fields, highlighting the need to model the coupling between local terrain features and atmospheric conditions at scales relevant to animal movement. Methods In this work, we propose a new empirical model for estimating the orographic updraft field. The model is developed using computational fluid dynamics simulations of canonical atmospheric conditions over moderately complex terrain. To isolate buoyancy and thermal effects, and focus on terrain-induced effects, we use only simulations of a neutrally stratified atmosphere to develop the model. The model, which we name Engineering Vertical Velocity Estimator (EVVE), is simple to implement and is a function of the underlying terrain elevation map, the desired height above ground level (AGL), and wind conditions at a reference height (80 m). We validate the model with data from the Alaiz mountain (Spain) field campaign. Results Compared to observations, the proposed improved model estimates the updrafts at 120 m AGL with a mean error of 0.11 m/s ( $$\sigma =0.28$$ σ = 0.28 m/s), compared to 0.85 m/s ( $$\sigma =0.58$$ σ = 0.58 m/s) for its baseline. For typical land-based wind turbine hub heights of 80 m AGL, the proposed model has a mean error of 0.04 m/s ( $$\sigma =0.25$$ σ = 0.25 m/s), compared to baseline 0.54 m/s ( $$\sigma =0.45$$ σ = 0.45 m/s) estimations. We illustrate an application of the model in movement ecology by comparing simulated tracks and presence maps of golden eagles ( Aquila chrysaetos ) moving across two distinct landscapes. The tracks and presence maps are obtained using a simple heuristic-based movement model, with the updraft field given by the proposed model and a wind vector-based estimation approach that is currently in wide use in movement ecology studies of raptors and other soaring birds. Conclusions We highlight that movement model results can be sensitive to the underlying orographic updraft model, especially in studies of fine-scale movements in regions of complex topography. We suggest adopting the proposed model rather than the wind vector estimation method for studies of soaring bird movements.

17 WIND ENERGY↗

DOE OSTI2026

Developing New Beam Constraints to Improve Sterile Neutrino Searches at NOvA

NOvA is a long-baseline neutrino oscillation experiment utilizing two functionally identical liquid scintillator detectors located 14 mrad off-axis from the NuMI beam at Fermilab. The position of the Near Detector (ND), 1 km from the neutrino production target, means it sees a range of off-axis angles (12.6-17.2 mrad), resulting in a noticeable energy gradient across the detector face. In addition, the neutrino production points vary widely along the decay pipe, which leads to a distribution of true neutrino baselines (L), and limits the sensitivity of our current sterile neutrino searches. This work presents a new "baseline reconstruction" technique designed to exploit our off-axis position to triangulate the precise origin of neutrinos within the NuMI beamline. By projecting the total momentum vector of Charged Current (CC) and Neutral Current (NC) interactions back to the central beam axis, we estimate the individual baseline for each event. We describe the technique, evaluate its performance, and demonstrate how this can improve the sensitivity of future searches for sterile neutrinos through reduced flux systematic uncertainties and improved L/E resolution.

Burns, Jessica [Cincinnati U.]↗

DOE OSTI2024

Nucleon Isovector Axial Form Factors

We present results for the isovector axial vector form factors obtained using thirteen 2 + 1 + 1-flavor highly improved staggered quark (HISQ) ensembles generated by the MILC collaboration. The calculation of nucleon two- and three-point correlation functions has been done using Wilson-clover fermions. In the analysis of these data, we quantify the sensitivity of the results to strategies used for removing excited state contamination and invoke the partially conserved axial current relation between the form factors to choose between them. Our data driven analysis includes removing contributions from multihadron $Nπ$ states that make significant contributions. Our final results are $g_A$ = 1.292(53) stat (24) sys for the axial charge; $g_S$ = 1.085(50) stat (103) sys and $g_T$ = 0.991(21) stat (10) sys for the scalar and tensor charges; $\langle{r^2_A}\rangle$ = 0.439(56) sta t(34) sys fm 2 for the mean squared axial charge radius, $g^*_P$ = 9.03(47) stat (42) sys for the induced pseudoscalar charge; and $g_{πNN}$ = 14.14(81) stat (85) sys for the pion-nucleon coupling. We also provide a parametrization of the axial form factor $G_A(Q^2)$ over the range 0 ≤ $Q^2$ ≤ 1 GeV 2 for use in phenomenology and a comparison with other lattice determinations. We find that the various lattice data agree within 10% but are significantly different from the extraction of $G_A(Q^2)$ from the $ν$-deuterium scattering data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Nucleon isovector axial form factors

We present results for the isovector axial vector form factors obtained using thirteen 2 + 1 + 1 -flavor highly improved staggered quark (HISQ) ensembles generated by the MILC collaboration. The calculation of nucleon two- and three-point correlation functions has been done using Wilson-clover fermions. In the analysis of these data, we quantify the sensitivity of the results to strategies used for removing excited state contamination and invoke the partially conserved axial current relation between the form factors to choose between them. Our data driven analysis includes removing contributions from multihadron N π states that make significant contributions. Our final results are g A = 1.292 ( 53 ) stat ( 24 ) sys for the axial charge; g S = 1.085 ( 50 ) stat ( 103 ) sys and g T = 0.991 ( 21 ) stat ( 10 ) sys for the scalar and tensor charges; ⟨ r A 2 ⟩ = 0.439 ( 56 ) stat ( 34 ) sys fm 2 for the mean squared axial charge radius, g P * = 9.03 ( 47 ) stat ( 42 ) sys for the induced pseudoscalar charge; and g π N N = 14.14 ( 81 ) stat ( 85 ) sys for the pion-nucleon coupling. We also provide a parametrization of the axial form factor G A ( Q 2 ) over the range 0 ≤ Q 2 ≤ 1 GeV 2 for use in phenomenology and a comparison with other lattice determinations. We find that the various lattice data agree within 10% but are significantly different from the extraction of G A ( Q 2 ) from the ν -deuterium scattering data. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Post-hazard condition assessment of nuclear piping-equipment systems: Novel approach to feature extraction and deep learning

Over the past decade, the use of artificial intelligence techniques in the field of health-monitoring has gained significant interest, especially for structures such as building and bridges. However, applications to industrial systems such as equipment-piping systems in nuclear plants have not been explored. Here, in this paper, it is shown that the existing techniques developed for buildings and bridges cannot be extended directly to equipment-piping systems as the response of such systems is governed by multiple localized modes unlike that in buildings and bridges. This paper proposes a new approach that consists of three key aspects: (i) a novel vector of degradation-sensitive features extracted from measured data, (ii) using a deep Artificial Neural Network (ANN) for diagnosis of degradation location and degradation severity, and (iii) consideration of uncertainty in degradation severity when training the ANN. Degradation in piping-equipment systems can occur due to flow-accelerated erosion and corrosion. These locations can potentially exhibit damage such as localized yielding or initiation of cracking due to an external event such as an earthquake. Moreover, such locations can at times go undetected by current inspection techniques. Therefore, a robust framework is needed for detection of degradation after a seismic event. This manuscript proposes a proof-of-concept framework, which utilizes data collected from sensors to generate a deep ANN database for predicting degraded locations and severity in a piping-equipment system. Degradation severity is classified as minor, moderate, and severe. In the suggested methodology, a novel vector of degradation-sensitive features is extracted from the sensor data to train the ANN. A simple piping-equipment system is selected to demonstrate feature extraction as a means to simplify pattern recognition, explore the design and parameters of an ANN, and develop a sensor placement strategy. The effectiveness of the proposed framework is demonstrated on a realistic primary safety system of a two-loop nuclear reactor. It is shown that the proposed post-hazard condition assessment framework is able to detect degraded locations along with the severity levels, including minor degradation, with considerably higher accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Machine Learning with Gradient-Based Optimization of Nuclear Waste Vitrification with Uncertainties and Constraints

Gekko is an optimization suite in Python that solves optimization problems involving mixed-integer, nonlinear, and differential equations. The purpose of this study is to integrate common Machine Learning (ML) algorithms such as Gaussian Process Regression (GPR), support vector regression (SVR), and artificial neural network (ANN) models into Gekko to solve data based optimization problems. Uncertainty quantification (UQ) is used alongside ML for better decision making. These methods include ensemble methods, model-specific methods, conformal predictions, and the delta method. An optimization problem involving nuclear waste vitrification is presented to demonstrate the benefit of ML in this field. ML models are compared against the current partial quadratic mixture (PQM) model in an optimization problem in Gekko. GPR with conformal uncertainty was chosen as the best substitute model as it had a lower mean squared error of 0.0025 compared to 0.018 and more confidently predicted a higher waste loading of 37.5 wt% compared to 34 wt%. The example problem shows that these tools can be used in similar industry settings where easier use and better performance is needed over classical approaches. Future works with these tools include expanding them with other regression models and UQ methods, and exploration into other optimization problems or dynamic control.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

Data-Driven Transition Path Analysis Yields a Statistical Understanding of Sudden Stratospheric Warming Events in an Idealized Model

Abstract Atmospheric regime transitions are highly impactful as drivers of extreme weather events, but pose two formidable modeling challenges: predicting the next event (weather forecasting) and characterizing the statistics of events of a given severity (the risk climatology). Each event has a different duration and spatial structure, making it hard to define an objective “average event.” We argue here that transition path theory (TPT), a stochastic process framework, is an appropriate tool for the task. We demonstrate TPT’s capacities on a wave–mean flow model of sudden stratospheric warmings (SSWs) developed by Holton and Mass, which is idealized enough for transparent TPT analysis but complex enough to demonstrate computational scalability. Whereas a recent article (Finkel et al. 2021) studied near-term SSW predictability, the present article uses TPT to link predictability to long-term SSW frequency. This requires not only forecasting forward in time from an initial condition, but also backward in time to assess the probability of the initial conditions themselves. TPT enables one to condition the dynamics on the regime transition occurring, and thus visualize its physical drivers with a vector field called the reactive current . The reactive current shows that before an SSW, dissipation and stochastic forcing drive a slow decay of vortex strength at lower altitudes. The response of upper-level winds is late and sudden, occurring only after the transition is almost complete from a probabilistic point of view. This case study demonstrates that TPT quantities, visualized in a space of physically meaningful variables, can help one understand the dynamics of regime transitions.

Meteorology & Atmospheric Sciences↗

DOE OSTI2024

Light-driven nanoscale vectorial currents

Controlled charge flows are fundamental to many areas of science and technology, serving as carriers of energy and information, as probes of material properties and dynamics and as a means of revealing or even inducing broken symmetries. Emerging methods for light-based current control offer particularly promising routes beyond the speed and adaptability limitations of conventional voltage-driven systems. However, optical generation and manipulation of currents at nanometre spatial scales remains a basic challenge and a crucial step towards scalable optoelectronic systems for microelectronics and information science. Here we introduce vectorial optoelectronic metasurfaces in which ultrafast light pulses induce local directional charge flows around symmetry-broken plasmonic nanostructures, with tunable responses and arbitrary patterning down to subdiffractive nanometre scales. Local symmetries and vectorial currents are revealed by polarization-dependent and wavelength-sensitive electrical readout and terahertz (THz) emission, whereas spatially tailored global currents are demonstrated in the direct generation of elusive broadband THz vector beams 17 . We show that, in graphene, a detailed interplay between electrodynamic, thermodynamic and hydrodynamic degrees of freedom gives rise to rapidly evolving nanoscale driving forces and charge flows under the extremely spatially and temporally localized excitation. These results set the stage for versatile patterning and optical control over nanoscale currents in materials diagnostics, THz spectroscopies, nanomagnetism and ultrafast information processing.

36 MATERIALS SCIENCE↗