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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A direct numerical simulation of Jet A flame kernel quenching

The safe operation of aeronautical engines requires an understanding of flame ignition, propagation and extinction. In this study, direct numerical simulations are performed using a 29 species reduced chemical mechanism for jet fuel surrogate Jet A to understand the flame quenching process. Here, initially laminar spherical flames of varying sizes and equivalence ratios are subject to an identical periodic domain of decaying and isotropic high intensity turbulence with a turbulent Reynolds number of 2400. All cases become quenched, except for the larger kernel with lower Karlovitz number. An analysis of the flame structure shows broadened preheat zone, flame shortening on the product side, differential species diffusion and partial fuel pyrolysis in the fresh mixture. Two extinction mechanisms are identified arising from flame shortening and high flame stretch. Flame shortening occurs due to turbulence-chemistry interactions that resemble the flame–flame interaction in a laminar counterflow reactant-to-reactant configuration, which contorts and breaks up the ignition kernel. Flame stretch is a local effect that attenuates the heat release rate and causes the flame to retreat towards the product mixtures, similar to what has been observed for reactant-to-product laminar counterflow flames. Chemical explosive mode analysis was also performed to quantify the flame structure and local combustion mode. The diffusion–reaction balance in pinched-off flame islands favors extinction of these smaller structures, while auto-ignition modes are observed within the flame kernel after fresh mixture is engulfed and preheated in the product kernel. Statistics of the density-weighted displacement speed conditional on local combustion mode indicates strong correlation between the local extinction mode and negative displacement speed. The local balance between diffusion and reaction ultimately determines the propensity for local extinction in both laminar and turbulent flames, the extent of which has an impact on global flame propagation.

42 ENGINEERING↗
DOE OSTI2023

MBX: A many-body energy and force calculator for data-driven many-body simulations

Many-Body eXpansion (MBX) is a C++ library that implements many-body potential energy functions (PEFs) within the “many-body energy” (MB-nrg) formalism. MB-nrg PEFs integrate an underlying polarizable model with explicit machine-learned representations of many-body interactions to achieve chemical accuracy from the gas to the condensed phases. MBX can be employed either as a stand-alone package or as an energy/force engine that can be integrated with generic software for molecular dynamics and Monte Carlo simulations. MBX is parallelized internally using Open Multi-Processing and can utilize Message Passing Interface when available in interfaced molecular simulation software. In this study, MBX enables classical and quantum molecular simulations with MB-nrg PEFs, as well as hybrid simulations that combine conventional force fields and MB-nrg PEFs, for diverse systems ranging from small gas-phase clusters to aqueous solutions and molecular fluids to biomolecular systems and metal-organic frameworks.

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

Simulations of multicomponent relativistic thermalization

Multicomponent relativistic fluids have been studied for decades. However, simulating the dynamics of the particles and fluids in such a mixture has been a challenge due to the fact that such simulations are computationally expensive in three spatial dimensions. Here, we report on the development and application of a multidimensional relativistic Monte Carlo code to explore the thermalization process in a relativistic multicomponent environment in a computationally inexpensive way. As an illustration we simulate the fully relativistic three-dimensional Brownian-motion-like solution to the thermalization of a high-mass particle (proton) in a bath of relativistic low-mass particles (electrons). Here, we follow the thermalization and ultimate equilibrium distribution of the Brownian-like particle as can happen in the cosmic plasma during big-bang nucleosynthesis. We also simulate the thermalization of energetic particles injected into the plasma as can occur, for example, by the decay of massive unstable particles during the big bang.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

WRF-simulated low-level jets over Iowa: characterization and sensitivity studies

Abstract. Output from 6 months of high-resolution simulations with the Weather Research and Forecasting (WRF) model are analyzed to characterize local low-level jets (LLJs) over Iowa for winter and spring in the contemporary climate. Low-level jets affect rotor plane aerodynamic loading, turbine structural loading and turbine performance, and thus accurate characterization and identification are pertinent. Analyses using a detection algorithm wherein the wind speed above and below the jet maximum must be below 80 % of the jet wind speed within a vertical window of approximately 20–530 m a.g.l. (above ground level) indicate the presence of an LLJ in at least one of the 14 700 4 km×4 km grid cells over Iowa on 98 % of nights. Nocturnal LLJs are most frequently associated with stable stratification and low turbulent kinetic energy (TKE) and hence are more frequent during the winter months. The spatiotemporal mean LLJ maximum (jet core) wind speed is 9.55 m s−1, and the mean height is 182 m. Locations of high LLJ frequency and duration across the state are seasonally varying, with a mean duration of 3.5 h. The highest frequency occurs in the topographically complex northwest of the state in winter and in the flatter northeast of the state in spring. Sensitivity of LLJ characteristics to the (i) LLJ definition and (ii) vertical resolution at which the WRF output is sampled is examined. LLJ definitions commonly used in the literature are considered in the first sensitivity analysis. These sensitivity analyses indicate that LLJ characteristics are highly variable with definition. Use of different definitions identifies both different frequencies of LLJs and different LLJ events. Further, when the model output is down-sampled to lower vertical resolution, the mean jet core wind speed height decreases, but spatial distributions of regions of high frequency and duration are conserved. Implementation of a polynomial interpolation to extrapolate down-sampled output to full-resolution results in reduced sensitivity of LLJ characteristics to down-sampling.

17 WIND ENERGY↗
DOE OSTI2023

Effects of Lower Troposphere Vertical Mixing on Simulated Clouds and Precipitation Over the Amazon During the Wet Season

Planetary boundary layer (PBL) schemes parameterize unresolved turbulent mixing within the PBL and free troposphere (FT). Previous studies reported that precipitation simulation over the Amazon in South America is quite sensitive to PBL schemes and the exact relationship between the turbulent mixing and precipitation processes is, however, not disentangled. In this study, regional climate simulations over the Amazon in January–February 2019 are examined at process level to understand the precipitation sensitivity to PBL scheme. The focus is on two PBL schemes, the Yonsei University (YSU) scheme, and the asymmetric convective model v2 (ACM2) scheme, which show the largest difference in the simulated precipitation. During daytime, while the FT clouds simulated by YSU dissipate, clouds simulated by ACM2 maintain because of enhanced moisture supply due to the enhanced vertical moisture relay transport process: (a) vertical mixing within PBL transports surface moisture to the PBL top, and (b) FT mixing feeds the moisture into the FT cloud deck. Due to the thick cloud deck over Amazon simulated by ACM2, surface radiative heating is reduced and consequently the convective available potential energy is reduced. As a result, precipitation is weaker from ACM2. Two key parameters dictating the vertical mixing are identified, p , an exponent determining boundary layer mixing and λ , a scale dictating FT mixing. Sensitivity simulations with altered p , λ , and other treatments within YSU and ACM2 confirm the precipitation sensitivity. The FT mixing in the presence of clouds appears most critical to explain the sensitivity between YSU and ACM2.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A pair of early- and late-forming galaxy cluster samples: A novel way of studying halo assembly bias assisted by a constrained simulation

The halo assembly bias, a phenomenon referring to dependencies of the large-scale bias of a dark matter halo other than its mass, is a fundamental property of the standard cosmological model. First discovered in 2005 from the Millennium Run simulation, it has been proven very difficult to be detected observationally, with only a few convincing claims of detection so far. The main obstacle lies in finding an accurate proxy of the halo formation time. In this study, by utilizing a constrained simulation that can faithfully reproduce the observed structures larger than 2 Mpc in the local universe, for a sample of 634 massive clusters at z ≤ 0.12, we found their counterpart halos in the simulation and used the mass growth history of the matched halos to estimate the formation time of the observed clusters. This allowed us to construct a pair of early- and late-forming clusters, with a similar mass as measured via weak gravitational lensing, and large-scale biases differing at the ≈3σ level, suggestive of the signature of assembly bias, which is further corroborated by the properties of cluster galaxies, including the brightest cluster galaxy and the spatial distribution and number of member galaxies. Our study paves a way to further detect assembly bias based on cluster samples constructed purely on observed quantities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Studying Open Quantum Systems Relevant to Chemistry on a Trapped-Ion Quantum Simulator (Final Technical Report)

This project advances the trapped-ion quantum simulator as a versatile platform for studying open quantum system phenomena. We aim to contribute to the emerging quantum simulation toolkits and enable simulation of nanoscale energy processes. Trapped-ion platforms offer unique capabilities: their vibrational motion can be precisely manipulated, measured, and coherently coupled to auxiliary qubits. The vibrational mode can function both as a highly sensitive quantum sensor and a programmable environment bath. Using this platform, we achieved three major outcomes. First, we demonstrated using the vibrational mode as an ultrasensitive probe for testing fundamental physics, including possible nonlinear quantum mechanics effects. Second, we established that these modes can act as controllable baths in which tunable noise and loss can enhance or modify energy-transfer dynamics, providing the experimental preparation toward studying mechanisms relevant to chemical reactions and light-harvesting systems. Third, by introducing controllable nonlinear gain and loss, we showed theoretically how simulations using trapped ions can model vibrationally-assisted energy transport in a non‐Hermitian quantum system comprising a chromophore dimer weakly coupled to a vibrational mode. Exploring the non‐Hermitian dynamics of the whole system including vibrations, we found that energy transfer accompanied by absorption of phonons from a vibrational mode can be significantly enhanced near an exceptional point. This theoretical work on simulation of energy transfer processes in driven non‐Hermitian quantum systems revealed an interesting novel path to study open quantum systems dynamics under conditions of gain and loss. We then further explored the benefits of controllable gain and loss with an experimental realization of quantum analogs of nonlinear oscillators, namely, the van der Pol oscillator. Here we observed mutual synchronization mediated by collective dissipation between two oscillators. In parallel, we explored related quantum networking protocols using the same trapped-ion platform, developing fast, high-fidelity schemes for distributing entanglement. Together, these achievements show that trapped-ion vibrational modes provide a highly programmable and high-fidelity platform for investigating complex dissipative quantum behavior, while enabling new approaches to remote quantum sensing, energy science, and nonlinear quantum dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Modelling Beam Loss Within the LANSCE Proton Storage Ring

Several upgrades are being considered for the proton storage ring (PSR) at the Los Alamos Neutron Science Center (LANSCE) to reduce beam loss and thereby reduce the cooldown period of the PSR following a beam run. First, we have considered an increased beam pipe diameter would reduce beam loss due to beam scraping caused by misalignments and tuning errors. However, this would require increased pole-to-pole gap height within the dipoles and quadrupoles, which would change their effective length and alter their fringe fields. The effect of different magnet gaps on the beam optical parameters and on beam loss was studied using the simulation codes MAD-X and PyORBIT. Second, we are developing a detailed particle tracking model within the framework of the simulation package General Particle Tracer (GPT) for the PSR H - stripping system. This model will be used to study the effect of the stripper foil parameters (position, composition, areal density, depth, etc.) on first-turn losses, where most of the observed beam loss and emittance growth occurs due to foil scattering, foil stripping, and Lorentz stripping. The model implements C++ custom elements created to model foil scattering, foil stripping, and Lorentz stripping, as these are not built-in features of GPT. The preliminary results of the MAD-X and PyORBIT simulations, as well as the GPT simulation model, is described in this tech note.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Adoption of image-driven machine learning for microstructure characterization and materials design: A Perspective

Microstructure characterization enables the development of structure-processing-property relationships critical to several research areas within the broad field of materials science, from alloy design to the assessment of corrosion resistance, and failure analysis. Conventional approaches to material characterization have relied on either qualitative inference by the human ex-pert or software applications that can extract high-level features from images, such as boundary segmentation, average grain diameter, etc. Such approaches rely heavily on subject matter expert user intervention and knowledge of what phases or more generally, what microstructural features, are of interest. The recent surge in the adoption of machine learning techniques to address problems in materials engineering has brought with it an increased interest and application of Image Driven Machine Learning (IDML) approaches. In this work, we review the applications of IDML to the field of materials characterization. A canonical hierarchy of stages is defined, which when put sequentially together completes an IDML study: problem definition, dataset building, model selection and training, model evaluation, and integration with existing instrumentation or simulation workflow. The studies reviewed in this work are analyzed from the perspective of each of these stages. Such a review permits agranular assessment of the field, for example the impact of IDML on materials characterization at the nanoscale, the size of a typical dataset required to train a semantic segmentation model on electron microscopy images, ubiquitousness of transfer learning in the domain, etc. Finally, we discuss the importance of interpretability and explainability in the field of IDML for materials characterization, and provide an overview of two emerging techniques in the field: semantic segmentation and generative adversarial networks.

Baskaran, Arun↗
DOE OSTI2024

Wind-induced torsion of parabolic trough collectors in operation

Parabolic trough collector (PTC) systems, a type of concentrating solar power (CSP), use parabolic mirrors to reflect the sun's rays toward an absorber tube to heat the fluid inside. PTCs track the sun throughout the day and are sensitive to angular misalignment, which reduces optical performance. Torsion of the PTC, exacerbated by wind loading, contributes to angular misalignment, which has been quantified by previous studies using numerical simulations and experimental tests. However, angular misalignment due to torsion and wind loading has not yet been studied in an operational plant, which can experience more extreme and variable conditions. This study characterizes the angular misalignment due to torsion at three rows of the Nevada Solar One CSP plant and investigates the influence of wind loading. These findings reveal substantial torsion along the PTC support structure that depends on wind conditions and orientation. Strong winds perpendicular to the PTCs increase the median angular displacement of the outermost row by up to 18 mrad and increase the standard deviation by up to 8 mrad when compared to weak wind conditions. These findings can inform enhanced PTC design, controls and modeling that will improve the performance, reliability and bankability of future CSP plants.

14 SOLAR ENERGY↗
DOE OSTI2019

Simulation and Experimental Results of Interdigital Capacitor (IDC) Sensors to Monitor Insulation Degradation of Cables

The efficacy of using interdigital capacitor (IDC) sensors to monitor cable insulation aging was studied through electromagnetic simulations and experimental measurements. Simulation results on both planar and conformal dielectric materials indicated that IDCs on thin substrates and with conducting backplanes provide superior performance in terms of higher capacitance. An airgap study showed that even a very small airgap between the sensor and the cable’s insulation results in significant reduction in capacitance; therefore, airgaps should be minimized. Finally, cables that underwent accelerated aging exposure in an oven were tested using fabricated IDC sensors. Measurement results showed a 12.5% increase in capacitance when the IDC was pressed against the insulation versus a 91.3% increase in capacitance when the IDC was pressed against the jackets of the cables for cable specimens ranging from unaged to aged 35 days.

Cable NDE, Interdigital Capacitance↗
DOE OSTI2021

OpenAWSEM with Open3SPN2: A fast, flexible, and accessible framework for large-scale coarse-grained biomolecular simulations

We present OpenAWSEM and Open3SPN2, new cross-compatible implementations of coarse-grained models for protein (AWSEM) and DNA (3SPN2) molecular dynamics simulations within the OpenMM framework. These new implementations retain the chemical accuracy and intrinsic efficiency of the original models while adding GPU acceleration and the ease of forcefield modification provided by OpenMM’s Custom Forces software framework. By utilizing GPUs, we achieve around a 30-fold speedup in protein and protein-DNA simulations over the existing LAMMPS-based implementations running on a single CPU core. We showcase the benefits of OpenMM’s Custom Forces framework by devising and implementing two new potentials that allow us to address important aspects of protein folding and structure prediction and by testing the ability of the combined OpenAWSEM and Open3SPN2 to model protein-DNA binding. The first potential is used to describe the changes in effective interactions that occur as a protein becomes partially buried in a membrane. We also introduced an interaction to describe proteins with multiple disulfide bonds. Using simple pairwise disulfide bonding terms results in unphysical clustering of cysteine residues, posing a problem when simulating the folding of proteins with many cysteines. We now can computationally reproduce Anfinsen’s early Nobel prize winning experiments by using OpenMM’s Custom Forces framework to introduce a multi-body disulfide bonding term that prevents unphysical clustering. Our protein-DNA simulations show that the binding landscape is funneled towards structures that are quite similar to those found using experiments. In summary, this paper provides a simulation tool for the molecular biophysics community that is both easy to use and sufficiently efficient to simulate large proteins and large protein-DNA systems that are central to many cellular processes. These codes should facilitate the interplay between molecular simulations and cellular studies, which have been hampered by the large mismatch between the time and length scales accessible to molecular simulations and those relevant to cell biology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Evaluating and Improving Scale–Awareness of a Convective Parameterization Closure Using Cloud–Resolving Model Simulations of Convection

This study aims to evaluate a revised closure of the Zhang-McFarlane (ZM) scheme for deep convection and to further improve its scale awareness. Output from two cloud-resolving model (CRM) simulations of both midlatitude organized (MC3E) and tropical unorganized (TWP-ICE) convection are used for the evaluation. By averaging the CRM output over different subdomain sizes from 64 km to 4 km, convection-associated fields are obtained to represent the fields at different horizontal resolutions of global climate models (GCMs). The convection parameterization closure is then tested using these fields. Results show that the cloud base mass flux M b determined from the closure has a relatively good relationship with CRM-simulated convection at coarse resolutions for both MC3E and TWP-ICE, but their correlation degrades at higher resolutions, especially when approaching the gray zone scale under 10 km. The fluctuation of CAPE consumption rate per unit cloud base mass flux for different convective events is found to be an important factor affecting the correlation. When averaged over the base domain, the cloud base mass flux determined by the closure is found to increase with decreasing subdomain size, implying that the closure is not scale-aware. To improve its scale-awareness, a modified closure is designed to alleviate the overprediction of convection at high GCM resolutions. Furthermore, results show that the modified closure has the scale-aware ability to suppress the overprediction problem at high GCM resolutions for both MC3E and TWP-ICE.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Challenges for density functional theory in simulating metal–metal singlet bonding: A case study of dimerized VO2

VO2 is renowned for its electric transition from an insulating monoclinic (M1) phase, characterized by V–V dimerized structures, to a metallic rutile (R) phase above 340 K. This transition is accompanied by a magnetic change: the M1 phase exhibits a non-magnetic spin-singlet state, while the R phase exhibits a state with local magnetic moments. Simultaneous simulation of the structural, electric, and magnetic properties of this compound is of fundamental importance, but the M1 phase alone has posed a significant challenge to the density functional theory (DFT). In this study, we show none of the commonly used DFT functionals, including those combined with on-site Hubbard U to treat 3d electrons better, can accurately predict the V–V dimer length. The spin-restricted method tends to overestimate the strength of the V–V bonds, resulting in a small V–V bond length. Conversely, the spin-symmetry-breaking method exhibits the opposite trends. Each of these two bond-calculation methods underscores one of the two contentious mechanisms, i.e., Peierls lattice distortion or Mott localization due to electron–electron repulsion, involved in the metal–insulator transition in VO2. To elucidate the challenges encountered in DFT, we also employ an effective Hamiltonian that integrates one-dimensional magnetic sites, thereby revealing the inherent difficulties linked with the DFT computations.

Chemistry↗
DOE OSTI2022

Hanford Low-Activity Waste Glass Minor Component Concentration Boundary Expansion

This study investigated the effects of increasing concentrations of minor components on simulated Hanford low-activity waste (LAW) glass properties to allow for LAW glass composition envelope expansion. A test matrix of ten glasses was generated by spiking three previously studied simulated Hanford LAW glass compositions with increasing concentrations of minor component mixes up to a total concentration of 2 wt%. In this study, the ten glasses were fabricated, and key properties were measured including crystal formation after canister centerline cooling (CCC), crystallinity as a function of temperature, density, viscosity, electrical conductivity, glass durability using the Vapor Hydration Test (VHT), product consistency using the Product Consistency Test (PCT), and sulfur solubility. The data from these tests are presented in this report.

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

Calibration of Al7075 with Plate-Puncture Predictions

The following details calibration of a material model for Al7075-T6511. This aluminum alloy is commonly used across a host of engineering applications. Owing to its widespread prevalence, there is great benefit in improving simulation predictions for this alloy. In the present effort, a calibration is performed of its elastic-plastic response accounting for both rate and temperature dependence. The calibration is informed by a series of tests that include specimens of different geometries tested at different rates and temperatures. All specimens are derived from the same barstock, 3.5 inches in diameter. The fitted model itself uses an anisotropic, Hill yield surface coupled with a Johnson-Cook hardening model. Failure predictions are had by means of a modified Wilkins failure criterion. Following calibration of the material model, a validation exercise is performed against platepuncture experiments. These experiments include multiple probe shapes, probe diameters, and plate thicknesses. The puncture experiments are replicated in simulation with mesh studies performed to assess uncertainty. Key quantities of interest, notably the absorbed energy up to failure, are compared between simulation and experiment providing a means to assess the suitability of the calibration in puncture simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Surrogate-driven Variance-based Sensitivity Analysis of Thermal Storage Tanks in Integrated Energy Systems

Sensitivity analysis and uncertainty quantification are essential steps for enhancing the accuracy of computational models by identifying and mitigating uncertainties. This study focuses on these steps for the Thermal Energy Delivery System at Idaho National Laboratory, specifically targeting the thermocline tank. Using a Modelica/Dymola simulation model, the study perturbed various design parameters and boundary conditions, including shape factor, porosity, outlet temperature, inlet mass flow rate, and system pressure, to predict and quantify uncertainty in the tank’s ax- ial temperature. A dataset of over 1,000 simulations was generated, and surrogate models were developed using the pyMAISE (Michigan Artificial Intelligence Standard Environment) library, which is an Automatic Machine Learning library for nuclear engineering applications. The optimal model, a feedforward neural network with two hidden layers, achieved an R2 score above 0.99 and a mean absolute error below 1 Kelvin. Sensitivity analyses using Sobol indices and Fourier amplitude sensitivity testing methods on this surrogate model revealed that the inlet mass flow rate at initial timestamps and porosity significantly impacts predicted temperatures across all sensors and time steps.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗