Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model assignment”

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 145 records · Page 8

In Silico Guidance for In Vitro Androgen and Glucocorticoid Receptor ToxCast Assays

Molecular initiating events (MIEs) are key events in adverse outcome pathways (AOPs) that link molecular chemistry to target biology. As they are based in chemistry, these interactions are excellent targets for computational chemistry approaches to in silico modelling. In this work, we aim to link ligand chemical structure to MIEs for androgen receptor (AR) and glucocorticoid receptor (GR) binding using ToxCast data. This has been done using an automated computational algorithm to perform maximal common substructure searches on chemical binders for each target from the ToxCast dataset. The models developed show a high level of accuracy, correctly assigning 87.20% of AR binders and 96.81% of GR binders in a 25% test set using holdout cross-validation. The 2D structural alerts developed can be used as in silico models to predict these MIEs, and as guidance for in vitro ToxCast assays to confirm hits. These models can target such experimental work, reducing the number of assays to be performed to gain required toxicological insight. Development of these models has also allowed some structural alerts to be identified as predictors for agonist or antagonist behavior at the receptor target. This work represents a first step in using computational methods to guide and target experimental approaches.

Allen, Timothy H.↗

CoRE MOF DB: A curated experimental metal-organic framework database with machine-learned properties for integrated material-process screening

Here, we present an updated version of the Computation-Ready, Experimental (CoRE) Metal-Organic Framework (MOF) database, which includes a curated set of computation-ready MOF crystal structures designed for high-throughput computational materials discovery. Data collection and curation procedures were improved from the previous version to enable more frequent updates in the future. Machine-learning-predicted properties, such as stability metrics and heat capacities, are included in the dataset to streamline screening activities. An updated version of MOFid was developed to provide detailed information on metal nodes, organic linkers, and topologies of an MOF structure. DDEC6 partial atomic charges of MOFs were assigned based on a machine-learning model. Gibbs ensemble Monte Carlo simulations were used to classify the hydrophobicity of MOFs. The finalized dataset was subsequently used to perform integrated material-process screening for various carbon-capture conditions using high-fidelity temperature-swing adsorption (TSA) simulations. Our workflow identified multiple MOF candidates that are predicted to outperform CALF-20 for these applications.

CoRE MOF database↗

Behavioral state resource selection in invasive wild pigs in the Southeastern United States

Elucidating correlations between wild pig ( Sus scrofa ) behavior and landscape attributes can aid in the advancement of management strategies for controlling populations. Using GPS data from 49 wild pigs in the southeastern U.S., we used hidden Markov models to define movement path characteristics and assign behaviors (e.g., resting, foraging, travelling). We then explored the connection between these behaviors and resource selection for both sexes between two distinct seasons based on forage availability (i.e., low forage, high forage). Females demonstrated a crepuscular activity pattern in the high-forage season and a variable pattern in the low-forage season, while males exhibited nocturnal activity patterns across both seasons. Wild pigs selected for bottomland hardwoods and dense canopy cover in all behavioral states in both seasons. Males selected for diversity in vegetation types while foraging in the low-forage season compared to the high-forage season and demonstrated an increased use of linear anthropogenic features across seasons while traveling. Wild pigs can establish populations and home ranges in an array of landscapes, but our results demonstrate male and female pigs exhibit clear differences in movement behavior and there are key resources associated with common behaviors that can be targeted to improve the efficiency of management programs.

54 ENVIRONMENTAL SCIENCES↗

SMBH seeds from dissipative dark matter

The existence of supermassive black holes (SMBHs) with masses greater than ~ 10 9 M ⊙ at high redshift ( z ≳ 7) is difficult to accommodate in standard astrophysical scenarios. We study the possibility that (nearly) totally dissipative self-interacting dark matter (tdSIDM)–in rare, high density dark matter fluctuations in the early Universe — produces SMBH seeds through catastrophic collapse. We use a semi-analytic model, tested and calibrated by a series of N-body simulations of isolated dark matter halos, to compute the collapse criteria and timescale of tdSIDM halos, where dark matter loses nearly all of its kinetic energy in a single collision in the center-of-momentum frame. Applying this model to halo merger trees, we empirically assign SMBH seeds to halos and trace the formation and evolution history of SMBHs. We make predictions for the quasar luminosity function, the M BH -σ v * relation, and cosmic SMBH mass density at high redshift and compare them to observations. We find that a dissipative dark matter interaction cross-section of σ/ m ~ 0.05 cm 2 /g is sufficient to produce the SMBHs observed in the early Universe while remaining consistent with ordinary SMBHs in the late Universe.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimal reconstruction of baryon acoustic oscillations for DESI 2024

Baryon acoustic oscillations (BAO) provide a robust standard ruler to measure the expansion history of the Universe through galaxy clustering. Density-field reconstruction is now a widely adopted procedure for increasing the precision and accuracy of the BAO detection. With the goal of finding the optimal reconstruction settings to be used in the DESI 2024 galaxy BAO analysis, we assess the sensitivity of the post-reconstruction BAO constraints to different choices in our analysis configuration, performing tests on blinded data from the first year of DESI observations (DR1), as well as on mocks that mimic the expected clustering and selection properties of the DESI DR1 target samples. Overall, we find that BAO constraints remain robust against multiple aspects in the reconstruction process, including the choice of smoothing scale, treatment of redshift-space distortions, fiber assignment incompleteness, and parameterizations of the BAO model. We also present a series of tests that DESI followed in order to assess the maturity of the end-to-end galaxy BAO pipeline before the unblinding of the large-scale structure catalogs.

79 ASTRONOMY AND ASTROPHYSICS↗

Hydrostratigraphic Region 4 Model - with HSUPc

Farnsworth Unit (FWU) CO2EOR Eclipse compositional model: Hydrostratigraphic Region 4 model uses the Morrow2 (UNOCAL) relative permeability curve with capillary pressure assigned heterogeneously by hydrostratigraphic unit

Capillary Pressure↗

Extending Symmetry-Preserving Attention Networks (SPANet) for Jet Assignment in Fully Hadronic \(t\bar{t}\) Events

Fully hadronic \(t\bar{t}\) reconstruction requires assigning reconstructed jets to the two top-quark decay branches, \(t\bar{t}\to(Wb)(Wb)\to(qqb)(qqb)\). This is a combinatorially large, symmetry-rich, and frequently underconstrained set-assignment problem due to detector effects and partial reconstructability. SPANet provides a strong symmetry-preserving baseline for this task, but its standard inference eagerly commits to a single hypothesis and does not explicitly treat branch-level reconstructability decisions. We present a two-stage extension of the SPANet pipeline for fully hadronic \(t\bar{t}\) events: (i) SPANet is modified into a proposal model that generates a shortlist of candidate jet assignments; (ii) a custom dual-head set transformer is trained to re-rank these candidates and identify reconstructible branches. Our extension gives modest improvements in candidate-selection efficiency, with SPANet achieving \(69.9\%\) and our extension achieving \(71.1\%\). Finally, an oracle study of the generated shortlist shows a substantial upper bound of \(88.5\%\) on the evaluated fully reconstructable subset, suggesting that the shortlist contains significant residual information that is not fully exploited.

Lisondodi Tada, Mateo [Puerto Rico U., Mayaguez] (↗

Computational Model for Microballistic Perforation of Multilayer Graphene

The peridynamic theory of solid mechanics is applied to the continuum modeling of the impact of small, high-velocity silica spheres on multilayer graphene targets. The model treats the laminate as a brittle elastic membrane. The material model includes separate failure criteria for the initial rupture of the membrane and for propagating cracks. Material variability is incorporated by assigning random variations in elastic properties within Voronoi cells. The computational model is shown to reproduce the primary aspects of the response observed in experiments, including the growth of a family of radial cracks from the point of impact.

36 MATERIALS SCIENCE↗

Predicting molecular dipole moments by combining atomic partial charges and atomic dipoles

The molecular dipole moment ( μ ) is a central quantity in chemistry. It is essential in predicting infrared and sum-frequency generation spectra as well as induction and long-range electrostatic interactions. Furthermore, it can be extracted directly—via the ground state electron density—from high-level quantum mechanical calculations, making it an ideal target for machine learning (ML). Here, we choose to represent this quantity with a physically inspired ML model that captures two distinct physical effects: local atomic polarization is captured within the symmetry-adapted Gaussian process regression framework which assigns a (vector) dipole moment to each atom, while the movement of charge across the entire molecule is captured by assigning a partial (scalar) charge to each atom. The resulting “MuML” models are fitted together to reproduce molecular μ computed using high-level coupled-cluster theory and density functional theory (DFT) on the QM7b dataset, achieving more accurate results due to the physics-based combination of these complementary terms. The combined model shows excellent transferability when applied to a showcase dataset of larger and more complex molecules, approaching the accuracy of DFT at a small fraction of the computational cost. We also demonstrate that the uncertainty in the predictions can be estimated reliably using a calibrated committee model. The ultimate performance of the models—and the optimal weighting of their combination—depends, however, on the details of the system at hand, with the scalar model being clearly superior when describing large molecules whose dipole is almost entirely generated by charge separation. These observations point to the importance of simultaneously accounting for the local and non-local effects that contribute to μ ; furthermore, they define a challenging task to benchmark future models, particularly those aimed at the description of condensed phases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling turbulent impurity transport in the SOL of DIII-D with a reduced model

A novel impurity transport model that approximates SOL turbulence as a fluctuating poloidal electric field is shown to be an acceptable replacement for the traditional approach of assigning an arbitrary radial diffusion coefficient to the impurity ions. The model is implemented in the DIVIMP impurity transport code and applied to an L-Mode tungsten divertor experiment on DIII-D. The poloidal electric field is represented as fluctuating between ±1000 V m –1 based on previous measurements. The resulting intermittent v r = E θ × B T transport causes ions to transport both into the core as well as into the far-SOL. Simultaneous agreement with estimates of the W density just inside the separatrix as well as in the far-SOL is obtained (n W ~ 10 14 m –3 and n W ~ 10 12 m –3 , respectively). Prompt re-deposition of the W ions was necessary to obtain agreement (f redep ~ 99%). Here, we conclude that simulating impurity transport using a physics-based approximation for turbulence in the SOL, versus arbitrarily assigning diffusion coefficients, may enable better reactor scale predictions of core impurity contamination.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamic Temporal Graph Sequence Data for Resilience-Oriented Distribution Network Reconfiguration

This dataset comprises temporal dynamic graph sequences generated from power grid simulations focused on grid reconfiguration to enhance resilience. The simulations model failure propagation under varying conditions, with nodes assigned distinct failure probabilities. For each time step, the dataset captures the evolution of node states (functional or failed) and features critical to grid operations, such as pv_output, load_profile, load_dispatch, dg_output, loss, and voltage. Node types include sources, normal loads, and nodes with specific equipment like PVs, micro turbines, or shunt capacitors. The dataset is structured to support the training of dynamic graph neural networks, facilitating research on node feature prediction and edge dynamics under failure scenarios. Three distinct configurations are included, providing a robust foundation for modeling power grid resilience.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Peridynamic Model for Microballistic Perforation of Multilayer Graphene

The peridynamic theory of solid mechanics is applied to the continuum modeling of the impact of small, high-velocity silica spheres on multilayer graphene targets. The model treats the laminate as a brittle elastic membrane. The material model includes separate failure criteria for the initial rupture of the membrane and for propagating cracks. Material variability is incorporated by assigning random variations in elastic properties within Voronoi cells. Here, the computational model is shown to reproduce the primary aspects of the response observed in experiments, including the growth of a family of radial cracks from the point of impact.

42 ENGINEERING↗

Impact of the plastic deformation microstructure in metals on the kinetics of recrystallization: A phase-field study

The sensitivity of recrystallization kinetics in metals to the heterogeneity of microstructure and deformation history is a widely accepted experimental fact. However, most of the available recrystallization models employ either a mean field approach or use grain-averaged parameters, and thus neglecting the mesoscopic heterogeneity induced by prior deformation. Here in the present study, we investigate the impact of deformation-induced dislocation (subgrain) structure on the kinetics of recrystallization in metals using the phase-field approach. The primary focus here is upon the role of dislocation cell boundaries. The free energy formulation of the phase-field model accounts for the heterogeneity of the microstructure by assigning localized energy to the resulting dislocation microstructure realizations generated from experimental data. These microstructure realizations are created using the universal scaling laws for the spacing and the misorientation angles of both the geometrically necessary and incidental dislocation boundaries. The resulting free energy is used into an Allen-Cahn based model of recrystallization kinetics, which are solved using the finite element method. The solutions thus obtained shed light on the critical role of the spatial heterogeneity of deformation in the non-smooth growth of recrystallization nuclei and on the final grain structure. The results showed that, in agreement with experiment, the morphology of recrystallization front exhibits protrusions and retrusions. By resolving the subgrain structure, the presented algorithm paves the way for developing predictive kinetic models that fully account for the deformed state of recrystallizing metals.

36 MATERIALS SCIENCE↗

Porous Flow Modeling of Axial Gas Redistribution in Fragmented LWR Fuel Rods using MOOSE

Understanding how gas axially redistributes within fragmented fuel pellets is crucial for predicting the behavior of Light Water Reactor (LWR) fuel rods, particularly during transient and accidental scenarios. The time scale of this phenomenon plays a fundamental role in determining the progression and hazard of a Loss Of Coolant Accident (LOCA), especially when high burn-up fuel in a severe state of fragmentation is involved. Here, this study presents a Computational Fluid Dynamics (CFD) model developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE) to predict the time-scale of plenum depressurization in Light-Water Reactor (LWR) fuel rods driven by axial gas transport through fragmented pellets. The model examines the effects of incorporating non-linearities in the friction term by comparing the results with experimental data. The experiment employed surrogate fuel rods containing pellets subjected to mechanical and/or thermal loadings to simulate various severity of cracking, and aimed at studying the influence of fuel conditions on axial gas redistribution. The results of this analysis indicate that under certain flow regime conditions - determined by the value of an equivalent Reynolds number - accounting for the non-linear friction term in Navier-Stokes equations guarantees better predictions for the time-scale of plenum depressurization. Also, the model enabled the simulation of the pressure decay by assigning distinct permeability values to each pellet instead of a single uniform value. Multiple simulations were run across all possible pellet position combinations, having each pellet assigned with values of permeability extracted from the experimental data. This allows to quantify the impact of the considering various non-uniform distributions of permeability on the dynamics of axial gas redistribution. The present work findings enhance the understanding of axial gas transport, and provide valuable insights for the integration of a model for predicting the axial gas redistribution during a LOCA scenario into the BISON fuel performance code.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling the interfacial failure and resulting mechanical properties of z-pinned additively manufactured composites

The mechanical properties of Fused Filament Fabrication (FFF) parts are limited in the build direction, albeit advantages such as design flexibility and in-house customization. Here, in this paper, a novel z-pinning approach that enhances mechanical properties in the build direction by depositing material across multiple layers within the part was investigated through numerical simulations and validated experimentally. A finite element model for z-pinned composite structures was developed by assigning fiber-orientation dependent material properties obtained using a micromechanics approach to beads and pins. The bead-to-bead and pin-to-pin adhesions within the z-pinned structures were modeled using a cohesive traction separation law. The properties of cohesive elements for carbon fiber-reinforced polylactic acid (CF-PLA) z-pinned composites were calibrated using tensile experiments. The elastic modulus and tensile strength of CF-PLA z-pinned composites in the build direction were predicted with the developed numerical model. The numerical investigation on various geometrical parameters revealed that the largest pin volume increases the stiffness and tensile strength by 40% and thus, has the greatest influence on the mechanical properties of the z-pinned composites. The effect of z-pin geometrical parameters on the mechanical properties was summarized to aid in the design of z-pinned additively manufactured composite structures.

36 MATERIALS SCIENCE↗

A Spatiotemporal-Aware Weighting Scheme for Improving Climate Model Ensemble Predictions

Multimodel ensembling has been widely used to improve climate model predictions, and the improvement strongly depends on the ensembling scheme. In this work, we propose a Bayesian neural network (BNN) ensembling method, which combines climate models within a Bayesian model averaging framework, to improve the predictive capability of model ensembles. Our proposed BNN approach calculates spatiotemporally varying model weights and biases by leveraging individual models' simulation skill, calibrates the ensemble prediction against observations by considering observation data uncertainty, and quantifies epistemic uncertainty when extrapolating to new conditions. More importantly, the BNN method provides interpretability about which climate model contributes more to the ensemble prediction at which locations and times. Thus, beyond its predictive capability, the method also brings insights and understanding of the models to guide further model and data development. In this study, we design experiments using an ensemble of CMIP6 climate model simulations to illustrate the BNN ensembling method's capability with respect to prediction accuracy, interpretability, and uncertainty quantification (UQ). We demonstrate that BNN can correctly assign larger weights to the regions and seasons where the individual model fits the observation better. Moreover, its offered interpretability is consistent with our understanding of localized climate model performance. Additionally, BNN shows an increasing uncertainty when the prediction is farther away from the period with constrained data, which appropriately reflects our trustworthiness of the models in the changing climate.

54 ENVIRONMENTAL SCIENCES↗

Planar Defect Layers Template a High-Pressure InBi Polymorph

The short- and long-range order of III–V materials under high pressure has long been the subject of debate, with advancements in structural characterization leading to significant revisions to the accepted structural models. Despite these revisions, previous high-pressure structural assignments in the In–Bi system include the site-disordered β-Sn structure type, a structure type demonstrated to be nonexistent in analogous III–V systems. While X-ray diffraction is consistent with site disordering in InBi at high pressure, cluster expansion calculations indicate that disordering requires temperatures above 3000 K. Here, we propose InBi as a model material for studying unique high-pressure planar defects due to its highly anisotropic stress-dependent properties and structure. Specifically, we identify two sets of planar defects that mimic the diffraction pattern of a site disordered β-Sn structure type and are compatible with the calculated disorder barrier. We derive these defects by symmetry relations over crystallographic transitions. Density functional theory calculations of the proposed defects suggest that these defects are stabilized by diminishing interlayer separations with pressure. Further, we find that one of the proposed defects closely resembles a bulk high-pressure phase of InBi, InBi-ϵ, and we assert that the proposed defects order upon heating, acting as a template for InBi-ϵ growth. The proposed defects and their electronic structure provide a basis for the trend of superconducting critical temperature with increasing pressure. These methods for identifying defects are generalizable to other materials with reports of site disorder at high pressure, prompting a broader search for related high-pressure defects.

36 MATERIALS SCIENCE↗

Emulation of seismic-phase traveltimes with machine learning

SUMMARY We present a machine learning (ML) method for emulating seismic-phase traveltimes that are computed using a global-scale 3-D earth model and physics-based ray tracing. Accurate traveltime predictions based on 3-D earth models are known to reduce the bias of event location estimates, increase our ability to assign phase labels to seismic detections and associate detections to events. However, practical use of 3-D models is challenged by slow computational speed and the unwieldiness of pre-computed lookup tables that are often large and have prescribed computational grids. In this work, we train a ML emulator using pre-computed traveltimes, resulting in a compact and computationally fast way to approximate traveltimes that are based on a 3-D earth model. Our model is trained using approximately 850 million P-wave traveltimes that are based on the global LLNL-G3D-JPS model, which was developed for more accurate event location. The training-set consists of traveltimes between 10 393 global seismic stations and randomly sampled event locations that provide a prescribed, distance-dependent geographic sample density for each station. Prediction accuracy is dependent on event-station distance and whether the station was included in the training set. For stations included in the training set the mean absolute deviation (MAD) of the difference between traveltimes computed using ray tracing through the 3-D model and the ML emulator for local, regional, and teleseismic distances are 0.090, 0.125 and 0.121 s, respectively. For tested station locations not included in the training set, MAD values for the three distance ranges increase to 0.173, 0.219 and 0.210 s, respectively. Empirical traveltime residuals for a global reference data are indistinguishable when ML emulation or the 3-D model is used to compute traveltimes. This result holds regardless of whether the recording station is used in ML training or not.

58 GEOSCIENCES↗