Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SIMULATOR”

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 559 records · Page 31

Simulating Pyrocumulonimbus Clouds Using a Multiscale Wildfire Simulation Framework

Pyrocumulonimbus (pyroCb) clouds, driven by extreme fires under favorable meteorological conditions, can inject smoke into the stratosphere at magnitudes comparable to those of moderate volcanic eruptions, potentially altering the global radiative balance and atmospheric composition. However, simulating pyroCb is particularly challenging in Earth system models. Using the Energy Exascale Earth System Model (E3SM), we developed a novel global multiscale framework to model pyroCb events in California, which includes a high‐resolution fire radiative power time series, a one‐dimensional plume‐rise parameterization, a fire‐induced vertical water vapor transport scheme, and a surface wildfire sensible heat flux representation. Our simulation successfully reproduces many pyroCb features, including cloud height, spatiotemporal evolution, and convective intensity in comparison with satellite and ground‐based observations. Sensitivity experiments show that realistic pyroCb simulation depends on vertical water vapor transport. These advances provide a basis for future exploration of pyroCb impacts at regional and global scales within climate models.

E3SM↗

A virtual squeeze on chemistry: Accelerated reaction simulations [Nanoreactor Simulations]

We describe how molecular simulations have the potential to give valuable insights into experimental results, but can be limited by the time- and length-scales they can simulate. Now, reactive chemistry can be driven through a novel simulation approach, which could have ramifications for many research areas, including astrobiology and the origins of life.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The HalfDome multi-survey cosmological simulations: N-body simulations

Upcoming cosmological surveys have the potential to reach groundbreaking discoveries on multiple fronts, including the neutrino mass, dark energy, and inflation. Most of the key science goals require the joint analysis of datasets from multiple surveys to break parameter degeneracies and calibrate systematics. To realize such analyses, a large set of mock simulations that realistically model correlated observables is required. In this paper we present the N-body component of the HalfDome cosmological simulations, designed for the joint analysis of Stage-IV cosmological surveys, such as Rubin LSST, Euclid, SPHEREx, Roman, DESI, PFS, Simons Observatory, CMB-S4, and LiteBIRD. Our 300TB initial data release includes full-sky lightcones and halo catalogs between z = 0–4 for 11 fixed cosmology realizations, as well as an additional run with local primordial non-Gaussianity (f NL = 20). The simulations evolve 61443 particles in a 3.75 h -1 Gpc box, reaching a minimum halo mass of ∼6 × 1012 h -1 M ⊙ and maximum scale of k ∼ h Mpc-1. Our data is publicly available: instructions to access the data and plans for future data releases can be found at https://halfdomesims.github.io.

Bayer, Adrian E↗

Simulation of argon response and light detection in the DarkSide-50 dual phase TPC

A Geant4-based Monte Carlo package named G4DS has been developed to simulate the response of DarkSide-50, an experiment operating since 2013 at LNGS, designed to detect WIMP interactions in liquid argon. In the process of WIMP searches, DarkSide-50 has achieved two fundamental milestones: the rejection of electron recoil background with a power of ~10 7 , using the pulse shape discrimination technique, and the measurement of the residual 39 Ar contamination in underground argon, ~3 orders of magnitude lower with respect to atmospheric argon. These results rely on the accurate simulation of the detector response to the liquid argon scintillation, its ionization, and electron-ion recombination processes. This work provides a complete overview of the DarkSide Monte Carlo and of its performance, with a particular focus on PARIS, the custom-made liquid argon response model.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Simulating image coaddition with the Nancy Grace Roman Space Telescope – II. Analysis of the simulated images and implications for weak lensing

One challenge for applying current weak lensing analysis tools to the Nancy Grace Roman Space Telescope is that individual images will be undersampled. Our companion paper presented an initial application of Imcom – an algorithm that builds an optimal mapping from input to output pixels to reconstruct a fully sampled combined image – on the Roman image simulations. In this paper, we measure the output noise power spectra, identify the sources of the major features in the power spectra, and show that simple analytic models that ignore sampling effects underestimate the power spectra of the coadded noise images. We compute the moments of both idealized injected stars and fully simulated stars in the coadded images, and their one- and two-point statistics. We show that the idealized injected stars have root-mean-square ellipticity errors (1–6) × 10 -4 per component depending on the band; the correlation functions are ≥2 orders of magnitude below requirements, indicating that the image combination step itself is using a small fraction of the overall Roman second moment error budget, although the fourth moments are larger and warrant further investigation. The stars in the simulated sky images, which include blending and chromaticity effects, have correlation functions near the requirement level (and below the requirement level in a wide-band image constructed by stacking all four filters). We evaluate the noise-induced biases in the ellipticities of injected stars, and explain the resulting trends with an analytical model. We conclude by enumerating the next steps in developing an image coaddition pipeline for Roman.

79 ASTRONOMY AND ASTROPHYSICS↗

Interpretation of proton radiography experiments of hohlraums with three-dimensional simulations [Interpretation of proton radiography experiments of hohlraums with 3D simulations]

Proton radiography experiments of laser-irradiated hohlraums performed at the OMEGA laser facility are analyzed using three-dimensional (3D) hydrodynamic simulations coupled to a proton trajectography package. Experiments with three different laser irradiation patterns were performed, and each produced a distinct proton image. By comparing these results with synthetic proton images obtained by sending protons through plasma profiles in the hohlraum obtained from 3D radiation hydrodynamic simulations, it is found that the simulated images agree favorably with the experimental images when electric fields, due to the electron pressure gradients that arise from 3D structures occurring during plasma expansion, are included. Furthermore, these comparisons provide quantitative estimates of the electric field present inside the hohlraums.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Time-dependent wave front propagation simulation of a hard x-ray split-and-delay unit: Towards a measurement of the temporal coherence properties of x-ray free electron lasers

For the European x-ray free electron laser (XFEL) a split-and-delay unit based on geometrical wavefront beam splitting and multilayer mirrors is built which covers the range of photon energies from 5 keV up to 20 keV. Maximum delays between Δτ = ±2.5 ps at hν=20 keV and up to Δτ = ±23 ps at hν = 5 keV will be possible. Time-dependent wave-optics simulations have been performed by means of Synchrotron Radiation Workshop software for XFEL pulses at hν = 5 keV. The XFEL radiation was simulated using results of time-dependent simulations applying the self-amplified spontaneous emission code FAST. Main features of the optical layout, including diffraction on the beam splitter edge and optics imperfections measured with a nanometer optic component measuring machine slope measuring profiler, were taken into account. The impact of these effects on the characterization of the temporal properties of XFEL pulses is analyzed. An approach based on fast Fourier transformation allows for the evaluation of the temporal coherence despite large wavefront distortions caused by the optics imperfections. In this manner, the fringes resulting from time-dependent two-beam interference can be filtered and evaluated yielding a coherence time of τ c = 0.187 fs (HWHM) for real, nonperfect mirrors, while for ideal mirrors a coherence time of τ c = 0.191 fs (HWHM) is expected.

36 MATERIALS SCIENCE↗

Hardware-based Advanced Electromagnetic Transient Simulation for A Large-Scale PV Plant in Real Time Digital Simulator

Power electronics-based resources, such as high-voltage direct current (HVdc) substations, photovoltaic (PV) plants, wind plants, electric vehicle charging stations, and energy storage systems, are increasingly being integrated within the power grid. Recently, multiple reports have emphasized the necessity for high-fidelity electromagnetic transient (EMT) simulations of these large-scale power electronics-based resources to accurately understand their behavior in power grids. However, performing hardware-based EMT simulations with high-fidelity models for such large power electronics systems is challenging due to the small time-step requirements and the involvement of a large number of states. This paper presents the implementation of hardware-based high-fidelity EMT dynamic model of a large-scale PV plant, accomplished through custom model development using the specific-C language in real-time digital simulator hardware (RTDS) and software (RSCAD).

Choi, Jongchan↗

In-Situ Feature Extraction of Large Scale Combustion Simulations Using Segmented Merge Trees, In: SC '14: Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis

The ever increasing amount of data generated by scientific simulations coupled with system I/O constraints are fueling a need for in-situ analysis techniques. Of particular interest are approaches that produce reduced data representations while maintaining the ability to redefine, extract, and study features in a post-process to obtain scientific insights. This paper presents two variants of in-situ feature extraction techniques using segmented merge trees, which encode a wide range of threshold based features. The first approach is a fast, low communication cost technique that generates an exact solution but has limited scalability. The second is a scalable, local approximation that nevertheless is guaranteed to correctly extract all features up to a predefined size. We demonstrate both variants using some of the largest combustion simulations available on leadership class supercomputers. Our approach allows state-of-the-art, feature-based analysis to be performed in-situ at significantly higher frequency than currently possible and with negligible impact on the overall simulation runtime.

Computer Science↗

Modeling of Cube Array Roughness: RANS, Large Eddy Simulation, and Direct Numerical Simulation

Abstract Flow over arrays of cubes is an extensively studied model problem for rough wall turbulent boundary layers. While considerable research has been performed in computationally investigating these topologies using direct numerical simulation (DNS) and large eddy simulation (LES), the ability of sublayer-resolved Reynolds-averaged Navier–Stokes (RANS) to predict the bulk flow phenomena of these systems is relatively unexplored, especially at low and high packing densities. Here, RANS simulations are conducted on six different packing densities of cubes in aligned and staggered configurations. The packing densities investigated span from what would classically be defined as isolated, up to those in the d-type roughness regime, filling in the gap in the present literature. Three different sublayer-resolved turbulence closure models were tested for each case: a low Reynolds number k–ϵ model, the Menter k–ω SST model, and a full Reynolds stress model. Comparisons of the velocity fields, secondary flow features, and drag coefficients are made between the RANS results and existing LES and DNS results. There is a significant degree of variability in the performance of the various RANS models across all comparison metrics. However, the Reynolds stress model demonstrated the best accuracy in terms of the mean velocity profile as well as drag partition across the range of packing densities.

Engineering↗

Perfect crystal propagator for physical optics simulations with Synchrotron Radiation Workshop

Until now, a treatment of dynamical diffraction from perfect crystals has been missing in the "Synchrotron Radiation Workshop" (SRW) wavefront propagation computer code despite the widespread use of crystals on X-ray synchrotron beamlines. Now a special propagator" module for calculating dynamical diffraction from a perfect crystal in the Bragg case has been written in C++, integrated into the SRW C/C++ library and made available for simulations using the Python interface of SRW. The propagator performs local processing of the frequency-domain electric field in the angular representation. A 2-D Fast Fourier Transform is used for changing the field representation from/to the coordinate representation before and after applying the crystal propagator. This ensures seamless integration of the new propagator with the existing functionalities of the SRW package, allows compatibility with existing propagators for other optical elements, and enables the simulation of complex beamlines transporting partially coherent X-rays. The code has been benchmarked by comparison with predictions made by plane-wave and spherical-wave dynamical diffraction theory. Test simulations for a selection of X-ray synchrotron beamlines are also shown.

47 OTHER INSTRUMENTATION↗

Crystallized Silicon Nanostructures - Experimental Characterization and Atomistic Simulations

We have synthesized silicon nanocrystalline structures from thermal annealing of thin film amorphous silicon-based multilayers. The annealing procedure that was carried out in vacuum at temperatures up to 1100 °C is integrated in a X-ray diffraction (XRD) setup for real-time monitoring of the formation phases of the nanostructures. The microstructure of the crystallized films is investigated through experimental measurements combined with atomistic simulations of realistic nanocrystalline silicon (nc-Si) models. The multilayers consisting of uniformly alternating thicknesses of hydrogenated amorphous silicon and silicon oxide (SiO2) were deposited by plasma enhanced chemical vapor deposition on crystalline silicon and Corning glass substrates. The crystallized structure consisting of nc-Si structures embedded in an amorphous matrix were further characterized through XRD, Raman spectroscopy, and Fourier transform infrared measurements. We are able to show the different stages of nanostructure formation and how the sizes and the crystallized mass fraction can be controlled in our experimental synthesis. The crystallized silicon structures with large crystalline filling fractions exceeding 50% have been simulated with a robust classical molecular dynamics technique. The crystalline filling fractions and structural order of nc-Si obtained from this simulation are compared with our Raman and XRD measurements.

36 MATERIALS SCIENCE↗

Preparing an Incompressible-Flow Fluid Dynamics Code for Exascale-Class Wind Energy Simulations

The U.S. Department of Energy has identified exascale-class wind farm simulation as critical to wind energy scientific discovery. A primary objective of the ExaWind project is to build high-performance, predictive computational fluid dynamics (CFD) tools that satisfy these modeling needs. GPU accelerators will serve as the computational thoroughbreds of next-generation, exascale-class supercomputers. Here, we report on our efforts in preparing the ExaWind unstructured mesh solver, Nalu-Wind, for exascale-class machines. For computing at this scale, a simple port of the incompressible-flow algorithms to GPUs is insufficient. To achieve high performance, one needs novel algorithms that are application aware, memory efficient, and optimized for the latest-generation GPU devices. The result of our efforts are unstructured-mesh simulations of wind turbines that can effectively leverage thousands of GPUs. In particular, we demonstrate a first-of-its-kind, incompressible-flow simulation using Algebraic Multigrid solvers that strong scales to more than 4000 GPUs on the Summit supercomputer.

algebraic multigrid↗

ParaEMT (NLR's (formerly known as NREL's) EMT Simulator: An Open Source, and HPC-Compatible Large-Scale Power System Electro-Magnetic Transient (EMT) Simulator) [SWR-22-16]

This open-source tool aims to provide a powerful and flexible platform for simulating electromagnetic transients (EMT) of large-scale inverter-based resource (IBR)-dominated power systems. The purpose of making it open source is to foster a community of EMT simulation and assist in exploring new EMT algorithms and applying advanced computational techniques to EMT simulation.

Tan, Jin↗

The Large-Eddy Simulation (LES) Atmospheric Radiation Measurement (ARM) Symbiotic Simulation and Observation (LASSO) Activity for Continental Shallow Convection

The U.S. Department of Energy’s Atmospheric Radiation Measurement (ARM) user facility recently initiated the Large-Eddy Simulation (LES) ARM Symbiotic Simulation and Observation (LASSO) activity focused on shallow convection at ARM’s Southern Great Plains (SGP) atmospheric observatory in Oklahoma. LASSO is designed to overcome an oft-shared difficulty of bridging the gap from point-based measurements to scales relevant for model parameterization development, and it provides an approach to add value to observations through modeling. LASSO is envisioned to be useful to modelers, theoreticians, and observationalists needing information relevant to cloud processes. LASSO does so by combining a suite of observations, LES inputs and outputs, diagnostics, and skill scores into data bundles that are freely available, and by simplifying user access to the data to speed scientific inquiry. The combination of relevant observations with observationally constrained LES output provides detail that gives context to the observations by showing physically consistent connections between processes based on the simulated state. A unique approach for LASSO is the generation of a library of cases for days with shallow convection combined with an ensemble of LES for each case. The library enables researchers to move beyond the single-case-study approach typical of LES research. The ensemble members are produced using a selection of different large-scale forcing sources and spatial scales. Since large-scale forcing is one of the most uncertain aspects of generating the LES, the ensemble informs users about potential uncertainty for each date and increases the probability of having an accurate forcing for each case.

54 ENVIRONMENTAL SCIENCES↗