Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transient Modeling”

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,045 records · Page 58

Development of the TREAT M2 experiment as a transient benchmark

A transient benchmark based on the Transient Reactor Test Facility (TREAT) M2 Calibration experiment (M2-CAL) is under development. TREAT, at Idaho National Laboratory, is a graphite moderated air-cooled research reactor which has been used extensively for fuel material testing under extreme and accident conditions. Accurate benchmark models are a beneficial component in the operations, experimental planning, and development of TREAT. In this work, we present a transient benchmark model for the M2-CAL experiment core loading. The benchmark model incorporates the coupling between the Monte Carlo Code SERPENT and the computational fluid dynamics code OpenFOAM to capture the temperature feedback mechanism. The M2-CAL transient 2580 was simulated in this work. A pre-transient analysis was performed to determine the optimum core composition and conditions before the beginning of the transient. The analysis was validated against historic TREAT kinetic measurements and the worth of the transient rod T-2. The axial power distribution in the flux wire for the M2-CAL experiment was determined and contrasted with the experiment. The model was then used to simulate the M2-CAL transient 2580 experiment based on the reported pre-transient and transient conditions. Several transient observables were calculated and compared to the experiment. The model shows good agreement with the experimental power traces, period, and average increase of power at the power ramp (time > 7.8 s). Inverse point kinetics analysis was performed during the period of the power ramp for the model and the experiment based on that, the temperature feedback component was isolated and contrasted with the experimental feedback. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Hot-spot tectonics of Eistla Regio, Venus: Results from Magellan images and Pioneer Venus gravity

Eistla Regio (ER) is a broad, low, discontinuous topographic rise striking roughly EW at low northern latitudes of Venus. Some 2000 x 7000 km in dimensions, it is the third largest rise in planform on Venus after Aphrodite Terra and Beta Phoebe Regiones. These rises are the key physiographic elements in a hot spot model of global tectonics including transient plume behavior. Since ER is the first such rise viewed by Magellan and the latitude is very favorable for Pioneer Venus gravity studies, some of the predictions of a time dependent hot spot model are tested. Western ER is defined as the rise including Gula and Sif Mons and central ER as that including Sappho Patera. Superior conjunction prevented Magellan from returning data on eastern ER (Pavlova) during the first mapping cycle. It is concluded that the western and central portions of ER, while part of the same broad topographic rise and tectonic framework, have distinctly different surface ages and gravity signatures. The western rise, including Gula and Sif Mons, is the expression of deep seated uplift with volcanism limited to the individual large shields. The eastern portion has been widely resurfaced more recently by thermal anomalies in the mantle.

Grimm, Robert E.↗

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↗

Transient periodicities in X-ray-active red dwarfs - First results from Mount Cuba and interpretation with an oscillating loop model

Results from a program at the Mount Cuba Astronomical Observatory to optically monitor transient periodicities in flare stars are reported. The data are analyzed for periodicities by means of a modified autocorrelation approach. A randomization technique is used to assess the statistical significance of periods. In AD Leonis, variability with amplitude 0.4 percent and a period of 4.1 min was detected during one run on March 6, 1991: the probability that this periodicity is due to chance is 10 exp -5. Further observations obtained on AD Leo within 1-2 hr on the same night showed no significant variability (0.01). It is shown that optical periodicities of the sort detected may be interpreted as arising from oscillations in coronal loops.

Mullan, D. J.↗

Development of a three-dimensional APOLLO3 neutrons deterministic scheme for the CABRI reactor

CABRI is an experimental reactor to study the fuel behavior during reactivity injection transients. These transients being highly multiphysics, the development of suitable modeling and simulation tools to simulate them is important for the optimization of the tests and the control of the experimental conditions. This paper focuses on the development of an APOLLO3 deterministic core calculation tool dedicated to the CABRI transient analysis. It represents the first stage of the incremental process for the implementation of a multiphysics time-dependent modeling of the CABRI transient. The neutron calculation scheme is based on a classical two-step approach. The first step consists of a 281-energy group calculation flux with the TDT-MOC (Method Of Characteristics) solver for cross-section space and energy (23 groups) collapsing for the CABRI different assembly clusters. The bias on a 2D core neutron calculation due to the self-shielding calculation and collapsing on a restricted pattern are investigated thanks to a comparison with a direct full 2D calculation on a quarter of core. The second step relies on a pin-resolved transport 3D transport core calculation with the SN solver MINARET. A progressive numerical validation process is followed to quantify the calculation biases on reactivity and reaction rates at each step using reference calculations with the stochastic code TRIPOLI4. The next development stage toward a multiphysics scheme will be the implementation of the 3D-kinetics equation resolution and the coupling with a core thermal-hydraulics model. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Boundary-layer transition detection with infrared imaging emphasizing cryogenic applications

This paper reviews the technique of boundary-layer transition detection using infrared (IR) imaging, emphasizing cryogenic wind-tunnel testing. With the exception of the low-temperature effects on the IR radiation, the discussion is relevant to conventional wind-tunnel and flight testing as well. At low temperatures, IR imaging encounters a reduction in the radiated energy throughout the IR spectrum, combined with a shift to longer wavelengths of the bulk of the radiation. This radiation behavior affects the minimum resolvable temperature difference (MRTD) of the IR imaging system because of its fixed wave band sensitivity. In the absence of commercial long wavelength IR imaging systems, operating at wavelengths longer than 13 micron, some measures can be taken to alleviate the problem caused by the MRTD limitation. The thermal signature of transition can be enhanced by allowing a small and controlled temperature increase of the wind-tunnel flow that induces a transient heat transfer to the model. This action temporarily reveals the model area under the turbulent regime through its higher heating rate compared with the laminar regime. The contrast between the areas exposed to the two regimes can be enhanced by subtraction of thermograms (the equilibrium thermogram from the transient thermogram). Further visual improvement can be obtained through shade stretching or binary shading.

Gartenberg, Ehud↗

Probabilistic Prediction of Geomagnetic Storms and the K p index

Geomagnetic activity is often described using summary indices to summarize the likelihood of space weather impacts, as well as when parameterizing space weather models. The geomagnetic index K p in particular, is widely used for these purposes. Current state-of-the-art forecast models provide deterministic K p predictions using a variety of methods – including empirically-derived functions, physics-based models, and neural networks – but do not provide uncertainty estimates associated with the forecast. This paper provides a sample methodology to generate a 3-hour-ahead K p prediction with uncertainty bounds and from this provide a probabilistic geomagnetic storm forecast. Specifically, we have used a two-layered architecture to separately predict storm (K p ≥ 5 – ) and non-storm cases. As solar wind-driven models are limited in their ability to predict the onset of transient-driven activity we also introduce a model variant using solar X-ray flux to assess whether simple models including proxies for solar activity can improve the predictions of geomagnetic storm activity with lead times longer than the L1-to-Earth propagation time. By comparing the performance of these models we show that including operationally-available information about solar irradiance enhances the ability of predictive models to capture the onset of geomagnetic storms and that this can be achieved while also enabling probabilistic forecasts.

79 ASTRONOMY AND ASTROPHYSICS↗

Ligand-Controlled Energetics and Charge Transfer in Pure and Doped Nanocrystals

The research in this award period initially focused on the spectroscopy and dynamics of spherical CdSe nanocrystals and CdSe nanoplatelets. In the later part of the grant period we turned our attention to InP-based nanocrystals. It has long been believed that the transient absorption signal from approximately spherical CdSe nanocrystals is dominated by conduction band state filling. However, this long-held belief has recently been challenged, based on femtosecond absorption measurements. However, the recent studies challenging the conventional wisdom do not account of finite rates of spin-lattice relaxation. Our work showed that this is a crucial error and that the recent results are misinterpreted, that is, the previous prevailing wisdom is correct. Another study focused on CdSe nanocrystal photochemistry used transient absorption (TA) spectroscopy to determine the spatial extents of CdSe nanoplatelet (NPL) excitons. Our work shows that the spatial extents of the excitons in the NPLs are far less than the physical dimension of the NPL. Using a model developed to understand the transient absorption spectroscopy, we obtain an average excitonic area of 21.2 ± 2.5 nm 2 , independent of the nanoplatelet size. Our work on InP/ZnSe and InP/ZnS core/shell nanocrystals shows that when there is a small lattice mismatch (InP-ZnSe, 3.5%) a coherent core-shell interface is obtained. In contrast, the InP-ZnS lattice mismatch is much larger, 8.3%. In this case, the experimental results showed best agreement with calculations in which lattice strain is ignored, indicating that the interfaces in InP/ZnS nanocrystals are largely incoherent.

14 SOLAR ENERGY↗

Development, validation, and verification of multi-pass thermo-mechanical welding simulations using the open-source MOOSE framework: NeT TG4 benchmark weldment

This study develops and validates a sequentially coupled thermo-mechanical welding simulation for the three-pass 316L stainless steel NeT TG4 benchmark weldment using the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) and the Nuclear Engineering Material model Library (NEML). A diffused ellipsoidal heat source was calibrated against thermocouple data and weld macrographs to accurately model the fusion zone geometry and transient thermal fields. Material hardening is represented using the Lemaitre-Chaboche mixed isotropic-kinematic hardening model, while four annealing models - no annealing, single-stage at 1050 °C and 1300 °C, and two-stage at 800 °C/1300 °C - were implemented to assess the impact of annealing models on the accuracy of the predicted welding-induced plasticity, distortions, and residual stresses. The predictions were validated against experimental measurements and benchmarked against results from commercial software, demonstrating that thermo-mechanical MOOSE welding simulations achieve comparable accuracy with enhanced computational efficiency. This work highlights the potential of using open-source finite element frameworks like MOOSE for advanced manufacturing simulations.

Ji, Wendy [Australian Nuclear Science and Technolo↗

Computational Thermal Hydraulics of a High-Performance Low-Enriched-Uranium Annular Target for HFIR Irradiation

Molybdenum-99 has historically been generated via isolation from fissioned highly enriched uranium (HEU) targets. Here, this isotope is in high demand due to its daily use across the world in radiopharmaceutical medical procedures. The primary objective of this work was to design and analyze an experimental target assembly containing one low-enriched-uranium (LEU) annular target for irradiation at the High Flux Isotope Reactor (HFIR). Efforts included incorporating spatially dependent energy sources from neutron and gamma interactions, quantifying thermal contact conductance at material interfaces, performing grid-independent studies, comparing turbulence models, and simulating various steady-state and transient scenarios relevant for irradiation qualification and eventual insertion. These models provide velocity, pressure, and temperature distributions in both space and time. Such results enable the selection of an appropriate irradiation location, fission rate density, and flow-limiting orifice size and demonstrate compliance with HFIR safety requirements such that insertion into the reactor can be approved. This analysis shows that across all scenarios, wetted surface temperatures remain below the coolant saturation temperature with no net vapor formation in the coolant. In every scenario, all components stay below 30% of the aluminum 6061 melting temperature. Computational fluid dynamics and system-level models predict peak target temperatures that agree within 4%, though the predicted axial location of the peak differs by about 10% of the heated length due to differences in flow development length. These results de-risk the irradiation of LEU (annular targets) and strengthen a domestic, HEU-independent 99 Mo supply by providing important fuel performance data to form the foundation for a robust licensing basis.

Molybdenum-99↗

Phase change in liquid face seals. II - Isothermal and adiabatic bounds with real fluids

Analytical studies of phase change effects in parallel and tapered liquid face seals are presented. An isothermal and adiabatic model of low Reynolds number flow are considered by numerical integration of the descriptive equations for a real fluid, and its thermodynamic properties are calculated for each step, using a computer program for the steam tables or fluid thermodynamic properties. It was shown that for low leakage rate the isothermal model is more accurate and for high leakage rates the adiabatic model is more accurate; that both models yield the same conclusions regarding stability; and that the transient of collapse is described by the adiabatic model which predicts a catastrophic collapse and then either failure or explosive return to a larger film thickness value. Finally, it is shown that converging seals may become unstable and the mass leakage rate is reduced significantly below the all liquid value when boiling occurs.

Hughes, W. F.↗

Characteristics of three-dimensional quasi-geostrophic transient eddy propagation in the vicinity of a simulated Gulf Stream

Output of an eddy-resolving model of the North Atlantic is diagnosed in the vicinity of the Gulf Stream (GS), using the three-dimensional quasi-geostrophic transient wave activity flux diagnostic developed by Plumb [1986]. The model North Atlantic shows large transient wave activity along the model GS and along the edge of the Mediterranean Outflow. The fluxes show advection of wave activity by the GS from the separation point and radiation of wave activity from the GS in the top 1000 m.

boundary↗

Transient liquid-crystal technique used to produce high-resolution convective heat-transfer-coefficient maps

In this transient technique the preheated isothermal model wall simulates the classic one-dimensional, semi-infinite wall heat transfer conduction problem. By knowing the temperature of the air flowing through the model, the initial temperature of the model wall, and the surface cooling rate measured at any location with time (using the fast-response liquid-crystal patterns recorded on video tape), the heat transfer coefficient can be calculated for the color isothermal pattern produced. Although the test was run transiently, the heat transfer coefficients are for the steady-state case. The upstream thermal boundary condition was considered to be isothermal. This transient liquid-crystal heat-transfer technique was used in a transient air tunnel in which a square-inlet, 3-to-1 exit transition duct was placed. The duct was preheated prior to allowing room temperature air to be suddenly drawn through it. The resulting isothermal contours on the duct surfaces were revealed using a surface coating of thermochromic liquid crystals that display distinctive colors at particular temperatures. A video record was made of the temperature and time data for all points on the duct surfaces during each test. The duct surfaces were uniformly heated using two heating systems: the first was an automatic temperature-controlled heater blanket completely surrounding the test duct like an oven, and the second was an internal hot-air loop through the inside of the test duct. The hot-air loop path was confined inside the test duct by insulated heat dams located at the inlet and exit ends of the test duct. A recirculating fan moved hot air into the duct inlet, through the duct, out of the duct exit, through the oven, and back to the duct inlet. The temperature nonuniformity of the test duct model wall was held very small. Test results are reported for two inlet Reynolds numbers of 200,000 and 1,150,000 (based on the square-inlet hydraulic diameter) and two free-stream turbulence intensities of about 1 percent, which is typical of wind tunnels, and up to 20 percent (using a grid), which is typical of real engine conditions.

Hippensteele, Steven A.↗

Dynamics modeling of molten salt reactor with reduced and expanded representations of delayed neutron precursors

Molten salt reactors (MSRs) present unique challenges in dynamic behavior due to the mobility of their fuel. In these reactors, delayed neutron precursors (DNPs) drift with the fuel circulation through the primary loop. As a result, a fraction of DNPs decays outside the core, effectively reducing the available delayed neutron population for reactivity control. Consequently, precise modeling of the distribution and behavior of DNPs is critical for accurate reactor dynamics simulations. In this study, the System Dynamics Analysis Tool (SDAT) was used to simulate a thermal-spectrum MSR under steady-state conditions and following transients. The effects of using reduced and expanded representations of DNPs with fewer or more groups than the conventional 6-group model were investigated. Their impact on the simulated distribution of precursors in the primary loop, reactivity loss value, and reactor response to transients was analyzed. Simulation results showed that reduced models lead to the loss of the actual DNPs distribution data, resulting in less accurate estimates of reactivity loss. Reactor power predictions using these reduced models showed significant deviations compared to those using the conventional 6-group model in transient simulations. Expanded models offered a more accurate representation of the distribution of DNPs and reactivity loss estimates. Reactor power predictions using expanded models showed minimal deviation from the conventional 6-group model during the simulated transients.

analysis↗

Dynamic Modeling of Solar Dynamic Components and Systems

The purpose of this grant was to support NASA in modeling efforts to predict the transient dynamic and thermodynamic response of the space station solar dynamic power generation system. In order to meet the initial schedule requirement of providing results in time to support installation of the system as part of the initial phase of space station, early efforts were executed with alacrity and often in parallel. Initially, methods to predict the transient response of a Rankine as well as a Brayton cycle were developed. Review of preliminary design concepts led NASA to select a regenerative gas-turbine cycle using a helium-xenon mixture as the working fluid and, from that point forward, the modeling effort focused exclusively on that system. Although initial project planning called for a three year period of performance, revised NASA schedules moved system installation to later and later phases of station deployment. Eventually, NASA selected to halt development of the solar dynamic power generation system for space station and to reduce support for this project to two-thirds of the original level.

Hochstein, John I.↗