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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,693 records · Page 94

All sources NASA NTRS DOE OSTI

DOE OSTI2022

PUSHing Core-collapse Supernovae to Explosions in Spherical Symmetry. V. Equation of State Dependency of Explosion Properties, Nucleosynthesis Yields, and Compact Remnants

Abstract In this fifth paper of the series, we use the parameterized, spherically symmetric explosion method PUSH to investigate the impact of eight different nuclear equations of state (EOS). We present and discuss the explosion properties and the detailed nucleosynthesis yields, and predict the remnant (neutron star or black hole) for all our simulations. For this, we perform two sets of simulations. First, a complete study of nonrotating stars from 11 to 40 M ⊙ at three different metallicities using the SFHo EOS; and, second, a suite of simulations for four progenitors (16 M ⊙ at three metallicities and 25 M ⊙ at solar metallicity) for eight different nuclear EOS. We compare our predicted explosion energies and yields to observed supernovae and to the metal-poor star HD 84937. We find EOS-dependent differences in the explosion properties and the nucleosynthesis yields. However, when comparing to observations, these differences are not large enough to rule out any EOS considered in this work.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Evaluating Thermostats' Deadbands Using HVAC Hardware-In-the-Loop Experiment for Advanced Control Strategies: Preprint

Smart thermostats have gained significant popularity due to their potential for optimizing energy consumption and enhanced user control while ensuring occupants' comfort. The deadband, also referred to as temperature differential, is defined as the temperature difference between the desired setpoint and upper threshold or lower threshold for the HVAC equipment to turn on. It is a key factor influencing energy efficiency and user satisfaction. This paper presents a comparative analysis of the deadbands of five different smart thermostats, tested with a heat pump, aiming to identify variations in their deadband settings and implications for energy management. The experimental study was conducted using a HVAC hardware-in-theloop (HIL) system that integrates smart thermostats with physical HVAC equipment in a simulated house environment. The study explores the trade-offs between energy efficiency and occupant comfort and highlights how different thermostats participating in demand response event cycle differently based on their deadband settings. The findings offer valuable insights into how selecting the right thermostat or configuring smart thermostat with appropriate deadband settings can be leveraged to enhance demand response capabilities, shift loads effectively and improve operational flexibility in HVAC systems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2021

A generic method for interfacing IEDs using low voltage interfaces to real-time simulators with hardware in the loop

This study presents an interface method (IM) to interconnect real-time simulators (RTSs) without amplifiers for different intelligent electronic devices (IEDs), such as protective relays, smart meters, and other devices. It is a significant topic for the power systems protection community because there are many protection devices with different requirements. These protection devices have pinouts that measure bus voltages/currents and breaker pole states and generate trip/close pulse signals. Most of the emulation test beds with RTSs and HIL need amplifiers to connect IEDs from different product manufacturers because the low-voltage interface levels are not available in some instruction manuals. Additionally, commercial amplifiers have cutoff frequency limitations and have an excessive cost when several IEDs are wired onto test beds. Then, the validation of protection and control systems with RTSs and IEDs in the loop from different vendors without amplifiers results in an unfeasible alternative. In this paper, the IM was based on using an interface box, power source, and RTS to identify IED inputs/outputs. Here, the sequence to recognize IED pinouts and calculate the current/voltage scaling factors is described. The IM was validated satisfactorily with RTS and a protection device (with unknown pinouts) in the loop.

42 ENGINEERING↗

DOE OSTI2021

Novel Data Driven Noise Emulation Framework using Deep Neural Network for Generating Synthetic PMU Measurements

Sensors play a critical role in supporting day-to-day grid operations and they are essential to operator’s decision-making process. Furthermore, sensors and sensor behaviors need to be emulated with grid simulations to perform modeling studies and to design cutting edge power systems applications. Ensuring the accurate behavior of these applications requires accurate emulation of sensors and pertinent signals. However, most grid simulators and modeling tools assume either zero error scenarios or simplistic noise models that may not always correlate to real-world sensors. To address the above issue, this work presents an initial study on the noise characteristics of phasor measurement units (PMUs), along with models for recreating their unique noise signatures. The proposed methods (both analytical and machine-learning-based) provide a substantial increase in a sensor’s model fidelity, a feature that can be leveraged by an end-user application to yield more accurate system representations. The proposed methods were then applied to micro PMU data from the EPFL microgrid campus to extract sensor noise profiles. This data was used to train a deep learning model, which was tested to emulate the noise characteristics present in actual signals. Based on the observed results and the employed data-driven methodology, the proposed methods may be adapted to replicate the behavior of other grid sensors and power new applications capable of detecting sensor degradation and eventual device failures in near real-time.

PMU, noise emulation, synthetic measurements, deep↗

DOE OSTI2022

Investigation of the role of hydrogen molecules in 1D simulation of divertor detachment

The role of neutral and charged hydrogenic molecules in detached regimes of tokamak plasmas is investigated using a simplified 1D numerical model. Using MAST Upgrade like conditions, simulations are implemented to study the rollover of target flux Γ in upstream density scan. It is found that if H 2 and H 2 + are considered in simulations a lower target temperature and a larger upstream density will be required to trigger divertor detachment under the same input power and particle flux, and the critical detachment threshold (the critical ratio of upstream static pressure to the power entering the recycling region) is found to be $^{ρ_{up}}_{ρ_{recl}}$ ~ 8.1 NMW -1 at rollover. Molecule–plasma interactions are found to be as crucial as atom–plasma interactions during divertor detachment, both of which account for the majority of plasma momentum loss in the cases studied here. Further analysis of the momentum loss decomposition shows molecule-plasma elastic collisions dominate molecule-plasma interactions, while molecular charge exchange cannot effectively reduce plasma momentum. In terms of H alpha emission, a strong rise of H alpha signal is found to be due to molecular excitation channels when the upstream density further increases after rollover.

SD1D↗

DOE OSTI2022

Assessing free tropospheric quasi-equilibrium for different GCM resolutions using a cloud-resolving model simulation of tropical convection

Abstract This study examines the free-tropospheric quasi-equilibrium at different global climate model (GCM) resolutions using the simulation of tropical convection by a cloud-resolving model during the Tropical Western Pacific International Cloud Experiment. The simulated dynamic and thermodynamic fields within the model domain are averaged over subdomains of different sizes equivalent to different GCM resolutions. These coarse-grained fields are then used to compute CAPE and its change with time, and their relationships with simulated convection. Results show that CAPE change with time is controlled predominantly by variations of thermodynamic properties in the planetary boundary layer for all subdomain sizes ranging from 64 to 4 km. Lag correlation analysis shows that CAPE generation by the free-tropospheric dynamical advection (dCAPE ls ) leads convective precipitation but is in phase with convective mass flux at 600 mb and 500 mb vertical velocity for all subdomain sizes. However, the correlation coefficients and regression slopes decrease as the subdomain size decreases for subdomain sizes smaller than 16 km. This is probably due to increased randomness of convection and more scale-dependence of the relationships when the subdomain size reaches the grey zone. By examining the sensitivity of the relationships of convection with dCAPE ls to temporal scales in different subdomain size, it shows that the quasi-equilibrium between dCAPE ls and convection holds well for timescales of 30 min or longer at all subdomain sizes. These results suggest that the free tropospheric quasi-equilibrium assumption may still be useable even for GCM resolutions in the grey zone.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Probabilistic constrained Bayesian inversion for transpiration cooling

Abstract To enable safe operations in applications such as rocket combustion chambers, the materials require cooling to avoid material damage. Here, transpiration cooling is a promising cooling technique. Numerous studies investigate possibilities to simulate and evaluate the complex cooling mechanism. One naturally arising question is the amount of coolant required to ensure a safe operation. To study this, we introduce an approach that determines the posterior probability distribution of the Reynolds number using an inverse problem and constraining the maximum temperature of the system under parameter uncertainties. Mathematically, this chance inequality constraint is dealt with by a generalized polynomial chaos expansion of the system. The posterior distribution will be evaluated by different Markov chain Monte Carlo based methods. A novel method for the constrained case is proposed and tested among others on two‐dimensional transpiration cooling models.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

ezAlign: A Tool for Converting Coarse-Grained Molecular Dynamics Structures to Atomistic Resolution for Multiscale Modeling

Soft condensed matter is challenging to study due to the vast time and length scales that are necessary to accurately represent complex systems and capture their underlying physics. Multiscale simulations are necessary to study processes that have disparate time and/or length scales, which abound throughout biology and other complex systems. Herein we present ezAlign, an open-source software for converting coarse-grained molecular dynamics structures to atomistic representation, allowing multiscale modeling of biomolecular systems. The ezAlign v1.1 software package is publicly available for download at github.com/LLNL/ezAlign. Its underlying methodology is based on a simple alignment of an atomistic template molecule, followed by position-restraint energy minimization, which forces the atomistic molecule to adopt a conformation consistent with the coarse-grained molecule. The molecules are then combined, solvated, minimized, and equilibrated with position restraints. Validation of the process was conducted on a pure POPC membrane and compared with other popular methods to construct atomistic membranes. Additional examples, including surfactant self-assembly, membrane proteins, and more complex bacterial and human plasma membrane models, are also presented. By providing these examples, parameter files, code, and an easy-to-follow recipe to add new molecules, this work will aid future multiscale modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Simulation of Reverse Osmosis Membrane Compaction Using Material Point Method (MPM)

Access to fresh drinking water in the future requires effective management of industrial wastewater and water purification from available resources. In this study, we present the simulation methodology and the analysis of a Reverse Osmosis (RO) membrane under various pressure conditions using the Material Point Method (MPM). In contrast to other numerical methods, MPM solves the continuum governing equations on material points in a Lagrangian framework. The method does not require a grid connecting the material points hence making it suitable to simulate large deformations during membrane compaction. The time integration is carried out using the explicit Euler method, while the spatial discretization is performed using linear or cubic-spline shape functions. A series of planar images containing detailed pore structures obtained from X-ray tomography experiments is converted to a three-dimensional collection of material points to simulate the membrane. Compressive loads are applied to the top layer of the membrane to simulate the experiments. The membrane deformation and pore size distribution before and after load application are reported and compared with the experimental measurements. The presentation discusses the numerical methods used, the performance of the solver on high-performance computing machines, and the results of membrane compaction in detail.

ENGINEERING,MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Expanded verification and validation studies of hypersonic aerodynamics with multiple physics-fidelity models

Hypersonic aerothermodynamics is an important domain of modern multiphysics simulation. The Multi-Fidelity Toolkit is a simulation tool being developed at Sandia National Laboratories to predict aerodynamic properties for compressible flows from a range of physics fidelities and computational speeds. These models include the Reynolds-averaged Navier–Stokes (RANS) equations, the Euler equations with momentum-energy integral technique (MEIT), and modified Newtonian aerodynamics with flat-plate boundary layer (MNA+FPBL) equations, and they can be invoked independently or coupled with hierarchical Kriging to interpolate between high-fidelity simulations using lower-fidelity data. However, as with any new simulation capability, verification and validation are necessary to gather credibility evidence. This work describes formal code- and solution-verification activities, as well as model validation with uncertainty considerations. Code verification activities on the MNA+FPBL model build on previous work by focusing on the viscous portion of the model. Viscous quantities of interest are compared against those from an analytical solution for flat-plate, inclined-plate, and cone geometries. The code verification methodology for the MEIT model is also presented. Test setup and results of code verification tests on the laminar and turbulent models within MEIT are shown. Solution-verification activities include grid-refinement studies on simulations that model the HIFiRE-1 wind tunnel experiments. These experiments are used for validation of all model fidelities. A thorough validation comparison with prediction error and uncertainty is also presented. Three additional HIFiRE-1 experimental runs are simulated in this study, and the solution verification and validation work examines the effects of the associated parameter changes on model performance. Finally, a study is presented that compares the computational costs and fidelities from each of the different models.

42 ENGINEERING↗

DOE OSTI2022

Direct Numerical Simulations of Turbulent Flow over Misaligned Traveling Waves

Direct numerical simulations of turbulent flow over prescribed traveling waves in a half-channel flow setup subject to a combined streamwise and spanwise pressure gradient are performed at a friction Reynolds number of 180. The simulations undertaken in this study consider both flow-aligned (including opposing and following) and misaligned waves as we attempt to quantify the effect of wave misalignment on turbulence statistics, the pressure drag force, and the wave attenuation rate. For the simulations, we consider three characteristic wave age values corresponding to slow-, intermediate-, and fast-moving waves. Wave misalignment is taken into account by applying a spanwise pressure gradient vertical to the traveling waves, which results in a three-dimensional turbulent flow field above the moving waves. Key flow quantities such as the mean velocity, velocity variances, and momentum fluxes are found to vary with the wave parameters, confirming the findings of previous studies.

direct numerical simulation↗

DOE OSTI2024

Interpreting the spatial distribution of soil properties with a physically-based distributed hydrological model

Digital soil maps are commonly data-driven as the development of physically-based models for soil mapping is difficult due to the complexity of soils. However, physically-based hydrologic models have been successful in simulating water dynamics. Since water movement is a major driver of pedogenesis, the physical rules that govern water movement might help explain and predict the spatial variation of soil properties. Here, we demonstrate the novel use of a physically-based, distributed hydrologic model to inform the spatial distribution of soil properties. The Distributed Hydrology Soil Vegetation Model (DHSVM) was utilized to simulate soil moisture content (SM) and water table depth (WTD) in two hillslope catchments under pasture and forest management wherein hydrologic model outputs were then compared with soil properties measured in situ. SM sensors and wells were installed in both catchments to validate simulations of soil water movement via Nash-Sutcliffe Efficiency (E). In-situ observations were made at 87 sites within both catchments to study the connection between simulated water movement (SM and WTD) and observed soil properties, namely the depth and thickness of the argillic (Bt), fragic (Btx), and C horizons, and the depth of redoximorphic features. The simulated time series of SM and WTD were also clustered per season using Dynamic Time Warping (DTW), which identified similarity among time series at varying timescales. Model validation suggested that simulations of surficial SM (0–20 cm) were reasonable (E = 0.45), however, simulated subsurface SM (45–60 cm) and WTD were not sufficiently accurate. The thickness of Btx horizons were spatially grouped into different populations by SM clusters from every season except spring. For the other properties, only SM dynamics of specific seasons grouped into significantly different populations, suggesting that the explanatory power of simulated water movement varies seasonally and was greater during winter. Here, we show clusters of simulated SM separated soil properties into statistically different populations, showing that hydrologic models could inform areas that followed different water dynamics related to pedogenic trajectories and related biogeochemical processes not necessarily simulated by the model. As such, physically-based modeling of water dynamics can, therefore, inform and advance digital soil mapping by linking water movement patterns stemming from hydrologic model outputs to spatial patterns of soil properties and pedogenesis.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Transition to turbulence in viscoelastic channel flow of dilute polymer solutions

The transition to turbulence in a plane Poiseuille flow of dilute polymer solutions is studied by direct numerical simulations of a finitely extensible nonlinear elastic fluid with the Peterlin closure. The range of Reynolds number ($Re$)$2000 \le Re \le 5000$is studied but with the same level of elasticity in viscoelastic flows. The evolution of a finite-amplitude perturbation and its effects on the transition dynamics are investigated. A viscoelastic flow begins transition at an earlier time than its Newtonian counterparts, but the transition time appears to be insensitive to polymer concentration in the dilute and semi-dilute regimes studied. Increasing polymer concentration, however, decreases the maximum attainable energy growth during the transition process. The critical or minimum perturbation amplitude required to trigger transition is computed. Interestingly, both Newtonian and viscoelastic flows follow almost the same power-law scaling of$Re^\gamma$with the critical exponent$\gamma \approx -1.25$, which is in close agreement with previous studies. However, a shift downward is observed for viscoelastic flow, suggesting that smaller perturbation amplitudes are required for the transition. A mechanism of the early transition is investigated by the evolution of wall-normal and spanwise velocity fluctuations and flow structure. The early growth of these fluctuations and the formation of quasi-streamwise vortices around low-speed streaks are promoted by polymers, hence causing an early transition. These vortical structures are found to support the critical exponent$\gamma \approx -1.25$. Once the transition process is completed, polymers play a role in dampening the wall-normal and spanwise velocity fluctuations and vortices to attain a drag-reduced state in viscoelastic turbulent flows.

Mechanics↗

DOE OSTI2020

Investigation of the Effect of Temperature on the Structure of SARS-CoV-2 Spike Protein by Molecular Dynamics Simulations

Statistical and epidemiological data imply temperature sensitivity of the SARS-CoV2 coronavirus. However, the molecular level understanding of the virus structure at different temperature is still not clear. Spike protein is the outermost structural protein of the SARS-CoV-2 virus which interacts with the Angiotensin Converting Enzyme 2 (ACE2), a human receptor, and enters the respiratory system. In this study, we performed an all atom molecular dynamics simulation to study the effect of temperature on the structure of the Spike protein. After 200 ns of simulation at different temperatures, we came across some interesting phenomena exhibited by the protein. We found that the solvent exposed domain of Spike protein, namely S1, is more mobile than the transmembrane domain, S2. Structural studies implied the presence of several charged residues on the surface of N-terminal Domain of S1 which are optimally oriented at 10– 30°C. Bioinformatics analyses indicated that it is capable of binding to other human receptors and should not be disregarded. Additionally, we found that receptor binding motif (RBM), present on the receptor binding domain (RBD) of S1, begins to close around temperature of 40°C and attains a completely closed conformation at 50°C. We also found that the presence of glycan moieties did not influence the observed protein dynamics. Nevertheless, the closed conformation disables its ability to bind to ACE2, due to the burying of its receptor binding residues. Our results clearly show that there are active and inactive states of the protein at different temperatures. This would not only prove beneficial for understanding the fundamental nature of the virus, but would be also useful in the development of vaccines and therapeutics.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Water Structure and Dynamics Near the Surfaces of Silicalite-1

Zeolites are crystalline microporous aluminosilicates that are commonly used as industrial sorbents and membranes. Solvent structuring and dynamics near zeolite crystal surfaces are thought to influence interfacial transport and molecular recognition processes, as well as fundamental aspects of their crystallization from solution. Here, in this work, we use molecular dynamics (MD) simulations to investigate the behavior of interfacial water near the exposed (010), (100), and (101) crystal faces of silicalite-1, one of the most widely studied zeolites. The MD simulations reveal that water’s translational and orientational order is strongly influenced by the distinct corrugations, pore apertures, and functional group distributions presented on each surface. Specifically, we observe two distinct hydration layers near each surface. Water molecules in the contact layer are hydrogen bonded to two or three exposed surface silanol groups. The relative populations of water molecules that are doubly and triply hydrogen bonded to the surface strongly depend on the distribution of the exposed silanols. The interactions with exposed surface silanols also influence water’s local orientational order and distribution across each surface. We also show that the structuring of the solvent near the different faces of silicalite-1 strongly impacts the interfacial dynamics. The translational and orientational relaxation dynamics of water are slowest in the contact layers and correlate with structural ordering near each face.

42 ENGINEERING↗

DOE OSTI2019

Computational characterization of the secondary droplets formed during the impingement of a train of ethanol drops

This paper uniquely characterizes the secondary droplets formed during the impingement of a train of ethanol drops, using three-dimensional direct numerical simulations performed under conditions studied experimentally by Yarin and Weiss. Our numerical results have been previously validated against experimental data demonstrating the ability to accurately capture the splashing dynamics. In this work, the predictive ability of the model is leveraged to gain further insight into secondary droplet formation. We present a robust post-processing algorithm, which scrutinizes the liquid volume fraction field in the volume-of-fluid method and quantifies the number, volume and velocity of secondary droplets. The high-resolution computational simulations enable secondary droplet characterization within close proximity of the impingement point at small length and time scales, which is extremely challenging to achieve experimentally. By studying the temporal evolution of secondary droplet formation, direct connections are made between liquid structures seen in the simulation and the instantaneous distribution of secondary droplets, leading to detailed insight into the instability-driven breakup process of lamellae. Time-averaged secondary droplet characteristics are also studied to describe the global distribution of secondary droplets. Such analysis is vital to understanding fuel drop impingement in direct injection engines, facilitating the development of highly accurate spray–wall interaction models for use in Lagrangian solvers.

42 ENGINEERING↗

DOE OSTI2025

Sensitivity of Regional WRF‐Chem Air Quality and Weather Simulations to Biomass‐Burning Emission Data Sets: A Case Study of the Impact of Canadian Wildfire on the US°

This study focuses on the period from June 26 to 29, 2023, when record‐breaking Canadian wildfires severely impacted air quality in the Midwest United States. Using the Weather Research and Forecasting Model with Chemistry (WRF‐Chem) and four biomass‐burning data sets (Fire Inventory from NCAR version 1, Fire Inventory from NCAR version 2.5, Quick Fire Emissions Data set [QFED], and Regional ABI‐VIIRS Emission), we analyzed aerosol transport from Canada to the US and assessed the model's accuracy in predicting PM 2.5 , O 3 , CO and aerosol weather feedback. Model simulations were compared with ground‐based and remote sensing observations as well as field measurements from the Community Research on Climate and Urban Science (CROCUS) project. Our findings show that the movement of a low‐pressure system from the Great Lakes to the Atlantic, combined with the high‐pressure system over the Atlantic, caused the transport of aerosols from Canadian wildfires to the US. Results show WRF‐Chem significantly underestimated key atmospheric components: aerosol optical depth (AOD) by over 50%, PM 2.5 by 65%–90% and peak O 3 concentrations by 50%–55% across four biomass burning data sets. Additionally, CO and NO 2 concentrations were underpredicted. The substantial underestimation of PM 2.5 led to an overestimation of temperature by up to 3.6 °C primarily due to excessive downward shortwave radiation, which resulted from the underestimation of direct aerosol effects and an increase in sensible heat flux. Among the biomass‐burning data sets, QFED produced the most accurate AOD and PM 2.5 predictions due to improved wildfire emission estimates, leading to a 1.0 to 1.5 °C reduction in temperature overestimation during the daytime. These findings underscore the need for improving wildfire emission estimates for trace gases and aerosols to enhance air quality and weather feedback predictions.

WRF-chem model↗

DOE OSTI2023

The Relationship between Age, Metallicity, and Abundances for Disk Stars in a Simulated Milky Way

Observations of the Milky Way's low-α disk show that several element abundances correlate with age at fixed metallicity, with unique slopes and small scatters around the age–[X/Fe] relations. In this study, we turn to simulations to explore the age–[X/Fe] relations for the elements C, N, O, Mg, Si, S, and Ca that are traced in a FIRE-2 cosmological zoom-in simulation of a Milky Way–like galaxy, m12i, and understand what physical conditions give rise to the observed age–[X/Fe] trends. We first explore the distributions of mono-age populations in their birth and current locations, [Fe/H], and [X/Fe], and find evidence for inside-out radial growth for stars with ages <7 Gyr. We then examine the age–[X/Fe] relations across m12i's disk and find that the direction of the trends agrees with observations, apart from C, O, and Ca, with remarkably small intrinsic scatters, σ int (0.01 - 0.04 dex). This σ int measured in the simulations is also metallicity dependent, with σ int ≈ 0.025 dex at [Fe/H] = -0.25 dex versus σ int ≈ 0.015 dex at [Fe/H] = 0 dex, and a similar metallicity dependence is seen in the GALAH survey for the elements in common. Additionally, we find that σ int is higher in the inner galaxy, where stars are older and formed in less chemically homogeneous environments. The age–[X/Fe] relations and the small scatter around them indicate that simulations capture similar chemical enrichment variance as observed in the Milky Way, arising from stars sharing similar element abundances at a given birth place and time.

79 ASTRONOMY AND ASTROPHYSICS↗