Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modeling and analyzing parasitic parameters in high frequency converters

This research focuses on electromagnetic interference (EMI) / electromagnetic compatibility (EMC) design and analysis in power electronics systems. To limit the EMI under the standards, different methods and strategies are investigated. Parasitic parameters of high frequency (HF) transformer are analyzed using a novel analytical method, finite element method (FEM), and experimental measurements for different structures and windings arrangements. Also, the magnetic field, electric field, electric displacement, and electric potential distribution are simulated and analyzed. Moreover, a high voltage system is considered and analyzed to improve the EMC. The EMI propagation paths are analyzed. The EMI noise level of the system is obtained and compared to the IEC61800-3 standard. To improve the EMC, the parasitic parameters of the transformer, as the main path of EMI circulation, are analyzed and optimized to block the propagation. Furthermore, the geometry structure of the HF transformer is optimized to lower the parasitics in the system. Three pareto-optimal techniques are investigated for the optimization. The models and results are verified by 3D-FEM and experimental results for several given scenarios. Furthermore, the EMC modeling and conducted EMI analysis are developed for a system including an AC-DC-DC power supply (rectifier and dual active bridge (DAB) converter). Moreover, the common mode (CM) EMI noise propagation through the system is discussed and the noise sources and effect of components on the noise are analyzed. Additionally, the CM impedance of different parts of the system and the noise levels are discussed. Finally, EMI attenuation techniques were applied to the system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications

Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring. We have deployed strainmeters and then conducted injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters (30 to 40m depth) was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 n/d over a few days. The signal at two strainmeters at shallow depth were consistent, although the magnitude of the horizontal strains were different reflecting the different radial directions from the well. The signal at a deep strainmeter deployed at reservoir depth was much different, with tensile vertical strains and compressive horizontal strains. These data can be interpreted by inverting poroelastic forward models developed using numerical and analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. This model is fast and can be inverted to estimate reservoir stiffness and geometry. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. The proxy model is periodically updated and refined using the finite element model to ensure accuracy. This approach significantly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. We have shown that the strain tensor in the caprock is sensitive to pressure in the reservoir, boundaries in the reservoir, and pressure in the caprock caused by leaks. These results indicate that coupling strain tensor data with inversion has the potential to help evaluate reservoirs during initial characterization, and to monitor them during the CO2 injection and storage process.

Murdoch, Larry↗
DOE OSTI2022

Determination of transition metal ions in fossil fuel associated wastewaters using chelation ion chromatography

Here, this study outlines the development and subsequent validation of a method using chelation ion chromatography (CIC) pretreatment followed by traditional ion chromatography (IC) and post column UV/vis detection to measure transition metals in fossil fuel wastewaters, such as oil & gas (O&G) brines and coal mine drainage (CMD) waters. Measurement of transition metals is often an important characterization step in the research of environmental and energy systems. IC represents one way to measure these metals with the advantages of being versatile, simple and relatively low cost compared to other analytical methods. However, high concentrations of alkali and alkaline earth metals present in fossil fuel wastewaters will decrease IC detectability of transition metals in these waters. In this study, a CIC method was developed for the analysis of transition metal ions (Fe 3+ , Cu 2+ , Ni 2+ , Zn 2+ , Co 2+ , Mn 2+ , and Fe 2+ ) in fossil fuel associated wastewaters such as Appalachian CMD and O&G wastewaters from the Permian and Bakken shale basins in the United States. CIC system incorporated an on-line chelator column (e.g., the MetPac CC-1) with high selectivity for transition metals over alkali and alkaline earth metals for salt matrix removal prior to transition metal separation and detection. Additional method developments also included acidifying all samples to 2% v/v HCl and using gradient elution rather than isocratic. The recoverability of transition metals in simple salt solutions commonly found in CMD and brine samples (e.g. NaCl, Na 2 SO 4 , CaCl 2 ) using CIC was evaluated and compared to that using traditional IC. Our results found that the CIC system significantly improved transition metal recoveries for samples in 10,000 mg/L CaCl 2 matrix, reaching 87%-108% recovery for all analytes, as opposed to 2-323% recovery in traditional IC. The limits of detection in this study achieved 10.09 – 161.2 μg/L, comparable to reported values in similar IC studies. The developed method was also verified with certified water samples, resulting in 89% - 111% recoveries in samples with higher analyte concentrations (i.e. >4x the LoDs). The developed method achieved 87% -112% recoveries for most analytes in CMD samples and 72%-138% recoveries for Bakken shale samples, relative to ICP-MS values. Overall, the current IC method can be a very good screening tool for fast and cheap analysis for transition metals at mg/L level, to facilitate selection of samples for more detailed ICP-MS analysis.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Optimal Realization of Yang–Baxter Gate on Quantum Computers

Quantum computers provide a promising method to study the dynamics of many-body systems beyond classical simulation. On the other hand, the analytical methods developed and results obtained from the integrable systems provide deep insights on the many-body system. Quantum simulation of the integrable system not only provides a valid benchmark for quantum computers but is also the first step in studying integrable-breaking systems. The building block for the simulation of an integrable system is the Yang–Baxter gate. It is vital to know how to optimally realize the Yang–Baxter gates on quantum computers. Based on the geometric picture of the Yang–Baxter gates, the optimal realizations of two types of Yang–Baxter gates with a minimal number of controlled NOT (CNOT) or gates are presented. It is also shown how to systematically realize the Yang–Baxter gates via the pulse control. The different realizations on IBM quantum computers are tested and compared. It is found that the pulse realizations of the Yang–Baxter gates always have a higher gate fidelity compared to the optimal CNOT or realizations. On the basis of the above optimal realizations, the simulation of the Yang–Baxter equation on quantum computers is demonstrated. Finally, these results provide a guideline and standard for further experimental studies based on the Yang–Baxter gate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Prediction of tissue optical properties using the Monte Carlo modeling of photon transport in turbid media and integrating spheres

Monte Carlo methods are an established technique for simulating light transport in biological tissue. Integrating spheres make experimental measurements of the reflectance and transmittance of a sample straightforward and inexpensive. This work presents an extension to existing Monte Carlo photon transport methods to simulate integrating sphere experiments. Crosstalk between spheres in dual-sphere experiments is accounted for in the method. Analytical models, previous works on Monte Carlo photon transport, and experimental measurements of a synthetic tissue phantom validate this method. We present two approaches for using this method to back-calculate the optical properties of samples. Experimental and simulation uncertainties are propagated through both methods. Both back-calculation methods find the optical properties of a sample accurately and precisely. Our model is implemented in standard Python 3 and CUDA C++ [J. Nickolls, I. Buck, M. Garland, and K. Skadron, ACM Queue 6 , 40 ( 2008 ) ] and is publicly available in Code 1.

Cook, Patrick D. (ORCID:0000000279345428)↗
DOE OSTI2024

CLimate Impact: Determining Etiology thRough pAthways (CLDERA)

Climate impacts have broad economic, health, political, and national security ramifications. Societally relevant impacts are typically farther downstream, are the product of multiple interacting processes, and can arise over small regions and timeframes because their sources are short-term and localized. Short-term forcings (as can be seen in volcanic eruptions, climatic tipping points (e.g., the collapse of rainforests or the disappearance of sea ice), or in increasingly plausible climate interventions) fundamentally possess low signal-to-noise and could benefit from accounting for the multiple conditional processes through which a downstream impact arises. Under the Grand Challenge LDRD CLDERA (CLimate impacts: Discovering Etiology thRough pAthways), we have developed tools to enable downstream impact attribution from geographically and temporally localized source forcings in the climate. CLDERA developed methods that can distinguish how a localized source drives the climate system to respond with particular impacts. The how is embodied in pathways – the spatio-temporally evolving chain of physical processes that connects a source to a series of increasingly distant impacts. Novel analytic methods in pursuit of downstream impact attribution were developed and demonstrated on simulations and observations of the 1991 eruption of Mt. Pinatubo in the Philippines. As described within this report we have • developed stratospheric expertise and aerosol modeling capabilities in E3SM, • created original methods to detect and model pathways from source-to-impact, and • advanced climate attribution through novel methods, cases, and approaches. Further, CLDERA developed a tiered verification process consisting of controlled datasets to prototype, verify, and refine the original method development. CLDERA increased Sandia’s footprint in the climate analytics community and developed new climate collaborations whilst also creating a cadre of climate analysts at Sandia. The products from CLDERA have been extensive with a total of 9 journal articles published, 12 articles submitted and under review, and an additional 8 articles in preparation. We have produced 1750 simulated years and developed 9 code-bases. This report details these accomplishments and serves as a summary of the work completed during the CLDERA Grand Challenge.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Geometric Interpretation of the Cluster Location Problem Part I: Theory

We present a new framing of the seismic location problem using principles drawn from differential geometry. Our interpretation relies upon the common assumption that travel times observed across a network are continuous, differentiable functions of source location. In consequence, travel‐time functions constitute a differentiable map between the source region and a Riemannian manifold. The manifold is said to be the image of the source region embedded in a generally high‐dimension travel‐time vector space. A cluster of events in the source region has an image of discrete points on the manifold, that, except in the simplest cases, cannot be viewed directly. However, it is possible to project the image of a cluster into a tangent space of the manifold for direct visualization. The projection operator can be computed directly from the data without a velocity model, but produces a distorted rendering of the cluster geometry. With a model we can predict the distortions and correct them to estimate cluster geometry. We develop these points with the simplest possible example, one for which direct visualization of the manifold is possible, using the example as an introduction to the relevant concepts from differential geometry in a familiar setting. The tangent space, a local linearization of the manifold, plays a key role. We develop a metric to estimate the limits of linearization, that is, to determine when the curvature of the manifold invalidates the linear assumption. We also examine the interplay of model error, inadequate network geometry, and pick error. We then generalize our results from the simple case to the general case of 3D source regions observed by general networks. Although we do suggest a new “project and correct” method for location, we do not develop it into a practical algorithm. In conclusion, our intention rather is to highlight new analytical methods grounded in differential geometry.

East Pacific Ocean Islands↗
DOE OSTI2021

Data Science Techniques, Assumptions, and Challenges in Alloy Clustering and Property Prediction

Data analytics methods have been increasingly applied to understanding materials chemistry, processing due to the manufacturing approach, and uni-axial and cyclic property relationships in the highly complex space of alloy design. There are several benefits to applying data analytics to this space, including the ability to manage non-linearities in the responses of the alloy attributes and the resulting mechanical properties. However, key difficulties in applying and understanding the results of data analytics include the often lack of reported assumptions and data processing steps necessary to improve interpretation and reproducibility in derived results. In this work, the methods used to generate clustering and correlation analyses for experimental 9% Cr ferritic-martensitic steel data were investigated and the resulting implications for mechanical property predictions were assessed. This work uses principal component analysis, partitioning around medoids, t-SNE, and k-means clustering to investigate trends in composition, processing and microstructure information with creep and tensile properties, building on work done previously using a smaller version of the same dataset. The initial assumptions, preprocessing steps and methods are investigated and outlined in order to depict the fine level of detail required to convey the steps taken to process data and produce analytical results. Here, the variations in the resulting analyses are explored due to the influence of new and more varied data.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Use of machine learning to analyze chemistry card sort tasks

Education researchers are deeply interested in understanding the way students organize their knowledge. Card sort tasks, which require students to group concepts, are one mechanism to infer a student’s organizational strategy. However, the limited resolution of card sort tasks means they necessarily miss some of the nuance in a student’s strategy. Here in this work, we propose new machine learning strategies that leverage a potentially richer source of student thinking: free-form written language justifications associated with student sorts. Using data from a university chemistry card sort task, we use vectorized representations of language and unsupervised learning techniques to generate qualitatively interpretable clusters, which can provide unique insight in how students organize their knowledge. We compared these to machine learning analysis of the students’ sorts themselves. Machine learning-generated clusters revealed different organizational strategies than those built into the task; for example, sorts by difficulty or even discipline. There were also many more categories generated by machine learning for what we would identify as more novice-like sorts and justifications than originally built into the task, suggesting students’ organizational strategies converge when they become more expert-like. Finally, we learned that categories generated by machine learning for students’ justifications did not always match the categories for their sorts, and these cases highlight the need for future research on students’ organizational strategies, both manually and aided by machine learning. In sum, the use of machine learning to analyze results from a card sort task has helped us gain a more nuanced understanding of students’ expertise, and demonstrates a promising tool to add to existing analytic methods for card sorts.

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

Analysis of semivolatile organics in liquid radioactive residue using sorptive stir bar and solvent back-extraction

Current methods for semivolatiles analysis in radioactive samples can produce large volumes of radioactive solvent residue. A method utilizing stir-bar sorptive extraction has been explored in this work for its applicability to radioactive waste samples. This low solvent analytical method may accelerate remediation, minimize hazardous solvent waste, and reduce exposure risk to workers. Organic compounds (polyaromatic hydrocarbons, chlorinated aromatics, and phenolic compounds) were chosen as surrogates for common Liquid Waste System (LWS) contaminants at Savannah River Site (Aiken, SC). Stir-bar extraction parameters (extraction time, matrix modification, and effective pH range) and solvent back extraction parameters (solvent type, volume, and extraction time) were optimized experimentally for the chosen compounds. Affinity of the stir-bar extraction polymer to radionuclides Cs-137 and Am-241 was observed to determine radionuclide concentration effects. The stir-bar method achieved mean recovery of 100 ± 0.7% (1σ), relative to 114 ± 7% using solvent extraction, while reducing weekly method hands-on time by 93.4% and solvent volume consumption by 99.3%. Sensitivity was improved by 378% in simulated tank waste and 278% in real-world LWS matrix, relative to solvent extraction. This work has produced a safe and optimized method for the low solvent analysis of organics in legacy radioactive tank waste by stir-bar sorptive extraction.

GC-MS↗
DOE OSTI2023

PeakDecoder enables machine learning-based metabolite annotation and accurate profiling in multidimensional mass spectrometry measurements

Multidimensional measurements using state-of-the-art separations and mass spectrometry provide advantages in untargeted metabolomics analyses for studying biological and environmental bio-chemical processes. However, the lack of rapid analytical methods and robust algorithms for these heterogeneous data has limited its application. Here, we develop and evaluate a sensitive and high-throughput analytical and computational workflow to enable accurate metabolite profiling. Our workflow combines liquid chromatography, ion mobility spectrometry and data-independent acquisition mass spectrometry with PeakDecoder, a machine learning-based algorithm that learns to distinguish true co-elution and co-mobility from raw data and calculates metabolite identification error rates. We apply PeakDecoder for metabolite profiling of various engineered strains of Aspergillus pseudoterreus, Aspergillus niger, Pseudomonas putida and Rhodosporidium toruloides. Results, validated manually and against selected reaction monitoring and gas-chromatography platforms, show that 2683 features could be confidently annotated and quantified across 116 microbial sample runs using a library built from 64 standards.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Real-Time Screening for Uranium Enrichment by Paper Spray Ionization Mass Spectrometry for Field Applications

A rapid isotope ratio screening technique for uranium enrichment is demonstrated by utilizing paper spray ionization (PSI) high-resolution mass spectrometry (HRMS). Measurements were conducted using an ambient ionization mass spectrometer coupled to a custom-made JEOL PSI attachment apparatus, with de minimus sample preparation requirements. The current method detection limit for individual isotopes (e.g., 235 U and 238 U) is approximately 50 pg, with subsequent optimization expected to further improve U isotopic sensitivity. The PSI analytical method described herein can support rapid analysis (both in-field and in-lab screening) of isotopes-of-interest, as demonstrated by empirical differentiation of depleted uranium (DU) and low enriched uranium (LEU) analytical aliquots. Furthermore, this analytical workflow holds promise for applications in nuclear forensics, international nuclear safeguards, and nonproliferation missions.

Anions↗
DOE OSTI2025

INTEGRATION OF DATA ANALYTICS WITH SYSTEM HEALTH PROGRAMS

Industry equipment reliability and asset management programs are essential elements that help ensure the safe and economical operation of nuclear power plants. The effectiveness of these programs is addressed in several industry developed and regulatory programs. However, these programs have proven to be labor intensive and expensive. There is an opportunity to significantly enhance the collection, analysis, and use of this information to provide more cost-effective plant operation. Additionally, there is an acute industry need to leverage advanced technology to reduce costs and improve operational effectiveness. The goal of this paper is to provide effective and efficient analytical methods and tools to support risk-informed decisions for the equipment reliability and asset management programs at nuclear power plants. This is accomplished by creating a direct bridge between component health/lifecycle data and decision making (e.g., maintenance scheduling and project prioritization). Here we are supporting typical system engineer decisions regarding maintenance activity scheduling and component ageing management. This is performed in a risk-informed context where herein the term “risk” is broadly constructed to include both plant reliability and economics. This framework combines data analytics tools to analyze equipment reliability data with risk-informed methods designed to support system engineer decisions (e.g., maintenance and replacement schedules, optimal maintenance posture) in a customizable workflow. A challenge is that the structure of this workflow strongly depends on the decision that needs to be made, the type of data available, and the constraints that need to be considered. Current methods are designed to provide specific answers to specific problems; however, these methods might prove to be inadequate even when problem settings slightly change (e.g., different types of requirements, additional dependencies between system reliability and economics). We tackled this challenge by designing framework in a flexible and modular fashion such that the user can assemble and customize his/her own workflow that integrates SSC economic lifecycle models (e.g., maintenance and replacement costs), system reliability models, and optimization methods.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Surface functionality and formation mechanisms of carbon and graphene quantum dots

Graphene quantum dots (GQDs) and Carbon dots (C-dots) have been widely studied in recent years due to their structural and optoelectrical properties. These properties have prompted the exploration of the role of these carbon-based materials in many potential applications. This includes solar cells, photodetectors, bioimaging, sensors, batteries and drug delivery. These properties and applications of GQDs and C-dots are highly dependent on their size, shape and surface functionality. In this work, GQDs and C-dots mixtures have been synthesized by an inexpensive wet chemical method by varying the synthesis temperature from 85°, 100° to 115 °C. The surface functionalities of the synthesized carbon-based materials were investigated by several analytical methods. We discovered a higher degree of oxidation at higher temperatures. The mechanism of formation of different sized GQDs and C-dots with different functionalities have been explained with the help of XPS and NEXAFS analysis. The influence of size and surface functionalities on the optical properties of these nanomaterials is analyzed by UV–Vis and PL spectroscopic techniques. Furthermore, this study demonstrates that physicochemical properties of GQDs and C-dots can be controlled by changing the synthesis temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Updates to Relevance Vector Machine: Multiclass Classification, Variable Selection, and Proof-of-Concept Application to Safeguards Fresh Fuel Verification using List-Mode Neutron Collar Data

To expand the capabilities of safeguards authorities to verify the integrity of fresh fuel assemblies, Oak Ridge National Laboratory has retrofit the existing electronics of the JCC-71 uranium neutron coincidence collar, which contains 18 3 He neutron detectors and an external 241 AmLi(α, n) neutron interrogation source arranged to surround a fresh nuclear fuel assembly. The new electronics system allows analysts to record list-mode neutron multiplicity data in addition to the singles and doubles rates that are currently measured. Based on previous proof-of-concept research, analysis of these new data will identify off-normal fuel configurations in an assembly and characterize or localize the specific partial fuel defects. The purpose of this report it to document the analysis algorithm development and then to demonstrate its capability for the safeguards verification of fresh fuel assemblies using list mode neutron collar data. To analyze the complex list-mode data collected with the upgraded uranium neutron collar, multivariate classification algorithms are being developed using a novel classification method, the relevance vector machine. This approach may be applied to multiclass problems to estimate the probability that test data belongs to one of many possible classes of data. In addition, our method identifies the most useful variables/channels for making predictions, which illuminates the basis for the model’s predictions, and this interpretability is largely unique among data analytics methods. Variable selection occurs during model training and parameter tuning and does not need any external hyperparameter tuning routines. Finally, we apply the modified relevance vector machine to a simulated dataset of list-mode neutron collar data generated with the radiation transport code MCNP. The method can correctly identify off-normal fuel configurations, categorize the data according to four fuel defect scenarios, and rank the channels in the data according to prediction utility. For nuclear safeguards applications, it is concluded that this method has the potential to increase the sensitivity and reliability to detect missing fuel rods from a standard 17 x 17 Pressurized Water Reactor (PWR) fresh fuel assembly. Within this analysis, “off-normal” (i.e., missing fuel rods) were correctly classified in 17 simulated test scenarios with one quarter (25%) of the fresh fuel rods missing using a training data set of 58 simulated measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Gas-phase ion-molecule interactions in a collision reaction cell with triple quadrupole-inductively coupled plasma mass spectrometry: Investigations with N 2 O as the reaction gas

Nitrous oxide (N 2 O) was used as a reaction gas to investigate the gas phase ion-molecule interactions using the Agilent 8900 QQQ-ICP-MS. A multi-element standard containing 45 elements with masses ranging from 9 to 208 u was measured in the presence and absence of N 2 O. The main product ion species observed were oxides and nitrides. Comparison of the N 2 O reaction results with similar measurements conducted with O 2 revealed that N 2 O was more effective at forming oxides in general: the elements Cd and Pb were shown to produce oxides with N 2 O where the reaction did not occur with O 2 . Nitrous oxide was also shown to produce a significant amount of nitride species in a few cases. The general reactivity was shown to be consistent with density functional theory (DFT)-predicted reaction enthalpies, such that all predicted exothermic reactions produced product ions at levels at least 1% of the unreacted ion. Our results show that reaction enthalpy is a reasonable predictor of reactivity with N2O on the timescales of the interactions in non-thermal ICP-MS/MS systems. Our work demonstrates the utility of two relatively new platforms (commercial elemental ICP-MS/MS and EMSL Arrows interface to the NWChem program suite), which allows for the study of a large number of elements within a short period. While DFT with the basis sets utilized here is not the most accurate computational method, it is also not computationally expensive and is shown to be suitable for predicting gas phase reactivity in the QQQ-ICP-MS for the majority of ions studied. Here, the ease and rapidity of data collection and DFT calculations has the potential to be very impactful for the identification of targeted reaction chemistries to be leveraged for analytical method development, such as for the inline separation of isobaric interferences from analytes of interest.

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

The ABACUS cosmological N -body code

We present abacus, a fast and accurate cosmological N-body code based on a new method for calculating the gravitational potential from a static multipole mesh. The method analytically separates the near- and far-field forces, reducing the former to direct 1/r 2 summation and the latter to a discrete convolution over multipoles. The method achieves 70 million particle updates per second per node of the Summit supercomputer, while maintaining a median fractional force error of 10 -5 . We express the simulation time-step as an event-driven ‘pipeline’, incorporating asynchronous events such as completion of co-processor work, input/output, and network communication. abacus has been used to produce the largest suite of N-body simulations to date, the abacussummit suite of 60 trillion particles, incorporating on-the-fly halo finding. abacus enables the production of mock catalogues of the volume and resolution required by the coming generation of cosmological surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Invariant surface elastic properties in FCC metals and their correlation to bulk properties revealed by machine learning methods

In this work, we present a combination of machine-learned models that predicts the surface elastic properties of general free surfaces in face-centered cubic (FCC) metals. These models are built by combining a semi-analytical method based on atomistic simulations to calculate surface properties with the artificial neural network (ANN) method or the boosted regression tree (BRT) method. The latter is also used to link bulk properties and surface orientation to surface properties. The surface elastic properties are represented by their invariants considering plane elasticity within a polar method. The resulting models are shown to accurately predict the surface elastic properties of seven pure FCC metals (Cu, Ni, Ag, Au, Al, Pd, Pt). The BRT model reveals the correlations between bulk and corresponding surface properties in terms of invariants, which can be used to guide the design of complex nano-sized particles, wires and films. Finally, by expressing the surface excess energy density as a function of surface elastic invariants, fast predictions of surface energy as a function of in-plane deformations can be made from these model constructs.

36 MATERIALS SCIENCE↗