Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transient dynamic analysis”

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 163 records · Page 9

Development and Porting of Nuclear Reactor Computational Models for the NRIC Virtual Test Bed in FY23

The U.S. DOE’s National Reactor Innovation Center (NRIC) hosts an open-source website and associated GitHub repository called the Virtual Test Bed (VTB) on which computational models for advanced reactors are documented and shared with the reactor community. In FY23 under the NRIC program, computational models for nuclear reactor analysis were developed and contributed to the VTB Github repository during FY23 by contributors at Argonne National Laboratory. Other models developed outside of NRIC were also ported to the VTB as part of this work, serving a broader mission to make computational reactor analysis models more widely available to the reactor community. The model development activities focused on demonstration of gas-cooled microreactor models including computational fluid dynamics simulations of flow through an industry-inspired air jacket design, and simulation of multiphysics transients for a gas-cooled microreactor assembly. Additionally, development of a molten salt reactor (led at Idaho National Laboratory) was supported through Argonne’s expertise in multigroup cross sections generation. Models for these two reactor types were targeted due to their relevancy to NRIC’s Demonstration of Microreactor Experiments (DOME) and Laboratory for Operation and Testing in the U.S. (LOTUS) physical test beds, which are slated to host microreactor and molten salt reactor experiments. The model porting activities consisted of developing detailed documentation for several physics models which originate from the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. We include only contributions from Argonne National Laboratory in this report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Improved Cross Section Generation Capability of Griffin in FY22

The Griffin code is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of deterministic steady-state transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication as well as transient solvers for spatial dynamics with the improved quasi-static method. Griffin uses cross section data in the ISOXML format generated from external deterministic or Monte Carlo cross section generation codes. In recent years, the MC 2 -3 modules have been added to Griffin for fast reactor cross section generation, and the self-shielding application programming interface (SSAPI) was implemented in the ISOXML module for thermal reactor cross section generation. The on-the-fly slowing down method and double-heterogeneity treatment have been implemented to SSAPI and verified against particulate fuel-bearing graphite-moderated thermal reactor problems with high accuracy. This year, work has been focused on improving the cross section generation capability of ISOXML and streamlining the cross section generation procedures. In addition, the form function data were added to ISOXML in order to support the pin power reconstruction capability that was newly implemented in Griffin in this fiscal year. To facilitate the cross section generation using MC 2 -3 and SSAPI in Griffin, the cross section generation workflows have been set up for both fast and thermal spectrum reactors. The MOOSE action system tool was devised for fast spectrum problems, and the MOOSE stochastic tool was adopted to the branch calculation procedure for thermal spectrum problems. Meanwhile, to ensure the accuracy of group-constants, the thermal up-scattering kernel calculator accounting for resonance scattering was implemented in ISOXML, demonstrating the accurate computation of a Doppler-broadened scattering kernel of any Legendre order within a reasonable timescale. Other aspects of ISOXML, such as deletion solver and data, documentation, ISOXML file management, and the interface for Mixture, were improved as well.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Transient uncertainty quantification and Global Sensitivity Analysis of the open-source Molten Chloride Reactor Experiment (MCRE) using GP-PCA surrogate models

Uncertainties in the thermophysical properties of molten salts impact both the steady-state and transient behavior of Molten Salt Reactors (MSRs). In this work, we aim to quantify the influence of such uncertainties on the transient operation of the Molten Chloride Reactor Experiment (MCRE), utilizing the open-source specifications provided for this reactor. Seven representative transient scenarios are considered. For each scenario, we evaluate the impact of thermophysical property uncertainties on four key multiphysics model output variables of interest (VoIs): maximum power density, maximum fuel temperature, maximum reflector temperature, and average fuel velocity magnitude. In addition, we perform a Global Sensitivity Analysis (GSA) by computing Sobol’ indices for the uncertain input parameters to determine their contribution to the variability of each VoI. Conducting GSA is computationally intensive due to the large number of required evaluations of the high-fidelity multiphysics model. To mitigate this cost, we develop a surrogate modeling framework that combines Gaussian Process (GP) regression with Principal Component Analysis (PCA), enabling efficient sample generation for the GSA. Our results show that for energy-related VoIs, thermal conductivity is the dominant contributor to uncertainty. In contrast, for flow-related VoIs, density and dynamic viscosity are the primary sources of uncertainty. The specific heat of the fuel salt was found to play a secondary role in the transient analyses.

42 - ENGINEERING↗

A multi-dimensional dataset for two-phase instability in low pressure natural circulation based on direct transient local measurement

Thermal hydraulic instability can challenge the control and the safety of nuclear reactors. This experimental research captures and investigates the instability in a low-pressure (1–10 atmospheres) natural circulation loop with an adiabatic chimney installed above the heated section. Direct transient local measurement of two-phase flow is performed at multiple axial locations along the flow channel with thermocouples and conductivity probes traversed by customized high-precision mechanisms. Based on the periodicity of oscillations, a post-processing method is proposed to align the data captured non-simultaneously to reconstruct the instability behavior in one representative period. The averaging operation over multiple oscillation cycles is performed to reduce the random errors. Cyclic phenomena are clearly demonstrated by the post-processed data, and rapid void fraction change induced by the flashing is observed to be significant in the driving mechanism for low pressure oscillations. Parametric effects of inlet subcooling, heat flux, inlet flow restriction and system pressure are studied within the experiment range, and qualitative results are explained and compared with existing literature. The ratio between the oscillation period and the travelling time of the fluid particle through the test section is compared with the theoretical value for the density wave oscillation. The new dataset contains directly measured transient multi-dimensional two-phase data along the flow channel for natural circulation under density wave oscillations. The generated dataset is valuable for the validation of system analysis codes and CFD codes in the simulation of natural circulation phenomena and dynamic two-phase flow.

42 ENGINEERING↗

Modeling Distributed Energy Resources for Analyzing Distribution Systems with High Renewable Penetration: Preprint

Increasing levels of inverter-based distributed energy resources (IBDERs) impact the legacy overcurrent distribution protection systems. The fault current injections of IBDERs are limited to between 1-2 p.u. of the rated current. The combination of the limited fault current and varying system load makes it increasingly difficult to reliably set the overcurrent protection. To address this challenge, researchers are proposing adaptive overcurrent protection mechanisms that can adapt to changing system conditions. Accurate modeling of IBDERs is required to develop reliable protection schemes. This paper presents a photovoltaics-based IBDER model in OpenDSS for adaptive overcurrent protection. The paper presents the validation of the proposed IBDER model against a detailed model in electromagnetic transient programs. Protection analysis of the Electric Power Research Institute J1 feeder with the proposed IBDER model is also presented.

adaptive↗

Rate Controlling in Low-Temperature Standard NH 3 -SCR: Implications from Operando EPR Spectroscopy and Reaction Kinetics

Here, a series of seven Cu/SSZ-13 catalysts with Si/Al = 6.7 are used to elucidate key rate-controlling factors during low-temperature standard ammonia-selective catalytic reduction (NH 3 -SCR), via a combination of SCR kinetics and operando electron paramagnetic resonance (EPR) spectroscopy. Strong Cu-loading-dependent kinetics, with Cu atomic efficiency increasing nearly by an order of magnitude, is found when per chabazite cage occupancy for Cu ion increases from ~0.04 to ~0.3. This is due mainly to the release of intercage Cu transfer constraints that facilitates the redox chemistry, as evidenced from detailed Arrhenius analysis. Operando EPR spectroscopy studies reveal strong connectivity between Cu-ion dynamics and SCR kinetics, based on which it is concluded that under low-temperature steady-state SCR, kinetically most relevant Cu species are those with the highest intercage mobility. Transient binuclear Cu species are mechanistically relevant species, but their splitting and cohabitation are indispensable for low-temperature kinetics.

catalytic reactions↗

Pump-induced motion of an interface between competing orders

In this work, we study the motion of an interface separating two regions with different electronic orders following a short-duration pump that drives the system out of equilibrium. Using a generalized Ginzburg-Landau approach and assuming that the main effect of the nonequilibrium drive is to transiently heat the system we address the question of the direction of interface motion—in other words, which ordered region expands and which contracts after the pump. Our analysis includes the effects of differences in free-energy landscape and in order parameter dynamics and identifies circumstances in which the drive may act to increase the volume associated with the subdominant order, for example when the subdominant order has a second-order free-energy landscape while the dominant order has a first-order one.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dynamic Validation of CNN-Based Surrogate Models for Inverter-Based Resources in Open-Source Solvers

Traditionally, distribution system planning has focused on steady-state analyses, with limited consideration of dynamic behavior. However, as large or medium-scale inverter-based resources (IBRs), particularly grid-following (GFL) inverters in commercial or industry buildings, become more prevalent, understanding their dynamic impact is essential for grid planning and operation. This article presents an innovative deep-learning (DL)-approach using convolutional neural networks technique to model the GFL inverters. Developed from real grid-tied commercial IBR transient data, these dynamic DL models overcome proprietary constraints by requiring minimal knowledge of internal converter physics while maintaining high accuracy and flexibility. To demonstrate their applicability, the models were incorporated into GridLAB-D, an open-source, three-phase distribution analysis tool. This integration enables dynamic simulations of large-scale distribution networks with high IBR penetration stability analysis. Rigorous testing and validation, aligned with industry standards, confirmed the reliability and efficiency of this approach, paving the way for enhanced planning and operational assessments of modern power systems.

Deep-learning↗

Mammalian complex III heme dynamics studied with pump-probe spectroscopy and red light illuminations

The electronic or molecular mechanisms which initiate photobiomodulation (PBM) in cells are not yet fully understood. Porcine complex III (C-III) of the electron transport chain was characterized with transient absorption spectroscopy (TAS). Then we applied our recently developed continuous wave laser coupled TAS procedure (CW-TAS) to investigate the effect of red light irradiances on the heme dynamics of C-III in its c 1 reduced state. The time constants associated with the axial ligand photodissociation were found to be 3.3 ± 0.3 ps for the oxidized state and 4.9 ± 0.4 ps for the c 1 reduced state. The analysis of the CW-TAS procedure yielded no significant changes in the C-III heme dynamics. We rule out the possibility of 635 nm CW light at 4.7 mW/cm 2 inducing a PBM effect on the heme dynamic of C-III, specifically with the photodissociation of its axial ligand.

59 BASIC BIOLOGICAL SCIENCES↗

Load Following Analysis of the Holos-Quad 10MWe Micro-Reactor with Plant Dynamics Code

This report presents the load following analysis of the “Holos-Quad” 10MWe micro-reactor design developed by HolosGen LLC. This analysis is focused on the Holos-Quad micro-reactor ability to match the changing grid demand at 10%/min rate. The control mechanisms for the plant are identified, simulated, and compared. Based on the comparison results, control strategy for load following of the Holos-Quad micro-reactor is developed. The control strategy and load following capabilities are demonstrated in a full-range down-and-up reactor power transient from 100% to 0% and back to 100% load at 10%/min rate. All calculations are carried out with Argonne’s Plant Dynamics Code.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The Development of a Generalized Riser Flow Regime Map Based Upon Higher Moment and Chaotic Statistics Using Electrical Capacitance Volume Tomography (ECVT)

Dynamic analyses have been applied to the temporal signals from an Electro Capacitance Volume Tomography instrument located near mid-height on the riser of an industrial-scale cold-flow circulating fluidized bed to characterize gas-solids flow behavior in the riser. Twelve capacitance electrodes surround the cylindrical riser over a height of 1.3 m. The instrument used a neural network deconvolution algorithm to determine the spatially resolved solids fraction recorded at 52 Hz. Experiments were carried out over a range of gas and solids flows in the transport regime using a Geldart Group B bed material, high density polyethylene with mean particle size of 880 μm. The radial solids distribution was found to vary from one-time step to the next between profiles typical of laminar and turbulent flow. The duration of time spent in each of these flow profiles depended upon the operating regime – dilute, core-annular, or fast fluidized bed. The chaotic structure of the temporal data was characterized using the three conventional approaches: the first 4 moments from the distribution of signal in time, system memory parameters from the autocorrelation function and the Hurst exponent, and analysis of the correlationentropy and correlation dimension of the attractor. These signal analysis techniques were used to clearly distinguish differences between different transport operating regimes. Specifically, it was experimentally observed that a riser transitions from core annular flow profile to dilute and dense regimes via increasing the frequency of short term transients to either dilute or dense flow profiles, respectively. A regime map was generated based upon these dynamics using solids flux and gas velocity axes. Fast fluidized, core annular, and dilute each exhibited different degree of dynamic characteristics typical of fluid dominated or particle compromising behavior. It should be noted that the magnitude for the different statistics was in the same range regardless of the regime, it was the radial profile for the statistic that changed and subsequently identified that there was a change in the regime. Finally, a reduced regime map was developed consisting of plotting the gas velocity normalized by the upper transport velocity versus the solids flux normalized by the saturation carrying capacity. The use of this reduced plot allowed the data from widely different conditions to be plotted and compared on the same<p>graph. Note that in many instances, some of the statistics identified the operating point as being in one regime while others indicated that it was in another indicating a transition region between dilute or core annular regimes and between the core annular and fast fluidization regimes. This now provides a tool that can be used to optimize process performance, identify changes in operating states, or replicate process dynamics during process scaling or changing operating parameters. </p>

Breault, Ronald↗

Heat load measurements for the PIP-II pHB650 cryomodule

This study presents a brief overview of the 1st and 2nd phases and an in-depth analysis of the 3rd phase heat load testing performed on the pHB650 (prototype High Beta 650 MHz) cryomodule at PIP2IT (PIP-II Injector Test Facility), with a focus on both the results and the methodological advancements that have improved testing efficiency and accuracy. A key challenge identified in the testing campaign is the higher-than-expected heat loads observed in the first PIP-II (Proton Improvement Plan II) prototype cryomodules (pSSR1 and pHB650) tested at PIP2IT. Elevated heat loads are concerning given the fixed capacity of the PIP-II cryoplant that is currently being installed at Fermilab. However, understanding the sources of these elevated heat loads offers a critical opportunity to implement effective heat load mitigations on upcoming PIP-II cryomodules to stay within the available capacity of the PIP-II cryoplant. The study includes a summary of test results, descriptions of measurement procedures, and key observations on parameters directly and indirectly related to heat load measurements. Direct observations include measured heat loads and the effectiveness of JT heat exchanger under varying conditions, while indirect observation analyze factors such as the temperature distribution on the two-phase pipe and relief piping under varying conditions. Thermal acoustic oscillations (TAO) were identified during testing, which was mitigated by replacing the original G10 stem with a stainless steel stem equipped with wipers for the cryomodule cooldown valve. A major innovation during pHB650 Phase 3 testing was the development of an automated Python script to streamline data acquisition, analysis, and reporting of heat load results. This script automatically retrieved data from ACNET (Accelerator Control Network), performed heat load calculations, and generated detailed reports featuring plots and tables. This advancement significantly reduced manual labor and enhanced the thoroughness of data analysis compared to earlier campaigns. The heat load test reports were promptly uploaded to the electronic logbook shortly after each test, enabling rapid feedback and collaboration between the SRF and cryogenic teams. The heat load measurements included various components: HTTS (high-temperature thermal shield), LTTS (low-temperature thermal shield), 2K isothermal and non-isothermal heat loads. Results were recorded both within the cryomodule and between the bayonet can supply and return. Measurements were conducted under different operating conditions such as "standard", "linac", and "simulated dynamic". Additionally, HTTS and LTTS heat loads were calculated in real time, allowing for the tracking of thermal stability and identification of changes during testing, both in steady-state and transient conditions. The results of this testing campaign not only provide valuable insights into the performance of the pHB650 cryomodule but also highlight best practices and lessons learned that will inform future cryomodule testing at PIP2IT. These include adopting automated tools for data analysis, refining real-time measurement capabilities, and emphasizing detailed pre-test planning. The framework established in this campaign aims to set an improved standard for cryomodule testing and heat load reporting in future cryomodule test campaigns.

Porwisiak, D. [Fermilab; Wroclaw Tech. U.]↗

Power System Recovery from Momentary Cessation with Transient Stability Improvement

Power system dynamics will be significantly changed by integrating wind farms and solar photovoltaic plants into power systems. This study investigates the effect of momentary cessation of inverter-based resources (IBRs) on transient stability and provides a recovery strategy for bulk IBRs in power systems. The theoretical analysis was initially carried out on a one-machine infinite-bus system, demonstrating the IBR impact in a critical group. The analysis was then expanded to a multimachine power system with IBRs using the single-machine equivalent method. The study found that IBRs in critical and noncritical groups exert contrasting effects on transient stability. Finally, a strategy for enhancing transient stability is proposed by controlling IBRs during power system recovery. The proposed strategy was verified by simulation on IEEE 9-bus and IEEE 39-bus power systems with the addition of IBRs. A 39-bus power system simulation demonstrates the scalability of the proposed method. Here, the proposed strategy provides effective and executable measures for improving system security in the presence of IBRs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Griffin Software Development Plan

Griffin is a MOOSE-based reactor physics application for advanced reactor multiphysics modeling and simulation. The application is developed in a consistent multiphysics environment with strong software quality assurance. Griffin inherited most of the capabilities of MAMMOTH/Rattlesnake and is adopting the capabilities from PROTEUS that are needed in the code. The toolset includes a variety of deterministic radiation transport solvers for fixed source, k-eigenvalue, ad-joint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. The code contains the cross-section preparation capabilities applicable to fast and thermal reactors, including TRISO-fueled reactors. Core management capabilities include core performance, fuel depletion and shuffling, equilibrium core calculation, pebble-bed reactor run-in and equilibrium core, molten-salt reactor delayed neutron precursor drift, and control rod and drum movement with cusping correction. This software development plan presents the current and future capabilities and features in Griffin for the design and analysis of non-light-water reactor systems in steady-state and transient conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Ultrashort XUV pulse absorption spectroscopy of partially oxidized cobalt nanoparticles

High-order harmonic generation (HHG) based transient extreme ultraviolet (XUV) absorption spectroscopy is an emerging technique to trace photoinduced charge carrier dynamics in condensed phase materials with femtosecond and even attosecond temporal resolution and elemental specificity. However, its application to nanoparticulate samples that are relevant, for example, for novel photocatalytic light harvesting concepts, has been limited. This is in part due to the challenge to produce residual-free samples on ultrathin, XUV-transparent substrates as well as a widespread understanding that sparsely distributed nanoparticles do not provide sufficient contrast for XUV absorption measurements. Here, we present static XUV absorption spectra of partially oxidized Co nanowire-structures with diameters of approximately 4.5 nm and lengths between 10 and 40 nm, recorded with an ultrashort pulse HHG light source. Nanoparticles are synthesized by the agglomeration of Co atoms inside superfluid helium droplets, followed by surface deposition and oxidation in ambient air. The method is uniquely suited for residual-free synthesis of transition metal nanowires and their deposition on ultrathin substrates. Analysis by high-resolution transmission electron microscopy reveals the formation of CoO nanowires with regions of unoxidized Co in their interior. The nanoparticle samples are investigated in an HHG-driven ultrafast XUV absorption setup. Despite the low surface coverage of only 23%, the recorded spectrum exhibits a distinct absorption feature at the Co M 2,3 (2p) edge near 60 eV with a peak height of about 40 mOD. Finally, the results support the feasibility of table-top ultrafast transient XUV absorption studies of photoinduced dynamics in transition metal oxide nanoparticles with sub-monolayer surface coverage.

74 ATOMIC AND MOLECULAR PHYSICS↗

Thermal-striping analysis methodology for sodium-cooled reactor design

Computational fluid dynamics simulations have been performed to study the applicability of engineering CFD methods for thermals striping analysis. Thermal striping is the fluctuating temperature profile in a solid caused by fluctuating fluid temperature, and the resulting fluctuating thermal stresses cause high cycle fatigue and eventual material failure. Here this study presents the methodology for thermal striping analysis that including a transient conjugate heat transfer RANS model and benchmark the results against a sodium triple jet experiment. The multiphysics calculations are implemented to analyze the thermal stresses in the solid domain caused by the coupled heat transfer between the fluid and solid domains. Validation data from a liquid sodium triple jet experiment include time-averaged temperature measurements and power spectra of the temperature signal. The numerical results agree well with these experimental measures, demonstrating key features such as the dominant frequency of temperature fluctuations. The applicability of a low-cost wall treatment method is demonstrated, enabling key computational savings. Finally, the performance of two finite element stress analysis software packages is compared, and the validity of the lower-cost method is confirmed. These results demonstrate the applicability of engineering methods for computational thermal striping calculations, enabling thermal striping estimations in large fluid systems such as the core of a Sodium-cooled nuclear reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗