Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulators”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

India National Gas Hydrate Program Expedition 02 Summary of Scientific Results: Numerical Simulation of Reservoir Response to Depressurization

The India National Gas Hydrate Program Expedition 02 (NGHP-02) discovered gas hydrate at high saturation in sand reservoirs at several sites in the deepwater Bay of Bengal. To assess the potential response of those deposits to scientific depressurization experiments, comprehensive geologic models were constructed to enable numerical simulation for two sites. Both sites (NGHP-02-09 and NGHP-02-16) feature thick sequences of thinly-interbedded reservoir and non-reservoir facies at sub-seafloor depths less than 300m and sub-sea depths of 2400m or more. These settings pose significant challenges to current modeling capabilities. First, the thinly-interbedded reservoir architecture complicates the determination of basic reservoir parameters from both log and core data due to measurement resolution issues. Secondly, the fine scale variation in sediment properties imparts great contrasts in key parameters over very short distances, creating high gradients at multiple scales and varying orientations that necessitate careful design of high-definition simulation grids. Thirdly, the deposits include internal sources of water, as well as a range of complex boundary conditions, including variable permeability within the overlying mud-rich “seals,” that complicate reservoir depressurization. Lastly, because of the unique combination of great water depth and relatively shallow sub-seafloor depth: models designed to maximize the dissociation rate impose large pressure drawdowns on relatively low-strength sediments. This condition renders the proper evaluation and integration of the geomechanical response to hydrate dissociation critical. In this report, we review the history of gas hydrate reservoir simulation, discuss methods for creating geologic input models, and summarize the key findings and implications of the collaborative NGHP-02 numerical simulation effort. Together, the studies confirm the viability of the modeled accumulations for scientific testing and identify key challenges related to the selection of specific test sites and the design of test wells.

Boswell, Ray↗
DOE OSTI2015

Gyrokinetic particle simulation of microturbulence for general magnetic geometry and experimental profiles

Developments in gyrokinetic particle simulation enable the gyrokinetic toroidal code (GTC) to simulate turbulent transport in tokamaks with realistic equilibrium profiles and plasma geometry, which is an essential step in the code–experiment validation process. These novel developments include numerical equilibrium representation using B-splines, a new Poisson solver based on finite difference using field-aligned mesh and magnetic flux coordinates, a new zonal flow solver for general geometry, and improvements on the conventional four-point gyroaverage with nonuniform background marker loading. The gyrokinetic Poisson equation is solved in the perpendicular plane instead of the poloidal plane. Exploiting these new features, GTC is able to simulate a typical DIII-D discharge with experimental magnetic geometry and profiles. The simulated turbulent heat diffusivity and its radial profile reflects good agreement with other gyrokinetic codes. The newly developed nonuniform loading method provides a modified radial transport profile to that of the conventional uniform loading method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Kinetic electron and ion instability of the lunar wake simulated at physical mass ratio

The solar wind wake behind the moon is studied with 1D electrostatic particle-in-cell (PIC) simulations using a physical ion to electron mass ratio (unlike prior investigations); the simulations also apply more generally to supersonic flow of dense magnetized plasma past non-magnetic objects. In this work, a hybrid electrostatic Boltzmann electron treatment is first used to investigate the ion stability in the absence of kinetic electron effects, showing that the ions are two-stream unstable for downstream wake distances (in lunar radii) greater than about three times the solar wind Mach number. Simulations with PIC electrons are then used to show that kinetic electron effects can lead to disruption of the ion beams at least three times closer to the moon than in the hybrid simulations. This disruption occurs as the result of a novel wake phenomenon: the non-linear growth of electron holes spawned from a narrow dimple in the electron velocity distribution. Most of the holes arising from the dimple are small and quickly leave the wake, approximately following the unperturbed electron phase-space trajectories, but some holes originating near the center of the wake remain and grow large enough to trigger disruption of the ion beams. Lastly, non-linear kinetic-electron effects are therefore essential to a comprehensive understanding of the 1D electrostatic stability of such wakes, and possible observational signatures in ARTEMIS data from the lunar wake are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Gyrokinetic particle simulation of a field reversed configuration

Gyrokinetic particle simulation of the field-reversed configuration (FRC) has been developed using the gyrokinetic toroidal code (GTC). The magnetohydrodynamic equilibrium is mapped from cylindrical coordinates to Boozer coordinates for the FRC core and scrape-off layer (SOL), respectively. A field-aligned mesh is constructed for solving self-consistent electric fields using a semi-spectral solver in a partial torus FRC geometry. This new simulation capability has been successfully verified and driftwave instability in the FRC has been studied using the gyrokinetic simulation for the first time. Initial GTC simulations find that in the FRC core, the ion-scale driftwave is stabilized by the large ion gyroradius. In the SOL, the driftwave is unstable on both ion and electron scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Simulation of stimulated Brillouin scattering and stimulated Raman scattering in shock ignition

We study stimulated Brillouin scattering (SBS) and stimulated Raman scattering (SRS) in shock ignition by comparing fluid and particle-in-cell (PIC) simulations. Under typical parameters for the OMEGA experiments [Theobald et al., Phys. Plasmas 19, 102706 (2012)], a series of 1D fluid simulations with laser intensities ranging between 2 × 10{sup 15} and 2 × 10{sup 16 }W/cm{sup 2} finds that SBS is the dominant instability, which increases significantly with the incident intensity. Strong pump depletion caused by SBS and SRS limits the transmitted intensity at the 0.17n{sub c} to be less than 3.5 × 10{sup 15 }W/cm{sup 2}. The PIC simulations show similar physics but with higher saturation levels for SBS and SRS convective modes and stronger pump depletion due to higher seed levels for the electromagnetic fields in PIC codes. Plasma flow profiles are found to be important in proper modeling of SBS and limiting its reflectivity in both the fluid and PIC simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Gyrokinetic particle simulation of beta-induced Alfven-acoustic eigenmode

The beta-induced Alfven-acoustic eigenmode (BAAE) in toroidal plasmas is verified and studied by global gyrokinetic particle simulations. When ion temperature is much lower than electron temperature, the existence of the weakly damped BAAE is verified in the simulations using initial perturbation, antenna excitation, and energetic particle excitation, respectively. When the ion temperature is comparable to the electron temperature, the unstable BAAE can be excited by realistic energetic particle density gradient, even though the stable BAAE (in the absence of energetic particles) is heavily damped by the thermal ions. In the simulations with reversed magnetic shear, BAAE frequency sweeping is observed and poloidal mode structure has a triangle shape with a poloidal direction similar to that observed in tokamak experiments. Here, the triangle shape changes the poloidal direction, and no frequency sweeping is found in the simulations with normal magnetic shear.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Simulations of the infrared, Raman, and 2D-IR photon echo spectra of water in nanoscale silica pores

Vibrational spectroscopy is frequently used to characterize nanoconfined liquids and probe the effect of the confining framework on the liquid structure and dynamics relative to the corresponding bulk fluid. However, it is still unclear what molecular-level information can be obtained from such measurements. In this paper, we address this question by using molecular dynamics (MD) simulations to reproduce the linear infrared (IR), Raman, and two-dimensional IR (2D-IR) photon echo spectra for water confined within hydrophilic (hydroxyl-terminated) silica mesopores. To simplify the spectra the OH stretching region of isotopically dilute HOD in D{sub 2}O is considered. An empirical mapping approach is used to obtain the OH vibrational frequencies, transition dipoles, and transition polarizabilities from the MD simulations. The simulated linear IR and Raman spectra are in good general agreement with measured spectra of water in mesoporous silica reported in the literature. The key effect of confinement on the water spectrum is a vibrational blueshift for OH groups that are closest to the pore interface. The blueshift can be attributed to the weaker hydrogen bonds (H-bonds) formed between the OH groups and silica oxygen acceptors. Non-Condon effects greatly diminish the contribution of these OH moieties to the linear IR spectrum, but these weaker H-bonds are readily apparent in the Raman spectrum. The 2D-IR spectra have not yet been measured and thus the present results represent a prediction. The simulated spectra indicates that it should be possible to probe the slower spectral diffusion of confined water compared to the bulk liquid by analysis of the 2D-IR spectra.

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

High-Fidelity Simulation of Flow-Induced Vibrations in Helical Steam Generators for Small Modular Reactors

Flow-induced vibration (FIV) is a widespread problem in energy systems as they rely on fluid movement for energy conversion. Vibrating structures may be damaged as fatigue or wear occur. Given the importance of reliable components in the nuclear industry, FIV has long been a major concern in the safety and operation of nuclear reactors. In particular, nuclear fuel rods and steam generators have been known to suffer from FIV and related failures. In this paper we discuss the use of the computational fluid dynamics code Nek5000 coupled to the structural code Diablo to simulate the flow in helical coil heat exchangers and associated FIV. In particular, one-way coupled calculations are performed, where pressure and tractions data are loaded into the structural model. The main focus of this paper is on validation of this capability. Fluid-only Nek5000 large eddy simulations are first compared against dedicated high-resolution experiments. Then, one-way coupled calculations are performed with Nek5000 and Diablo for two data sets that provide FIV data for validation. These calculations were aimed at simulating available legacy FIV experiments in helical steam generators in the turbulent buffeting regime. In this regime one-way coupling is judged sufficient since the pressure loads do not cause substantial displacements. It is also the most common source of vibration in helical steam generators at the low flows expected in integral pressurized water reactors. We discuss validation of two-way coupled experiments and benchmarks toward the simulation of fluid elastic instability. We briefly discuss the application of these methods to grid-to-rod fretting.

Flow-Induced Vibration↗
DOE OSTI2014

Fully coupled simulation of cosmic reionization. II. Recombinations, clumping factors, and the photon budget for reionization

We use a fully coupled cosmological simulation including dark matter dynamics, multispecies hydrodynamics, nonequilibrium chemical ionization, flux-limited diffusion radiation transport, and a parameterized model of star formation and feedback (thermal and radiative) to investigate the epoch of hydrogen reionization in detail. In this paper, the first of several application papers, we investigate the mechanics of reionization from stellar sources forming in high-z galaxies, the utility of various formulations for the gas clumping factor on accurately estimating the effective recombination time in the intergalactic medium (IGM), and the photon budget required to achieve reionization. We also test the accuracy of the static and time-dependent models of Madau et al. as predictors of reionization completion/maintenance. We simulate a WMAP7 ΛCDM cosmological model in a 20 comoving Mpc cube, resolved with 800{sup 3} uniform fluid cells and dark matter particles. By tuning our star formation recipe to approximately match the observed high-redshift star formation rate density and galaxy luminosity function, we have created a fully coupled radiation hydrodynamical realization of hydrogen reionization, which begins to ionize at z ≈ 10 and is completed at z ≈ 5.8 without further tuning. We find that roughly two ionizing photons per H atom are required to convert the neutral IGM to a highly ionized state. After reionization concludes, we find that the quantity n-dot{sub ion}×(1 Gyr)/n{sub H} is ∼9 at z = 5, in rough agreement with measurements of the ionizing emissivity by Becker and Bolton. The complicated events during reionization that lead to this number can be generally described as inside-out, but in reality, the narrative depends on the level of ionization of the gas one attributes as being ionized. We find that the formula for the ionizing photon production rate needed to maintain the IGM in an ionized state derived by Madau et al. should not be used to predict the epoch of reionization completion because it ignores history-dependent terms in the global ionization balance which are not ignorable. We find that the time-dependent model for the ionized volume fraction Q{sub HII} is more predictive, but overestimates the redshift of reionization completion z{sub reion} by Δz ≈ 1. We propose a revised formulation of the time-dependent model that agrees with our simulation to a few percent accuracy. Finally, we use our simulation data to measure the absorption of ionizing radiation due to circumgalactic gas resolved on our mesh and find f-bar{sub esc}(CGM)≈0.7.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

SIMULATIONS OF DISK GALAXIES WITH COSMIC RAY DRIVEN GALACTIC WINDS

We present results from high-resolution hydrodynamic simulations of isolated Small Magellanic Cloud (SMC)- and Milky-Way-sized galaxies that include a model for feedback from galactic cosmic rays (CRs). We find that CRs are naturally able to drive winds with mass loading factors of up to ∼10 in dwarf systems. The scaling of the mass loading factor with circular velocity between the two simulated systems is consistent with η∝v{sub circ}{sup 1-2} required to reproduce the faint end of the galaxy luminosity function. In addition, simulations with CR feedback reproduce both the normalization and the slope of the observed trend of wind velocity with galaxy circular velocity. We find that winds in simulations with CR feedback exhibit qualitatively different properties compared to supernova-driven winds, where most of acceleration happens violently in situ near star forming sites. The CR-driven winds are accelerated gently by the large-scale pressure gradient established by CRs diffusing from the star-forming galaxy disk out into the halo. The CR-driven winds also exhibit much cooler temperatures and, in the SMC-sized system, warm (T ∼ 10{sup 4} K) gas dominates the outflow. The prevalence of warm gas in such outflows may provide a clue as to the origin of ubiquitous warm gas in the gaseous halos of galaxies detected via absorption lines in quasar spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

THREE-DIMENSIONAL CORE-COLLAPSE SUPERNOVA SIMULATED USING A 15 M{sub ⊙} PROGENITOR

We have performed ab initio neutrino radiation hydrodynamics simulations in three and two spatial dimensions (3D and 2D) of core-collapse supernovae from the same 15 M{sub ☉} progenitor through 440 ms after core bounce. Both 3D and 2D models achieve explosions; however, the onset of explosion (shock revival) is delayed by ∼100 ms in 3D relative to the 2D counterpart and the growth of the diagnostic explosion energy is slower. This is consistent with previously reported 3D simulations utilizing iron-core progenitors with dense mantles. In the ∼100 ms before the onset of explosion, diagnostics of neutrino heating and turbulent kinetic energy favor earlier explosion in 2D. During the delay, the angular scale of convective plumes reaching the shock surface grows and explosion in 3D is ultimately lead by a single, large-angle plume, giving the expanding shock a directional orientation not dissimilar from those imposed by axial symmetry in 2D simulations. We posit that shock revival and explosion in the 3D simulation may be delayed until sufficiently large plumes form, whereas such plumes form more rapidly in 2D, permitting earlier explosions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Dark Energy Survey Year 3 results: Cosmology from combined galaxy clustering and lensing validation on cosmological simulations

In this work, we present a validation of the Dark Energy Survey Year 3 (DES Y3) 3×2-point analysis choices by testing them on Buzzard2.0, a new suite of cosmological simulations that is tailored for the testing and validation of combined galaxy clustering and weak-lensing analyses. We show that the buzzard2.0 simulations accurately reproduce many important aspects of the DES Y3 data, including photometric redshift and magnitude distributions, and the relevant set of two-point clustering and weak-lensing statistics. We then show that our model for the 3×2-point data vector is accurate enough to recover the true cosmology in simulated surveys assuming the true redshift distributions for our source and lens samples, demonstrating robustness to uncertainties in the modeling of the nonlinear matter power spectrum, nonlinear galaxy bias, and higher-order lensing corrections. Additionally, we demonstrate for the first time that our photometric redshift calibration methodology, including information from photometry, spectroscopy, clustering cross-correlations, and galaxy–galaxy lensing ratios, is accurate enough to recover the true cosmology in simulated surveys in the presence of realistic photometric redshift uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Finite Element Simulation and X-Ray Microdiffraction Study of Strain Partitioning in a Layered Nanocomposite

The depth-dependent strain partitioning across the interfaces in the growth direction of the NiAl/Cr(Mo) nanocomposite between the Cr and NiAl lamellae was directly measured experimentally and simulated using a finite element method (FEM). Depth-resolved X-ray microdiffraction demonstrated that in the as-grown state both Cr and NiAl lamellae grow along the 111 direction with the formation of as-grown distinct residual ~0.16% compressive strains for Cr lamellae and ~0.05% tensile strains for NiAl lamellae. Three-dimensional simulations were carried out using an implicit FEM. First simulation was designed to study residual strains in the composite due to cooling resulting in formation of crystals. Strains in the growth direction were computed and compared to those obtained from the microdiffraction experiments. Second simulation was conducted to understand the combined strains resulting from cooling and mechanical indentation of the composite. Numerical results in the growth direction of crystal were compared to experimental results confirming the experimentally observed trends.

Barabash, R. I.↗
DOE OSTI2017

Self-consistent Modeling of Reionization in Cosmological Hydrodynamical Simulations

The ultraviolet background (UVB) emitted by quasars and galaxies governs the ionization and thermal state of the intergalactic medium (IGM), regulates the formation of high-redshift galaxies, and is thus a key quantity for modeling cosmic reionization. The vast majority of cosmological hydrodynamical simulations implement the UVB via a set of spatially uniform photoionization and photoheating rates derived from UVB synthesis models. We show that simulations using canonical UVB rates reionize and, perhaps more importantly, spuriously heat the IGM, much earlier (z∼15) than they should. This problem arises because at z>6, where observational constraints are nonexistent, the UVB amplitude is far too high. We introduce a new methodology to remedy this issue, and we generate self-consistent photoionization and photoheating rates to model any chosen reionization history. Following this approach, we run a suite of hydrodynamical simulations of different reionization scenarios and explore the impact of the timing of reionization and its concomitant heat injection on the thermal state of the IGM. We present a comprehensive study of the pressure smoothing scale of IGM gas, illustrating its dependence on the details of both hydrogen and helium reionization, and argue that it plays a fundamental role in interpreting Lyα forest statistics and the thermal evolution of the IGM. The premature IGM heating we have uncovered implies that previous work has likely dramatically overestimated the impact of photoionization feedback on galaxy formation, which sets the minimum halo mass able to form stars at high redshifts. We make our new UVB photoionization and photoheating rates publicly available for use in future simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Distributed Macroscopic Traffic Simulation with Open Traffic Models

This paper presents OTM-MPI, an extension of the Open Traffic Models platform (OTM) for running macroscopic traffic simulations in high-performance computing environments. OTM-MPI represents the first open-source, distributed-memory, macroscopic simulation model developed for modern high performance parallel machines and large networks. Macroscopic simulations are appropriate for studying regional traffic scenarios when aggregate trends are of interest, rather than individual vehicle traces. They are also appropriate for studying the routing behavior of classes of vehicles, such as app-informed vehicles. The network partitioning was performed with METIS. Inter-process communication was done with MPI (message-passing interface). Results are provided for two networks: one realistic network which was obtained from Open Street Maps for Chattanooga, TN, and another larger synthetic grid network. The software recorded a speedups of 198x using 256 cores for Chattanooga, and 475x with 1,024 cores for the synthetic network.

macro-scopic traffic simulation↗
DOE OSTI2018

Simulated building energy demand biases resulting from the use of representative weather stations

Numerical building models are typically forced with weather data from a limited number of “representative cities” or weather stations representing different climate regions. The use of representative weather stations reduces computational costs, but often fails to capture spatial heterogeneity in weather that may be important for simulations aimed at understanding how building stocks respond to a changing climate. We quantify the potential reduction in bias from using an increasing number of weather stations over the western U.S. The approach is based on deriving temperature and load time series using incrementally more weather stations, ranging from 8 to roughly 150, to capture weather across different seasons. Using 8 stations, one from each climate zone, across the western U.S. results in an average absolute summertime temperature bias of 7.2°F with respect to a spatially-resolved gridded dataset. The mean absolute bias drops to 2.8°F using all available weather stations. Temperature biases of this magnitude could translate to absolute summertime mean simulated load biases as high as 13.8%, a significant error for capacity expansion planners who may use these types of simulations. Increasing the size of the domain over which biases are calculated reduces their magnitude as positive and negative biases may cancel out. Using 8 representative weather stations can lead to a 20-40% overestimation of peak building loads during both summer and winter. Using weather stations close to population centers reduces both mean and peak load biases. This approach could be used by others designing aggregate building simulations to understand the sensitivity to their choice of weather stations used to drive the models.

Burleyson, Casey D.↗
DOE OSTI2021

Higher-order particle representation for particle-in-cell simulations

In this paper we present an alternative approach to the representation of simulation particles for unstructured electrostatic and electromagnetic PIC simulations. In our modified PIC algorithm we represent particles as having a smooth shape function limited by some specified finite radius, r 0 . A unique feature of our approach is the representation of this shape by surrounding simulation particles with a set of virtual particles with delta shape, with fixed offsets and weights derived from Gaussian quadrature rules and the value of r 0 . As the virtual particles are purely computational, they provide the additional benefit of increasing the arithmetic intensity of traditionally memory bound particle kernels. The modified algorithm is implemented within Sandia National Laboratories' unstructured EMPIRE-PIC code, for electrostatic and electromagnetic simulations, using periodic boundary conditions. We show results for a representative set of benchmark problems, including electron orbit, a transverse electromagnetic wave propagating through a plasma, numerical heating, and a plasma slab expansion. In this work, good error reduction across all of the chosen problems is achieved as the particles are made progressively smoother, with the optimal particle radius appearing to be problem-dependent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗