Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Behavior modeling”

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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Geometric pumping and dephasing at topological phase transition

A measure-preserving formalism (MPF) is constructed and applied to spin/band models, which yield observations about pumping. It occurs at the topological phase transition, i.e., a Z 2 flip, suggesting that Z 2 can imply bulk effects. The model's asymptotic behavior is analytically solved via MPF. The pumping probability is geometric, fractional, and has a ceiling of $\frac{1}{2}$. Intriguingly, theorems are proved about occurrence conditions, which are linked to the system's dimension and the distinction between rational and irrational numbers. Finally, experimental detection is discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Steepest-descent algorithm for simulating plasma-wave caustics via metaplectic geometrical optics

The design and optimization of radiofrequency-wave systems for fusion applications is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at wave cutoffs and caustics. To accurately model the wave behavior in these regions, more advanced and computationally expensive “full-wave” simulations are typically used, but this is not strictly necessary. A new generalized formulation called metaplectic geometrical optics (MGO) has been proposed that reinstates GO near caustics. The MGO framework yields an integral representation of the wavefield that must be evaluated numerically in general. We present an algorithm for computing these integrals using Gauss-Freud quadrature along the steepest-descent contours. Benchmarking is performed on the standard Airy problem, for which the exact solution is known analytically. Furthermore, the numerical MGO solution provided by the new algorithm agrees remarkably well with the exact solution and significantly improves on previously derived analytical approximations of the MGO integral.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Prediction of defect properties in concentrated solid solutions using a Langmuir-like model

The alleged existence of sluggish diffusion in high-entropy alloys has drawn controversy. In high-entropy alloys and, in general, in all solids, transport properties are controlled by point defect concentration, which must be known before performing atomistic simulations to compute transport coefficients. In this work, we present a general Langmuir-like model for defect concentration in an arbitrarily complex solid solution and apply this model to generate expressions for concentrations of vacancies and small interstitial atoms. We then calculate the vacancy concentration as a function of temperature in the equiatomic CoNiCrFeMn and FeAl alloys with modified embedded-atom-method potentials for various chemical orderings, showing there is no clear correlation between vacancy thermodynamics and chemical ordering in the CoNiCrFeMn alloy, but clear systematic patterns for FeAl. We believe this is due to the high stability of disordered, random, and ordered intermetallic phases, respectively, in the CoNiCrFeMn and FeAl systems. Finally, this work provides future avenues to the prediction of thermal interstitials and vacancies in solid solutions, which is necessary for models of nonequilibrium behavior of solid solutions.

composition↗
DOE OSTI2025

When more data hurts: Optimizing data coverage while mitigating diversity-induced underfitting in an ultrafast machine-learned potential

Machine-learned interatomic potentials (MLIPs) are becoming an essential tool in materials modeling. However, optimizing the generation of training data used to parametrize the MLIPs remains a significant challenge. This is because MLIPs can fail when encountering local environments too different from those present in the training data. The difficulty of determining a priori the environments that will be encountered during molecular dynamics simulation necessitates diverse, high-quality training data. Here, this study investigates how training data diversity affects the performance of MLIPs using the Ultra-Fast force field (UF 3 ) to model amorphous silicon nitride. We employ expert and autonomously generated data to create the training data and fit four force field variants to subsets of the data. Our findings reveal a critical balance in training data diversity: insufficient diversity hinders generalization, while excessive diversity can exceed the MLIP's learning capacity, reducing simulation accuracy. Specifically, we found that the UF 3 variant trained on a subset of the training data, in which nitrogen-rich structures were removed, offered vastly better prediction and simulation accuracy than any other variant. By comparing these UF 3 variants, we highlight the nuanced requirements for creating accurate MLIPs, emphasizing the importance of application-specific training data to achieve optimal performance in modeling complex material behaviors.

ab initio molecular dynamics↗
DOE OSTI2024

Hybrid Modeling of Three-Phase Grid-Supporting Inverters for Dynamic Studies

Grid technologies connected by power electronic converter (PEC) interfaces continually implement grid support functions mandated by grid codes and standards. The transition to converter-based generation demands precise PEC models to assess system dynamics, which have been previously overlooked in conventional power systems. This study proposes a hybrid method for analyzing grid-connected three-phase PEC dynamics with the IEEE standard 1547-2018 Volt-VAr mode that combines physics and data-driven techniques. The physics model reflects the PEC’s internal behavior, whereas the data-driven modeling technique evaluates the grid-supporting capabilities of the smart PEC. The system identification approach is used to generate dynamic PEC models based on changing grid voltage and measured current injected into the grid by the PEC. In the Volt-VAr support mode, a detailed topological model including switches is utilized to compare the goodness-of-fit of the extracted hybrid dynamic model. The results demonstrate that the hybrid PEC model in the Volt-VAr mode accurately matches the dynamics with the topological model.

Subedi, Sunil↗
DOE OSTI2024

Pretest Modeling A Spent Nuclear Fuel Seismic Shake Test

The U.S. Department of Energy Spent Fuel and Waste Science and Technology (SFWST) program is planning to conduct a series of full-scale shake table tests to simulate hypothetical earthquake conditions and record the response of surrogate spent nuclear fuel (SNF) assemblies in a dry canister storage system mockup. The shake table motions will represent a range of hypothetical earthquake conditions at hypothetical locations in the continental U.S. to generally define the range of mechanical loads that SNF can be expected to experience during extended dry storage periods. This paper describes the pretest predictions made with LS-DYNA models of the mockup storage systems. The test will use two dry storage system configurations, a mockup vertical concrete cask (VCC) and a mockup horizontal storage module (HSM). The test will use a production-quality canister and basket. Within the canister will be four instrumented fuel assemblies with fuel rods containing surrogate mass and 28 instrumented dummy assemblies that are intended to match the mass and outer dimension of a fuel assembly. The finite element models include models of the VCC and HSM on the shake table to calculate the system level dynamic responses and separate single fuel assembly models to calculate stress and strain in fuel assembly components. Both types of models include nonlinear behavior like rod-to-rod contact and the ability for VCC’s to rock and slide. This paper presents the expected response of the VCC, HSM, and fuel assemblies to the shake table testing that is planned to start in April of 2024. The earthquake conditions represent seismic hazards in the 2,000-to-20,000-year return period range. The test data is expected to confirm the expectation that fuel rod cladding will remain intact, fuel assembly structural components like guide tubes will remain intact, and no significant VCC sliding or tipping will occur in the range of conditions to be tested with the shake table.

Klymyshyn, Nicholas A.↗
DOE OSTI2024

Multilevel Parareal Algorithm with Averaging for Oscillatory Problems

The present study is an extension of the work done by Peddle, Haut, and Wingate and Haut and Wingate, where a two-level Parareal method with mapping and averaging is examined. The method proposed in this paper is a multilevel Parareal method with arbitrarily many levels, which is not restricted to the two-level case. We give an asymptotic error estimate which reduces to the two-level estimate for the case when only two levels are considered. Introducing more than two levels has important consequences for the averaging procedure, as we choose separate averaging windows for each of the different levels, which is an additional new feature of the present study. The different averaging windows make the proposed method especially appropriate for nonlinear multiscale problems, because we can introduce a level for each intrinsic scale of the problem and adapt the averaging procedure such that we reproduce the behavior of the model on the particular scale resolved by the level. The method is applied to nonlinear differential equations. The nonlinearities can generate a range of frequencies in the problem. The computational cost of the new method is investigated and studied on several examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Simulating the Impact of Dynamic Rerouting on Metropolitan-scale Traffic Systems

The rapid introduction of mobile navigation aides that use real-time road network information to suggest alternate routes to drivers is making it more difficult for researchers and government transportation agencies to understand and predict the dynamics of congested transportation systems. Computer simulation is a key capability for these organizations to analyze hypothetical scenarios; however, the complexity of transportation systems makes it challenging for them to simulate very large geographical regions, such as multi-city metropolitan areas. In this article, we describe enhancements to the Mobiliti parallel traffic simulator to model dynamic rerouting behavior with the addition of vehicle controller actors and vehicle-to-controller reroute requests. The simulator is designed to support distributed-memory parallel execution using discrete event simulation and be scalable on high-performance computing platforms. We demonstrate the potential of the simulator by analyzing the impact of varying the population penetration rate of dynamic rerouting on the San Francisco Bay Area road network. Using high-performance parallel computing, we can simulate a day in the San Francisco Bay Area with 19 million vehicle trips with 50 percent dynamic rerouting penetration over a road network with 0.5 million nodes and 1 million links in less than three minutes. We present a sensitivity study on the dynamic rerouting parameters, discuss the simulator’s parallel scalability, and analyze system-level impacts of changing the dynamic rerouting penetration. Furthermore, we examine the varying effects on different functional classes and geographical regions and present a validation of the simulation results compared to real-world data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Computational Modeling of Heterogeneity of Stress, Charge, and Cyclic Damage in Composite Electrodes of Li-Ion Batteries

Charge heterogeneity is a prevalent feature in many electrochemical systems. In a commercial cathode of Li-ion batteries, the composite is hierarchically structured across multiple length scales including the sub-micron single-crystal primary-particle domains up to the macroscopic particle ensembles. The redox kinetics of charge transfer and mass transport strongly couples with mechanical stresses. This interplay catalyzes substantial heterogeneity in the charge (re)distribution, stresses, and mechanical damage in the composite electrode during charging and discharging. We assess the heterogeneous electrochemistry and mechanics in a LiNi x Mn y Co z O 2 (NMC) cathode using a fully coupled electro-chemo-mechanics model at the cell level. A microstructure-resolved model is constructed based on the synchrotron X-ray tomography data. We calculate the stress field in the composite and then quantitatively evaluate the kinetics of surface charge transfer and Li transport biased by mechanical stresses. We further model the cyclic behavior of the cell. The repetitive deformation of the active particles and the weakening of the interfacial strength cause gradual increase of the interfacial debonding. The mechanical damage impedes electron transfer, incurs more charge heterogeneity, and results in the capacity degradation in batteries over cycles.

25 ENERGY STORAGE↗
DOE OSTI2024

Generic Multi-Layer Perceptron Inference Accelerator on FPGA (vneuron) v1.0

We have designed and implemented a neural network inference compute engine (vneuron) that can be deployed in the fabric of any FPGA without using special hardware accelerator primitive. The "vneuron" is purely written in verilog, and supports scalable neural network structure with fully connected layers and ReLU activation ( Multi-Layer Perceptron architecture) with 16 bits of precision. We have demonstrated it on an Xilinx Artix 7 FPGA for a 16-input, 8-output MLP with 3 layer, 1600 parameters. It takes 40 DSP48E and 40 BRAM18, and takes 131 clock cycles for computing (1048 ns when clocked at 125MHz). We include PyTorch quantization from a given floating point model, and provide behavioral verification simulation in the disclosed software package.

Du, Qiang↗
DOE OSTI2024

Framework For Performing Time-dependent Pebble Bed Reactor Simulations

The present work details the creation of a high-fidelity Monte Carlo methodology for analyzing the run-in and subsequent approach to equilibrium for PBRs. The methodology entails a Python module wrapped around Serpent so as to perform neutronics calculations, move pebbles, refuel the core, and discharge pebbles, thereby modeling the explicit behavior of the PBR run-in. The code kugelpy is within the GitHub repository `pyrates`.

Stewart, RyanH. [Idaho National Laboratory (INL), ↗
DOE OSTI2021

A population data-driven workflow for COVID-19 modeling and learning

CityCOVID is a detailed agent-based model that represents the behaviors and social interactions of 2.7 million residents of Chicago as they move between and colocate in 1.2 million distinct places, including households, schools, workplaces, and hospitals, as determined by individual hourly activity schedules and dynamic behaviors such as isolating because of symptom onset. Disease progression dynamics incorporated within each agent track transitions between possible COVID-19 disease states, based on heterogeneous agent attributes, exposure through colocation, and effects of protective behaviors of individuals on viral transmissibility. Throughout the COVID-19 epidemic, CityCOVID model outputs have been provided to city, county, and state stakeholders in response to evolving decision-making priorities, while incorporating emerging information on SARS-CoV-2 epidemiology. Here we demonstrate our efforts in integrating our high-performance epidemiological simulation model with large-scale machine learning to develop a generalizable, flexible, and performant analytical platform for planning and crisis response.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

An Automatic Learning Framework for Smart Residential Communities

Predictability has been foundational to matching supply and demand in the day-to-day operation of the electric power system. Demand predictability is eroding because of the increased use of renewable energy resources and more sophisticated loads, such as electric vehicles and smart appliances. In this paper, an automatic software framework is described which can be used for load forecasting in smart communities. A time-varying clustering-based Markov chain approach is used to predict the energy consumption of residential buildings in a smart community. The training data is based on 1-minute meter data of occupied homes over one month. The data points are first clustered based on the energy consumption and the time of the day. Then, the original data is converted using the Centroids of the clusters. A time-varying Markov chain is subsequently trained to model the energy consumption behavior of residents for each home using the transformed data. The trained model is shown to successfully predict load in 5-minute intervals over a 24 hours period.

Zandi, Helia↗
DOE OSTI2020

FOSWEC v2: System identification and control testing at Hinsdale Wave Research Laboratory

This report summarizes the design and execution of a wave tank test of the floating oscillating surge wave energy converter (FOSWEC) in the O.H. Hinsdale Wave Research Laboratory Directional Wave Basin at Oregon State University. This device, which uses two "flaps" that pivot about a central platform when excited by waves, has a natural frequency within the range of the waves by which it is excited. The FOSWEC was originally considered to be a 1:33 scale device, however, for the current tests, no fixed relative scale is used (i.e., the WEC is considered to be scaled for the basin’s wave environment in which it operates). The primary goal of this test was to assess the degree to which previously developed modeling, experimentation, and control design methods could be applied to a broad range of wave energy converter designs. Testing was conducted to identify a dynamic model for the impedance and excitation behavior of the device. Using these models, a series of closed loop tests were conducted using a causal impedance matching controller. This report provides a brief description of the results, as well as a summary of the device and experimental design. The results show that the methods applied to this experimental device perform well and should be broadly applicable. The data collected during testing is compressed into FOSWEC.zip. Please refer to Appendix C (pages 61-63) of the test report for descriptions of each test ID corresponding to the compressed files.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Polymer Retention Evaluation in a Heavy Oil Sand for a Polymer Flooding Application on Alaska’s North Slope

For a polymer flooding field trial in a heavy oil reservoir on Alaska’s North Slope, polymer retention is a key parameter. Because of the economic impact of retention, this parameter was extensively studied using field core material and conditions. Several important findings were noted. Polymer retention based on effluent viscosity measurements can be overestimated unless the correct (non-linear) relation between polymer concentration and viscosity is used. Polymer degradation (either mechanical or oxidative) can also lead viscosity-based measurements to overestimate retention. Inaccessible pore volume (IAPV) can be overestimated if insufficient brine is flushed through the sand between polymer banks. Around 100 pore volumes (PV) of brine may be needed to displace mobile polymer to approach a true residual resistance factor and properly measure IAPV. Even for a sand pack with kwsor = 20 mD, IAPV was zero for HPAM with molecular weight (Mw) of 18 million g/mol. Fine-grain particles (<20 μm) strongly impacted polymer retention values. Native NB#1 sand with a significant component of particles <20 μm exhibited 290 μg/g, while the same sand exhibited 28 μg/g after these small particles were removed. Polymer retention did not necessarily correlate with mineral composition. The NB#1, NB#3, and OA sands had similar elemental and clay compositions, but the NB#1 sand exhibited ~10 times higher retention than the NB#3 sand. Polymer retention did not necessarily correlate with permeability. NB#1 sand exhibited much higher retention than OA sand, even though NB#1 sand was twice as permeable as OA sand. No evidence of chromatographic separation of HPAM molecular weights was found in our experiments. Although retention tended to be greater without a residual oil saturation (than at Sor), the effect was not strong. Aging a core (with high oil saturation) at 60°C reduced HPAM retention by a factor of two. As expected, under similar conditions, polymer retention was greater for a higher Mw HPAM (18 million g/mol) than for a lower Mw HPAM (10-12 million g/mol). In many cases with high polymer retention values (e.g., 240 μg/g), polymer arrival at the end of the core was relatively quick, but achieving the injected concentration occurred gradually over many pore volumes. This effect was not due chromatographic separation of polymer molecular weights. Results from modeling of this behavior was consistent with concentrationdependent polymer retention. The form assumed for the retention function in a simulator can have an important impact on the timing and magnitude of the oil response from a polymer flood. Field-based observations can underestimate polymer retention, depending on when the tracer and polymer concentrations were measured and the assumptions made about reservoir heterogeneity.

Wang, Dongmei↗
DOE OSTI2020

Updated Available Drawdowns for Big Hill SPR Caverns - Model Including the Caprock Fault

This report updates the estimated values of the baseline available drawdowns for the caverns at the Big Hill storage facility, and an updated table listing the available drawdowns. An updated finite element numerical analysis model, which included a fault in the caprock layers, was constructed and the daily data of actual wellhead pressures and oil-brine interfaces was used. The number of available drawdowns for each of the Big Hill SPR caverns is estimated using the new model. All caverns are predicted to have five available drawdowns remaining from a geomechanical perspective. BH-101 and 105 have a region of concern at the floor edge and/or sloping floor, where tensile and dilatant stresses are predicted to occur during each workover. The tensile state is predicted to occur because of the geometries of the edge and floor. Therefore, geomechanical examination for two caverns would be recommended after a drawdown leach. The well integrity of each cavern is not investigated in this report. The estimates for the number of baseline available drawdowns are subject to change in the future as the knowledge of physical phenomena at the sites, and the further development of the models of geomechanical behavior at the sites, evolve over time.

02 PETROLEUM↗
DOE OSTI2020

FOSWEC dynamics and controls test report

This report describes the testing of a model scale wave energy converter. This device, which uses two aps that pivot about a central platform when excited by waves, has a natural frequency within the range of the waves by which it is excited. The primary goal of this test was to assess the degree to which previously developed modeling, experimentation, and control design methods could be applied to a broad range of wave energy converter designs. Testing was conducted to identify a dynamic model for the impedance and excitation behavior of the device. Using these models, a series of closed loop tests were conducted using a causal impedance matching controller. This report provides a brief description of the results, as well as a summary of the device and ex- perimental design. The results show that the methods applied to this experimental device perform well and should be broadly applicable.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

K-Modules

K-Modules are capacitive energy modules designed to be stacked together to form a capacitor bank (aka a Marx bank). Discharging the stack produces a single, high voltage, high current pulse that when coupled to an X-ray (or electron) tube, produces energetic X-rays (or an intense electron beam, E-beam). The energetic X-rays are used in Flash X-ray Radiography (FXR), a technique to image objects moving at extremely high speed to measure their speed, shape and internal density profiles. The images help validate computer models for materials behavior under extreme pressure and temperature. FXR is a technology critical to dynamic experiments for the national security mission, and most users are DOE and DOD labs. Febetrons, the X-ray generating devices, are used at these sites and rely on capacitive energy modules for their operation. A stack of the newly developed K-modules will significantly extend the capabilities of these Febetrons

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗