Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D temperature 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 541 records · Page 30

Aircraft observations and cloud-microwave radiative model simulations of a deep convective system

A 3D cloud model is used to simulate the storm structure of the June 29, 1986 storm system observed during the Cooperative Huntsville Meteorological Experiment (Dodge et al., 1986), and the results are linked to a microwave radiative transfer model for simulating aircraft observations. The results of cloud and radiative model simulations are then compared with radar data and aircraft microwave observations, respectively. It is shown that the degree of agreement between model simulations and the aircraft microwave observations depends on the model-derived microphysics of cloud and precipitation. It is found that the calculated values of the brightness temperature (Tb) at 92 GHz are close to the observed temperatures, while the calculated Tb values at 18 and 37 GHz are lower than the observed Tb.

Yeh, H.-Y. M.↗

NASA Tech Briefs, December 2009

Topics include: A Deep Space Network Portable Radio Science Receiver; Detecting Phase Boundaries in Hard-Sphere Suspensions; Low-Complexity Lossless and Near-Lossless Data Compression Technique for Multispectral Imagery; Very-Long-Distance Remote Hearing and Vibrometry; Using GPS to Detect Imminent Tsunamis; Stream Flow Prediction by Remote Sensing and Genetic Programming; Pilotless Frame Synchronization Using LDPC Code Constraints; Radiometer on a Chip; Measuring Luminescence Lifetime With Help of a DSP; Modulation Based on Probability Density Functions; Ku Telemetry Modulator for Suborbital Vehicles; Photonic Links for High-Performance Arraying of Antennas; Reconfigurable, Bi-Directional Flexfet Level Shifter for Low-Power, Rad-Hard Integration; Hardware-Efficient Monitoring of I/O Signals; Video System for Viewing From a Remote or Windowless Cockpit; Spacesuit Data Display and Management System; IEEE 1394 Hub With Fault Containment; Compact, Miniature MMIC Receiver Modules for an MMIC Array Spectrograph; Waveguide Transition for Submillimeter-Wave MMICs; Magnetic-Field-Tunable Superconducting Rectifier; Bonded Invar Clip Removal Using Foil Heaters; Fabricating Radial Groove Gratings Using Projection Photolithography; Gratings Fabricated on Flat Surfaces and Reproduced on Non-Flat Substrates; Method for Measuring the Volume-Scattering Function of Water; Method of Heating a Foam-Based Catalyst Bed; Small Deflection Energy Analyzer for Energy and Angular Distributions; Polymeric Bladder for Storing Liquid Oxygen; Pyrotechnic Simulator/Stray-Voltage Detector; Inventions Utilizing Microfluidics and Colloidal Particles; RuO2 Thermometer for Ultra-Low Temperatures; Ultra-Compact, High-Resolution LADAR System for 3D Imaging; Dual-Channel Multi-Purpose Telescope; Objective Lens Optimized for Wavefront Delivery, Pupil Imaging, and Pupil Ghosting; CMOS Camera Array With Onboard Memory; Quickly Approximating the Distance Between Two Objects; Processing Images of Craters for Spacecraft Navigation; Adaptive Morphological Feature-Based Object Classifier for a Color Imaging System; Rover Slip Validation and Prediction Algorithm; Safety and Quality Training Simulator; Supply-Chain Optimization Template; Algorithm for Computing Particle/Surface Interactions; Cryogenic Pupil Alignment Test Architecture for Aberrated Pupil Images; and Thermal Transport Model for Heat Sink Design.

Source record↗

Cold Season Performance of the NU-WRF Regional Climate Model in the Great Lakes Region

As Earth’s largest collection of freshwater, the Laurentian Great Lakes have enormous ecological and socio-economic value. Their basin has become a regional hotspot of climatic and limnological change, potentially threatening its vital natural resources. Consequentially, there is a need to assess the current state of climate models regarding their performance across the Great Lakes region and develop the next generation of high-resolution regional climate models to address complex limnological processes and lake-atmosphere interactions. In response to this need, the current paper focuses on the generation and analysis of a 20-member ensemble of 3-km National Aeronautics and Space Administration (NASA)-Unified Weather Research and Forecasting (NU-WRF) simulations for the 2014-2015 cold season. The study aims to identify the model’s strengths and weaknesses; optimal configuration for the region; and the impacts of different physics parameterizations, coupling to a 1D lake model, time-variant lake-surface temperatures, and spectral nudging. Several key biases are identified in the cold-season simulations for the Great Lakes region, including an atmospheric cold bias that is amplified by coupling to a 1D lake model but diminished by applying the Community Atmosphere Model radiation scheme and Morrison microphysics scheme; an excess precipitation bias; anomalously early initiation of fall lake turnover and subsequent cold lake bias; excessive and overly persistent lake ice cover; and insufficient evaporation over Lakes Superior and Huron. The research team is currently addressing these key limitations by coupling NU-WRF to a 3D lake model in support of the next generation of regional climate models for the critical Great Lakes Basin.

Michael Notaro↗

Coaxial Ceramic Direct Ink Writing on Heterogenous and Rough Surfaces: Investigation of Core–Shell Interactions

In this work, coaxial conductor–ceramic direct ink writing enables the printing of sensitive or encapsulated materials onto heterogeneous and rough substrates. While encasing the core fluid within a stiff ceramic shell, continuity may be maintained, even while printing onto conventionally challenging substrates. Here, we report the development of a coaxial ceramic direct ink writing suite and explore coflow interrelationships based on microfluidic principles. A coaxial nozzle is designed to facilitate the coextrusion of an alumina shell, whereas indium–tin-oxide inks constitute the core. In this manner, a core–shell ceramic element may be printed onto rough substrates for future high-temperature applications. Colloidal inks are engineered to provide the required rheological and sintering performance. Moreover, flow simulations in conjunction with microfluidic coflow principles are used to explore the coaxial printing processing space, thus controlling the core–shell architectures. Physical modeling is further used to analyze core deformations and eccentricity. As a result, simulations are validated experimentally, and the analyses are used to deposit coaxial ceramic features onto heterogeneous, high-temperature ceramic substrates.

20 FOSSIL-FUELED POWER PLANTS↗

Prediction of Dynamical Impact of Changes in Stratospheric Ozone

Under this grant one paper by Lou et al and a second paper by Kindler et al is in journal. These papers both describe N2O simulations using UKMO and Goddard assimilated wind fields and comparisons of the results against CLAES N2O observations. The results of these studies indicate some of the difficulties of using the assimilated wind fields, and the vertical motions in particular, in simulating long term variations in trace gases in the stratosphere. On the other hand, qualitatively the results possess a number of features of the observations even on time scales longer than a month or two. More recently we have started to examine results obtained using NCAR models a 3D version of which also uses the UKMO assimilated wind fields. Calculations have already been made with their 2D model with emphasis on the seasonal cycle in ozone at high latitudes in the upper stratosphere. Simultaneously trends in stratospheric ozone have been studied in detail from SAGE and UARS observations. Moreover, observations of the trends since 1984 do not show a significant interhemispheric asymmetry in upper stratospheric ozone. Therefore any asymmetry in the trends must have occurred prior to the mid-eighties and would most likely have been related to interhemispheric differences in upper stratospheric temperature trends. Another activity has been to compile an ozone climatology from UARS and SAGE observations. This effort has been performed as part of a UARS team activity to assemble a climatology of all the UARS long-lived trace gases for 1992-1993.

Cunnold, Derek M.↗

E-Theatre

Mesoscale/severe storms and results from Terra Fritz Hasler (NASA/Goddard) will show the latest spectacular animations of the hurricanes and tropical storms: Floyd, Georges, Mitch, Bonnie etc. from Geostationary Operational Environmental Satellites (GOES) and Tropical Rain Measuring Mission (TRMM) supported by MMS 3D nested numerical model results. See dust storms in Africa and smoke plumes from fires in Mexico. See how High-Definition Television (HDTV) is revolutionizing entertainment and the way we communicate scientific results. Fly in from space and take a tour of Africa. Drop in on Albuquerque using 1-m Ikonos 'Spy Satellite' technology. See data from the new Earth Observing System (EOS) Terra satellite as well as Landsat 7, and SeaWiFS. See climate change in action with Global Land and Ocean productivity changes over the last 20 years. Remote sensing observations of Sea Surface Temperature (SST), El Nino, height, winds and ocean color from GOES, National Oceanic and Atmosphere Administration Advanced Very High Resolution Radiometer (AVHRR), Special Sensor Microwave Imager (SSMI), and SeaWiFS put in context with ten year LANL global ocean model. Compare the Gulf Stream observed by Terra with the simulations. Additional information is contained in the original extended abstract.

Hasler, Arthur F.↗

Tungsten–carbon surface evolution and erosion modeling for a small angle slot divertor in DIII-D

Abstract We modeled tungsten–carbon mixed surface evolution, sputtering erosion, and transport for the tungsten coated region of a small angle slot (SAS) divertor design for the DIII-D tokamak. This divertor concept aims to achieve a closed slot dissipative plasma to minimize heat load and surface erosion, and to study high- Z material performance. Our advanced simulations use coupled ITMC-DYN material mixing/response and 3D full kinetic REDEP/WBC erosion/redeposition code packages, with divertor plasma solution from the SOLPS-ITER package with 4 MW power input. The SAS design geometry and resulting in-slot plasma parameters cause significant differences in predicted sputter/transport from a conventional divertor. For 2% C/D incident plasma ratio, an equilibrium mixed C/W surface is attained at ∼30 s of discharge, from wall sputtered carbon transported to the 10 cm long tungsten divertor region. Tungsten remains exposed to the plasma, but the evolved surface composition varies with different C/D ratios. Tungsten is primarily sputtered from the mixed surface by impinging carbon ions in the +1 to +4 charge states, with some self-sputtering. Redeposition of sputtered tungsten to the divertor is significant, ∼80% along the higher plasma temperature attached plasma SAS entrance region, but this is less than the typically near-unity values for a conventional divertor. Plasma-incident carbon is highly backscattered (∼50%) from the mixed surface, with little redeposition (<10%); this helps maintain tungsten in the surface sputter zone. Carbon is mainly sputtered from the mixed surface by D + ions, also with low redeposition (∼10%–30%). Finally, the modeling shows non-zero but low sputtered tungsten current from the divertor to the core plasma direction. These results appear favorable for effective testing of a tungsten-containing SAS divertor in DIII-D, and extrapolation of mixed-material evolution/response findings to the analogous low- Z /high- Z , Be/W, ITER plasma facing system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the 3D printability of one-part moisture-curable polyurethanes via direct ink writing (DIW)

Direct ink writing (DIW) is an ambient temperature additive manufacturing (AM) method that accommodates many possible ink materials. Here, in this work, we demonstrate using a moisture-curable commercially available polyurethane (PU) sealant as an alternative ink for DIW. We discussed the fundamentals of PU chemistry and determined the best 3D printing parameters. Studies on rheological, thermogravimetric, spectroscopic characterization, and initial finite element analyses (FEA) showed properties expected from a performance sealant with high elongation and low modulus of a 3D-printed object. This affirms the flexibility of the DIW technique as an accessible AM method amenable for future materials development from commercial model formulations.

36 MATERIALS SCIENCE↗

Analysis of core asymmetries in inertial confinement fusion implosions using three-dimensional hot-spot reconstruction

Three-dimensional effects play a crucial role during the hot-spot formation in inertial confinement fusion (ICF) implosions. A data analysis technique for 3D hot-spot reconstruction from experimental observables has been developed to characterize the effects of low modes on 3D hot-spot formations. In nuclear measurements, the effective flow direction, governed by the maximum eigenvalue in the velocity variance of apparent ion temperatures, has been found to agree with the measured hot-spot flows for implosions dominated by mode ℓ = 1. Asymmetries in areal-density (ρR) measurements were found to be characterized by a unique cosine variation along the hot-spot flow axis. In x-ray images, a 3D hot-spot x-ray emission tomography method was developed to reconstruct the 3D hot-spot plasma emissivity using a generalized spherical-harmonic Gaussian function. The gradient-descent algorithm was used to optimize the mapping between the projections from the 3D hot-spot emission model and the measured x-ray images along multiple views. Furthermore, this work establishes a platform to analyze 3D low-mode core asymmetries in ICF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design and Analysis of an Integrated Additively Manufactured Test Article for Plasma-Facing Components

Current plasma-facing components (PFCs) used in helium-cooled divertor modules are complex structures with tungsten tile, steel sleeve components, and cartridges, all assembled in a helium-cooled multiple jet (HEMJ) structure. The goal of this project is to simplify the complex PFC design using additive manufacturing techniques to create a single integrated tungsten test article. Apart from the flexibility this opens up in exploring a wide array of geometries for the article, having a single integrated article significantly reduces the number of joints and parts in the article, thus reducing chances of leaks. A process called electron beam melting has shown to produce very high-density samples and unique geometries, enabling HEMJ or similar designs. To validate and optimize this novel design, the model underwent a series of computational fluid dynamics and finite element analysis simulations to replicate steady-state heat flux in the divertors. The simulations presented in this study consider a steady-state base heat flux of 5 MW/m 2 , with water serving as the coolant. Future research will explore the use of helium as a coolant, simulate edge-localized-mode conditions, and include experimental validation. Since 3D-printed tungsten is anisotropic, the build direction versus build plane of the article are taken into consideration for the test article strength. Because of the high operating temperatures and low ductility of tungsten, thermal creep and brittle fracture are important failure mechanisms to consider. In conclusion, the cap is evaluated with various flow velocities and nozzle diameters, and an optimal design choice is made for which this cap will survive the divertor conditions with a conservative safety margin.

Divertor↗

Thermodynamic balance of three-dimensional stratospheric winds derived from a data assimilation procedure

The NASA/Goddard three-dimensional chemistry and transport model is driven by winds from a stratospheric data assimilation system. Synoptic- and planetary-scale patterns, apparent in satellite observations of trace constituents, are successfully reproduced for seasonal integrations. As model integrations proceed, however, the quality of simulations decreases, and systematic differences between calculation and measurement appear. The differences are explained by examining the zonal-mean residual circulation. The vertical velocity w-bar (sup star) is calculated two ways: (1) from the diabatic heating rates and temperature tendency and (2) from the Eulerian vertical velocity and the horizontal eddy heat flux convergence. The results from these calculations differ substantially. Periodic insertion of observational data during the assimilation process continually shocks the general circulation model and produces these differences, which leads to an overestimate of the mean vertical heat and cconstituent transport. Such differences are expected to be general to all data assimilation products. This interpretation is corroborated by two-dimensional (2D) model calculations. When w-bar(sup star) is calculated from (2), the 2D ozone evolution is unrealistic and qualitatively similar to the 3D model simulation. The 2D ozone evolution is reasonable when w-bar (sup star) from (1).

Weaver, Clark J.↗

NASA Tech Briefs, May 2007

Topics include: Noise-Canceling Helmet Audio System; Program Analyzes Spacecraft/Ground Radio Links; Two-Way Communication Using RFID Equipment and Techniques; Six-Message Electromechanical Display System; Scanning Terahertz Heterodyne Imaging Systems; Master Clock and Time-Signal-Distribution System; Synchronous Phase-Resolving Flash Range Imaging; Integrated Radial Probe Transition From MMIC to Waveguide; Bar-Code System for a Microbiological Laboratory; MMIC Amplifier Produces Gain of 10 dB at 235 GHz; Mapping Nearby Terrain in 3D by Use of a Grid of Laser Spots; Digital Beam Deflectors Based Partly on Liquid Crystals; Narrow-Band WGM Optical Filters With Tunable FSRs; Better Finite-Element Analysis of Composite Shell Structures; Computing Spacecraft-Pointing Vectors for Limb Tracking; Enhanced Master Controller Unit Tester; Rover Graphical Simulator; Increasing Durability of Flame-Sprayed Strain Gauges; Multifunctional, High-Temperature Nanocomposites; Multilayer Impregnated Fibrous Thermal Insulation Tiles; Radiation-Shielding Polymer/Soil Composites; Film/Adhesive Processing Module for Fiber-Placement Processing of Composites; Fabrication of Submillimeter Axisymmetric Optical Components; Electrochemical Disposal of Hydrazines in Water; Statistical Model of Evaporating Multicomponent Fuel Drops; Resistively Heated SiC Nozzle for Generating Molecular Beams; Compact Packaging of Photonic Millimeter-Wave Receiver; Diffractive Combiner of Single-Mode Pump Laser-Diode Beams; Wide-Band, High-Quantum-Efficiency Photodetector; A Robustly Stabilizing Model Predictive Control Algorithm; Modeling Evaporation of Drops of Different Kerosenes; Development of Vapor-Phase Catalytic Ammonia Removal System; Several Developments in Space Tethers; Design Concept for a Nuclear Reactor-Powered Mars Rover; Formation-Initialization Algorithm for N Spacecraft; and DNSs of Multicomponent Gaseous and Drop-Laden Mixing Layers Achieving Transition to Turbulence.

Source record↗

Stability Analysis for HIFiRE Experiments

The HIFiRE-1 flight experiment provided a valuable database pertaining to boundary layer transition over a 7-degree half-angle, circular cone model from supersonic to hypersonic Mach numbers, and a range of Reynolds numbers and angles of attack. This paper reports selected findings from the ongoing computational analysis of the measured in-flight transition behavior. Transition during the ascent phase at nearly zero degree angle of attack is dominated by second mode instabilities except in the vicinity of the cone meridian where a roughness element was placed midway along the length of the cone. The growth of first mode instabilities is found to be weak at all trajectory points analyzed from the ascent phase. For times less than approximately 18.5 seconds into the flight, the peak amplification ratio for second mode disturbances is sufficiently small because of the lower Mach numbers at earlier times, so that the transition behavior inferred from the measurements is attributed to an unknown physical mechanism, potentially related to step discontinuities in surface height near the locations of a change in the surface material. Based on the time histories of temperature and/or heat flux at transducer locations within the aft portion of the cone, the onset of transition correlated with a linear N-factor, based on parabolized stability equations, of approximately 13.5. Due to the large angles of attack during the re-entry phase, crossflow instability may play a significant role in transition. Computations also indicate the presence of pronounced crossflow separation over a significant portion of the trajectory segment that is relevant to transition analysis. The transition behavior during this re-entry segment of HIFiRE-1 flight shares some common features with the predicted transition front along the elliptic cone shaped HIFiRE-5 flight article, which was designed to provide hypersonic transition data for a fully 3D geometric configuration. To compare and contrast the crossflow dominated transition over the HIFiRE-1 and HIFiRE-5 configurations, this paper also analyzes boundary layer instabilities over a subscale model of the HIFiRE-5 flight configuration that was tested in the Mach 6 quiet tunnel facility at Purdue University.

Li, Fei↗

Hardware in the Loop Testing of an Iodine-Fed Hall Thruster

CUBESATS are relatively new spacecraft platforms that are typically deployed from a launch vehicle as a secondary payload,1 providing low-cost access to space for a wide range of end-users. These satellites are comprised of building blocks having dimensions of 10x10x10 cm cu and a mass of 1.33 kg (a 1-U size). While providing low-cost access to space, a major operational limitation is the lack of a propulsion system that can fit within a CubeSat and is capable of executing high delta v maneuvers. This makes it difficult to use CubeSats on missions requiring certain types of maneuvers (i.e. formation flying, spacecraft rendezvous). Recently, work has been performed investigating the use of iodine as a propellant for Hall-effect thrusters (HETs) 2 that could subsequently be used to provide a high specific impulse path to CubeSat propulsion. Iodine stores as a dense solid at very low pressures, making it acceptable as a propellant on a secondary payload. It has exceptionally high ρIsp (density times specific impulse), making it an enabling technology for small satellite near-term applications and providing the potential for systems-level advantages over mid-term high power electric propulsion options. Iodine flow can also be thermally regulated, subliming at relatively low temperature ( less than100 C) to yield I2 vapor at or below 50 torr. At low power, the measured performance of an iodine-fed HET is very similar to that of a state-of-the-art xenon-fed thruster. Just as importantly, the current-voltage discharge characteristics of low power iodine-fed and xenon-fed thrusters are remarkably similar, potentially reducing development and qualifications costs by making it possible to use an already-qualified xenon-HET PPU in an iodine-fed system. Finally, a cold surface can be installed in a vacuum test chamber on which expended iodine propellant can deposit. In addition, the temperature doesn't have to be extremely cold to maintain a low vapor pressure in the vacuum chamber (it is under 10(exp -6) torr at -75 C), making it possible to 'cryopump' the propellant with lower-cost recirculating refrigerant-based systems as opposed to using liquid nitrogen or low temperature gaseous helium cryopanels. In the present paper, we describe testing performed using an iodine-fed 200 W Hall thruster mounted to a thrust stand and operated in conjunction with MSFCs Small Projects Rapid Integration and Test Environment (SPRITE) Portable Hardware In the Loop (PHIL) hardware. This work is performed in support of the iodine satellite (iSAT) project, which aims to fly a 200-W iodine-fed thruster on a 12-U CubeSat. The SPRITE PHIL hardware allows a given vehicle to do a checkout of its avionics algorithm by allowing it to monitor and feed data to simulated sensors and effectors in a digital environment. These data are then used to determine the attitude of the vehicle and a separate computer is used to interpret the data set and visualize it using a 3D graphical interface. The PHIL hardware allows the testing of the vehicles bus by providing 'real' hardware interfaces (in the case of this test a real RS422 bus) and specific components can be modeled to show their interactions with the avionics algorithm (e.g. a thruster model). For the iSAT project the PHIL is used to visualize the operating cycle of the thruster and the subsequent effect this thrusting has on the attitude of the satellite over a given period of time. The test is controlled using software running on an Andrews Space Cortex 160 flight computer. This computer is the current baseline for a full iSAT mission. While the test could be conducted with a lab computer and software, the team chose to exercise the propulsion system with a representative CubeSat-class computer. For purposes of this test, the "flight" software monitored the propulsion and PPU systems, controlled operation of the thruster, and provided thruster state data to the PHIL simulation. Commands to operate the thruster were initiated from an operator's workstation outside the vacuum chamber and passed through the Cortex 160 to exercise portions of the flight avionics. Two custom-designed pieces of electronics hardware have been designed to operate the propellant feed system. One piece of hardware is an auxiliary board that controls a latch valve, proportional flow control valves (PFCVs) and valve heaters as well as measuring pressures, temperatures and PFCV feedback voltage. An onboard FPGA provides a serial link for issuing commands and manages all lower level input-output functions. The other piece of hardware is a power distribution board, which accepts a standard bus voltage input and converts this voltage into all the different current-voltage types required to operate the auxiliary board. These electronics boards are located in the vacuum chamber near the thruster, exposing this hardware to both the vacuum and plasma environments they would encounter during a mission, with these components communicating to the flight computer through an RS-422 interface. The auxiliary board FPGA provides a 28V MOSFET switch circuit with a 20ms pulse to open or close the iodine propellant feed system latch valve. The FPGA provides a pulse width modulation (PWM) signal to a DC/DC boost converter to produce the 12-120V needed for control of the proportional flow control valve. There are eight MOSFET-switched heating circuits in the system. Heaters are 28V and located in the latch valve, PFCV, propellant tank and propellant feed lines. Both the latch valve and PFCV have thermistors built into them for temperature monitoring. There are also seven resistance temperature device (RTD) circuits on the auxiliary board that can be used to measure the propellant tank and feedline temperatures. The signals are conditioned and sent to an analog to digital converter (ADC), which is directly commanded and controlled by the FPGA.

Polzin, Kurt A.↗

Vidyut3d: A Gpu Accelerated Fluid Solver for Non-Equilibrium Plasmas on Adaptive Grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure twin electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate approximately 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

Sitaraman, Hariswaran↗

SAM Code Enhancements for Modeling of Liquid Metal-Cooled Fast Reactor Concepts

The SAM code is under development and supported by DOE-NE’s Nuclear Energy Advanced Modeling and Simulation program as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. These advanced reactor concepts incorporate novel and improved approaches to achieve safety and economic feasibility. This report summarizes two major efforts in addressing the modeling gaps in SAM for liquid-metal-cooled fast reactor (LMFR) applications, i.e. thermal mixing and stratification phenomena in large pools and corrosion-oxidation of components in flowing lead. A new one-dimensional model for thermal mixing and stratification effects in large pools and enclosures is developed and implemented. Thermal mixing and stratification occur when fluid enters a pool at a temperature different than the bulk fluid itself, a scenario often encountered during transients in pool-type systems. These phenomena are critical for the safety of reactors, impacting phenomena like natural circulation, which is essential for passive cooling. The improved model in SAM addresses limitations of state-of-the-art approaches by combining one-dimensional (1D) channels, representing the coolant jet flow, with lumped-parameter zero-dimensional (0D) pools, representing the rest of coolant in the tank. Energy exchange between the 1D jet and the 0D pools is based on heat transfer correlations calibrated against 3D simulations. It is verified that this model can handle various flow configurations, including hot jets in colder pools, cold jets in hotter pools, and the presence of features like ceilings, free surfaces, and obstacles. Additionally, validation against experimental data demonstrates the ability of the model to capture mixing and stratification effects in a wide range of conditions. The flexibility and improved accuracy of the new model make it a valuable tool for reactor safety analysis, allowing for the simulation of different geometries encountered in advanced reactors. A system-level corrosion modeling capability is developed and implemented in SAM to support Lead Fast Reactor (LFR) development. Although the initial focus of this capability will be on LFR application, this can later be expanded to include other liquid metals such as Lead-Bismuth Eutectic (LBE) and PbLi. This report summarizes the common corrosion mitigation strategies and outlines the progress on implementing and validating a corrosion-oxidation model in SAM. Verification and validation of the corrosion-oxidation portion of the model was performed using analytical solution and measured data from samples tested in the non-isothermal pumped lead loop at IPPE Obninsk. The iron transport and corrosion/precipitation portion of the model was assessed using an analytical model and measured corrosion depths from a natural convection lead loop experiment performed at CEA. It is demonstrated that the model implemented in SAM performed well in these assessments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Io's volcanic and sublimation atmospheres

Fully 3D axisymmetric gasdynamic equations simulating SO2 and H2S frost sublimation and SO2 dayside and nightside volcanic atmospheres on Io are numerically solved, using a time-explicit finite-volume formulation. Both the sublimation and volcanic atmospheres generate horizontal supersonic winds away from the subsolar point or the volcanic vent. While the sublimation atmosphere is primarily driven by horizontal pressure gradients determined by surface temperatures, the volcanic atmosphere is driven by pressure gradients that are determined by the source rate. Sublimation and condensation produce patterns of surface deposits which are characteristic of the two types of atmospheres. The volcanic model is quantitatively consistent with Voyager observations of ring deposits.

Moreno, Miguel A.↗