Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Randomized methods”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Modeling plasticity-mediated void growth at the single crystal scale: A physics-informed machine learning approach

Modeling the evolution of voids during plastic flow as well as their effects on plastic dissipation is critical for both component manufacturing and lifetime estimation purposes. To this end, we propose a rate-dependent constitutive model to homogenize the effects of semi-randomly distributed voids on single crystal plasticity whilst capturing void interaction and plastic anisotropy. Here, this present work focuses on the case of face centered cubic crystals to introduce an anisotropic gauge function applicable within the crystal plasticity formalism. The approach combines analytical methods to describe the micromechanics of the system in combination with symbolic regression to capture analytically intractable mechanisms from data. The hybrid framework uses a physics-informed genetic programming-based symbolic regression algorithm to solve a multiform optimization problem simultaneously producing a new gauge function and a new strain rate equation. This is also a multi-objective optimization problem with many competing objectives. A new search and selection step is introduced to the genetic algorithm that promotes convergence toward a global solution that better satisfies all the objectives. Overall, the symbolic equations produced leverage data-driven methods to achieve greater accuracy than comparable alternatives on an analytically intractable problem while maintaining model transparency.

36 MATERIALS SCIENCE↗
DOE OSTI2025

CANA v1.0.0: efficient quantification of canalization in automata networks

The biomolecular networks underpinning cell function exhibit canalization, or the buffering of fluctuations required to function in a noisy environment. We present a new major release of $\tt{CANA}$, v1.0.0, an open-source Python package for understanding canalization in automata network models, discrete dynamical systems in which activation of biomolecular entities (e.g. transcription of genes) is modeled as the activity of coupled automata. One understudied putative mechanism for canalization is the functional equivalence of biomolecular regulators (e.g. among the transcription factors for a gene). We study this mechanism using the theory of symmetry in discrete functions. We present a new exact method, $\tt{schematodes}$, for finding maximal symmetry groups among the inputs to discrete functions, and integrate it into $\tt{CANA}$. The $\tt{schematodes}$ method substantially outperforms the inexact method of previous $\tt{CANA}$ versions both in speed and accuracy. We apply $\tt{CANA}$ v1.0.0 to study symmetry in 74 experimentally supported automata network models from the Cell Collective (CC) repository. The symmetry distribution is significantly different in the CC than in random automata with the same in-degree (connectivity) and bias (average output) (Kolmogorov–Smirnov test, P ≪ .001). Its spread is much wider than in a null model (IQR 0.31 versus IQR 0.20 with equal medians), demonstrating that the CC is enriched in functions with extreme symmetry or asymmetry.

Boolean networks↗
DOE OSTI2024

Landscape of nuclear deformation softness with spherical quasiparticle random-phase approximation

We investigate the stability and softness of nuclei against quadrupole, octupole, and hexadecapole deformation. By applying the spherical Skyrme-force Hartree-Fock Bardeen-Cooper-Schrieffer quasiparticle random phase approximation, we diagnose ground-state deformation when imaginary solutions are obtained, i.e., the spherical ground state collapses. We also calculate the multipole polarizability in spherical nuclei with no collapse, as a measure of softness. This numerically light and theoretically sound method is found able to capture deformation patterns across the nuclide chart. As a result, the connection between the intrinsic shape of nuclei and the dynamics of their low-lying collective states is established and the role of shell structure is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Efficient, Predictive Tomography of Multi-Qubit Quantum Processors

After decades of R&D, quantum computers comprising more than 2 qubits are appearing. If this progress is to continue, the research community requires a capability for precise characterization (“tomography”) of these enlarged devices, which will enable benchmarking, improvement, and finally certification as mission-ready. As world leaders in characterization -- our gate set tomography (GST) method is the current state of the art – the project team is keenly aware that every existing protocol is either (1) catastrophically inefficient for more than 2 qubits, or (2) not rich enough to predict device behavior. GST scales poorly, while the popular randomized benchmarking technique only measures a single aggregated error probability. This project explored a new insight: that the combinatorial explosion plaguing standard GST could be avoided by using an ansatz of few-qubit interactions to build a complete, efficient model for multi-qubit errors. We developed this approach, prototyped it, and tested it on a cutting-edge quantum processor developed by Rigetti Quantum Computing (RQC), a US-based startup. We implemented our new models within Sandia’s PyGSTi open-source code, and tested them experimentally on the RQC device by probing crosstalk. We found two major results: first, our schema worked and is viable for further development; second, while the Rigetti device is indeed a “real” 8-qubit quantum processor, its behavior fluctuated significantly over time while we were experimenting with it and this drift made it difficult to fit our models of crosstalk to the data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Predicting Initial Trans-Membrane Pressure for Optimized Operations in UF Unit Using Random Forest

With the growing scarcity of freshwater, innovative process design mechanisms like Reverse Osmosis (RO) are increasingly gaining attention among water treatment utilities to address the rising demand. Ensuring reliable water production necessitates efficient resource utilization, minimizing downtime in (ultra-filtration) UF systems. Recent advancements in machine learning (ML) have enabled the development of accurate data-driven models for Model Predictive Control (MPC), often requiring minimal prior knowledge of underlying physical processes. In this study, we present predictive regression models based on Random Forest (RF) and Auto-Regressive (AR) approaches to forecast the initial Trans-Membrane Pressure (TMP) for each filtration cycle in data generated by Direct Potable Reuse (DPR) systems. The proposed RF-based model demonstrates superior performance compared to baseline methods, including historical mean, Last Observation Carried Forward (LOCF), and naïve AR models, across various forecasting horizons in terms of root mean square error (RMSE) metric. To evaluate how different classes of process variables contribute to TMP dynamics over time, we examine the feature importance of independent covariates across multiple forecast horizons. This analysis provides insight into the temporal relevance of operational and sensor-derived features, guiding control and monitoring strategies. Additionally, the impact of hyperparameter tuning on TMP prediction performance is studied for both direct and recursive RF modelling approaches across increasing forecast horizons. Accurate prediction of initial TMP is critical for optimizing RO operations, as it enables the development of robust modelling frameworks by accurately estimating membrane fouling trends, thereby enhancing process efficiency and long-term reliability. The demonstrated efficacy of the RF-based approach highlights its potential as a tool for real-time decision-making in water treatment systems, paving the way for advanced process optimization and sustainable water resource management.

Mukherjee, Subrata [ORNL] (ORCID:0000000309930338)↗
DOE OSTI2025

Predicting Initial Trans-Membrane Pressure for Optimized Operations in UF Unit Using Random Forest

With the growing scarcity of freshwater, innovative process design mechanisms like Ultra-filtration(UF) units are increasingly gaining attention among water treatment utilities to address the rising demand. Ensuring reliable water production necessitates efficient resource utilization, minimizing downtime in UF systems. Recent advancements in machine learning (ML) have enabled the development of accurate data-driven models for Model Predictive Control (MPC), often requiring minimal prior knowledge of underlying physical processes. In this study, we present predictive regression models based on Random Forest (RF) and Auto-Regressive (AR) approaches to forecast the initial Trans-Membrane Pressure (TMP) for each filtration cycle in data generated by Direct Potable Reuse (DPR) systems. The proposed RF-based model demonstrates superior performance compared to baseline methods, including historical mean, Last Observation Carried Forward (LOCF), and naïve AR models, across various forecasting horizons in terms of root mean square (RMSE) metric. Accurate prediction of initial TMP is critical for optimizing CCRO operations, as it enables the development of robust modelling frameworks that enhance process efficiency and reliability. The demonstrated efficacy of the RF-based approach highlights its potential as a tool for real-time decision-making in water treatment systems, paving the way for advanced process optimization and sustainable water resource management.

Mukherjee, Subrata [ORNL] (ORCID:0000000309930338)↗
DOE OSTI2025

Green solvent-processable poly(3-hexylthiophene) derivative maintains high hole mobility in diodes

The development of environmentally sustainable fabrication methods for organic solar cells (OSCs) is critical to enable their large-scale adoption. Conventional solution-processed OSCs often rely on halogenated, toxic solvents that pose health and environmental risks, limiting their scalability. This is because π-conjugated polymers tend to have a low solubility in non-halogenated solvents. A common strategy to enhance solubility in alternative solvents is through the incorporation of polar solubilizing groups, often on every repeat unit of at least one monomer. However, too many polar side chains tend to disrupt favorable morphology and decrease charge transport properties. In this study, we demonstrate that a minimal degree of side chain functionalization can improve green solvent processability while preserving electrical performance. We tested this hypothesis using our previously reported random copolymer derivative of poly(3-hexylthiophene) (P3HT) where ∼10 mol% of the side chains are 6-pentanoatehexyl side chains. We find that this low-level functionalization significantly broadens the solvents that can be used to process P3HT to less toxic solvents such as o-xylene and anisole. Furthermore, the copolymer processed from greener solvents showed high hole mobilities (∼10⁻ 3 cm 2 /Vs by the Space-Charge Limiting Current method in diodes), comparable to P3HT films cast from toxic chloroform. These findings suggest that minimal side-chain modification is a viable strategy for expanding green solvent compatibility while preserving electrical performance, paving the way toward more sustainable organic semiconductor designs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ensemble variational Fokker-Planck methods for data assimilation

Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Itˆo process. This work introduces the Variational Fokker-Planck (VFP) framework for data assimilation, a general approach that includes previously known particle flow filters as special cases. The McKean-Vlasov-Itˆo process that transforms particles is defined via an optimal drift that depends on the selected diffusion term. It is established that the underlying probability density - sampled by the ensemble of particles - converges to the Bayesian posterior probability density. For a finite number of particles the optimal drift contains a regularization term that nudges particles toward becoming independent random variables. Based on this analysis, we derive computationally-feasible approximate regularization approaches that penalize the mutual information between pairs of particles, and avoid particle collapse. Moreover, the diffusion plays a role akin to a particle rejuvenation approach that aims to alleviate particle collapse. The VFP framework is very flexible. Different assumptions on prior and intermediate probability distributions can be used to implement the optimal drift, and localization and covariance shrinkage can be applied to alleviate the curse of dimensionality. A robust implicit-explicit method is discussed for the efficient integration of stiff McKean- Vlasov-Itˆo processes. Here, the effectiveness of the VFP framework is demonstrated on three progressively more challenging test problems, namely the Lorenz ’63, Lorenz ’96 and the quasi-geostrophic equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A Fast Algorithm for Scanning Transmission Electron Microscopy Imaging and 4D-STEM Diffraction Simulations

Scanning transmission electron microscopy (STEM) is an extremely versatile method for studying materials on the atomic scale. Many STEM experiments are supported or validated with electron scattering simulations. However, using the conventional multislice algorithm to perform these simulations can require extremely large calculation times, particularly for experiments with millions of probe positions as each probe position must be simulated independently. Recently, the plane-wave reciprocal-space interpolated scattering matrix (PRISM) algorithm was developed to reduce calculation times for large STEM simulations. Here, we introduce a new method for STEM simulation: partitioning of the STEM probe into “beamlets,” given by a natural neighbor interpolation of the parent beams. This idea is compatible with PRISM simulations and can lead to even larger improvements in simulation time, as well requiring significantly less computer random access memory (RAM). We have performed various simulations to demonstrate the advantages and disadvantages of partitioned PRISM STEM simulations. We find that this new algorithm is particularly useful for 4D-STEM simulations of large fields of view. We also provide a reference implementation of the multislice, PRISM, and partitioned PRISM algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Isoconversional Kinetic Analysis of Oxygen Desorption for Sr 0.75 Ca 0.25 FeO 3-δ Perovskite in Dry and Steam-Based Environments

The desorption kinetics of Sr 0.75 Ca 0.25 FeO 3-δ perovskite material are examined in both dry and steam-based environments using redox gaseous products of a laboratory-scale fixed bed. First, desorption kinetics associated with the dry environment are proposed based on Friedman’s isoconversional method. It is found from the reconstructed reaction model that the desorption kinetics of this perovskite are controlled by a three-step mechanism, where hypothetically, the phase-boundary reaction rapidly prevails first, followed by diffusion and/or nucleation, and finally, unimolecular decay or random nucleation-controlled reaction. Subsequent analyses of the Arrhenius parameters indicated that both the early and later desorption stages, respectively, controlled by the phase boundary and unimolecular/random nucleation reactions, are associated with low energetic demand. In contrast, a high energetic penalty is required for the diffusion/nucleation reaction. Further, a satisfactory a priori verification of the proposed kinetics suggested that the three-step reaction mechanism reasonably describes the desorption of this perovskite in a dry environment. With the presence of steam in the desorption environment, oxygen production is substantially inhibited, but the hypothetical three-step mechanism is still found to control the desorption. Meanwhile, the primary effects of the steam on these mechanisms include the widening of the conversion extent associated with the diffusion/nucleation reaction, along with a higher energetic demand compared to the dry environment.

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

Machine learning-enabled discovery of ionic liquid–solvent electrolytes exhibiting high ionic conductivity

Ionic liquids (ILs), which are a class of materials with versatile nature and growing popularity, are facing impediments toward widespread usage as electrolytes due to various factors such as low ionic conductivity, high viscosity, high market price etc. One of the ways these limitations can be addressed is by mixing ILs with a molecular solvent. In a combinatorial sense, there exists an immense number of specific IL–solvent combinations. An exhaustive experimental or even simulation-based investigation of the chemical space spanned by such combinations can be extremely time-consuming, expensive, and nearly impossible. An alternative approach is to employ machine learning-based models developed from available databases. Although there exists prior literature that integrates machine learning to investigate mixtures of specific solvents with ILs, these models lack generalization necessitating development of a large number of ML models to handle various solvents. To remedy this shortcoming, as a part of designing green electrolytes with high ionic conductivity that can have potential applications in next-generation batteries and solar cells, this work aims to develop a unified machine learning model to predict ionic conductivity of any IL–solvent mixture system. In this regard, three models, namely, Random Forest, extreme gradient boosting (XGBoost), and artificial neural network (ANN) were formulated using the NIST ILThermo database. The dataset contained 549 unique ionic liquids from 16 cation families and 81 unique solvents, representing a total of 23 712 datapoints. SHAPLEY additive explanation (SHAP) method was used to assess the impact of various features on model prediction and their significance was compared with literature to gain physical insight about the model behavior. Finally, using the developed models, approximately 2.5 million IL–solvent mixtures at five different compositions were screened at room temperature. The high-throughput screening yielded nearly 19 000 IL–solvent mixtures for which ionic conductivity was found to exceed the ionic conductivity of conventional Li-ion battery electrolyte.

25 ENERGY STORAGE↗
DOE OSTI2024

A Novel Authentication Management for the Data Security of Smart Grid

Bidirectional wireless communication is employed in various smart grid components such as smart meters and control and monitoring applications where security is vital. The Trusted Third Party (TTP) and wireless connectivity between the smart meter and the third party in the key management-based encryption techniques for the smart grid are expected to be totally trustworthy and dependable. In a wired/wireless medium, however, a man-in-the-middle may seek to disrupt, monitor and manipulate the network, or simply execute a replay attack, revealing its vulnerability. Recognizing this, this study presents a novel authentication management (model) comprised of two layer security schema. The first layer implements an efficient novel encryption method for secure data exchange between meters and control center with the help of two partially trusted simple servers (constitutes the TTP). In this setting, one server handles the data encryption between the meter and control center/central database, and the other server administers the random sequence of data transmission. The second layer monitors and verifies exchanged data packets among smart meters. It detects abnormal packets from suspicious sources. To implement this node-to-node authentication, One class support vector machine algorithm is proposed which takes advantages of the location information as well as the data transmission history (node identification, packet size, and data transmission frequency). This schema secures data communication, and imposes a comprehensive privacy throughout the system without considerably extending the complexity of the conventional key management scheme.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Observing Profiles of Derived Kinematic Field Quantities Using a Network of Profiling Sites

Abstract Observations of thermodynamic and kinematic parameters associated with derivatives of the thermodynamics and wind fields, namely, advection, vorticity, divergence, and deformation, can be obtained by applying Green’s theorem to a network of observing sites. The five nodes that comprise the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) profiling network, spaced 50–80 km apart, are used to obtain measurements of these parameters over a finite region. To demonstrate the applicability of this technique at this location, it is first applied to gridded model output from the High-Resolution Rapid Refresh (HRRR) numerical weather prediction model, using profiles from the locations of ARM network sites, so that values calculated from this method can be directly compared to finite difference calculations. Good agreement is found between both approaches as well as between the model and values calculated from the observations. Uncertainties for the observations are obtained via a Monte Carlo process in which the profiles are randomly perturbed in accordance with their known error characteristics. The existing size of the ARM network is well suited to capturing these parameters, with strong correlations to model values and smaller uncertainties than a more closely spaced network, yet it is small enough that it avoids the tendency for advection to go to zero over a large area.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Localization and coherent imaging of hidden moving objects using laser speckle

Imaging and sensing of moving objects through opaque scattering media is a challenging but important problem in a variety of applications, including environmental sensing, biomedical imaging, and material inspection. We have previously demonstrated a technique to coherently image a moving object through thick, heavily scattering random media using correlations of speckle images as a function of the object’s spatial translation. Here, we demonstrate that this technique can be combined with localization to achieve imaging without prior knowledge of the object’s motion, greatly extending the application domain. This method is effective beyond the thin or weakly scattering regime and, rather than motion being deleterious, exploits the information available when the hidden object is moving, as could be the case in a cluttered terrestrial environment or through substantial levels of biological tissue scatter.

Hastings, Ryan L. (ORCID:0009000095977807)↗
DOE OSTI2024

Compiling Quantum Circuits for Dynamically Field-Programmable Neutral Atoms Array Processors

Dynamically field-programmable qubit arrays (DPQA) have recently emerged as a promising platform for quantum information processing. In DPQA, atomic qubits are selectively loaded into arrays of optical traps that can be reconfigured during the computation itself. Leveraging qubit transport and parallel, entangling quantum operations, different pairs of qubits, even those initially far away, can be entangled at different stages of the quantum program execution. Such reconfigurability and non-local connectivity present new challenges for compilation, especially in the layout synthesis step which places and routes the qubits and schedules the gates. In this paper, we consider a DPQA architecture that contains multiple arrays and supports 2D array movements, representing cutting-edge experimental platforms. Within this architecture, we discretize the state space and formulate layout synthesis as a satisfiability modulo theories problem, which can be solved by existing solvers optimally in terms of circuit depth. For a set of benchmark circuits generated by random graphs with complex connectivities, our compiler OLSQ-DPQA reduces the number of two-qubit entangling gates on small problem instances by 1.7x compared to optimal compilation results on a fixed planar architecture. To further improve scalability and practicality of the method, we introduce a greedy heuristic inspired by the iterative peeling approach in classical integrated circuit routing. Using a hybrid approach that combined the greedy and optimal methods, we demonstrate that our DPQA-based compiled circuits feature reduced scaling overhead compared to a grid fixed architecture, resulting in 5.1X less two-qubit gates for 90 qubit quantum circuits. These methods enable programmable, complex quantum circuits with neutral atom quantum computers, as well as informing both future compilers and future hardware choices.

Physics↗
DOE OSTI2022

Comparing the extraction performance of cyclodextrin-containing supramolecular deep eutectic solvents versus conventional deep eutectic solvents by headspace single drop microextraction

A headspace single drop microextraction (HS-SDME) method coupled with high performance liquid chromatography was developed to compare the extraction of eighteen aromatic organic pollutants from aqueous solutions using cyclodextrin-based supramolecular deep eutectic solvents (SUPRADESs) and alkylammonium halide-based conventional deep eutectic solvents (DESs). Different derivatives of beta-cyclodextrin (β-CD) were employed as hydrogen bond acceptors (HBA) in SUPRADESs and the extraction performance investigated. SUPRADES comprised of the 20 wt% native β-CD HBA provided the highest enrichment factors of analytes compared to SUPRADESs comprised of other derivatives of β-CD (random methylated β-cyclodextrin, heptakis(2,3,6-tri-O-methyl)-β-cyclodextrin, and 2-hydroxypropyl β-cyclodextrin). In addition, native β-CD and its derivatives were dissolved in the neat DESs and their effect on the extraction of analytes examined. Dissolution of 20 wt% native β-CD in the choline chloride ([Ch + ][Cl - ]):2Urea DES resulted in a significant increase in the extraction efficiencies of target analytes compared to the neat [Ch + ][Cl - ]:2Urea DES. Under optimum conditions, the extraction method required a solvent microdroplet of 6.5 μL, 1000 rpm stir rate, 30% (w/v) salt concentration, and a temperature of 40 °C. The tetrabutylammonium chloride: 2 lactic acid DES resulted in the highest enrichment factors while the [Ch + ][Cl - ]:2Urea DES had the lowest for most of the analytes among the evaluated solvents. The method provided limits of detection (LODs) down to 35 μg L -1 . Finally, the developed method was applied for the analysis of spiked tap and lake water, where relative recoveries ranging from 83.7% -119.7% and relative standard deviations lower than 19.2% were achieved.

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

Sub‐Diffraction Correlation of Quantum Emitters and Local Strain Fields in Strain‐Engineered WSe 2 Monolayers

Strain-engineering in atomically thin metal dichalcogenides is a useful method for realizing single-photon emitters (SPEs) for quantum technologies. Correlating SPE position with local strain topography is challenging due to localization inaccuracies from the diffraction limit. Currently, SPEs are assumed to be positioned at the highest strained location and are typically identified by randomly screening narrow-linewidth emitters, of which only a few are spectrally pure. In this work, hyperspectral quantum emitter localization microscopy is used to locate 33 SPEs in nanoparticle-strained WSe 2 monolayers with sub-diffraction-limit resolution (≈30 nm) and correlate their positions with the underlying strain field via image registration. In this system, spectrally pure emitters are not concentrated at the highest strain location due to spectral contamination; instead, isolable SPEs are distributed away from points of peak strain with an average displacement of 240 nm. These observations point toward a need for a change in the design rules for strain-engineered SPEs and constitute a key step toward realizing next-generation quantum optical architectures.

2D material↗
DOE OSTI2024

Engineering phonon transport through cation disorder in dimensionally constricted high entropy MXene

Designing materials with low thermal conductivity is a crucial objective for applications in thermal insulation and thermoelectrics. Traditional methods such as doping, mechanical strain and introducing defects in perfect crystals have been widely explored to impede the flow of heat. Here, this work introduces dimensional constriction and cationic disorder as novel avenues to manipulate lattice thermal conductivity (LTC). High entropy materials characterized by random distribution of multiple elements, creates a suitable environment for thermal insulation due to its configurational disorder and local lattice distortions. On the other hand, MXenes, derived from MAX-phase, have garnered considerable attention due to their unique structural attributes, leading to potential applications in catalysis and energy storage. Ti 2 AlC MAX-phase is examined to understand the impact of dimensional constriction on phonon transport of Ti 2 C with cationic disorder, i.e., (Ti 0.25 Nb 0.25 Cr 0.25 Ta 0.25 ) 2 C. The exponential reduction in LTC of HE-MXene is attributed to disorder scattering that significantly limits phonon mean free path (MFP) and relaxation time. The spread of mode-resolved LTC with MFP highlights the influence of disorder on phonon scattering. This work provides a systematic approach to engineer LTC through dimensional constriction and cationic disorder, laying the foundation for tailored materials with desired thermal properties.

2D materials↗