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

The homogeneous mixing of MgO and H 2 O at extreme conditions

Investigating water worlds presents a unique opportunity to understand the fundamental processes of planetary formation and evolution. One key aspect is characterizing the interactions between water and rock under the pressures and temperatures present within these worlds. Investigating the conditions for the homogeneous mixing of these materials is imperative to characterizing bulk properties and evolution of water-rich exoplanets. Here, in this study, we use density functional molecular dynamics simulations to study MgO-H 2 O mixtures at high pressure–temperature conditions where H 2 O occurs in solid, superionic, or liquid form. MgO, the representative rocky material, can be either solid or liquid. We start from 500 K at 120 GPa, increasing the temperature step by step up to 8000 K. By inspection, we determine the temperature at which MgO-HO homogeneously mix in our simulations. At 6000 K and 174 GPa is when we find the system to homogeneously mix. This heat-until-it-mixes approach provides us with an upper bound on the temperature for the mixing of MgO and H 2 O. We find that homogeneous mixing occurs at sufficiently low temperatures to be relevant for the collisional growth of a water-rich planet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Control coordination in inverter‐based microgrids using AoI‐based 5G schedulers

Abstract A coordinated set point automatic adjustment with correction enabled (C‐SPAACE) framework that uses 5G communication for real‐time control coordination between inverter‐based resources (IBR) in microgrids is proposed. Utilising slicing capability, 5G offers low‐latency communication to C‐SPAACE under normal conditions. However, given the multitude of power grid use cases, a certain 5G slice for C‐SPAACE may have access only to limited radio spectrum resources, which if not managed well, greatly undermines the communication needs of C‐SPAACE framework. Thus, optimally scheduling the available spectrum resources among IBRs in a sliced 5G network‐based C‐SPAACE framework becomes a critical problem. To address this issue, the authors utilise a novel age of information (AoI) metric and designs an AoI‐based 5G scheduler to provide low‐latency communication to C‐SPAACE. Following this, a co‐simulation environment is designed using PSCAD/EMTDC and Python to simulate a microgrid supported by 5G communication. Time‐domain simulation case studies are performed using the proposed co‐simulation environment to evaluate the performance of C‐SPAACE using 5G with both AoI‐based and other baseline (non‐AoI) schedulers.

Choudhury, Biplav↗
DOE OSTI2024

Analysis of single-mode Richtmyer–Meshkov instability using high-order incompressible vorticity-streamfunction and shock-capturing simulations

In this study, two- and three-dimensional simulation results obtained using a new high-order incompressible, variable-density vorticity–streamfunction (VS) method and data from previous ninth-order weighted essentially nonoscillatory (WENO) shock-capturing simulations are used to investigate the nonlinear dynamics of single-mode Richtmyer–Meshkov instability using a model of a Mach 1.3 air(acetone)/SF6 shock tube experiment. A comparison of the density fields from both simulations with the experimental images demonstrates very good agreement in the large-scale structure with both methods but differences in the small-scale structure. The WENO method captures the small-scale disordered structure observed in the experiment, while the VS method partially captures such structure and yields a strong rotating core. The perturbation amplitude growth from the simulations generally agrees well with the experiment. The simulation bubble and spike amplitudes agree well at early times. At later times, the WENO bubble amplitude is smaller than the VS amplitude and vice versa for the spike amplitude. The predictions of nonlinear single-mode instability growth models are shown to agree with the simulation amplitudes at early-to-intermediate times but underpredict the amplitudes at later times in the nonlinear regime. Visualizations of the mass fraction and enstrophy isosurfaces, velocity and vorticity fields, and baroclinic vorticity production and vortex stretching terms from the three-dimensional simulations indicate that, with the exception of the small-scale structure within the rollups, the VS and WENO results are in good agreement.

42 ENGINEERING↗
DOE OSTI2020

Detection of DoS Attacks Using ARFIMA Modeling of GOOSE Communication in IEC 61850 Substations

Integration of Information and Communication Technology (ICT) in modern smart grids (SGs) offers many advantages including the use of renewables and an effective way to protect, control and monitor the energy transmission and distribution. To reach an optimal operation of future energy systems, availability, integrity and confidentiality of data should be guaranteed. Research on the cyber-physical security of electrical substations based on IEC 61850 is still at an early stage. In the present work, we first model the network traffic data in electrical substations, then, we present a statistical Anomaly Detection (AD) method to detect Denial of Service (DoS) attacks against the Generic Object Oriented Substation Event (GOOSE) network communication. According to interpretations on the self-similarity and the Long-Range Dependency (LRD) of the data, an Auto-Regressive Fractionally Integrated Moving Average (ARFIMA) model was shown to describe well the GOOSE communication in the substation process network. Based on this ARFIMA-model and in view of cyber-physical security, an effective model-based AD method is developed and analyzed. Two variants of the statistical AD considering statistical hypothesis testing based on the Generalized Likelihood Ratio Test (GLRT) and the cumulative sum (CUSUM) are presented to detect flooding attacks that might affect the availability of the data. Our work presents a novel AD method, with two different variants, tailored to the specific features of the GOOSE traffic in IEC 61850 substations. The statistical AD is capable of detecting anomalies at unknown change times under the realistic assumption of unknown model parameters. The performance of both variants of the AD method is validated and assessed using data collected from a simulation case study. We perform several Monte-Carlo simulations under different noise variances. The detection delay is provided for each detector and it represents the number of discrete time samples after which an anomaly is detected. In fact, our statistical AD method with both variants (CUSUM and GLRT) has around half the false positive rate and a smaller detection delay when compared with two of the closest works found in the literature. Our AD approach based on the GLRT detector has the smallest false positive rate among all considered approaches. Whereas, our AD approach based on the CUSUM test has the lowest false negative rate thus the best detection rate. Depending on the requirements as well as the costs of false alarms or missed anomalies, both variants of our statistical detection method can be used and are further analyzed using composite detection metrics.

IEC 61850 electrical substations↗
DOE OSTI2020

Interactions Between the Amazonian Rainforest and Cumuli Clouds: A Large‐Eddy Simulation, High‐Resolution ECMWF, and Observational Intercomparison Study

Abstract The explicit coupling at meter and second scales of vegetation's responses to the atmospheric‐boundary layer dynamics drives a dynamic heterogeneity that influences canopy‐top fluxes and cloud formation. Focusing on a representative day during the Amazonian dry season, we investigate the diurnal cycle of energy, moisture and carbon dioxide at the canopy top, and the transition from clear to cloudy conditions. To this end, we compare results from a large‐eddy simulation technique, a high‐resolution global weather model, and a complete observational data set collected during the GoAmazon14/15 campaign. The overall model‐observation comparisons of radiation and canopy‐top fluxes, turbulence, and cloud dynamics are very satisfactory, with all the modeled variables lying within the standard deviation of the monthly aggregated observations. Our analysis indicates that the timing of the change in the daylight carbon exchange, from a sink to a source, remains uncertain and is probably related to the stomata closure caused by the increase in vapor pressure deficit during the afternoon. We demonstrate quantitatively that heat and moisture transport from the subcloud layer into the cloud layer are misrepresented by the global model, yielding low values of specific humidity and thermal instability above the cloud base. Finally, the numerical simulations and observational data are adequate settings for benchmarking more comprehensive studies of plant responses, microphysics, and radiation.

Vilà‐Guerau de Arellano, J.↗
DOE OSTI2026

Impact of vertical and seasonal variation in leaf traits on simulating soybean canopy photosynthesis via 1D and 3D modeling

Accurate modeling of photosynthesis is crucial for predicting crop productivity and quantifying the carbon cycle in agroecosystems. Leaf traits are essential inputs for modeling canopy photosynthesis. Yet, many existing models still use fixed plant functional type (PTF)-based values to parameterize leaf traits under a big-leaf or two-big-leaf assumption, neglecting their vertical profiles and seasonal changes. This simplification may introduce significant uncertainties in estimating gross primary productivity (GPP). In this study, we simulated soybean GPP and tested the effects of vertical and seasonal variation in three key leaf photosynthetic traits: the maximum carboxylation rate at 25 °C (Vcmax 25 ), leaf chlorophyll content (LCC), and leaf mass per area (LMA) in the 1D-SCOPE and 3D-Helios models. Weekly field measurements were conducted during the growing season of 2024 to support the simulation. We designed ten leaf trait parameterization schemes by incorporating different combinations of vertical profiles and seasonal changes, while assuming homogeneous canopy architecture in both models. Our results revealed that Vcmax 25 vertical and seasonal variation had the strongest influence on simulated GPP in both 1D and 3D models, while LCC and LMA effects were minimal. Particularly, the scheme with an empirically parameterized Vcmax 25 profile achieved comparable performance to the scheme with the measured Vcmax 25 profile. Both 1D-SCOPE and 3D-Helios accurately modeled GPP (SCOPE: R 2 = 0.87, Bias = 0.55 µmol m⁻² s⁻¹; Helios: R 2 = 0.9, Bias = 0.22 µmol m⁻² s⁻¹) under the most complex scheme, and their responses to vertical and seasonal variation in leaf traits were consistent, demonstrating the robustness of our findings. Based on our findings, we propose a scalable framework for parameterizing leaf traits to improve GPP simulations. This study contributes to improving the representation of leaf trait dynamics in canopy-level photosynthesis models, potentially enhancing our ability to predict crop productivity and understand agroecosystem carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Bead–Spring Simulation of Ionomer Melts—Studying the Effects of Chain-Length and Associating Group Fraction on Equilibrium Structure and Extensional Flow Behavior

Ionomers are associative polymers with diverse applications ranging from selective membranes and high-performance adhesives to abrasion- and chemical-resistant coatings, insulation layers, vacuum packaging, and foamed sheets. Within equilibrium melt, the ionic or associating groups are known to form thermally reversible, associative clusters whose presence can significantly affect the system’s mechanical, viscoelastic, and transport properties. It is, thus, of great interest to understand how to control such clusters’ size distribution, shape, and stability through the designed choice of polymer architecture and the ionic groups’ fraction, arrangement, and interaction strength. In this work, we represent linear associating polymers using a Kremer–Grest type bead–spring model and perform large-scale MD simulations to explore the effect of polymer chain-length (l) and fraction (f s ) of randomly placed associating groups on the size distribution and stability of formed clusters. We consider different chain-lengths (below and above entanglement), varying fractions of associating groups (represented by ‘sticky’ beads) between 5 and 20%, and a fixed sticky–sticky nonbond interaction strength of four times that between regular non-associating beads. For all melts containing associating groups the equilibrium structure factor S(q) displays a signature ionomer peak at low wave vector q whose intensity increases with increasing f s and l. The average cluster size Nc increases with f s . However, the effect of chain-length on Nc appears to be pronounced only at higher values of f s . Under extensional flows, the computed stress (and viscosity) is higher at higher f s and l regardless of strain rate. Beyond a critical strain rate, we observe fragmentation of the associative clusters, which has interesting effects on the stress/viscous response.

36 MATERIALS SCIENCE↗
DOE OSTI2024

How Does Vertical Wind Shear Influence Updraft Characteristics and Hydrometeor Distributions in Supercell Thunderstorms?

Abstract Vertical wind shear is known to affect supercell thunderstorms by displacing updraft hydrometeor mass downshear, thereby facilitating the storms’ longevity. Shear also impacts the size of supercell updrafts, with stronger shear leading to wider, less dilute, and stronger updrafts with likely greater hydrometeor production. To more clearly define the role of shear across different vertical layers on hydrometeor concentrations and displacements relative to supercell updrafts, a suite of idealized numerical model simulations of supercells was conducted. Shear magnitudes were systematically varied across the 0–1, 1–6, and 6–12 km AGL layers, while the thermodynamic environment was held fixed. Simulations show that as shear magnitude increases, especially from 1 to 6 km, updrafts become wider and less dilute with an increase in hydrometeor loading, along with an increase in the low-level precipitation area/rate and total precipitation accumulation. Even with greater updraft hydrometeor loading amid stronger shear, updrafts are more intense in stronger shear simulations due to larger thermal buoyancy owing to wider, less dilute updraft cores. Furthermore, downshear hydrometeor displacements are larger in environments with stronger 1–6-km shear. In contrast, there is relatively less sensitivity of hydrometeor concentrations and displacements to variations in either 0–1- or 6–12-km shear. Results are consistent across free tropospheric relative humidity sensitivity simulations, which show an increase in updraft size and hydrometeor mass with increasing free tropospheric relative humidity owing to a reduction in entrainment-driven dilution for wider updrafts in moister environments. Significance Statement Rotating thunderstorms, known as supercells, are able to persist for multiple hours. One common explanation is that large changes in wind speed and/or direction with height, or shear, transport rain/hail away from supercell updrafts, supporting their maintenance. The strong shear within supercell environments, however, may also lead to greater rail/hail amounts, thereby leading to weaker storms due to this extra mass of water/ice within updrafts. Furthermore, the impact of shear across different height layers on supercell rain/hail characteristics has not been thoroughly investigated. In this study, computer simulations of supercells were conducted to determine that shear occurring between 1 and 6 km above ground level has a large impact on rain/hail distribution in supercells and that stronger shear in this layer leads to wider/stronger supercells with greater rain/hail accumulations at the surface. Additionally, some of the extra mass of water/ice is transported farther away from updrafts due to the stronger environmental storm-relative winds.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Comprehensive Analysis of Model Parameter Uncertainty Influence on Evaluation of HVAC Operation to Mitigate Indoor Virus: A Case Study for an Office Building in a Cold and Dry Climate

Simulation-based studies of HVAC operation to mitigate indoor virus have been conducted to understand tradeoffs between indoor air quality (IAQ) and energy consumption. However, the influence of model parameter uncertainty in these studies has not been systematically quantified, which is critical when providing guidance to building operators. To address this gap, we identify 20 model parameters for a typical medium office building system in a cold and dry climate that can influence IAQ and energy consumption for indoor virus scenarios. Next, the distributions of the parameter values are estimated from literature and a set of simulation samples for three representative days is generated that simultaneously sample all of the parameters from their distributions. Two HVAC virus mitigation strategies are studied: increased filtration using MERV 13 filters and increased ventilation, via supply of 100% outdoor air into buildings. The model parameter uncertainty leads to significant variability in the results, particularly for the IAQ because of the highly uncertain virus generation rate. Use of 100% outdoor air can be beneficial for some uncertain scenarios on the hot day, but shows less IAQ improvement and/or significant energy increases on the other days. The virus removal efficiency and pressure drop of the HVAC filter, fan efficiency, and internal heat gain are the most important parameters to determining the tradeoffs of the two strategies. Our results demonstrate how this model parameter uncertainty analysis methodology can provide practical guidance to building operators.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Enhancing Neutrino Event Reconstruction with Pixel-Based 3D Readout for Liquid Argon Time Projection Chambers

In this paper we explore the potential improvements in neutrino event reconstruction that a 3D pixelated readout could offer over a 2D projective wire readout for liquid argon time projection chambers. We simulate and study events in two generic, idealized detector configurations for these two designs, classifying events in each sample with deep convolutional neural networks to compare the best 2D results to the best 3D results. In almost all cases we find that the 3D readout provides better reconstruction efficiency and purity than the 2D projective wire readout, with the advantages of 3D being particularly evident in more complex topologies, such as electron neutrino charged current events. We conclude that the use of a 3D pixelated detector could significantly enhance the reach and impact of future liquid argon TPC experiments physics program, such as DUNE.

Adams, C↗
DOE OSTI2021

Intermediate scattering functions of a rigid body monoclonal antibody protein in solution studied by dissipative particle dynamic simulation

In the past decade, there was increased research interest in studying internal motions of flexible proteins in solution using Neutron Spin Echo (NSE) as NSE can simultaneously probe the dynamics at the length and time scales comparable to protein domain motions. However, the collective intermediate scattering function (ISF) measured by NSE has the contributions from translational, rotational, and internal motions, which are rather complicated to be separated. Widely used NSE theories to interpret experimental data usually assume that the translational and rotational motions of a rigid particle are decoupled and independent to each other. To evaluate the accuracy of this approximation for monoclonal antibody (mAb) proteins in solution, dissipative particle dynamic computer simulation is used here to simulate a rigid-body mAb for up to about 200 ns. The total ISF together with the ISFs due to only the translational and rotational motions as well as their corresponding effective diffusion coefficients is calculated. The aforementioned approximation introduces appreciable errors to the calculated effective diffusion coefficients and the ISFs. For the effective diffusion coefficient, the error introduced by this approximation can be as large as about 10% even though the overall agreement is considered reasonable. Thus, we need to be cautious when interpreting the data with a small signal change. In addition, the accuracy of the calculated ISFs due to the finite computer simulation time is also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Influence of defects on the femtosecond laser damage resistance of multilayer dielectric gratings

Multilayer dielectric (MLD) gratings with high diffraction efficiency and a high laser-induced damage (LID) threshold for pulse compressors are key to scaling the peak and average power of chirped pulse amplification lasers. However, surface defects introduced by manufacturing, storage, and handling processes can reduce the LID resistance of MLD gratings and impact the laser output. The underlying mechanisms of such defect-initiated LID remain unclear, especially in the femtosecond regime. In this Letter, we model dynamic processes in interactions of a 20-fs near-infrared (NIR) laser pulse and a MLD grating design in the presence of cylindrically symmetrical nodules and particle contaminants and cracks at the surface. Utilizing a dynamic model based on a 2D finite difference in time domain (FDTD) field solver coupled with photoionization, electron collision, and refractive index modification, we study the simulation results for the damage site distribution initiated by defects of various types and sizes and its impact on the LID threshold of the grating design.

Zhang, Simin (ORCID:0000000174948230)↗
DOE OSTI2023

Enhancing strontium clock atom interferometry using quantum optimal control

Strontium clock atom interferometry is a promising new technology, with multiple experiments under development around the world to explore its potential for dark matter and gravitational wave detection. In these detectors, large momentum transfer using sequences of many laser pulses is necessary, and thus high-fidelity pulses are important since small errors become magnified. Quantum optimal control (QOC) is a framework for developing control pulse waveforms that achieve high fidelity and are robust against experimental imperfections. Resonant single-photon transitions using the narrow clock transition of strontium involve significantly different quantum dynamics than more established atom interferometry methods based on far-detuned two-photon Raman or Bragg transitions, which leads to new opportunities and challenges when applying QOC. Here, we study in simulation QOC pulses for strontium clock interferometry and demonstrate their advantage over basic square pulses (primitive pulses) and composite pulses in terms of robustness against multiple noise channels. In conclusion, this could improve the scale of large momentum transfer in Sr clock interferometers, paving the way to achieving these scientific goals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Surface enrichment dictates block copolymer orientation

Orientation of block copolymer (BCP) morphology in thin films is critical to applications as nanostructured coatings. Despite being well-studied, the ability to control BCP orientation across all possible block constituents remains challenging. Here, in this study, we deploy coarse-grained molecular dynamics simulations to study diblock copolymer ordering in thin films, focusing on chain makeup, substrate surface energy, and surface tension disparity between the two constituent blocks. We explore the multi-dimensional parameter space of ordering using a machine-learning approach, where an autonomous loop using a Gaussian process (GP) control algorithm iteratively selects high-value simulations to compute. The GP kernel was engineered to capture known symmetries. The trained GP model serves as both a complete map of system response, and a robust means of extracting material knowledge. We demonstrate that the vertical orientation of BCP phases depends on several counter-balancing energetic contributions, including entropic and enthalpic material enrichment at interfaces, distortion of morphological objects through the film depth, and of course interfacial energies. BCP lamellae are found more resistant to these effects, and thus more robustly form vertical orientations across a broad range of conditions; while BCP cylinders are found to be highly sensitive to surface tension disparity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Simulations with machine learning potentials identify the ion conduction mechanism mediating non-Arrhenius behavior in LGPS

Abstract Li 10 Ge(PS 6 ) 2 (LGPS) is a highly concentrated solid electrolyte, in which Coulombic repulsion between neighboring cations is hypothesized as the underlying reason for concerted ion hopping, a mechanism common among superionic conductors such as Li 7 La 3 Zr 2 O 12 (LLZO) and Li 1.3 Al 0.3 Ti 1.7 (PO 4 ) 3 (LATP). While first principles simulations using molecular dynamics (MD) provide insight into the Li + transport mechanism, historically, there has been a gap in the temperature ranges studied in simulations and experiments. Here, we used a neural network potential trained on density functional theory (DFT) simulations, to run up to 40-nanosecond long MD simulations at DFT-like accuracy to characterize the ion conduction mechanisms across a range of temperatures that includes previous simulations and experimental studies. We have confirmed a Li + sublattice phase transition in LGPS around 400 K, below which the ab -plane diffusivity D a b ∗ is drastically reduced. Concomitant with the sublattice phase transition near 400 K, there is less cation-cation (cross) correlation, as characterized by Haven ratios closer to 1, and the vibrations in the system are more harmonic at lower temperature. Intuitively, at high temperature, the collection of vibrational modes may be sufficient to drive concerted ion hops. However, near room temperature, the vibrational modes available may be insufficient to overcome electrostatic repulsion, thus resulting in less correlated ion motion and comparatively slower ion conduction. Such phenomena of a sublattice phase transition, below which concerted hopping plays a less significant role, may be extended to other highly concentrated solid electrolytes such as LLZO and LATP.

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

In situ observation of the magnetization configuration and reversal in cylindrical nanowires

Curvilinear magnetic structures often have unique magnetic behavior compared to their rectilinear counterparts. This is due to the unique curvilinear boundary conditions as well as the curvature induced Dzyaloshinskii–Moriya-like interaction and the curvature induced anisotropy. The effects of a curvilinear geometry are best studied in 3D structures, where the curvature can have a significant spatial extent. Of these 3D structures, the simplest structure to study is the cylindrical nanowire. Here, we have simulated the magnetization reversal in cylindrical NiFe nanowires and present in situ Lorentz TEM images to support the findings of the simulations. We studied the domain formation and reversal of nanowires with two distinct diameters that give rise to a different reversal behavior. We have, thus, found that the zero-field magnetization configuration in these wires can take on a double helix type of configuration. The reversal in these structures then proceeds through the winding and unwinding of these helical configurations rather than through domain wall propagation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗