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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Hydrodynamic Interactions in Topologically Linked Ring Polymers.

Despite decades of interdisciplinary research on topologically linked ring polymers, their dynamics remain largely unstudied. These systems represent a major scientific challenge as they are often subject to both topological and hydrodynamic interactions (HI), which render dynamical solutions either mathematically intractable or computationally prohibitive. Here we circumvent these limitations by preaveraging the HI of linked rings. We show that the symmetry of ring polymers leads to a hydrodynamic decoupling of ring dynamics. This decoupling is valid even for nonideal polymers and nonequilibrium conditions. Physically, our findings suggest that the effects of topology and HI are nearly independent and do not act cooperatively to influence polymer dynamics. We use this result to develop highly efficient Brownian dynamics algorithms that offer enormous performance improvements over conventional methods and apply these algorithms to simulate catenated ring polymers at equilibrium, confirming the independence of topological effects and HI. The methods developed here can be used to study and simulate large systems of linked rings with HI, including kinetoplast DNA, Olympic gels, and poly[n]catenanes.

Rauscher, Phillip M.↗
DOE OSTI2023

Aviation Mission Profile Performance Evaluation in a High Specific Power Surface Permanent Magnet Machine with Additively Manufactured Coils and Integrated Heat Pipes

The purpose of this paper is to accurately evaluate the performance and overall efficiency of an electrical machine driving a propeller (as a part of a hybrid or electric propulsion system) during the flight mission profile using electromagnetic-thermal simulations. The study focuses on a 250 kW, 5000 rpm, 1.25 kHz, 36/30 slot/pole surface permanent magnet (SPM) machine with additively manufactured (AM) hollow conductors and integrated heat pipes (HP). The machine drive parameters (frequency, and current) are tuned to maximize the efficiency for a given flight cycle, while providing the highest specific power (SP) during takeoff. Afterward using Multiphysics simulations, the transient thermal behavior of the machine is studied and temperature profile for different parts of the machine is extracted. Finally, steady-state temperatures during cruising are estimated and overall efficiency is re-evaluated.

additive manufacturing↗
DOE OSTI2021

Reduced-order modeling of near-field THMC coupled processes for nuclear waste repositories in shale

Performance assessment (PA) of geologic radioactive waste repositories requires three-dimensional simulation of highly nonlinear, thermo-hydro-mechanical-chemical (THMC), multiphase flow and transport processes across many kilometers and over tens to hundreds of thousands of years. Integrating the effects of a near-field geomechanical process (i.e. buffer swelling) into coupled THC simulations through reduced-order modeling, rather than through fully coupled geomechanics, can reduce the dimensionality of the problem and improve computational efficiency. In this study, PFLOTRAN simulations model a single waste package in a shale host rock repository, where re-saturation of a bentonite buffer causes the buffer to swell and exert stress on a highly fractured disturbed rock zone (DRZ). Three types of stress-dependent permeability functions (exponential, modified cubic, and Two-part Hooke’s law models) are implemented to describe mechanical characteristics of the system. Overall, our modeling study suggests that compressing fractures reduces DRZ permeability, which could influence the rate of radionuclide transport and exchange with corrosive species in host rock groundwater that could accelerate waste package degradation. Less permeable shale host rock delays buffer swelling, consequently retarding DRZ permeability reduction as well as chemical transport within the barrier system.

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

DEVELOPMENT OF PLASMA SIMULATION TOOL FOR THE ACCEL-ERATING CAVITIES

Plasma processing of superconducting radio frequency (SRF) cavities has shown an improvement in accelerating gradient by reducing the radiation due to field emission and multipacting. Plasma processing is a common tech-nique where the free oxygen produced by the plasma breaks down and removes hydrocarbons from surfaces. This increases the work function and reduces the second-ary emission coefficient. The hydrocarbon fragments of H2, CO, CO2, and H2O are removed from the system with the process gas which is flowing through the system. Here, we present COMSOL for the first time to simulate the plasma processing of an SRF cavity. In this work, we use Jefferson Lab's C75 SRF cavities design as our case study. Using simulation, we predict the condition of plasma ignition inside the SRF cavity. The simulation provides information about the optimal rf coupling to the cavity, mode for plasma ignition, choice of gas concen-tration, power, and pressure.

Senevirathne, I.↗
DOE OSTI2026

Dependence of wind-farm-induced gravity waves and wind farm performance on non-dimensional atmospheric parameters and simulation configuration

This large-eddy simulation (LES) study examines how wind-farm-induced atmospheric gravity waves (AGWs) and wind farm performance depend on non-dimensional atmospheric parameters and simulation configuration. A hypothetical aligned wind farm of actuator disks is simulated under neutral surface conditions, with a stable capping inversion and a mildly stable free atmosphere, to assess the effects of stratification beyond the atmospheric boundary layer (ABL) on ABL flow. Simulation set-ups fully resolving AGWs are validated to minimize spurious wave generation and reflection from the domain boundaries. The validated set-up is then used to analyze AGW types and characteristics, as well as stratification impacts under conventionally neutral boundary layer (CNBL) conditions. These conditions are governed by four non-dimensional parameters: the Froude numbers of the free atmosphere and capping inversion (Fr, Fr i ), and the aspect ratios of the ABL and wind farm (H̃ i , S h ). Simulation configurations that fully resolve AGWs – capturing at least one wavelength both horizontally and vertically – yield the most realistic stratification effects on ABL flow, whereas partial or unresolved configurations produce non-physical, channel-like behavior. A coherent description of the AGW phenomena is provided, highlighting the central role of capping inversion displacement in linking ABL fluctuations with AGWs. Trapped waves are confined within the capping inversion, while interfacial and internal waves aloft are identified as the AGW types most relevant to wind farm performance. The wavy inversion, analogous to an interfacial wave, forms converging and diverging zones that drive power fluctuations across the farm. The interfacial wavelength, measured over the wind farm, corresponds to one diverging, one converging, and one mildly diverging zone. As the interfacial wavelength decreases with Fr i , multiple convergence–divergence zones develop under sub-critical conditions (Fr i <1.0), while for super-critical conditions (Fr i > 1.0), the wavelength approaches the farm length. Wave amplitude increases with decreasing H̃ i (i.e., shallower capping inversions). Wind farm performance is most sensitive to H̃ i : shallow boundary layers increase blockage and reduce efficiency, while deeper layers enhance efficiency. Increasing Fr and Fr i mitigates blockage, and increasing S h mainly improves wake recovery. Although local power fluctuations arise from AGWs, overall wind farm efficiency remains nearly constant with Fr and Fr i , improving primarily with larger H̃ i and S h .

17 WIND ENERGY↗
DOE OSTI2024

Nucleosynthetic Analysis of Three-dimensional Core-collapse Supernova Simulations

Abstract We study in detail the ejecta conditions and theoretical nucleosynthetic results for 18 three-dimensional core-collapse supernova (CCSN) simulations done by F ornax . Most of the simulations are carried out to at least 3 s after bounce, which allows us to follow their longer-term behaviors. We find that multidimensional effects introduce many complexities into the ejecta conditions. We see a stochastic electron fraction evolution, complex peak temperature distributions and histories, and long-tail distributions of the time spent within nucleosynthetic temperature ranges. These all lead to substantial variation in CCSN nucleosynthetic yields and differences from 1D results. We discuss the production of lighter α -nuclei, radioactive isotopes, heavier elements, and a few isotopes of special interest. Comparing pre-CCSN and CCSN contributions, we find that a significant fraction of elements between roughly Si and Ge are generally produced in CCSNe. We find that 44 Ti exhibits an extended production timescale as compared to 56 Ni, which may explain its different distribution and higher than previously predicted abundances in supernova remnants such as Cas A and SN1987A. We also discuss the morphology of the ejected elements. This study highlights the high-level diversity of ejecta conditions and nucleosynthetic results in 3D CCSN simulations and emphasizes the need for additional long-term (∼10 s) 3D simulations to properly address such complexities.

Astronomy & Astrophysics↗
DOE OSTI2021

A Case Study on Pathogen Transport, Deposition, Evaporation and Transmission: Linking High-Fidelity Computational Fluid Dynamics Simulations to Probability of Infection

A high-fidelity, low-Mach computational fluid dynamics simulation tool that includes evaporating droplets and variable-density turbulent flow coupling is well-suited to ascertain transmission probability and supports risk mitigation methods development for airborne infectious diseases such as COVID-19. A multi-physics large-eddy simulation-based paradigm is used to explore droplet and aerosol pathogen transport from a synthetic cough emanating from a kneeling humanoid. For an outdoor configuration that mimics the recent open-space social distance strategy of San Francisco, maximum primary droplet deposition distances are shown to approach 8.1 m in a moderate wind configuration with the aerosol plume transported in excess of 15 m. In quiescent conditions, the aerosol plume extends to approximately 4 m before the emanating pulsed jet becomes neutrally buoyant. A dose–response model, which is based on previous SARS coronavirus (SARS-CoV) data, is exercised on the high-fidelity aerosol transport database to establish relative risk at eighteen virtual receptor probe locations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Chemical Feedback in the Self-Assembly and Function of Air–Liquid Interfaces: Insight into the Bottlenecks of CO 2 Direct Air Capture

As fossil fuels remain a major source of energy throughout the world, developing efficient negative emission technologies, such as direct air capture (DAC), which remove carbon dioxide (CO 2 ) from the air, becomes critical for mitigating climate change. Although all DAC processes involve CO 2 transport from air into a sorbent/solvent, through an air–solid or air–liquid interface, the fundamental roles the interfaces play in DAC remain poorly understood. In this work, we study the interfacial behavior of amino acid (AA) solvents used in DAC through a combination of vibrational sum frequency generation spectroscopy and molecular dynamics simulations. This study revealed that the absorption of atmospheric CO 2 has antagonistic effects on subsequent capture events that are driven by changes in bulk pH and specific ion effects that feedback on surface organization and interactions. Among the three AAs (leucine, valine, and phenylalanine) studied, we identify and separate behaviors from CO 2 loading, chemical changes, variations in pH, and specific ion effects that tune structural and chemical degrees of freedom at the air–aqueous interface. The fundamental mechanistic findings described here are anticipated to enable new approaches to DAC based on exploiting interfaces as a tool to address climate change.

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

Structure and dynamics of $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquids close to Earth's liquid core conditions

Using an artificial neural-network machine learning interatomic potential, we have performed molecular dynamics simulations to study the structure and dynamics of $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquid close to the Earth's liquid core conditions. Here, the simulation results reveal that the short-range structural order (SRO) in the $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquid is very strong. About 80% of the atoms are arranged in crystallinelike SRO motifs. In particular, ~ $70$% of Fe-centered clusters can be classified as either hexagonal-close-packed-like or icosahedrallike SRO motifs. The SRO clusters centered on Fe, Si, or O atoms are strongly intermixed and homogenously distributed throughout the liquid. The atomic structure of the liquid and the fractions of dominant SRO clusters are not sensitive to pressure/temperature used in the simulations except that the SRO of the O-centered clusters is enhanced close to inner core pressures. The O-diffusion coefficient is about two to three times larger than the Fe and Si ions and increases more rapidly in the deeper core regions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Inertia Estimation Under High Penetration of Inverter-based Resources

Many types of renewable energy sources, especially solar and wind generation, are connected to the electric grid through power-electronic-based interfaces (inverters). These inverter-based resources (IBRs) are mechanically decoupled from the grid, which reduces system inertia and thus may compromise its stability and reliability. In this study, we examine the impact of high penetration of IBRs on the power grid’s inertia. To achieve this, we intentionally introduce disturbances into a simulation case study, such as a step load change, to observe and record the system’s frequency responses. This study offers a nuanced understanding of how the integration of IBRs affects grid stability, and it provides essential guidance for future grid management and resilience strategies.

Yadav, Ajay↗
DOE OSTI2025

Decoupling thermal and irradiation effects on grain boundary segregation

Radiation-induced segregation (RIS) is most often measured by peak solute concentration at a boundary. However, this may give an incomplete picture of segregation quantity and phenomena. Radiation-induced and thermal segregation at grain boundaries was investigated in Fe-9.6 at.% Cr after 9 MeV Fe 3+ ion irradiation at 400 °C. The experimental results were compared to kinetic Monte Carlo (KMC) simulations. The study revealed that Cr enrichment (peak segregation) at the grain boundaries was comparable in both the irradiated and non-irradiated conditions, although irradiation resulted in broader segregation profiles in both experiments and simulations, indicating an overall increase in grain boundary segregation due to irradiation. This broadening is attributed to back diffusion into the grain interior. While it is an established phenomenon, this study offers a quantitative evaluation using experimental data and KMC modeling. Further, these results emphasize the importance of analyzing the entire segregation profile and decoupling the thermal and irradiation contributions to solute segregation.

Grain boundary segregation↗
DOE OSTI2025

Effects of high-pressure hydrogen exposure on filler-elastomer adhesion

Elastomers are known to gain enhanced mechanical properties through compounding with nanosized filler particles such as silica or carbon black. Filler dispersion and filler-polymer interfacial strength are key contributing factors to this improvement. The interfacial strength is critical to part lifetime in pressurized gas sealing applications such as O-rings, where weak binding between the filler particle and polymer matrix can lead to internal void structures. With the aim to build a fundamental understanding of precursors to pressurized hydrogen-induced failure in elastomers, we use all-atom molecular dynamics simulations to study the impact of hydrogen oversaturation on filler-polymer interaction strength. We systematically study the interface between a commonly used elastomer, ethylene-propylene-diene monomer (EPDM) and silica by varying gas concentration, crosslink density, and surface chemistry. Our simulations predict that decompression leads to a localization of excess gas near the interface. In conclusion, we demonstrate that this localized gas can weaken interfacial adhesion and quantify the interaction using thermodynamic approaches.

EPDM↗
DOE OSTI2022

Effect of void positioning on the detonation sensitivity of a heterogeneous energetic material

We show although it is well-established that voids profoundly influence the initiation and reaction behaviors of heterogeneous energetic materials such as polymer-bonded explosives (PBX) and propellants, there has been little study of how void location in different constituents in the microstructures of such materials affect the macroscale behavior. Here, we use three-dimensional (3D) mesoscale simulations to study how void placement within the reactive grains versus the polymer binder influences the shock-to-detonation transition (SDT) in a polymer-bonded explosive. The material studied here has a microstructure comprised of 75% PETN (pentaerythritol tetranitrate) grains and 25% HTPB (hydroxyl-terminated polybutadiene) polymer binder by volume. Porosities up to 10% in the form of spherical voids distributed in both the grains and polymer are considered. An Arrhenius reactive burn relation is used to model the chemical kinetics of the PETN grains under shock loading, thereby resolving the heterogeneous detonation behavior of the PBX. The influence of void location on the shock initiation sensitivity of the material is quantitatively ranked by comparing the predicted run distance to detonation (RDD) for each sample. The analysis includes inherent quantification of uncertainties arising from the stochastic variations in the microstructure morphologies and void distributions by using statistically equivalent microstructure sample sets (SEMSS), leading to probabilistic formulations for the RDD as a function of shock pressure. The calculations reveal that the location of voids in the composite microstructure significantly affects the RDD. Specifically, voids exclusively within the grains cause the PBX to be more sensitive (having shorter RDD) than voids in the polymer binder. Unique probabilistic relationships are derived to map the probability of observing RDD for each void location material case, allowing for prediction of initiation behavior anywhere in the shock pressure – RDD space. These findings agree with trends reported in the literature.

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

How negative feedback and the ambient environment limit the influence of recombination in common envelope evolution

ABSTRACT We perform 3D hydrodynamical simulations to study recombination and ionization during the common envelope (CE) phase of binary evolution, and develop techniques to track the ionic transitions in time and space. We simulate the interaction of a $2\, \mathrm{M_\odot }$ red giant branch primary and a $1\, \mathrm{M_\odot }$ companion modelled as a particle. We compare a run employing a tabulated equation of state (EOS) that accounts for ionization and recombination, with a run employing an ideal gas EOS. During the first half of the simulations, ∼15 per cent more mass is unbound in the tabulated EOS run due to the release of recombination energy, but by simulation end the difference has become negligible. We explain this as being a consequence of (i) the tabulated EOS run experiences a shallower inspiral and hence smaller orbital energy release at late times because recombination energy release expands the envelope and reduces drag, and (ii) collision and mixing between expanding envelope gas, ejecta and circumstellar ambient gas assists in unbinding the envelope, but does so less efficiently in the tabulated EOS run where some of the energy transferred to bound envelope gas is used for ionization. The rate of mass unbinding is approximately constant in the last half of the simulations and the orbital separation steadily decreases at late times. A simple linear extrapolation predicts a CE phase duration of ${\sim}2\, {\rm yr}$, after which the envelope would be unbound.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Projecting California Light-Duty Vehicle Attributes (2019-2035)

In this report, we estimated a set of future vehicle attribute scenarios for new light-duty vehicles for the California market using the Automotive Deployment Options Projection Tool (ADOPT). ADOPT starts simulations with all existing vehicle makes and models and endogenously creates new vehicle models over time based on assumed technology improvements and market conditions. For this study, we simulated four scenarios with varied technology, policy, and electric vehicle infrastructure assumptions. We aggregated simulation results into up to 30 vehicle classes (covering size and price classifications) for each powertrain, for six powertrains. We simulated model years 2019 to 2035. We present results for four vehicle attributes: vehicle acceleration, fuel economy (including electric and gasoline for plug-in hybrid electric vehicles [PHEVs]), vehicle range, and vehicle purchase price. We also present results showing the estimated number of vehicle models available in each vehicle class over time. In the Mid scenario, which contains conditions between our most conservative and most optimistic assumptions for emerging electric and hydrogen technologies, we observe improvements in acceleration, range (for battery-electric vehicles [BEVs]), and vehicle purchase price for electric and hydrogen powertrains. In some cases, we project that consumer preferences lead to trade-offs between vehicle attributes, such as reduced fuel economy in exchange for improved acceleration. Conventional vehicles show modest improvements in fuel economy in these scenarios due to the high numbers of BEV and PHEV sales, which reduce the improvements required in conventional vehicles to meet fleet fuel economy standards. In scenarios with advanced technology assumptions (including reduced battery price and improved energy density), we observe further improvements in BEV range and price relative to the improvements in the Mid scenario.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Collectivity in ultra-peripheral heavy-ion collisions

We present full (3+1 )D dynamical simulations to study collective behavior in ultra-peripheral nucleus-nucleus collisions (UPC) at the Large Hadron Collider (LHC) with the 3DGlauber+MUSIC+UrQMD framework [1, 2]. By extrapolating from asymmetric p+Pb collisions, we simulate a quasi-real photon γ* interacting with the Pb nucleus in an ultra-peripheral collision at the LHC, assuming strong final-state effects. We study the elliptic flow hierarchy between p+Pb and γ*+Pb collisions, which is dominated by the difference in longitudinal flow decorrelations. Our theoretical framework provides a quantitative tool to study collectivity in small asymmetric collision systems at current and future collider experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Experimental and theoretical study of energy transfer in a chromophore triad: What makes modeling dynamics successful?

Simulation of electronic dynamics in realistically large molecular systems is a demanding task that has not yet achieved the same level of quantitative prediction already realized for its static counterpart. This is particularly true for processes occurring beyond the Born–Oppenheimer regime. Non-adiabatic molecular dynamics (NAMD) simulations suffer from two convoluted sources of error: numerical algorithms for dynamics and electronic structure calculations. While the former has gained increasing attention, particularly addressing the validity of ad hoc methodologies, the effect of the latter remains relatively unexplored. Indeed, the required accuracy for electronic structure calculations to reach quantitative agreement with experiment in dynamics may be even more strict than that required for static simulations. In this study, we address this issue by modeling the electronic energy transfer in a donor–acceptor–donor (D–A–D) molecular light harvesting system using fewest switches surface hopping NAMD simulations. In the studied system, time-resolved experimental measurements deliver complete information on spectra and energy transfer rates. Subsequent modeling shows that the calculated electronic transition energies are “sufficiently good” to reproduce experimental spectra but produce over an order of magnitude error in simulated dynamical rates. We further perform simulations using artificially shifted energy gaps to investigate the complex relationship between transition energies and modeled dynamics to understand factors affecting non-radiative relaxation and energy transfer rates.

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

The DESI One-Percent Survey: exploring a generalized SHAM for multiple tracers with the UNIT simulation

We perform SubHalo Abundance Matching (SHAM) studies on UNIT simulations with {σ, V ceil , v smear }-SHAM and {σ, V ceil , $f$ sat }-SHAM. They are designed to reproduce the clustering on 5–30 h -1 Mpc of luminous red galaxies (LRGs), emission-line galaxies (ELGs), and quasi-stellar objects (QSOs) at 0.4 < z < 3.5 from DESI (Dark Energy Spectroscopic Instrument) One Percent Survey. V ceil is the incompleteness of the massive host (sub)haloes and is the key to the generalized SHAM. v smear models the clustering effect of redshift uncertainties, providing measurements consistent with those from repeat observations. A free satellite fraction $f$ sat is necessary to reproduce the clustering of ELGs. We find ELGs present a more complex galaxy–halo mass relation than LRGs reflected in their weak constraints on σ. LRGs, QSOs, and ELGs show increasing V ceil values, corresponding to the massive galaxy incompleteness of LRGs, the quenched star formation of ELGs and the quenched black hole accretion of QSOs. For LRGs, a Gaussian v smear presents a better profile for subsamples at redshift bins than a Lorentzian profile used for other tracers. The impact of the statistical redshift uncertainty on ELG clustering is negligible. The best-fitting satellite fraction for DESI ELGs is around 4 per cent, lower than previous estimations for ELGs. The mean halo mass log 10 ($\langle$M vir $\rangle$) in h -1 M ⊙ for LRGs, ELGs, and QSOs are 13.16 ± 0.01, 11.90 ± 0.06, and 12.66 ± 0.45, respectively. Our generalized SHAM algorithms facilitate the production of multitracer galaxy mocks for cosmological tests.

79 ASTRONOMY AND ASTROPHYSICS↗