Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transient 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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Development of Whole System Digital Twins for Advanced Reactors: Leveraging Graph Neural Networks and SAM Simulations

Here, in this work, we introduce a novel method to develop whole system digital twins (DTs) for advanced nuclear reactors. This method treats a complex reactor system as a heterogeneous graph: with the system components as different types of graph nodes and their physical interconnections as edges. Based on the heterogeneous graph, a graph neural network combining graph convolution and temporal node attention is developed as the DT, facilitating a comprehensive understanding of the system's dynamic behavior. By utilizing the System Analysis Module (SAM) code for simulating various operational transients, we develop a graph-based database that trains the DT. This DT is characterized by two primary functions: It can infer the entire system's status using sparse node information, and it can predict the progress of transients based on current and historical system information. Our approach is validated through case studies on the Experimental Breeder Reactor II (EBR-II) system and a generic Fluoride-salt-cooled High-temperature Reactor (gFHR), demonstrating the DT's accuracy in forecasting operational transients. The DT's rapid computation capabilities enhance its potential for supporting advanced reactor operations, offering benefits in intelligent simulation, autonomous control, and anomaly detection, paving the way for improved safety analysis and intelligent component health management for advanced reactor systems and reducing their operations and maintenance cost.

EBR-II↗
NASA NTRS1975

The sources of material comprising a mass ejection coronal transient

The origin of the material ejected during a white-light coronal transient on August 26 and 27, 1973, is investigated using simultaneous observations of a slowly ascending prominence and the more rapid accompanying coronal transient. Nearly identical simultaneous images of the prominence were obtained over a six-hour period in the H-alpha and He II (304 A) emission lines; contemporaneous Skylab coronograph observations over approximately 1.5 hours showed that the mass-ejection coronal transient rose above the ascending prominence. Based on analysis of these observations, it is concluded that: (1) the bulk of the ejected material originated in the lower corona, despite the lack of an observed depletion there; (2) the material in the transient was at coronal temperature and was visible in polarized radiation due to Thomson scattering of photospheric light by free electrons; and (3) the total event was far larger, more energetic, and longer lasting than would have been inferred from the prominence observations alone.

Hildner, E.↗
NASA NTRS2001

Stage Separation Performance Analysis Project

Stage separation process is an important phenomenon in multi-stage launch vehicle operation. The transient flowfield coupled with the multi-body systems is a challenging problem in design analysis. The thermodynamics environment with burning propellants during the upper-stage engine start in the separation processes adds to the complexity of the-entire system. Understanding the underlying flow physics and vehicle dynamics during stage separation is required in designing a multi-stage launch vehicle with good flight performance. A computational fluid dynamics model with the capability to coupling transient multi-body dynamics systems will be a useful tool for simulating the effects of transient flowfield, plume/jet heating and vehicle dynamics. A computational model using generalize mesh system will be used as the basis of this development. The multi-body dynamics system will be solved, by integrating a system of six-degree-of-freedom equations of motion with high accuracy. Multi-body mesh system and their interactions will be modeled using parallel computing algorithms. Adaptive mesh refinement method will also be employed to enhance solution accuracy in the transient process.

Chen, Yen-Sen↗
DOE OSTI2024

Performance Improvements of the Griffin Solvers in FY24

The Griffin code is a MOOSE-based reactor physics application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, we have made significant efforts to improve the performance of transport solver options and cross-section generation for the efficient use of Griffin in advanced reactor applications. For the HFEM-PN solver, the residual evaluations of HFEM kernels were optimized by utilizing the pre- computed averaged cross sections for individual elements. Numerical integration involving the evaluation of basis functions at quadrature points was bypassed by facilitating precomputed element mass matrices for response matrices. Red-black iterations were improved by introducing a new generalized minimum residual based solver. The memory usage of response matrix storage was significantly reduced by applying basis function rotations on interfaces and calculating volumetric odd-parity moments on the fly. Additionally, the adjoint flux and transient calculation capabilities of the HFEM-PN solver were successfully implemented and verified using the TWIGL benchmark problem. For the DFEM-SN solver, memory footprint and computation time were significantly reduced by not treating angular flux vectors as the MOOSE nonlinear system vectors. Specifically for IQS, scalar adjoint weighting was introduced to further eliminate angular adjoint flux storage in the MOOSE auxiliary system. It was demonstrated through the three-dimensional Advanced Burner Test Reactor core problem that the memory usage for transient calculations with the IQS method was reduced by over 7.5× compared to before the optimizations. For the self-shielding application programming interface, a new double-heterogeneity treatment method, named the Bell Function-Based Analytic Two-Region Slowing Down Method, was developed to efficiently flux-volume homogenize TRISO particles with the matrix. Additionally, optimizations were made to hyper- fine group (HFG) slowing down calculations by pretabulating collision probability coefficients and grouping isotopes, significantly reducing the computational time for calculating scattering sources per HFG. Lastly, the pin power reconstruction module was extended to account for temporal behavior in a microreactor analysis problem, specifically for a control drum transient. Verification tests for each of these improvements demonstrated significant performance enhancements and memory reduction.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1990

BATSE software for the analysis of the gamma ray burst spatial distribution

The Burst and Transient Source Experiment (BATSE) on the Gamma Ray Observatory (GRO) is designed to study astronomical gamma ray sources and to provide better positional, spectral, and time resolution about these objects than has previously been possible from one experiment. The procedure to be used in the analysis of the gamma ray burst spatial distribution is presented. Data is input from BATSE via the Gamma Ray Burst Catalog (listing individual burst positions, flux values, and associated errors) and the Sky Sensitivity Map (which summarizes observational selection effects in table format). A FORTRAN program generates Monte Carlo burst catalogs, which are models to be compared to the actual distribution. The Monte Carlo models are then filtered through the Sky Sensitivity Map so that they suffer from the same selection effects as the actual catalog data. Additionally, each burst position is converted into a probability distribution to mimic BATSE positional sensitivity. The Burst Catalog, Monte Carlo burst catalog, and Sky Sensitivity Map are then passed onto an IDL program that compares the catalogs for statistical significance. The Sky Sensitivity Map is used to estimate how often each sky area is observed above the minimum flux level in question. Each burst found in this sky area is then weighted according to the frequency with which this sky area is observed. The catalogs are then compared via tests of homogeneity (based on their radial distributions) and isotropy (based upon their angular distributions). The results of the statistical comparisons along with graphs and charts of the summaries, are output from the IDL program for study.

Hakkila, Jon↗
DOE OSTI2022

Transient Data Library of Solar Grid Integrated Distributed System

This submission contains an open-source library of transient events in distributed system with high solar PV. The library includes the collected data, related documents and scripts for loading the data. The data library is built for transient event detection and machine learning based analysis algorithm development. The data was collected via both field test and software simulation. The units for the data are included in the data file headers for each data series. A text editor or spreadsheet software, such as Excel, and Matlab is required to view the data.

algorithms↗
NASA NTRS2014

Timing Studies of X Persei and the Discovery of Its Transient Quasi-periodic Oscillation Feature

We present a timing analysis of X Persei (X Per) using observations made between 1998 and 2010 with the Proportional Counter Array (PCA) onboard the Rossi X-ray Timing Explorer (RXTE) and with the INTEGRAL Soft Gamma-Ray Imager (ISGRI). All pulse arrival times obtained from the RXTE-PCA observations are phase-connected and a timing solution is obtained using these arrival times. We update the long-term pulse frequency history of the source by measuring its pulse frequencies using RXTE-PCA and ISGRI data. From the RXTEPCA data, the relation between the frequency derivative and X-ray flux suggests accretion via the companion's stellar wind. However, the detection of a transient quasi-periodic oscillation feature, peaking at approximately 0.2 Hz, suggests the existence of an accretion disc. We find that doublebreak models fit the average power spectra well, which suggests that the source has at least two different accretion flow components dominating the overall flow. From the power spectrum of frequency derivatives, we measure a power-law index of approximately − 1, which implies that, on short time-scales, disc accretion dominates over noise, while on time-scales longer than the viscous time-scales, the noise dominates. From pulse profiles, we find a correlation between the pulse fraction and the count rate of the source.

accretion↗
NASA NTRS1999

Turbomachinery Design Quality Checks to Avoid Friction Induced Structural Failure

A unique configuration of the P&W SSME Alternate Fuel Turbopump turbine disk/blade assembly, combined with a severe thermal environment, resulted in several structural anomalies that were driven by frictional contact forces. Understanding the mechanics of these problems provides new quality checks for future turbo machinery designs. During development testing in 1997 of the SSME alternate fuel turbopump at Stennis Space Center, several potentially serious problems surfaced with the turbine disk/blade assembly that had not been experienced in extensive earlier testing. Changes to the operational thermal environment were noted based on analytical prediction of modifications that affected performance and on stationary thermal measurements adjacent to the rotor assembly. A detailed structural investigation was required to reveal the mechanism of distress induced by the change. The turbine disk experienced cracking in several locations due to increased thermal gradient induced stress during start and shutdown transients. This was easily predictable using standard analysis procedures and expected once the thermal environment was characterized. What was not expected was the curling of a piston ring used for blade axial retention in the disk, indentation of the axial face of the blade attachment by a spacer separating the first and second stage blades, and most significantly, galling and cracking of the blade root attachment that could have resulted in blade release. Past experience, in gas turbine environments, set a precedent of never relying on friction for help and to evaluate it only in specific instances where it was obvious that it would degrade capability. In each of the three cases above, friction proved to be a determining factor that pushed the components into an unsatisfactory mode of operation. The higher than expected temperatures and rapid thermal transients combined with friction to move beyond past experience. The turbine disk/blade assembly configuration contributed to the potential for these problems to occur by limiting the radial deflection from thermals and centrifugal loading. The cooled solid bore configuration was chosen to improve rotordynamic stability by limiting the length of rotor overhang while still protecting the roller bearing by maintaining zero slope under the inner race. During a start transient, the rim area of the disk heats rapidly and expands axially and circumferentially and requires corresponding radial and axial growth of the disk to maintain relative positioning of the disk, blades, spacers and retainer rings. The stiffness, large thermal mass, and bore cooling flow combine to severely limit the disk rim radial growth which results in the potential for relative movement between these parts. Friction then becomes a player in the determination of component stress.

Moore, Jerry H.↗
DOE OSTI2022

Nothing to see here: failed supernovae are faint or rare

ABSTRACT The absence of Type IIP core-collapse supernovae arising from progenitors above 17 solar masses suggests the existence of another evolutionary path by which massive stars end their lives. The direct collapse of a stellar core to a black hole without the production of a bright, explosive transient is expected to produce a long-lived, dim, red transient known as a failed supernova. Despite the detection of a number of candidates for disappearing massive stars in recent years, conclusive observational evidence for failed supernovae remains elusive. A custom-built pipeline designed for the detection of faint transients is used to re-analyse 10 yr of observations of 231 nearby galaxies from the PTF/ZTF surveys. This analysis recovers known supernovae, and yields a number of interesting transients. However, none of these are consistent with a failed supernova. Through Monte Carlo tests the recovery efficiency of our pipeline is quantified. By assuming failed supernovae occur as a Poissonian process with zero detections in the data set, 95 per cent upper limits to the rate of failed supernovae are calculated as a function of failed supernova absolute magnitude. We estimate failed supernovae to be less than 0.61, 0.33, 0.26, or 0.23 of the core-collapse SN rate for absolute magnitudes of −11, −12, −13, and −14, respectively. Finally, we show that if they exist, the Vera C. Rubin Observatory will find 1.7–3.7 failed SNe per year for an absolute bolometric luminosity of ∼6 × 1039 erg s−1 out to distances of 33–43 Mpc, depending on their assumed spectral energy distribution.

Byrne, R. A. (ORCID:0000000177661737)↗
DOE OSTI2024

Embedding Neural Thermal Scattering (NeTS) Modules in SERPENT for Higher Fidelity Advanced Reactor Analysis

When a neutron born in fission thermalizes to the order of $k$ $B$ $T$, it’s de-Broglie wavelength and energy approach the order of inter-atomic spacing and elementary lattice oscillations, respectively. $S$($a,β,t$) or the scattering law, uuantify these temperature-dependent crystallographic contributions to total cross section (or reaction rate). In a Monte Carlo analysis, cumulative distribution functions (CDFs) of $S$($a,β,t$) are loaded to memory from “A Compact ENDF” (ACE) files for stochastically selecting thermal scattered neutron trajectories. In this work, novel neural thermal scattering (NeTS) modules for $S$($a,β,t$) CDFs are designed, trained, serialized and embedded within SERPENT using Python’s limited C-API for on-the-fly deployment of crystalline graphite $S$($a,β,t$) sampling. Torchscript tracing and Numba just-in-time (JIT) compilation streamline neural inference on NVIDIA GPUs with CUDA libraries. Demonstrations of bare sphere thermalization of fast and thermal sources show excellent agreement between embedded NeTS in SERPENT and MCNP. With an explicit model of the reactor, NeTS can predict on-the-fly changes in TREAT neutron spectra as a function of local temperature, which can serve to improve transient and accident predictions in a multiphysics analysis framework. This framework can be further extended to account on-the-fly for changes in local graphitic microstructure to scattering cross sections, and outlines a novel coupling of modern machine learning with state-of-the-art reactor physics methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Galvanostatic Intermittent Titration Technique Reinvented: Part I. A Critical Review

The galvanostatic intermittent titration technique (GITT), introduced in 1977 by Weppner and Huggins, provided a readily accessible means to measuring the chemical diffusion coefficient of electrochemical electrode materials. The method continues to be widely used today, but the reported diffusivity values are highly inconsistent, ranging as much as four orders of magnitude for some Li layered oxide compositions. Even qualitative trends of diffusivity are inconsistent, suggesting significant flaws in the implementation of the method. Other variants of the GITT method also suffer from similar inconsistency problems. Here we identify numerous sources of significant error including composition-dependent reaction overpotentials, mathematical flaws in the relaxation analysis methods, finite-size and non-planar geometry effects, inter-particle inhomogeneity issues, early transient effects, and surface area uncertainties. We propose a simple relaxation analysis scheme using the time variable t relax + τ - t relax , where t relax is relaxation time and τ is the galvanostatic pulse duration. We also propose to use dense diffusion-limited samples to isolate the bulk-diffusion process in the time domain. Chemical diffusivity can be extracted much more reliably with this improved implementation of the GITT method.

25 ENERGY STORAGE↗
DOE OSTI2024

Advanced grid-forming (GFM) inverter controls, modeling and system impact study for inverter dominated grids

Our main goals in this project are to understand: 1) the stability of the power system with very high penetration of PV or other inverter-based resources (IBR), 2) the interactions among IBRs and traditional generators, and 3) the impact of the GFM inverters on bulk power system. To achieve these goals, we took a holistic research and development approach. First, an impedance-based large system modeling and stability analysis tool was developed to aid the large IBR system small-signal stability analysis. Second, consensus-based control and transient overload ride through were developed to control and coordinate multiple GFM inverters in a large system to maintain transient stability during big grid events. Third, high fidelity electromagnetic transient (EMT) IBR models were developed and integrated into a large reduced WECC model to study the system impact from IBRs. Lastly, extensive validation tests of the GFM control which was developed in this project and implemented in two 1MW commercial solar inverters were performed using different validation platforms, e.g., control hardware-in-the-loop (CHIL) platform, power hardware- in-the-loop (PHIL) platform, and a real PV plant.

14 SOLAR ENERGY↗
NASA NTRS2005

Unbalance Response Prediction for Accelerating Rotors With Load-Dependent Nonlinear Bearing Stiffness

Rolling-element bearing forces vary nonlinearly with bearing deflection. Thus an accurate rotordynamic analysis requires that bearing forces corresponding to the actual bearing deflection be utilized. Previous papers have explored the transient effect of suddenly applied imbalance and the steady-state unbalance response, using bearing forces calculated by the rolling-element bearing analysis code COBRA-AHS. The present work considers the acceleration of a rotor through one or more critical speeds. The rotordynamic analysis showed that for rapid acceleration rates the maximum response amplitude may be considerably less than predicted by steady-state analysis. Above the critical speed, transient vibration at the rotor natural frequency occurs, similar to that predicted for a Jeffcott rotor with constant-stiffness bearings. A moderate amount of damping will markedly reduce the vibration amplitude, but this damping is not inherent in ball bearings.

Fleming, David P.↗
DOE OSTI2022

Dicyanobenzothiadiazole (DCBT) Organic Dye as a Visible Light Absorbing Strong Photoinduced Oxidant with a 16 Microsecond Long–Lived Excited State

Strong photoinduced oxidants are important to organic synthesis and solar energy conversion, to chemical fuels or electric. For these applications, visible light absorption is important to solar energy conversion and long–lived excited states are needed to drive catalysis. With respect to these desirable qualities, a series of five 5,6–dicyano[2,1,3]benzothiadiazole (DCBT) dyes are examined as organic chromophores that can serve as strong photooxidants in catalytic systems. The series utilizes a DCBT core with aryl groups on the periphery with varying electron donation strengths relative to the core. The dyes are studied via both steady–state and transient absorption and emission studies. Additionally, computational analysis, voltammetry, crystallography, and absorption spectroelectrochemistry are also used to better understand the behavior of these dyes. Ultimately, a strong photooxidant is arrived at with an exceptionally long excited state lifetime for an organic chromophore of 16 µs. Finally, the long–lived excited state photosensitizer is well–suited for use in catalysis, and visible light driven photosensitized water oxidation is demonstrated using a water–soluble photosensitizer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1982

Thermal capabilities and graphical output of PAFEC

The program for automatic finite element calculations (PAFEC) developed for the computation of heat transfer in structures is discussed. The passive and interactive graphics capabilities of the PAFEC system are presented and future developments are outlined. This finite element system is shown to have significant thermal capabilities in support of its general structures. The program offers the following types of analysis: interactive graphics; user defined program control steps; steady state, transient heat transfer; boundary element methods; linear static, stress and displacements; modes and frequencies calculations; direct dynamic time intergration; frequency response analysis; elastohydrodynamic lubrication; large deflection analysis, buckling; creep and plasticity analysis; and substructures analysis.

Akin, J. E.↗
DOE OSTI2025

Transient Modeling and Simulation of a Generic Stable Salt Reactor

A SAM system-level model of a generic stable salt reactor has been developed to investigate thermal-hydraulic behavior and safety performance under steady and transient conditions. The model integrates information generated from a reactor physics analysis using PROTEUS and PERSENT, and a computation fluid dynamics (CFD) analysis using STAR-CCM+. A loose, iterative coupling scheme between PROTEUS and SAM is implemented to calculate the equilibrium power and temperature distributions in the steady-state critical core condition. The converged steady-state model is then used in PERSENT to calculate the four reactivity feedback temperature coefficients (Doppler, fuel density, coolant density, and core radial expansion) and kinetic parameters that are needed in SAM to model the temperature feedback effects in transient simulations. Within the fully enclosed liquid fuel pins, natural convection is the dominant heat transfer mechanism. The STAR-CCM+ model of the fuel pin considers conjugate heat transfer from the liquid fuel salt to the pin cladding and external reactor coolant. The CFD results of the axial and radial temperature profiles are used to empirically determine an effective fuel salt thermal conductivity in the SAM fuel pin model so that the temperatures predicted by the SAM model match as closely as possible the CFD results. In the central region of the fuel pin, the effective thermal conductivity is as high as similar to 60 times the physical fuel salt thermal conductivity. The whole-plant SAM model is then used to simulate an unprotected station blackout transient. The results of this simulation showed that the large negative fuel axial expansion reactivity feedback reduces fission power to similar to 2.4% nominal power. The core is cooled by natural circulation, which removes heat in the core to the emergency heat removal system, and ultimately, to the ambient. However, peak fuel salt and cladding temperatures can potentially reach as high as 1500 K, albeit briefly, if the shutdown mechanism fails to operate.

stable salt reactor; transient simulations; system↗
NASA NTRS2024

Data-Driven Method for Nonlinear Modal Interaction Identification in Fighter Aircraft

The purpose of this paper is to report a new data-driven analysis methodology for identifying nonlinear modal interactions and its application to aeroelastic wind-tunnel data from a generic fighter configuration, which reveals global modal interactions in this system caused by multiple resonance capture events initiated by local nonlinearities. It is shown that during supersonic testing at the Transonic Dynamics Tunnel, vibro-impacts at the embedded external fuel tank shakers on the inboard underwing stations of the generic fighter introduce strong nonlinearity that gives rise to nonlinear modal interactions. Specifically, time–frequency analysis of the acceleration data demonstrates that 6∶5 and 7∶6 transient resonance captures occur across the aircraft exclusively during this vibro-impact forcing window. Further analysis suggests that these resonance captures are dependent on Mach number. These results show experimentally that local vibro-impacts can cause global nonlinear effects in an aeroelastic system and that these effects can be identified using the new methodology herein. The data-driven analysis scheme employed shows applicability to full-aircraft flight test data for identifying nonlinear modal interactions without knowledge of the precise nonlinear contributors.

nonlinear dynamics↗