Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transient Simulation”

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,765 records · Page 98

A Eulerian three-fluid flow framework for the study of fuel dispersal behavior under loss-of-coolant accident conditions

As the nuclear industry explores new fuel designs to accommodate increased burnup, studying fuel behavior during loss-of-coolant accidents is essential for ensuring the safe operation of light-water reactors. Here, this study presents a computational model to simulate the intricate three-phase flow dynamics of fuel dispersal following cladding breach, where solid fuel particles and carrier fission gases interact with the surrounding liquid or gas medium. The model utilizes a Eulerian three-fluid framework to depict the phase interaction within interpenetrating continua, treating fragmented fuel as granular material. Closure relationships for bulk-flow properties of the solid phase and detailed interfacial exchange terms in high-temperature, high-speed gas-solid-liquid flow scenarios are proposed. The model was validated using adiabatic experimental results in the literature, focusing on characterizing fuel dispersal behavior during loss-of-coolant accidents conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched probability density functions reported in experimental studies. Key highlights of this study include the theoretical description of a coupled three-phase flow with phase change, interpretation of the high-pressure boundary conditions and insights into the transient behavior of gas-solid-liquid phase dynamics during a simulated high-speed dispersal event.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Turbofan compressor dynamics during afterburner transients

The effects of afterburner light-off and shut-down transients on the compressor stability are investigated. The reported experimental results are based on detailed high-response pressure and temperature measurements on the TF30-P-3 turbofan engine. The tests were performed in an altitude test chamber simulating high-altitude engine operation. It is shown that during both types of transients, flow breaks down in the forward part of the fan-bypass duct. At a sufficiently low engine inlet pressure this resulted in a compressor stall. Complete flow breakdown within the compressor was preceded by a rotating stall. At some locations in the compressor, rotating stall cells initially extended only through part of the blade span. For the shut-down transient the time between first and last detected occurrence of rotating stall is related to the flow Reynolds number. An attempt was made to deduce the number and speed of propagation of rotating stall cells.

Kurkov, A. P.↗

Carbon Dioxide Absorption Modeling for Off-Gas Treatment in the Nuclear Fuel Cycle

The absorption of carbon dioxide is an important process in many practical applications such as reduction of greenhouse gases, separation and purification processes in the chemical and petroleum industries, and capture of radioactive isotopes in the nuclear fuel cycle The goal of this research is to develop a dynamic model to simulate CO 2 absorption by using different alkanolamines as absorption solvents. The model is based upon transient mass and energy balances for the chemical species commonly present in CO 2 gas-liquid absorption. A computer code has been written to implement the proposed model. Simulation results are discussed. The reported model simulates well the response to dynamic changes in input conditions. The proposed model can be used to optimize and control the separation of carbon-14 in the form of CO 2 in the nuclear industry.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Status of Multiple Channel Fuel Performance Capabilities Within the SAS4A/SASSYS-1 Safety Analysis Software

SAS4A/SASSYS-1 (SAS) is a fast-running simulation tool used to perform deterministic analysis of anticipated events as well as design basis and beyond design basis accidents for advanced liquid-metal-cooled nuclear reactors. It is a critical element of safety analysis capabilities for the U.S. Department of Energy and is utilized within industry to perform the transient safety analyses required to support the licensing of Liquid Metal-cooled Fast Reactors (LMFRs). Although SAS is exceptionally fast for most transient scenarios, fuel performance calculations, along with the associated pre-transient characterization of the fuel pin, may be required for transient scenarios where fuel pin failure is hypothesized. Both the pre-transient characterization and the transient fuel performance calculation are necessary to properly quantify margins to potential fuel failure and assess the time spent potentially exceeding such margins during events. While safety analysis calculations with fuel performance models provide a more detailed characterization of the reactor during a transient, the pre-transient characterization can be time-consuming and computationally expensive. Often, large numbers of fuel pins have been exposed to similar pre-transient irradiation conditions. Similarly, the same pre-transient fuel characterization may be applicable to numerous transient conditions. This provides an opportunity to optimize the SAS computational framework such that pre-transient fuel characterization can be shared across multiple channels (fuel pins) and across multiple simulations, thus dramatically reducing overall computational costs. This report summarizes progress toward enhancing the SAS computational framework to support shared, multiple channel fuel performance characterizations intended to significantly reduce computational costs. Preliminary testing has shown that the computational time saved by using the pre-transient sharing capability is approximately equal to the time it takes to perform the pre-transient characterization.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Real-time Hardware-in-the-Loop Evaluation of a Partially Turboelectric Propulsion Control Design

In support of aviation fuel burn and emission reduction goals, NASA is pursuing high-payoff research investments that promise to transform aviation. This includes investments in Electrified Aircraft Propulsion (EAP). Multiple technology challenges must be addressed to unlock the full potential of EAP. This includes addressing challenges related to propulsion controls, which will be vital for ensuring efficient coordinated operation of EAP subsystems. This paper presents results from real-time hardware-in-the-loop testing of a control design for a single aisle partially-turboelectric aircraft propulsion concept conducted at the NASA Electric Aircraft Testbed (NEAT) facility. The control system under test is designed for a propulsion concept consisting of two wing-mounted turbofan engines that produce thrust and generate electrical power to drive a boundary layer ingesting tailfan propulsor via an electrical motor. An integrated control strategy is applied to ensure coordinated operation of the turbofan and tailfan subsystems during steady-state and transient operation throughout the flight envelope. The NEAT test of this integrated control design consists of a partially hardware-in-the-loop, partially simulated configuration. A subscale representation of the electrical system design is implemented in hardware and mechanically coupled to electric machines that emulate turbomachinery and propulsor shaft dynamics. The hardware configuration is then operated under the control of a real-time computer application that runs a simulation of the propulsion system and the developed control logic. The NEAT facility test campaign includes a series of experiments that subject the control design to throttle transients conducted throughout the flight envelope and full-flight mission profiles. Testing under simulated performance degradation is also conducted to evaluate control design robustness. This includes constant and abrupt changes in degradation levels. Results from the hardware-in-the-loop test are presented and shown to be in good agreement with pre-test simulation predictions demonstrating the efficacy of the integrated control design approach.

Electrified Aircraft Propulsion↗

Real-time Hardware-in-the-Loop Evaluation of a Partially Turboelectric Propulsion Control Design

In support of aviation fuel burn and emission reduction goals, NASA is pursuing high-payoff research investments that promise to transform aviation. This includes investments in Electrified Aircraft Propulsion (EAP). Multiple technology challenges must be addressed to unlock the full potential of EAP. This includes addressing challenges related to propulsion controls, which will be vital for ensuring efficient coordinated operation of EAP subsystems. This paper presents results from real-time hardware-in-the-loop testing of a control design for a single aisle partially-turboelectric aircraft propulsion concept conducted at the NASA Electric Aircraft Testbed (NEAT) facility. The control system under test is designed for a propulsion concept consisting of two wing-mounted turbofan engines that produce thrust and generate electrical power to drive a boundary layer ingesting tailfan propulsor via an electrical motor. An integrated control strategy is applied to ensure coordinated operation of the turbofan and tailfan subsystems during steady-state and transient operation throughout the flight envelope. The NEAT test of this integrated control design consists of a partially hardware-in-the-loop, partially simulated configuration. A subscale representation of the electrical system design is implemented in hardware and mechanically coupled to electric machines that emulate turbomachinery and propulsor shaft dynamics. The hardware configuration is then operated under the control of a real-time computer application that runs a simulation of the propulsion system and the developed control logic. The NEAT facility test campaign includes a series of experiments that subject the control design to throttle transients conducted throughout the flight envelope and full-flight mission profiles. Testing under simulated performance degradation is also conducted to evaluate control design robustness. This includes constant and abrupt changes in degradation levels. Results from the hardware-in-the-loop test are presented and shown to be in good agreement with pre-test simulation predictions demonstrating the efficacy of the integrated control design approach.

Electrified Aircraft Propulsion↗

H Infinity Control of Magnetic Bearings to Ensure Both System and External Periodic Disturbance Robustness

Both self-excited and forced disturbances often lead to severe rotor vibrations in a magnetic bearing systems with long slender shafts. This problem has been studied using the H-infinity method, and stability with good robustness can be achieved for the linearized model of a magnetic bearing when small transient disturbances are applied. In this paper, the H-infinity control method for self-excited and forced disturbances is first reviewed. It is then applied to the control of a magnetic bearing rotor system. In modelling the system, the shaft is first discretized into 18 finite elements and then three levels of condensation are applied. This leads to a system with three masses and three compliant elements which can be described by six state variable coordinates. Simulation of the resultant system design has been performed at speeds up to 10,000 rpm. Disturbances in terms of different initial displacements, initial impulses, and external periodic inputs have been imposed. The simulation results show that good stability can be achieved under these different transient disturbances using the proposed controller while at the same time reducing the sensitivity to external periodic disturbances.

Jiang, Yuhong↗

Data trustworthiness signatures for nuclear reactor dynamics simulation

With the increased reliance on digitization in industrial control systems, the need for effective monitoring techniques has risen dramatically. Specifically, there is now a growing concern about the so-called false data injection (FDI) attacks. These attacks aim to alter the raw sensors’ data to cause malicious outcomes. Any serious FDI algorithm is based on an intimate knowledge of the system and its associated physics models, which renders conventional outlier/anomaly detection techniques almost obsolete in the face of such attacks. Thus, a critical need has emerged to develop a new class of defense methods that are capable of detecting FDI attacks under the assumption that the attacker has a strong familiarity with the system and its physics modeling. This class of defense methods are denoted by model-based defenses which are premised on the assumption that the attacker, while having a good understanding of the system, does not have full privileged access to all proprietary data and historical records of operation. However, (s)he is assumed to be capable of learning system behavior using self-learning techniques during an initial lie-in-wait period. To defend against this scenario, we propose a new model-based randomized window algorithm that searches time-series data for signatures that can serve as classifiers between normal and FDI scenarios. The classifiers are based on the correlations between the dominant degrees of freedom (DOFs) and the less-dominant DOFs (expected to be very sensitive to the system details that are unknown to the attacker). For demonstration, RELAP5 models are employed to calculate representative nuclear reactor behavior during a number of transient scenarios. Finally, falsified data are injected into the RELAP5-simulated behavior, and the proposed signature-identification algorithm is employed to detect the injected data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Large-volume flux closure during plasmoid-mediated reconnection in coaxial helicity injection

A large-volume flux closure during transient coaxial helicity injection (CHI) in NSTX-U is demonstrated through resistive magnetohydrodynamics (MHD) simulations. Several major improvements, including the improved positioning of the divertor poloidal field coils, are projected to improve the CHI start-up phase in NSTX-U. Simulations in the NSTX-U configuration with constant in time coil currents show that with strong flux shaping the injected open field lines (injector flux) rapidly reconnect and form large volume of closed flux surfaces. This is achieved by driving parallel current in the injector flux coil and oppositely directed currents in the flux shaping coils to form a narrow injector flux footprint and push the injector flux into the vessel. As the helicity and plasma are injected into the device, the oppositely directed field lines in the injector region are forced to reconnect through a local Sweet–Parker type reconnection, or to spontaneously reconnect when the elongated current sheet becomes MHD unstable to form plasmoids. In these simulations for the first time, it is found that the closed flux is over 70% of the initial injector flux used to initiate the discharge. These results could work well for the application of transient CHI in devices that employ super conducting coils to generate and sustain the plasma equilibrium.

Current Drive↗

Water Impact Test and Simulation of a Composite Energy Absorbing Fuselage Section

In March 2002, a 25-ft/s vertical drop test of a composite fuselage section was conducted onto water. The purpose of the test was to obtain experimental data characterizing the structural response of the fuselage section during water impact for comparison with two previous drop tests that were performed onto a rigid surface and soft soil. For the drop test, the fuselage section was configured with ten 100-lb. lead masses, five per side, that were attached to seat rails mounted to the floor. The fuselage section was raised to a height of 10-ft. and dropped vertically into a 15-ft. diameter pool filled to a depth of 3.5-ft. with water. Approximately 70 channels of data were collected during the drop test at a 10-kHz sampling rate. The test data were used to validate crash simulations of the water impact that were developed using the nonlinear, explicit transient dynamic codes, MSC.Dytran and LS-DYNA. The fuselage structure was modeled using shell and solid elements with a Lagrangian mesh, and the water was modeled with both Eulerian and Lagrangian techniques. The fluid-structure interactions were executed using the fast general coupling in MSC.Dytran and the Arbitrary Lagrange-Euler (ALE) coupling in LS-DYNA. Additionally, the smooth particle hydrodynamics (SPH) meshless Lagrangian technique was used in LS-DYNA to represent the fluid. The simulation results were correlated with the test data to validate the modeling approach. Additional simulation studies were performed to determine how changes in mesh density, mesh uniformity, fluid viscosity, and failure strain influence the test-analysis correlation.

Fasanella, Edwin L.↗

Plant Physiology Increases the Magnitude and Spread of the Transient Climate Response to CO 2 in CMIP6 Earth System Models

Increasing concentrations of CO 2 in the atmosphere influence climate both through CO 2 ’s role as a greenhouse gas and through its impact on plants. Plants respond to atmospheric CO 2 concentrations in several ways that can alter surface energy and water fluxes and thus surface climate, including changes in stomatal conductance, water use, and canopy leaf area. These plant physiological responses are already embedded in most Earth system models, and a robust literature demonstrates that they can affect global-scale temperature. However, the physiological contribution to transient warming has yet to be assessed systematically in Earth system models. Here this gap is addressed using carbon cycle simulations from phases 5 and 6 of the Coupled Model Intercomparison Project (CMIP) to isolate the radiative and physiological contributions to the transient climate response (TCR), which is defined as the change in globally averaged near-surface air temperature during the 20-yr window centered on the time of CO 2 doubling relative to preindustrial CO 2 concentrations. In CMIP6 models, the physiological effect contributes 0.12°C ( σ : 0.09°C; range: 0.02°–0.29°C) of warming to the TCR, corresponding to 6.1% of the full TCR ( σ : 3.8%; range: 1.4%–13.9%). Moreover, variation in the physiological contribution to the TCR across models contributes disproportionately more to the intermodel spread of TCR estimates than it does to the mean. The largest contribution of plant physiology to CO 2 -forced warming—and the intermodel spread in warming—occurs over land, especially in forested regions.

54 ENVIRONMENTAL SCIENCES↗

Improved Mars Helicopter Aerodynamic Rotor Model for Comprehensive Analyses

The Mars Helicopter is part of the NASA Mars 2020 rover mission scheduled to launch in July of 2020. Its goal is to demonstrate the viability and potential of heavier-than-air vehicles in the Martian atmosphere. Ultimately, it aims to bridge the resolution gap between orbiters and the rover as well as allow access to otherwise inaccessible regions. The low density of the Martian atmosphere and the relatively small-scale rotor result in very low Reynolds number flows. The low density and low Reynolds numbers reduce the lifting force and lifting efficiency, respectively. This paper describes the generation of the improved Mars Helicopter aerodynamic rotor model. The goal is to generate a performance model for the Mars Helicopter rotor using a free wake analysis, since this has a low computational cost for design. The improvements in the analysis are two-fold and are expanded on from two prior publications. First, the fidelity of the simulations is increased by performing higher-order two-dimensional time-accurate OVERFLOW simulations allowing for higher accuracy aerodynamic coefficients and a better understanding of the boundary layer behavior as well as its transient features. Second, a version of the model is generated to duplicate the exact testing conditions in the 25-ft. diameter Space Simulator at the Jet Propulsion Laboratory, which allows for better correlation of rotor performance figures. Previous work correlated performance with that test, but did not consider the higher temperatures in the experiment compared to those of the Martian atmosphere. The higher temperatures in the experiment are expected to give conservative performance estimates, as they give rise to an increase in speed of sound and decrease in observed Reynolds numbers.

Analyses↗

Improved Mars Helicopter Aerodynamic Rotor Model for Comprehensive Analyses

The Mars Helicopter is part of the NASA Mars 2020 rover mission scheduled to launch in July of 2020. Its goal is to demonstrate the viability and potential of heavier-than-air vehicles in the Martian atmosphere. Ultimately, it aims to bridge the resolution gap between orbiters and the rover as well as allow access to otherwise inaccessible regions. The low density of the Martian atmosphere and the relatively small-scale rotor result in very low Reynolds number flows. The low density and low Reynolds numbers reduce the lifting force and lifting efficiency, respectively. This paper describes the generation of the improved Mars Helicopter aerodynamic rotor model. The goal is to generate a performance model for the Mars Helicopter rotor using a free wake analysis, since this has a low computational cost for design. The improvements in the analysis are two-fold and are expanded on from two prior publications. First, the fidelity of the simulations is increased by performing higher-order two-dimensional time-accurate OVERFLOW simulations allowing for higher accuracy aerodynamic coefficients and a better understanding of the boundary layer behavior as well as its transient features. Second, a version of the model is generated to duplicate the exact testing conditions in the 25–ft. diameter Space Simulator at the Jet Propulsion Laboratory, which allows for better correlation of rotor performance figures. Previous work correlated performance with that test, but did not consider the higher temperatures in the experiment compared to those of the Martian atmosphere. The higher temperatures in the experiment are expected to give conservative performance estimates, as they give rise to an increase in speed of sound and decrease in observed Reynolds numbers.

Grip, Håvard F.↗

DNB evaluation of PWR steam-line break event using Westinghouse and CASL coupled code systems

ANCKVIPRE and the VERA-CS code systems have been successfully applied to predicting PWR core response with respect to departure from nucleate boiling (DNB) at the limiting time step of a postulated hot zero power (HZP) main steam-line break (MSLB) event with and without offsite power available. ANCKVIPRE links the ANC code with the VIPRE-W code. ANC is Westinghouse's multi-dimensional nodal diffusion code for all nuclear core design calculations. It predicts core reactivity, assembly power, rod power, detector thimble flux, and other relevant core characteristics. VIPRE-W is the Westinghouse version of a thermalhydraulic subchannel code developed for light water reactor core design applications. The reactor core simulation was based on input of the core boundary conditions from RETRAN system transient and STAR-CCM+ CFD core inlet distribution calculations. To account for the asymmetric power distribution in the reactor core due to the broken steam pipe in one loop and the stuck RCCA, both ANCKVIPRE and the VERA-CS models were set up for the whole reactor core. The advanced modeling and simulations confirmed that the HZP MSLB high-flow case with offsite power available is more limiting with respect to the DNB acceptance criterion, consistent with the limiting case analyzed in the 4-loop plant Safety Analysis Report (SAR). The results also demonstrated that the currently NRC-approved Westinghouse coupled code system such as ANCKVIPRE remains conservative.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Large-volume flux closure during plasmoid-mediated reconnection in coaxial helicity injection

A large-volume flux closure during transient coaxial helicity injection (CHI) in NSTX-U is demonstrated through resistive magnetohydrodynamics (MHD) simulations. Several major improvements, including the improved positioning of the divertor poloidal field coils, are projected to improve the CHI start-up phase in NSTX-U. Simulations in the NSTX-U configuration with constant in time coil currents show that with strong flux shaping the injected open field lines (injector flux) rapidly reconnect and form large volume of closed flux surfaces. This is achieved by driving parallel current in the injector flux coil and oppositely directed currents in the flux shaping coils to form a narrow injector flux footprint and push the injector flux into the vessel. As the helicity and plasma are injected into the device, the oppositely directed field lines in the injector region are forced to reconnect through a local Sweet-Parker type reconnection, or to spontaneously reconnect when the elongated current sheet becomes MHD unstable to form plasmoids. In these simulations for the first time, it is found that the closed flux is over 70% of the initial injector flux used to initiate the discharge. Furthermore, these results could work well for the application of transient CHI in devices that employ super conducting coils to generate and sustain the plasma equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transient variation of martian ground-atmosphere thermal boundary layer structure.

Results of a numerical simulation of the diurnal redistribution of temperature by radiative and molecular-conductive processes in the Martian soil-atmosphere system. An attempt is made to assess the importance of atmospheric molecular conduction near the surface and to estimate the characteristic depth of the diurnal temperature wave. The computational results are found to indicate a dual structure in the diurnal temperature wave propagation pattern, with a diffusive-type wave in the lowest 150 m superimposed on a radiatively induced disturbance with a characteristic scale of 1.8 km. Atmospheric molecular thermal conduction typically accounts for about 15% of the total heating/cooling in the lowest 25 m. Thermal conduction in both the soil and atmosphere appears to be an important factor in the thermal coupling of these subsystems. A free-convection regime in the conduction layer is predicted by the model for about five hours of the Martian day.

Pallmann, A. J.↗

Transient, hypervelocity flow in an axisymmetric nozzle

The performance of an axisymmetric nozzle designed to produce uniform parallel flow with a nominal Mach number of 8 is examined. The nozzle employs a flow decoupling approach and is calibrated by measuring Pitot pressures near the exit plane for a number of operating conditions. The detailed starting process of the nozzle is simulated numerically, and a quasi-one-dimensional method is used for investigating the transient wave structure and the nozzle starting time over a number of tunnel operating conditions. An axisymmetric Navier-Stokes code is utilized in order to illustrate the multidimensional and viscous aspects of the starting process. Although the steady performance observed is found adequate, settling times are much longer than predicted by the classical quasi-one-dimensional models.

Jacobs, P. A.↗

Pulse-by-pulse transient thermal deformation in crystal optics under high-repetition-rate FEL

Time-domain modeling of the thermal deformation of crystal optics can help define acceptable operational ranges across the pulse-energy repetition-rate phase space. In this paper, we have studied the transient thermal deformation of a water-cooled diamond crystal for a cavity-based X-ray free-electron laser (CBXFEL), either an X-ray free-electron laser oscillator (XFELO) or a regenerative amplifier X-ray free-electron laser (RAFEL), by numerical simulations including finite-element analysis and advanced data processing. Pulse-by-pulse transient thermal deformation of a 50 µm-thick diamond crystal has been performed with X-ray pulse repetition rates between 50 kHz and 1 MHz. Results for temperature and thermal deformation have been compared with the results of transient analysis using a continuous wave (CW) power loading. Temperature and thermal deformation results from pulse-by-pulse transient analysis vary with time about the results for the CW case for the same average power. The variation amplitude increases with pulse energy and decreases with repetition rate. When the repetition rate increases to infinity, both temperature and thermal deformation converge to the results for the CW case. Two critical time scales for the operation of crystal optics in a CBXFEL are (1) first-turn time, i.e. the time for the XFEL pulse to complete the first turn around the cavity so that the crystal sees the recirculated XFEL pulse, and (2) period-end time, i.e. the time that the next electron bunch arrives for the amplification, so that the crystal outcouples the amplified FEL power. For the same average power, simulation results show that the crystal thermal deformation seen by the XFEL beam decreases with repetition rate at the first-turn time of a 300 m-long cavity and increases with repetition rate at the period-end time. For the wavefront preservation requirement of the crystal optics, a pulse-energy versus repetition-rate phase space has been established. The upper bounds of the pulse energy at both first-turn and period-end times decreases with repetition rate, especially at the period-end time. The upper bound of the thermal deformation of the crystal at the period-end time for any repetition frequency can be estimated from the CW case. For a water-cooled diamond crystal of dimension 5 mm × 5 mm × 0.05 mm, the time to reach a quasi steady-state is about 50 ms for temperature and 50 µs for thermal deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗