Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Characterization Study of TestBed Infrastructure Performance in a Distributed Simulation Environment: Baseline Analysis

Characterization of the performance of Air Traffic Management Exploration (ATM-X) TestBed integration environment has been investigated and documented for one system configuration for progressively increasing traffic. Several statistical parameters were used to assess the performance of the TestBed distributed system such as mean, standard deviation, skewness, and kurtosis of latency, and update rate for aircraft state messages that are transmitted through the simulated system under investigation. It is necessary to assess the performance characteristics of distributed systems in terms of the indicated statistical parameters mentioned above. It is critical to verify the system performance with respect to a researcher’s required system performance. Computer host specifications are documented in terms of Central Processing Unit (CPU) clock speed and core count. Transmission Control Protocol/ Internet Protocol (TCP/IP) message protocol was used for data transmission. The system network topology also contributes to the latency and update rate variations from the one imposed by the data source. The motivation for selecting the TestBed infrastructure as the focus of this study can be attributed to the number of services and capabilities it provides that help simplify the process of preparing and conducting a simulation. These capabilities include an easy to use GUI for simulation configuration, access to TestBed library by the end-user of other simulation software components, a modular adapter paradigm that allows simple connectivity of external software to TestBed, connectivity with other simulation laboratories, and a Software Development Kit (SDK) for quicker development. Two types of traffic generators, Air Traffic Generator (ATG) and Multi Aircraft Control System (MACS) were used to generate messages that were injected into the TestBed distributed environment. Eight different air traffic scenarios with progressively increasing loads were generated for each air traffic simulator. The corresponding air traffic loads between the two simulators had an identical number of aircraft per scenario, but different flight plans. It was observed that the performance of MACS degraded for air traffic scenarios containing more than 200 aircraft (37.5 KB/s nominal throughput). However, ATG performed adequately under all tested air traffic loads up to 1200 aircraft (225. KB/s nominal throughput). The tests show that MACS exhibits better latency performance with smaller aircraft loads when compared to ATG. The tests also show that the TestBed infrastructure successfully transmits 1200 aircraft without significant degradation of its performance. From the latency trends for both MACS and ATG, it is clear that as aircraft load increases, the latency in the system increases as well as its standard deviation. Likewise, the trends for the update data rate for both MACS and ATG show that as the aircraft load increases, so does the standard deviation and mean of the update rates which can be attributed to the performance of MACS and ATG applications. The analysis of the results of this study have proven that the overall system performance is dependent on the individual performance of each system component that is connected to TestBed, which subsequently propagates into the system. All TestBed characterization tests were conducted in SimLabs at NASA Ames Research Center in November 2019. This study addresses the need for a baseline TestBed characterization, and the results will serve as a reference for more complex simulation systems.

Air Traffic Management simulations↗
DOE OSTI2022

Diurnal cycle of precipitation over global monsoon systems in CMIP6 simulations

Accurate simulation of the diurnal cycle of precipitation over different climate regimes, including global monsoon systems, is an ongoing challenge for current state-of-the-art climate models. Here, in this study, the diurnal cycles of summertime precipitation from 21 CMIP6 historical simulations are evaluated with the NASA IMERG data via the diurnal harmonic amplitude and phase from the Fourier analysis in seven regional monsoon systems: Northern and Southern Africa monsoons, South Asia monsoon, East Asia monsoon, Australia monsoon, and North and South America monsoons. Particularly, we focus on four subregions including ocean, coastal ocean, land, and coastal land considering the different underlying physical processes. It is found that the CMIP6 historical simulations perform fairly well on reproducing the diurnal cycle of precipitation over ocean and coastal ocean, while the simulated diurnal harmonic phase and amplitude over land and coastal land show larger biases and model spread compared to biases over ocean and coastal ocean. The simulation of precipitation diurnal cycle is substantially improved in the multi-model mean of both CMIP6 historical and CMIP6 AMIP simulations compared with that of CMIP5 historical simulations. But the model bias of precipitation diurnal cycle is in general similar between CMIP6 AMIP and CMIP6 historical simulations, which suggests that the impact of the interactive ocean on the simulated diurnal precipitation averaged over a large domain is generally minor. The influence of model resolution on the diurnal cycle is shown in mixed directions across various subregions, and significant improvement from the low- to high-resolution models is only notable over coastal ocean.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Myna: Connecting powder bed fusion build data to simulation tools for digital twin applications

Additive manufacturing (AM), as a digital process, can generate a detailed digital thread linking a part’s design and manufacturing to its operational performance. As AM systems advance, an increasing amount of process data is stored in manufacturing databases. In principle, this data can be utilized by simulation-based digital twin approaches, such as real-time process control and asynchronous post-processing guidance. However, few tools currently exist for systematically integrating digital thread data with computational tools. Here, in this study, we propose a software package, called Myna, for connecting data from powder bed fusion processes to simulation tools. The utility of such a platform is demonstrated using build data from the Oak Ridge National Laboratory Manufacturing Demonstration Facility “Peregrine v2023-10” public dataset to automatically configure and run 54 semi-analytical 3DThesis melt pool simulations, 78 numerical Additive FOAM melt pool simulations, and 3 ExaCA microstructure simulations. The simulated, spatially registered microstructures are then compared directly with electron backscatter diffraction characterization of the corresponding as-built part locations. The resulting simulated microstructure showed variation as a function of process parameters, particularly stripe width; however, the experimental data had little variation between the microstructure texture and grain size resulting from different processing conditions. Analysis of the discrepancies suggest that it is possible a two-phase ferritic-austenitic solidification model is needed to accurately predict grain size and texture for certain stainless steel 316L feedstock compositions under powder bed fusion conditions, providing direction for future research. As illustrated here, due to the number and complexity of the simulations involved in AM process-structure–property predictions, automated methods to connect process data and simulations will remain necessary tools for testing hypotheses and implementing digital twin applications.

Knapp, Gerald L. [Oak Ridge National Laboratory (O↗
DOE OSTI2017

SQERTSS: Dynamic rank based throttling of transition probabilities in kinetic Monte Carlo simulations

Lattice based Kinetic Monte Carlo (KMC) simulations offer a powerful simulation technique for investigating large reaction networks while retaining spatial configuration information, unlike ordinary differential equations. However, large chemical reaction networks can contain reaction processes with rates spanning multiple orders of magnitude. This can lead to the problem of “KMC stiffness” (similar to stiffness in differential equations), where the computational expense has the potential to be overwhelmed by very short time-steps during KMC simulations, with the simulation spending an inordinate amount of KMC steps / cpu-time simulating fast frivolous processes (FFPs) without progressing the system (reaction network). In order to achieve simulation times that are experimentally relevant or desired for predictions, a dynamic throttling algorithm involving separation of the processes into speed-ranks based on event frequencies has been designed and implemented with the intent of decreasing the probability of FFP events, and increasing the probability of slow process events -- allowing rate limiting events to become more likely to be observed in KMC simulations. This Staggered Quasi-Equilibrium Rank-based Throttling for Steady-state (SQERTSS) algorithm designed for use in achieving and simulating steady-state conditions in KMC simulations. Lastly, as shown in this work, the SQERTSS algorithm also works for transient conditions: the correct configuration space and final state will still be achieved if the required assumptions are not violated, with the caveat that the sizes of the time-steps may be distorted during the transient period.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Active- and transfer-learning applied to microscale-macroscale coupling to simulate viscoelastic flows

Active- and transfer-learning are applied to microscale dynamics of polymer flows for the multiscale discovery of effective constitutive approximations required in viscoelastic flow simulation. The result is macroscopic rheology directly connected to a microstructural model. Micro and macroscale simulations are adaptively coupled by means of Gaussian process regression (GPR) to run the expensive microscale computations only as necessary. This multiscale method is demonstrated with flows of a polymer solution as a model system. At the microscale level dissipative particle dynamics (DPD) is employed to model the fluid as a suspension of bead-spring micro-structures subjected to steady shear flow. The results yield the non-Newtonian viscosity and the first normal stress difference at strain rates as training data used in a GPR model. DPD parameters are calibrated with respect to experimental data for a real polymer solution. Compliance with these data requires adjustment of the DPD model's cutoff radius, which then becomes a function of the second invariant of the strain rate tensor. The FENE-P model is chosen for the macroscale description using the spectral element method (SEM) to simulate channel flow and flow past a circular cylinder. The DPD results at the lowest possible shear strain rate yield an estimate of the zero-shear rate viscosity, which allows the initiation of the macroscale flow by SEM as a Newtonian fluid. The resulting strain-rate field is surveyed to determine additional shear strain rate sampling points for the DPD system. This new information allows an initial fitting of parameters of the constitutive equation followed by new SEM simulations at the macroscale. Additionally, guided by active-learning GPR to select new sampling points, this process continues until convergence is achieved. The effectiveness of this new simulation paradigm for viscoelastic flows is tested with different macroscale operating conditions. The effective closure learned in the channel simulation is then transferred directly to the flow past a circular cylinder at low Reynolds number, where the results show that only two additional DPD simulations are required to achieve a satisfactory constitutive model. With an increase of the Reynolds number, the active-learning scheme automatically detects the inaccuracy of the learned constitutive model, and initiates additional DPD simulations for the extra data needed to once again close the microscale-macroscale coupled system. This new paradigm of active- and transfer-learning for multiscale modeling is readily applicable to other microscale-macroscale coupled simulations of complex fluids and other materials. Furthermore, the coupling between microscale and macroscale solvers can be seamlessly implemented with our open source multiscale universal interface (MUI) library.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Validation of SOLPS-ITER and EDGE2D-EIRENE simulations for H, D, and T JET ITER-like wall low-confinement mode plasmas

Both experiments and simulations with SOLPS-ITER and EDGE2D-EIRENE show that the onset of detachment for the low-field side (LFS) divertor – defined here as the line-averaged upstream density (〈n e 〉 edge ) at which the plasma flux to the LFS target (I LFS-plate ) starts to decrease with increasing 〈n e 〉 edge – is independent of the isotope mass. However, there are three major simulation-experiment discrepancies: (i) the absolute values of I LFS-plate and the electron density (n e ) in the LFS divertor at the onset of detachment are significantly lower in simulations, i.e., approximately a factor of 2 for I LFS-plate and a factor of 3-4 for n e ; (ii) the degree of detachment – defined here as the difference between I LFS-plate at the onset of detachment and at an 〈n e 〉 edge value close to the density limit – is smaller in simulations compared to experiments; and (iii) the experimentally observed larger degree of detachment for D and T plasmas compared to H plasmas cannot be clearly distinguished from the simulation results. There are strong indications that discrepancy (i) is to a large extent caused by neglecting Lyman-opacity effects in our simulations. The simulations predict a similar net volumetric recombination source for all isotopes due to the fact that molecule-activated recombination (MAR) compensates for the reduced electron–ion recombination (EIR) for H, whereas MAR is negligible for D and T. This similar net volumetric recombination source for all isotopes leads to an isotope-independent degree of detachment in simulations. An analysis of the Balmer-α and Balmer-$γ$ emission confirms the underestimate of MAR in simulations (especially for D and T) for the JET metallic wall, which was previously observed for devices with a carbon wall. The underestimate of MAR is an important cause for discrepancy (ii) and the fact that there is a stronger underestimate of MAR for D and T than for H explains discrepancy (iii). Extending the plasma grid to the vessel wall increases I LFS-plate and n e at the onset of detachment by 25%, and the EIR source increases by 80% in detached conditions. Hence, while the extended grid results are closer to the experimental observations, the previously described qualitative discrepancies still persist.

Balmer emission↗
DOE OSTI2018

Adaptive ensemble simulations of biomolecules

Recent advances in both theory and computational power have created opportunities to simulate biomolecular processes more efficiently using adaptive ensemble simulations. Ensemble simulations are now widely used to compute a number of individual simulation trajectories and analyze statistics across them. Adaptive ensemble simulations offer a further level of sophistication and flexibility by enabling highlevel algorithms to control simulations-based on intermediate results. We review some of the adaptive ensemble algorithms and software infrastructure currently in use and outline where the complexities of implementing adaptive simulation have limited algorithmic innovation to date. In conclusion, we describe an adaptive ensemble API to overcome some of these barriers and more flexibly and simply express adaptive simulation algorithms to help realize the power of this type of simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Molecular Simulations of CO 2 and H 2 Solubility, CO 2 Diffusivity, and Solvent Viscosity at 298 K for 27 Commercially Available Physical Solvents

CO 2 and H 2 solubilities, CO 2 /H 2 solubility selectivities, CO 2 diffusivities, and solvent viscosities in 27 commercially available physical solvents at 298 K were calculated from molecular simulations using the CHARMM36 all-atom force field for most solvents, and the simulation results were compared with available experimental data in this report. The van der Waals radius parameters for solvents were slightly tuned to reproduce the experimental solvent density. The simulated CO 2 solubilities are comparable with the experimental data, with an average absolute difference of 28%. For the homologous compounds containing the –(OCH 2 CH 2 )– repeat unit, both simulated and experimental data show that CO 2 solubility decreases when the number of repeat units is increased; CO 2 solubilities in these homologous compounds exhibit almost a perfect positive linear correlation with the solvent free-volume fractions. The simulated H 2 solubilities and CO 2 /H 2 solubility selectivities are also comparable with the experimental data, with differences of 22% and 17%, respectively. The H 2 solubilities in all solvents researched here correlate very well with the solvent free-volume fractions, exhibiting a positive linear correlation coefficient of 0.84. Additionally, simulations show that CO 2 solubility decreases when the temperature is increased. In contrast, H 2 solubility increases at elevated temperature, which is partly due to the increased solvent free-volume fraction at elevated temperature. Lastly, although the viscosity difference tends to be large (30%–246%) between simulation and experiment, both simulated and experimental data exhibit a similar solvent viscosity trend. Furthermore, simulations show that CO 2 diffusivities in solvents are very strongly correlated with the solvent viscosities and the relationship between them is given by D$_{CO_{2}}$ = (2.6 ± 0.3) × 10 –9 /η solvent 0.59 ± 0.03.

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

Simulated Dust Transport in the Convective Boundary Layer

Soil dusts are an important source of aerosol in agricultural regions and can affect the Earth's radiation budget through the modification of cloud properties, and in particular, through their ability to act as ice nucleating particles. In order to impact cloud properties, agricultural soil dusts need to be transported from the point of emission to cloud–relevant altitudes. Vertical transport within the planetary boundary layer is strongly controlled by turbulence and is challenging to represent accurately in regional and global models. Large–eddy simulations (LES) are run at resolutions capable of resolving most of the turbulent energy directly and can thus better simulate vertical transport. In this study, we leverage the LES ARM Symbiotic Simulation and Observation (LASSO) large–eddy simulations to simulate vertical transport of agricultural dust within the turbulent boundary layer using a modified version of the stochastic dispersion model FLEXPART–WRF. We find that the modified model is better capable of simulating particle transport due to turbulence, and that particle size was the greatest factor in determining particle lifetime. Individual meteorology and particle density had intermediate effects upon particle transport, while release height had little effect upon simulation results. Finally, we utilize a quasi–single column model (QSCM) approach to determine how our results compare to a parameterized treatment of turbulence. The QSCM simulations led to greater tracer transport out of the boundary layer, with ramifications for any studies utilizing a Lagrangian stochastic model to understand tracer dispersion. Furthermore, these results highlight the importance of accurately simulating turbulence for understanding particle transport.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Enhancing traffic safety analysis with digital twin technology: integrating vehicle dynamics and environmental factors into microscopic traffic simulation

Traffic safety is a critical concern in transportation engineering and urban planning. Traditional traffic safety analysis requires trained observers to collect data in the field, which is time-consuming, labor-intensive, and sometimes inaccurate. In recent years, microscopic traffic simulation, which simulates individual vehicles’ movements within a transportation network, have been utilized to study traffic safety. However, microscopic traffic simulation only focuses on traffic-related factors, such as traffic volume, traffic signals, and lane configurations, neglecting vehicle dynamics and environment-related factors like weather and lighting conditions, which can significantly impact traffic safety. In light of this, this paper explores the application of digital twin technology in traffic safety analysis, integrating vehicle simulators, which consider vehicle dynamics and environmental factors, and microscopic traffic simulators, which simulate the operations of traffic flow, for enhanced safety evaluations. Various scenarios, including different weather conditions and visibility levels, are simulated using a digital twin of a road segment in Tuscaloosa, Alabama. The simulations employ Surrogate Safety Measures (SSMs) like Time to Collision (TTC) and Deceleration Rate to Avoid a Crash (DRAC) to assess safety under varying conditions. The results demonstrate that traffic digital twin can identify potential safety issues that traditional microscopic simulation cannot, providing insights for improving traffic control strategies and transportation infrastructure to enhance traffic safety.

Civil engineering↗
DOE OSTI2020

Fully compressible simulations of waves and core convection in main-sequence stars

Context: Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and requires rather high artificial viscosity and thermal diffusivity for numerical stability. A direct comparison with observations is thus limited. Aims: In this study, we explore the capabilities of our compressible multidimensional Seven-League Hydro (SLH) code to simulate stellar oscillations. Methods: We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic 3 M ⊙ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results: Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with the observed spectra of massive stars. Conclusion: The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows for the coupling of gravity and pressure waves in the outer convective envelopes of evolved stars to be studied in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Evaluating Connected and Automated Vehicles in Co-Simulation Environment of Traffic Microsimulation and Vehicle Dynamics

Connected and automated vehicles (CAVs) have the potential to improve many aspects of the current transportation systems such as safety, mobility, and energy efficiency. In order to evaluate the benefits and impacts of a CAV, the CAV control algorithm is typically implemented on vehicles simulated in a traffic microsimulation environment. However, traffic microsimulation usually lacks detailed vehicle and powertrain dynamics, making it challenging to fully understand how a CAV control algorithm will perform and respond on an actual vehicle. Whether the same benefits measured in the simulation will also be observed in real-world remains an open question. One potential approach to fill in this gap is to conduct a co-simulation of traffic microsimulation with detailed vehicle and powertrain dynamics models, often developed in MATLAB Simulink. However, current microsimulation tools such as VISSIM and SUMO do not have a ready-to-use interface for co-simulation with vehicle dynamics and Simulink. Also, even if such an interface exists, it will be tool-specific, making it challenging to shift from one tool to another or test CAV controls in different tools. There are needs for tool-agnostic co-simulation as different microsimulation tools have their pros and cons, and researchers often need to use different tools based on the purposes of the simulation, project needs, and applications. In this work, Flexible Interface for X-in-the-loop Simulation (FIXS) is developed that can support the co-simulation of microsimulation, CAV control algorithm, and vehicle dynamics model in Simulink. Enabled by the FIXS, the benefit and performance of a CAV control algorithm can be better understood with the consideration of vehicle responses and dynamics. The connection to VISSIM and SUMO is handled internally by the interface, and users can easily switch tools by changing a configuration file. The co-simulation capability is demonstrated for a VISSIM eco-approach and departure CAV scenario and a SUMO cooperative merging scenario for both a passenger CAV and a class 8 heavy-duty connected and automated trucks.

Shao, Yunli↗
DOE OSTI2018

Equilibration and analysis of first-principles molecular dynamics simulations of water

First-principles molecular dynamics (FPMD) simulations based on density functional theory are becoming increasingly popular for the description of liquids. In view of the high computational cost of these simulations, the choice of an appropriate equilibration protocol is critical. In this study, we assess two methods of estimation of equilibration times using a large dataset of first-principles molecular dynamics simulations of water. The Gelman-Rubin potential scale reduction factor [A. Gelman and D. B. Rubin, Stat. Sci. 7, 457 (1992)] and the marginal standard error rule heuristic proposed by White [Simulation 69, 323 (1997)] are evaluated on a set of 32 independent 64-molecule simulations of 58 ps each, amounting to a combined cumulative time of 1.85 ns. The availability of multiple independent simulations also allows for an estimation of the variance of averaged quantities, both within MD runs and between runs. We analyze atomic trajectories, focusing on correlations of the Kohn-Sham energy, pair correlation functions, number of hydrogen bonds, and diffusion coefficient. The observed variability across samples provides a measure of the uncertainty associated with these quantities, thus facilitating meaningful comparisons of different approximations used in the simulations. We find that the computed diffusion coefficient and average number of hydrogen bonds are affected by a significant uncertainty in spite of the large size of the dataset used. Lastly, a comparison with classical simulations using the TIP4P/2005 model confirms that the variability of the diffusivity is also observed after long equilibration times. Complete atomic trajectories and simulation output files are available online for further analysis.

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

CFD Simulation of Helium Flow Loop Test Section

A helium flow loop is being assembled at Oak Ridge National Laboratory to analyze heat transfer enhancement for systems such as blanket and divertor components. To efficiently identify optimum geometries for heat transfer enhancement in these applications, simulation work is performed to optimize test section designs that are built and tested in the helium flow loop that operates at 4 MPa and a mass flow rate of 100 g/s. Different ribbed geometries that examine rib shape, rib height, rib orientation, rib spacing, and three dimensional orientation are modeled and simulated in STAR-CCM+ to compare their ability to remove heat and mitigate pressure drop. Following the simulations, models are selected and manufactured for the helium flow loop tests. Simulations initially focus on a hydrodynamic study to determine the appropriate mesh and physics models and then add a heat flux to analyze the heat transfer abilities of the models. The simulations are run in steady state and use a Reynolds-averaged Navier-Stokes k-ε turbulence model. The helium is modeled as an ideal gas. The simulation explores models of geometries that enhance the heat transfer and decrease pressure drop with an overall goal of increasing fluid collision with the wall. Enhanced geometries are simulated to select appropriate designs for manufacturing, and preliminary experimental results are used to validate the simulations. Furthermore, the factors that are being analyzed in the comparison between the experimental and the simulated results include matching thermocouple temperatures, pressure drop, roughness, and fluid velocity.

42 ENGINEERING↗
DOE OSTI2026

Numerical-heating effects in atmospheric pressure streamer discharges simulated with a PIC code

Artificial heating in plasma simulations is a well-known phenomenon which occurs when, among other things, the Debye length is poorly resolved by the simulation mesh. Here, in this work, the degree to which numerical-heating occurs during a simulation of a nanosecond atmospheric pressure streamer discharge is examined. The streamer is simulated using a two-dimensional finite-element, particle-in-cell code Empire, which uses direct simulation Monte Carlo for binary particle interactions. Initially, an estimate of the numerical-heating rate applied to Empire is performed using a simple plasma model. Second, a positive atmospheric pressure streamer discharge simulation is performed to study the effects of numerical heating on plasma density, electron temperature, and streamer velocity. The nominal Debye length is approximately 1 μm and the amount of numerical heating introduced in the simulation is varied by using mesh sizes ranging from 2 μm to 20 μm. A measurable numerical heating quantity is proposed that can be used to estimate the appropriate element size and quantify the numerical-heating that can be expected over the simulation time for an atmospheric pressure streamer. In conclusion while Δx/λ D violations can be an issue it is not likely to be an issue with streamer discharges that are temporally short and occur in environments where collision frequencies are high. This result validates the rationale of grid size choices for a large amount of previously published works where Δx/λ D violation was not clearly addressed. Primary finding of this work is that numerical heating is of minor concern for plasma simulations where electron–neutral collisions are numerous such that multiple collisions can occur within a single plasma period.

Nikic, Dejan [University of New Mexico, Albuquerqu↗
DOE OSTI2020

Domain-specific compilers for dynamic simulations of quantum materials on quantum computers

Abstract Simulation of the dynamics of quantum materials is emerging as a promising scientific application for noisy intermediate-scale quantum (NISQ) computers. Due to their high gate-error rates and short decoherence times, however, NISQ computers can only produce high-fidelity results for those quantum circuits smaller than some given circuit size. Dynamic simulations, therefore, pose a challenge as current algorithms produce circuits that grow in size with each subsequent time-step of the simulation. This underscores the crucial role of quantum circuit compilers to produce executable quantum circuits of minimal size, thereby maximizing the range of physical phenomena that can be studied within the NISQ fidelity budget. Here, we present two domain-specific (DS) quantum circuit compilers for the Rigetti and IBM quantum computers, specifically designed to compile circuits simulating dynamics under a special class of time-dependent Hamiltonians. The compilers outperform state-of-the-art general-purpose compilers in terms of circuit size reduction by around 25%–30% as well as wall-clock compilation time by around 40% (dependent on system size and simulation time-step). Drawing on heuristic techniques commonly used in artificial intelligence, both compilers scale well with simulation time-step and system size. Code for both compilers is open-source and packaged into a full-stack quantum simulation software with tutorials included for ease of use for future researchers wishing to perform dynamic simulations of quantum materials on quantum computers. As our DS compilers provide significant improvements in both compilation time and simulation fidelity, they provide a building block for accelerating progress toward physical quantum supremacy.

Physics↗
DOE OSTI2024

Efficient end-to-end simulation of time-dependent coherent X-ray scattering experiments

Physical optics simulations for beamlines and experiments allow users to test experiment feasibility and optimize beamline settings ahead of beam time in order to optimize valuable beam time at synchrotron light sources like NSLS-II. Further, such simulations also help to develop and test experimental data processing methods and software in advance. The Synchrotron Radiation Workshop ( SRW ) software package supports such complex simulations. We demonstrate how recent developments in SRW significantly improve the efficiency of physical optics simulations, such as end-to-end simulations of time-dependent X-ray photon correlation spectroscopy experiments with partially coherent undulator radiation (UR). The molecular dynamics simulation code LAMMPS was chosen to model the sample: a solution of silica nanoparticles in water at room temperature. Real-space distributions of nanoparticles produced by LAMMPS were imported into SRW and used to simulate scattering patterns of partially coherent hard X-ray UR from such a sample at the detector. The partially coherent UR illuminating the sample can be represented by a set of orthogonal coherent modes obtained by simulation of emission and propagation of this radiation through the coherent hard X-ray (CHX) scattering beamline followed by a coherent-mode decomposition. GPU acceleration is added for several key functions of SRW used in propagation from sample to detector, further improving the speed of the calculations. The accuracy of this simulation is benchmarked by comparison with experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2019

DeepDriveMD: Deep-Learning Driven Adaptive Molecular Simulations for Protein Folding

Simulations of biological macromolecules are important in understanding the physical basis of complex processes such as protein folding. However, even with increasing computational capacity and specialized architectures, the ability to simulate protein folding at atomistic scales still remains challenging. This stems from the dual aspects of high dimensionality of protein conformational landscapes, and the inability of atomistic molecular dynamics (MD) simulations to sufficiently sample these landscapes to observe folding events. Machine learning/deep learning (ML/DL) techniques, when combined with atomistic MD simulations offer the opportunity to potentially overcome these limitations by: (1) effectively reducing the dimensionality of MD simulations to automatically build latent representations that correspond to biophysically relevant reaction coordinates (RCs), and (2) driving MD simulations to automatically sample potentially novel conformational states based on these RCs. We examine how coupling DL approaches with MD simulations can lead to effective approaches to fold small proteins on supercomputers. In particular, we study the computational costs and effectiveness of scaling DL-coupled MD workflows implemented using RADICAL-Cybertools in folding two prototypical systems, namely Fs-peptide and the fast-folding variant of the villin head piece protein. We demonstrate that a DL-coupled MD workflow is able to effectively learn latent representations and drive adaptive simulations. Compared to traditional MD-based approaches, our approach achieves an effective performance gain in sampling the folded states by at least 2.3x. Together, our study provides quantitative basis to understand how coupling DL approaches to MD simulations, can lead to effective performance gains and reduced times to solution on supercomputing resources.

Lee, Hyungro↗