Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transient 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 847 records · Page 47

Efficient Free Triplet Generation Follows Singlet Fission in Diketopyrrolopyrrole Polymorphs with Goldilocks Coupling

Microcrystal electron diffraction, grazing incidence wide-angle scattering, and UV–vis spectroscopy were used to determine the unit-cell structure and the relative composition of dimethylated diketopyrrolopyrrole H- and J-polymorphs within thin films subjected to vapor solvent annealing (VSA) for different times. The electronic structure and excited-state deactivation pathways of the different polymorphs were examined by transient absorption spectroscopy, conductive probe atomic force microscopy, and molecular modeling. We find that VSA initially converts amorphous films into mixtures of H- and J-polymorphs and promotes further conversion from H to J with longer VSA times. Here, though both polymorphs exhibit efficient SF to form coupled triplets, free triplet yields are higher in J-polymorph films compared to mixed films because coupling in J-aggregates is lower and, in turn, more favorable for triplet decoupling.

36 MATERIALS SCIENCE↗

Climate forcing controls on carbon terrestrial fluxes during shale weathering

Climate influences near-surface biogeochemical processes and thereby determines the partitioning of carbon dioxide (CO 2 ) in shale, and yet the controls on carbon (C) weathering fluxes remain poorly constrained. Using a dataset that characterizes biogeochemical responses to climate forcing in shale regolith, we implement a numerical model that describes the effects of water infiltration events, gas exchange, and temperature fluctuations on soil respiration and mineral weathering at a seasonal timescale. Our modeling approach allows us to quantitatively disentangle the controls of transient climate forcing and biogeochemical mechanisms on C partitioning. We find that ~3% of soil CO 2 (1.02 mol C/m 2 /y) is exported to the subsurface during large infiltration events. Here, net atmospheric CO 2 drawdown primarily occurs during spring snowmelt, governs the aqueous C exports (61%), and exceeds the CO 2 flux generated by pyrite and petrogenic organic matter oxidation (~0.2 mol C/m 2 /y). We show that shale CO 2 consumption results from the temporal coupling between soil microbial respiration and carbonate weathering. This coupling is driven by the impacts of hydrologic fluctuations on fresh organic matter availability and CO 2 transport to the weathering front. Diffusion-limited transport of gases under transient hydrological conditions exerts an important control on CO 2 (g) egress patterns and thus must be considered when inferring soil CO 2 drawdown from the gas phase composition. Our findings emphasize the importance of seasonal climate forcing in shaping the net contribution of shale weathering to terrestrial C fluxes and suggest that warmer conditions could reduce the potential for shale weathering to act as a CO 2 sink.

58 GEOSCIENCES↗

The glass-ceiling convective regime and the origin and diversity of coronae on Venus

Venus and Earth are rocky planets of roughly the same size and bulk density, yet their surface volcanic and tectonic features appear substantially different. On Venus, the coexistence of large volcanic highlands—interpreted as the surface expression of long-lived mantle plumes—alongside coronae, smaller features thought to be caused by transient thermal diapirs, remains enigmatic. Using two-dimensional numerical models of mantle convection with sharp and broad mineral phase transitions for pyrolite, we show that both scales of upwellings can be generated in a stagnant lid planet with an interior temperature 250 to 400 K warmer than Earth’s. The smaller plumes originate from a ~600 km deep internal layer that exists as a consequence of the different sequence of mineral phase transitions that occur in warmer mantles less processed and differentiated by partial melting and volcanism. Future models that include melting will provide further tests of our hypothesis.

Science & Technology - Other Topics↗

Species Transport Framework Development in SAM for System-Level Tritium Source Term Analysis

The SAM code is under development as a modern system-level modeling and simulation tool for advanced non–light water reactor safety analyses, with recent efforts to add capabilities to evaluate radiological source term risks in these novel reactor concepts. By leveraging the established system-level multiphysics thermal-hydraulic models in SAM, a framework for tightly coupled species transport modeling has been integrated into the code for engineering-scale source term evaluation. This species transport framework was first applied to the simulation of tritium, which is a well-known source term in conventional light water reactors. Tritium poses a unique risk in salt-cooled reactors, especially those with lithium-bearing salts such as the fluoride salt–cooled high-temperature reactor (FHR) concept, as tritium is generated in the salt coolant in significant quantities due to neutron interactions. A compounding factor is the increased mobility of tritium at high temperatures, which is able to permeate through metals while also potentially being retained in graphite pebbles and structures. Engineering-scale models for the tritium transport pathways in a FHR have been developed using the new species transport framework in SAM. The capabilities are assessed through analytical verification problems and validated with data from a graphite retention experiment. In conclusion, the system-level model is demonstrated by performing an initial estimate of baseline tritium generation and flows in a generic reference SAM FHR model, setting a foundation for future studies of source term transient analysis with the potential for further multiscale and multiphysics integration.

SAM↗

Spontaneous skyrmion conformal lattice and transverse motion during dc and ac compression

Abstract We use atomistic-based simulations to investigate the behavior of ferromagnetic skyrmions being continuously compressed against a rigid wall under dc and ac drives. The compressed skyrmions can be annihilated close to the wall and form a conformal crystal with both a size and a density gradient, making it distinct from conformal crystals observed previously for superconducting vortices and colloidal particles. For both dc and ac driving, the skyrmions can move transverse to the compression direction due to a combination of density and size gradients. Forces in the compression direction are converted by the Magnus force into transverse motion. Under ac driving, the amount of skyrmion annihilation is reduced and we find a skyrmion Magnus ratchet pump. We also observe shear banding in which skyrmions near the wall move up to twice as fast as skyrmions further from the wall. When we vary the magnitude of the applied drive, we find a critical current above which the skyrmions are completely annihilated during a time scale that depends on the magnitude of the drive. By varying the magnetic parameters, we find that the transverse motion is strongly dependent on the skyrmion size. Smaller skyrmions are more rigid, which interferes with the size gradient and destroys the transverse motion. We also confirm the role of the size gradient by comparing our atomistic simulations with a particle-based model, where we find that the transverse motion is only transient. Our results are relevant for applications where skyrmions encounter repulsive magnetic walls, domain walls, or interfaces.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Time-dependent SOLPS-ITER simulations of the tokamak plasma boundary for model predictive control using SINDy *

Abstract Time-dependent SOLPS-ITER simulations have been used to identify reduced models with the sparse identification of nonlinear dynamics (SINDy) method and develop model-predictive control of the boundary plasma state using main ion gas puff actuation. A series of gas actuation sequences are input into SOLPS-ITER to produce a dynamic response in upstream and divertor plasma quantities. The SINDy method is applied to identify reduced linear and nonlinear models for the electron density at the outboard midplane n e , s e p O M P and the electron temperature at the outer divertor T e , s e p d i v . Note that T e , s e p d i v is not necessarily the peak value of T e along the divertor. The identified reduced models are interpretable by construction (i.e. not black box), and have the form of coupled ordinary differential equations. Despite significant noise in T e , s e p d i v , the reduced models can be used to predict the response over a range of actuation levels to a maximum deviation of 0.5% in n e , s e p O M P and 5%–10% in T e , s e p d i v for the cases considered. Model retraining using time history data triggered by a preset error threshold is also demonstrated. A model predictive control strategy for nonlinear models is developed and used to perform feedback control of a SOLPS-ITER simulation to produce a setpoint trajectory in n e , s e p O M P using the integrated plasma simulator framework. The developed techniques are general and can be applied to time-dependent data from other boundary simulations or experimental data. Ongoing work is extending the approach to model identification and control for divertor detachment, which will present transient nonlinear behavior from impurity seeding, including realistic latency and synthetic diagnostic signals derived from the full SOLPS-ITER output.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Detecting anomalous SRF cavity behavior with unsupervised learning

We present an unsupervised learning framework for detecting anomalous superconducting radio-frequency (SRF) cavity behavior at the Continuous Electron Beam Accelerator Facility (CEBAF), emphasizing its initial performance and effectiveness. Key to the system’s success was the development of data acquisition systems (DAQs) that capture fast-sampled, information-rich signals, essential for detecting transient effects. The approach involves creating daily cavity-specific models using principal component analysis to handle variations in rf signal behavior and mitigate performance degradation from data drift. This unsupervised method eliminates the need for expensive labeling by continuously updating models with recent data. Deployed and operational for 3 months before a scheduled shutdown, the system successfully identified several issues with DAQ signals, confirming its effectiveness. Despite access to only a fraction of CEBAF’s SRF cavity signals, the framework efficiently detected several instances requiring intervention, demonstrating a significant improvement over traditional, labor-intensive methods of manual plot inspection. Published by the American Physical Society 2025

43 PARTICLE ACCELERATORS↗

High-Performance Computing Based EMT Simulation: Power Grid with IBRs

Electromagnetic transient (EMT) simulation of power grids with high-fidelity models of inverter-based resources (IBRs) is time-consuming and difficult to scale. The necessity for high-fidelity models of IBRs that incorporate the dynamics of individual inverters within IBRs has been showcased in recent studies. These studies focused on events with partial power reduction in each IBR during a transmission line fault in the power grid. These types of events have been documented in multiple North American Electric Reliability Council (NERC) reports in the past decade. It is imperative then to find solutions to speed-up EMT simulations and scale the size of the region with IBRs studied in EMT simulations. In this paper, a combination of numerical simulation algorithms with high-performance computing techniques are employed in discretization and linear solvers employed in the proposed RE-INTEGRATE EMT simulation platform for power grid with IBRs. For ease of scalability, modular and object-oriented programming is used as these techniques are implemented. Additionally, automation software is developed to convert legacy software codes to the proposed RE-INTEGRATE EMT simulation platform. Thereafter, this platform is evaluated on multi-core central processing units (CPUs). Finally, scale-up tests are performed to showcase the scalability that is possible.

Marthi, Phani Ratna Vanamali [ORNL] (ORCID:0000000↗

Oscillation Risks of Grid-Following and Grid-Forming Inverter-Based Resources in Series-Compensated Networks

Here, this paper investigates the dynamic behavior of a grid-connected inverter-based resource (IBR) when connected radially to a series compensated line. Potential interactions between the series compensation and the IBR have been identified for both types: grid-following (GFL) or grid-forming (GFM). The study begins with electromagnetic transient (EMT) simulations to demonstrate stability issues. Subsequently, nonlinear analytical models are formulated in the dq frame, validated against the EMT simulation, and leveraged to assess eigenvalues and participation factors. Influencing factors of the dominant oscillation modes have been identified. The analysis results show that series compensation may make a mode associated with the synchronization unit unstable. Furthermore, customized feedback systems are built for the synchronizing loop. Series compensation can increase the sensitivity of the voltage phase angle towards the synchronizing angle, and introduce phase lag in the real power response towards the synchronizing angle. These factors may cause interactions with the phase-locked loop in GFL-IBR systems and with power-based synchronization in GFM-IBR systems, potentially leading to instability.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Distributed Small-Signal Stability Conditions for Inverter-Based Unbalanced Microgrids

The proliferation of inverter-based generation and advanced sensing, controls, and communication infrastructure have facilitated accelerated deployment of microgrids. A coordinated network of microgrids can maintain reliable power delivery to critical facilities during extreme events. Low-inertia offered by the power-electronics–interfaced energy resources, however, can present significant challenges to ensuring stable operation of the microgrids. In this work, distributed small-signal stability conditions for inverter-based microgrids are developed that involve the droop-controller parameters and the network parameters (e.g. line impedances, loads). The distributed closed-form parametric stability conditions derived in this paper can be verified in a computationally efficient manner, facilitating reliable design and operations of networks of microgrids. Dynamic phasor models have been used to capture the effects of electromagnetic transients. Furthermore, numerical results are presented, along with PSCAD simulations, to validate the analytical stability conditions. Effects of design choices, such as the conductor types, and inverter sizes, on the small-signal stability of inverter-based microgrids are investigated to derive useful engineering insights.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Xyce(™) Parallel Electronic Simulator v.7.5

The Xyce Parallel Electronic Simulator simulates electronic circuit behavior in DC, AC, HB, MPDE and transient mode using standard analog (DAE) and/or device (PDE) device models including several age and radiation aware devices. It supports a variety of computing platforms (both serial and parallel) computers. Lastly, it uses a variety of modern solution algorithms dynamic parallel load-balancing and iterative solvers.! ! Xyce is primarily used to simulate the voltage and current behavior of a circuit network (a network of electronic devices connected via a conductive network). As a tool, it is mainly used for the design and analysis of electronic circuits.! ! Kirchoff's conservation laws are enforced over a network using modified nodal analysis. This results in a set of differential algebraic equations (DAEs). The resulting nonlinear problem is solved iteratively using a fully coupled Newton method, which in turn results in a linear system that is solved by either a standard sparse-direct solver or iteratively using Trilinos linear solver packages, also developed at Sandia National Laboratories.

Source record↗

Transient superconductivity at nano- and meso-scales

High-T c superconductivity is surpassed by few, if any, other unsolved problems in contemporary physics in terms of its richness, complexity, impact on other fields, and potential technological importance in energy applications. Recent discoveries reveal that the highest transition temperatures emerge under extraordinary experimental conditions such as the strong perturbation of a single atomic layer material by the underlying substrate, application of ~200 GPa pressure or ~MV/cm static electric fields, or irradiation by intense femtosecond optical fields. Research in superconductivity under extreme experimental conditions is challenging and requires the development of novel experimental methods and new theoretical tools suitable to tackle the problem. Our team has approached the challenges of transient superconductivity under intense optical fields by rethinking the types of experimental observables that are best suited to elucidate the physics of light-induced effects. This exercise has identified the need for: i) ultra-fast spatio-temporal probes enabling access to plasmonic properties and also nano-scale inhomogeneities that are ubiquitous in unconventional superconductors; ii) meso-scale structures and hybrid meta-surfaces imperative for strong enhancement of optical fields and iii) theoretical approaches suitable to explain and predict the response of superconductors under ultra-fast photo-excitation. Co-PIs have already made major advances with developing capabilities i)-iii). Therefore, our team is poised to make significant progress in transient superconductivity by exploring novel spatio-temporal effects that currently remain unattainable. Co-investigators will search for light-induced superconductivity in the high-T c cuprates to build upon spectacular recent results still lacking independent confirmations. Novel data acquisition methods combined with nano-plasmonic imaging will allow our team to test the hypothesis of photo-induced superconducting pairing. Spatio-temporal experiments will provide insights into the interplay between superconductivity and competing orders in the cuprates. Co-PIs will also investigate strong light-matter interaction and superconductivity in FeSe monolayers. Finally, we will study spatio-temporal electrodynamics of planar Josephson junctions. This latter direction will allow our team to thoroughly characterize one of the most fundamental examples of inhomogeneity in all of unconventional superconductivity. The proposed experiments in combination with theoretical analysis will provide information that is difficult to obtain using alternative methods. Basov and Averitt will carry out pump-probe spectroscopy and nano-imaging experiments. Averitt and Hone will design and fabricate state-of-the-art meta-surfaces. Theoretical and computational studies of ultrafast transient phenomena will be carried out by Millis and Fogler. Modeling of time-resolved nano-optical effects will be done by Fogler. The bulk of the requested budget will be used to support graduate students at Columbia & UCSD that will be co-supervised by co-PIs. The proposed program is transformative, since it will enable an entire suite of experiments previously either impossible or technically implausible. It will deliver critically important insights not only into mechanisms of unconventional superconductivity but also to many other correlated quantum materials.

36 MATERIALS SCIENCE↗

Shallow-to-Deep Convective Transition in the Amazon (Final Report)

The frequent and extensive deep cloud formations over the Amazon Basin provide a vast area of atmospheric heating in the Tropics, forcing large-scale global circulations on daily to secular time scales. Motivated by the importance of these systems to the global climate system and the well-documented seasonal cycle in rainfall over the Amazon Basin, this study focused on improving understanding of the processes governing these systems on seasonal time scales and their impact on the typical growth of convection from early afternoon through evening, termed the shallow-to-deep transition. This study is based on observations from the Department of Energy (DOE) ARM Mobile Facility (AMF) near Manacapuru, Brazil collected during the Observations and Modeling of the Green Ocean Amazon (GOAmazon) field campaign from 1 January 2014 to 30 November 2015, together with total column water vapor (CWV) and surface meteorological observations from a dense Global Navigation Satellite System (GNSS) network within 50 km of the AMF site provided by the CHUVA Program collected in coordination with the GOAmazon campaign. The specific objectives of this study were to address the following questions: 1) What environmental conditions support the shallow-to-deep convective transition over the Central Amazon? 2) How does the relative importance of these conditions compare during the wet and dry seasons? 3) How are these conditions modulated by basin-scale transient forcing? 4) How well are these processes captured by regional model simulations convective resolving scale? Our main hypothesis was that intra-seasonal atmospheric disturbances were a primary driver of rainfall during the wet season in the Central Amazon, while localized surface heating of the land was the primary driver of enhanced rainfall during the dry season. Overall, our findings supported our hypothesis, with the convectively active phase of atmospheric Kelvin waves strongly favoring the development of larger, more organized rain systems during the wet season, while localized isolated deep convection dominated the types of rain events seen during the dry season. These results have implications for the predictability of rainfall over the Central Amazon and for its representation in climate and weather models.

54 ENVIRONMENTAL SCIENCES↗

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↗

Updated SAM Model for the Molten Salt Reactor Experiment (MSRE)

The development of reference standard problems based on prototypical reactor designs is of particular importance to verify the adequacy of computer codes and evaluation models for specific reactor types. To support the multiphysics coupled simulation of molten-salt-fueled reactor (MSR) using SAM and Griffin computer codes for safety and licensing analysis, much efforts have been put into enhancing code capabilities and developing reference models for the MSR primary loop in SAM. In this work, a previously developed Molten Salt Reactor Experiment (MSRE) primary loop model was updated to include a two-dimensional (2-D) core region and external core components in one-dimension (1-D) or zero-dimension (0-D). To ensure accurate feedback calculation in multi-physics simulations, the delayed neutron precursor tracking model and solid graphite model were added in the SAM model. In addition, the 2-D and 1-D domains are tightly coupled using the recently developed single-solve approach in SAM. The updated model has been tested under both steady-state and transient scenarios to demonstrate its potential for the multi-physics simulation of MSRE with coupled SAM and Griffin.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Model Specification of Droop-Controlled, Grid-Forming Inverters (REGFM_A1)

This document describes a positive-sequence model of droop-controlled, grid-forming (GFM) inverter-based resources (IBRs). It can be considered as an initial model for evaluating the impacts of GFM IBRs on the transients and dynamics of transmission systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Report on DIC Cyclic Dryout Testing of ATF Cladding Materials

Boiling water reactors can experience anticipated operational occurrences, where an expected deviation from normal operating conditions is witnessed for short time scales. One example of such transients arises from cladding dryout, where stable boiling becomes unstable, resulting in cladding heat up before re-wetting occurs. To understand the impact of such transients on fuel cladding lifetime and performance, cyclic dryout conditions were imposed upon fresh fuel cladding in the absence of irradiation. Pressurized cladding tube segments were exposed to rapid temperature increases to intermediate temperatures below a cladding rupture threshold and deformation was monitored in situ via 3-D digital image correlation. The impact of cycle frequency on cladding deformation was assessed. A clear understanding of elastic and plastic strain responses were developed for Zr alloy deformation under temperature transients, and these results can be used to support material model development and validate fuel performance codes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Study of the Protection Improvements for a Weak Grid Area With High Inverter-Based Resources (IBRs)

This project designs enhanced protection scheme for the real-world weak grid area with a high penetration of IBRs. As the existing protection schemes are originally designed for traditional synchronous machines, we first evaluate if the protection scheme will continue to operate reliably in systems with high levels of IBRs. Hardware relays are tested using a controller-hardware-in-the-loop setup. PSCAD electromagnetic transient simulation with IBR original equipment manufacturer black-box models is used to perform fault studies and generate COMTRADE data, which are replayed by a real-time digital simulator (RTDS) to feed input to the hardware relays. Three scenarios are analyzed: normal operation, an N-1 contingency, and an IBR-only scenario. The evaluation results reveal the following: 1) the protection scheme remains reliable under normal conditions and N-1 contingencies and 2) in IBR-only scenarios, differential protection (87L) continues to operate reliably, whereas local protection elements, such as distance and directional elements, fail because of the lack of regulated negative sequence current contributed by IBRs. Enhanced protection is designed to address the challenge of lack of negative sequence current from IBRs, including increased restraining factors a2 and k2 to block 32Q or using V instead QV ORDER for ground faults, enhanced mho distance element with voltage and phase angle supervision for L-L faults. The efficacy of enhanced protection logic is validated and proven to work reliably. Additionally, IEEE Std. 2800-2022 negative sequence current compliant GFL and GFM IBRs from another vendor are tested and proven to work reliably without need for enhanced logic. Therefore, this work provides valuable decision-making for utilities facing protection system challenges due to IBRs, either designing enhanced protection scheme or requesting their IBRs being IEEE Std. 2800-2022 compliant to produce regulated negative sequence current for protection relay to make correct decision.

24 POWER TRANSMISSION AND DISTRIBUTION↗