Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transient dynamics”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Demonstrating the Transient System Impact of Cyber-Physical Events Through Scalable Transmission and Distribution (T&D) Co-Simulation

Modern power systems become more vulnerable to cyber threats due to their growing interconnectivity, interdependence, and complexity. Widespread deployment of distributed energy resources (DERs) further expands the threat landscape to the grid edge, where fewer cybersecurity protections exist. In this article, a systematic cyber-physical events demonstration, enabled by an integrated transmission, distribution, and communication co-simulation framework, is performed. It analyzes cyber risks to power grid under DER-enabled automatic generation control from different angles. Unlike existing works, the simulation captures millisecond-to-minutes frequency and voltage transient dynamics at a cross-region system scale.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Evaluation of Steam Cycle Upgrades to Improve the Competitiveness of U.S. Coal Power Plants (Final Scientific / Technical Report)

Increasing the competitiveness of the existing pulverized-coal utility fleet in the United States may be achieved by decreasing heat rate, via increases in steam cycle efficiency through upgraded steam temperatures and use of latest technology available in steam turbine and blading design. The average net plant efficiency of the US coal-fired fleet is approximately 33% (HHV). Plant efficiency increases to approximately 41.4% (HHV) at 1,350°F (732°C) steam temperature. However, achieving these Advanced Ultra-Super Critical (AUSC) steam conditions requires the use of advanced high-temperature materials. While there has been a significant amount of DOE-funded materials R&D, most of the related design work has focused on new (greenfield) units, rather than on opportunities to retrofit this advanced technology to the existing utility fleet. If technology, based upon the advanced materials required for AUSC steam conditions, may be applied to the existing fleet, using an economically viable retrofit, a higher capacity factor can be expected as a result of the increased plant competitiveness. The Electric Power Research Institute (EPRI) was awarded a project by the US Department of Energy to examine the technical and economic feasibility of a series of steam cycle upgrades to the two most prevalent types of U.S. coal power plants: 2,400 psig (16.6 MPa) subcritical and 3,500 psig (24.1 MPa) supercritical pulverized coal units. The nine upgrade options that were originally being considered included increasing the main and reheat steam temperatures from 1,000°F (538°C) to 1,100°, 1,200°, and 1,350°F (593°C, 649°C, and 732°C) while holding the steam pressures constant at their original design values, and cases where just the main steam or reheat steam temperatures were increased. The objective was to minimize the modifications required to the existing power plant while still providing a significant improvement in heat rate. The upgrade options assumed that the boiler enclosure envelope remained unchanged from each base case, and that all applicable OEM design guidelines for normal commercial units were imposed. For the highest temperature supercritical case, an option of using a low-pressure molten salt loop to transfer heat from the furnace to the steam was examined. The first major task of the work scope was designed to examine the technical feasibility of various upgrade options, while the subsequent work determined economic viability of the technically feasible upgrade options. Prior to evaluating the effect of these increased temperatures, a “base case” model of a subcritical and supercritical PC boiler was created, which was used for comparative purposes. Upgrade options were evaluated at full-load, part-load and dynamic transient conditions. Once the technical feasibility of each upgrade option was evaluated, the economic value of the heat rate improvement of each feasible option was determined by detailed modeling of unit dispatch in several regional power markets. The dispatch model was used to estimate the amount of revenue from power sales the upgraded unit would receive in comparison to a non-upgraded version of the same power plant. As a parallel task to the dispatch analysis, the capital cost of implementing the upgrades was estimated. The capital cost estimates were then compared to the increased revenue estimated by the dispatch modeling to determine the economic attractiveness of each upgrade option. Several upgrade options were determined to be technically feasible. The net present value (NPV) of the costs for steam cycle upgrades considered in this study ranged from approximately $\$$111 to $\$$130 million. The economic modeling results show that the unit dispatch changes resulting from steam cycle upgrades are relatively small, due largely to heat rate (and operating cost) changes being relatively small. Additionally, the cost of each upgrade exceeds the net revenue increases associated with the upgrade case. Note that the breakeven values are higher for subcritical retrofits, but the capital costs for the subcritical upgrades are also slightly higher. In typical new pulverized coal plants, fuel accounts for approximately 25% of the cost of electricity (COE), while capital costs represent around 50% of the COE. Therefore, in order to improve the heat rate by 4% one can only afford to increase the capital cost by 2%, at the same cost of electricity. The conclusion of this study is that without a cost for emitting CO 2 , it will be difficult to pay for significant efficiency improvements on plants firing low cost coals.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Impact of Integration Scheme on Performance of Anisotropic Plasticity Models

Given the prevalent role of metals in a variety of industries, schemes to integrate corresponding constitutive models in finite element applications have long been studied. A number of formulations have been developed to accomplish this task; each with their own advantages and costs. Often the focus has been on ensuring the accuracy and numerical stability of these algorithms to enable robust integration. While important, emphasis on these performance metrics may often come at the cost of computational expense potentially neglecting the needs of individual problems. In the current work, the performance of two of the most common integration methods for anisotropic plasticity -- the convex cutting plane (CCP) and closest point projection (CPP) -- across a variety of metrics is assessed; including accuracy and cost. A variety of problems are considered ranging from single elements to large representative simulations including both implicit quasistatic and explicit transient dynamic type responses. The relative performance of each scheme in the different instances is presented with an eye towards guidance on when the different algorithms may be beneficial.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Sierra/SolidMechanics 5.6 User's Manual

Sierra/SolidMechanics (Sierra/SM) is a three-dimensional solid mechanics code with a versatile element library, nonlinear material models, large deformation capabilities, and contact. It is built on the SIERRA Framework. SIERRA provides a data management framework in a parallel computing environment that allows the addition of capabilities in a modular fashion. Contact capabilities are parallel and scalable. This document provides information about the functionality in Sierra/SM and the command structure required to access this functionality in a user input file. This document is divided into chapters based primarily on functionality. For example, the command structure related to the use of various element types is grouped in one chapter; descriptions of material models are grouped in another chapter. Sierra/SM provides both explicit transient dynamics and implicit quasistatics and dynamics capabilities. Both the explicit and implicit modules are highly scalable in a parallel computing environment. In the past, the explicit and implicit capabilities were provided by two separate codes, known as Presto and Adagio, respectively. These capabilities have been consolidated into a single code. The executable is named Adagio, but it provides the full suite of solid mechanics capabilities, for both implicit and explicit. The Presto executable has been disabled as a consequence of this consolidation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Stabilized Hyperfoam Modeling of the General Plastics EF4003 (3 PCF) Flexible Foam

Constitutive model parameterizations for the General Plastics EF4003 low density 3 pound per cubic foot are needed for design and qualification purposes in normal and abnormal mechanical simulations. The material is expected to be deformed in two ways: first during preloading, and second under impact conditions of the system (transient dynamic). All analyses are to be performed at room temperature. The goal is to provide the analysis community a robust constitutive model parameterization to represent the compression behavior of the EF4003 foam from small deformations up to massive compressive deformations when the foam is densifying. It is worth noting the EF4003 exhibits anisotropy in its stress-strain behavior between the rise and transverse directions (See figure 2.8c-d) as well as plateau behavior that is very likely to cause material stability issues, due to the buckling transition, (and has historically done so) when using Sandia’s current workhorse models for flexible foams, Hyperfoam and Flex Foam. A Stability-informed Hyperfoam parameterization procedure is developed and executed to calibrate a hyperfoam model for the EF4003 room temperature, transversely loaded data. A rise orientation parameterization was not attempted due to localization in the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Embedded, Real-Time, and Distributed Traveling Wave Fault Location Method Using Graph Convolutional Neural Networks

This work proposes and develops an implementation of a fault location method to provide a fast and resilient protection scheme for power distribution systems. The method analyzes the transient dynamics of traveling waves (TWs) to generate features using the discrete wavelet transform (DWT), which are then used to train several graph convolutional network (GCN) models. Faults are simulated in the IEEE 34-node system, which is divided into three protection zones (PZs). The goal is to identify the PZ in which the fault occurs. The GCN models create a distributed protection scheme, as all nodes are able to retrieve a prediction. Given that message-passing between nodes occurs both during training and in the execution of the model, the resiliency of such schemes to communication losses was analyzed and demonstrated. One of the models, which only uses voltage measurements, was implemented on a Texas Instruments F28379D development board. The execution times were monitored to assess the speed of the protection scheme. It is shown that the proposed method can be executed in approximately a millisecond, which is comparable to existing TW protection in the transmission system. For experimental purposes, a DWT-based detection method is employed. A design of a setup to playback TWs using two development boards is also addressed.

Jiménez-Aparicio, Miguel (ORCID:000000016864461X)↗
DOE OSTI2025

XRISM Forecast for the Coma Cluster: Stormy, with a Steep Power Spectrum

The XRISM Resolve microcalorimeter array measured the velocities of hot intracluster gas at two positions in the Coma galaxy cluster: ${3}^{{\prime} }\times {3}^{{\prime} }$ squares at the center and at 6$^{\prime} $ (170 kpc) to the south. We find the line-of-sight velocity dispersions in those regions to be σ z = 208 ± 12 km s −1 and 202 ± 24 km s −1 , respectively. The central value corresponds to a 3D Mach number of M = 0.24 ± 0.015 and a ratio of the kinetic pressure of small-scale motions to thermal pressure in the intracluster plasma of only 3.1% ± 0.4%, at the lower end of predictions from cosmological simulations for merging clusters like Coma, and similar to that observed in the cool core of the relaxed cluster A2029. Meanwhile, the gas in both regions exhibits high line-of-sight velocity differences from the mean velocity of the cluster galaxies, Δv z = 450 ± 15 km s −1 and 730 ± 30 km s −1 , respectively. A small contribution from an additional gas velocity component, consistent with the cluster optical mean, is detected along a sight line near the cluster center. The combination of the observed velocity dispersions and bulk velocities is not described by a Kolmogorov velocity power spectrum of steady-state turbulence; instead, the data imply a much steeper effective slope (i.e., relatively more power at larger linear scales). This may indicate either a very large dissipation scale, resulting in the suppression of small-scale motions, or a transient dynamic state of the cluster, where large-scale gas flows generated by an ongoing merger have not yet cascaded down to small scales.

coma cluster↗
DOE OSTI2025

A new biogeochemical modelling framework (FLaMe-v1.0) for lake methane emissions on the regional scale: development and application to the European domain

This study presents a new physical-biogeochemical modelling framework for simulating lake methane (CH 4 ) emissions at regional scales. The new model, FLaMe-v1.0 (Fluxes of Lake Methane), rests on an innovative, computationally efficient lake clustering approach that enables the simulation of CH 4 emissions across a large number of lakes. Building on the Canadian Small Lake Model (CSLM) that simulates the lake physics, we develop a suite of biogeochemical modules to simulate transient dynamics of organic Carbon (C), Oxygen (O 2 ), and CH 4 . We first test the performance of FLaMe-v1.0 by analyzing physical and biogeochemical processes in two theoretical lakes with characteristics that can be considered representative for many lakes (an oligotrophic, deep lake driven by cold climate versus a eutrophic, shallow lake driven by warm climate). Next, we evaluate the model by comparing simulated and observed timeseries of CH 4 emissions in four well-surveyed lakes. We then apply FLaMe-v1.0 at the European scale to evaluate simulated diffusive and ebullitive lake CH 4 fluxes against in-situ measurements in both boreal and central European regions. Finally, we provide a first assessment of the spatio-temporal variability in CH 4 emissions from European lakes with a surface area comprised between 0.1–1000 km 2 (n= 108 407, total area = 1.33 × 105 km 2 ), indicating a total emission of 0.97 ± 0.23 Tg CH 4 yr −1 , with the uncertainty constrained by combining FLaMe-v1.0 and machine learning techniques. Moreover, 30 % and 70 % of these CH 4 emissions are through diffusive and ebullitive pathways, respectively. Annually averaged CH 4 emission rates per unit lake area during 2010–2016 have a South-to-North decreasing gradient, resulting in a mean over the European domain as 7.39 g CH 4 m −2 yr −1 . Our simulations reveal a strong seasonality (with ice-blocking effects accounted for) in European lake CH 4 emissions, with nearly ten times higher emissions during late summer than during winter. This pronounced seasonal variation highlights the importance of accounting for the sub-annual variability in CH 4 emissions to accurately constrain regional CH 4 budgets. In the future, FLaMe-v1.0 could be embedded into Earth System Models to investigate the feedback between climate warming and global lake CH 4 emissions.

Maisonnier, Manon [Free Univ. of Brussels (Belgium↗
DOE OSTI2020

Model-Agnostic Algorithm for Real-Time Attack Identification in Power Grid using Koopman Modes

Malicious activities on measurements from sensors like Phasor Measurement Units (PMUs) can mislead the control center operator into taking wrong control actions resulting in disruption of operation, financial losses, and equipment damage. In particular, false data attacks initiated during power systems transients caused due to abrupt changes in load and generation can fool the conventional model-based detection methods relying on thresholds comparison to trigger an anomaly. In this paper, we propose a Koopman mode decomposition (KMD) based algorithm to detect and identify false data attacks in real-time. The Koopman modes (KMs) are capable of capturing the nonlinear modes of oscillation in the transient dynamics of the power networks and reveal the spatial embedding of both natural and anomalous modes of oscillations in the sensor measurements. The Koopman-based spatio-temporal nonlinear modal analysis is used to filter out the false data injected by an attacker. The performance of the algorithm is illustrated on the IEEE 68-bus test system using synthetic attack scenarios generated on GridSTAGE, a recently developed multivariate spatio-temporal data generation framework for simulation of adversarial scenarios in cyber-physical power systems.

Nandanoori, Sai Pushpak↗
DOE OSTI2021

Integrated transmission electron microscope

An integrated transmission electron microscope comprising multiple electron sources for tuned beams of ultrafast, scanning probe, and parallel illumination in varied beam energies can be alternated within sub-microseconds onto a sample with dynamic ‘transient state’ processes to acquire atomic-scale structural/chemical data with site specificity. The various electron sources and condenser optics enable high-resolution imaging, high-temporal resolution imaging, and chemical imaging, using fast-switching magnets to direct the different electron beams onto a single maneuverable objective pole piece where the sample resides. Such multimodal in situ characterization tools housed in a single microscope have the potential to revolutionize materials science.

Jungjohann, Katherine L.↗
DOE OSTI2023

Efficiency of charge transfer in changing the dissociation dynamics of OD + transients formed after the photo-fragmentation of D 2 O

Here, we present an investigation of the relaxation dynamics of deuterated water molecules after direct photo-double ionization at 61 eV. We focus on the very rare D + + O + + D reaction channel in which the sequential fragmentation mechanisms were found to dominate the dynamics. Aided by theory, the state-selective formation and breakup of the transient OD + (a 1 Δ, b 1 Σ + ) is traced, and the most likely dissociation path—OD + : a 1 Δ or b 1 Σ + → A 3 Π → X 3 Σ – → B 3 Σ – —involving a combination of spin–orbit and non-adiabatic charge transfer transitions is determined. The multi-step transition probability of this complex transition sequence in the intermediate fragment ion is directly evaluated as a function of the energy of the transient OD + above its lowest dissociation limit from the measured ratio of the D + + O + + D and competing D + + D + + O sequential fragmentation channels, which are measured simultaneously. Our coupled-channel time-dependent dynamics calculations reproduce the general trends of these multi-state relative transition rates toward the three-body fragmentation channels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Dynamic processes in transient phases during self-assembly of organic semiconductor thin films

Understanding and manipulating crystallization processes has been an important challenge for solution-processed organic thin films, both for fundamental studies and for fabricating thin films with near-intrinsic charge transport properties. We report an in situ X-ray scattering study of the crystallization of 2-decyl-7-phenyl-[1]benzothieno[3,2-b][1]benzothiophene (Ph-BTBT-C 10 ) during its deposition from solution. At temperatures modestly below the smectic-E/crystalline phase boundary, the crystallization goes through a transient liquid crystal state before reaching the final stable crystalline phase. Significant dynamics occur in the first few seconds of the transition, which are observed through fluctuations in the X-ray scattering intensity, and are correlated with the time interval that the transient thin film coexists with the evaporating solvent. The transition to the stable crystalline phase takes minutes or even hours under these conditions, which may be a result of the asymmetry of the molecule. Transient phases are of potential interest for applications, since they can act as a route to self-assembly of organic thin films. However, our observations show that the long-lived monolayer-stacked intermediate state does not act as a template for the bilayer-stacked crystalline phase. Rather, the grain structure is replaced through nucleation, where the nucleation free-energy barrier is related to a potential barrier that prevents molecules to flip their long axis by 180°.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Transient Relativistic Fluid Dynamics in a General Hydrodynamic Frame

In this work, we propose a new theory of second-order viscous relativistic hydrodynamics which does not impose any frame conditions on the choice of the hydrodynamic variables. It differs from Mueller-Israel-Stewart theory by including additional transient degrees of freedom, and its first-order truncation reduces to Bemfica-Disconzi-Noronha-Kovtun theory. Conditions for causality and stability are explicitly given in the conformal regime. As an illustrative example, we consider Bjorken flow solutions to our equations and identify variables which make a hydrodynamic attractor manifest.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A Data Set for Intercomparing the Transient Behavior of Dynamical Model-Based Subseasonal to Decadal Climate Predictions

Climate predictions using coupled models in different time scales, from intraseasonal to decadal, are usually affected by initial shocks, drifts, and biases, which reduce the prediction skill. These arise from inconsistencies between different components of the coupled models and from the tendency of the model state to evolve from the prescribed initial conditions toward its own climatology over the course of the prediction. Aiming to provide tools and further insight into the mechanisms responsible for initial shocks, drifts, and biases, this paper presents a novel data set developed within the Long Range Forecast Transient Intercomparison Project, LRFTIP. This data set has been constructed by averaging hindcasts over available prediction years and ensemble members to form a hindcast climatology, that is a function of spatial variables and lead time, and thus results in a useful tool for characterizing and assessing the evolution of errors as well as the physical mechanisms responsible for them. A discussion on such errors at the different time scales is provided along with plausible ways forward in the field of climate predictions

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Influence of Nd(II) and Nd(III) Ions on the Speciation and Kinetic Dynamics of Radiolytic Transients in Molten LiCl-KCl Eutectic Salt Mixtures

In the pyrochemical reprocessing of used nuclear fuel, understanding the effects of ionizing radiation on the chemical behavior of fission products in molten chloride eutectic salt is crucial. Here, this study investigates the impact of radiation on the chemical and kinetic speciation of Nd ions in molten LiCl-KCl eutectic, which are challenging to separate due to their multivalent oxidation states and comproportionation reactions. Electron pulse radiolysis techniques were used to determine transient absorption spectra, chemical kinetics, and activation parameters for the reaction of Nd ions with the primary products of molten LiCl-KCl eutectic radiolysis: the solvated electron (e S - ) and dichlorine radical anion (Cl 2 •- ). Initially, Nd(III) reacted rapidly with e S - , forming Nd(II) with a second-order rate coefficient of k = (4.54 ± 0.07) × 10 10 M -1 s -1 at 673 K. The resulting Nd(II) ions then reacted slowly with Cl 2 •- , regenerating Nd(III) with an estimated rate of k = (1.72 ± 0.04) × 10 10 M -1 s -1 at 673 K. Additionally, the comproportionation reaction of Nd(III) and Nd(s) to form Nd(II) was monitored chemically and electrochemically. This work suggests that the comproportionation reaction of Nd(III) with Nd(s) to form Nd(II) is unstable in molten LiCl-KCl eutectic.

36 - MATERIALS SCIENCE↗
DOE OSTI2023

Topology induced crossover between Langevin, subdiffusion, and Brownian diffusion regimes in supercooled water

Despite extensive studies of supercooled water, it remains challenging to understand its peculiar dynamic anomalous properties. In this work, we integrated full atomistic simulations of supercooled water over the temperature range of room temperature to 200 K using a quantum-mechanics-based polarizable force field with the dressed dynamics method that couples fast collision events and slow reorganization dynamics of hydrogen-bond networks. Our analysis unveils the salient multiscale features in the transient relaxation dynamics of supercooled water. A classical Langevin behavior dominates at fast timescales, while long-time relaxations unveil two different activation barriers in two temperature regions: below and above 230 K. The modulation of the entropy spectrum by temperature is elucidated in terms of a three-state model underlined by the complexity of the water dynamics associated with a topological transition of a strong hydrogen-bond network. This state-dependent network topology is quantitatively characterized by power-law exponents of inverse network connectivity from 200 to 298 K. Finally, this work provides valuable guidance for further studies on the transient relaxation dynamics of supercooled water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗