Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic 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 469 records · Page 26

Surrogate Modeling of Nonlinear Dynamic Systems: A Comparative Study

Surrogate models play a vital role in overcoming the computational challenge in designing and analyzing nonlinear dynamic systems, especially in the presence of uncertainty. This paper presents a comparative study of different surrogate modeling techniques for nonlinear dynamic systems. Four surrogate modeling methods, namely, Gaussian process (GP) regression, a long short-term memory (LSTM) network, a convolutional neural network (CNN) with LSTM (CNN-LSTM), and a CNN with bidirectional LSTM (CNN-BLSTM), are studied and compared. All these model types can predict the future behavior of dynamic systems over long periods based on training data from relatively short periods. The multi-dimensional inputs of surrogate models are organized in a nonlinear autoregressive exogenous model (NARX) scheme to enable recursive prediction over long periods, where current predictions replace inputs from the previous time window. Three numerical examples, including one mathematical example and two nonlinear engineering analysis models, are used to compare the performance of the four surrogate modeling techniques. The results show that the GP-NARX surrogate model tends to have more stable performance than the other three deep learning (DL)-based methods for the three particular examples studied. The tuning effort of GP-NARX is also much lower than its deep learning-based counterparts.

42 ENGINEERING↗

Quantifying heavy quark transport coefficients with an improved transport model

We report that the heavy-flavor transport coefficients contain important information on the strong interaction at finite temperatures. The extraction of these numbers from experimental data requires dynamical modeling of heavy-flavor transport that is coupled to realistic medium evolution. Furthermore, meaningful extractions necessitate both a faithful implementation of the physical inputs to be tested and the quantification of model uncertainty. For these purposes, we have developed a partonic transport model LIDO. It has an improved treatment of in-medium parton bremsstrahlung, which has been calibrated to theoretical calculations in a simple medium to reduce modeling uncertainty. Regarding the interaction between heavy quark and the medium, few-body perturbative scatterings are applied to large-momentum transfer (q) processes, while a diffusion equation models the dynamics of small-q processes. Such a separation restricts the explicit use of medium quasi-particles to large-q processes only. Another advantage is that deviations from the leading-order probe-medium coupling can be parametrized as an additional contribution to the diffusion constant. The heavy quark transport coefficients are then extracted with uncertainty estimation from a Bayesian analysis including both the RHIC and the LHC data. The results are found to be consistent with earlier extraction of the light-quark transport coefficients at high momentum and be comparable with lattice calculations of the heavy-flavor diffusion constant in the static limit at low momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Scalable modular dynamic molten salt reactor system model with decay heat

Nodal dynamic models allow rapid reduced-order simulations of complex systems. This paper presents a publicly available system dynamic model of a molten salt reactor system with a topology of the Molten Salt Reactor Experiment (MSRE), in MATLAB/Simulink, building upon a model previously published and validated with MSRE data. It adds scaling of nominal power, dynamic representation of decay heat generation, a generic decay heat removal system, improved documentation, modularity, and test cases. One hour of system transient runs in less than 10 min on a laptop. This model is intended for first order engineering calculations, developing insights into system transient behavior, equipment sizing, parameter sensitivity studies. The model contains an easy to use toolkit which can be customized and expanded to describe any molten salt reactor system. Finally, scenarios involving off-normal transients are presented. Relation between shutdown reactivity insertions, decay heat removal, core recriticality, and safe shutdown conditions are discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Conversion degree and heat transfer in the cold cap and their effect on glass production rate in an electric melter

A predictive model of melt rate in waste glass vitrification operations is needed to inform melter operations during normal and off-normal operations. This paper describes the development of a model of the cold cap (the reacting melter feed floating on molten glass in a glass melter) that couples heat transfer with the feed-to-glass conversion kinetics. The model was applied to four melter feeds designed for high-level and low-activity nuclear waste feeds using the material properties, either measured or estimated, to obtain temperature and conversion distribution within the cold cap. The cold cap model, when coupled with a computational fluid dynamics model of a Joule-heated glass melter, allows the prediction of the glass production rate and power consumption. The results show reasonable agreement with the melting rates measured during pilot-scale melter tests.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Modeling of the Coral Microbiome: the Influence of Temperature and Microbial Network

Host-associated microbial communities are shaped by extrinsic and intrinsic factors to the holobiont organism. Environmental factors and microbe-microbe interactions act simultaneously on the microbial community structure, making the microbiome dynamics challenging to predict. The coral microbiome is essential to the health of coral reefs and sensitive to environmental changes. Here, we develop a dynamic model to determine the microbial community structure associated with the surface mucus layer (SML) of corals using temperature as an extrinsic factor and microbial network as an intrinsic factor. The model was validated by comparing the predicted relative abundances of microbial taxa to the relative abundances of microbial taxa from the sample data. The SML microbiome from Pseudodiploria strigosa was collected across reef zones in Bermuda, where inner and outer reefs are exposed to distinct thermal profiles. A shotgun metagenomics approach was used to describe the taxonomic composition and the microbial network of the coral SML microbiome. By simulating the annual temperature fluctuations at each reef zone, the model output is statistically identical to the observed data. The model was further applied to six scenarios that combined different profiles of temperature and microbial network to investigate the influence of each of these two factors on the model accuracy. The SML microbiome was best predicted by model scenarios with the temperature profile that was closest to the local thermal environment, regardless of the microbial network profile. Our model shows that the SML microbiome of P. strigosa in Bermuda is primarily structured by seasonal fluctuations in temperature at a reef scale, while the microbial network is a secondary driver. Coral microbiome dysbiosis (i.e., shifts in the microbial community structure or complete loss of microbial symbionts) caused by environmental changes is a key player in the decline of coral health worldwide. Multiple factors in the water column and the surrounding biological community influence the dynamics of the coral microbiome. However, by including only temperature as an external factor, our model proved to be successful in describing the microbial community associated with the surface mucus layer (SML) of the coral P. strigosa. The dynamic model developed and validated in this study is a potential tool to predict the coral microbiome under different temperature conditions.

59 BASIC BIOLOGICAL SCIENCES↗

Estimated attribution of the RFS program on soybean biodiesel in the U.S. using the bioenergy scenario model

Biofuels are supported by various governmental policies in the U.S. and globally as an alternative transportation fuel for environmental, geopolitical, and economic reasons. Much debate surrounds the effectiveness of these policies as well as the overall net environmental effect of increased biofuels use. In the U.S., recent studies have shown that the Renewable Fuels Standard (RFS) Program, overall, may not have been the leading driver of the ethanol industry from 2005 to 2020, contrary to common perception. Similar scrutiny has not been applied to biodiesel. Here, this study uses the Bioenergy Scenario Model, a well-vetted system dynamics model, to retrospectively apportion historical biodiesel production between the RFS Program and other potentially influential drivers, such as the economics of biodiesel vs. diesel, the Biodiesel Tax Credit (BTC), California's Low Carbon Fuel Standard, and other factors. From 2002 to 2020 about 36% of biodiesel production can be attributed to the RFS Program, 35% to the BTC, and the rest to other factors. Thus, the overall effect of the RFS Program appears much larger on biodiesel than on corn ethanol. The finding that the same policy may have very disparate effects on different biofuels helps inform the design of future policies on biofuels.

09 BIOMASS FUELS↗

Impact of grain boundary and surface diffusion on predicted fission gas bubble behavior and release in UO 2 fuel

In this work, we quantify the impact of grain boundary (GB) and surface diffusion on fission gas bubble evolution and fission gas release in UO 2 nuclear fuel using simulations with a hybrid phase field/cluster dynamics model. Here, we begin with a comprehensive literature review of uranium vacancy and xenon atom diffusivity in UO 2 through the bulk, along GBs, and along surfaces. In our model we represent fast GB and surface diffusion using a heterogeneous diffusivity that is a function of the order parameters that represent bubbles and grains. We find that the GB diffusivity directly impacts the rate of gas release via GB transport, and that the GB diffusivity is likely below 104 times the lower value from Olander and van Uffelen. We also find that the surface diffusivity impacts bubble coalescence and mobility, and that the bubble surface diffusivity is likely below 10 -4 times the value from Zhou and Olander.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Effect of injected flux and current temporal phasing on self-organization in the HIT-SI3 experiment

The HIT-SI3 device at the University of Washington uses three oscillating inductive helicity injectors to form and sustain spheromak plasma equilibria. By adjusting the temporal phase of the injector waveforms with respect to each other, the toroidal spectrum of the imposed perturbations can be controlled. Using a recently implemented GPU-based control system, the available mode spectra were explored experimentally by scanning the space of relative injector phasing. In this space, significant variation in the toroidal mode spectrum ($n$ = 1, 2, 3) of the perturbations was observed. Additionally, variation in characteristics of driven equilibria was also observed, including a ≈ $30$% range in toroidal current gain ($I$ $Φ$ / $I$ $Inj$ ). Experimental results are compared with both a composite-equilibria and nonlinear dynamic model, including extended MHD simulations using the NIMROD code and composite Taylor state equilibria computed using the PSI-Tet code. In conclusion, qualitative agreement is seen with the nonlinear models, but not with composite-equilibria models, suggesting the use of nonlinear models to better capture observed plasma dynamics and provide predictive use for future experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamic life-cycle carbon analysis for fast pyrolysis biofuel produced from pine residues: implications of carbon temporal effects

Abstract Background Woody biomass has been considered as a promising feedstock for biofuel production via thermochemical conversion technologies such as fast pyrolysis. Extensive Life Cycle Assessment studies have been completed to evaluate the carbon intensity of woody biomass-derived biofuels via fast pyrolysis. However, most studies assumed that woody biomass such as forest residues is a carbon–neutral feedstock like annual crops, despite a distinctive timeframe it takes to grow woody biomass. Besides, few studies have investigated the impacts of forest dynamics and the temporal effects of carbon on the overall carbon intensity of woody-derived biofuels. This study addressed such gaps by developing a life-cycle carbon analysis framework integrating dynamic modeling for forest and biorefinery systems with a time-based discounted Global Warming Potential (GWP) method developed in this work. The framework analyzed dynamic carbon and energy flows of a supply chain for biofuel production from pine residues via fast pyrolysis. Results The mean carbon intensity of biofuel given by Monte Carlo simulation across three pine growth cases ranges from 40.8–41.2 g CO 2 e MJ −1 (static method) to 51.0–65.2 g CO 2 e MJ −1 (using the time-based discounted GWP method) when combusting biochar for energy recovery. If biochar is utilized as soil amendment, the carbon intensity reduces to 19.0–19.7 g CO 2 e MJ −1 (static method) and 29.6–43.4 g CO 2 e MJ −1 in the time-based method. Forest growth and yields (controlled by forest management strategies) show more significant impacts on biofuel carbon intensity when the temporal effect of carbon is taken into consideration. Variation in forest operations and management (e.g., energy consumption of thinning and harvesting), on the other hand, has little impact on the biofuel carbon intensity. Conclusions The carbon temporal effect, particularly the time lag of carbon sequestration during pine growth, has direct impacts on the carbon intensity of biofuels produced from pine residues from a stand-level pine growth and management point of view. The carbon implications are also significantly impacted by the assumptions of biochar end-of-life cases and forest management strategies.

09 BIOMASS FUELS↗

Rapidly ramp cryogenic air separation unit without loss of O2 product purity—application for low-carbon fossil-fuel plants

The rapid integration of intermittent renewable sources into the electricity grid is driving the need for low-carbon, fossil-fuel power plant capable of rapid ramping. Therefore, a cryogenic air separation unit (ASU) as part of low-carbon, fossil-fuel power plant should be capable of rapid ramping. However, highly integrated and nonlinear processes of the ASU would significantly restrict this rapid ramping. To overcome this fundamental issue, we study the basic dynamic process of a state-of-the-art double-column ASU. A flow-driven dynamic model was established in Aspen Plus Dynamics, assuming perfect flow controls, to capture the basic dynamics of ramping ASU. We find the vapor-liquid counter-current flow structure in the low-pressure column is critical to the air separation when rapidly ramping the ASU. This flow structure is established based on a complicated heat integration process. It is simplified as an apparent counter-current heat transfer process in this work, which greatly reduces the complexity for studying the dynamics of ASU. For keeping this basic flow structure, the heat integration is maintained as its heat duty follows the ASU load, through several basic feed-forward and feed-back controllers on the critical stream flowrates. Moreover, we find a fundamental issue of mismatched dynamics in the low-pressure column, resulting in a significant loss of O2 product purity when rapidly ramping down the ASU and an obvious fluctuation of O2 product purity. Based on these explorations, we propose a basic control method for rapidly ramping the ASU without loss of O2 product purity. Results show the ASU basic dynamic process successfully ramps at a rate up to 10%/min (40-100% load) while maintaining the O2 product purity within 95.2-95.6 mol%.

Cheng, Mao↗

MRCI Subtask 2.4/2.5: Regional/Subregional Analysis and Risk Assessment Final Technical Summary Report

The objective of the Midwest Regional Carbon Initiative (MRCI) project is to implement a collaborative Regional Initiative (RI) to accelerate the deployment of carbon capture, and storage (CCS) in the Midwest-Northeastern quadrant of the United States. This report is a Technical Summary report describing work performed on Tasks 2.4 (Conducting Regional/Subregional Analysis) and 2.5 (Assessing and Managing Risk) during the MRCI project. In Task 2.4, detailed numerical reservoir simulation models were developed for selected carbon storage (CS) systems identified under Task 2.1 (Battelle, 2021a) in the MRCI study area. The objective of Task 2.4 is to demonstrate a dynamic modeling methodology for evaluating the suitability of the (selected) CS systems in the MRCI region for hosting a commercial-scale storage project. In this study, CO2 injectivity was evaluated for different CS systems with an annual injection rate of 1 million metric tonnes (MMT) of CO2 considered as the minimum requirement for a commercial scale project. The objective of Task 2.5 is to assess key risks associated with storing CO2 in the different CS systems across the MRCI region and to demonstrate a method(s) for assessing these risks that may be used by developers of future CO2 storage projects in the region. The risk analysis was limited to evaluating two types of leakage risks (i.e., wellbore leakage, flow across unfractured caprock) at three modeled sites considered in Task 2.4.

MRCI,Report,Summary,Technical Challenges,dynamic m↗

Developing of Quaternary Pumped Storage Hydropower for Dynamic Studies

Quaternary pumped storage hydropower (Q-PSH) technology, as one of the new advanced-PSH technology, has been developed by taking advantage of Conventional-PSH (C-PSH) and Adjustable Speed-PSH (AS-PSH). By combining adjustable-speed pump unit and conventional hydropower turbine unit in the quaternary configuration, Q-PSH has the more competitive capability of providing fast power support in the future high renewable penetrated power system. Acting as energy storage (ES), Q-PSH provides promising power supply to deal with the uncertainty and variability from renewable energy generation. This paper focuses on the dynamic modeling of Q-PSH technology employing full-converter machine and the impact of Q-PSH on the frequency response in a system. The detailed model of Q-PSH is developed and implemented in the IEEE 14-bus system based on GE Positive Sequence Load Flow (PSLF) platform, which captures the dynamic of multiple operation modes, especially hydraulic short-circuit (HSC) operation mode. Several cases are set up to reveal the advantages of Q-PSH technology when power electronic based renewable energy generation is deployed in the system. Sensitivity studies of the controller in pump governor show the impact of parameters in pump response performance. The comparison case illustrates the impact of frequency response provided by the Q-PSH in the system.

dynamic modeling↗

Shadow molecular dynamics for flexible multipole models

Shadow molecular dynamics provide an efficient and stable atomistic simulation framework for flexible charge models with long-range electrostatic interactions. Shadow molecular dynamics simulations are driven by approximate “shadow” Born–Oppenheimer potentials for which the exact charges and forces are directly accessible without relying on costly (and approximate) iterative solvers. While previous implementations have been limited to atomic monopole charge distributions, we extend this approach to flexible multipole models. We derive detailed expressions for the shadow energy functions, potentials, and force terms, explicitly incorporating monopole–monopole, dipole–monopole, and dipole–dipole interactions. In our formulation, both atomic monopoles and atomic dipoles are treated as extended dynamical variables alongside the propagation of the nuclear degrees of freedom. We demonstrate that introducing the additional dipole degrees of freedom preserves the stability and accuracy previously seen in monopole-only shadow molecular dynamics simulations. In addition, we present a shadow molecular dynamics scheme where the monopole charges are held fixed while the dipoles remain flexible. Our extended shadow dynamics provide a framework for stable, computationally efficient, and versatile molecular dynamics simulations involving long-range interactions between flexible multipoles. This is of particular current interest in combination with machine-learned interatomic potentials, including long-range electrostatic interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Droplet bioprinting of acellular and cell-laden structures at high-resolutions

Advances in digital light projection(DLP) based (bio) printers have made printing of intricate structures at high resolution possible using a wide range of photosensitive bioinks. A typical setup of a DLP bioprinter includes a vat or reservoir filled with liquid bioink, which presents challenges in terms of cost associated with bioink synthesis, high waste, and gravity-induced cell settling, contaminations, or variation in bioink viscosity during the printing process. Here, we report a vat-free, low-volume, waste-free droplet bioprinting method capable of rapidly printing 3D soft structures at high resolution using model bioinks and model cells. A multiphase many-body dissipative particle dynamics model was developed to simulate the dynamic process of droplet-based DLP printing and elucidate the roles of surface wettability and bioink viscosity. Process variables such as light intensity, photo-initiator concentration, and bioink formulations were optimized to print 3D soft structures (∼0.4–3 kPa) with a typical layer thickness of 50 µm, an XY resolution of 38 ± 1.5 μm and Z resolution of 237 ± 5.4 µm. To demonstrate its versatility, droplet bioprinting was used to print a range of acellular 3D structures such as a lattice cube, a Mayan pyramid, a heart-shaped structure, and a microfluidic chip with endothelialized channels. Droplet bioprinting, performed using model C3H/10T1/2 cells, exhibited high viability (90%) and cell spreading. Additionally, microfluidic devices with internal channel networks lined with endothelial cells showed robust monolayer formation while osteoblast-laden constructs showed mineral deposition upon osteogenic induction. Overall, droplet bioprinting could be a low-cost, no-waste, easy-to-use, method to make customized bioprinted constructs for a range of biomedical applications.

DLP↗

BSM (Bioenergy Scenario Model) 2023 FKA: Biomass Scenario Model [SWR-09-09]

The U.S. Department of Energy's (DOE's) Bioenergy Technologies Office and the National Renewable Energy Laboratory (NREL) developed the BSM (Bioenergy Scenario model) to explore the development of a U.S. biofuels industry. The BSM is a system dynamics model built on the STELLA software platform. The model represents the dynamic interactions of the major sectors of the biofuels industry—feedstock production, feedstock logistics, biomass to biofuels conversion, and biofuels end use, including fuels inventory, dispensing, distribution, fuel use, and the vehicle fleet. The BSM represents contextual aspects of the developing biofuels industry, including investment in new biomass to biofuel conversion technologies, competition from petroleum fuels, vehicle demand for biofuels, and various government policies, using all of these to simulate the development of the industry. The purpose of the BSM is to generate and explore plausible scenarios for the evolution of a biofuels industry in the United States, and as a high-level system model it is not designed for precise, quantitative forecasting. Instead, it is best used to (1) analyze and evaluate alternate policies; (2) generate scenarios; (3) identify high-impact levers and bottlenecks to system evolution; and (4) seed focused discussion among policymakers, analysts, and stakeholders.

Bush, Brian↗

New Nonreactive Force Field for Accurate Molecular Dynamics Simulations of TATB at Extreme Conditions

Insensitive high explosives based on TATB (1,3,5-triamino-2,4,6-trinitrobenzene) are needed in applications when safety is of paramount importance, but the basic material properties that give rise to its insensitivity are not fully understood. Molecular dynamics modeling using empirical force fields (FFs) has been the main route to characterize many complicated dynamical properties of TATB single crystal, but these FFs have not been comprehensively tested at extreme conditions typical of detonation. We collect a benchmark data set of (quasi)static TATB physical properties as determined by experiments and electronic structure calculations and apply this data set to validate four existing TATB FFs along with a new TATB FF that we develop here and denote as the CEA-LLNL-Missouri (CLM) FF. Benchmark data include vibrational spectra, the TATB crystal temperature–pressure–volume equation of state and lattice parameters, properties of TATB crystal polymorphs and transitions to the gaseous and liquid states, dimer energy landscapes, the pressure-dependent elastic tensor, and the energy landscape for inelastic deformation via sliding of TATB crystal layers. As a general assessment, we find that the two existing nonreactive FFs are more accurate in describing TATB’s physical properties compared to the two variants of the ReaxFF reactive FF considered. The new CLM FF is found to consistently yield similar or better agreement with experiments and electronic structure theory than any of the existing FF models, and it presents a distinct improvement in accurately modeling TATB elasticity and equation of state. So this work is expected to help improve the accuracy of FF-based modeling of complicated dynamic responses that ultimately govern the safety and performance characteristics of this material.

36 MATERIALS SCIENCE↗

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗

Network-Wide Traffic Signal Control Using Bilinear System Modeling and Adaptive Optimization

This study proposes a new multi-input multi-output optimal bilinear signal control method in which a bilinear dynamic model approximation is used to capture the nonlinear dynamics of the urban traffic networks. With signal green time splits as the control input and traffic delay changes as the output for each intersections in the network, a bilinear system model was developed, which, on the basis of linear system modeling, takes interactions among traffic delays and signal timing splits into consideration. Based on the bilinear system modeling framework, we conducted two steps in each time interval to derive traffic control strategies: (1) we used the normalized least-squared algorithm to estimate system parameters; and (2) we solved an online optimization problem to obtain the updated traffic control inputs for the signal timing that minimizes future traffic delays. We evaluated the proposed method in a microscopic traffic simulation environment (VISSIM) with a 35-intersection network of Bellevue city in Washington. Two different traffic demand patterns: (1) normal traffic demands; and (2) time-varying traffic demands were simulated to compare the performance of different control strategies. Experimental results show that (1) the proposed bilinear system model can better describe traffic system dynamics than linear-model based methods, such as our previously developed linear-quadratic regulator control; and (2) the proposed method outperforms the state-of-the-art signal control strategies, namely the max-pressure and the self-organizing traffic light control methods. We have also shown that the proposed method is applicable to all other possible network layouts and signal controller phasing structures.

42 ENGINEERING↗