Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer experiments”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Defect-induced phonon-resonant scattering and its influence on thermal transport of irradiated thorium-dioxide

Thermal transport in proton irradiated thorium-dioxide (ThO 2 ) is investigated. Using a combination of experiments and first-principles computational framework, the role of lattice defects on thermal conductivity is analyzed. A resonant-phonon scattering mechanism beyond the traditionally considered Rayleigh scattering is found to significantly influence low-temperature thermal transport in the presence of irradiation-induced point defects. The existence of localized phonon modes associated with irradiation-induced defects is suggested by the inability of the first-principles based thermal conductivity model—which considers only three-phonon interactions and phonon-defects scattering using the Tamura formalism—to predict the experimental results, unless a resonant scattering mechanism is included. The emergence of additional peaks in the Raman spectra in the proximity of phonon-resonant frequency provides further evidence for the existence of localized modes. Coupled with a microstructure evolution model, this analysis enables more accurate analysis for contrasting the contributions of different phonon scattering mechanisms across all irradiation doses and temperatures.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Interface Tracking Investigation of Geometric Effects on the Bubbly Flow in PWR Subchannels

We report that absorbing heat from the fuel rod surface, water as coolant can undergo subcooled boiling within a pressurized water reactor (PWR) fuel rod bundle. Because of the buoyancy effect, the vapor bubbles generated will then rise along and interact with the subchannel geometries. Reliable prediction of bubble behavior is of immense importance to ensure safe and stable reactor operation. However, given a complex engineering system like a nuclear reactor, it is very challenging (if not impossible) to conduct high-resolution measurements to study bubbly flows under reactor operation conditions. The lack of a fundamental two-phase-flow database is hindering the development of accurate two-phase-flow models required in more advanced reactor designs. In response to this challenge, first-principles–based numerical simulations are emerging as an attractive alternative to produce a complementary data source along with experiments. Leveraged by the unprecedented computing power offered by state-of-the-art supercomputers, direct numerical simulation (DNS), coupled with interface tracking methods, is becoming a practical tool to investigate some of the most challenging engineering flow problems. In the presented research, turbulent bubbly flow is simulated via DNS in single PWR subchannel geometries with auxiliary structures (e.g., supporting spacer grid and mixing vanes). The geometric effects these structures exert on the bubbly flow are studied with both a conventional time-averaging approach and a novel dynamic bubble tracking method. Lastly, the new insights obtained will help inform better two-phase models that can contribute to safer and more efficient nuclear reactor systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

Engineering the CernVM-Filesystem as a High Bandwidth Distributed Filesystem for Auxiliary Physics Data

A common use pattern in the computing models of particle physics experiments is running many distributed applications that read from a shared set of data files. We refer to this data is auxiliary data, to distinguish it from (a) event data from the detector (which tends to be different for every job), and (b) conditions data about the detector (which tends to be the same for each job in a batch of jobs). Relatively speaking, conditions data also tends to be relatively small per job where both event data and auxiliary data are larger per job. Unlike event data, auxiliary data comes from a limited working set of shared files. Since there is spatial locality of the auxiliary data access, the use case appears to be identical to that of the CernVM- Filesystem (CVMFS). However, we show that distributing auxiliary data through CVMFS causes the existing CVMFS infrastructure to perform poorly. We utilize a CVMFS client feature called 'alien cache' to cache data on existing local high-bandwidth data servers that were engineered for storing event data. This cache is shared between the worker nodes at a site and replaces caching CVMFS files on both the worker node local disks and on the site's local squids. We have tested this alien cache with the dCache NFSv4.1 interface, Lustre, and the Hadoop Distributed File System (HDFS) FUSE interface, and measured performance. In addition, we use high-bandwidth data servers at central sites to perform the CVMFS Stratum 1 function instead of the low-bandwidth web servers deployed for the CVMFS software distribution function. We have tested this using the dCache HTTP interface. As a result, we have a design for an end-to-end high-bandwidth distributed caching read-only filesystem, using existing client software already widely deployed to grid worker nodes and existing file servers already widely installed at grid sites. Files are published in a central place and are soon available on demand throughout the grid and cached locally on the site with a convenient POSIX interface. This paper discusses the details of the architecture and reports performance measurements.

Dykstra, D.↗
DOE OSTI2023

Modeling the sequential dissociative double ionization of O 2 by ultrashort intense infrared laser pulses

A density matrix approach for sequential double ionization (DM-SDI) of molecules has been developed recently and was applied to the N 2 molecule. In this article, we extended the DM-SDI model to O 2 , which is a more complicated system to model than N 2 , due to its electronic structures and spin-orbit and laser couplings in the manifold of doubly charged states. We obtained a good agreement on the kinetic energy release spectrum of O + + O + from previous experiments. Thanks to the low computational cost of the model, we explored the mechanism behind the ionization and dissociation dynamics as well as the effects of lasers on the spectrum. Furthermore, this work will pave the way to model sequential dissociative double ionization of larger molecules and to probe molecular dynamics by measuring kinetic energy release spectra from this process.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Recovering noise-free quantum observables

We introduce a technique for recovering noise-free observables in noisy quantum systems by combining the results of many slightly different experiments. Our approach is applicable to a variety of quantum systems, but we illustrate it with applications to quantum information and quantum sensing. The approach corresponds to repeating the same quantum evolution many times with known variations on the underlying systems' error properties, e.g., the spontaneous emission and dephasing times $T_1$ and $T_2$. As opposed to standard quantum error correction methods, which have an overhead in the number of qubits (many physical qubits must be added for each logical qubit), our method has only an overhead in the number of evaluations, allowing the overhead to, in principle, be hidden via parallelization. We show that the effective spontaneous emission $T_1$ and dephasing $T_2$ times can be increased using this method in both simulation and experiments on an actual quantum computer. We also show how to correct more complicated entangled states and how Ramsey fringes relevant to quantum sensing can be significantly extended in time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Electronic structure of intrinsic defects in c -gallium nitride: Density functional theory study without the jellium approximation

Here, we report the first nonjellium, systematic, density functional theory (DFT) study of intrinsic and extrinsic defects and defect levels in zinc-blende (cubic) gallium nitride. We use the local moment counter charge (LMCC) method, the standard Perdew-Becke-Ernzerhoff (PBE) exchange-correlation potential, and two pseudopotentials, where the Ga 3$\textit{d}$ orbitals are either in the core ($d^0$) or explicitly in the valence set ($d^{10}$). We studied 64, 216, 512, and 1000 atom supercells, and demonstrated convergence to the infinite limit, crucial for delineating deep from shallow states near band edges, and for demonstrating the elimination of finite cell-size errors. Contrary to common claims, we find that exact exchange is not required to obtain defect levels across the experimental band gap. As was true in silicon, silicon carbide, and gallium arsenide, the extremal LMCC defect levels of the aggregate of defects yield an effective LMCC defect band gap that is within 10% of the experimental gap (3.3 eV) for both pseudopotentials. We demonstrate that the gallium vacancy is more complicated than previously reported. There is dramatic metastability–a nearest-neighbor nitrogen atom shifts into the gallium site, forming an antisite, nitrogen vacancy pair, which is more stable than the simple vacancy for positive charge states. Our assessment of the $d^0$ and $d^{10}$ pseudopotentials yields minimal differences in defect structures and defect levels. The better agreement of the $d^0$ lattice constant with experiment suggests that the more computationally economical $d^0$ pseudopotentials are sufficient to achieve the fidelity possible within the physical accuracy of DFT, and thereby enable calculations in larger supercells necessary to demonstrate convergence with respect to finite size supercell errors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Examples of symmetry-preserving truncations in tensor field theory

We consider the tensor formulation of the nonlinear O(2) sigma model and its gauged version (the compact Abelian Higgs model), on a D-dimensional cubic lattice, and show that tensorial truncations are compatible with the general identities derived from the symmetries of these models. This means that the universal properties of these models can be reproduced with highly simplified formulations desirable for implementations with quantum computers or for quantum simulations experiments. We discuss the extensions to global non-Abelian symmetries, discrete symmetries and pure gauge Abelian models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Solving key challenges in collider physics with foundation models

Foundation models are neural networks that are capable of simultaneously solving many problems. Large language foundation models like ChatGPT have revolutionized many aspects of daily life, but their impact for science is not yet clear. In this paper, we use a new foundation model for hadronic jets to solve three key challenges in collider physics. In particular, we show how experiments can (1) save significant computing power when developing reconstruction algorithms, (2) perform a complete uncertainty quantification for high-dimensional measurements, and (3) search for new physics with model agnostic methods using low-level inputs. In each case, there are significant computational or methodological challenges with current methods that limit the science potential of deep learning algorithms. By solving each problem, we take jet foundation models beyond proof-of-principle studies and into the toolkit of practitioners.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Optimally tuned starting point for single-shot GW calculations of solids

The dependence of ab initio many-body perturbation theory within the $\textit{GW}$ approximation on the eigensystem used in calculating quasiparticle corrections limits this method's predictive power. Here, we investigate the accuracy of the recently developed Wannier-localized optimally tuned screened range-separated hybrid (WOT-SRSH) functional as a generalized Kohn-Sham starting point for single-shot $\textit{GW}$ ($G_0W_0$) calculations for a range of semiconductors and insulators. Comparison to calculations based on well-established functionals, namely, PBE, PBE0, and HSE, as well as to self-consistent GW schemes and to experiment, shows that band gaps computed via $G_0W_0$@WOT-SRSH have a level of precision and accuracy that is comparable to that of more advanced methods such as quasiparticle self-consistent $\textit{GW}$ and eigenvalue self-consistent $\textit{GW}$. In this work, we also find that $G_0W_0$@WOT-SRSH improves the description of states deeper in the valence band manifold. Finally, we show that $G_0W_0$@WOT-SRSH significantly reduces the sensitivity of computed band gaps to ambiguities in the underlying WOT-SRSH tuning procedure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Resonant neutron scattering lengths

Unlike most of the periodic table, many rare-earth elements display considerable resonant scattering for thermal neutrons. Although this property is accompanied by strong neutron absorption, modern high-intensity neutron sources make diffraction experiments possible with these elements. Computation of scattering intensities is accomplished by fitting the variation in resonant scattering lengths ( b 0 , b ′ and b ′′) to a semi-empirical Breit–Wigner formalism, which can be evaluated over the range of neutron energies useful for diffraction, typically E = 10–600 meV; λ = 0.4–2.8 Å (with good extrapolation to longer wavelengths).

Von Dreele, Robert B. (ORCID:0000000253390634)↗
DOE OSTI2016

Kernels for scalable data analysis in science: Towards an architecture-portable future

In this paper, we pose and address some of the unique challenges in the analysis of scientific Big Data on supercomputing platforms. Our approach identifies, implements and scales numerical kernels that are critical to the instantiation of theory-inspired analytic workflows on modern computing architectures. We present the benefits of scalable kernels towards constructing algorithms such as principal component analysis and non-negative matrix factorization on an image-analysis use case at the Oak Ridge Leadership Computing Facility (OLCF). Based on experience with the use-case, we conclude that piecing scalable analytic kernels into user-defined analytic workflows are a flexible, modular and agile way to enable architecture-portable productivity for the data-intensive sciences.

Sukumar, Sreenivas R.↗
DOE OSTI2025

Chemical Reactor Network Modeling of Ammonia Rich-Quench-Lean Combustion Using a Partially Stirred Reactor Approach

Ammonia is a promising alternative to hydrogen with high energy density and favorable storage and transport characteristics. However, low flammability and a propensity for high nitrogen oxide (NO x ) emissions make direct utilization challenging. Recently, two-stage rich-quench-lean (RQL) combustion strategies have shown promise in achieving low NO x emissions with ammonia. In this approach, the rich stage serves to oxidize a portion of the fuel while thermally decomposing as much of the remaining ammonia as possible, generating hydrogen. In the second (lean) stage, air is rapidly introduced, burning out the hydrogen and residual ammonia. Two-stage RQL combustion of ammonia has been investigated in the open literature both experimentally and numerically. In general, idealized chemical reactor network (CRN) models predict NO x concentrations below those of 2D/3D computational fluid dynamics models and experiments. The primary drivers of these discrepancies may be largely attributed to finite rate mixing nonadiabatic operation. The typical CRN model is comprised of a perfectly-stirred-reactor (PSR), followed by a plug-flow-reactor (PFR), meant to represent the flame, and postflame zones, respectively. In the two-stage RQL approach two PSR-PFR networks are arranged sequentially, corresponding to the rich and lean stages, with secondary air injection in between. In the authors' past work, this arrangement has demonstrated the significant sensitivity of exit NO x to the rich stage equivalence ratio, while the amount of secondary air injection was shown to be less critical. In this paper, the CRN model is extended to (1) include the impacts of heat loss and (2) utilize a partially-stirred-reactor (PaSR) approach to study the impacts of mixing on emissions performance. Varying amounts of heat loss are applied to the rich relaxation zone to understand emissions performance and changes to optimization of equivalence ratio and residence time. Premixed and nonpremixed configurations are considered in the rich stage PaSR, with varying degrees of mixing intensity to study the interaction between mixing, transport, and kinetic timescales. Critically, the impact of mixing between hot products and secondary air injection is studied to understand practical injector needs. Results show unburnt ammonia leaving the rich stage as a primary contributor to NO x emissions – driven both by increased heat loss and reduced mixing rates. Furthermore, heat losses have been shown to create conditions that are conducive to increased N 2 O formation in the lean stage. In conclusion, the results of this study will be considered in the context of developing optimized two-stage RQL combustors for ammonia.

Combustion↗
DOE OSTI2016

Assessing task performance in FFDM, DBT, and synthetic mammography using uniform and anthropomorphic physical phantoms

Purpose The purpose of this study is to quantify the differences in detectability between full field digital mammography (FFDM), digital breast tomosynthesis (DBT), and synthetic mammography (SM) for challenging, low contrast signals, in the context of both a uniform and an anthropomorphic, textured phantom. Methods Images of the phantoms were acquired using a Hologic Selenia Dimensions system. Images were taken at 50%, 100%, and 200% of the dose delivered under automatic exposure control (AEC). Low‐contrast disks, created using an inkjet printer with iodine‐doped ink, were inserted into the phantom. The disks varied in diameter from 210 to 630 μ m, and in local contrast from 1.1% to 2.8% in regular increments. Human observers located the disks in a 4 alternative forced choice experiment. Proportion correct (PC) was computed as the number of correct localizations out of the total number of tries. Results Overall, scores from FFDM and DBT were consistently greater than scores from SM. At an exposure corresponding to the AEC setting, mean PC scores for the largest disks with the uniform phantom were 0.80 for FFDM, 0.83 for DBT, and 0.66 for SM, with the same rank ordering at other doses. Scores were similar but lower for the nonuniform background. At an exposure twice the AEC setting, however, the difference between uniform and nonuniform scores was most pronounced for DBT alone. Differences between scores for FFDM and SM were statistically significant, while those between FFDM and DBT were not. Scores were used to compute the minimum contrast level needed to reach 62.5% detection rate. The minimum contrast for SM was 36%–81% higher compared to FFDM or DBT, in either background. Conclusions This study shows that an anthropomorphic phantom and lesions inserts may be used to conduct a reader study. Detectability was significantly lower for synthetic mammography than for FFDM or DBT, for all conditions. Additionally, observer performance was consistently lower for the anthropomorphic phantom, indicating the greater challenge due to anatomical background. Because of this, it may be important to use realistic phantoms in observer studies in order to draw conclusions that are more clinically relevant.

Ikejimba, Lynda C.↗
DOE OSTI2015

Tunneling spectroscopy of a phosphorus impurity atom on the Ge(111)-(2 × 1) surface

We numerically model the Ge(111)-(2 × 1) surface electronic properties in the vicinity of a P donor impurity atom located near the surface. We find a notable increase in the surface local density of states (LDOS) around the surface dopant near the bottom of the empty surface state band π*, which we call a split state due to its limited spatial extent and energetic position inside the band gap. We show that despite the well-established bulk donor impurity energy level position at the very bottom of the conduction band, a surface donor impurity on the Ge(111)-(2 × 1) surface might produce an energy level below the Fermi energy, depending on the impurity atom local environment. It is demonstrated that the impurity located in subsurface atomic layers is visible in a scanning tunneling microscope (STM) experiment on the Ge(111)-(2 × 1) surface. The quasi-1D character of the impurity image, observed in STM experiments, is confirmed by our computer simulations with a note that a few π-bonded dimer rows may be affected by the presence of the impurity atom. We elaborate a model that allows classifying atoms on the experimental low-temperature STM image. We show the presence of spatial oscillations of the LDOS by the density-functional theory method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

IRNet: A General Purpose Deep Risidual Regression Framework for Materials Discovery

Materials discovery is crucial for making scientific advances in many domains. Collections of data from experiments and first-principle computations have spurred interest in applying machine learning methods to create predictive models capable of mapping from composition and crystal structures to materials properties. Generally, these are regression problems with the input being a 1D vector composed of numerical attributes representing the material composition and/or crystal structure. While neural networks consisting of fully connected layers have been applied to such problems, their performance often suffers from the vanishing gradient problem when network depth is increased. Hence, predictive modeling for such tasks has been mainly limited to traditional machine learning techniques such as Random Forest. In this paper, we study and propose design principles for building deep regression networks composed of fully connected layers with numerical vectors as input. We introduce a novel deep regression network with individual residual learning, IRNet, that places shortcut connections after each layer so that each layer learns the residual mapping between its output and input. We use the problem of learning properties of inorganic materials from numerical attributes derived from material composition and/or crystal structure to compare IRNet’s performance against that of other machine learning techniques. Using multiple datasets from the Open Quantum Materials Database (OQMD) and Materials Project for training and evaluation, we show that IRNet provides significantly better prediction performance than the state-of-the-art machine learning approaches currently used by domain scientists. We also show that IRNet’s use of individual residual learning leads to better convergence during the training phase than when shortcut connections are between multi-layer stacks while maintaining the same number of parameters.

deep learning↗
DOE OSTI2017

What Do Structures Tell Us About Chemokine Receptor Function and Antagonism?

Chemokines and their cell surface G protein–coupled receptors are critical for cell migration, not only in many fundamental biological processes but also in inflammatory diseases and cancer. Recent X-ray structures of two chemokines complexed with full-length receptors provided unprecedented insight into the atomic details of chemokine recognition and receptor activation, and computational modeling informed by new experiments leverages these insights to gain understanding of many more receptor:chemokine pairs. In parallel, chemokine receptor structures with small molecules reveal the complicated and diverse structural foundations of small molecule antagonism and allostery, highlight the inherent physicochemical challenges of receptor:chemokine interfaces, and suggest novel epitopes that can be exploited to overcome these challenges. The structures and models promote unique understanding of chemokine receptor biology, including the interpretation of two decades of experimental studies, and will undoubtedly assist future drug discovery endeavors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Baseline Assessment of TREAT for Modeling and Analysis Needs

TREAT is an air-cooled, graphite moderated, thermal, heterogeneous test facility designed to evaluate reactor fuels and structural materials under conditions simulating various types of nuclear excursions and transient undercooling situations that could occur in a nuclear reactor. After 21 years in a standby mode, TREAT is being re-activated to revive transient testing capabilities. Given the time elapsed and the concurrent loss of operating experience, current generation and advanced computational methods are being applied to begin TREAT modeling and simulation prior to renewed at-power operations. Such methods have limited value in predicting the behavior of TREAT without proper validation. Hence, the U.S. DOE has developed a number of programs to support development of benchmarks for both critical and transient operations. Extensive effort has been expended at INL to collect detailed descriptions, drawings and specifications for all aspects of TREAT, and to resolve conflicting data found through this process. This report provides a collection of these data, with updated figures that are significantly more readable than historic drawings and illustrations, compositions, and dimensions based on the best available sources. This document is not nor should it be considered to be a benchmark report. Rather, it is intended to provide one-stop shopping, to the extent possible, for other work that seeks to prepare detailed, accurate models of the core and its components. Given the nature of the variety of historic documents available and the loss of institutional memory, the only completely accurate database of TREAT data is TREAT itself. Unfortunately, disassembly of TREAT for inspection, assay, and measurement is highly unlikely. Hence the data provided herein is intended serve as a best-estimate substitute.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

SciDAC GSEP: Gyrokinetic Simulation of Energetic Particle Turbulence and Transport

Energetic particle (EP) confinement is a key physics issue for burning plasma experiment ITER, the crucial next step in the quest for clean and abundant energy, since ignition relies on self-heating by energetic fusion products (α-particles). Due to the strong coupling of EP with burning thermal plasmas, plasma confinement property in the ignition regime is one of the most uncertain factors when extrapolating from existing fusion devices to the ITER tokamak. EP population in current tokamaks are mostly produced by auxiliary heating such as neutral beam injection (NBI) and radio frequency (RF) heating. Remarkable progress in developing comprehensive EP simulation codes and understanding basic EP physics has been made by two concurrent SciDAC EP projects GSEP funded by the Department of Energy (DOE) Office of Fusion Energy Science (OFES), which have successfully established gyrokinetic turbulence simulation as a necessary paradigm shift for studying the EP confinement in burning plasmas. Verification and validation have rapidly advanced through close collaborations between simulation, theory, and experiment. Furthermore, productive collaborations with computational scientists have enabled EP simulation codes to effectively utilize current petascale computers and emerging exascale computers. We review here key physics progress in the GSEP projects regarding verification and validation of gyrokinetic simulations, nonlinear EP physics, EP coupling with thermal plasmas, and reduced EP transport models. Advances in high performance computing through collaborations with computational scientists that enable these large scale electromagnetic simulations are also highlighted. These results have been widely disseminated in numerous peer-reviewed publications including many Phys. Rev. Lett. papers and many invited presentations at prominent fusion conferences such as the biennial International Atomic Energy Agency (IAEA) Fusion Energy Conference and the annual meeting of the American Physics Society, Division of Plasma Physics (APS-DPP).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗