Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transients (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 739 records · Page 41

Modeling of Transient Flow Mixing of Streams Injected into a Mixing Chamber

Ignition is recognized as one the critical drivers in the reliability of multiple-start rocket engines. Residual combustion products from previous engine operation can condense on valves and related structures thereby creating difficulties for subsequent starting procedures. Alternative ignition methods that require fewer valves can mitigate the valve reliability problem, but require improved understanding of the spatial and temporal propellant distribution in the pre-ignition chamber. Current design tools based mainly on one-dimensional analysis and empirical models cannot predict local details of the injection and ignition processes. The goal of this work is to evaluate the capability of the modern computational fluid dynamics (CFD) tools in predicting the transient flow mixing in pre-ignition environment by comparing the results with the experimental data. This study is a part of a program to improve analytical methods and methodologies to analyze reliability and durability of combustion devices. In the present paper we describe a series of detailed computational simulations of the unsteady mixing events as the cold propellants are first introduced into the chamber as a first step in providing this necessary environmental description. The present computational modeling represents a complement to parallel experimental simulations' and includes comparisons with experimental results from that effort. A large number of rocket engine ignition studies has been previously reported. Here we limit our discussion to the work discussed in Refs. 2, 3 and 4 which is both similar to and different from the present approach. The similarities arise from the fact that both efforts involve detailed experimental/computational simulations of the ignition problem. The differences arise from the underlying philosophy of the two endeavors. The approach in Refs. 2 to 4 is a classical ignition study in which the focus is on the response of a propellant mixture to an ignition source, with emphasis on the level of energy needed for ignition and the ensuing flame propagation issues. Our focus in the present paper is on identifying the unsteady mixing processes that provide the propellant mixture in which the ignition source is to be placed. In particular, we wish to characterize the spatial and temporal mixture distribution with a view toward identifying preferred spatial and temporal locations for the ignition source. As such, the present work is limited to cold flow (pre-ignition) conditions

Voytovych, Dmytro M.↗

Design of Experiments for Dynamic Test Runs in Solvent-Based CO 2 Capture Pilot Plants

Test runs in the pilot plants consume significant resources, and therefore, the learning from test runs should be maximized. Test runs conducted in the pilot plants are often steady state. It takes several hours for reaching steady-state in the pilot plants, and thus, the duration of the test runs needs to be long even for collecting few steady-state data points. On the other hand, a large number of measurements can be collected through dynamic test runs in a short span of time. This paper presents a systematic design of dynamic experiments (DoDEs) for identifiability of model parameters, which is achieved by persistently exciting the inputs signals. A pseudorandom binary sequence (PRBS) is designed as the input signal for DoDE due to its efficiency in obtaining sufficient spectral content. However, due to the long sequence size of the PRBS signal, a Schroeder-phase input signal, which is a multisine signal, is also designed. Tests for both types of signals are run in the Pilot Solvent Test Unit (PSTU) at the National Carbon Capture Center in Wilsonville, Alabama. The transient data are used to solve dynamic data reconciliation and parameter estimation problem. The estimated parameters are found to be not only superior to those estimated from using data collected from hundreds of steady-state test runs in a nonreactive (air–water) system, but the parameters could be estimated by using the dynamic data collected for about 24 h from the pilot plant for the MEA-H 2 O–CO 2 system.

CO2 capture↗

Enhanced simulation software for rocket turbopump, turbulent, annular liquid seals

One of the main objectives of this work is to develop a new dynamic analysis for liquid annular seals with arbitrary profile and to analyze a general distorted interstage seal of the space shuttle main engine high pressure oxygen turbopump (SSME-ATD-HPOTP). The dynamic analysis developed is based on a method originally proposed by Nelson and Nguyen. A simpler scheme based on cubic splines is found to be computationally more efficient and has better convergence properties at higher eccentricities. The first order solution of the original analysis is modified by including a more exact solution that takes into account the variation of perturbed variables along the circumference. A new set of equations for dynamic analysis are derived based on this more general model. A unified solution procedure that is valid for both Moody's and Hirs' friction models is presented. Dynamic analysis is developed for three different models: constant properties, variable properties, and thermal effects with variable properties. Arbitrarily varying seal profiles in both axial and circumferential directions are considered. An example case of an elliptical seal with varying degrees of axial curvature is analyzed in detail. A case study based on predicted clearances of an interstage seal of the SSME-ATD-HPOTP is presented. Dynamic coefficients based on external specified load are introduced to analyze seals that support a preload. The other objective of this work is to study the effect of large rotor displacements of SSME-ATD-HPOTP on the dynamics of the annular seal and the resulting transient motion. One task is to identify the magnitude of motion of the rotor about the centered position and establish limits of effectiveness of using current linear models. This task is accomplished by solving the bulk flow model seal governing equations directly for transient seal forces for any given type of motion, including motion with large eccentricities. Based on the above study, an equivalence is established between linearized coefficients based transient motion and the same motion as predicted by the original governing equations. An innovative method is developed to model nonlinearities in an annular seal based on dynamic coefficients computed at various static eccentricities. This method is thoroughly tested for various types of transient motion using bulk flow model results as a benchmark.

Padavala, Satya↗

Alfven Waves in the Solar Wind, Magnetosheath, and Outer Magnetosphere

Alfven waves Propagating outward from the Sun are ubiquitous in the solar wind and play a major role in the solar wind-magnetosphere interaction. The passage of the waves generally occurs in the form of a series of discrete steepened discontinuities, each of which results in an abrupt change in the interplanetary magnetic field direction. Some orientations of the magnetic field permit particles energized at the Earth's bow shock to gain access to the foreshock region immediately upstream from the Earth's bow shock. The thermal pressure associated with these particles can greatly perturb solar wind plasma and magnetic field parameters shortly prior to their interaction with the Earth's bow shock and magnetosphere. The corresponding dynamic pressure variations batter the magnetosphere, driving magnetopause motion and transient compressions of the magnetospheric magnetic field. Alfven waves transmit information concerning the dynamic pressure variations applied to the magnetosphere to the ionosphere, where they generate the traveling convection vortices (TCVs) seen in high-latitude ground magnetograms. Finally, the sense of Alfvenic perturbations transmitted into the magnetosheath reverses across local noon because magnetosheath magnetic field lines drape against the magnetopause. The corresponding change in velocity perturbations must apply a weak torque to the Earth's magnetosphere.

Sibeck, D. G.↗

Tracking carrier and exciton dynamics in mixed-cation lead mixed-halide perovskite thin films

Mixed-cation lead mixed-halide perovskites simultaneously possess structural stability and high power conversion efficiency. A thorough study of both carrier and exciton dynamics is needed to understand the photophysical properties that underpin its superior photovoltaic performance. By utilizing a broadband transient absorption spectroscopy, we observe the carrier and exciton dynamics in a FA 0.85 Cs 0.15 Pb(I 0.97 Br 0.03 ) 3 (FCPIB) perovskite by simultaneously resolving the carrier and exciton contribution to the transient change of the absorption spectra, from which the carrier density and exciton oscillator strength can be determined. Our data reveal a quick and significant conversion of the photogenerated carriers to excitons, on top of the usual carrier recombination process. Moreover, the decay of carrier density shows a change of kinetics from a second-order recombination at high pump fluence to a third-order recombination at low pump fluence. Our analysis utilizes band anharmonicity, presents an independent determination of electronic temperature and quasi-Fermi energy, and reveals an interesting interplay among the processes of carrier cooling, exciton formation/decay and carrier recombination, all as a function of time after photoexcitation. Our work demonstrates the use of pump fluence as a knob to tune the relative populations of carriers and excitons in halide perovskite materials.

36 MATERIALS SCIENCE↗

Multiphysics Coupling of PROTEUS-NODAL and SAM for Molten Salt Reactor Simulation

The purpose of this report is to document the new features and capabilities that have been implemented in the PROTEUS-NODAL code and the coupling procedure of PROTEUS-NODAL with SAM to analyze MSRs. This report is organized as follows. Section 2 presents the variational nodal P 1 solver for cylindrical geometries and its verification test results. Section 3 discusses the transient analysis capability implemented in the P 1 solvers of PROTEUS-NODAL. The solution methods for the coupled system of TFSP and precursor drift equations, the coarse mesh finite difference (CMFD) acceleration, the kinetics parameter evaluation method, and the thermal feedback are discussed. The verification test results against the SP 3 solutions are also discussed. Section 4 describes the coupling scheme of PROTEUS-NODAL and SAM under the MOOSE framework for MSR applications. The developed MOOSE applications named “TreeFrog” and “TreeKangaroo” are discussed. Verification test results for the coupled system are also discussed by comparing the MSFR steady state solution with that obtained from a manually coupled PROTEUS-NODAL and ANSYS CFX computational fluid dynamics (CFD) calculations and the MSFR transient analysis results with the open literature solutions and the PROTEUS-NODAL solutions obtained with the standalone thermal-hydraulics solver. Section 5 provides conclusions and future works.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Developing a PSCAD Model of the Reduced 240-Bus WECC Test System: Preprint

Fast-acting controls of inverter-based resources (IBR) introduce unprecedented dynamics to bulk power systems and island power grids, especially the high-frequency dynamics ranging from $>5$ Hz to hundreds of Hz. These oscillatory stability problems cannot be captured by the phasor domain simulations, e.g., in PSS/E, Powerworld and TSAT, due to insufficient modeling. Instead, they should be studied by electromagnetic transient (EMT) simulations by using proper EMT dynamic models. This paper develops and validates the PSCAD model of the reduced 240-bus Western Electricity Coordinating Council (WECC) test system. The PSCAD model will be released to the public, which can be used by both academia and power industry practitioners for studying the fast dynamic problems without any Critical Energy/Electric Infrastructure Information (CEII).

240-bus WECC test system↗

Theory of Ion-Mediated Segmental Localization, Activated Structural Relaxation, and the Glass Transition in Polymerized Ionic Liquids

Here we combine polymer integral equation theory for structural correlations with a force-level theory for activated segmental relaxation to construct a microscopic theory for how ions modify the onset of segmental localization, the polymer alpha relaxation time, and glass transition temperature, T g , in dry polymerized ionic liquids (PolyILs). Our most important findings concerning the behavior of T g are (i) a weak dependence on anion–cation Coulomb attraction strength due to strong decoupling of ion and segmental dynamics, (ii) strong reduction with increasing mobile ion size due to both plasticization effects and weakening of ion-modified monomer caging, and (iii) increase with backbone rigidity and intrachain dynamic cooperativity. We quantitatively estimate T g variations as ion size and polymer static and dynamic rigidity are widely varied. The experimental data of 19 PolyILs are organized into a specific pattern consistent with the theoretical predictions with glass transition temperatures varying over a very wide range (~250 K) as observed experimentally. Based on the theoretical results, suggestions are made concerning the search for high conductivity Li- or Na-based PolyILs based on further lowering of T g via new polymer synthesis. Calculations of the alpha time as a function of packing fraction reveal the dynamic fragility increases for larger mobile ions, a more flexible backbone, and/or a smaller degree of intrachain dynamic cooperativity, as a consequence of the increasing importance of collective elasticity in determining the activation barrier. An analysis of transient segmental localization and a crossover to activated dynamics reveals similar, but not identical, trends as predicted at the laboratory vitrification timescale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tandem bulk oxygen diffusion and surface reactions in reducible metal oxides control redox cycle dynamics

The interplay between bulk oxygen diffusion and surface reactions in reducible metal oxides is key in heterogeneous catalysts, but direct measurements of oxygen mobility, transient kinetics, and in situ spectroscopies have been lacking. Here, we reveal complex dynamic behavior of ceria-zirconia by H 2 using transient kinetics via mass spectrometry and in situ Raman and near-ambient pressure x-ray photoelectron spectroscopies. Molecular dynamics simulations with a machine learning potential delineate competitive oxygen diffusion mechanisms, with an optimal mobility at intermediate reductions. We expose a compensation between vacancy availability and lattice distortion at intermediate to high reductions and Frenkel defects at low reductions, underscoring a potential deficiency of 16 O/ 18 O exchange experiments in deducing oxygen mobility. Vacancies in proximity require electron localization on Ce atoms further away. The continuous replenishment of surface oxygen results in a varying reduction rate, with H 2 dissociation being the rate-limiting step. Multiscale transient simulations, consistent with experiments, indicate catalysts of potentially spatially varying oxidation states. The approach is broadly applicable to reducible oxide materials.

36 MATERIALS SCIENCE↗

GeN-ROM—An OpenFOAM®-based multiphysics reduced-order modeling framework for the analysis of Molten Salt Reactors

This work presents a projection-based multiphysics Model Order Reduction (MOR) framework for the analysis of nuclear systems and its application to parametric simulations of Molten Salt Reactors (MSR). The framework, named GeN-ROM, is developed using OpenFOAM® and employs a Proper Orthogonal Decomposition aided Reduced-Basis technique (POD-RB). It can be used to reduce steady-state and transient multiphysics problems involving parametric fluid dynamics, heat exchange, and neutronics phenomena. For the treatment of structural elements in the hydraulic systems, a porous medium approach has been adopted. The reduction process is data-driven and snapshot information is extracted via POD to learn the solution manifold and to build global spatial basis functions. At the data collection phase, GeN-ROM makes use of the solvers available in GeN-Foam, a similarly OpenFOAM®-based multiphysics framework developed for the analysis of nuclear reactors. The global bases are used both to approximate the solution fields and to project the full-order equations onto lower-dimensional subspaces, thus considerably reducing the number of unknowns in a numerical system. This reduction leads to significant computational speedups, which is ideal for multi-query applications such as uncertainty quantification or design optimization. The developed tool has been tested using a 2D multiphysics model of the Molten Salt Fast Reactor (MSFR) with steady-state and transient scenarios, with speedups on the order of 10 – 10 5 .

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Analysis of the transient compressible vapor flow in heat pipe

The transient compressible one-dimensional vapor flow dynamics in a heat pipe is modeled. The numerical results are obtained by using the implicit non-iterative Beam-Warming finite difference method. The model is tested for simulated heat pipe vapor flow and actual flow in cylindrical heat pipes. A good comparison of the present transient results for the simulated heat pipe vapor flow with the previous results of a two-dimensional numerical model is achieved and the steady state results are in agreement with the existing experimental data. The transient behavior of the vapor flow under subsonic, sonic, and supersonic speeds and high mass flow rates are successfully predicted. The one-dimensional model also describes the vapor flow dynamics in cylindrical heat pipes at high temperatures.

Jang, Jong Hoon↗

Survey of Dynamic Mode Decomposition Methods

Dynamic mode decomposition (DMD) is a data-driven reduced order modeling (ROM) technique used for dynamic systems. The widely adopted algorithm was first introduced and demonstrated on fluid flow data by Schmid. In recent years, various other fields, such as nuclear engineering, have begun to adopt this method. For example, DMD has been used for estimating α-eigenvalues, as an ROM for pulsed neutron problems, for predicting isotopic composition in burnup calculations, as acceleration techniques for iterative methods, and in capturing dynamic behaviors in molten salt reactor transients. This report seeks to demonstrate the capabilities and limits of the standard DMD algorithm, and identify problem spaces where variants may be better suited. The primary variant this report considers is Multi-Resolution DMD (mrDMD). Because this serves as a survey, synthetically produced data is used in lieu of simulation results. The remainder of this report will go into detail on the DMD theory, outline the standard DMD and mrDMD algorithms, present test cases highlighting the applicability of each, and finally present a discussion on how to determine the best suited algorithm for a given problem. All calculations performed in this report are carried out using the open source DMD library, PyDMD.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Drop impact tester development for spent nuclear fuel vibration integrity study

The accumulated damage of SNF rod is defined as the integrated damage induced by the SNF assembly harmonic vibrations and SNF assembly & cask system contact interactions induced transient shocks, experienced by SNF rods during normal conditions of transportation (NCT). The CIRFT technology has provided means to resolve/understand the SNF rod harmonic vibration reliability issues associated with maximum axial clad tubing stress induced by SNF inertia related dynamic bending loading. However, the SNF rod integrity in association with the transient shock loads during NCT through dynamic contact interactions among SNF assembly cannot be investigated systematically through CIRFT technology alone. Thus, the proposed drop impact tester development become a viable tool for SNF system transient shock effect investigation. Furthermore, the transient shocks induced loading mode is normally perpendicular to the SNF rod axial orientation, i.e., in a transverse orientation to the SNF rod. To proceed with such loading mode investigation, a new device, “SNF Drop Impact Tester” that provides the transverse impact shock load onto a SNF rod through a free-drop projectile was developed in FY20 for SNF System Vibration Reliability Investigation; where the SNF system contact impact intensity can be calibrated accordingly with projectile weight and projectile head geometry. This newly developed research tool can further provide a detailed understanding about the effect of dynamic contact-interaction loadings; and combined harmonic vibration and transient shock loading modes on the fatigue damage evolution of the HBU SNF under NCT.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Analysis of the transient compressible vapor flow in heat pipes

The transient compressible one-dimensional vapor flow dynamics in a heat pipe is modeled. The numerical results are obtained by using the implicit non-iterative Beam-Warming finite difference method. The model is tested for simulated heat pipe vapor flow and actual vapor flow in cylindrical heat pipes. A good comparison of the present transient results for the simulated heat pipe vapor flow with the previous results of a two-dimensional numerical model is achieved and the steady state results are in agreement with the existing experimental data. The transient behavior of the vapor flow under subsonic, sonic, and supersonic speeds and high mass flow rates are successfully predicted. The one-dimensional model also describes the vapor flow dynamics in cylindrical heat pipes at high temperatures.

Jang, J. H.↗

Time Resolved Visualization of Liquid Jet Interaction With H2-Air Rotating Detonations Using MHz Rate Diesel PLIF

The characterization of the dynamic response of liquid jets to transient detonation wave passage is critical for optimization and modeling of liquid fueled rotating detonation combustors. In this work, a rotating detonation combustor (RDC) is operated on hydrogen and air to sustain stable detonation waves that acts as a detonation driver and interact in a one-way coupled manner with a single liquid fuel jet that propagates into the combustion chamber with cycle periods of ~ 250 μs. Diesel is utilized as a realistic fuel surrogate with higher aromatic compounds to enable fluorescence excitation, using the 355 nm thirdharmonic output of a burst-mode Nd:YAG laser, imaged at rates up to 1 MHz. By optimizing the technique to accommodate orders of magnitude variations in the fuel density throughout the injection process, the PLIF data enable quantitative measurements including the refill time, the relative recovery between liquid and gaseous jets and jet trajectory across various momentum flux ratio. As the passage of the detonation wave imparts significant changes in the momentum flux ratio, the qualitative liquid break-up process and spatial distribution varies significantly in time. As the injection system recovers ~ 70% of the cycle period and return to a quasi-steady position and allow comparisons with theoretical jet trajectories. These data, enabled by ultra-high-speed PLIF imaging, represent some of the first detailed measurements for quantifying the dynamic response and recovery of liquid jets exposed to periodic detonations in an operating RDC.

Propulsion↗

Time Resolved Visualization of A Liquid Jet in an RDE Using MHz Rate Diesel PLIF

The characterization of the dynamic response of liquid jets to transient detonation wave passage is critical for optimization and modeling of liquid fueled rotating detonation combustors. In this work, a rotating detonation combustor (RDC) is operated on hydrogen and air to sustain stable detonation waves that acts as a detonation driver and interact in a one-way coupled manner with a single liquid fuel jet that propagates into the combustion chamber with cycle periods of ~ 250 μs. Diesel is utilized as a realistic fuel surrogate with higher aromatic compounds to enable fluorescence excitation, using the 355 nm third harmonic output of a burst-mode Nd:YAG laser, imaged at rates up to 1 MHz. By optimizing the technique to accommodate orders of magnitude variations in the fuel density throughout the injection process, the PLIF data enable quantitative measurements including the refill time, the relative recovery between liquid and gaseous jets and jet trajectory across various momentum flux ratio. As the passage of the detonation wave imparts significant changes in the momentum flux ratio, the qualitative liquid break-up process and spatial distribution varies significantly in time. As the injection system recovers ~ 70% of the cycle period and return to a quasi-steady position and allow comparisons with theoretical jet trajectories. These data, enabled by ultra-high-speed PLIF imaging, represent some of the first detailed measurements for quantifying the dynamic response and recovery of liquid jets exposed to periodic detonations in an operating RDC.

Propulsion↗