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,189 records · Page 66

Incorporating two-dimensional conduction modeling techniques into an energy simulation program: The EnergyPlus radiant system example

Due to the inherent complexity of multi-dimensional solutions for transient conduction in a multilayered building element, most heat balance-based algorithms assume that conduction occurs only in one dimension—through the building element between the inside and outside environments. Usually, since the length and width of the elements are much greater than the thickness of the element, this is a valid assumption. However, there are some cases such as hydronic radiant slabs where this assumption may not lead to an accurate estimate of the energy consumption of the system. In a hydronic radiant slab, water tubes are embedded in a slab, and hot or chilled water is circulated to meet any heating or cooling loads. The energy consumption of such a system will depend highly on the temperature of the water being sent to the slab (i.e., produced by the boiler or chiller). The temperature of the water that must be sent to the slab will depend on the spacing of the tubes, the multilayer construction of the building element, and the boundary conditions to which the construction is exposed. Thus, to accurately predict the energy consumption of such a system, a transient two-dimensional heat conduction model must be integrated with a heat balance-based simulation. This paper demonstrates the potential loss of accuracy associated with the one-dimensional radiant system models by providing a comparison of one- and two-dimensional solutions using EnergyPlus.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A numerical and experimental investigation of electrothermal aircraft deicing

Experimental data were obtained for the deicing characteristics of a stationary UH-1H helicopter blade which had been fitted with an electrothermal deicer assembly. The tests were run in the NASA Lewis Icing Research Tunnel, and yielded transient temperature responses for the substrate, heater and abrasion shield at selected positions around the blade. The data at the abrasion shield-ice interface clearly documented when melting, shedding or refreezing occurred. Comparisons were made between the experimental data and a one-dimensional numerical model. The agreement was generally very good, with the simulations being shown to be capable of predicting the transient temperature responses along with phase change and ice shedding. At many blade positions, the model was capable of accurately simulating the thermal response of the electrothermal deicer assembly.

Leffel, K.↗

JPS heater and sensor lightning qualification

Simulated lightning strike testing of the Redesigned Solid Rocket Motor (RSRM) field joint protection system heater assembly was performed at Thiokol Corp., Wendover Lightning Facility. Testing consisted of subjecting the lightning evaluation test article to simulated lightning strikes and evaluating the effects of heater cable transients on cables within the systems tunnel. The maximum short circuit current coupled onto a United Space Boosters, Inc. operational flight cable within the systems tunnel, induced by transients from all cables external to the systems tunnel, was 92 amperes. The maximum open-circuit voltage coupled was 316 volts. The maximum short circuit current coupled onto a United Space Boosters, Inc. operational flight cable within the systems tunnel, induced by heater power cable transients only, was 2.7 amperes; the maximum open-circuit voltage coupled was 39 volts. All heater power cable induced coupling was due to simulated lightning discharges only, no heater operating power was applied during the test. The results showed that, for a worst-case lightning discharge, the heater power cable is responsible for a 3.9 decibel increase in voltage coupling to operational flight cables within the systems tunnel. Testing also showed that current and voltage levels coupled onto cables within the systems tunnel are partially dependant on the relative locations of the cables within the systems tunnel.

Cook, M.↗

Global Positioning System Simulator Field Operational Procedures

Global Positioning System (GPS) simulation is an important activity in the development or qualification of GPS signal receivers for space flight. Because a GPS simulator is a critical resource it is highly desirable to develop a set of field operational procedures to supplement the basic procedures provided by most simulator vendors. Validated field procedures allow better utilization of the GPS simulator in the development of new test scenarios and simulation operations. These procedures expedite simulation scenario development while resulting in scenarios that are more representative of the true design, as well as enabling construction of more complex simulations than previously possible, for example, spacecraft maneuvers. One difficulty in the development of a simulation scenario is specifying various modes of test vehicle motion and associated maneuvers requiring that a user specify some (but not all) of a few closely related simulation parameters. Currently this can only be done by trial and error. A stand-alone procedure that implements the simulator maneuver motion equations and solves for the motion profile transient times, jerk and acceleration would be of considerable value. Another procedure would permit the specification of some configuration parameters that would determine the simulated GPS signal composition. The resulting signal navigation message, for example, would force the receiver under test to use only the intended C-code component of the simulated GPS signal. A representative class of GPS simulation-related field operational procedures is described in this paper. These procedures were developed and used in support of GPS integration and testing for many successful spacecraft missions such as SAC-A, EO-1, AMSAT, VCL, SeaStar, sounding rockets, and by using the industry standard Spirent Global Simulation Systems Incorporated (GSSI) STR series simulators.

Kizhner, Semion↗

Computer program for analysis of high speed, single row, angular contact, spherical roller bearing, SASHBEAN. Volume 1: User's guide

The computer program SASHBEAN (Sikorsky Aircraft Spherical Roller High Speed Bearing Analysis) analyzes and predicts the operating characteristics of a Single Row, Angular Contact, Spherical Roller Bearing (SRACSRB). The program runs on an IBM or IBM compatible personal computer, and for a given set of input data analyzes the bearing design for it's ring deflections (axial and radial), roller deflections, contact areas and stresses, induced axial thrust, rolling element and cage rotation speeds, lubrication parameters, fatigue lives, and amount of heat generated in the bearing. The dynamic loading of rollers due to centrifugal forces and gyroscopic moments, which becomes quite significant at high speeds, is fully considered in this analysis. For a known application and it's parameters, the program is also capable of performing steady-state and time-transient thermal analyses of the bearing system. The steady-state analysis capability allows the user to estimate the expected steady-state temperature map in and around the bearing under normal operating conditions. On the other hand, the transient analysis feature provides the user a means to simulate the 'lost lubricant' condition and predict a time-temperature history of various critical points in the system. The bearing's 'time-to-failure' estimate may also be made from this (transient) analysis by considering the bearing as failed when a certain temperature limit is reached in the bearing components. The program is fully interactive and allows the user to get started and access most of its features with a minimal of training. For the most part, the program is menu driven, and adequate help messages were provided to guide a new user through various menu options and data input screens. All input data, both for mechanical and thermal analyses, are read through graphical input screens, thereby eliminating any need of a separate text editor/word processor to edit/create data files. Provision is also available to select and view the contents of output files on the monitor screen if no paper printouts are required. A separate volume (Volume-2) of this documentation describes, in detail, the underlying mathematical formulations, assumptions, and solution algorithms of this program.

Aggarwal, Arun K.↗

Open Source Modeling for Power Systems Dynamics

Time-domain simulations for power systems (both phasor and electromagnetic transients) are mostly dominated by commercial tools, such as PSS/E or PSCAD. However, with the increasing penetration of inverter-based resources and changes in model paradigm, new open-source tools are looking into become a possible alternative for a range of stakeholders. In this presentation we will focus on how to improve communication with industry and other stakeholders to properly showcase the value of open source simulators. Participants will identify current barriers facing the adoption of open-source time-domain tools, and will discuss how to improve the bench-marking processes across different tools, while maintaining software modularity to include new features and models as new challenges continue to happen in energy systems.

ENERGY PLANNING, POLICY, AND ECONOMY,MATHEMATICS A↗

Ring Model Development and Validation for Prismatic HTGR Core Thermal-hydraulics and Safety Analysis

Because of the complex core geometry, prismatic high temperature gas-cooled reactors (prismatic HTGRs) often exhibit complex thermal fluid behaviors during both normal operating and transient conditions. Most HTGR designs rely on passive safety system for decay heat removal, such as the reactor cavity cooling system (RCCS). During postulated accidents like Pressurized Conduction Cooldown (PCC) event, the decay heat is first radially transferred from the core region to the reactor vessel outer surface, then to the RCCS cooling panels. The peak fuel temperature is controlled by heat transfer mechanisms with two distinctive characteristic length scales, i.e., the core-wise effective heat conduction and local heat conduction in the fuel pellet scale. As both length scales are essential to determine the fuel temperature, from the modeling perspective, computer codes must be able to capture heat transfer in both scales. This is challenging for both computational fluid dynamics (CFD) tools because of extremely large amount of computation resources required, and for system analysis codes because of the challenge to model the complex core geometry. Under the support of DOE-NE’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, efforts have been pursued to support HTGR technology development and its modeling and simulation needs. There is a particular need for advanced modeling and simulation tools to predict thermal-fluid behavior during safety-related transients. In our previous studies, the ring model was adopted in SAM and further developed to simulate the normal operating condition and a PCC event using the MHTGR-350 design of General Atomics as the reference design. This current work represents a continuation of these previous efforts, and the focus is to critically review and examine simplifications and assumptions made to develop the ring model, and to perform code validation using experimental data from an integral-effect test facility, the High Temperature Test Facility (HTTF) at the Oregon State University. In this study, the test PG-27 from the HTTF test suite was selected for code validation purpose. The test PG-27 is a transient test designed to simulate the PCC event of the MHTGR design. Very good agreements between SAM prediction and experimental measurements were found in both coolant and solid structure temperatures during the transient.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Numerical simulations of microflare evolution in the solar transition region and corona

Several observers report transient ultraviolet brightenings, often referred to as microflares, in the solar atmosphere. In this paper, the results are presented of a series of one-dimensional numerical simulations examining possible relationships between microflares and the generation of dynamical chromospheric and transition region features. Low-energy and medium-energy microflares eject long-lived cool, dense gas plugs into the corona, with the gas plug traversing the loop apex in the medium energy case. In the case of high-energy microflares, the gas plug is rapidly heated to the temperature of the surrounding corona, and the results resemble the dynamics occurring in standard solar flare thick-target electron beam models.

Sterling, Alphonse C.↗

Possible evidence of excitonic condensation in a topological insulator

The transient excitonic condensate is a nonequilibrium electron-hole Bardeen-Cooper-Schrieffer state in a photoexcited semiconductor and semimetal, where electron-hole pairs undergo a phase transition and condense into a single coherent quantum state. Despite numerous experimental works to realize the predicted excitonic condensation phase, experimental evidence still remains elusive. This is largely due to the absence of direct measurements of a material's transient momentum-dependent electronic structure and the excitonic state in the condensation regime. Here, using time and angle-resolved photoemission spectroscopy, we find direct evidence of a transient excitonic condensate in the spin-polarized spatially indirect excitonic topological states in Bi2Te3. Accompanying the formation of the excitonic topological states by photoexcitation, we reveal a splitting of the hole's and electron's quasi-equilibrium chemical potential followed by the band flattening and backbending of the transient topological surface state. Moreover, within the same momentum range, we report a reshaping of the bulk valence band in the form of a Mexican-hat-like Bogoliubov dispersion-hallmarks of the excitonic condensation, followed by the opening of an energy gap at the Fermi level. The fluence and temperature dependence of these renormalization effects are reminiscent of excitonic condensation within Bardeen-Cooper-Schrieffer (BCS)-like behavior. These results, together with theoretical simulation, point to the possible formation of a transient excitonic condensate and provide opportunities to manipulate topologically protected Bose condensates with light.

ARPES↗

Evaluation of Test/Analysis Correlation Methods for Crash Applications

A project has been initiated to improve crash test and analysis correlation. The work in this paper concentrated on the test and simulation results for a fuselage section. Two drop tests of the section were conducted. The first test was designed to excite the linear structural response for comparison with finite element modal analysis results. The second test was designed to provide data for correlation with crash simulations. An MSC.Dytran model was developed to generate nonlinear transient dynamic results. Following minor modifications, the same model was executed in MSC.Nastran to generate modal analysis results. The results presented in this paper concentrate on evaluation of correlation methodologies for crash test data and finite element simulation results.

Lyle, Karen H.↗

Transient fuel performance analysis for the preliminary fuel concept of general atomics fast modular reactor

This study investigates the transient fuel performance of General Atomics Fast Modular Reactor (GA-FMR) during accident scenarios, focusing on the behavior of its innovative fuel system that combines high-assay low enriched uranium dioxide (HALEUO2) fuel with SiGA® ceramic matrix composite silicon carbide cladding. The preliminary fuel design’s response was analyzed during reactivity-initiated accidents (RIA) and loss of coolant accidents (LOCA) using BISON fuel performance analysis code, which included both the diffusion enhanced and BISON-FASTGRASS coupled UO 2 models. The RIA analysis demonstrated that effective reactivity control reduced fuel temperature, though with transient fission gas release resulting in additional tensile stress state on the cladding. LOCA simulations revealed differing predictions between the two models: the BISON UO 2 model showed more transient fission gas release but minimal pellet expansion, while the BISON-FASTGRASS UO 2 model predicted less pronounced fission gas release but more fuel swelling and thermal expansion, potentially leading to pellet-cladding mechanical interaction. Here, these findings highlight critical areas for fuel design optimization and identify knowledge gaps requiring further experimental and computational investigation to advance GA-FMR fuel development.

Lee, Soon K. [Argonne National Laboratory (ANL), A↗

Dynamic characteristics of a 20 kHz resonant power system - Fault identification and fault recovery

A detailed simulation of a dc inductor resonant driver and receiver is used to demonstrate the transient characteristics of a 20 kHz resonant power system during fault and overload conditions. The simulated system consists of a dc inductor resonant inverter (driver), a 50-meter transmission cable, and a dc inductor resonant receiver load. Of particular interest are the driver and receiver performance during fault and overload conditions and on the recovery characteristics following removal of the fault. The information gained from these studies sets the stage for further work in fault identification and autonomous power system control.

Wasynczuk, O.↗

Development of a Computational Framework for Multiphysics Multiphase Species Tracking using NEAMS Tools

This report implements a high-fidelity multiphysics modeling framework using the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program tools to track isotopic species in Molten Salt Reactors (MSRs), with a specific focus on the 91-depletion chain within the Molten Salt Reactor Experiment (MSRE). The model integrates neutronics, thermal-hydraulics, depletion, and thermochemistry to simulate the production, transport, and phase transitions of isotopes under steady-state and transient conditions. The main findings reveal that isotopes such as bromine-91 largely remain in the liquid phase, while others, including krypton-91and yttrium-91, transition to the gas phase, significantly influencing the reactor’s radiological source term. The study also shows that during transients, like a reactivity insertion transient, rapid void formation and the expansion of the liquid-gas interface led to substantial transfers of dissolved isotopes into the gas phase, altering isotope distribution and largely increasing the source term in the off-gas system. Additionally, the research highlights that short-lived isotopes dominate the initial off-gas response during transients, while longer-lived isotopes determine the equilibrium state, underscoring the necessity of dynamic simulations for accurate species tracking and reactor safety analysis. The developed methodology will be applied in the future to the tracking of a larger number of species and introduce other species tracking mechanisms, such as deposition and plating.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A three-dimensional transport model for the middle atmosphere

In this paper we describe fundamental properties of an 'off-line' three-dimensional transport model, that is, a model which uses prescribed rather than predicted winds. The model is currently used primarily for studying problems of the middle atmosphere because we have not (yet) incorporated a formulation for the convective transport of trace species, a prerequisite for many tropospheric problems. The off-line model is simpler and less expensive than a model which predicts the wind and mass evolution (an 'on-line' model), but it is more complex than the two-dimensional (2-D) zonally averaged transport models often used in the study of chemistry and transport in the middle atmosphere. It thus serves as a model of intermediate complexity and can fill a useful niche for the study of transport and chemistry. We compare simulations of four tracers, released in the lower stratosphere, in both the on- and off-line models to document the difference resulting from differences in modeling the same problem with this intermediate model. These differences identify the price to be paid in going to a cheaper and simpler calculation. The off-line model transports a tracer in three dimensions. For this reason, it requires fewer approximations than 2-D transport model, which must parameterize the effects of mixing by transient and zonally asymmetric wind features. We compare simulations of the off-line model with simulations of a 2-D model for two problems. First, we compare 2-D and three-dimensional (3-D) models by simulating the emission of an NO(x)-like tracer by a fleet of high-speed aircraft. The off-line model is then used to simulate the transport of C-14 and to contrast its simulation properties to that of the host of 2-D models which participated in an identical simulation in a recent NASA model intercomparison. The off-line model is shown to be somewhat sensitive to the sampling strategy for off-line winds. Simulations with daily averaged winds are in very good qualitative agreement but are less diffusive than when driven with instantaneous winds sampled at half-hour intervals. Simulations with the off-line and 2-D models are quite similar in the middle and upper stratosphere but behave quite differently in the lower stratosphere, where the 3-D model has a substantially more vigorous circulation. The off-line model is quite realistic in its simulation of C-14. While there are still systematic differences between the 3-D calculation and the observations, the differences seem to be substantially reduced when compared with the body of 2-D simulations documented in the above mentioned NASA intercomparison, particularly at 31 deg N.

Rasch, Philip J.↗

Thermal Testing and Model Correlation of the Magnetospheric Multiscale (MMS) Observatories

International Conference on Envronmental Systems (ICES), Seattle WA NCTS 20964-15. The Magnetospheric Multiscale (MMS) mission is a Solar Terrestrial Probes mission comprising four identically instrumented spacecraft that will use Earths magnetosphere as a laboratory tostudy the microphysics of three fundamental plasma processes: magnetic reconnection, energetic particle acceleration, and turbulence. This paper presents the complete thermal balance (TB) test performed on the first of four observatories to go through thermal vacuum (TV) and the minibalance testing that was performed on the subsequent observatories to provide a comparison of all four. The TV and TB tests were conducted in a thermal vacuum chamber at the Naval Research Laboratory (NRL) in Washington, D.C. with the vacuum level higher than 1.3 x 10-4 Pa (10-6 torr)and the surrounding temperature achieving -180 C. Three TB test cases were performed that included hot operational science, cold operational science and a cold survival case. In addition to the three balance cases a two hour eclipse and a four hour eclipse simulation was performed during the TV test to provide additional transient data points that represent the orbit in eclipse (or Earth's shadow) The goal was to perform testing such that the flight orbital environments could be simulated as closely as possible. A thermal model correlation between the thermal analysis and the test results was completed. Over 400 1-Wire temperature sensors, 200 thermocouples and 125 flight thermistor temperature sensors recorded data during TV and TB testing. These temperatureversus time profiles and their agreements with the analytical results obtained using Thermal Desktop and SINDAFLUINT are discussed. The model correlation for the thermal mathematical model (TMM) is conducted based on the numerical analysis results and the test data. The philosophy of model correlation was to correlate the model to within 3 C of the test data using the standard deviation and mean deviation error calculation. Individual temperature error goal is to be within 5 C and the heater power goal is to be within 5 of test data. The results of the model correlation are discussed and the effect of some material and interface parameters on the temperature profiles are presented.

Thermal↗

Thermal Testing and Model Correlation of the Magnetospheric Multiscale (MMS) Observatories

The Magnetospheric Multiscale (MMS) mission is a Solar Terrestrial Probes mission comprising four identically instrumented spacecraft that will use Earth's magnetosphere as a laboratory to study the microphysics of three fundamental plasma processes: magnetic reconnection, energetic particle acceleration, and turbulence. This paper presents the complete thermal balance (TB) test performed on the first of four observatories to go through thermal vacuum (TV) and the minibalance testing that was performed on the subsequent observatories to provide a comparison of all four. The TV and TB tests were conducted in a thermal vacuum chamber at the Naval Research Laboratory (NRL) in Washington, D.C. with the vacuum level higher than 1.3 x 10 (sup -4) pascals (10 (sup -6) torr) and the surrounding temperature achieving -180 degrees Centigrade. Three TB test cases were performed that included hot operational science, cold operational science and a cold survival case. In addition to the three balance cases a two hour eclipse and a four hour eclipse simulation was performed during the TV test to provide additional transient data points that represent the orbit in eclipse (or Earth's shadow) The goal was to perform testing such that the flight orbital environments could be simulated as closely as possible. A thermal model correlation between the thermal analysis and the test results was completed. Over 400 1-Wire temperature sensors, 200 thermocouples and 125 flight thermistor temperature sensors recorded data during TV and TB testing. These temperature versus time profiles and their agreements with the analytical results obtained using Thermal Desktop and SINDA/FLUINT are discussed. The model correlation for the thermal mathematical model (TMM) is conducted based on the numerical analysis results and the test data. The philosophy of model correlation was to correlate the model to within 3 degrees Centigrade of the test data using the standard deviation and mean deviation error calculation. Individual temperature error goal is to be within 5 degrees Centigrade and the heater power goal is to be within 5 percent of test data. The results of the model correlation are discussed and the effect of some material and interface parameters on the temperature profiles are presented.

Correlation↗

Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry

A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modelling of environmentally induced discharges in geosynchronous satellites

The NASCAP computer code was used to compute the charging and discharging characteristics of a typical communications satellite in geosynchronous orbit. For the case of a severe substorm satellite surface differential charging in sunlight was found to be substantially less than that required to produce discharges in ground simulation studies. A discharge process was postulated involving discharges triggered at edges (or imperfection) followed by discharges to space. The characteristics of such discharges was parametrically varied to evaluate the possible effects on the satellite. Results indicated that discharge characteristics inferred from satellite monitors could be caused by predicted space discharges, that single cell discharges to space can reduce surface potentially over entire satellite, and that low density electron trajectory computations indicate that discharge generated electrons do not return to the satellite by long trajectories. Current transients predicted do not agree with available ground simulation results indicating that additional work must be done both analytically and experimentally to understand and fully explain these discrepancies.

Stevens, N. J.↗