Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forward modeling”

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 685 records · Page 38

Roadmap on multiscale materials modeling

Modeling and simulation is transforming modern materials science, becoming an important tool for the discovery of new materials and material phenomena, for gaining insight into the processes that govern materials behavior, and, increasingly, for quantitative predictions that can be used as part of a design tool in full partnership with experimental synthesis and characterization. Modeling and simulation is the essential bridge from good science to good engineering, spanning from fundamental understanding of materials behavior to deliberate design of new materials technologies leveraging new properties and processes. This Roadmap presents a broad overview of the extensive impact computational modeling has had in materials science in the past few decades, and offers focused perspectives on where the path forward lies as this rapidly expanding field evolves to meet the challenges of the next few decades. The Roadmap offers perspectives on advances within disciplines as diverse as phase field methods to model mesoscale behavior and molecular dynamics methods to deduce the fundamental atomic-scale dynamical processes governing materials response, to the challenges involved in the interdisciplinary research that tackles complex materials problems where the governing phenomena span different scales of materials behavior requiring multiscale approaches. The shift from understanding fundamental materials behavior to development of quantitative approaches to explain and predict experimental observations requires advances in the methods and practice in simulations for reproducibility and reliability, and interacting with a computational ecosystem that integrates new theory development, innovative applications, and an increasingly integrated software and computational infrastructure that takes advantage of the increasingly powerful computational methods and computing hardware.

36 MATERIALS SCIENCE↗

Chlorine Isotope Separations using Thermal Diffusion

As of the close of 2023 a thermal diffusion isotope separations (TDIS) apparatus was constructed, successful shakedown testing was achieved and enrichments of isotopic concentrations relative to natural abundance 35/37 Cl were collected. Further, the model for these enrichment experiments drove the timing for sampling and other critical, extrapolated functions that are discussed herein. Our decision to move forward in 2024 with the installation of a larger set of separation tubes and associated equipment was based on the successful development of the predictive model. The Chlorine Isotopes Project Team at PNNL is prepared to claim that the installation and hence its separative power is restricted only by the spatial limitations of the laboratory, and this at present appears to be the limiting feature to first pass high enriched 37 Cl.

07 ISOTOPE AND RADIATION SOURCES↗

Measurement of prompt D + and ${D}_s^{+}$ production in p Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The production of prompt D + and $D^+_s$ mesons is studied in proton-lead collisions at a centre-of-mass energy of $\sqrt{^{S}NN}$ = 5.02 TeV. The data sample corresponding to an integrated luminosity of (1.58 ± 0.02)nb –1 is collected by the LHCb experiment at the LHC. The differential production cross-sections are measured using D + and $D^+_s$ candidates with transverse momentum in the range of 0 < p T < 14 GeV/c and rapidities in the ranges of 1.5 < y* < 4.0 and –5.0 < y* < –2.5 in the nucleon-nucleon centre-of-mass system. For both particles, the nuclear modification factor and the forward-backward production ratio are determined. These results are compared with theoretical models that include initial-state nuclear effects. In addition, measurements of the cross-section ratios between D + , $D^+_s$ and D 0 mesons are presented, providing a baseline for studying the charm hadronization in lead-lead collisions at LHC energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Core-to-core X-ray emission spectra from Wannier based multiplet ligand field theory

Recent advances using Density Functional Theory (DFT) to augment Multiplet Ligand Field Theory (MLFT) have led to ab-initio calculations of many formerly empirical parameters. This development makes MLFT more predictive instead of interpretive, thus improving its value for studies of highly correlated 3 d , 4 d , and ƒ-electron systems. Synchrotron time is always at a premium, and tools that provide predictive capabilities have clear value when it comes to planning studies. Here, in this work, we develop a DFT + MLFT based approach for core-to-core Kα x-ray emission spectra (XES) and evaluate its performance for a range of transition metal systems. We find good agreement between theory and experiment, as well as the ability to capture key spectral trends related to spin and oxidation state. We also discuss limitations of the model in the context of the remaining free parameters and suggest directions forward.

Ab initio theoretical treatment↗

AI-accelerated protein-ligand docking for SARS-CoV-2 is 100-fold faster with no significant change in detection

Protein-ligand docking is a computational method for identifying drug leads. The method is capable of narrowing a vast library of compounds down to a tractable size for downstream simulation or experimental testing and is widely used in drug discovery. While there has been progress in accelerating scoring of compounds with artificial intelligence, few works have bridged these successes back to the virtual screening community in terms of utility and forward-looking development. We demonstrate the power of high-speed ML models by scoring 1 billion molecules in under a day (50 k predictions per GPU seconds). We showcase a workflow for docking utilizing surrogate AI-based models as a pre-filter to a standard docking workflow. Our workflow is ten times faster at screening a library of compounds than the standard technique, with an error rate less than 0.01% of detecting the underlying best scoring 0.1% of compounds. Our analysis of the speedup explains that another order of magnitude speedup must come from model accuracy rather than computing speed. In order to drive another order of magnitude of acceleration, we share a benchmark dataset consisting of 200 million 3D complex structures and 2D structure scores across a consistent set of 13 million “in-stock” molecules over 15 receptors, or binding sites, across the SARS-CoV-2 proteome. We believe this is strong evidence for the community to begin focusing on improving the accuracy of surrogate models to improve the ability to screen massive compound libraries 100 × or even 1000 × faster than current techniques and reduce missing top hits. The technique outlined aims to be a fast drop-in replacement for docking for screening billion-scale molecular libraries.

59 BASIC BIOLOGICAL SCIENCES↗

Publishing Environmental Assessment and Management Science: Crossing the Hurdles

Benefits accrue to scientists, resource managers, companies, and policymakers when environmental scientists publish in peer-reviewed journals. However, environmental scientists and practitioners face challenges, including the sometimes low value placed on journal articles, institutional vested interests in outcomes, and the changing priorities of employers and project sponsors. Confidentiality agreements can also lead scientists to assume publication is not an option. Case studies may be viewed by potential authors as too routine for peer-reviewed journals. On the basis of 30 years of experience, we suggest that publishing hurdles can be overcome and that environmental scientists have a range of options. The topics of manuscripts can include not only results from case studies and perspectives based on them but also byproducts of assessments, including definitions, plans, monitoring methods and models, and decision frameworks. Environmental scientists have unique opportunities to move science forward with their practical knowledge if they can move across the institutional, logistical, data-related, and content-related hurdles.

54 ENVIRONMENTAL SCIENCES↗

Decorrelation of participant and spectator angular momenta in heavy-ion collisions

High-energy heavy-ion collisions contain enormous angular momentum, |$\vec{𝐽}$|, which is 𝑂⁡(10 3 –10 6 ⁢ℏ) in the range of collision energy, $\sqrt{s_{NN}}$, spanned experimentally by the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). Here, a fraction of $\vec{𝐽}$ is transferred to the overlapping collision region, which is indispensable for measuring observables such as vorticity-driven hadron spin alignment with $\hat{𝐽}$. Experiments estimate the orientation of $\hat{𝐽}$ of the participant nucleons within the collision overlap region, $\hat{𝐽}$ part , by using that of the forward- and backward-going spectating nucleons $\hat{𝐽}$ spec . Using two models, we study the decorrelation between $\hat{𝐽}$ part and $\hat{𝐽}$ spec , driven both by angular-momentum conservation and event-by-event fluctuations, as well as by the decorrelation between the orientation of the elliptic overlap region and the $\hat{𝐽}$ part . $\sqrt{s_{NN}}$-dependent decorrelation is observed in both of these cases and is large enough to be an important corrective factor used when experimentally observing phenomena driven by $\vec{𝐽}$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Workbooks for Cambium 2024 Data

These workbooks contain a subset of cost and emissions data from the 2024 Cambium datasets, with levelization calculations to assist users in translating Cambium’s year-over-year values to a representative value for a user-specified project timeline. These workbooks provide modeled data for 18 GEA regions covering the contiguous United States, projected forward through 2050. Mappings of these regions to ZIP codes and counties is given within this workbook in the corresponding tabs. For the full Cambium 2024 data sets, see the Cambium 2024 project on NREL's Scenario Viewer. For more details on input assumptions and methodology see the associated report: Cambium 2024 Scenario Descriptions and Documentation. Users are advised to review section 4 of the report, which discusses limitations and caveats of the data. This data is planned to be updated annually. Information on the latest versions can be found on the NREL energy analysis page on Cambium.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

60 years of science in ICF: from conception to scientific breakeven on the National Ignition Facility

The recent achievements of a burning plasma, fusion ignition, and scientific energy gain with deuterium-tritium (DT) fuel at Lawrence Livermore National Laboratory’s National Ignition Facility (NIF) represents a major milestone in the development of inertial confinement fusion (ICF) and all of fusion research. In these experiments, fuel pressures well in excess of hundreds of GBars were achieved in the compressed fuel, and robust alpha heating of the fuel, far in excess of the energy provided by the implosion, were demonstrated for the first time. These achievements occurred 60 years after the inception of ICF and the first laser demonstration, and were made possible by more than five decades of research at laser facilities around the world. Advances in laser technology both in wavelength and precision, motivated by improved understanding of laser-plasma interaction physics and the demands of targets; improvements in target fabrication inspired by the need to control and minimize hydrodynamic instabilities in the implosion; and multi-dimensional simulations and diagnostics have been critical to this achievement. This paper will summarize the scientific and technical advances, the surprises, and the challenges that had to be overcome to achieve these goals.

fusion↗

High concentrations of ice: Investigations using polarimetric radar observations combined with in situ measurements and cloud modeling (Final Report)

This DOE-funded joint project had the over-arching goal of understanding the reasons for observed size distributions of ice particles in cold clouds. Its approach involves the use of cloud models and field observations by radar and aircraft of real storms. The project addresses three major research objectives: (1) Utilize a novel polarimetric radar technique to retrieve size distributions and amounts of ice from radar data collected during the previous DOE ARM field campaigns; (2) Evaluate the results of ice microphysical retrievals using available in situ aircraft measurements and other remote sensors (e.g., vertically pointing cloud radars and wind profilers); (3) Improve the treatment of microphysical processes leading to cloud glaciation and ice multiplication in numerical cloud models. The study was performed by three research organizations: University of Oklahoma, The Hebrew University of Jerusalem, Israel, and Lund University, Sweden. The Israel and Swedish partners were funded by the sub-awards from the University of Oklahoma. Cloud modeling studies have been performed by the research teams at The Hebrew University of Jerusalem (HUJ) and Lund University to identify the origins of high concentrations of cloud ice in areas of high ice water content (HIWC). The Hebrew University Cloud Model (HUCM) with full spectral bin microphysics and the Lund University aerosol-cloud (AC) model with a hybrid bin / bulk microphysics scheme complementing HUCM were utilized for simulations. Both research teams had particular focus on secondary ice production (SIP) as one of the possible sources of enhanced ice concentration. The HUJ group suggested a novel concept of ice multiplication during droplet freezing. It is assumed that splintering and droplet fragmentation during droplet freezing takes place because of dendritic growth within a supercooled drop. The resulting simulations of SIP generated small ice in concentrations exceeding hundreds per liter similar to what was observed in the HIWC regions of the tropical storms. The Lund team explored the SIP mechanisms such as breakup of ice particles due to ice-ice collisions and ice sublimation that are expected to dominate the continental storms. They also quantified the impact of homogeneous nucleation of cloud droplets on the total number concentration of ice at very low temperatures near the tops of the clouds. Additionally, the Lund AC model is able to simulate the effect of aerosols of various types (including biological) on the cloud life cycle and the corresponding ice production. The HUCM / AC model was used to simulate one of the “golden” cases of the DOE MC3E campaign on 20 May 2011. The model output was converted into the fields of polarimetric radar variables using the polarimetric radar forward operator developed at the University of Oklahoma and compared with radar observations and in situ microphysical measurements onboard research aircraft. It was demonstrated that specific differential phase KDP is the best radar parameter to identify the HIWC areas and quantify the corresponding ice parameters. For the first time, the shape of the vertical profile of KDP was realistically reproduced by the cloud model with the KDP maximum in the dendritic growth layer (DGL) centered at the -15°C isotherm. The University of Oklahoma team has developed a methodology for polarimetric radar retrievals of such microphysical parameters of ice as ice water content (IWC), mean volume diameter (Dm), and total number concentration (Nt) of ice particles. These retrievals have been validated using in situ aircraft measurements during 6 field campaigns and proved to be quite robust and reliable. This allowed to build the first climatology of the vertical profiles of polarimetric radar variables and retrieved microphysical parameters for the three types of weather systems: continental MCSs, maritime MCSs, and tropical cyclones / hurricanes (Hu and Ryzhkov 2022). The data were collected by a multitude of the WSR-88D radars in 13 continental and 10 maritime MCSs and 11 landfalling hurricanes. The HIWC areas were identified within the examined storms and the corresponding “HIWC statistics” was compared with the “background” one without HIWC. An overarching conclusion of the study is that maritime tropical storms (MCSs and hurricanes) are characterized by smaller size ice in higher concentration compared to the continental MCSs. High ice water content in the HIWC areas is primarily caused by a strong jump in a number concentration of ice particles rather than the increase of their size compared to the “background” environment. This may point to the homogeneous nucleation of excessive amounts of supercooled droplets and / or secondary ice production as the possible origins of HIWC. Such a climatology provides a good observational reference for the modelers to evaluate the performance of their models. As an example, the in-depth analysis of the 20 May 2011 MC3E case shows that the advanced cloud models developed in the course of this study still tend to underestimate the number concentration of ice in the HIWC areas although they succeed in reproducing realistically looking vertical profiles of IWC and Nt. The results of the project research are summarized in 13 journal papers.

54 ENVIRONMENTAL SCIENCES↗

Machine learning the spectral function of a hole in a quantum antiferromagnet

Understanding charge motion in a background of interacting quantum spins is a fundamental problem in quantum many-body physics. The most extensively studied model for this problem is the so-called t-t'-t''-J model, where the determination of the parameter t' in the context of cuprate superconductors is challenging. Here we present a theoretical study of the spectral functions of a mobile hole in the t-t'-t''-J model using two machine-learning techniques: K-nearest neighbor regression (KNN) and a feed-forward neural network (FFNN). We employ the self-consistent Born approximation to generate a dataset of about 1.3 x 10 5 spectral functions. Here we show that, for the forward problem, both methods allow for the accurate and efficient prediction of spectral functions, allowing, e.g., rapid searches through parameter space. Furthermore, we find that for the inverse problem (inferring Hamiltonian parameters from spectra), the FFNN can, but the KNN cannot, accurately predict the model parameters using merely the density of states. Our results suggest that it may be possible to use deep-learning methods to predict materials parameters from experimentally measured spectral functions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Parametric Study of SANSMIC Input and Resulting Impact on Predicted Cavern Geometry and Leaching Efficiency

The Sandia Solution Mining Code (SANSMIC) has been used for many years to examine the development of salt cavern geometry, both in a confirmatory manner with comparisons made to real-world sonar data and in a predictive manner when updated sonar data are not available. SANSMIC models require some modeling choices in order to incorporate real-world data. Key modeling choices include the vertical resolution of cavern geometry to implement, as well as how to incorporate daily raw water injection data into the SANSMIC model. This report documents five studies that address the impact of the modeling choices on the predicted cavern geometries and calculated leaching efficiencies. In most cases, hypothetical cylindrical initial cavern geometries are used to provide a common baseline against which to test the systematic variation of input variables including cavern radius, oil-brine-interface (OBI) depth, vertical cell size, raw water injection rate, raw water injection duration, workover time, and number of leaching stages. The use of smaller cell sizes is recommended moving forward to provide a better one-to-one relationship between sonar data and the modeled cavern. A new methodology for incorporating raw water injection data is also recommended, in order to more closely model real-world injection and workover times. Overall, the systematic studies performed here have increased our confidence in previous SANSMIC model results, as well future use of the code for predicting leaching effects on cavern geometries. Some minor changes to modeling choices are recommended, which can easily be applied with the version of SANSMIC currently under development.

42 ENGINEERING↗

FACET: a new long-lived particle detector in the very forward region of the CMS experiment

We describe a proposal to add a set of very forward detectors to the CMS experiment for the high-luminosity era of the Large Hadron Collider to search for beyond the standard model long-lived particles, such as dark photons, heavy neutral leptons, axion-like particles, and dark Higgs bosons. The proposed subsystem is called FACET for Forward-Aperture CMS ExTension, and will be sensitive to any particles that can penetrate at least 50 m of magnetized iron and decay in an 18 m long, 1 m diameter vacuum pipe. The decay products will be measured in detectors using identical technology to the planned CMS Phase-2 upgrade.

Particle and Resonance Production↗

Instantiation of HCML Demonstrating Bayesian Predictive Modeling for Attentional Control

The research team developed models of Attentional Control (AC) that are unique to existing modeling approaches in the literature. The goal was to enable the research team to (1) make predictions about AC and human performance in real-world scenarios and (2) to make predictions about individual characteristics based on human data. First, the team developed a proof-of-concept approach for representing an experimental design and human subjects data in a Bayesian model, then demonstrated an ability to draw inferences about conditions of interest relevant to real-world scenarios. Ultimately, this effort was successful, and we were able to make reasonable (meaning supported by behavioral data) inferences about conditions of interest to develop a risk model for AC (where risk is defined as a mismatch between AC and attentional demand). The team additionally defined a path forward for a human-constrained machine learning (HCML) approach to make predictions about an individual's state based on performance data. The effort represents a successful first step in both modeling efforts and serves as a basis for future work activities. Numerous opportunities for future work have been defined.

97 MATHEMATICS AND COMPUTING↗

Dark Energy Survey Year 6 Results: Point-Spread Function Modeling

We present the point-spread function (PSF) modeling for weak lensing shear measurement using the full six years of the Dark Energy Survey (DES Y6) data. We review the PSF estimation procedure using the PIFF (PSFs In the Full FOV) software package and describe the key improvements made to PIFF and modeling diagnostics since the DES year three (Y3) analysis: (i) use of external Gaia and infrared photometry catalogs to ensure higher purity of the stellar sample used for model fitting, (ii) addition of color-dependent PSF modeling, the first for any weak lensing analysis, and (iii) inclusion of model diagnostics inspecting fourth-order moments, which can bias weak lensing measurements to a similar degree as second-order modeling errors. Through a comprehensive set of diagnostic tests, we demonstrate the improved accuracy of the Y6 models evident in significantly smaller systematic errors than those of the Y3 analysis, in which all g band data were excluded due to insufficiently accurate PSF models. For the Y6 weak lensing analysis, we include g band photometry data in addition to the riz bands, providing a fourth band for photometric redshift estimation. Looking forward to the next generation of wide-field surveys, we describe several ongoing improvements to PIFF, which will be the default PSF modeling software for weak lensing analyses for the Vera C. Rubin Observatory’s Legacy Survey of Space and Time.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Application of unsupervised deep learning to image segmentation and in-situ contact angle measurements in a CO 2 -water-rock system

Rock surface wettability is a critical property that regulates multiphase flows in porous media, which can be quantified using the surface contact angle (CA). X-ray micro-computed tomography (μCT) provides an effective approach to in-situ measurements of surface CAs. However, the CA measurement accuracy depends significantly on the quality of CT image segmentation, which is the clustering of CT pixels into separate phases. Inspired by this, we developed a deep learning (DL)-based CA measurement workflow. Motivated by the recent tremendous progress in unsupervised learning techniques and aiming to avoid expensive manual data annotations, an unsupervised DL pipeline for CT image segmentation was proposed and implemented, which includes unsupervised model training and post-processing. The unsupervised model training was driven by a novel loss function constrained with feature similarity and spatial continuity and implemented by iterative forward and backward paths; the former clustered the pixel-wise feature vectors extracted by convolution neural networks, whereas the latter updated the parameters using gradient descent. An over-segmentation strategy was adopted for model training. The post-processing steps based on agglomerative hierarchical clustering (AHC) were implemented to further merge the over-segmented model output to the desired cluster number, which is intended to improve the efficiency of image segmentation. The developed unsupervised DL pipeline was compared with other commonly-used image segmentation methods using pixel-wise and physics-based evaluation metrics on a synthetic raw-image dataset, which had a known ground truth. The unsupervised DL pipeline showed the best performance. Next, the segmented images were input to an automatic CA measurement tool, and the results were validated by comparisons with manual measurements. The CA values from the manual and automatic measurements showed similar distributions and statistical properties. The automatic measurement demonstrated a wider spectrum because of the much larger number of measurement data points. The primary novelty of the unsupervised DL pipeline developed in this study lies in the novel loss function and the over-segmentation strategy associated with AHC post-processing. Finally, the workflow has been proven an efficient tool for pore-scale wettability characterization, which has a wide range of applications in fundamental studies of multiphase flows in natural porous media, which have critical implications to geological carbon sequestration, hydrocarbon energy recovery, and contaminant transport in groundwater.

42 ENGINEERING↗

Markovian transport processes in a heterogeneous, variably saturated watershed: A multi-domain spatial Markov model

Transport processes in watersheds remain challenging to describe accurately with upscaled models because of their intricate complexities at multiple scales that can lead to multi-modal, non-Fickian breakthrough curves. One of the recent advances in solute transport modeling has been the Lagrangian spatial Markov model (SMM), which describes transport using velocity increments and models the correlated transitions between them. The method has been applied previously to saturated porous media and fractured media, but these successful applications suggest it may be useful for describing watershed scale processes as well. An existing, 3-D, heterogeneous variably saturated model of a headwaters catchment and classical random walk particle tracking were used to generate a forward simulation and the resulting transport was analyzed with a spatial Markov framework. Existing SMM methods were used to analyze the transport paths through the complex, unconfined domain and we found that the vadose zone plays an important role in the correlation structure. As such, a new kind of SMM is proposed, termed a multi-domain SMM (MD-SMM), that uses separate correlation models for auto-transitions in the vadose zone and the saturated zone, and for the cross-transitions between them. The MD-SMM improves representation of transport relative to the single domain model and can naturally delineate times when a particle is in the vadose zone or the saturated zone, adding an additional degree of realism to the upscaled model. Lastly, we show that the transition probability matrix for transport in the saturated portion of this unconfined system is approximately Markovian, which is an important validation of the proposed framework.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of dual-weighted residual and machine learning error estimation for projection-based reduced-order models of steady partial differential equations

Projection-based reduced-order models (pROMs) show great promise as a means to accelerate many-query applications such as forward error propagation, solving inverse problems, and design optimization. In order to deploy pROMs in the context of high-consequence decision making, accurate error estimates are required to determine the region(s) of applicability in the parameter space. The following paper considers the dual-weighted residual (DWR) error estimate for pROMs and compares it to another promising pROM error estimate, machine learned error models (MLEM). Here, we show how DWR can be applied to ROMs and then evaluate DWR on two partial differential equations (PDEs): a two-dimensional linear convection–reaction–diffusion equation, and a three-dimensional static hyper-elastic beam. It is shown that DWR is able to estimate errors for pROMs extrapolating outside of their training set while MLEM is best suited for pROMs used to interpolate within the pROM training set.

42 ENGINEERING↗