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

The Aemulus Project. I. Numerical Simulations for Precision Cosmology

The rapidly growing statistical precision of galaxy surveys has led to a need for ever more precise predictions of the observables used to constrain cosmological and galaxy formation models. The primary avenue through which such predictions will be obtained is suites of numerical simulations. These simulations must span the relevant model parameter spaces, be large enough to obtain the precision demanded by upcoming data, and be thoroughly validated in order to ensure accuracy. In this paper, we present one such suite of simulations, forming the basis for the AEMULUS Project, a collaboration devoted to precision emulation of galaxy survey observables. We have run a set of 75 (1.05 h(-1) Gpc)(3) simulations with mass resolution and force softening of 3.51 x 10(10) (Omega(m)/0.3) h(-1) M-circle dot and 20 h(-1) kpc, respectively, in 47 different wCDM cosmologies spanning the range of parameter space allowed by the combination of recent cosmic microwave background, baryon acoustic oscillation, and Type Ia supernova results. We present convergence tests of several observables including spherical overdensity halo mass functions, galaxy projected correlation functions, galaxy clustering in redshift space, and matter and halo correlation functions and power spectra. We show that these statistics are converged to 1% (2%) or to the sample variance of the statistic, whichever is larger, for halos with more than 500 (200) particles, respectively, and scales of r > 200 h(-1) kpc. in real space or k similar to 3 h Mpc(-1). in harmonic space for z <= 1. We find that the dominant source of uncertainty comes from varying the particle loading of the simulations. This leads to large systematic errors for statistics using halos with fewer than 200 particles and scales smaller than k similar to 4 h Mpc(-1). We provide the halo catalogs and snapshots detailed in this work to the community at. https://AemulusProject.github.io.

large-scale structure of universe↗
DOE OSTI2017

A Comparison of Grid-based and SPH Binary Mass-transfer and Merger Simulations

There is currently a great amount of interest in the outcomes and astrophysical implications of mergers of double degenerate binaries. In a commonly adopted approximation, the components of such binaries are represented by polytropes with an index of n = 3/2. We present detailed comparisons of stellar mass-transfer and merger simulations of polytropic binaries that have been carried out using two very different numerical algorithms—a finite-volume “grid” code and a smoothed-particle hydrodynamics (SPH) code. We find that there is agreement in both the ultimate outcomes of the evolutions and the intermediate stages if the initial conditions for each code are chosen to match as closely as possible. We find that even with closely matching initial setups, the time it takes to reach a concordant evolution differs between the two codes because the initial depth of contact cannot be matched exactly. There is a general tendency for SPH to yield higher mass transfer rates and faster evolution to the final outcome. We also present comparisons of simulations calculated from two different energy equations: in one series, we assume a polytropic equation of state and in the other series an ideal gas equation of state. In the latter series of simulations, an atmosphere forms around the accretor, which can exchange angular momentum and cause a more rapid loss of orbital angular momentum. In the simulations presented here, the effect of the ideal equation of state is to de-stabilize the binary in both SPH and grid simulations, but the effect is more pronounced in the grid code.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Application of Artificial Neural Networks to Reliable Nuclear Data for Nonproliferation Modeling and Simulation

Detection and identification of special nuclear materials (SNMs) are an essential part of the US nonproliferation effort. Modern cutting-edge SNM detection methodologies rely more and more on modeling and simulation techniques. Experiments with radiological samples in realistic configurations, is the ultimate tool that establishes the minimum detection limits of SNMs in a host of different geometries. Modern modeling and simulation approaches have the potential to significantly reduce the number of experiments with radioactive sources needed to determine these detection limits and reduce the financial barrier to SNM detection. Unreliable nuclear data is one of the principal causes of uncertainty in modeling and simulating nuclear systems. In particular, nuclear cross sections introduce a significant uncertainty in the nuclear data. The goal of this research is to develop a methodology that will autonomously extract the correct nuclear resonance characteristics of experimental data in a reliable way, a task previously left to expert judgement. Accurate nuclear data will in turn allow contemporary modeling and simulation to become far more reliable, de-escalating the extent of experimental testing. Modeling and simulation techniques reduce the use and distribution of radiological sources, while at the same time increase the reliability of the currently used methods for the detection and identification of SNMs.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2017

The laboratory environmental algae pond simulator (LEAPS) photobioreactor: Validation using outdoor pond cultures of Chlorella sorokiniana and Nannochloropsis salina

A bench-scale photobioreactor system, the Laboratory Environmental Algae Pond Simulator (LEAPS), was designed and constructed to simulate outdoor pond cultivation for a wide range of geographical locations and seasons. The LEAPS consists of six well-mixed glass column photobioreactors submerged in a temperature controlled waterbath (-2 °C to >60 °C) and illuminated from above by a programmable multicolor LED lighting system (0 to 2,500 µmol/m 2 -sec). Measured incident light intensities and water temperatures deviated from the respective light and temperature setpoints on average only 2.3% and 0.9%, demonstrating accurate simulation of light and temperature conditions measured in outdoor ponds. In order to determine whether microalgae strains cultured in the LEAPS exhibit the same linear phase biomass productivity as in outdoor ponds, Chlorella sorokiniana and Nannochloropsis salina were cultured in the LEAPS bioreactors using light and temperature scripts measured previously in the respective outdoor pond studies. For Chlorella sorokiniana, the summer season biomass productivity in the LEAPS was 6.6% and 11.3% lower than in the respective outdoor ponds in Rimrock, Arizona, and Delhi, California; however, these differences were not statistically significant. For Nannochloropsis salina, the winter season biomass productivity in the LEAPS was 9.2% higher than in the respective outdoor ponds in Tucson, Arizona, a difference shown to be not statistically significant. Potential reasons for the positive or negative divergence in LEAPS performance, relative to outdoor ponds, are discussed. To demonstrate the utility of the LEAPS in predicting productivity, two other strains – Scenedesmus obliquus and Stichococcus minor – were evaluated using the summer season script for Rimrock, Arizona. For both strains, the productivity was around 11.6 g/m 2 -day at 25 cm culture depth. In conclusion, the LEAPS is an accurate pond simulator and thus offers a reliable, fast, and cost-effective way to screen microalgae strains and operating conditions for high biomass productivity and co-product yields.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Systematic Validation of Ensemble Cloud-Process Simulations using Polarimetric Radar Observations and Simulator over the NASA Wallops Flight Facility

The BiLateral Operational Storm-Scale Observation and Modeling (BLOSSOM) project was initiated to establish a long-term supersite to improve understanding of cloud physical states and processes as well as to support satellite and climate model programs over the WFF site via a bilateral approach of storm-scale observations and process modeling. This study highlights a noble systematic validation framework of the BLOSSOM ensemble cloud-process simulations through mixed-phase, light-rain, and deep-convective precipitation cases. The framework consists of creating a domain-shifted ensemble of large-scale forcing datasets, and configuring and performing cloud-process simulations with three different bulk microphysics schemes. Validation uses NASA S-band dual-POLarimetric radar (NPOL) observations in the form of statistical composites and skill scores via a polarimetric radar simulator and newly developed CfRad Data tool (CfRAD). While the simulations capture the overall structures of the reflectivity composites, polarimetric signals are still poorly simulated, mainly due to a lack of representation of ice microphysics diversity in shapes, orientation distributions, and their complex mixtures. Despite the limitation, this new ensemble-based validation framework demonstrates that 1) no particular forcing or microphysics scheme outperforms the rest and 2) the skill scores of coarse- and fine-resolution ensemble simulations with different domain-shifted forcing and microphysics schemes are highly correlated with each other with no clear improvement. On the other hand, this suggests that coarse-resolution ensemble simulations are relevant for selecting the best meteorological forcing and microphysics scheme before conducting computationally demanding large eddy simulations (LESs) in support of aircraft and satellite instrument development as well as cloud-precipitation-convection parameterizations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Refinement of the experimental dynamic structure factor for liquid para-hydrogen and ortho-deuterium using semi-classical quantum simulation

The dynamic structure factor of liquid para-hydrogen and ortho-deuterium in corresponding thermodynamic states (T = 20.0 K, n = 21.24 nm{sup −3}) and (T = 23.0 K, n = 24.61 nm{sup −3}), respectively, has been computed by both the Feynman-Kleinert linearized path-integral (FK-LPI) and Ring-Polymer Molecular Dynamics (RPMD) methods and compared with Inelastic X Ray Scattering spectra. The combined use of computational and experimental methods enabled us to reduce experimental uncertainties in the determination of the true sample spectrum. Furthermore, the refined experimental spectrum of para-hydrogen and ortho-deuterium is consistently reproduced by both FK-LPI and RPMD results at momentum transfers lower than 12.8 nm{sup −1}. At larger momentum transfers the FK-LPI results agree with experiment much better for ortho-deuterium than for para-hydrogen. More specifically we found that for k ∼ 20.0 nm{sup −1} para-hydrogen provides a test case for improved approximations to quantum dynamics.

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

Initial performances of first undulator-based hard x-ray beamlines of NSLS-II compared to simulations

Commissioning of the first X-ray beamlines of NSLS-II included detailed measurements of spectral and spatial distributions of the radiation at different locations of the beamlines, from front-ends to sample positions. Additionally, comparison of some of these measurement results with high-accuracy calculations of synchrotron (undulator) emission and wavefront propagation through X-ray transport optics, performed using SRW code, is presented.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

Simulated performance of the optical Thomson scattering diagnostic designed for the National Ignition Facility

An optical Thomson scattering diagnostic has been designed for the National Ignition Facility to characterize under-dense plasmas. We report on the design of the system and the expected performance for different target configurations. The diagnostic is designed to spatially and temporally resolve the Thomson scattered light from laser driven targets. The diagnostic will collect scattered light from a 50 × 50 × 200 μm volume. The optical design allows operation with different probe laser wavelengths. A deep-UV probe beam (λ{sub 0} = 210 nm) will be used to Thomson scatter from electron plasma densities of ∼5 × 10{sup 20} cm{sup −3} while a 3ω probe will be used for plasma densities of ∼1 × 10{sup 19} cm{sup −3}. The diagnostic package contains two spectrometers: the first to resolve Thomson scattering from ion acoustic wave fluctuations and the second to resolve scattering from electron plasma wave fluctuations. Expected signal levels relative to background will be presented for typical target configurations (hohlraums and a planar foil).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

A 2.5D Computational Method to Simulate Cylindrical Fluidized Beds

In this paper, the limitations of axisymmetric and Cartesian two-dimensional (2D) simulations of cylindrical gas-solid fluidized beds are discussed. A new method has been proposed to carry out pseudo-two-dimensional (2.5D) simulations of a cylindrical fluidized bed by appropriately combining computational domains of Cartesian 2D and axisymmetric simulations. The proposed method was implemented in the open-source code MFIX and applied to the simulation of a lab-scale bubbling fluidized bed with necessary sensitivity study. After a careful grid study to ensure the numerical results are grid independent, detailed comparisons of the flow hydrodynamics were presented against axisymmetric and Cartesian 2D simulations. Furthermore, the 2.5D simulation results have been compared to the three-dimensional (3D) simulation for evaluation. This new approach yields better agreement with the 3D simulation results than with axisymmetric and Cartesian 2D simulations.

Li, Tingwen↗
DOE OSTI2014

Benchmarking all-atom simulations using hydrogen exchange

Molecular dynamics simulations have recently become capable of observing multiple protein unfolding and refolding events in a single-millisecond–long trajectory. This major advance produces atomic-level information with nanosecond resolution, a feat unmatched by experimental methods. Such simulations are being extensively analyzed to assess their description of protein folding, yet the results remain difficult to validate experimentally. We apply a combination of hydrogen exchange, NMR, and other techniques to test the simulations with a resolution of single H-bonds. Several significant discrepancies between the simulations and experimental data were uncovered for regions of the energy surface outside of the native basin. This comparison yields suggestions for improving the force fields and provides a general method for experimentally validating folding simulations.

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

Simulation-based inference for parameter estimation of complex watershed simulators

High-resolution, spatially distributed process-based (PB) simulators are widely employed in the study of complex catchment processes and their responses to a changing climate. However, calibrating these PB simulators using observed data remains a significant challenge due to several persistent issues, including the following: (1) intractability stemming from the computational demands and complex responses of simulators, which renders infeasible calculation of the conditional probability of parameters and data, and (2) uncertainty stemming from the choice of simplified representations of complex natural hydrologic processes. Here, we demonstrate how simulation-based inference (SBI) can help address both of these challenges with respect to parameter estimation. SBI uses a learned mapping between the parameter space and observed data to estimate parameters for the generation of calibrated simulations. To demonstrate the potential of SBI in hydrologic modeling, we conduct a set of synthetic experiments to infer two common physical parameters – Manning's coefficient and hydraulic conductivity – using a representation of a snowmelt-dominated catchment in Colorado, USA. We introduce novel deep-learning (DL) components to the SBI approach, including an “emulator” as a surrogate for the PB simulator to rapidly explore parameter responses. We also employ a density-based neural network to represent the joint probability of parameters and data without strong assumptions about its functional form. While addressing intractability, we also show that, if the simulator does not represent the system under study well enough, SBI can yield unreliable parameter estimates. Approaches to adopting the SBI framework for cases in which multiple simulator(s) may be adequate are introduced using a performance-weighting approach. The synthetic experiments presented here test the performance of SBI, using the relationship between the surrogate and PB simulators as a proxy for the real case.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Parameter estimation for chaotic systems using a hybrid adaptive cuckoo search with simulated annealing algorithm

This paper introduces a novel hybrid optimization algorithm to establish the parameters of chaotic systems. In order to deal with the weaknesses of the traditional cuckoo search algorithm, the proposed adaptive cuckoo search with simulated annealing algorithm is presented, which incorporates the adaptive parameters adjusting operation and the simulated annealing operation in the cuckoo search algorithm. Normally, the parameters of the cuckoo search algorithm are kept constant that may result in decreasing the efficiency of the algorithm. For the purpose of balancing and enhancing the accuracy and convergence rate of the cuckoo search algorithm, the adaptive operation is presented to tune the parameters properly. Besides, the local search capability of cuckoo search algorithm is relatively weak that may decrease the quality of optimization. So the simulated annealing operation is merged into the cuckoo search algorithm to enhance the local search ability and improve the accuracy and reliability of the results. The functionality of the proposed hybrid algorithm is investigated through the Lorenz chaotic system under the noiseless and noise condition, respectively. The numerical results demonstrate that the method can estimate parameters efficiently and accurately in the noiseless and noise condition. Finally, the results are compared with the traditional cuckoo search algorithm, genetic algorithm, and particle swarm optimization algorithm. Simulation results demonstrate the effectiveness and superior performance of the proposed algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Graphics processing unit accelerated three-dimensional model for the simulation of pulsed low-temperature plasmas

A 3-dimensional particle-in-cell/Monte Carlo collision simulation that is fully implemented on a graphics processing unit (GPU) is described and used to determine low-temperature plasma characteristics at high reduced electric field, E/n, in nitrogen gas. Details of implementation on the GPU using the NVIDIA Compute Unified Device Architecture framework are discussed with respect to efficient code execution. The software is capable of tracking around 10 × 10{sup 6} particles with dynamic weighting and a total mesh size larger than 10{sup 8} cells. Verification of the simulation is performed by comparing the electron energy distribution function and plasma transport parameters to known Boltzmann Equation (BE) solvers. Under the assumption of a uniform electric field and neglecting the build-up of positive ion space charge, the simulation agrees well with the BE solvers. The model is utilized to calculate plasma characteristics of a pulsed, parallel plate discharge. A photoionization model provides the simulation with additional electrons after the initial seeded electron density has drifted towards the anode. Comparison of the performance benefits between the GPU-implementation versus a CPU-implementation is considered, and a speed-up factor of 13 for a 3D relaxation Poisson solver is obtained. Furthermore, a factor 60 speed-up is realized for parallelization of the electron processes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Development of a Wind Plant Large-Eddy Simulation with Measurement-Driven Atmospheric Inflow: Preprint

This paper details the development of an aeroelastic wind plant model with large-eddy simulation (LES). The chosen LES solver is the Simulator for Wind Farm Applications (SOWFA) based on the OpenFOAM framework, coupled to NREL's comprehensive aeroelastic analysis tool, FAST. An atmospheric boundary layer (ABL) precursor simulation was constructed based on assessments of meteorological tower, lidar, and radar data over a 3-hour window. This precursor was tuned to the specific atmospheric conditions that occurred both prior to and during the measurement campaign, enabling capture of a night-to-day transition in the turbulent ABL. In the absence of height-varying temperature measurements, spatially averaged radar data were sufficient to characterize the atmospheric stability of the wind plant in terms of the shear profile, and near-ground temperature sensors provided a reasonable estimate of the ground heating rate describing the morning transition. A full aeroelastic simulation was then performed for a subset of turbines within the wind plant, driven by the precursor. Analysis of two turbines within the array, one directly waked by the other, demonstrated good agreement with measured time-averaged loads.

17 WIND ENERGY↗
DOE OSTI2019

Comparison of Wind Farm Layout Optimization Results Using a Simple Wake Model and Gradient-Based Optimization to Large-Eddy Simulations: Preprint

Comparison of wind farm layout optimization results with simulations or wind tunnel tests is important because the models used during optimization use simplifying assumptions that can alter the design space. Some of the characteristics of the simple models may negatively influence the resulting layouts. In this paper we perform wind farm layout optimization and compare the resulting improvements to large eddy simulation (LES) results. We begin by describing the models used, including changes specific for use with gradient-based optimization. We then compare our version of the models with the original versions and large eddy simulation LES) data from previous publications. Using the models described, we performed gradient-based wind farm layout optimization, using exact gradients, and compare the results to LES simulations. The models predicted an improvement in annual energy production (AEP) of 7.4%, while the LES reported an AEP improvement of 9:9%. We conclude that the improvements found by optimizing with the simple models are not just an artifact of the models, but are real improvements.

17 WIND ENERGY↗
DOE OSTI2014

Evaluation of convection‐permitting model simulations of cloud populations associated with the Madden‐Julian Oscillation using data collected during the AMIE/DYNAMO field campaign

Abstract Regional convection‐permitting model simulations of cloud populations observed during the 2011 Atmospheric Radiation Measurement (ARM) Madden‐Julian Oscillation Investigation Experiment/Dynamics of the Madden‐Julian Oscillation Experiment (AMIE/DYNAMO) field campaign are evaluated against ground‐based radar and ship‐based observations. Sensitivity of model simulated reflectivity, surface rain rate, and cold pool statistics to variations of raindrop breakup/self‐collection parameters in four state‐of‐the‐art two‐moment bulk microphysics schemes in the Weather Research and Forecasting (WRF) model is examined. The model simulations generally overestimate reflectivity from large and deep convective cells, and underestimate stratiform rain and the frequency of cold pools. In the sensitivity experiments, introduction of more aggressive raindrop breakup or decreasing the self‐collection efficiency increases the cold pool occurrence frequency in all of the simulations, and slightly reduces the reflectivity and precipitation statistics bias in some schemes but has little effect on the overall mean surface precipitation. Both the radar observations and model simulations of cloud populations show an approximate power law relationship between convective echo‐top height and equivalent convective cell radius.

Hagos, Samson↗
DOE OSTI2017

PFLOTRAN-E4D: A parallel open source PFLOTRAN module for simulating time-lapse electrical resistivity data

Time-lapse electrical resistivity tomography (ERT) is finding increased application for remotely monitoring processes occurring in the near subsurface in three-dimensions (i.e. 4D monitoring). However, there are few codes capable of simulating the evolution of subsurface resistivity and corresponding tomographic measurements arising from a particular process, particularly in parallel and with an open source license. Herein we describe and demonstrate an electrical resistivity tomography module for the PFLOTRAN subsurface simulation code, named PFLOTRAN-E4D. The PFLOTRAN-E4D module operates in parallel using a dedicated set of compute cores in a master-slave configuration. At each time step, the master processes receives subsurface states from PFLOTRAN, converts those states to bulk electrical conductivity, and instructs the slave processes to simulate a tomographic data set. The resulting multi-physics simulation capability enables accurate feasibility studies for ERT imaging, the identification of the ERT signatures that are unique to a given process, and facilitates the joint inversion of ERT data with hydrogeological data for subsurface characterization. PFLOTRAN-E4D is demonstrated herein using a field study of stage-driven groundwater/river water interaction ERT monitoring along the Columbia River, Washington, USA. Results demonstrate the complex nature of changes subsurface electrical conductivity, in both the saturated and unsaturated zones, arising from water table changes and from river water intrusion into the aquifer. The results also demonstrate the sensitivity of surface based ERT measurements to those changes over time. PFLOTRAN-E4D is available with the PFLOTRAN development version with an open-source license at https://bitbucket.org/pflotran/pflotran-dev .

Environmental Molecular Sciences Laboratory↗