Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulator study”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Comprehensive Simulation Study to Evaluate Future Vehicle Energy and Cost Reduction Potential

Under the umbrella of EERE’s Office of Sustainable Transportation, the U.S. Department of Energy’s (DOE) Vehicle Technologies Office (VTO) and Hydrogen and Fuel Cell Technologies Office (HFTO) seek to develop sustainable, affordable, and efficient technologies for transportation of goods and people. Translating investments in advanced transportation component technologies and powertrains to estimate the potential for vehicle-level fuel savings is critical to understanding DOE’s impact and success in this mission For this study, Argonne National Laboratory (Argonne) simulated technologies funded by VTO and HFTO for light-duty vehicles across the following: Powertrain configurations (conventional, power-split hybrid electric vehicle, extended-range electric vehicle, battery electric drive, and fuel-cell vehicles); Vehicle classes (compact car, mid-size car, small sport utility vehicle [SUV], mid-size SUV, and pickup truck); Fuels (gasoline, diesel, natural gas, hydrogen, and battery electricity). We assessed each technology for five different timeframes: laboratory years 2015 (reference), 2020, 2025, 2030, and 2045. We assumed a delay of 5 years between laboratory year and model year (i.e., the year the technology is introduced into production). Finally, we included uncertainties for both technology performance and cost by considering two cases (note that these cases are not intended as predictions of future performance): Low case , aligned with DOE technology manager estimates of expected original equipment manufacturer (OEM) improvements based on business as usual regulatory and market environments; High case , aligned with aggressive technology advancements based on research and development (R&D) targets developed through support by VTO and HFTO. We estimated the energy and cost impact of different technologies using Autonomie (Argonne undated), a state-of-the-art vehicle system simulation tool developed by Argonne and used to assess the energy consumption, performance, and cost of multiple advanced vehicle technologies. The tool comprises a complete set of vehicle models to assess impacts across a wide range of classes (from light- to heavy-duty), powertrain configurations (from conventional to hybrid electric vehicles [HEVs], fuel cell electric vehicles [FCEVs], plug-in hybrid electric vehicles [PHEVs], and battery electric vehicles [BEVs]), components, and control strategies, including vehicle-level and component-level controls developed and calibrated using dynamometer test data. Autonomie has been used to support a wide range of studies: analyzing various component technologies, sizing powertrain components to meet different vehicle requirements, comparing the benefits of powertrain configurations, optimizing both heuristic and route-based vehicle energy control, and predicting transportation energy use when paired with a traffic modeling tool such as POLARIS. This report documents the assumptions made and the vehicle-level energy consumption benefits and associated technology costs estimated for various types of light-duty vehicles. Details regarding vehicle assumptions and simulation results are available in the spreadsheets accompanying this report.

08 HYDROGEN↗
DOE OSTI2023

Crosslinking and Gelation of Polymer Brushes and Free Polymer Chains in a Confined Space during Controlled Radical Polymerization-A Computer Simulation Study

The confinement effect on crosslinking and gelation during controlled/living copolymerization of a monovinyl and a divinyl monomer was studied by Monte Carlo computer simulation using the dynamic lattice liquid method. The simulated reactions took place in plane-parallel slits of various widths for various grafting densities (GDs) and various initial initiator/crosslinker ratios ([X] 0 /[Ini] 0 ). Monomers and crosslinkers were initially randomly distributed in the slit, but initiators were either anchored to one solid plane (polymer brushes) or randomly distributed (free chains). For brushes, the gel point (GP) depended on GD and [X] 0 /[Ini] 0 . The higher these parameters, the lower was the gel point. The gel point was also affected by the slit width. For free, nongrafted chains, the dependence of GP on slit width was negligible. For brushes, a nonuniform density of crosslinks was found, the highest density was close to the wall. This asymmetry was high at the gel point and decreased at higher conversions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mechanism of monovalent and divalent ion mobility in Nafion membrane: An atomistic simulation study

Polymer cation exchange membranes (CEMs) are widely used in water treatment processes. The fundamental factors that control the separation and selectivity of cations with different valences in CEMs are not fully understood. In this study, we use atomistic molecular dynamics simulations to investigate the underlying molecular mechanisms that control the mobility of cations with different valences in Nafion membranes. Our results indicate structural differences in binding of monovalent and divalent cations, which gives rise to differences in ion diffusion in Nafion. Monovalent cations are shown to be “territorially” bound, i.e., bound while partially hydrated, to the fixed charge groups whereas divalent cations are shown to be “site” bound, i.e., bound while fully dehydrated, to the charge groups on the polymer. This difference in binding structure gives rise to differences in transport characteristics of cations in Nafion.

Chemistry↗
DOE OSTI2024

Large-Eddy Simulation Study of Flow and Combustion Dynamics in a Full-Scale Hydrogen–Air Rotating Detonation Combustor-Stator Integrated System

In the present work, a first-of-its-kind three-dimensional (3D) large-eddy simulation (LES) study is conducted to numerically investigate the combustion dynamics as well as aero-thermal phenomena in a full-scale nonpremixed hydrogen–air rotating detonation engine (RDE) (with a diverging-shaped lower-end wall), when integrated with nozzle guide vanes (NGV) acting as the turbine stator. The wall-modeled LES framework incorporates hydrogen–air detailed chemical kinetics and adaptive mesh refinement (AMR). A comparative analysis is carried out for two operating conditions with different fuel/air mass flow rates but global equivalence ratio of unity, and considering RDE configurations without and with stator. The LES model is validated against available experimental data for the low mass flux condition with respect to detonation wave speed/height, wave dynamics, and axial static pressure distribution. Numerical results indicate significant deflagrative combustion occurring in the fill region near the inner wall due to formation of recirculation zones in the injection near-field driven by the backward facing step. The leading detonation wave is found to be trailed by an azimuthal reflected-shock combustion (ARSC) wave, consistent with experimental observations, which consumes unburned vitiated reactants that leak through the main detonation wave. The main detonation wave characteristics, such as detonation wave speed/height and combustion efficiency, do not change appreciably with the presence of NGV. A novel combustion diagnostic technique based on chemical explosive mode analysis (CEMA) is employed to quantify the fraction of heat release occurring in the detonative mode versus deflagrative mode for the simulated conditions. The exit flow is found to be nearly fully subsonic and supersonic for the low and high mass flux conditions, respectively. Further analysis of the exit flow profiles shows that the presence of NGV renders the flow more axial and significantly impacts the exit Mach number and total pressure, while the total temperature shows negligible change. In addition, the low mass flux operating point, despite exhibiting more deflagrative losses within the combustor, yields overall lower pressure drop from plenum to exhaust, which is mainly attributed to lower pressure drop across the injectors. Lastly, the rotating detonation engine-nozzle guide vanes (RDE-NGV) configuration exhibits higher total pressure loss compared to rotating detonation engine (RDE) without stator across both the mass flux conditions. In conclusion, this study extends the state-of-the-art in numerical modeling of pressure gain combustion (PGC) systems by demonstrating high-fidelity 3D reacting LES of full-scale RDE-NGV systems relevant to RDE-turbine integration for stationary power generation.

Combustion↗
DOE OSTI2023

Simulation Study of Monoenergetic Photon Source LiDAR System Performance for Emergency Response Applications

The following report summarizes the results of a Monte Carlo simulation and modeling study of the Lawrence Berkeley National Laboratory’s (LBL) Monoenergetic Photon Source (MPS) as part of a high-energy LiDAR system for use in emergency response applications. A set of image quality indicators (IQI) and nonsensitive phantoms or generalized surrogates for search targets have been identified for use in the planned measurement campaign at the Berkeley Lab Laser Accelerator (BELLA) Center toward the end FY2023. These items were modeled and the time-of-flight (TOF) spectral responses following irradiation were calculated for a set of varying operating conditions. The goal of this study was to both predict the LiDAR system performance as well as determine if minor modifications to the IQI and phantom designs were required before fabrication. It was found that the IQI and target designs both overlap and extend towards the limit of performance of the system and are thus suitable for the final project measurement campaign.

43 - PARTICLE ACCELERATORS↗
DOE OSTI2021

Kinetic simulation study of magnetized collisionless shock formation on a terawatt laser system

Perpendicular, magnetized, collisionless shocks in hydrogen and neon plasmas are studied with 2-D particle-in-cell simulations for parameters accessible to experiments on OMEGA EP. The simulations are performed with realistic ion-electron mass ratios by which the relative importance of different micro-instabilities can be accurately captured. The dispersion relation of the modified two-stream instability (MTSI), the main dissipation mechanism for these shocks, is used to find suitable parameters for upcoming experiments. Simulations show that magnetized collisionless shocks can be readily formed within a few tenths of an ion gyro-period in both hydrogen and neon gases, with a background magnetic field of 50 T, achievable using the magneto-inertial fusion electrical discharge system (MIFEDS). Here, a portion of the ions are reflected to the upstream region and accelerated in both shock normal and tangential directions, indicating the formation of a supercritical shock. Shock front reformation is seen in longer time 1-D simulations. The results show that the formation time and width of these shocks are determined by MTSI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

3-D Geological Modeling for Numerical Flow Simulation Studies of Gas Hydrate Reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. Here, this study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-stepping grid were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a three-dimensional (3-D) trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

58 GEOSCIENCES↗
DOE OSTI2023

Graphic contrastive learning analyses of discontinuous molecular dynamics simulations: Study of protein folding upon adsorption

A comprehensive understanding of the interfacial behaviors of biomolecules holds great significance in the development of biomaterials and biosensing technologies. In this work, we used discontinuous molecular dynamics (DMD) simulations and graphic contrastive learning analysis to study the adsorption of ubiquitin protein on a graphene surface. Our high-throughput DMD simulations can explore the whole protein adsorption process including the protein structural evolution with sufficient accuracy. Contrastive learning was employed to train a protein contact map feature extractor aiming at generating contact map feature vectors. Subsequently, these features were grouped using the k-means clustering algorithm to identify the protein structural transition stages throughout the adsorption process. The machine learning analysis can illustrate the dynamics of protein structural changes, including the pathway and the rate-limiting step. Our study indicated that the protein–graphene surface hydrophobic interactions and the π–π stacking were crucial to the seven-stage adsorption process. Upon adsorption, the secondary structure and tertiary structure of ubiquitin disintegrated. The unfolding stages obtained by contrastive learning-based algorithm were not only consistent with the detailed analyses of protein structures but also provided more hidden information about the transition states and pathway of protein adsorption process and structural dynamics. Our combination of efficient DMD simulations and machine learning analysis could be a valuable approach to studying the interfacial behaviors of biomolecules.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Two-stage membrane-based process utilizing highly CO 2 -selective membranes for cost and energy efficient carbon capture from coal flue gas: A process simulation study

Membrane technology for CO 2 capture has become an attractive strategy due to its cost and energy efficiency and low materials costs. In the past decade, membrane-based process designs for post-combustion power plant CO 2 capture have been developed, utilizing existing highly CO 2 -permeable membranes with relatively low CO 2 /N 2 selectivity (<50), and have obtained reasonably economic carbon capture. However, few membrane-based process designs were proposed for moderate to highly CO 2 -selective membranes (CO 2 /N 2 selectivity of 50–300, and >300, respectively), which have vastly emerged in recent years, such as various facilitated transport membranes (FTMs). Herein, we proposed a two-stage membrane-base process design targeting economic carbon capture from coal-fired flue gas. This process design features the utilization of highly CO 2 -selective membranes for one-stage CO 2 enrichment to 95% dry-base purity in the first stage and recycle of the remaining CO 2 by a highly CO 2 -permeable membrane in the second stage, in order to achieve economic CO 2 capture with 90% capture rate and >95% CO 2 product purity. Through an integration-iteration membrane model and the Aspen Plus process simulation, a sensitivity study of operating pressures (feed and permeate pressures) and membrane properties (CO 2 permeance and CO 2 /N 2 selectivity) was conducted. Critical CO 2 /N 2 selectivity of 300–400 was found for the highly CO 2 -selctive membranes to meet the demand for cost and energy efficient results. The lowest possible membrane area of 4.8 × 10 5 m 2 and fractional energy of 19.3% were obtained, which is comparable to or even more attractive than reported membrane-based process designs. Here, this work provides a new membrane process design option for highly CO 2 -selective membranes and gives insights on the influence of membrane performance and operation condition.

42 ENGINEERING↗
DOE OSTI2022

Simulation study of chemo-mechanical impacts of CO 2 injection in morrow b sandstone reservoir

Here, this study investigates the impacts of geomechanical and geochemical changes on carbon storage in a partially depleted oil reservoir, using results from four different coupled simulation models. Models were used to examine the relative importance of storage mechanisms, and how changing reservoir parameters might affect these mechanisms through time. The study uses data from a Morrowan sandstone reservoir in the Farnsworth Unit (FWU), Ochiltree County, Texas which is currently undergoing CO 2 enhanced oil recovery (EOR). Partially depleted oil reservoirs such as the FWU offer attractive carbon utilization and/or storage targets because of existing infrastructure and economic benefits from incremental oil recovery as well as tax credits. However, prediction of storage capacity or long-term fluid migration in these fields can be difficult because of the wide variation in formation fluids and operational histories that may have undergone. CO 2 injection can cause complex geomechanical and geochemical responses in a reservoir as a result of interplay between dynamic changes in pore pressure, reservoir temperature, fluid composition, and interactions between formation fluids, CO 2 , and reservoir rock. Thus, multiple coupled numerical simulation models must be developed and used to more precisely understand what CO 2 storage mechanisms are most significant, as well as the long-term fate of the stored CO 2 . Our study used results from hydrodynamic, coupled hydro-geomechanical, coupled hydro-geochemical, and coupled hydro-geomechanical-geochemical models to examine how changes in geomechanical and geochemical properties can impact the injectivity or storage capacity of CO 2 . Models simulated historical field operations and then forward-modeled a water-alternate-gas (WAG) operation for 20 years, followed by a 1000-year post-injection monitoring. The work demonstrates that in this specific reservoir, geomechanical impacts appear to be more significant than any geochemical processes, which showed minimal impact on reservoir properties. Furthermore, these simulation scenarios corroborated other studies showing more than 90% of CO 2 injected are being stored and will remain stored in the reservoir. The study was able to draw on a rich set of data from this currently operational project for the work. Because the Morrow is a widespread oil-bearing reservoir in the midcontinent of the United States and is fairly typical of many sandstone reservoirs elsewhere, the study has significance beyond the boundaries of our field area.

58 GEOSCIENCES↗
DOE OSTI2023

Microstructure-process relationships in monolithic U-10Mo fuel foil single-pass rolling: A parametric simulation study

In this work, a previously validated coupling of Kinetic Monte Carlo (KMC) Potts Model and finite element method (FEM) simulations was implemented to investigate the effects of microstructural features in as-cast and homogenized monolithic U-10Mo foils on the emergent microstructure after rolling and reheating. Parameters that could potentially affect recrystallization behavior of the rolled U-10Mo foils were considered: grain size distribution, uranium carbide (UC) size distribution, UC volume fraction, spatial distribution of UC, and rolling reduction magnitude. Grain structure and the magnitude of rolling reduction have the strongest influence on recrystallization kinetics and the fabricated grain size distribution. The UC distribution had only a weak effect on the recrystallization kinetics and final microstructures. While particle-stimulated nucleation (PSN) occurred in simulation more frequently as grain size increased, its incidence did not appear to considerably affect the recrystallization kinetics or grain size distribution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Organomercury Measurements from Sludge Batch 10 Simulant Studies

Researchers at the Savannah River National Laboratory have recently performed testing to evaluate the Sludge Batch 10 flowsheet using simulated sludge waste. In the course of this testing several samples were taken to determine the concentration of organomercury species resulting from sludge batch processing. These samples were submitted to the Savannah River National Laboratory Sensing and Metrology department quantitation using a mercury analyzer. The signal amplifier used to perform organomercury quantitation in the mercury analyzer experienced a malfunction in the course of analysis, leading to uncertainty of organomercury concentrations observed. Seventy-three organomercury samples from six project submissions to the Savannah River National Laboratory Sensing and Metrology department have been critically reviewed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Simulation study of fast ion losses associated with the rotating n = 1 resonant magnetic perturbations in KSTAR

In this work we describe a simulation framework for testing the fast ion loss mechanism associated with the experimentally observed beam ion losses when an externally applied toroidally rotating perturbed magnetic field is used to control edge localized modes in the KSTAR tokamak. The simulations reproduce the key qualitative features of neutral beam injection (NBI) ion detection by a fast ion loss detector. The NBI ion losses in the simulation mainly occur for passing particles due to orbit stochastization, which is caused by orbit resonance with the 3D field perturbations. The relative toroidal angle of NBI ion deposition to the orbit island determines the radial path of the particles between confinement or loss. The fast ion loss quantity depends on the toroidal angle rotation of the 3D field with respect to the beam deposition position. The averaged transport of toroidal canonical angular momentum reveals that linear and nonlinear resonance of the NBI ions within the 3D field are the main factors determining fast ion transport and loss.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Impact of refrigerant undercharge faults on building indoor conditions and HVAC system operation in residential Buildings: A simulation study

This study investigates the impact of refrigerant undercharge on indoor temperature and HVAC system performance in residential buildings. Simulation models for typical residential buildings in Orlando, FL and Indianapolis, IN were developed using the ResStock database. A refrigerant undercharge fault model was then applied to the simulations with varying levels of fault intensity. The paper offers an extensive analysis, revealing that variations in supply air temperature, equipment runtime, and cooling energy consumption due to the level of refrigerant undercharge faults are notably significant on a summer representative day. Similarly, on a winter representative day, changes in supply air temperature and runtime are significant as well as changes in supplemental heat energy consumption. We find that occupants may remain oblivious to these faults during the cooling season, particularly when the HVAC system is oversized; in that case, supply air temperature data could help detect a fault. Another challenge is that during the heating season, when the supplemental heater operates, it is difficult to identify a refrigerant undercharge fault using only indoor and supply air temperature data. Finally, this study finds that supply air temperature, equipment runtime, and supplemental heater energy consumption data can help in detecting refrigerant undercharge faults.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

A direct numerical simulation study for confined non-isothermal jet impingement at moderate nozzle-to-plate distances: Capturing jet-to-ambient density effects

A direct numerical simulation (DNS) campaign is deployed for a series of confined downward oriented, non-isothermal turbulent impinging jet configurations. A baseline Reynolds number of 9960 is obtained through a precursor DNS pipe flow simulation (Re τ = 505). Three jet temperature configurations (confinement height to nozzle diameter of three) enter a cylindrical domain that share ambient and impingement plate temperatures (298.15 K). The range of jet temperatures are crafted such that the ratio of inlet to ambient density varies from unity to 0.52, showcasing the effect of density disparity on flow characteristics such as core collapse, radial mixing of momentum and energy, near-wall stagnation behavior, wall-jet profiles, and large-scale vortical structures. Surface quantities provided include mean radial heat flux and wall-shear stress profiles, and heat flux histograms at select radial stations. Results showcase increased radial normal stresses for higher temperature jets that support increased mixing, resulting in large-scale recirculation structures that are smaller, while retaining similar normalized radial wall profiles for shear stress, heat flux and pressure. Furthermore, radial plots for wall shear stress and Nusselt number showcase strong radial decay as compared to previous configurations that share similar jet and ambient temperatures. For the 373.15 K case, a Gaussian-like histogram for heat fluxes at the impingement plate transitions to a log-normal profile as radial distances increase. In contrast, the 573.15 K configuration displays a bi-modal heat flux characteristic at the impingement plate, and in similar manner to the moderate temperature counterpart, transitions to a log-normal profile at larger radial distances.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Wind farm response to mesoscale-driven coastal low level jets: a multiscale large eddy simulation study

Realistic atmospheric turbulence–wind farm interactions during coastal low-level jet (LLJ) events are captured using high-fidelity, mesoscale-driven large eddy simulations (LES) to understand wind turbine loads, wakes and overall performance. The simulation has been carried out using the ExaWind aeroelastic solver, AMR-Wind. The simulations have been compared against a baseline unstable case matching the wind speed, wind direction and TI at hub-height location. Results indicate that the LLJ has negative impacts on the turbine hub and tower loads, and opens up potential avenues for design load mitigation strategies.

17 WIND ENERGY↗