Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer simulation”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

First Principles Nonadiabatic Excited-State Molecular Dynamics in NWChem

Computational simulation of non-adiabatic molecular dynamics is an indispensable tool for understanding complex photoinduced processes such as internal conversion, energy transfer, charge separation, and spatial localization of excitons, to name a few. Here, we report an implementation of the fewest-switches surface hopping algorithm in the NWChem computational chemistry program. Here the surface hopping method is combined with linear-response time-dependent density functional theory calculations of adiabatic excited state potential energy surfaces. To treat quantum transitions between arbitrary electronic Born{Oppenheimer states, we have implemented both numerical and analytical differentiation schemes for derivative non-adiabatic couplings. A numerical approach for the time-derivative non-adiabatic couplings together with an analytical method for calculating non-adiabatic coupling vectors is an efficient combination for surface hopping approaches. Additionally, electronic decoherence schemes and a state reassigned unavoided crossings algorithm are also implemented to improve the accuracy of the simulated dynamics and to handle trivial unavoided crossings. We apply our code to study the ultrafast decay of photoexcited benzene, including a detailed analysis of the potential energy surface, population decay time scales, and vibrational coordinates coupled to the excitation dynamics. The development presented in this work is a baseline for future implementations of more sophisticated frameworks for simulating of non-adiabatic molecular dynamics in NWChem.

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

A scalable matrix-free spectral element approach for unsteady PDE constrained optimization using PETSc/TAO

In this work, we provide a new approach for the efficient matrix-free application of the transpose of the Jacobian for the spectral element method for the adjoint-based solution of partial differential equation (PDE) constrained optimization. This results in optimizations of nonlinear PDEs using explicit integrators where the integration of the adjoint problem is not more expensive than the forward simulation. Solving PDE constrained optimization problems entails combining expertise from multiple areas, including simulation, computation of derivatives, and optimization. The Portable, Extensible Toolkit for Scientific computation (PETSc) together with its companion package, the Toolkit for Advanced Optimization (TAO), is an integrated numerical software library that contains an algorithmic/software stack for solving linear systems, nonlinear systems, ordinary differential equations, differential algebraic equations, and large-scale optimization problems and, as such, is an ideal tool for performing PDE-constrained optimization. This paper describes an efficient approach in which the software stack provided by PETSc/TAO can be used for large-scale nonlinear time-dependent problems. Time integration can involve a range of high-order methods, both implicit and explicit. The PDE-constrained optimization algorithm used is gradient-based and seamlessly integrated with the simulation of the physical problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Cement sensors with acoustic bandgaps using carbon nanotubes

Cement is widely used in wellbores to stabilize the steel casing used in wellbore operations for oil and gas production, enhanced geothermal systems and carbon sequestration, and to limit fluid movement between sub-surface strata. Flaws such as microcracks in wellbore cement can lead to leakage along the wellbore compromising wellbore integrity. There is an increasing need for methods to monitor cement crack propagation in wellbore environments. In this study, we develop and report the first cementitious sensors capable of exhibiting high frequency acoustic bandgaps (ABGs) using carbon nanotubes (CNTs). Computational simulations of a sensor unit cell are used to design cement-multi walled carbon nanotubes (MWCNTs) sensors that show a wide bandgap. When the cement-MWCNTs sensors is embedded in cement specimens, bandgaps were measured experimentally under 300 kHz and under 600 kHz, consistent with the computationally predicted bandgaps in the range of 290–360 kHz, 410–460 kHz and 515–585 kHz. These bandgap features were absent in homogeneous cement specimens. X-ray tomographic reconstructions showed microscopic debonding at cement-MWCNTs sensor interface. Frequency response analysis of a three-dimensional computational model indicated a shift of frequency of minimum transmission due to the interface debonding, but no perturbation of bandgap response was observed. Here, the cement-MWCNTs sensors developed in this study show the potential of a packed CNT inclusion material in cementitious matrix to create ABGs in a cement matrix.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dispersion-enhanced sequential batch sampling for adaptive contour estimation

In computer simulation and optimal design, sequential batch sampling offers an appealing way to iteratively stipulate optimal sampling points based upon existing selections and efficiently construct surrogate modeling. Nonetheless, the issue of near duplicates poses tremendous quandary for sequential learning. It refers to the situation that selected critical points cluster together in each sampling batch, which are individually but not collectively informative towards the optimal design. Near duplicates severely diminish the computational efficiency as they barely contribute extra information towards update of the surrogate. To address this issue, we impose a dispersion criterion on concurrent selection of sampling points, which essentially forces a sparse distribution of critical points in each batch, and demonstrate the effectiveness of this approach in adaptive contour estimation. Specifically, we adopt Gaussian process surrogate to emulate the simulator, acquire variance reduction of the critical region from new sampling points as a dispersion criterion, and combine it with the modified expected improvement (EI) function for critical batch selection. The critical region here is the proximity of the contour of interest. This proposed approach is vindicated in numerical examples of a two-dimensional four-branch function, a four-dimensional function with a disjoint contour of interest and a time-delay dynamic system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Uncertainty Quantification in Atomistic Modeling of Metals and Its Effect on Mesoscale and Continuum Modeling: A Review

The design of next-generation alloys through the integrated computational materials engineering (ICME) approach relies on multiscale computer simulations to provide thermodynamic properties when experiments are difficult to conduct. Atomistic methods such as density functional theory (DFT) and molecular dynamics (MD) have been successful in predicting properties of never before studied compounds or phases. However, uncertainty quantification (UQ) of DFT and MD results is rarely reported due to computational and UQ methodology challenges. Over the past decade, studies that mitigate this gap have emerged. These advances are reviewed in the context of thermodynamic modeling and information exchange with mesoscale methods such as the phase-field method (PFM) and calculation of phase diagrams (CALPHAD). The importance of UQ is illustrated using properties of metals, with aluminum as an example, and highlighting deterministic, frequentist, and Bayesian methodologies. Finally, challenges facing routine uncertainty quantification and an outlook on addressing them are also presented.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Comparing Emulation Methods for Computer Models with High Dimensional Output

This Master’s Capstone project will present a comparison of statistical models for computer simulation studies. The four models included in this comparison study were chosen for both their proven and diverse methodologies. We will present a case study on hurricane flood data in the Delaware bay which highlights the strengths and weaknesses of each model when applied to a very large spatial field. As computers have gotten faster, we have become interested in modeling increasingly large spatial fields in both size and resolution. Statistical algorithms that are able to efficiently handle these fields have never been more important. We therefore find this comparison to be extremely topical.

42 ENGINEERING↗
DOE OSTI2022

Monte Carlo simulations demonstrating physics of equivalency of gamma, electron-beam, and X-ray for radiation sterilization

We report the sterilization of medical devices using the gamma rays from the decay of cobalt-60 has accumulated decades of experience of the performance of the materials and devices that are irradiated. The use of radiation using electron beams and x-rays has much less experience and this leads to questions of equivalency between these three technologies. Computer simulations were conducted to model the relevant physical processes of the interactions of each of the three forms of radiation in order to compare the spectra of electron energies at energies below 500 keV. It is predominantly the electrons below this threshold that produce the sterilization dose. No difference in energy spectra was seen between the three types of initial radiation. It is concluded that there is no energy dependent difference between gamma, e-beam, and x-ray for radiation sterilization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Ablation loading of barium ions into a surface-electrode trap

Trapped-ion quantum information processing may benefit from qubits encoded in isotopes that are practically available in only small quantities, e.g., due to low natural abundance or radioactivity. Laser ablation provides a method of controllably liberating neutral atoms or ions from low-volume targets, but energetic ablation products can be difficult to confine in the small ion-electrode distance, micron-scale microfabricated traps amenable to high-speed, high-fidelity manipulation of ion arrays. Here, we investigate ablation-based ion loading into surface-electrode traps of different sizes to test a model describing ion loading probability as a function of effective trap volume and other trap parameters. We characterize loading of ablated barium from a metallic source in two cryogenic surface-electrode traps with 730 and 50 μm ion-electrode distances. Furthermore, our loading rate agrees with a predictive analytical model, providing insight for the confinement of limited-quantity species of interest for quantum computing, simulation, and sensing.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Physics-Based Optical Neuromorphic Classification

Typical approaches to classify scenes from light convert the light field to electrons to perform the computation in the digital electronic domain. This conversion and downstream computational analysis require significant power and time. Diffractive neural networks have recently emerged as unique systems to classify optical fields at lower energy and high speeds. Previous work has shown that a single layer of diffractive metamaterial can achieve high performance on classification tasks. In analogy with electronic neural networks, it is anticipated that multilayer diffractive systems would provide better performance, but the fundamental reasons for the potential improvement have not been established. In this work, we present extensive computational simulations of two - layer diffractive neural networks and show that they can achieve high performance with fewer diffractive features than single layer systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Vibrational Signatures of Electronic Properties in Renewable-Energy Catalysis

The objective of this research program is to discover and develop new approaches for ab initio computational simulations of molecular motion. Specifically, this program aims to decipher the connection between the underlying electronic structure and the accordant molecular vibrations in reactive ions and radicals for renewable-energy purposes. Recent experimental progress in ion sources and optical spectroscopies has unearthed considerable new inner-sphere detail for these complexes, but the connection between these spectral signatures and mechanistic information often remains elusive, to the continued frustration of experimentalists. For this purpose, new anharmonic vibrational frequency methods, along with a publicly deployed software package, will be developed. Working closely with committed experimental collaborators, this conceptual and computational framework will be used to explain the results of new spectroscopy experiments, focusing specifically on the inner-shell mechanisms of renewable-energy catalysis. The oxidation half of catalytic water-splitting chemistry will be a central focus, along with fundamental studies of the manner in which strong ions and radicals activate solvent as a chemical species. The resulting products of the research program will include openly available software and algorithms for the ab initio simulation of challenging vibrational spectra, as well as critical mechanistic insight into energy-focused catalytic processes that are opaque to other existing analytical techniques.

42 ENGINEERING↗
DOE OSTI2020

Solvent–Solvent Correlations across Graphene: The Effect of Image Charges

Wetting experiments show pure graphene to be weakly hydrophilic, but its contact angle (CA) also reflects the character of the supporting material. Measurements and Molecular Dynamics simulations on suspended and supported graphene often reveal a CA reduction due to the presence of the supporting substrate. A similar reduction is consistently observed when graphene is wetted from both sides. The effect has been attributed to transparency to molecular interactions across the graphene sheet, however, the possibility of substrate-induced graphene polarization has also been considered. Computer simulations of CA on graphene have so far been determined by ignoring the material’s conducting properties. We improve the graphene model by incorporating its conductivity according to the Constant Applied Potential Molecular Dynamics. Using this method, we compare the wettabilities of suspended graphene and graphene supported by water by measuring the CA of cylindrical water drops on the sheets. The inclusion of graphene conductivity and concomitant polarization effects lead to a lower CA on suspended graphene but the CA reduction is significantly bigger when the sheets are also wetted from the opposite side. The stronger adhesion is accompanied by a profound change in the correlations among water molecules across the sheet. While partial charges on water molecules interacting across an insulator sheet attract charges of the opposite sign, apparent attraction among like charges is manifested across the conducting graphene. The change is associated with graphene polarization, as the image charges inside the conductor attract equally signed partial charges of water molecules on both sides of the sheet. Additionally, by using a non-polar liquid (diiodomethane), we affirm a detectable wetting translucency when liquid-liquid forces are dominated by dispersive interactions. Furthermore, our findings are important for predictive modeling toward a variety of applications including sensors, fuel cell membranes, water filtration, and graphene-based electrode materials in high-performance supercapacitors.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Investigations of water/oxide interfaces by molecular dynamics simulations

Water/oxide interfaces are ubiquitous on earth and show significant influence on many chemical processes. For example, understanding water and solute adsorption as well as catalytic water splitting can help build better fuel cells and solar cells to overcome our looming energy crisis; the interaction between biomolecules and water/oxide interfaces is one hypothesis to explain the origin of life. However, knowledge in this area is still limited due to the difficulty of studying water/solid interfaces. As a result, research using increasingly sophisticated experimental techniques and computational simulations has been carried out in recent years. Although it is difficult for experimental techniques to provide detailed microscopic structural information, molecular dynamics (MD) simulations have satisfactory performance. In this report, we discuss classical and ab initio MD simulations of water/oxide interfaces. Generally, we are interested in the following questions: How do solid surfaces perturb interfacial water structure? How do interfacial water molecules and adsorbed solutes affect solid surfaces and how do interfacial environments affect solvent and solute behavior? Finally, we discuss progress in the application of neural network potential based MD simulations, which offer a promising future because this approach has already enabled ab initio level accuracy for very large systems and long trajectories.

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

Beam Breakup Simulations for the Scorpius Flash-Radiography Accelerator

Beam breakup (BBU) in linear induction accelerators (LIAs) used for flash radiography is problematic because the high-frequency beam motion can blur the source spot, thereby degrading resolution. Amplification of BBU depends directly on details of cell design and is suppressed in operation by the solenoidal magnetic field focusing on the electron beam. Therefore, much effort has gone into design of the cells and magnetic focusing for the new Scorpius radiography LIA. Here, in this article, we use computer simulations to demonstrate that BBU in Scorpius should be no more than in present radiography LIAs at the Los Alamos National Laboratory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Special Analyses for the Hanford Integrated Disposal Facility Performance Assessment - 20102

In 2014, the Department of Energy (DOE) Office of River Protection and its contractors began to develop a performance assessment for the near-surface disposal of low-level and mixed low-level waste at the Hanford Site's Integrated Disposal Facility (IDF). The IDF is a doubly-lined landfill that was constructed between 2004 and 2006 to be the disposal facility for the vitrified low-activity waste that will be produced at the Waste Treatment and Immobilization Plant (WTP). IDF is also expected to receive solid secondary waste produced at the WTP and other solid wastes from site activities. The IDF has been in a preoperational state awaiting authorization from DOE and a RCRA permit modification from the State of Washington Department of Ecology to receive waste. Both the Disposal Authorization Statement and permit modification require a performance assessment demonstrating that the system of engineered and natural features will limit releases of radionuclides and hazardous chemicals from the IDF and be protective of human health and the environment. The simulated duration is 10,000 years. Based on the analyses presented in the 2017 Integrated Disposal Facility Performance Assessment, DOE issued a conditional Operating Disposal Authorization Statement for the IDF in June 2018. The long-term performance of the IDF to be protective of human health and the environment was evaluated under the requirements of DOE Order 435.1, Radioactive Waste Management. Computer simulations were performed to evaluate whether or not the IDF would comply with DOE requirements. In the time that has passed since the performance assessment was approved by DOE, new information has been discovered that had not been considered in the performance assessment. Since this new information has not been evaluated, the potential impact of the changes have not been taken into consideration in DoE's disposal authorization. DOE and its contractors follow a change control process to screen and, when necessary, evaluate new information that could potentially impact the conclusions of the completed performance assessment. This paper will describe the change control process and provide two examples of evaluations performed following the change control process. The first example evaluates a new waste form for liquid secondary waste that was not evaluated in the performance assessment. In the performance assessment, liquid secondary waste was assumed to be solidified with grout. A new recommendation to dispose of the liquid secondary waste after drying it to a powder was evaluated. The second example evaluates inventory implications from changes to the flow sheet that estimates the feed composition to the low-activity waste vitrification facility. The changes result in higher strontium concentrations in the vitrified waste stream. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Fluid modeling of low-temperature plasmas

Fluid models are essential for understanding and predicting low-temperature plasma (LTP) behavior in various scientific and industrial settings. This paper provides an introductory tutorial on fluid modeling of LTPs, covering model formulation, implementation, and computational simulations. The tutorial focuses on five main components of the formulation of LTP fluid models: fluid flow, energy, chemistry, electromagnetism, and material properties, as well as in essential aspects of model implementations, including multiscale phenomena, multiphysics coupling, and numerical convergence. Designed for students and early-career researchers, this work offers a practical foundation for developing and using fluid models, from in-house computational codes to commercial software, bridging fundamental theory with real-world applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Electron–Ion Temperature Ratio in Astrophysical Shocks

Collisionless shock waves in supernova remnants and the solar wind heat electrons less effectively than they heat ions, as is predicted by kinetic simulations. However, the values of T e /T p inferred from the Hα profiles of supernova remnant shocks behave differently as a function of Mach number or Alfvén Mach number than what is measured in the solar wind or predicted by simulations. Here we determine T e /T p for supernova remnant shocks using Hα profiles, shock speeds from proper motions, and electron temperatures from X-ray spectra. We also improve the estimates of sound speed and Alfvén speed used to determine Mach numbers. We find that the Hα determinations are robust and that the discrepancies among supernova remnant shocks, solar wind shocks, and computer-simulated shocks remain. We discuss some possible contributing factors, including shock precursors, turbulence, and varying preshock conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Seamless Wireless Communication Platform for Internet of Things Applications

The rapid growth of the Internet of Things (IoT) devices resulted in the proliferation of wireless technologies to cater to their increasing data rate requirements and support multiple applications. However, such ever-increasing wireless technologies present numerous challenges such as incompatible wireless standards, increased energy consumption, and insecure communication. The traditional gateways proposed in the literature suffers from limitations such as computational complexity, resource requirements, increased cost, and device size. We vision an era of seamless wireless communication to alleviate the aforementioned challenges in IoT applications. through three inter-dependent functionalities namely detection and identification of wireless technologies, energy-efficient transmit power control, and secure end-to-end communication. To prove the concept, a new gateway is proposed to achieve these three functionalities with only physical layer measurements so that the different communication protocols in the higher layers can be avoided. Novel schemes are conceptualized for resource-limited seamless IoT applications. Moreover, the conceptual seamless IoT platform is validated through software-based computer simulation and software-defined radio-based testbed implementation. Finally, the preliminary analysis demonstrates that the proposed platform has great potential in advancing seamless IoT applications.

97 MATHEMATICS AND COMPUTING↗