Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CALCULATION METHODS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

On the Crystal Structure of Colloidally Prepared Metastable Ag2Se Nanocrystals

Structural polymorphism is known for many bulk materials; however, on the nanoscale metastable polymorphs tend to form more readily than in the bulk, and with more structural variety. One such metastable polymorph observed for colloidal Ag2Se nanocrystals has traditionally been referred to as the “tetragonal” phase of Ag2Se. While there are reports on the chemistry and properties of this metastable polymorph, its crystal structure, and therefore electronic structure, has yet to be determined. We report that an anti-PbCl2-like structure type (space group P21/n) accurately describes the powder X-ray diffraction and X-ray total scattering patterns of colloidal Ag2Se nanocrystals prepared by several different methods. Density functional theory (DFT) calculations indicate that the anti-PbCl2-like Ag2Se polymorph is a dynamically stable, narrow-band gap semiconductor. DFT results reveal a dense theoretical Ag2Se phase space with many low-energy polymorphs, which helps explain the large number of polymorphs reported in the literature. Analysis and calculation data are stored in the zip archive. The `ag2se-calcs.aiida.` contains the provenance of the calculations and can be imported into an AiiDA database instance. The antiPbCl2like_Ag2Se_laboratory.cif file is the Reitveld refined Ag2Se structure starting from the PbCl2 structure.

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

Experimental and theoretical study of the soot-forming tendencies of furans as potential biofuels

Sustainably-produced biofuels, including furans derived from lignocellulosic biomass, are being studied as potential next-generation fuels. Biofuels have unique structures that are absent in conventional fuels and they potentially emit fewer particulates than conventional fuels because of the presence of oxygen atoms. In this work, we focus on quantification and measurement of sooting tendency. A set of furans was selected and their sooting tendencies were characterized by Yield Sooting Index (YSI), which was obtained by measuring the line-of-sight spectral radiance (LSSR) of a coflow laminar diffusion flame doped with 1000 ppm of each test fuel. The measured YSIs are lower than traditional gasoline fuels and benzenoid aromatics (e.g., YSIbenzene = 100). We also performed numerical simulations with a detailed chemical kinetic mechanism to analyze the flame chemistry of 2-methylfuran (2MF). 4→5 H-atom migration in the furan ring causes 31.0% of the total decomposition and is the dominant fuel consumption pathway leading to soot. DFT calculations via the G4 composite method were carried out to obtain activation energies (Ea) and bond dissociation energies (BDE). The results showed that the 4→5 H-atom migration in the furan ring has the second lowest Ea and subsequently forms a 1,2-diene-substituted ketone derivative as a stable species. The ring-opening products of the 4→5 H-atom migration are propyne + a linear ketene for 2-n-alkylfurans, and 1-butyne + a methyl-substituted ketene for 2-t-butylfuran. These results can explain the trends of the measured YSIs of 2-alkylfurans.

09 BIOMASS FUELS↗
DOE OSTI2022

Verification and Validation of the Stakeholder Tool for Assessing Radioactive Transportation (START) – 22323

The U.S. Department of Energy (DOE) Office of Integrated Waste Management is planning for the eventual transportation, storage, and disposal of spent nuclear fuel (SNF) and high-level radioactive waste (HLW) from nuclear power plant and DOE sites. The Stakeholder Tool for Assessing Radioactive Transportation (START) is a web-based, geospatial decision-support tool developed for evaluating routing options and other aspects of transporting SNF and HLW, covering rail, truck, barge, and intermodal infrastructure and operations in the continental United States. The verification and validation (V&V) process is intended to independently assess START to provide confidence in the ability of START to accurately provide intended results. The V&V process checks the START tool using a variety of methods, ranging from independent hand calculations to comparison of START performance and results to those of other codes. The V&V activity was conducted independently from the START development team with opportunities to provide feedback and collaborate throughout the process. The V&V analyzed attributes of transportation routes produced by START, including route distance and both population and population density captured within buffer zones around routes. Population in the buffer zone, population density in the buffer zone, and route distance were all identified as crucial outputs of the START code and were subject to V&V tasks. Some of the improvements identified through the V&V process were standardizing the underlying population data in START, changing the projection of the population raster data, and changes to the methodology used for population density to improve its applicability for expected users. This collaboration also led to suggested improvements to some of the underlying shape file segments within START. A formal series of test routes went through the V&V process from START version 3.2.1 following the implementation of the recommended improvements to START. The V&V of the buffer zone population and distance reported in START were achieved by utilizing a custom workflow developed in QGIS. This platform was chosen because it is fully independent of START. The buffer zone population and route distances in START showed excellent agreement to independent V&V test results. Over 200 route test cases were run; in all cases, the percent difference in the population within the buffer zone was less than +/- 5%, and the majority of cases were below +/- 1%. These cases were also duplicated in ArcMap, which is Environmental Systems Research Institute’s (ESRI) GIS desktop application. The ArcMap results also show good agreement; however, they cannot be considered independent of START because START also uses ESRI software. The V&V tasks are currently addressing population density and estimated radiation dose in START as well as implementing further automation of the V&V process to allow larger suites of test routes to be run more easily and quickly by less experienced GIS users. These improvements to the V&V process will make it a valuable tool for the development team and the stakeholders so that each new iteration of START can easily be checked for quality and consistency. This paper and presentation will describe the V&V process, findings and implications, and ongoing V&V activities.

START, DOE, Transporation, dose, GIS↗
DOE OSTI2021

Detection of Outliers in LiDAR Data Acquired by Multiple Platforms over Sorghum and Maize

High-resolution point cloud data acquired with a laser scanner from any platform contain random noise and outliers. Therefore, outlier detection in LiDAR data is often necessary prior to analysis. Applications in agriculture are particularly challenging, as there is typically no prior knowledge of the statistical distribution of points, plant complexity, and local point densities, which are crop-dependent. The goals of this study were first to investigate approaches to minimize the impact of outliers on LiDAR acquired over agricultural row crops, and specifically for sorghum and maize breeding experiments, by an unmanned aerial vehicle (UAV) and a wheel-based ground platform; second, to evaluate the impact of existing outliers in the datasets on leaf area index (LAI) prediction using LiDAR data. Two methods were investigated to detect and remove the outliers from the plant datasets. The first was based on surface fitting to noisy point cloud data via normal and curvature estimation in a local neighborhood. The second utilized the PointCleanNet deep learning framework. Both methods were applied to individual plants and field-based datasets. To evaluate the method, an F-score was calculated for synthetic data in the controlled conditions, and LAI, the variable being predicted, was computed both before and after outlier removal for both scenarios. Results indicate that the deep learning method for outlier detection is more robust than the geometric approach to changes in point densities, level of noise, and shapes. The prediction of LAI was also improved for the wheel-based vehicle data based on the coefficient of determination (R2) and the root mean squared error (RMSE) of the residuals before and after the removal of outliers.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Projection-to-Born-improved subtractions at NNLO

While the current frontier in fixed-order precision for collider observables is N 3 LO, important steps are necessary to consolidate NNLO cross-section predictions with improved stability and efficiency. Slicing methods have been successfully applied to obtain NNLO and N 3 LO predictions, but have shown poor performance in the presence of fiducial cuts due to large kinematical power corrections. In this paper we implement Projection-to-Born-improved q T (P2B q T ) and jettiness (P2Bτ 0 ) subtractions for a large class of color singlet processes in MCFM. This method allows for the efficient evaluation of fiducial power corrections in any non-local subtraction scheme using a Projection-to-Born subtraction. We demonstrate the significant numerical improvements of this method based on fiducial Drell-Yan and Higgs cross-sections. Moreover, with fiducial power corrections removed via this method, the leading-logarithmic power corrections that have only been calculated without fiducial cuts can be included, further improving the calculations. For di-photon production with photon isolation, we devise a novel method in combination with P2B-improved subtractions, which we name P2B γ τ 0 , and P2B γ q T for the two subtraction schemes, respectively. This method allows the inclusion of both fiducial power corrections due to kinematic cuts on the photons and a set of isolation power corrections in the fragmentation channel where a quark may enter the isolation cone. We find significant improvements in the convergence of NNLO di-photon cross-sections with photon isolation cuts, demonstrating that it is possible to achieve a stable and efficient calculation of di-photon cross-sections using slicing methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Accelerating Thermochemical Equilibrium Calculations for Nuclear Reactor Applications

Thermochemical properties play a key role in modeling and simulation of several key phenomena in nuclear reactors. There has been an increasing interest in incorporating CALPHAD-based formulations in multiphysics simulations including for Molten Salt Reactors where knowledge of phase evolution of the salt and the chemical potentials of various elements are of utmost importance in source term analyses and redox control. However, the size of such simulations is often limited by the high computational cost of full thermodynamic equilibrium calculations. This work discusses the current efforts aimed at accelerating thermochemical equilibrium calculations for multiphysics simulations performed using the open-source finite element / finite volume code Multiphysics Object Oriented Simulation Environment (MOOSE) [1]. While several methods have been proposed for accelerating phase equilibrium calculations [2], most focus on relatively small systems and often rely on a- priori knowledge of the state-space of the system. Nuclear materials, however, are often multi-component systems owing to the evolution of composition under irradiation and an approach based on a-priori mapping of phase diagram is often not enough. This work is aimed at demonstrating an on-the-fly surrogate modeling framework that uses active learning to reduce the number of full equilibrium calculations that must be performed. By combining with efficient coupling approaches, the surrogate framework helps in reducing the computational cost of thermodynamic equilibrium informed multiphysics simulations of nuclear materials. The performance is benchmarked against full coupling with the thermochemistry library Thermochimica [3]. This work uses a machine learning based approach for constructing surrogate models to predict the stable phases in a multicomponent system. The surrogates were constructed using neural networks and Gaussian process classification. In this work, we compare the relative performance of the two methods. We also demonstrate the use of caching previous calculations by interpolating the values from nearest neighbors. References [1] Lindsay, A.D., et al. "2.0 – MOOSE: Enabling massively parallel multiphysics simulation", SoftwareX, 20 (2022): 101202. [2] Roos, W.A. and Zietsman J.H. "Accelerating complex chemical equilibrium calculations – A Review", Calphad, 77 (2022): 102380. [3] Piro, M.H.A., et al. "The thermochemistry library Thermochimica", Computational Materials Science, 67 (2013): 266-272.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Center for Tokamak Transient Simulations (RPI Unstructured Mesh Developments for FES SciDAC4 Partnerships) (Final Report)

The goal of this project was the development of unstructured mesh technologies for fusion simulation codes” for FES SciDAC partnerships and to integrate those technologies into the simulation workflows of those partnerships. Specific developments were executed in support of the following FES SciDAC4 partnerships: Partnership Center for High‐fidelity Boundary Plasma Simulation (HBPS), Center for Integrated Simulation of Fusion Relevant RF (RF‐SciDAC), Center for Plasma Surface Interactions: Predicting the Performance and Impact of Dynamic PFC Surfaces (PSI2), and Center for Tokamak Transient Simulations (CTTS). The key unstructured mesh development areas addressed in this project include (i) methods to most effectively perform PIC calculations on unstructured meshes; (ii) creating meshes for fusion systems accounting for any desired level of geometric complexity and providing physics aware mesh configurations, (iii) adapting unstructured meshes to control the discretization errors, (iv) executing unstructured mesh calculations on GPU accelerated systems, (v) supporting physics/application‐specific PIC operations including surface/wall interactions of particles, (vi) providing infrastructure tools to support the interactions of solvers with unstructured meshes, and (vii) providing advanced methods for coupling plasma physics codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Robust foil regenerator flow loss and heat transfer tests under oscillating flow condition

A regenerator works as an energy storage component in a Stirling engine. Thus, it is one of the most critical components that makes the Stirling engine highly efficient. Traditionally, the regenerators are made of wire media such as woven screens and random fibers. These wire regenerators experience considerably high flow losses through the regenerator due to the flow separation and the resulting vortices. The high flow losses associated with the flow separation and vortices result in a lower engine efficiency. This research is to test a new robust foil regenerator to investigate its feasibility for free-piston Stirling engine applications. The robust foil regenerator is configured with parallel flow channels and stiffening ribs. Comparing to complex flows patterns in woven screen and random fiber regenerators, the flows in the parallel plates are simple and the flow losses are considerably small, and the heat transfer coefficient is also relatively small. To test the regenerator experimentally, a regenerator test rig was manufactured. To verify the accuracy of the test rig, a woven screen regenerator was made and tested prior to the robust foil regenerator test. The friction factor correlation of the woven screen regenerator shows a good agreement with the published regenerator correlations. The robust foil regenerator tests were conducted under the oscillating flow condition with different Reynolds numbers. The flow losses were measured with dynamic pressure transducers and the friction coefficient was derived from measured pressure drop. The measured temperature difference over the regenerator was used to calculate the heat transfer coefficient and Nusselt number. The figure of merit with two methods of NPH/NTU ratio and j/f ratio was calculated to evaluate the suitability of the robust foil regenerator. The results show that the robust foil regenerator has considerably higher figure of merit and therefore, better performance than other types of regenerator.

42 ENGINEERING↗
DOE OSTI2021

Adsorption and Surface Diffusion of Metals on α-Al2O3 for Advanced Manufacturing Applications

Abstract The adsorption and diffusion of Mo and Nb adatoms on the $$\alpha$$ α -Al 2 O 3 (0001) surface were explored using density functional theory-based methods. Adsorption energies of Mo and Nb adatoms at minima sites on the surface were determined. Mo and Nb adatoms prefer to adsorb to the same locations on the surface, and larger adsorption energies calculated for Nb compared to Mo indicate that Nb adatom-surface interactions are stronger than Mo. Using minima adsorption sites as initial and final locations for surface diffusion, energy barriers for diffusion were calculated using the nudged elastic band method. Overall, Mo and Nb follow roughly the same diffusion paths. The diffusion pre-factors for Mo and Nb are similar; however, Mo diffusion has a lower energy barrier and thus a larger diffusion coefficient compared with Nb. These results provide insight into the role of surface diffusion of Mo and Nb adatoms during advanced manufacturing processes.

Biaggne, Austin (ORCID:0000000328074960)↗
DOE OSTI2023

Measurement of Beam Energy in the Fermilab's Linac taken at the Transfer Line

Linac is the first machine in the accelerator chain at Fermilab where H$^{-}$ ions are accelerated from 35 keV to 401.5\,MeV and then injected into a synchronton known as Booster where they are stripped of their electrons to become protons. One of the tools used during tuning of the Linac extraction energy is two beam pickups known as Griffin Detectors. Our goal is to control the output energy using machine learning techniques to increase the reliability and quality of the beam delivered from Linac. The first step is to understand the data from the diagnostics to develop reliable and accurate energy measurement, and control methods before implementing machine learning techniques. Two methods of energy measurement were studied, and their results are compared. The first method was the time of flight measurement using Beam Position Monitors that provide beam phase measurement. The second method used the relation between beam transverse positions and dispersion values to calculate momentum variation. The results of these two measurement methods are found to be consistent.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Vacancy-driven variations in the phonon density of states of fast neutron irradiated nuclear graphite

Research examining radiation damage and its effects in graphite began in the 1940’s during the development of moderated nuclear reactors. Interest in this topic is expanding because of emerging applications associated with fullerenes and carbon nanostructures as well as its long-standing use as a fission moderator. In this work, we report the measurements of the full phonon density of states of irradiated nuclear graphite using inelastic neutron scattering experiments at room temperature. The samples were previously exposed to different levels of neutron doses and irradiation temperatures. The phonon density of states of perfect and defected graphite supercells including different configurations of vacancies and interstitials were calculated using first-principles direct method. Despite high neutron irradiation doses, the scattering measurements show that the induced damage is localized, and the layered structure is preserved. In conclusion, a comparison of the measured phonon densities of states of irradiated samples with those calculated of defected supercells indicates that the main changes observed are mainly attributed to formation of vacancies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tree tensor network hierarchical equations of motion based on time-dependent variational principle for efficient open quantum dynamics in structured thermal environments

In this work, we introduce an efficient method, TTN-HEOM, for exactly calculating the open quantum dynamics for driven quantum systems interacting with highly structured bosonic baths by combining the tree tensor network (TTN) decomposition scheme with the bexcitonic generalization of the numerically exact hierarchical equations of motion (HEOM). The method yields a series of quantum master equations for all core tensors in the TTN that efficiently and accurately capture the open quantum dynamics for non-Markovian environments to all orders in the system–bath interaction. These master equations are constructed based on the time-dependent Dirac–Frenkel variational principle, which isolates the optimal dynamics for the core tensors given the TTN ansatz. The dynamics converges to the HEOM when increasing the rank of the core tensors, a limit in which the TTN ansatz becomes exact. We introduce TENSO, tensor equations for non-Markovian structured open systems, as a general-purpose Python code to propagate the TTN-HEOM dynamics. We implement three general propagators for the coupled master equations: two fixed-rank methods that require a constant memory footprint during the dynamics and one adaptive-rank method with a variable memory footprint controlled by the target level of computational error. We exemplify the utility of these methods by simulating a two-level system coupled to a structured bath containing one Drude–Lorentz component and eight Brownian oscillators, which is beyond what can presently be computed using the standard HEOM. Our results show that the TTN-HEOM is capable of simulating both dephasing and relaxation dynamics of driven quantum systems interacting with structured baths, even those of chemical complexity, with an affordable computational cost.

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

Hydrodynamic Coupling to a Homogenized Radiation Transport Method based on Young Measures

Resolving radiation transport fields subject to opacity profiles with strong, oscillatory line structure while potentially falling under intermediate optical depth conditions presents a numerical challenge in radiation transport modeling. The Young measure-based homogenization technique formulated by Haut et al. (2017) was investigated as a candidate method for resolving radiation fields under such conditions more accurately. The method was compared against frequently-utilized mean opacity methods as the Rosseland and Planck formulations. In this work, all methods were tested through radiation slab calculations separately comprised of aluminum, copper, and krypton, each for different thermodynamic conditions. Following these offline radiation slab calculations, demonstrations shifted towards the SCEPTRE radiation transport code and, subsequently, the multiphysics ALEGRA code for approximately-coupled radiation-material simulations. Throughout all the simulations shown in this study, for a fixed computational cost, the homogenized method was observed to be more accurate than any of the solutions determined through traditional mean opacity approaches.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Convergent Protocols for Computing Protein–Ligand Interaction Energies Using Fragment-Based Quantum Chemistry

Fragment-based quantum chemistry methods offer a way to sidestep the steep nonlinear scaling of electronic structure calculations so that large molecular systems can be investigated using high-level methods. Here, we use fragmentation to compute protein–ligand interaction energies in systems with several thousand atoms, using a new software platform for managing fragment-based calculations that implements a screened many-body expansion. Convergence tests using a minimal-basis semiempirical method (HF-3c) indicate that two-body calculations, with single-residue fragments and simple hydrogen caps, are sufficient to reproduce interaction energies obtained using conventional supramolecular electronic structure calculations, to within 1 kcal/mol at about 1% of the computational cost. We also demonstrate that the HF-3c results are illustrative of trends obtained with density functional theory in basis sets up to augmented quadruple-ζ quality. Strategic deployment of fragmentation facilitates the use of converged biomolecular model systems alongside high-quality electronic structure methods and basis sets, bringing ab initio quantum chemistry to systems of hitherto unimaginable size. This will be useful for generation of high-quality training data for machine learning applications.

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

Persistent Protein Motions in a Rugged Energy Landscape Revealed by Normal Mode Ensemble Analysis

Proteins are allosteric machines that couple motions at distinct, often distant, sites to control biological function. Low-frequency structural vibrations are a mechanism of this long-distance connection and are often used computationally to predict correlations, but experimentally identifying the vibrations associated with specific motions has proved challenging. Spectroscopy is an ideal tool to explore these excitations, but measurements have been largely unable to identify important frequency bands. The result is at odds with some previous calculations and raises the question what methods could successfully characterize protein structural vibrations. Here we show the lack of spectral structure arises in part from the variations in protein structure as the protein samples the energy landscape. However, by averaging over the energy landscape as sampled using an aggregate 18.5 μs of all-atom molecular dynamics simulation of hen egg white lysozyme and normal-mode analyses, we find vibrations with large overlap with functional displacements are surprisingly concentrated in narrow frequency bands. These bands are not apparent in either the ensemble averaged vibrational density of states or isotropic absorption. However, in the case of the ensemble averaged anisotropic absorption, there is persistent spectral structure and overlap between this structure and the functional displacement frequency bands. We systematically lay out heuristics for calculating the spectra robustly, including the need for statistical sampling of the protein and inclusion of adequate water in the spectral calculation. The results show the congested spectrum of these complex molecules obscures important frequency bands associated with function and reveal a method to overcome this congestion by combining structurally sensitive spectroscopy with robust normal mode ensemble analysis.

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

Non-Hermitian quantum mechanics approach for extracting and emulating continuum physics based on bound-state-like calculations: Detailed description

Here, this work applies a reduced basis method to study the continuum physics of a finite quantum system—either few or many-body. Specifically, I develop reduced-order models, or emulators, for the underlying inhomogeneous Schrödinger equation and train the emulators against the equation's bound-state-like solutions at complex energies. The emulators rapidly and accurately interpolate and extrapolate the matrix elements of the Hamiltonian resolvent operator (Green's function) across a parameter space that includes both complex energy and other real-valued physical inputs in the Schrödinger equation. The spectra, discretized and compressed as the result of emulation, and the associated resolvent matrix elements (or amplitudes), have the defining characteristics of non-Hermitian quantum mechanics calculations, featuring complex eigenenergies with negative imaginary parts and branch cuts moved below the real axis in the complex energy plane. Therefore, one now has a method that extracts continuum physics from bound-state-like calculations and emulates those extractions in the input parameter space. Building on a prior Letter [Zhang, Phys. Rev. Lett. 135, 242501 (2025)], this article provides the full theoretical details, a comprehensive analysis of the method's performance, and a brief discussion of how it can be coupled with existing continuum approaches to perform emulations in their input parameter spaces.

ab initio calculations↗
DOE OSTI2021

Elaborated thermochemical treatment of HF, CO, N 2 , and H 2 O: Insight into HEAT and its extensions

Here, empirical, highly accurate non-relativistic electronic total atomization energies (eTAEs) are established by combining experimental or computationally converged treatments of the nuclear motion and relativistic contributions with the total atomization energies of HF, CO, N 2 , and H 2 O obtained from the Active Thermochemical Tables. These eTAEs, which have estimated (2σ) uncertainties of less than 10 cm -1 (0.12 kJ mol -1 ), form the basis for an analysis of high-level ab initio quantum chemical calculations that aim at reproducing these eTAEs for the title molecules. The results are then employed to analyze the performance of the high-accuracy extrapolated ab initio thermochemistry, or High-Accuracy Extrapolated Ab Initio Thermochemistry (HEAT), family of theoretical methods. The method known as HEAT-345(Q), in particular, is found to benefit from fortuitous error cancellation between its treatment of the zero-point energy, extrapolation errors in the Hartree-Fock and coupled cluster contributions, neglect of post-(T) core-correlation, and the basis-set error involved in higher-level correlation corrections. In addition to shedding light on a longstanding curiosity of the HEAT protocol—where the cheapest HEAT-345(Q) performs comparably to the theoretically more complete HEAT-456QP procedure—this study lays the foundation for extended HEAT variants that offer substantial improvements in accuracy relative to the established approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗