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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Nonlinear fluid simulation of particle and heat fluxes during burst of ELMs on DIII-D with BOUT++ code [Fluid Simulation of Particle and Heat Fluxes during Burst of ELMs on DIID with BOUT++ code]

In order to study the distribution and evolution of the transient particle and heat fluxes during edge-localized mode (ELM) bursts, a BOUT++ six-field two-fluid model based on the Braginskii equations with non-ideal physics effects is used to simulate pedestal collapse in divertor geometry. We used the profiles from the DIII-D H-mode discharge #144382 with fast target heat flux measurements as the initial conditions for the simulations. Moreover, a flux-limited parallel thermal conduction is used with three values of the flux-limiting coefficient ${{\alpha}_{j}}$ , free streaming model with ${{\alpha}_{j}}=1$ , sheath-limit with ${{\alpha}_{j}}=0.05$ , and one value in between. The studies show that a 20 times increase in ${{\alpha}_{j}}$ leads to ~6 times increase in the heat flux amplitude to both the inner and outer targets, and the widths of the fluxes are also expanded. The sheath-limit model of flux-limiting coefficient is found to be the most appropriate one, which shows ELM sizes close to the measurements. The evolution of the density profile during the burst of ELMs of DIII-D discharge #144382 is simulated, and the collapse in width and depth of ${{n}_{\text{e}}}$ are reproduced at different time steps. The growing process of the profiles for the heat flux at divertor targets during the burst of ELMs measured by IRTV (infrared television) is also reproduced by this model. The widths of heat fluxes towards targets are a little narrower, and the peak amplitudes are twice the measurements possibly due to the lack of a model of divertor radiation which can effectively reduce the heat fluxes. The magnetic flutter combined with parallel thermal conduction is found to be able to increase the total heat loss by around 33% since the magnetic flutter terms provide the additional conductive heat transport in the radial direction. Finally, the heat flux profile at both the inner and outer targets is obviously broadened by magnetic flutter. The lobe structures near the X-point at LFS are both broadened and elongated due to the magnetic flutter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

SARS-CoV2 Protein-Ligand Simulation Dataset: Layer 1 (Simulation Initial Conditions and Parameters)

A set of 24 protein structures/complexes from the SARS CoV-2 proteome and inputs prepared for simulation using the CHARMM36m forcefield in PDB and gromacs formats. Each system contains a pdb and gromacs top and related input files necessary for running a temperature replica-exchange simulation. In addition, we also include: charmm PSF files (generated from the gromacs topology), a list of temperatures at which replica-exchange simulations were done (tempRamp), example gromacs run input mdp files, and initial minimized structures where available (minimized.pdb).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Multiscale thermal hydraulic coupling methods for boiling water reactor simulation

During the last two years, the Virtual Environment for Reactor Applications (VERA) has been extended to simulate boiling water reactors (BWRs). The thermal hydraulic effects present in BWRs are far more complex than those in pressurized water reactors (PWRs). Therefore, the runtime is significantly increase and convergence behavior worse compared to PWRs. Most of the additional runtime is spent during the first few coupled iterations, when the power shape is still rapidly evolving, dramatically affecting the thermal hydraulics (TH). To alleviate the increased computational expense, a multiscale TH coupling approach was developed in VERA in which a highly efficient, simplified TH model is solved for several coupled iterations until the power shape is partially converged. The simplified solution is used to inform the assembly-wise flow distribution in the high-fidelity TH solver CTF, reducing the amount of work required to properly balance the pressure drop in each channel and reducing runtime for the couple calculation. Because the final TH calculations are performed with CTF, there is ultimately no impact on the accuracy of the converged solution. This paper presents the details of this multiscale TH coupling approach, along with results for a 4x4 array of GE-14 fuel bundles and whole-core coupled simulations of the Hatch core. These two cases show that the multiscale approach dramatically reduces the runtime of coupled BWR simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

An HPC-Based Hydrothermal Finite Element Simulator for Modeling Underground Geothermal Behavior with Example Simulations on The Treasure Island and UC Berkeley Campus

This submission contains the source code of the Hydrothermal Finite Element Simulator used for the Treasure Island and UC Berkeley campus geothermal simulation. It contains a report that summarizes the development and validation of this Hydrothermal Finite Element Simulator, with a case study on Treasure Island site. It also contains a report that investigates the feasibility of upgrading the existing campus energy delivery system at UC Berkeley to a fifth-generation district heating and cooling system that includes geothermal heat/cold storage.

15 GEOTHERMAL ENERGY↗
DOE OSTI2019

Modeling and Simulation - Introducing hardware-in-the-loop capabilities to the Human Systems Simulation Laboratory

The DOE Office of Nuclear Energy (DOE-NE) cybersecurity program would like to expand the use of the Human Systems Simulation Laboratory to include capabilities for cybersecurity research and development. Specifically, the introduction of hardware-in-the-loop (HITL) simulation of nuclear power plant architecture into the HSSL will enable researchers and nuclear utilities the ability to identify impacts generated by sophisticated hostile cyber adversaries to identify the consequences of simulated potential attack vectors, and to develop solutions for detecting or preventing these cyber attacks.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Production Through Simulation: Using Simulation Technologies to Create and Evaluate Nuclear Fuel Fabrication Facility Designs

This project explores how physics-based and discrete-event simulation technologies can be jointly employed to model a suite of nuclear fuel fabrication facility designs with variable radiological environments, optimize their design to lower operational costs subject to a set of design constraints, and evaluate their performance. Initially, three variants of U-20Pu-10Zr metallic nuclear fuel in three specific geometric configurations are modeled and studied across more than 50 physics simulations in nuclear physics software packages SCALE and Monte Carlo N-Particle Transport (MCNP6.2). From these tests, values for effective dose rate and criticality are recorded for each of the nine alloy-geometry configurations. These values are then incorporated into discrete-event models of fuel fabrication in the simulation program ExtendSim Pro 10 as process attributes. This allows one to track the dose to facility personnel as fuel is fabricated, dually enabling us to evaluate the safety of a given facility design in terms of absorbed dose per worker per year and to craft operational guidelines that adhere to federal and local safety regulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Large-eddy simulation of an atmospheric bore and associated gravity wave effects on wind farm performance in the southern Great Plains

Gravity waves are a common occurrence in the atmosphere, with a variety of generation mechanisms. Their impact on wind farms has only recently gained attention, with most studies focused on wind farm-induced gravity waves. In this study, the interaction between a wind farm and gravity waves generated by an atmospheric bore event is assessed using multiscale large-eddy simulations. The atmospheric bore is created by a thunderstorm downdraft from a nocturnal mesoscale convective system (MCS). The associated gravity waves impact the wind resource and power production at a nearby wind farm during the American Wake Experiment (AWAKEN) in the US southern Great Plains. A two-domain nested setup (Δx=300 and 20 m) is used in the Weather Research and Forecasting (WRF) model, forced with data from the High-Resolution Rapid Refresh model, to capture both the formation of the bore and its interaction with individual wind turbines. The MCS is resolved on the large outer domain, where the structure of the bore and the associated gravity waves are found to be especially sensitive to parameterized microphysics processes. On the finer inner domain, gravity wave interactions with individual wind turbines are resolved; wake dynamics are captured using a generalized actuator disk parameterization in WRF. The gravity waves are found to have a strong effect on the atmosphere above the wind farm; however, the effect of the waves is more nuanced closer to the surface where there is additional turbulence, both ambient and wake-generated. Notably, the gravity waves modulate the mesoscale environment by weakening and dissipating the preexisting low-level jet, which reduces hub-height wind speed and hence the simulated power output, which is confirmed by the observed supervisory control and data acquisition (SCADA) power data. Additionally, the gravity waves induce local wind direction variations correlated with fluctuations in pressure, which lead to fluctuations in the simulated power output as various turbines within the farm are subjected to waking from nearby turbines.

17 WIND ENERGY↗
NASA NTRS1994

Direct simulation Monte Carlo and Navier-Stokes simulations of blunt body wake flows

Numerical results obtained with direct simulation Monte Carlo and Navier-Stokes methods are presented for a Mach-20 nitrogen flow about a 70-deg blunted cone. The flow conditions simulated are those that can be obtained in existing low-density hypersonic wind tunnels. Three sets of flow conditions are considered with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is on the wake structure: how the wake structure changes as a function of rare faction, what the afterbody levels of heating are, and to what limits the continuum models are realistic as rarefunction in the wake is progressively increased. Calculations are made with and without an afterbody sting. Results for the afterbody sting are emphasized in anticipation of an experimental study for the current flow conditions and model configuration. The Navier-Stokes calculations were made with and without slip boundary conditions. Comparisons of the results obtained with the two simulation methodologies are made for both flowfield structure and surface quantities.

WAKE STRUCTURE↗
DOE OSTI2017

Investigating the dependence of SCM simulated precipitation and clouds on the spatial scale of large-scale forcing at SGP [Investigating the scale dependence of SCM simulated precipitation and cloud by using gridded forcing data at SGP]

Large-scale forcing data, such as vertical velocity and advective tendencies, are required to drive single-column models (SCMs), cloud-resolving models, and large-eddy simulations. Previous studies suggest that some errors of these model simulations could be attributed to the lack of spatial variability in the specified domain-mean large-scale forcing. This study investigates the spatial variability of the forcing and explores its impact on SCM simulated precipitation and clouds. A gridded large-scale forcing data during the March 2000 Cloud Intensive Operational Period at the Atmospheric Radiation Measurement program's Southern Great Plains site is used for analysis and to drive the single-column version of the Community Atmospheric Model Version 5 (SCAM5). When the gridded forcing data show large spatial variability, such as during a frontal passage, SCAM5 with the domain-mean forcing is not able to capture the convective systems that are partly located in the domain or that only occupy part of the domain. This problem has been largely reduced by using the gridded forcing data, which allows running SCAM5 in each subcolumn and then averaging the results within the domain. This is because the subcolumns have a better chance to capture the timing of the frontal propagation and the small-scale systems. As a result, other potential uses of the gridded forcing data, such as understanding and testing scale-aware parameterizations, are also discussed.

58 GEOSCIENCES↗
DOE OSTI2015

Simulating and evaluating best management practices for integrated landscape management scenarios in biofuel feedstock production: Evaluating Best Management Practices for Biofuel Feedstock Production

Sound crop and land management strategies can maintain land productivity and improve the environmental sustainability of agricultural crop and feedstock production. This study evaluates the improvement of water sustainability through an integrated landscaping management strategy, where landscaping design, land management operations, crop systems, and agricultural best management practices (BMPs) play equal roles. The strategy was applied to the watershed of the South Fork Iowa River in Iowa, with a focus on implementing riparian buffers and converting low productivity land to provide cellulosic biomass while benefiting soil and water quality. The Soil and Water Assessment Tool (SWAT) was employed to simulate the impact of integrated landscape design on nutrients, suspended sediments, and flow on the watershed and subbasin scales. First, the study evaluated the representation of buffer strip as a vegetative barrier and as a riparian buffer using trapping efficiency and area ratio methods in SWAT. For the riparian buffer, the area ratio method tends to be more conservative, especially in nitrate loadings, while the trapping efficiency method generates more optimistic results. The differences between the two methods increase with buffer width. The two methods may not be comparable for the field-scale vegetative barrier simulation because of limitations in model spatial resolution. Landscape scenarios were developed to quantify water quality under (1) current land use, (2) partial land conversion to switchgrass, and (3) riparian buffer implementation. Results show that when low productivity land (15.2% of total watershed land area) is converted to grow switchgrass, suspended sediment, total nitrogen, total phosphorus, and nitrate loadings are reduced by 69.3%, 55.5%, 46.1%, and 13.4%, respectively, in the watershed surface streams. The reduction was less extensive when riparian buffer strips (30 m or 50 m) were applied to the stream network at 1.4% of total land area in the watershed. At the subbasin level, the degree of nutrient and suspended sediment reduction varies extensively, ranging from a few percent up to 55%. Results indicate that effective landscape design on current agricultural land can potentially bring marked improvements in water quality and soil erosion control while producing food and fuel feedstock in the South Fork Iowa River watershed. The concept can be integrated with other watershed management programs to improve sustainability of land, water, and the ecosystem.

bioenergy production↗
DOE OSTI2021

Adaptive downstream tensorial eddy viscosity for hybrid large‐eddy simulations‐Reynolds‐averaged‐Navier–Stokes simulations

An adaptive method is developed to improve the accuracy of eddy‐viscosity Reynolds‐averaged‐Navier–Stokes (RANS) model in hybrid large‐eddy simulations (LES)‐RANS simulations by using available upstream LES results. The method first gets the tensorial eddy viscosity from the upstream LES solution at the LES‐RANS interface and then uses that information to improve the downstream RANS model by invoking the weak‐equilibrium assumption. The proposed method was evaluated via two test problems—flow in a channel and over a periodic hill. Results obtained show the proposed approach to increase the accuracy and stability of hybrid LES‐RANS simulations. Since the modification of the downstream RANS model is based on the tensorial eddy viscosity from the upstream LES solution, the method is adaptive to the problem being studied.

Zhang, Wanjia↗
DOE OSTI2014

Fast Multiresolution Reads of Massive Simulation Datasets

Today’s massively parallel simulation codes can produce output ranging up to many terabytes of data. Utilizing this data to support scientific inquiry requires analysis and visualization, yet the sheer size of the data makes it cumbersome or impossible to read without computational resources similar to the original simulation. We identify two broad classes of problems for reading data and present effective solutions for both. The first class of data reads depends on user requirements and available resources. Tasks such as visualization and user-guided analysis may be accomplished using only a subset of variables with a restricted spatial extent at a reduced resolution. The other class of reads requires full resolution multivariate data to be loaded, for example to restart a simulation. We show that utilizing the hierarchical multiresolution IDX data format enables scalable and efficient serial and parallel read access on a variety of hardware from supercomputers down to portable devices. Here, we demonstrate interactive view-dependent visualization and analysis of massive scientific datasets using low-power commodity hardware, and we compare read performance with other parallel file formats for both full and partial resolution data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

OpenRBC: A Fast Simulator of Red Blood Cells at Protein Resolution

We present OpenRBC, a coarse-grained molecular dynamics code, which is capable of performing an unprecedented in silico experiment—simulating an entire mammal red blood cell lipid bilayer and cytoskeleton as modeled by multiple millions of mesoscopic particles—using a single shared memory commodity workstation. To achieve this, we invented an adaptive spatial-searching algorithm to accelerate the computation of short-range pairwise interactions in an extremely sparse three-dimensional space. The algorithm is based on a Voronoi partitioning of the point cloud of coarse-grained particles, and is continuously updated over the course of the simulation. The algorithm enables the construction of the key spatial searching data structure in our code, i.e., a lattice-free cell list, with a time and space cost linearly proportional to the number of particles in the system. The position and the shape of the cells also adapt automatically to the local density and curvature. The code implements OpenMP parallelization and scales to hundreds of hardware threads. It outperforms a legacy simulator by almost an order of magnitude in time-to-solution and >40 times in problem size, thus providing, to our knowledge, a new platform for probing the biomechanics of red blood cells.

Biophysics↗