Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation model”

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 415 records · Page 23

Interpolation Methods to Enable Fast and Accurate EMT Simulation of PV Inverters

In this paper, interpolation methods to enable fast simulation of photovoltaic (PV) inverters are proposed. The interpolation methods are based on the event-driven interpolating method that can be introduced in the simulation model to run the model at a higher timestep while maintaining accuracy. The proposed methods are tested on a neutral point clamp (NPC) converter and a two-level converter (2LC) that are typically used in large-scale PV plants. The interpolation methods are integrated into the NPC and 2LC simulation models in PSCAD/Fortran and RSCAD C platforms. These models are tested at different simulation timesteps in off-line and real-time simulation platforms PSCAD and RSCAD respectively.

Marthi, Phani Ratna Vanamali↗

Tropical Cyclone Precipitation Response to Surface Warming in Aquaplanet Simulations With Uniform Thermal Forcing

While many modeling studies have attempted to estimate how tropical cyclone (TC) precipitation is impacted by climate change, the multitude of analysis techniques and methodologies have resulted in varying conclusions. Simplified models may be able to help overcome this problem. Radiative-convective equilibrium (RCE) model simulations have been used in various configurations to study fundamental aspects of Earth's climate. While many RCE modeling studies have focused on TC genesis, intensification, and size, limited work has been done using RCE to study TC precipitation. Here, in this study, the response of TC precipitation to sea surface temperature (SST) change is analyzed in global Community Atmosphere Model (CAM) aquaplanet simulations run with Radiative-Convective Equilibrium Model Intercomparison Project protocols, with the addition of planetary rotation. We expect that the insight gained about how TC precipitation responds to SST warming will help predict how TCs in the real world respond to climate change. In the CAM RCE simulations, the warmer SST simulations have less TCs on average, but the TCs tend to be larger in outer size and more intense. As simulation SST increases, more extreme precipitation rates occur within TCs, and more of the TC precipitation comes from these extreme rates. For extreme (99th percentile) TC precipitation, SST, and TC intensity increases dominate the 8.6% per K increase, while TC outer size changes have little impact. For accumulated TC precipitation, SST, and TC intensity contributions are still the majority, but TC outer size changes also contribute to the 6.6% per K increase.

54 ENVIRONMENTAL SCIENCES↗

CERESMIP: a climate modeling protocol to investigate recent trends in the Earth's Energy Imbalance

The Clouds and the Earth's Radiant Energy System (CERES) project has now produced over two decades of observed data on the Earth's Energy Imbalance (EEI) and has revealed substantive trends in both the reflected shortwave and outgoing longwave top-of-atmosphere radiation components. Available climate model simulations suggest that these trends are incompatible with purely internal variability, but that the full magnitude and breakdown of the trends are outside of the model ranges. Unfortunately, the Coupled Model Intercomparison Project (Phase 6) (CMIP6) protocol only uses observed forcings to 2014 (and Shared Socioeconomic Pathways (SSP) projections thereafter), and furthermore, many of the ‘observed' drivers have been updated substantially since the CMIP6 inputs were defined. Most notably, the sea surface temperature (SST) estimates have been revised and now show up to 50% greater trends since 1979, particularly in the southern hemisphere. Additionally, estimates of short-lived aerosol and gas-phase emissions have been substantially updated. These revisions will likely have material impacts on the model-simulated EEI. We therefore propose a new, relatively low-cost, model intercomparison, CERESMIP, that would target the CERES period (2000-present), with updated forcings to at least the end of 2021. The focus will be on atmosphere-only simulations, using updated SST, forcings and emissions from 1990 to 2021. The key metrics of interest will be the EEI and atmospheric feedbacks, and so the analysis will benefit from output from satellite cloud observation simulators. The Tier 1 request would consist only of an ensemble of AMIP-style simulations, while the Tier 2 request would encompass uncertainties in the applied forcing, atmospheric composition, single and all-but-one forcing responses. We present some preliminary results and invite participation from a wide group of models.

54 ENVIRONMENTAL SCIENCES↗

SECURED: Simulator-Enhanced Control and Understanding of Reactor systems for cyber-Event Defense

The study discusses a learning approach for analyzing cyber-events in reactor systems using integrated hardware and personal computer simulator models. Key points include the rise in cyber-attacks and their sophistication in industrial control systems (ICS), the necessity for awareness, understanding, resource allocation, and preparation to combat these threats, and the digital transformation of old and new nuclear plants, increasing their exposure to cyber threats. It highlights the cyber vulnerabilities of advanced reactor systems, which rely on digital instrumentation and control for operations and safety functions, making them susceptible to cyber-attacks. The approach involves demonstrating reactor system plant ICS cyber-attacks under various operational conditions utilizing tools like simulator models and hardware-based kits. A strategic solution approach tailored to critical infrastructure is emphasized, along with community engagement for public and government support, adopting effective learning approaches, and the preparation for anticipated future challenges. The presentation concludes with a call to action to address challenges, leverage opportunities, and advance through lesson learning in cybersecurity for nuclear energy systems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Optimizing Patient-Specific Medication Regimen Policies Using Wearable Sensors in Parkinson’s Disease

Effective treatment of Parkinson’s disease (PD) is a continual challenge for healthcare providers, and providers can benefit from leveraging emerging technologies to supplement traditional clinic care. We develop a data-driven reinforcement learning (RL) framework to optimize PD medication regimens through wearable sensors. We leverage a data set of n = 26 PD patients who wore wrist-mounted movement trackers for two separate six-day periods. Using these data, we first build and validate a simulation model of how individual patients’ movement symptoms respond to medication administration. We then pair this simulation model with an on-policy RL algorithm that recommends optimal medication types, timing, and dosages during the day while incorporating human-in-the-loop considerations on medication administration. The results show that the RL-prescribed medication regimens outperform physicians’ medication regimens, despite physicians having access to the same data as the RL agent. To validate our results, we assess our wearable-based RL medication regimens using n = 399 PD patients from the Parkinson’s Progression Markers Initiative data set. We show that the wearable-based RL medication regimens would lead to significant symptom improvement for these patients, even more so than training RL policies directly from this data set. In doing so, we show that RL models from even small data sets of wearable data can offer novel, generalizable clinical insights and medication strategies, which may outperform those derived from larger data sets without wearable data.

60 APPLIED LIFE SCIENCES↗

Energy Saving Quantification on Ductless Heat Pump (DHP) in Existing Homes

In residential retrofit applications, ductless mini-split heat pumps (DHP) are often reported to have high-energy savings potential, depending on the system they are supplementing or replacing. However, recently, there have been a number of utility studies and analyses indicating these energy savings are not being achieved when the existing system is left in place for backup heating or air conditioning. The Pacific Northwest National Laboratory conducted a three-phase project to help determine which control strategies would have the most energy savings impact in various climate zones around the United States. The first phase developed the standalone simulation model for the PNNL Lab Homes and investigated the energy-saving potential from different control strategies and HVAC system configurations. The second phase focused on conducting experiments in the PNNL Lab Homes, which tested the most promising solutions that were modeled in the first phase. The third phase used the field data to calibrate the simulation model and then extrapolated the results to different climate locations and different building sizes. This report focuses on the third phase of the study, including five major parts: model calibration, parametric model setup, results, a sensitivity analysis of air leakage rate, and conclusions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Reducing warpage in a hybrid large-scale additive manufacturing and compression molding process

In recent years, a hybrid manufacturing process, developed by combining extrusion-based large-scale additive manufacturing (AM) and compression molding (CM) techniques, has shown promising outcomes for producing structurally functional parts. The process can be used with both short fiber-reinforced composites and neat polymers and hence, even multi-material parts can be manufactured easily. This process offers the advantages of structural enhancement by having a desired fiber orientation using a large-scale AM process, as well as rapid manufacturing capability using a CM process. In the large-scale AM process, the alignment of fibers in the deposition direction enables significant improvement in the mechanical properties of the manufactured parts. However, the anisotropy resulting from the directional arrangement of fibers also introduces challenges related to warpage in the produced parts. This study aims to identify the causes of warpage and propose strategies to mitigate it. The research involves the use of preforms manufactured through the large-scale AM, which are then combined with the neat resin for CM manufacturing. A finite element-based numerical simulation model is developed, employing a sequentially coupled thermomechanical approach. Through a parametric study using the simulation models, optimization of printing direction and preform geometry is performed to minimize warpage. This contributes to the advancement and wider adoption of AM/CM hybrid manufacturing to produce structurally functional parts.

Jo, Eonyeon↗

Avista CEF2 Shared Energy Economy: Modeling and Simulation

This research project supported Avista’s Shared Energy Economy (SEE) Project in Spokane’s University District (UDistrict), as funded by the WA Department of Commerce Clean Energy Fund (CEF) Phase 2. The SEE project comprises a microgrid with two campus buildings, two solar photovoltaic generators, and two batteries. The microgrid can operate grid-connected or off-grid. This research project included modeling the buildings and other microgrid components, simulating a transactive energy system, and supporting Avista’s testing of different use cases for the WA Department of Commerce. The main research finding is that buildings can be modeled successfully in a data-driven approach.

14 SOLAR ENERGY↗

Fluid Dynamics of Axial Turbomachinery: Blade- and Stage-Level Simulations and Models

The current generation of axial turbomachines is the culmination of decades of experience, and detailed understanding of the underlying flow physics has been a key factor for achieving high efficiency and reliability. Driven by advances in numerical methods and relentless growth in computing power, computational fluid dynamics has increasingly provided insights into the rich fluid dynamics involved and how it relates to loss generation. This article presents some of the complex flow phenomena occurring in bladed components of gas turbines and illustrates how simulations have contributed to their understanding and the challenges they pose for modeling. The interaction of key aerodynamic features with deterministic unsteadiness, caused by multiple blade rows, and stochastic unsteadiness, i.e., turbulence, is discussed. High-fidelity simulations of increasingly realistic configurations and models improved with help of machine learning promise to further grow turbomachinery performance and reliability and, thus, help fluid mechanics research have a greater industrial impact.

Mechanics↗

Train small, model big: Scalable physics simulators via reduced order modeling and domain decomposition

Numerous cutting-edge scientific technologies originate at the laboratory scale, but transitioning them to practical industry applications is a formidable challenge. Traditional pilot projects at intermediate scales are costly and time-consuming. An alternative, the pilot-scale model, relies on high-fidelity numerical simulations, but even these simulations can be computationally prohibitive at larger scales. To overcome these limitations, we propose a scalable, physics-constrained reduced order model (ROM) method. The ROM identifies critical physics modes from small-scale unit components, projecting governing equations onto these modes to create a reduced model that retains essential physics details. We also employ Discontinuous Galerkin Domain Decomposition (DG-DD) to apply ROM to unit components and interfaces, enabling the construction of large-scale global systems without data at such large scales. Here this method is demonstrated on the Poisson and Stokes flow equations, showing that it can solve equations about 15–40 times faster with only ~1% relative error. Furthermore, ROM takes one order of magnitude less memory than the full order model, enabling larger scale predictions at a given memory limitation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hardware In the Loop for Demand Flexibility (HIL4DF) v1.0

The software package in question is a collection of simulation models in the Modelica language, representing a variety of mechanical system designs and envelope conditions related to LBL's FLEXLAB facility. The collection of models also features multiple controls sequences that can be simulated with the FLEXLAB model to simulate different demand flexibility scenarios. Additionally, this package will feature datasets from 3 experimental tests, used for calibration, validation and comparison against the Modelica models, this includes weather data that can be used to replicate different scenarios in simulation across the same weather conditions experienced in real experiments. Given FLEXLAB high level of instrumentation and available data, the models are calibrated across multiple measurement points, and thus results from the extension of this model to other climate zones or control sequences, would provide high level of confidence.

Huang, Weiping↗

Precipitation Characteristics and Future Changes Over the Southern Slope of Tibetan Plateau Simulated by a High-Resolution Global Nonhydrostatic Model

Current climate model resolution cannot accurately describe the complex topography over the Tibetan Plateau, which limits our understanding of past and future precipitation over this region. This study investigates the daily precipitation characteristics and its future changes over the Tibetan Plateau, especially for the southern slope of Tibetan Plateau (SSTP), by 14-km Nonhydrostatic ICosahedral Atmospheric Model (NICAM) with explicitly calculated convection for historical and future 30-years period. By comparing with the satellite Global Precipitation Measurement (GPM), NICAM well reproduces the historical precipitation spatial pattern, seasonal cycle, and the extreme precipitation belt over SSTP, but overestimates precipitation amount by~35%. It is found that heavy precipitation probability decreases as elevation becomes higher while the light precipitation generally shows the opposite. For the precipitation changes during June to September from 1979–2008 to 2075–2104, NICAM predicts that mean precipitation will decrease over low-level SSTP but increase over high-level SSTP. Nonprecipitation and heavy precipitation probability will increase while light precipitation probability will decrease in the future over SSTP. The extreme precipitation probability and intensity will increase ~50%/°C and ~8%/°C over SSTP, and this increase is more obvious as elevation becomes higher. The robust increase of extreme precipitation along the SSTP topography is unique and has not been identified by the climate model simulations before. The strong meridional gradient of specific humidity over SSTP is found to be further enhanced under global warming, and this gradient enhancement is suggested to be responsible for the increase in the extreme precipitation over SSTP.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deciphering the R Ψ m correlator in search for the chiral magnetic effect in relativistic heavy ion collisions

The chiral magnetic effect (CME) is extensively studied in heavy-ion collisions at the BNL Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). An azimuthal correlator called R Ψ m was proposed to measure the CME. By observing the same R Ψ 2 and R Ψ 3 (convex) distributions from A Multi-Phase Transport (ampt) model, by contrasting data and model as well as large and small systems and by event shape engineering (ESE), a recent preprint (arXiv:2006.04251v1) from STAR suggests that the R Ψ m observable is sensitive to the CME signal and relatively insensitive to backgrounds, and their Au + Au data are inconsistent with known background contributions. Our ampt results contradict those from a recent preprint by STAR. Our toy model simulations demonstrate that R Ψ 2 is sensitive to both the CME signal and physics backgrounds. Toy model simulations of isobar collisions show similar centrality dependence and magnitudes for the relative R Ψ 2 strengths as well as the relative Δ γ strengths. We conclude that R Ψ 2 and the inclusive Δ γ are essentially the same.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A reduced-model (nSOLT) simulation of neutral recycling effects on plasma turbulence in the divertor region of MAST-U

The 2D scrape-off-layer turbulence code (nSOLT), which includes 1D kinetic neutral–plasma interactions, is applied to study effects of neutral recycling on plasma turbulence for parameters illustrative of the MAST-U divertor region. Neutral recycling is modeled by injecting a fraction of the parallel plasma flux to the divertor back into the simulation domain as a source of Franck–Condon-distributed neutrals. Stationary sources, concentrated at the magnetic separatrix, model plasma streaming into the divertor region from the upstream scrape-off-layer and sustain plasma turbulence absent neutral recycling. Starting from one such no-neutrals equilibrium, we initiate recycling in a numerical experiment designed to diagnose and identify the effects of various neutral–plasma interactions on the divertor plasma, divertor turbulence, and plasma exhaust. The onset of recycling triggers an initial burst of enhanced cross field plasma transport that is quelled by ionization cooling and charge–exchange (CX) friction, with growing neutral pressure, leading to a quiescent, turbulence-free state. Diagnosis of this transient burst reveals that (1) the sudden increase in plasma density due to ionization dominates the onset of the burst, (2) electron cooling due to ionization increases collisionality and disconnects blob filaments from the sheath, and (3) CX friction drives tripole polarization of a blob that can dominate the curvature-driven dipole polarization, leading to the stagnation of blob propagation and reduced radial turbulent transport. It is shown that CX friction is negligible compared to sheath physics in determining equilibrium mean flow shearing rates, for parameters considered herein (specifically a short connection length to the divertor target), while it can significantly reduce interchange-instability growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DEPLOYING FAST CHARGING INFRASTRUCTURE FOR ELECTRIC VEHICLES IN URBAN NETWORKS: AN ACTIVITY-BASED APPROACH

This paper explores an important problem under the domain of network modeling, the optimal configuration of charging infrastructure for electric vehicles (EVs) in urban networks considering EV users' daily activities and charging behavior. This study proposes a charging behavior simulation model considering different initial state of charge (SOC), travel distance, availability of home chargers, and the daily schedule of trips for each traveler. The proposed charging behavior simulation model examines the complete chain of trips for EV users as well as the interdependency of trips traveled by each driver. The problem of finding the optimum charging configuration is then formulated as a mixed-integer nonlinear programming problem that considers the dynamics of travel time and travel distance, the interdependency of trips made by each driver, limited range of EVs, remaining battery capacity for recharging, waiting time in queue, and detour to access a charging station. This problem is solved using a metaheuristic approach for a large-scale case network. A series of examples are presented to demonstrate the model efficacy and explore the impact of energy consumption on the final SOC and the optimum charging infrastructure.

Chain of Trips↗

Predicting Deformation of Parts from HP Printers due to Sintering with Modeling and Simulation (CRADA NFE-18-07363)

This work proposes to leverage Oak Ridge National Laboratory’s (ORNL) sintering simulation software developed at the Manufacturing Demonstration Facility to 1) investigate the deformation of sintered parts and 2) to develop fast surrogate-based optimization methods for HP Inc. This includes the development of ORNL’s software on HP systems, integration and testing with HP part geometries and material properties, the development of surrogate models based on the results of the simulation, and modifications to the software based on feedback from HP. The expected impacts of this work are a drastic reduction in print costs, an increase in printer throughput, and the identification of novel, possibly exotic materials for 3D printing. This work will also decrease the time to market for HP’s print technology by decreasing the amount of resources and type dedicated to physical testing.

42 ENGINEERING↗

Modeling and simulation of earthquake soil structure interaction excited by inclined seismic waves

Presented is an application of wave potential formulation (WPF) together with domain reduction method (DRM) to modeling earthquake soil structure interaction (ESSI) behavior in horizontally layered ground under inclined incident seismic waves. Wave potential formulation is used to develop a spatially varying, inclined seismic wave field from incident Primary (P) and Secondary (S) waves that propagate through layered ground. Developed seismic wave field is then used to develop effective forces for Domain Reduction Method that are then used for analyzing ESSI response of a soil structure system. Developed methodology, called WPF-DRM, is verified using analytic solution for a free field response of layered ground subjected to inclined incident waves. Developed WPF-DRM methodology is illustrated through analysis of an ESSI response of a deeply embedded structure, a small modular reactor (SMR) subjected to incident S wave polarized in vertical plane (SV) with variation in inclinations and frequencies. Presented example highlights the influences of incident wave inclination and frequency on ESSI response of analyzed SMR.

58 GEOSCIENCES↗

First Estimation of Model Parameters for Neutrino-Induced Nucleon Knockout Using Simulation-Based Inference

To enable an accurate determination of oscillation parameters, accelerator-based neutrino experiments require detailed simulations of nuclear interaction physics in the GeV regime. While substantial effort from both theory and experiment is currently being invested to improve the fidelity of these simulations, their present deficiencies typically oblige experimental collaborations to resort to empirical tuning of simulation model parameters. As the precision requirements of the field continue to become more stringent, machine learning techniques may provide a powerful means of handling corresponding growth in the complexity of future neutrino interaction model tuning exercises. To study the suitability of simulation-based inference (SBI) for this physics application, in this paper we revisit a tuned configuration of the GENIE neutrino event generator that was originally developed by the MicroBooNE collaboration. Despite closely reproducing the adopted values of four physics parameters when confronted with the tuned cross-section predictions as input, we find that our trained SBI algorithm prefers modestly different values (within MicroBooNE's assigned uncertainties) and achieves slightly better goodness-of-fit when inference is run on the experimental data set originally used by MicroBooNE. We also find that our trained algorithm can create a fair approximation of an alternative neutrino scattering simulation, NuWro, that shares only a subset of its physics model parameters with GENIE.

Tame-Narvaez, Karla [Fermilab] (ORCID:000000022249↗