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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Spatio-Temporal Surrogates for Interaction of a Jet with High Explosives: Part I - Analysis with a Small Sample Size

Computer simulations, especially of complex phenomena, can be expensive, requiring high-performance computing resources. Often, to understand a phenomenon, multiple simulations are run, each with a different set of simulation input parameters. These data are then used to create an interpolant, or surrogate, relating the simulation outputs to the corresponding inputs. When the inputs and outputs are scalars, a simple machine learning model can suffice. However, when the simulation outputs are vector valued, available at locations in two or three spatial dimensions, often with a temporal component, creating a surrogate is more challenging. In this report, we use a two-dimensional problem of a jet interacting with high explosives to understand how we can build high-quality surrogates. The characteristics of our data set are unique - the vector-valued outputs from each simulation are available at over two million spatial locations; each simulation is run for a relatively small number of time steps; the size of the computational domain varies with each simulation; and resource constraints limit the number of simulations we can run. We show how we analyze these extremely large data-sets, set the parameters for the algorithms used in the analysis, and use simple ways to improve the accuracy of the spatio-temporal surrogates without substantially increasing the number of simulations required.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Prediction challenge: First principles simulation of the ultrafast electron diffraction spectrum of cyclobutanone

Computer simulation has long been an essential partner of ultrafast experiments, allowing the assignment of microscopic mechanistic detail to low-dimensional spectroscopic data. However, the ability of theory to make a priori predictions of ultrafast experimental results is relatively untested. Herein, as a part of a community challenge, we attempt to predict the signal of an upcoming ultrafast photochemical experiment using state-of-the-art theory in the context of preexisting experimental data. Specifically, we employ ab initio Ehrenfest with collapse to a block mixed quantum–classical simulations to describe the real-time evolution of the electrons and nuclei of cyclobutanone following excitation to the 3s Rydberg state. The gas-phase ultrafast electron diffraction (GUED) signal is simulated for direct comparison to an upcoming experiment at the Stanford Linear Accelerator Laboratory. Following initial ring-opening, dissociation via two distinct channels is observed: the C3 dissociation channel, producing cyclopropane and CO, and the C2 channel, producing CH2CO and C2H4. Direct calculations of the GUED signal indicate how the ring-opened intermediate, the C2 products, and the C3 products can be discriminated in the GUED signal. We also report an a priori analysis of anticipated errors in our predictions: without knowledge of the experimental result, which features of the spectrum do we feel confident we have predicted correctly, and which might we have wrong?

Chemistry↗
DOE OSTI2023

Ground Motion Models (GMMs) Improvements Using Earthquake Simulations on High Performance Computers

A computationally efficient simulation platform was developed that can provide representative synthetic ground motions from crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using 3D modeling and high-performance computing. The main objective was to use synthetic ground motion to provide constrains to refinements of exiting ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-5Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. As part of a strategy for selecting a reginal velocity model and validation of developed rupture modeling technique, ground motions from the moment magnitude Mw5.0 November 7, 2016, Cushing Oklahoma, and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes were simulated. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model. Saikia’s model demonstrated better performance in modelling high frequency regional wave propagation for CEUS region. The proposed 3D model includes lateral variations added to the 1D background model using the stochastic scheme of Pitarka and Mellors. Comparisons of the simulations with recordings of both earthquakes demonstrated the reliability of our deterministic simulation approach while emphasizing the importance of including small-scale variability in the regional velocity model needed to reproduce the observed high-frequency wave scattering effects. As part of validation analysis, comparisons with different GMMs for a Mw6.5 earthquake in the CESUS region resulted in a very good match between the simulated and empirical ground motion models. Initial investigations of within-event and between-event ground motion variabilities for Mw6.5 scenario earthquakes on a strike-slip fault, suggest that they are strongly related to spatial slip and slip rate variations, average rupture velocity, rupture area and rupture initiation location. For certain scenarios we found that the ground motion variability observed at near-fault distances (< 5 km) also persists at longer distances. Regardless of the rupture scenario, the simulated ground motion tends to fully saturate at short distances and for all periods. The near-fault saturation has to do with the attenuation of waves propagating along the fault and local rupture radiation pattern that also contribute to stronger ground motion variation at such distances. Analysis of effects of rupture initiation location suggest that the peak ground motion (PGV) and spectral acceleration (SA) can be quite variable due to rupture directivity effects. Such effects are stronger at periods longer than 1s. The effect of the 1D velocity models and surface topography on simulated ground motion were investigated by comparing three component synthetic seismograms computed at selected sites. Effect of surface topography was considered using the ratio between spectral accelerations simulated for two 1D models with flat surface topography and realistic model with surface topography. Overall, the topography slightly amplifies (by ~30%) the ground motion amplitude in the frequency range 1-3Hz. The effect of topography is more visible in the surface and coda waves portion of the seismograms.

58 GEOSCIENCES↗
DOE OSTI2022

A characterization of plasma properties of a heterogeneous magnetized low pressure discharge column

An approach is presented for characterizing heterogeneous magnetized plasma discharge tubes through the scattering of electromagnetic plane waves. Here, we formulate the analytical problem of electromagnetic scattering from a gyrotropic plasma column. The scattering accounts for the heterogeneous composition of the cylindrical discharge plasma and facilitates determining its propensity for gyrotropic scattering, particularly when electron collisional damping may be prevalent. The analytical results are validated using computational simulations. Scattered fields from the magnetized plasma are measured experimentally, and, by comparing the analytical and experimental results, the unknown parameters of the discharge, i.e., characteristic plasma and electron collisional damping frequencies, are determined. The technique is relatively straight-forward to use and removes the need for commercial computational electromagnetic simulations when experimental data on scattering characteristics of such cylindrical discharge plasmas are available.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ten questions concerning agent-based modeling of occupant behavior for energy and environmental performance of buildings

We report the complexity of occupant behavior is one of the major contributors to uncertainty in building performance simulation. Agent-based modeling (ABM), a computational simulation technique, has gained attention in the occupant modeling field due to its capability and flexibility to capture the heterogeneity and dynamics of human behavior and the emergent effects. While multiple efforts in the past decade have demonstrated the usefulness of the ABM approach for simulating occupants and their impacts on building performance, several crucial matters in the ABM research still remain unexplored. This paper presents ten questions that highlight the most important issues regarding ABM research and applications for occupant behavior in the context of building performance simulation. The questions and answers aim to provide insights into current and future ABM research, and more importantly to inspire new significant questions from young researchers in the field. This research is part of the IEA EBC Annex 79 project, occupant-centric building design and operation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Correlated Anion Disorder in Heteroanionic Cubic TiOF 2

Resolving anion configurations in heteroanionic materials is crucial for understanding and controlling their properties. For anion-disordered oxyfluorides, conventional Bragg diffraction cannot fully resolve the anionic structure, necessitating alternative structure determination methods. We have investigated the anionic structure of anion-disordered cubic (ReO 3 -type) TiOF 2 using X-ray pair distribution function (PDF), 19 F MAS NMR analysis, density functional theory (DFT), cluster expansion modeling, and genetic-algorithm structure prediction. Our computational data predict short-range anion ordering in TiOF 2 , characterized by predominant cis-[O 2 F 4 ] titanium coordination, resulting in correlated anion disorder at longer ranges. To validate our predictions, we generated partially disordered supercells using genetic-algorithm structure prediction and computed simulated X-ray PDF data and 19 F MAS NMR spectra, which we compared directly to experimental data. To construct our simulated 19 F NMR spectra, we derived new transformation functions for mapping calculated magnetic shieldings to predicted magnetic chemical shifts in titanium (oxy)fluorides, obtained by fitting DFT-calculated magnetic shieldings to previously published experimental chemical shift data for TiF 4 . We find good agreement between our simulated and experimental data, which supports our computationally predicted structural model and demonstrates the effectiveness of complementary experimental and computational techniques in resolving anionic structure in anion-disordered oxyfluorides. From additional DFT calculations, we predict that increasing anion disorder makes lithium intercalation more favorable by, on average, up to 2 eV, highlighting the significant effect of variations in short-range order on the intercalation properties of anion-disordered materials.

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

Space Charge Simulations of High Intensity Proton Beams in the AGS Booster

Computer simulation studies have been performed to understand the beam behavior and to explore intensity limitations of proton beams in the AGS Booster at higher beam intensities. During the 100 GeV polarized proton operations of RHIC Run 2024, sPHENIX operated in modes with a crossing angle at collisions in order to mitigate beam-beam effects. Three different running modes were employed: (a) sPHENIX operated with a negative (-2 mrad) crossing angle, and STAR operated with 0 mrad. Both experiments were brought into collisions at the start of the store. (b) sPHENIX was brought into collisions with 0 mrad first. Then STAR was brought into collisions after the beam-beam parameter from sPHENIX reduced to below $10 \times 10^{-3}$. (c) sPHENIX operated with a positive (+1.5 mrad) crossing angle, and STAR operated with 0 mrad. Both experiments were brought into collisions at the start of the store. The collisions with a crossing angle of up to $\pm 2$ mrad, as in running modes (a) and (c), lead to large Piwinski angle in the new sPHENIX detector, which reduces luminosity if other parameters are unchanged. There are two ways to compensate the reduction in luminosity: squeeze $\beta^{*}$ if there is sufficient dynamic aperture, or increase the injected beam intensity. The first part of polarized proton operation during RHIC Run 2024 was dedicated to increasing the intensity. Different configurations were tested with crossing angle and lattice adjustments on RHIC. At the same time, new injector configurations were developed and tested in an effort to push for both higher intensity and better quality of the beam injected into RHIC. When the beam intensity is increased, space charge is a concern particularly in the lower energy stages of acceleration, such as during the injection and the early part of the Booster cycle, which could become a dominant effect in limiting the intensity of the beam that can be delivered to RHIC.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Hands-on with OWL: the Oak–Ridge Wang–Landau Monte Carlo software suite

The Oak–Ridge Wang–Landau (OWL) package is an open-source scientific software specialized for large-scale, Monte Carlo simulations for the study of materials properties at finite temperature. In this paper, we discuss the main features and capabilities of OWL, followed by detailed descriptions of building and running the code. The readers will be guided through the usage and functionality of the code with a few hands-on examples. This paper is based on a tutorial on OWL given at the 32nd Center for Simulational Physics Workshop on Recent Developments in Computer Simulation Studies in Condensed Matter Physics.

Li, Ying Wai↗
DOE OSTI2020

On the Stark Effect of the O I 777-nm Triplet in Plasma and Laser Fields

The O I 777-nm triplet transition is often used for plasma density diagnostics. It is also employed in nonlinear optics setups for producing quasi-comb structures when pumped by a near-resonant laser field. Here, we apply computer simulations to situations of the radiating atom subjected to the plasma microfields, laser fields, and both perturbations together. Our results, in particular, resolve a controversy related to the spectral line anomalously broadened in some laser-produced plasmas. The importance of using time-dependent density matrix is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Simulating the Impact of Dynamic Rerouting on Metropolitan-scale Traffic Systems

The rapid introduction of mobile navigation aides that use real-time road network information to suggest alternate routes to drivers is making it more difficult for researchers and government transportation agencies to understand and predict the dynamics of congested transportation systems. Computer simulation is a key capability for these organizations to analyze hypothetical scenarios; however, the complexity of transportation systems makes it challenging for them to simulate very large geographical regions, such as multi-city metropolitan areas. In this article, we describe enhancements to the Mobiliti parallel traffic simulator to model dynamic rerouting behavior with the addition of vehicle controller actors and vehicle-to-controller reroute requests. The simulator is designed to support distributed-memory parallel execution using discrete event simulation and be scalable on high-performance computing platforms. We demonstrate the potential of the simulator by analyzing the impact of varying the population penetration rate of dynamic rerouting on the San Francisco Bay Area road network. Using high-performance parallel computing, we can simulate a day in the San Francisco Bay Area with 19 million vehicle trips with 50 percent dynamic rerouting penetration over a road network with 0.5 million nodes and 1 million links in less than three minutes. We present a sensitivity study on the dynamic rerouting parameters, discuss the simulator’s parallel scalability, and analyze system-level impacts of changing the dynamic rerouting penetration. Furthermore, we examine the varying effects on different functional classes and geographical regions and present a validation of the simulation results compared to real-world data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The microscopic mechanism of bulk melting of ice

We study the initial stages of homogeneous melting of a hexagonal ice crystal at coexistence and at moderate superheating. Our trajectory-based computer simulation approach provides a comprehensive picture of the events that lead to melting, from the initial accumulation of 5+7 defects, via the formation of L-D and interstitial-vacancy pairs, to the formation of a liquid nucleus. Of the different types of defects that we observe to be involved in melting, a particular kind of 5+7 type defect (type 5) plays a prominent role as it often forms prior to the formation of the initial liquid nucleus and close to the site where the nucleus forms. Hence, like other solids, ice homogeneously melts via the prior accumulation of defects.

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

Crowding-Induced Spatial Organization of Gene Expression in Cell-Sized Vesicles

Cell-free protein synthesis is an important tool for studying gene expression and harnessing it for applications. In cells, gene expression is regulated in part by the spatial organization of transcription and translation. Unfortunately, current cell-free approaches are unable to control the organization of molecular components needed for gene expression, which limits the ability to probe and utilize its effects. In this work, we show, using complementary computational and experimental approaches, that macromolecular crowding can be used to control the spatial organization and translational efficiency of gene expression in cell-sized vesicles. Computer simulations and imaging experiments reveal that, as crowding is increased, DNA plasmids become localized at the inner surface of vesicles. Ribosomes, in contrast, remain uniformly distributed, demonstrating that crowding can be used to differentially organize components of gene expression. We further carried out cell-free protein synthesis reactions in cell-sized vesicles and quantified mRNA and protein abundance. At sufficiently high levels of crowding, we observed localization of mRNA near vesicle surfaces, a decrease in translational efficiency and protein abundance, and anomalous scaling of protein abundance as a function of vesicle size. These results are consistent with high levels of crowding causing altered spatial organization and slower diffusion. Our work demonstrates a straightforward way to control the organization of gene expression in cell-sized vesicles and provides insight into the spatial regulation of gene expression in cells.

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

On the characterization of biases arising from methods and approximations used for sensitivity analyses

The preliminary assessment of the performance and safety of advanced reactors, as well as the identification of research and development needs, rely on computer simulations due to the lack of operational experience. Those simulations should be based on verified and validated computational tools, calculation schemes, and nuclear data libraries, and should be accompanied of a proper characterization of the involved uncertainties, providing confidence in the computational-based predictions. Thus, sensitivity and uncertainty studies, together with integral experiments, play an essential role in that process. Sensitivities are then a critical element and guidelines about how to produce accurate enough sensitivities depending on the subsequent analyses to be performed are of interest. The present study addresses the sensitivities biases arising from the use of different nuclear data libraries, computational methods and the assumption of modelling simplifications. Not only sensitivities for multiplication factor but also for safety-relevant reactivity responses have been analyzed. In the study, propagated uncertainties are also included for the considered parameters to illustrate the impact of the sensitivity profiles obtained for each case on the final uncertainty. Then, relevant recommendations are given for sensitivity analyses of highly complex systems. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

An analysis of Potential Benefits of Adding Diurnal Thermal Storage in Geothermal Heat Pump Systems

Geothermal heat pumps (GHP) are highly efficient for space heating and cooling, utilizing the subsurface of the ground as seasonal energy storage through borehole heat exchangers (i.e., heat stored in summer is used to provide heating in winter, and cold stored in winter is used to provide cooling in summer). On the other hand, diurnal thermal energy storage can actively be charged and discharged for load shifting. Conventional thermal energy storage uses tanks filled with thermal storage media (such as phase change materials) and heat exchangers. It occupies building floor space or land areas outside the building, which may not be available especially in residential buildings. This paper will introduce a novel configuration that can utilize existing GHP system for both diurnal and seasonal thermal storage. This system can actively charge the ground on daily basis by injecting heat or cold in borehole heat exchangers. A simulation model was created for the novel configuration. Computer simulations were conducted using the model for the new GHP system serving a typical single-family house in three U.S. cities representing hot, warm, and cold climates. Based on simulation results, the potential benefits of load shifting and downsize of the expensive ground heat exchangers used by the GHP system were assessed.

Anees, Fady [ORNL]↗
DOE OSTI2024

R-Adaptivity to Enable Compression of Elementary Computations in Extreme-Scale Finite Element Simulators

Modern computing systems are capable of exascale calculations, which are revolutionizing the development and application of high-fidelity numerical models in computational science and engineering. While these systems continue to grow in processing power, the available system memory has not increased commensurately, and electrical power consumption continues to grow. A predominant approach to limit the memory usage in large-scale applications is to exploit the abundant processing power and continually recompute many low-level simulation quantities, rather than storing them. However, this approach can adversely impact the throughput of the simulation and diminish the benefits of modern computing architectures. We present three novel contributions to reduce the memory burden while maintaining, and sometimes improving, performance in simulations based on finite element discretizations. The first contribution develops dictionary-based data compression schemes that detect and exploit the structure of the discretization, due to redundancies across the finite element mesh. While these schemes are shown to reduce memory requirements by more than 99% on meshes with large numbers of identical mesh cells, there are applications where this structure does not exist. The second contribution leverages a recently developed augmented Lagrangian optimization algorithm to enable r-adaptivity for meshes with the goal of enhancing the redundancies in the mesh. The third contribution extends these methods to patch-based linear solvers and preconditioners by compressing local matrices. Numerical results demonstrate the effectiveness of the proposed methods to detect, enhance and exploit mesh structure on a suite of examples inspired by large-scale applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Constant-depth circuits for dynamic simulations of materials on quantum computers

Abstract Dynamic simulation of materials is a promising application for near-term quantum computers. Current algorithms for Hamiltonian simulation, however, produce circuits that grow in depth with increasing simulation time, limiting feasible simulations to short-time dynamics. Here, we present a method for generating circuits that are constant in depth with increasing simulation time for a specific subset of one-dimensional (1D) materials Hamiltonians, thereby enabling simulations out to arbitrarily long times. Furthermore, by removing the effective limit on the number of feasibly simulatable time-steps, the constant-depth circuits enable Trotter error to be made negligibly small by allowing simulations to be broken into arbitrarily many time-steps. For an N -spin system, the constant-depth circuit contains only $\mathcal {O}(N^{2})$ O ( N 2 ) CNOT gates. Such compact circuits enable us to successfully execute long-time dynamic simulation of ubiquitous models, such as the transverse field Ising and XY models, on current quantum hardware for systems of up to 5 qubits without the need for complex error mitigation techniques. Aside from enabling long-time dynamic simulations with minimal Trotter error for a specific subset of 1D Hamiltonians, our constant-depth circuits can advance materials simulations on quantum computers more broadly in a number of indirect ways.

Materials simulation↗
DOE OSTI2026

Interpolation of computed gamma-ray detector response functions

Gamma-ray spectra measured by traditional detectors contain features that result from a combination of the effects of detector materials/geometry, the incident gamma-ray energy, and the angle of entry. The features, such as the full-energy photopeak, Compton continuum, annihilation peak, and escape peaks, are governed by simple relationships depending on incident energy and have been known for a long time. Monte Carlo computer simulations of gamma rays interacting with a detector will show these features, and with a resolution function applied, the results should look similar to real measurements. The traditional approach to creating a detector response function requires many separate simulations of monoenergetic gamma rays striking the detector. This paper presents a new approach to developing computed detector response functions. The new approach involves a much smaller number of monoenergetic gamma-ray simulations and uses interpolation to quickly generate the responses of gamma rays that were not simulated. During the interpolation process, the underlying physics equations are used to accurately compute the response of a given energy gamma ray from the small set of simulations. Such work enables accelerated generation of synthetic radiation detector data.

Detector response↗
DOE OSTI2022

Dataset of Simulated Intracardiac Transmembrane Voltage Recordings and ECG Signals

The dataset consists of pairs of computationally simulated intracardiac transmembrane voltage recordings and ECG signals. In total, 16140 organ-level simulations were conducted to create this dataset, using a range of cardiac geometries and physiological parameters. Simulations were performed at LLNL's Lassen supercomputer, concurrently utilizing 4 GPUs and 40 CPU cores. Each simulation produced pairs of 500ms-by-10 ECG signals and 500ms-by-75 transmembrane voltage signals. For convenience, the signals are concatenated and saved as matrices. Each of these matrices is then stored as a numpy array. See the documentation for further details.

Landajuela, Mikel↗