Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model types”

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 235 records · Page 13

DeePMD-kit v2: A software package for deep potential models

DeePMD-kit is a powerful open-source software package that facilitates molecular dynamics simulations using machine learning potentials known as Deep Potential (DP) models. This package, which was released in 2017, has been widely used in the fields of physics, chemistry, biology, and material science for studying atomistic systems. The current version of DeePMD-kit offers numerous advanced features, such as DeepPot-SE, attention-based and hybrid descriptors, the ability to fit tensile properties, type embedding, model deviation, DP-range correction, DP long range, graphics processing unit support for customized operators, model compression, non-von Neumann molecular dynamics, and improved usability, including documentation, compiled binary packages, graphical user interfaces, and application programming interfaces. This article presents an overview of the current major version of the DeePMD-kit package, highlighting its features and technical details. Additionally, this article presents a comprehensive procedure for conducting molecular dynamics as a representative application, benchmarks the accuracy and efficiency of different models, and discusses ongoing developments.

97 MATHEMATICS AND COMPUTING↗

Deep learning inversion of gravity data for detection of CO 2 plumes in overlying aquifers

In this work, we developed an effective U-Net based deep learning (DL) model for inversion of surface gravity data on a rectangular grid to predict 2-D high-resolution subsurface CO 2 distribution along a vertical cross-section due to CO 2 leakage through a wellbore within a deep CO 2 storage reservoir. We used synthetic data to model two types of CO 2 leakage scenarios: one CO 2 plume in a shallow aquifer (single plume case), and two plumes present at different depths (double plume case). The 3-D synthetic plume samples were created by sampling among predetermined CO 2 plume depths, saturations, and volumes. The corresponding surface gravity data on a rectangular grid were generated by a 3-D forward model. The U-Net model detected 72% of single-plume samples, and one or both plumes in 75% of double-plume samples. Most of the undetected single plumes have small gravity field strengths below the typical noise level of 5 μGal. This model generated reproducible, reliable predictions with acceptable errors and demonstrated improved spatial resolution over the conventional least-squares inversion. In contrast to the conventional least-squares inversion, which often overestimates the size of its target and underestimates its density, this U-Net model accurately delineated the boundary of a target. Furthermore, this DL inversion detected deep, small, or low saturation CO 2 plumes that are often more difficult to resolve with conventional gravity inversion methods. We note the limitations of this feasibility study, including the use of synthetic data with regular CO 2 plume shapes, and the prediction of a 2-D plume cross-section rather than the full 3-D plume, as well, we recognize the lower detection fraction for double-plume scenarios. Nevertheless, this study demonstrates that DL gravity inversion is a promising and potentially superior method to conventional least-squares inversion. Our U-Net based deep learning inversion approach may be adapted for inversion of other types of geophysical data. DL inversion can facilitate near real-time monitoring of geologic carbon sequestration to provide site operators with prompt information about subsurface CO 2 distribution for risk management and mitigation.

58 GEOSCIENCES↗

Phase-field framework with constraints and its applications to ductile fracture in polycrystals and fatigue

Abstract Modeling of ductile fracture in polycrystalline structures is challenging, since it requires integrated modeling of cracks, crystal plasticity, and grains. Here we extend the typical phase-field framework to the situations with constraints on the order parameters, and formulate two types of phase-field models on ductile fracture. The Type-I model incorporates three sets of order parameters, which describe the distributions of cracks, plastic strain, and grains, respectively. Crystal plasticity is employed within grain interiors accommodated by J 2 plasticity at grain boundaries. The applications of the Type-I model to single crystals and bicrystals demonstrate the influences of grain orientations and grain boundaries on crack growth. In the Type-II model, J 2 plasticity is assumed for the whole system and grain structures are neglected. Taking advantage of the efficiency of the fast Fourier transform, our Type-II model is employed to study low cycle fatigue. Crack closure and striation-like patterning of plastic strain are observed in the simulations. Crack growth rate is analyzed as a function of the J -integral, and the simulated fatigue life as a function of plastic strain agrees with the Coffin–Manson relation without a priori assumption.

36 MATERIALS SCIENCE↗

Role of Oxidizing Conditions in the Dispersion of Supported Platinum Nanoparticles Explored by Ab Initio Modeling

Achieving fine control over the dispersion of supported platinum nanoparticles (Pt) is a promising avenue to enhancing their catalytic activity and selectivity. Experimental observations suggest that exposing ceria-supported Pt nanoparticles to O 2 at 500 °C promotes their dispersion into smaller particles and eventually single atoms. The exact role of oxygen in this process is not yet well understood. Past density functional theory studies of ceria-supported Pt have typically narrowed their scope to single atoms and clusters of a few atoms. Herein, we combine several approaches and types of models in a consistent atomistic framework to evaluate the relative stability of ceria-supported Pt as a function of the degree of oxidation of Pt and of the particle size, ranging from single atoms to nanoparticles 1.5 nm of diameter. Finally, we find that the largest nanoparticles remain the thermodynamically most stable species on the lowest energy facet of ceria even under oxidizing conditions, suggesting that stronger adsorption sites are required to stabilize smaller clusters and single atoms and promote oxidative dispersion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

PyTorch Implementation of Log-Additive Convolutional Neural Networks

This code is a collection of python code that defines, trains, and tests Log-Additive Convolutional Neural Networks. The model components and training routine are based on the PyTorch python library. The code implements the Log-Additive Convolutional Neural Networks as described in Pagendam et al. 2023. In addition to the Log-Additive Convolutional Neural Networks, this library also defines the Log-Normal Density loss function as described in Pagendam et al. 2023. Code from this paper is not publicly available, so the Pytorch implementation of this type of model is unique to this library.

Callis, Skylar↗

Data Requirements for Application of Risk-Based Dynamic Contingency Analysis to Evaluate Hurricane Impact to Electrical Infrastructure in Puerto Rico

This paper presents a risk-based dynamic contingency analysis framework that was used to evaluate the hurricane impact to electrical infrastructure in Puerto Rico. PNNL developed a scalable risk-based framework for identifying high-voltage transmission resilience improvements by classifying and prioritizing high-risk power grid contingencies (system failures) under hurricane impact. The risk-based framework is founded on grid outage definitions with their associated probabilities of occurrence from hurricane events, in combination with an impact assessment derived from detailed dynamic cascading failure analysis. This paper focuses on a discussion around data requirements for transmission resilience planning for hurricane events, derived from the development of the risk-based framework and its application to Puerto Rico. This paper launches an important first step in encouraging the engineering community and power system industry to move towards establishing resilience planning as a routine practice. Since actual results for Puerto Rico contain sensitive information, sample simulation results will be used to illustrate the data requirements and risk-based dynamic cascading framework on the Puerto Rico power grid, as well as demonstrate the potential for such a simulation framework. The paper includes a discussion on the lessons learned, importance and need for improved datasets that are not usually considered in traditional power system planning. The paper will also elaborate on how the scalable simulation framework and datasets might be expanded to larger footprints and leveraged for modelling other types of natural disasters.

DCAT, Puerto Rico, hurricane, Power System Stabili↗

Hydrogen Bond Network Disruption by Hydration Layers in Water Solutions with Salt and Hydrogen-Bonding Polymers (PEO)

A mean field theory model describing the interaction of ion hydration layers with the network of hydrogen bonds of both water and the nonionic polymer poly(ethylene oxide) (PEO) is presented. The predictions of the model for types and statistics of hydrogen bonds, the number of water molecules bound to PEO, or their dependence on temperature are successfully verified from all-atom simulations at different NaCl and PEO concentrations. Furthermore, our simulations show that the binding of cations to PEO increases monotonically with salt concentration, in agreement with recent experimental results, through a mechanism in which the sum of the number of bound water and cations is independent of salt concentration. Here, the model introduced is general and can describe any salt or hydrogen-bond-forming polymer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

HydraGNN v4.0

The new version of HydraGNN v4.0 provides additional core capabilities, such as: Inclusion of multi-body atomistic cluster expansion MACE, polarizable atom interaction neural network PAINN, and equivariant principal neighborhood aggregation (PNAEq) among the message passing layers supported -Inclusion of graph transformers to directly model long-range interactions between nodes that are distant in the graph topology Integration of graph transformers with message passing layers by combining the graph embedding generated by the two mechanisms, which allows for an improved expressivity of the HydraGNN architecture Improved re-implementation of multi-task learning (MTL) to allow its use for stabilized training across imbalanced, multi-source, multi-fidelity data Introduction of multi-task parallelism, a newly proposed type of model parallelism specifically for MTL architectures, which allows to dispatch different output decoding heads to different GPU devices Integration of multi-task parallelism with pre-existing distributed data parallelism to enable a 2D parallelization for distributed training Improved portability of the distributed training across Intel GPUs, which has been testes on ALCF exascale supercomputer Aurora Inclusion of 2-level fine-grained energy profilers portable across NVIDIA, AMD, and Intel GPUs to monitor the power and energy consumption associated with different functions executed by the HydraGNN code during data pre-load and training Restructuring of previous examples and inclusion of new sets of examples to illustrate the download, preprocess, and training of HydraGNN models on new large-scale open-source datasets for atomistic materials modeling (e.g., Alexandria, Transition1x, OMat24, OMol25)

Lupo Pasini, Massimiliano [Oak Ridge National Labo↗

Assessing Cognitive Impacts of Errors from Machine Learning and Deep Learning Models: Final Report

Due to their recent increases in performance, machine learning and deep learning models are being increasingly adopted across many domains for visual processing tasks. One such domain is international nuclear safeguards, which seeks to verify the peaceful use of commercial nuclear energy across the globe. Despite recent impressive performance results from machine learning and deep learning algorithms, there is always at least some small level of error. Given the significant consequences of international nuclear safeguards conclusions, we sought to characterize how incorrect responses from a machine or deep learning-assisted visual search task would cognitively impact users. We found that not only do some types of model errors have larger negative impacts on human performance than other errors, the scale of those impacts change depending on the accuracy of the model with which they are presented and they persist in scenarios of evenly distributed errors and single-error presentations. Further, we found that experiments conducted using a common visual search dataset from the psychology community has similar implications to a safeguards- relevant dataset of images containing hyperboloid cooling towers when the cooling tower images are presented to expert participants. While novice performance was considerably different (and worse) on the cooling tower task, we saw increased novice reliance on the most challenging cooling tower images compared to experts. These findings are relevant not just to the cognitive science community, but also for developers of machine and deep learning that will be implemented in multiple domains. For safeguards, this research provides key insights into how machine and deep learning projects should be implemented considering their special requirements that information not be missed.

97 MATHEMATICS AND COMPUTING↗

Application of artificial intelligence methods in the international roughness index prediction of rigid and composite pavements: a systematic review

The International Roughness Index (IRI) is a widely adopted metric for quantifying pavement roughness, directly influencing vehicle safety, ride comfort, and overall roadway performance. In recent years, the use of Machine Learning (ML) models for IRI prediction has gained momentum, with the goal of improving the allocation of maintenance and rehabilitation resources by enabling accurate assessments of pavement conditions. Most prior reviews, however, have concentrated on flexible pavements, leaving a notable gap regarding rigid and composite pavements. To address this gap, the present study conducts a systematic review of Artificial Intelligence (AI) methods applied to IRI prediction for rigid and composite pavements. Literature published between 2004 and 2025 is synthesized to highlight prevailing trends, methodological contributions, and directions for future research. Particular attention is given to the types of models employed, the datasets used for training and validation, and the role of input variables and data-processing strategies. Across the included studies, ensemble learning methods (especially gradient boosting variants such as XGBoost), artificial neural networks, and hybrid architectures frequently achieved high predictive skill, with several models reporting test-set coefficients of determination approaching 0.9–0.96, indicating strong potential for capturing the influence of traffic, pavement structure, and climatic factors. Since these results are obtained from heterogeneous datasets and evaluation protocols, they are interpreted qualitatively rather than as strict cross-study rankings. Analysis of input variables revealed that pavement age and initial IRI were included in 91% (21 of 23) and 78% (18 of 23) of studies, respectively. Climatic variables such as the freezing index appeared in 57% (13 of 23), while traffic-related factors were considered in 65% (15 of 23). The findings underscore the importance of standardized, high-quality datasets, such as those from the Long-Term Pavement Performance (LTPP) program, along with data consistency, model interpretability, computational efficiency, and replicability in enhancing IRI prediction. Future research should focus on incorporating input variable selection techniques to identify the most influential predictors, thereby improving accuracy and robustness. Integrating these approaches with advanced non-linear data-driven models, coupled with robust hyperparameter optimization, holds considerable promise for strengthening the reliability of IRI prediction and supporting resilient pavement management strategies.

42 ENGINEERING↗

Probing bottom-associated production of a TeV scale scalar decaying to a top quark and dark matter at the LHC

Abstract A minimal non-thermal dark matter model that can explain both the existence of dark matter and the baryon asymmetry in the universe is studied. It requires two color-triplet, iso-singlet scalars with$$ \mathcal{O}\left(\textrm{TeV}\right) $$ O TeV masses and a singlet Majorana fermion with a mass of$$ \mathcal{O}\left(\textrm{GeV}\right) $$ O GeV . The fermion becomes stable and can play the role of the dark matter candidate. We consider the fermion to interact with a top quark via the exchange of QCD-charged scalar fields coupled dominantly to third generation fermions. The signature of a single top quark production associated with a bottom quark and large missing transverse momentum opens up the possibility to search for this type of model at the LHC in a way complementary to existing monotop searches.

Physics↗

Comparison on the properties of ITZs in fly ash-based geopolymer and Portland cement concretes with equivalent flowability

Highlights: • Modelled ITZs can effectively avoid the effects caused by different factors and test locations of ITZ. • ITZ properties of PC and geopolymer concretes are compared based on equivalent flowability of binders. • Geopolymer ITZ exhibits stronger bonding than that of PC, and higher mechanical properties than geopolymer paste. • Strong ITZs in geopolymer concrete are attributed to the layer of gel-rich paste with dense microstructures. • Elastic modulus and hardness properties of ITZs in PC concrete are increased due to the CH crystals. This paper aims to compare the properties of interfacial transition zones (ITZs) in Portland cement (PC) concrete and geopolymer concrete. Portland cement and geopolymer pastes were designed with the equivalent flowability to provide similar mix and casting condition of ITZs. Two types of modelled ITZs were prepared to facilitate the nanoindentations across ITZs, microstructural characterization, and comparison on the properties of ITZs with less influential factors. The results showed that the interfacial bonding of ITZs between geopolymer matrix and aggregate is relatively stronger than the counterpart in the PC concrete. There is a high amount of crystalline hydration products in the ITZs of PC concrete, but a layer gel-rich paste with denser microstructures in the ITZs of geopolymer concrete. Additionally, the interface morphology and nanoindentation analysis indicate that the property of ITZs in the modelled geopolymer concrete is not poorer than that of the corresponding geopolymer paste.

36 MATERIALS SCIENCE↗

Involute Working Group – FSI Analysis of Fuel Plates Using Finite Volume and Finite Element Methods

The three involute plate research reactors RHF, HFIR, and FRM II have expressed an interest in using computational software to carry their steady-state safety analysis. Since these tools represent a significant departure from the methods used currently (one-dimensional), the acceptability of the new approach by regulators requires thorough verification and validation of these tools. Therefore, Argonne National Laboratory and the three involute-plate reactors formed an informal group called the Involute Working Group aiming at qualifying computational tools to perform steady-state safety analysis. The present report focuses on a comparison of finite volume and finite element methods to model solids in fluid-structure interaction problems with the goal to estimate the coolant flow-induced fuel plate deflections obtained with the two methods. The finite volume method will be obsoleted in STARCCM+ by the end of 2021, nevertheless, this evaluation is important because the method was used by ANL researchers to model the response of the fuel plates, despite its drawbacks, which are discussed in the report. It was essential to check how those estimates compare to the results obtained with the finite element method that is considered superior for structural analysis. Various geometries, i.e., flat, cylindrical and circle-involute fuel plates, as well as coolant flow speed, were considered. The comparison shows that, independently of the plate geometry, the finite volume method significantly underestimates the deflection as compared to finite element method for coarser meshes. When the discretization is developed as a result of a mesh sensitivity study using finite element method, the result obtained using finite volume method can be a few times smaller than the corresponding finite element method solution. A code-to-code comparison , between STAR-CCM+ and LS-DYNA was included in the analysis. Within the LS-DYNA models, two types of finite element formulations were used: solid and shell finite elements. Mesh sensitivity study showed that both approaches converge to a similar value that was obtained with STAR-CCM+ finite element solver. The evaluation of the computational solvers was extended by adding two benchmark cases from the STAR-CCM+ Verification Suite and presented in the Appendix A. The selected cases are: (1) bending of a cantilever beam under external load, and (2) cylindrical shell deformation analysis, known in the literature as ‘Scordelis-Lo roof’. The problems were solved with finite volume, and finite element methods, and the results confirmed the previously discussed findings. The analysis shows that the finite element solver is superior to the finite volume solver in terms of representation of model geometry and estimating the structural behavior of fuel plates. Depending on the ratio of the load to the flexibility of the plate, the finite volume solver can greatly under- or overestimate the structural response if a very carefully selected mesh is not used.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Binding of Phosphate Species to Ca 2+ and Mg 2+ in Aqueous Solution

Phosphate derivatives and their interaction with metal cations are involved in many important biological phenomena, so an accurate characterization of the phosphate-metal interaction is necessary to properly understand the role of phosphate-metal contacts in mediating biological function. Herein, we improved the standard 12-6 Lennard-Jones (LJ) potential via the usage of the 12-6-4 LJ model, which incorporates ion-induced dipole interactions. Via parameter scanning, we fine-tuned the 12-6-4 LJ polarizability values to obtain accurate absolute binding free energies for the phosphate anions H 2 PO 4 – , HPO 4 2– , PO 4 3– coordinating with Ca 2+ and Mg 2+ . First, we modified the phosphate 12-6-4 LJ parameters to reproduce the solvation free energies of the series of phosphate anions using the thermodynamic integration (TI) method. Then, using the potential mean force (PMF) method, the polarizability of the metal-phosphate interaction was obtained. We show that the free energy profiles of phosphate ions coordinated to Ca 2+ and Mg 2+ generally show similar trends at longer metal–phosphate distances, while the absolute binding energy values increased with deprotonation. The resulting parameters demonstrate the flexibility of the 12-6-4 LJ-type nonbonded model and its usefulness in accurately describing cation–anion interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Leptonic cascade decays of a heavy Higgs boson through vectorlike leptons at the LHC

We demonstrate the potential of fully leptonic cascade decays of a heavy neutral Higgs boson through vectorlike leptons as a simultaneous probe for extended Higgs sectors and extra matter particles at the LHC. The processes we explore are unique in that their event topologies lead to di-boson-like leptonic final states with a lepton pair which does not reconstruct the mass of a gauge boson. By recasting existing $2\ell$ + $E^{\text{miss}}_{R}$ and $3/4\ell$ searches channels using run2 data from the LHC we obtain model independent bounds on the masses of heavy scalars and vectorlike leptons and use these results to explore future prospects at the HL-LHC. Our results can be directly applied to any kind of new physics scenarios sharing the final states and the event topology. For concreteness, we apply our results to a benchmark scenario: a two Higgs doublet model type-II augmented with vectorlike leptons. Remarkably, even with current data the sensitivity of our analysis shows a reach for masses of a heavy neutral Higgs and vectorlike leptons up to 2 TeV and 1.5 TeV, respectively. Even for low tan β ≳ 1, the analysis retains sensitivity to heavy Higgs masses slightly above 1 TeV. The future sensitivities at the HL-LHC extend the reach for heavy Higgses and new leptons to 2.7 TeV and 2 TeV, respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Machine learning-driven descriptions of protein dynamics at solid-liquid interfaces

This chapter has described how ML has enabled quantitative analysis of HS-AFM data to discover the physical phenomena governing protein dynamics and ordering at solid-liquid interfaces. The research detailed in this chapter modeled the rotation models of protein nanorods, the discovery of which would otherwise not be possible. By tracking the trajectories of individual protein rods from frame to frame, it was possible to model Brownian type motion and behaviors and Levy-flight dynamics that had not previously been shown. We also described the application of the Python package AtomAI, which has been developed specifically to analyze and extract physical phenomena, providing exemplar code for training an ensemble of deep neural networks to produce the semantic segmentation of AFM data and functions for encoding and decoding local environments. We last described a combinatorial approach to analyze very noisy data with a densely covered substrate where the emergence of order for the protein liquid crystals could be elucidated. By combining the methods from Case 1 and 2, it was possible to obtain the center of mass and angle for each rod in the images and track the assembly of the rods over time into a 2D liquid crystal array on the surface of mica.

protein dynamics, solid-liquid interfaces, atomic ↗

High-temperature chromium diffusion in austenitic stainless steel: Ab initio molecular dynamics simulations

Chromium self-diffusion through stainless steel (SS) matrix and along grain boundaries is an important mechanism controlling SS structural materials corrosion. For this work, Cr diffusion in austenitic SS was simulated using canonical ab initio molecular dynamics with realistic models of type-316 SS bulk, with and without Cr vacancies, and a low-energy Σ3 twin boundary typically observed at active corrosion sites. Cr self-diffusion coefficients at 750 and 850 °C calculated using Einstein’s diffusion equation are 4.2 × 10 –6 and 8.1 × 10 –6 Å 2 ps –1 in pristine bulk, 3.8 × 10 –3 and 5.5 × 10 –3 Å 2 ps –1 in bulk including Cr vacancies, and 9.5 × 10 –2 and 1.0 × 10 –1 Å 2 ps –1 at a Σ3[1 1 1]60° twin boundary.

36 MATERIALS SCIENCE↗

Charging needs for electric semi-trailer trucks

Battery-electric vehicles provide a pathway to decarbonize heavy-duty trucking, but the market for heavy-duty battery-electric semi-trailer trucks is nascent, and specific charging requirements remain uncertain. We leverage large-scale vehicle telematics data (>205 million miles of driving) to estimate the charging behaviors and infrastructure requirements for U.S. battery-electric semi-trailer trucks within three operating segments: local, regional, and long-haul. We model two types of charging - mid-shift (fast) and off-shift (slow) - and show that off-shift charging at speeds compatible with current light-duty charging infrastructure (i.e., =350 kW) can supply 35% to 77% of total energy demand for local and regional trucks with =300-mile range. Megawatt-level speeds are required for mid-shift charging, which make up 44% to 57% of energy demand for long-haul trucks with =500-mile range. However, demand shifts from mid-shift to off-shift charging as the range for battery-electric trucks increases and when off-shift charging is widely available. Finally, we observe geographic trends in charging demand, finding that local trucks have greater demand within urban areas, whereas long-haul trucks have more demand along rural interstate corridors. As the range for battery-electric trucks increases, we show that charging demand shifts from rural to urban locations due to observed vehicle dwell tendencies.

33 ADVANCED PROPULSION SYSTEMS↗