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,567 records · Page 87

Quantum Simulation of Polarized Light-Induced Electron Transfer with a Trapped-Ion Qutrit System

Electron transfer within and between molecules is crucial in chemistry, biochemistry, and energy science. Here this study describes a quantum simulation method that explores the influence of light polarization on electron transfer between two molecules. By implementing precise and coherent control among the quantum states of trapped atomic ions, we can induce quantum dynamics that mimic the electron-transfer dynamics in molecules. We use three-level systems (qutrits), rather than traditional two-level systems (qubits), to enhance the simulation efficiency and realize high-fidelity simulations of electron-transfer dynamics. We treat the quantum interference between the electron coupling pathways from a donor with two degenerate excited states to an acceptor and analyze the transfer efficiency. We also examine the potential error sources that enter the quantum simulations. The trapped-ion systems have favorable scalings with system size compared to those of classical computers, promising access to richer electron-transfer simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Designing reinforcement learning algorithms for building HVAC control: From experimental observation to simulation comparisons

Advanced supervisory-level control with reinforcement learning (RL) is regarded as a promising solution for HVAC systems to minimize energy consumption while maintaining thermal comfort and indoor air quality. However, most RL applications were conducted in the simulation environment rather than real-world HVAC systems. This paper developed a value-based RL controller termed Deep Q-Network (DQN) for a typical central HVAC system and evaluated its performance in a building test facility. By comparing DQN with a rule-based controller, the study not only demonstrated the cases where DQN could properly maintain indoor comfort but also discussed possible reasons why DQN failed in some other situations. Recognizing the limitations of value-based RL algorithms from the experimental tests, a simulation study was conducted to compare DQN with an alternative RL approach, an actor–critic algorithm termed Deep Deterministic Policy Gradient (DDPG). In scenarios with a relatively large action space, DDPG outperformed DQN by requiring fewer computational resources and achieving better thermal comfort, lower energy consumption, and more stable control actions. The findings suggest that the ability of DDPG to handle continuous control variables more effectively allows for faster convergence in training and more precise control in practice, which enhances the overall efficiency and reliability of the HVAC system.

Guo, Fangzhou↗

Universal Solvent Viscosity Reduction via Hydrogen Bonding Disruptors

Liquid Ion Solutions LLC (DBA RoCo Global) in partnership with Carnegie Mellon University and Carbon Capture Scientific LLC, has performed lab-scale development and evaluation of novel additives that lower the viscosity of water-lean amine solvents for post-combustion carbon dioxide capture. This project focuses on developing additives that minimize the formation of long-range electrostatic and hydrogen bonding (HB) networks, decreasing the solvent viscosity, improving diffusion, and improving the process economics. The project objectives included: 1) performing computer simulation to understand the molecular interactions of the additive molecules in water-lean CO 2 capture solvents, 2) design and synthesis of HB disruptors additives, 3) performance testing with additive molecules on model amine solvents, and 4) demonstration of the effectiveness of the optimized additives in the presence of synthetic flue gas. To meet the abovementioned objectives, the project team utilized a holistic approach that combines molecular simulation, experimental testing, and economic analysis studies. The project team developed ab initio molecular model and then perform computer simulation to develop relationship between hydrogen bonding, viscosity, and performed quantitative analysis of additive on the viscosity of the solvent. The team completed computational comparative study on a range of organic functional groups such as ethers, esters, cyclic carbonates, alkanes, and ammonium salts for their effect on viscosity gaining key insights into molecular interactions and the impact of various functional groups and molecular shapes on viscosity. Assisted with molecular simulation insights, the project team conducted additive synthesis and testing, including a proof-of-concept study, additive screening, optimization, and synthetic flue gas testing. The experimental proof-of-concept study proved that the hydrogen bonding acceptors result in significant decrease of viscosities. Detailed additive screening (exploring various functionalities and molecular structures) has been performed. Several promising additives showed excellent reduction in viscosity (30-41%) at 5% additive loading, and over 50% viscosity reduction at 10% additive loading for the model solvents. The team also performed complex screening studies on additive loadings and mixing effect among additives using the design of experiments. Based on multiple screening experiments, one additive-solvent candidate was down-selected for synthetic flue gas testing. A 100-hour continuous absorption/desorption study was conducted under simulated flue gas using a lab-scale continuous capture and separation system. No degradation (for both solvent and additive) was observed based on the GC results of the solvent samples collected from the continuous study. The team conducted preliminary engineering analyses and cost-benefit analyses to quantify the potential economic benefits of the additive approach for solvent viscosity reduction. Based on the experimental data, CO 2 capture cost savings from the capital and operating cost savings are estimated at $\$$4.7/tonne and $\$$0.3/tonne CO 2 captured, respectively. Considering the additive cost, the net benefit is estimated to be between $\$$4.32~$\$$4.86/tonne CO 2 captured.

20 FOSSIL-FUELED POWER PLANTS↗

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗

Potential impact of work from home jobs on residential energy bills: A case study in phoenix, AZ, USA

Nearly one-third of U.S. households face challenges paying energy bills. During the day, many residents have routine access to cooled environments provided by others—employers, shopping centers, and other public buildings. The COVID-19 pandemic, however, has significantly shifted the cost burden of air conditioning in hot cities. Specifically, during the pandemic, many companies either laid off employees or put them in work-from-home (WFH) assignments. Large tech companies are already promoting WFH as a long-term option for their employees, even after the pandemic. This change in the nature of the workforce might reduce daily travel expenses for workers, but could also significantly increase residential energy bills, particularly during summer in very hot climates. Here, this study uses building energy simulations to quantify the potential residential energy bill penalties resulting from WFH for typical residences in Phoenix. Four building archetypes are used in this study to represent variations in building vintage, occupancy, and characteristics. The results show that, for some single-family residences in Phoenix, WFH can increase annual energy bills by more than $1100 (up to a 70% increase). The study also demonstrates that building performance enhancement retrofit measures have the potential to reduce this WFH energy bill penalty for the existing buildings substantially.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

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

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗