Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Functional priors”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Experimental and Numerical Investigation of Flow Boiling in Additive Manufactured Foam Structures With Vapor Pathways

Abstract The unique properties of metal foams make them potential candidates for a range of applications, including microsystem thermal management. Using additive manufacturing to create foam-type structures can improve upon prior thermal solutions by eliminating thermal interface materials and allowing for customization/local control of parameters. In the present investigation, flow boiling in additive-manufactured metal foams is investigated both experimentally and numerically. Two test samples, one with uniform structure and the other with pathways for vapor removal, are compared both experimentally and numerically. A conjugate computational fluid dynamics and heat transfer (CFD-HT) model utilizing a three-dimensional volume of fluid (VOF) model with accompanying evaporation/condensation model provided in-depth visualization of the boiling flow phenomena. The experiments generated the thermohydraulic performance over a range of heat fluxes, demonstrating that the sample incorporating dedicated vapor pathways performed better in both pressure and heat transfer performance metrics compared to the uniform foam. Additionally, negative system-level effects (i.e., hydraulic oscillations) were shown to be abated using the vapor removal structures. The numerical model yielded further insight into the factors contributing to the improved performance. Results indicated the pathways functioned as vapor removal channels, allowing the generated vapor to vent from the foam structure into the lanes. Further computational investigations demonstrated changes in flow regimes, where the addition of vapor channels caused the flow to change from churn to annular. Bubble behavior unique to the vapor pathway structure was studied, showing stagnant regions that eject vapor into the channel.

Engineering↗
DOE OSTI2023

Trimming and Decontamination of Metagenomic Data can Significantly Impact Assembly and Binning Metrics, Phylogenomic and Functional Analysis

Background: Investigators using metagenomic sequencing to study microbiomes often trim and decontaminate reads without knowing their effect on downstream analyses. Objective: This study was designed to evaluate the impacts JGI trimming and decontamination procedures have on assembly and binning metrics, placement of MAGs into species trees, and functional profiles of MAGs extracted from complex rhizosphere metagenomes, as well as how more aggressive trimming impacts these binning metrics. Methods: Twenty-three Miscanthus x giganteus rhizosphere metagenomes were subjected to different combinations and thresholds of force, kmer, and quality trimming and decontamination using BBDuk. Reads were assembled and binned in KBase. Phylogenomic and statistical analyses were applied to evaluate the effects of trimming and decontamination on downstream analyses. Results: We found that JGI trimmed and decontaminated reads had significant impacts on assembly and binning metrics compared to raw reads, including significantly higher total contig counts, more contigs greater than 10k bp in length, and larger total lengths of raw assemblies compared to QC assemblies, and 2.0% lower average contamination of QC MAGs compared to raw MAGs. We also found that differences in the placement of MAGs in species trees increased with decreasing completeness and contamination thresholds. Furthermore, aggressive trimming (Q20) was found to significantly reduce MAG counts. Conclusion: Trimming and decontamination of metagenomics reads prior to assembly can change an investigator’s answer to the questions, “Who is there and what are they doing?” However, mild trimming and decontamination of metagenomic reads with high-quality scores are recommended for removing sample processing and sequencing artifacts.

Whitham, Jason M.↗
DOE OSTI2020

Tunable Perovskite-Derived Bismuth Halides: Cs 3 Bi 2 (Cl 1– x I x ) 9

Bismuth-based perovskites are of interest as safer alternatives to lead-based optoelectronic materials. Prior studies have reported on the compounds Cs 3 Bi 2 Cl 9 , Cs 3 Bi 2 I 9 , and Cs 3 Bi 2 Cl 3 I 6 . Here we examine a range of compounds of the formula Cs 3 Bi 2 (Cl 1–x I x ) 9 , where x takes values from 0.09 to 0.52. Powder and single-crystal X-ray diffraction were used to determine that all of these compounds adopt the layered vacancy-ordered perovskite structure observed for Cs 3 Bi 2 Cl 3 I 6 , which is also the high-temperature phase of Cs 3 Bi 2 Cl 9 . We find that, even with very small iodine incorporation, the structure is switched to that of Cs 3 Bi 2 Cl 3 I 6 , with I atoms displaying a distinct preference for the capping sites on the BiX 6 octahedra. Optical absorption spectroscopy was employed to study the evolution of optical properties of these materials, and this is complemented by density functional theory electronic structure calculations. Finally, three main absorption features were observed for these compounds, and with increasing x, the lowest-energy features are red-shifted.

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

Genetic Basis of Chromate Adaptation and the Role of the Pre-existing Genetic Divergence during an Experimental Evolution Study with Desulfovibrio vulgaris Populations

Hexavalent chromium [Cr(VI)] is a common environmental pollutant. However, little is known about the genetic basis of microbial evolution under Cr(VI) stress and the influence of the prior evolution histories on the subsequent evolution under Cr(VI) stress. In this study, Desulfovibrio vulgaris Hildenborough (DvH), a model sulfate-reducing bacterium, was experimentally evolved for 600 generations. By evolving the replicate populations of three genetically diverse DvH clones, including ancestor (AN, without prior experimental evolution history), non-stress-evolved EC3-10, and salt stress-evolved ES9-11, the contributions of adaptation, chance, and pre-existing genetic divergence to the evolution under Cr(VI) stress were able to be dissected. Significantly decreased lag phases under Cr(VI) stress were observed in most evolved populations, while increased Cr(VI) reduction rates were primarily observed in populations evolved from EC3-10 and ES9-11. The pre-existing genetic divergence in the starting clones showed strong influences on the changes in lag phases, growth rates, and Cr(VI) reduction rates. Additionally, the genomic mutation spectra in populations evolved from different starting clones were significantly different. A total of 14 newly mutated genes obtained mutations in at least two evolved populations, suggesting their importance in Cr(VI) adaptation. An in-frame deletion mutation of one of these genes, the chromate transporter gene DVU0426, demonstrated that it played an important role in Cr(VI) tolerance. Overall, our study identified potential key functional genes for Cr(VI) tolerance and demonstrated the important role of pre-existing genetic divergence in evolution under Cr(VI) stress conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Multilayer Network Clarifies Prevailing Water Consumption Telecouplings in the United States

Virtual water flows are used to map the indirect water consumption connections implied by the supply chain of a city, region, or country. This information can be used to manage supply chains to achieve environmental policy objectives and mitigate environmental risks to critical supply chains. A limitation of prior work is that these flows are typically analyzed using monolayer networks, which ignores crucial intersectoral or interlayer couplings. Here, in this study, we use a multilayer network to account for such couplings when analyzing blue virtual water flows in the United States. Our multilayer network consists of 115 different regions (nodes), covering the entire conterminous United States; 41 coupled economic sectors (layers); and ∼2 × 10 7 possible links. To analyze the multilayer network, we focus on three fundamental network properties: topological connectivity, mesoscale structure, and node centrality. The network has a high connectivity, with each node being on average connected to roughly 2/3 of the network's nodes. Interlayer flows are a major driver of connectivity, representing ∼54% of all the network's connections. Five different groups of tightly connected nodes (communities) characterize the network. Each community represents a preferred spatial mode of long-range virtual water interaction within the United States. We find that large (populous) cities have a stronger influence than small ones on network functioning because they attract and recirculate more virtual water through their supply chains. Our results also highlight differences between the multilayer and monolayer virtual water network, which overall show that the former provides a more realistic representation of virtual water flows.

Garcia Gonzalez, Susana↗
DOE OSTI2021

Cesium Ion Exchange Testing Using Crystalline Silicotitanate with Hanford Tank Waste 241-AP-107

At the time of this testing, the Low-Activity Waste Pretreatment System (LAWPS) was to provide for the initial production of immobilized low-activity waste by feeding Hanford tank supernate from tank farms to the Hanford Tank Waste Treatment and Immobilization Plant (WTP) Low-Activity Waste (LAW) Facility for immobilization. Washington River Protection Solutions requested that Hanford tank waste collected from tank 241-AP-107 (hereafter called AP-107) be processed using conceived pretreatment steps (suspended solids removal by filtration, Cs removal by ion exchange) then vitrified. A small-scale test platform to demonstrate the solids filtration, Cs removal, and LAW vitrification was constructed and installed at Pacific Northwest National Laboratory. Bench-scale ion exchange testing with approximately 9 L of AP-107 supernate was conducted using crystalline silicotitanate (CST) ion exchange media. The IONSIV R9140-B CST was provided by Honeywell UOP, LLC in 2018 (Batch 2081000057). The ion exchange media was first tested with simulant and was previously described. This report describes the Cs ion exchange batch contact and column test results with the AP-107 tank waste. Batch contact testing helps to evaluate CST performance on tank waste supernate prior to processing it in the ion exchange columns. Batch contacts were performed with the waste at four Cs concentrations at a phase ratio of 200 (liquid volume to exchanger mass) with AP-107. The distribution coefficient (K d ) at the equilibrium condition of 8.57 µg Cs/mL (AP-107 feed condition) was determined to be 669 mL AP-107/g CST. With a CST bed density of 1.00 g/mL, this K d corresponded to a predicted 50% Cs breakthrough of 669 bed volumes (BVs). The Cs load capacity at the equilibrium feed condition was determined to be 7.5 mg Cs/g dry CST. The column testing was prototypic to the intended LAWPS operations in a lead-lag column format, although on a small-scale basis with 10-mL CST beds. The feed was processed downflow through the lead column and then through the lag column at ~2.2 BV/h. Loading continued until the lag column reached the WTP waste acceptance criteria (WAC) for receiving supernatant waste for vitrification (a function of the Na and 137 Cs concentrations). For AP-107, the WAC is 0.114% of the influent 137 Cs concentration; this required a Cs decontamination factor of 876. The Cs effluent from the lag column reached the WAC after processing ~410 BVs. To keep the subsequent product effluent below the WAC, a replacement lag column was prepared, the lead column was removed from service (after processing a total of 471 BVs), the lag column was put into the lead column position, and the replacement lag column was installed. Feed processing continued and after another ~290 BVs the Cs effluent from the lag column again exceeded the WAC. In both cases, the lead columns only reached 25% Cs breakthrough before removal. Although 50% Cs breakthrough was not reached, this value was estimated and averaged based on extrapolation of the loading curves (640 BVs) and agreed within 4% of the predicted 50% Cs breakthrough from batch contact test results (669 BVs). Table ES.1 summarizes the observed column performance and relevant Cs loading characteristics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Non-intrusive reduced order modeling of natural convection in porous media using convolutional autoencoders: Comparison with linear subspace techniques

Natural convection in porous media is a highly nonlinear multiphysical problem relevant to many engineering applications (e.g., the process of CO 2 sequestration). Here, we extend and present a non-intrusive reduced order model of natural convection in porous media employing deep convolutional autoencoders for the compression and reconstruction and either radial basis function (RBF) interpolation or artificial neural networks (ANNs) for mapping parameters of partial differential equations (PDEs) on the corresponding nonlinear manifolds. To benchmark our approach, we also describe linear compression and reconstruction processes relying on proper orthogonal decomposition (POD) and ANNs. Further, we present comprehensive comparisons among different models through three benchmark problems. The reduced order models, linear and nonlinear approaches, are much faster than the finite element model, obtaining a maximum speed-up of 7 × 10 6 because our framework is not bound by the Courant–Friedrichs–Lewy condition; hence, it could deliver quantities of interest at any given time contrary to the finite element model. Our model’s accuracy still lies within a relative error of 7% in the worst-case scenario. We illustrate that, in specific settings, the nonlinear approach outperforms its linear counterpart and vice versa. We hypothesize that a visual comparison between principal component analysis (PCA) and t-Distributed Stochastic Neighbor Embedding (t-SNE) could indicate which method will perform better prior to employing any specific compression strategy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Flexoelectricity and New Phenomena

It is easy to miss the scientific implications of our recent work on Triboelectricity. Everyone knows that rubbing and contact can produce static electricity; less appreciated is that the thermodynamic driver has been an open question since static electricity was first observed by Thales of Miletus around 585 BC. People think they understand it, for instance one common explanation that can still be found in the current literature is that differences in the work function drives charge transfer, often called the Volta-Helmholtz hypothesis. As summarized in 1967 by Harper, this fails to explain many experimental observations, for instance that charging can occur when two pieces of the same material are rubbed against each other. We are the first to place triboelectricity on a solid foundation rooted in quantum mechanics – the flexoelectric effect. We were able to explain a significant number of previously unexplained phenomena: 1) Bipolar tribocurrents associated with stick-slip, due to the change in sign of the strain gradients. 2) A one-third power scaling of tribocurrents with indentation force. 3) Tribocharging when two identical materials are used, these being due to local variations in the asperities so there are usually local potential differences. 4) Inhomogeneous charging of insulators, related to the statistical nature of asperities. 5) An experimentally observed reversal in the sign of charge transfer for negative and positive curvature, which is related to a change in the sign of the strain gradient. Exploiting our DOE prior funded work on flexoelectricity, we obtained semi-quantitative matching to existing experimental measurements of the surface charge in triboelectric experiments. The work has been well received in the literature. The work has been the focus of a number of popular science press articles, and also formed the basis for a Podcast for children 6-10 “The Rise and Fall of Static Man” posted in December 2019 by NPR as part of their “Wow in the World” series. I was also briefly interviewed by the Chicago PBS station in January 2020. This work are significant for a wide range of energy applications; to quote from an independent source: "Triboelectric power has plenty of potential, says Wenzhuo Wu, an assistant professor of engineering at Purdue. If the basics of static electricity are better understood, we could maximize the efficiency of wind or wave power generators, Wu says. The body's own movement could be used to power internal medical devices. Imagine being able to create a roof shaped to harness the power of a raindrop — the friction of the rain passing over the surface — to generate triboelectricity, powering the building below it." This is the start of new science, some of which we already partially understand such as the role of band bending in charge transfer. We need to understand charge transfer combining elasticity, quantum mechanics, band bending and defect states. These directly involve several of the DOE Grand Challenges: "How do we control material processes at the level of electrons? How do remarkable properties of matter emerge from complex correlations of the atomic or electronic constituents and how can we control these properties? How do we characterize and control matter away—especially very far away— from equilibrium?" I will argue that this work truly falls into the class of disruptive science; it is not just a simple extension, linear science. Not everyone will accept the approach. Since we explain far more about triboelectricity than anyone before, the preponderance of evidence supports the model. The feedback I have received is that many agree with the work, to quote: "The model makes sense, says Michael McAlpine, a professor of engineering at the University of Minnesota. "It's such a simple explanation, I was surprised I didn't put my finger on that," McAlpine says." The proposal received strong reviews. It was also publicized on the Department of Energy Web Site.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Reactions of hydrazine with the amidogen radical and atomic hydrogen

The rate coefficient k 1 for NH 2 + N 2 H 4 was measured to be (5.4 ± 0.4) × 10 –14 cm 3 molecule –1 s –1 at 296 K. NH 2 was generated by pulsed laser photolysis of NH 3 at 193 nm, and monitored as a function of time by pulsed laser-induced fluorescence excited at 570.3 nm under pseudo-first order conditions in the presence of excess N 2 H 4 in an Ar bath gas. This reaction was also investigated computationally, with geometries and scaled frequencies obtained with M06-2X/6-311+G(2df,2p) theory, and single-point energies from CCSD(T)-F12b/cc-pVTZ-F12 theory, plus a term to correct approximately for electron correlation through CCSDT(Q). Three connected transition states are involved and rate constants were obtained via Multistructural Improved Canonical Variational Transition State Theory with Small Curvature Tunneling. Combination of experiment and theory leads to a recommended rate coefficient for hydrogen abstraction of k 1 = 6.3 × 10 –23 T 3.44 exp(+289 K/T) cm 3 molecule –1 s –1 . The minor channel for H + N 2 H 4 forming NH 2 + NH 3 was characterized computationally as well, to yield 5.0 × 10 –19 T 2.07 exp(-4032 K/T) cm 3 molecule –1 s –1 . Finally, these results are compared to several discordant prior estimates, and are employed in an overall mechanism to compare with measurements of half-lives of hydrazine in a shock tube.

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

Phospholipase D and phosphatidylinositol-4-phosphate 5-kinase 1 are involved in the regulation of oligodendrocyte morphological differentiation

Highlights: • PLD and PIP5K1 are involved in the regulation of oligodendroglial cell morphological differentiation. • Inhibition of PLD or PIP5K1 decreases differentiation or myelin marker expression. • PLD and PIP5K1 regulate differentiation of oligodendrocyte precursor cells. • Inhibition of PLD or PIP5K1 decreases differentiation marker expression of oligodendrocyte precursor cells. Oligodendroglial cells (oligodendrocytes) differentiate to form the myelin that wraps neuronal axons in the central nervous system (CNS). This myelin sheath supports the propagation of saltatory conduction and protects axons from physical stresses. When oligodendrocytes do not normally differentiate to myelinate axons, their key functions as oligodendrocytes in the CNS are severely impaired. The molecular mechanics that control differentiation still remain to be clarified. Arf6 belongs to the small GTPase family and is known to be a positive regulator of oligodendrocyte differentiation. Here, we show that the phospholipase D (PLD) and phosphatidylinositol-4-phosphate 5-kinase 1 (PIP5K1) molecules, the major effectors of Arf6, are involved in the regulation of oligodendrocyte differentiation. Knockdown of PLD1 or PIP5K type 1γ (PIP5K1C) by their respective specific siRNAs in mouse oligodendroglial FBD-102b cells inhibited morphological differentiation into structures bearing myelin-like processes; this finding is consistent with the concurrent changes in expression of differentiation and myelin marker proteins. Treatment with VU0155069 or UNC3230, specific inhibitors of PLD and PIP5K1, respectively, blunted morphological differentiation and decreased expression of myelin and differentiation marker proteins. Similar results have been obtained in studies using primary oligodendrocytes. These results suggest that the major Arf6 effector molecules PLD and PIP5K1 are among the molecules involved in the regulation of morphological differentiation in oligodendrocytes prior to myelination.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Formulation and solution approach for calibrating activity-based travel demand model-system via microsimulation

This study addresses the problem of calibrating utility-maximizing nested logit activity-based travel demand model-systems. After estimation, it is common practice to use aggregate measurements to calibrate the estimated model-system’s parameters prior to their application in transportation planning, policy making, and operations. However, calibration of activity-based model-systems has received much less attention. Existing calibration approaches are myopic heuristics in the sense that they do not consider the fundamental inter-dependencies among choice-models and do not have a systematic way to adjust model parameters. Also, other purely simulation-based approaches do not perform well in large-scale applications. In this study, we focus on utility-maximizing nested logit activity-based model-systems and calibrating aggregate statistics such as activity shares, mode shares, time-dependent & mode-specific OD flows, and time-dependent & mode-specific sensor counts. We formulate the calibration problem as a simulation-based optimization problem and propose a stochastic gradient-based solution procedure to solve it. The solution procedure relies on microsimulation to calculate expectations of the aggregate statistics of interest to the calibration problem. Additionally, we derive approximate analytical expressions for the gradient of the objective function —that are evaluated through microsimulation on mini-batches of the population. The proposed solution procedure is sensitive to the fundamental structure of the activity-based model-system and is non-myopic in considering the dependencies across its model components. The formulated optimization problem is non-convex, highly nonlinear, and potentially has multiple-minima. Lastly, we show —through a real-world application— that the proposed solution procedure outperforms other state-of-the-art purely simulation-based optimization approaches in terms of computational efficiency, stability, and convergence. We also compare various gradient-based solution algorithms to determine the best algorithm to update the parameters. This work has the potential to facilitate wider and easier application of activity-based model-systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Exploring the Future Energy Value of Long-Duration Energy Storage

Long-duration energy storage is commonly viewed as a key technology for providing flexibility to the grid and broader energy systems over a multidecadal time frame. However, prior work has typically used present-day grid infrastructures to characterize the relationship between the duration and arbitrage value of storage in electricity markets. This study leverages established National Renewable Energy Laboratory grid planning and operations tools, analysis, and data to execute a price-taker model of an energy storage system for several 8760 h price series representative of current and future contiguous United States grid infrastructures with varying shares of variable renewable energy (VRE). We find that the total value of energy storage typically increases with VRE shares, but any increase in the relative value of longer storage durations over time depends on the region and grid mix. Some regions see incremental value increasing notably, up to 20–40 h in 2050, while others do not. The negative effect of lower roundtrip efficiency on value is also found to be scenario-dependent, with the energy value in higher VRE scenarios being less sensitive to roundtrip efficiency and more supportive of longer storage durations. Long-duration storage value and deployment potential are a function of evolving electricity sector infrastructure, markets, and policy, making it critical to consistently revisit potential long-duration storage contributions to the grid.

14 SOLAR ENERGY↗
DOE OSTI2020

Technical Gaps in Hanford High-Level Waste Solids Settling Behavior and Settling Time Evaluation for Direct Feed High-Level Waste (DFHLW) Operations

Settling of high-level waste (HLW) solids in process vessels is a key conceptual process step in providing HLW feed directly to the Hanford Waste Treatment and Immobilization Plant (WTP) HLW Vitrification Facility. Direct Feed High-Level Waste (DFHLW) is a potential approach to initiating HLW vitrification prior to completing of the WTP Pretreatment Facility. Settling would be used with subsequent supernatant decant to concentrate HLW feed. To support planning for DFHLW, Washington River Protection Solutions (WRPS) requested support from the Pacific Northwest National Laboratory to evaluate the current data set available to predict the time needed for HLW solids to settle, to identify gaps in the understanding and predictive capability of HLW solids waste settling times, and to provide scoping estimates of the potential settling time. Eight technical gaps were identified for predicting settling times and characteristics of the formed sediment layers including: Gap 1: In-Tank Settling Rates Faster than Settling of Laboratory Samples, Gap 2: Effect of Sludge Leaching/Washing on Predicted Settling Times, Gap 3: Predicting Waste Settling from Waste Chemistry (Waste Type), Gap 4: Predicting Waste Settling from Particle Size and Density Distributions (PSDDs), Gap 5: Insufficient Laboratory and In-Tank Settling Data to Represent Hanford Waste, Gap 6: Methods for Real-Time, In-Tank Tracking of Settling, Gap 7: Prediction of Sediment Erosion Resistance as a Function of Settling Time, and Gap 8: Prediction of Sediment Solids Content as a Function of Settling Time. In addition to the data gaps, an overarching observation of the settling rate and settled layer data is the significant variation in behavior. At similar solids concentrations, settling rates can vary by as much as 3 orders of magnitude depending on the source waste tank, and significantly different settling rates are noted between laboratory and in situ tests for the same waste tank. The range of average solids concentration in existing HLW sediment, which may have been quiescent for decades, can vary from less than 7 wt% to greater than 74 wt% solids. The shear strengths (or yield stresses) measured on laboratory samples range from less than 27 Pa to greater than 6400 Pa. These variations can challenge process planning for the application of a settle/decant process for DFHLW. This report describes the significance of the gaps to the settle/decant process and presents uncertainties by way of examples. Potential technical approaches for resolving these gaps are described and the estimated difficulty in resolving these gaps is evaluated. Based on the significance of the gap and the difficulty of resolution, recommendations are made to address specific gaps. Scoping estimates of the potential settling times for DFHLW solids have been made based on the existing data set with its associated gaps. Depending on the process vessel depth and final sediment concentration, substantial fractions of the scoping estimate results for settling times for characterized HLW exceed the 2-week period that has been previously assumed for process planning. There is also significant disparity, potentially greater than a factor of 5000 difference, in the estimated settling times depending on process vessel depth and final sediment solid concentration. This variation in results underscores the significance of the identified gaps and uncertainties with respect to process planning for utilizing settle/decant operations for DFHLW.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Innovating the next generation of commercial smart building software

Nearly 30% of commercial building energy use is wasted due to equipment faults and HVAC controls problems. The result is increased emissions, compromised comfort and productivity, and less reliable coordination of building power needs with a clean grid. The energy impact alone represents $17 billion in potential savings. Today’s smart building software provides a robust solution to address these operational deficiencies. Energy management and information systems (EMIS) are saving up to 9% on average, with two-year paybacks. They are being incorporated into energy management processes, commissioning services, and utility programs. As effective as they are, two barriers prevent even deeper benefits; limited personnel to fix problems once they are identified, and the expense and time to manually implement changes in control systems. In partnership with the research community, the EMIS industry is developing new capabilities to overcome these barriers. Moving beyond siloed products for either fault detection and diagnostics, or optimal control, these new capabilities empower users to not only automatically identify faults, but also to push corrective action, and control improvements to their buildings. In this paper, several areas for enhancements are documented: ‘one-time’ correction of faults such as setpoints, schedules, and economizer lockouts; short-term active testing for automated proportional integral derivative (PID) loop tuning and functional testing; and continuous supervisory control for demand flexibility and year-round efficiency. Results are presented from a pair of partner implementations out of a dozen providers integrating these enhancements into their products, including field tests from across the country, and insights into operator acceptance and integration into operations and maintenance practices.

Casillas, Armando↗
DOE OSTI2022

Assessing models of force-dependent unbinding rates via infrequent metadynamics

Protein–ligand interactions are crucial for a wide range of physiological processes. Many cellular functions result in these non-covalent “bonds” being mechanically strained, and this can be integral to proper cellular function. Broadly, two classes of force dependence have been observed—slip bonds, where the unbinding rate increases, and catch bonds, where the unbinding rate decreases. Despite much theoretical work, we cannot predict for which protein–ligand pairs, pulling coordinates, and forces a particular rate dependence will appear. Here, we assess the ability of MD simulations combined with enhanced sampling techniques to probe the force dependence of unbinding rates. We show that the infrequent metadynamics technique correctly produces both catch and slip bonding kinetics for model potentials. We then apply it to the well-studied case of a buckyball in a hydrophobic cavity, which appears to exhibit an ideal slip bond. Finally, we compute the force-dependent unbinding rate of biotin–streptavidin. Here, the complex nature of the unbinding process causes the infrequent metadynamics method to begin to break down due to the presence of unbinding intermediates, despite the use of a previously optimized sampling coordinate. Allowing for this limitation, a combination of kinetic and free energy computations predicts an overall slip bond for larger forces consistent with prior experimental results although there are substantial deviations at small forces that require further investigation. Furthermore, this work demonstrates the promise of predicting force-dependent unbinding rates using enhanced sampling MD techniques while also revealing the methodological barriers that must be overcome to tackle more complex targets in the future.

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

Precursor Analysis Report: Industroyer Targeting Ukraine Electric Power Transport Utility (Ukrenergo) 2016

The Industroyer Targeting Ukraine Electric Power Transport Utility (Ukrenergo) 2016 Precursor Analysis Report leverages publicly available information about the December 2016 cyber attack against the Ukrainian Ukrenergo electric transmission utility and catalogs anomalous observables for each technique employed in the attack. This analysis is based upon the methodology of the Cybersecurity for the Operational Technology Environment (CyOTE) program. Industroyer is a modular malware framework designed to deploy several Industrial Control System (ICS) protocol-specific attack payloads to disrupt electricity distribution. Adversaries deployed Industroyer within the target network on a Microsoft Windows endpoint capable of directly manipulating or communicating with ICS. Industroyer abuses the functionality of a targeted ICS’s legitimate control system to achieve its intended impact. Adversaries likely first gained access to Ukrenergo enterprise networks in early 2016 after a successful spearphishing campaign against organizations in the electric power sector. Adversaries then began capturing credentials beginning on 1 December 2016. This allowed access to the ICS environment at the Pivnichna electric transmission substation outside Kyiv through a device dual-homed on the Information Technology (IT) and ICS networks. Adversaries conducted discovery, targeting, and access to this device using information and previously captured credentials from compromised enterprise IT machines. Finally, the adversaries deployed and launched the Industroyer malware just before midnight on 17 December. By midnight, Ukrenergo had lost control of a targeted substation, resulting in electric power outages for over an hour in the city of Kyiv and the Kyiv region. Researchers and analysts identified 31 unique techniques (used in a sequence of 33 steps) utilized during the attack with a total of 846 observables using MITRE ATT&CK® for Industrial Control Systems. The CyOTE program assesses observables accompanying techniques used prior to the triggering event to identify opportunities to detect malicious activity. If observables accompanying the attack techniques are perceived and investigated prior to the triggering event, earlier comprehension of malicious activity can take place. Twenty-nine of the identified techniques used during the Industroyer cyber attack were precursors to the triggering event. Analysis identified 548 observables associated with these precursor techniques, 353 of which were assessed to have an increased likelihood of being perceived in the 300 days preceding the triggering event. The response and comprehension time could have been reduced if the observables had been identified earlier. The information gathered in this report contributes to a library of observables tied to a repository of artifacts, data sources, and technique detection references for practitioners and developers to support the comprehension of indicators of attack. Asset owners and operators can use these products if they experience similar observables or to prepare for comparable scenarios.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2023

Cybersecurity Platform and Certification Framework Development for Extreme Fast Charging (XFC)-Integrated Charging Ecosystem (Final Project Report)

This report summarizes a pioneering effort in Electric Vehicle charging infrastructure ecosystem cybersecurity requirements, assessment methodologies, functional verification, as well as embodiment of the key technologies in the form of hardware and software tools being made available to the public. EPRI led a team of experts, as well as a stakeholder coalition encompassing all key actors in the EV charging infrastructure ecosystem that includes eXtreme Fast Charging (XFC) equipment (defined as 200kW or above). EV charging infrastructure in the United States is a patchwork of networks that have continued to grow organically and have been designed to serve the charging needs of the EV owners, who are their customers. In doing so, each network provider, as well as their connected entities such as the cloud Electric Vehicle Service Providers or EVSPs, utility back office, utility AMI networks, payment networks, as well as Original Equipment Manufacturer (EV manufacturer) telematics networks, have designed systems that may work well individually, but no single entity is responsible for the entire ecosystem to be secure in terms of data exchange. Furthermore, there is no uniformity in how each actor has implemented the cybersecurity requirements since no system-wide cybersecurity requirements existed prior to this project. The final project report describes the technical approach guided by the EV charging infrastructure cybersecurity working group, convened specifically for this project. The technical approach included definition of requirements at the ecosystem level, treated as a ‘system of systems’, and then passed down to individual systems (EVSE, EV, cloud EVSP, utility, and the payment networks), followed by developing the cybersecurity risk and vulnerability assessment methods, that were later applied to real-world cyber-physical systems at EPRI, ANL, and NREL laboratories, to validate both the process and the results. Finally, in a spotlight over the most vulnerable equipment, which is the EV charge station (AC or DC), the team developed a multi-layer cybersecurity implementation in the embedded domain embodied by the open-source Secure Network Interface Card (SNIC) demonstrating the various ways in which the infrastructure can be secured protecting against the identified attack surfaces. Finally, the entire process of EV infrastructure cybersecurity assessment was encapsulated in the Electric Vehicle Charging Cybersecurity Management (EVC2M) online GUI-based tool, that is expected to be released to the public. The report presents the objectives, the technical approach, the key results, and recommendations for future work.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

A new coupled multiphase flow–finite strain deformation–fault slip framework for induced seismicity

Production of hydrocarbons and water from subsurface reservoirs are known to cause permanent deformation of the reservoir and seismicity along faults both of which are detrimental to sustainable development of natural resources. Most of the prior studies on understanding fluid flow-induced plasticity and seismicity have focused on one or the other phenomenon due to the numerical difficulty associated with simultaneous modeling of the two failure phenomena because stress and deformation evolve non-linearly in both plasticity and seismicity. However, in reservoirs undergoing long-term production, plastic failure can alter the stress paths of points on a fault such that the onset, location, and magnitude of actual seismic events can no longer be predicted by a poroelastic simulation due to inaccurate stress and deformation history. We present a computational framework for coupled multiphase flow, finite strain poroplastic deformation, and dynamic fault slip and use it to understand the impact of plastic deformation on the onset, location, and magnitude of induced fault slip events. We evaluate the impact of plasticity on reservoir pressure, deformation, induced stress, and fault slip by comparing infinitesimal strain elastic and finite strain poro-elastoplastic models. For real-world applications, we consider different scenarios where the reservoir is either mechanically weaker or stronger than the caprock. We analyze the stress evolution as a function of the change in reservoir pressure to understand the role of contrast in reservoir and caprock elastic moduli on geomechanical stability. The results show that the poroplastic reservoir exhibits larger vertical deformation and delayed slip than the poroelastic reservoir after the same amount of oil production. For the same amount of pressure drop, a reservoir with a smaller modulus than the caprock displays a larger vertical displacement and an earlier onset of both plastic failure and fault slip. For a reservoir with a larger modulus than the caprock, vertical displacement is larger on the reservoir top boundary and smaller on the ground surface, and a higher pressure drop is needed to induce plastic failure and fault slip.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗