Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation, Modeling, Performance, and Visualization”

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.

574 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Scalable Contour Tree Computation by Data Parallel Peak Pruning

As data sets grow to exascale, automated data analysis and visualization are increasingly important, to intermediate human understanding and to reduce demands on disk storage via in situ analysis. Trends in architecture of high performance computing systems necessitate analysis algorithms to make effective use of combinations of massively multicore and distributed systems. One of the principal analytic tools is the contour tree, which analyses relationships between contours to identify features of more than local importance. Unfortunately, the predominant algorithms for computing the contour tree are explicitly serial, and founded on serial metaphors, which has limited the scalability of this form of analysis. While there is some work on distributed contour tree computation, and separately on hybrid GPU-CPU computation, there is no efficient algorithm with strong formal guarantees on performance allied with fast practical performance. Furthermore, we report the first shared SMP algorithm for fully parallel contour tree computation, with formal guarantees of O(lg V lg t) parallel steps and O(V lg V) work for data with V samples and t contour tree supernodes, and implementations with more than 30× parallel speed up on both CPU using TBB and GPU using Thrust and up 70× speed up compared to the serial sweep and merge algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Subtask 1.1 – Advanced Characterization of Unconventional Oil and Gas Reservoirs to Enhance CO2 Storage Resource Estimates

The Energy & Environmental Research Center (EERC) developed methods and approaches for advanced characterization of unconventional oil and gas reservoirs to enhance CO2 storage resource estimates using samples collected from organic-rich shales of the Bakken Formation. The project was funded under U.S. Department of Energy (DOE) Cooperative Agreement No. DE-FE0024233, with nonfederal cost share provided by Hitachi High Technologies America. The project participants included the EERC, Hitachi High Technologies America, the team at the National Energy Technology Laboratory’s (NETL’s) computed tomography (CT) scanning laboratory in Morgantown, West Virginia, and several NETL staff who have worked on CO2 storage resource estimation. Investigations on the use of subsurface geologic formations for CO2 storage have been ongoing for many years. In the last decade, improvements in directional drilling technology and hydraulic fracturing of organic-rich shale and other tight rock formations have opened new opportunities for CO2 storage, with the potential cobenefit of enhancing gas and oil recovery. Gas-bearing, organic-rich shales have been of particular interest as a carbon storage target because their ability to store large volumes of CO2 is implied by their ability to store methane; however, estimating the CO2 storage resource of organic-rich shales and other unconventional reservoirs is challenging because of the lack of detailed geologic and petrophysical data that are needed to develop improved volumetric equations. According to NETL’s method for estimating regional-scale CO2 storage resources, the key factors that control CO2 storage in geologic formations include void space within induced fractures, natural fractures and matrix pores, as well as the occurrence and distribution of organic matter (OM) and clay minerals. A key goal of this effort was to use field emission scanning electron microscopy (FESEM) analysis of fine-grain unconventional reservoir rocks coupled with advanced image analysis and machine learning techniques to quantify the occurrence and distribution of porosity, OM, and clay minerals. By better understanding the potential migration pathways of CO2 within tight rocks and quantifying the occurrence of sorptive components (OM and clays) along potential flow pathways, methods to estimate the CO2 storage resource of such rocks can be refined. The Bakken petroleum system (Bakken) was used as the source of rock samples for this effort because it includes organic-rich shales with contrasting degrees of thermal maturity as well as tight, nonshale clastic rocks with varying clay mineral content. To better understand the key components of interest in tight rock samples, the EERC worked closely with Hitachi to optimize the FESEM settings necessary to resolve the sample mineralogy, porosity and OM content at micro- to nanoscale resolution. The EERC also worked closely with Hitachi to apply and improve its advanced mineral identification and characterization system (AMICS) software so that it is better suited for mineral identification at the high resolution necessary to resolve features of interest in fine-grained rocks such as organic-rich shales. The EERC also applied advanced image analysis software packages, including Ilastik, to classify and quantify the key components of interest within the samples using FESEM imagery. These advanced techniques allowed for quantification and differentiation of the void space within fractures, pores and OM, as well as a technique to determine and quantify the key mineral and OM components surrounding pore spaces. Being able to better quantify the void space and the potential accessibility of CO2 to sorptive minerals or OM within tight rocks allows for improved estimation of the CO2 storage resource potential of organic-rich shales and other unconventional formations. The information derived from the image analysis techniques was used to refine NETL’s method for estimating regional-scale CO2 storage resources in collaboration with staff from NETL. An additional goal of this work was to better understand rates of CO2 permeation, flow pathways, and sorption in organic-rich shales and nonshale reservoir rocks and to evaluate potential changes in CO2 flow over time. The EERC conducted long-duration (several weeks) CO2 flow-through testing at reservoir conditions using Bakken shale and nonshale samples to assess CO2 permeation rates and retention/sorption within core plugs to represent bulk matrix behavior. Conventional isothermal adsorption testing over a range of pressure conditions was also performed on pulverized Bakken samples to better assess potential maximum sorptive capacity of the samples. The results of this testing illustrated the higher sorptive capacity of the Bakken shales when compared to the sorptive capacity of the nonshale reservoir rock, which contains very little OM when compared to the Bakken shales. The results also illustrate the differences in permeation behavior as a function of fracture flow versus matrix flow. In an attempt to visualize CO2 permeation and flow behavior in Bakken shale samples over time, the EERC collaborated with NETL’s CT Lab in Morgantown, West Virginia. A variety of tests were performed at the NETL lab on unfractured and fractured samples of the Bakken shale. The unfractured Bakken shale samples did not result in CO2 migration within the samples over the testing duration; thus imagery of CO2 migration pathways within the matrix of Bakken samples was not obtained. The results of the flow-through testing within fractured samples illustrate changes in CO2 permeation rates over time, with a trend of decreasing CO2 flow over time. The results of this effort suggest that FESEM analysis coupled with advanced image analysis and machine learning techniques are effective in characterizing and quantifying the mineralogy, porosity, fractures, and OM content of 2-D imagery at the high resolution needed to identify features of interest in tight rocks. In addition, the image analysis techniques allow for identification and quantification of different porosity types as well as the key mineral components and/or OM that CO2 would contact in the reservoir as it migrates through pores and fracture networks. This is key information needed to better estimate the CO2 storage resource potential of organic-rich shales and other tight reservoirs. Additional work is needed to determine the sorptive capacity of the individual sample components (OM and clay minerals) at reservoir conditions so that this approach and the refined CO2 storage resource equation can be applied to other shale and unconventional reservoirs. In addition, subsequent work should be performed to better understand the number of images needed to use stereological approaches to estimate the volume of sample components using 2-D imagery, to explore machine learning techniques to differentiate between induced vs. natural fractures, and to estimate the connectivity of pores and fracture networks using 2-D imagery. Being able to derive this information at the microscale would facilitate scaling up of the sample characterization data used for reservoir modeling and simulation to estimate the CO2 storage resource potential of tight formations.

03 NATURAL GAS↗
NASA NTRS2016

Crew Exploration Vehicle (CEV) (Orion) Occupant Protection

The purpose of this study was to determine the similarity between the response of the THUMS model and the Hybrid III Anthropometric Test Device (ATD) given existing Wright-Patterson (WP) sled tests. There were four tests selected for this comparison with frontal, spinal, rear, and lateral loading. The THUMS was placed in a sled configuration that replicated the WP configuration and the recorded seat acceleration for each test was applied to model seat. Once the modeling simulations were complete, they were compared to the WP results using two methods. The first was a visual inspection of the sled test videos compared to the THUMS d3plot files. This comparison resulted in an assessment of the overall kinematics of the two results. The other comparison was a comparison of the plotted data recorded for both tests. The metrics selected for comparison were seat acceleration, belt forces, head acceleration and chest acceleration. These metrics were recorded in all WP tests and were outputs of the THUMS model. Once the comparison of the THUMS to the WP tests was complete, the THUMS model output was also examined for possible injuries in these scenarios. These outputs included metrics for injury risk to the head, neck, thorax, lumbar spine and lower extremities. The metrics to evaluate head response were peak head acceleration, HIC15, and HIC36. For the neck, N (sub ij) was calculated. The thorax response was evaluated with peak chest acceleration, the Combined Thoracic Index (CTI), sternal deflection, chest deflection, and chest acceleration- 3 ms clip. The lumbar spine response was evaluated with lumbar spine force. Finally the lower extremity response was evaluated by femur and tibia force. The results of the simulation comparisons indicate the THUMS model had a similar response to the Hybrid III dummy given the same input. The primary difference seen between the two was a more flexible response of the THUMS compared to the Hybrid III. This flexibility was most pronounced in the neck flexion, shoulder deflection and chest deflection. Due to the flexibility of the THUMS, the resulting head and chest accelerations tended to lag the Hybrid III acceleration trace and have a lower peak value. The results of the injury metric comparison identified possible injury trends between simulations. Risk of head injury was highest for the lateral simulations. The risk of chest injury was highest for the rear impact. However, neck injury risk was approximately the same for all simulations. The injury metric value for lumbar spine force was highest for the spinal impact. The leg forces were highest for the rear and lateral impacts. The results of this comparison indicate the THUMS model performs in a similar manner as the Hybrid III ATD. The differences in the responses of model and the ATD are primarily due to the flexibility of the THUMS. This flexibility of the THUMS would be a more human like response. Based on the similarity between the two models, the THUMS should be used in further testing to assess risk of injury to the occupant.

Currie-Gregg, Nancy J.↗
NASA NTRS2000

Atomistic Method Applied to Computational Modeling of Surface Alloys

The formation of surface alloys is a growing research field that, in terms of the surface structure of multicomponent systems, defines the frontier both for experimental and theoretical techniques. Because of the impact that the formation of surface alloys has on surface properties, researchers need reliable methods to predict new surface alloys and to help interpret unknown structures. The structure of surface alloys and when, and even if, they form are largely unpredictable from the known properties of the participating elements. No unified theory or model to date can infer surface alloy structures from the constituents properties or their bulk alloy characteristics. In spite of these severe limitations, a growing catalogue of such systems has been developed during the last decade, and only recently are global theories being advanced to fully understand the phenomenon. None of the methods used in other areas of surface science can properly model even the already known cases. Aware of these limitations, the Computational Materials Group at the NASA Glenn Research Center at Lewis Field has developed a useful, computationally economical, and physically sound methodology to enable the systematic study of surface alloy formation in metals. This tool has been tested successfully on several known systems for which hard experimental evidence exists and has been used to predict ternary surface alloy formation (results to be published: Garces, J.E.; Bozzolo, G.; and Mosca, H.: Atomistic Modeling of Pd/Cu(100) Surface Alloy Formation. Surf. Sci., 2000 (in press); Mosca, H.; Garces J.E.; and Bozzolo, G.: Surface Ternary Alloys of (Cu,Au)/Ni(110). (Accepted for publication in Surf. Sci., 2000.); and Garces, J.E.; Bozzolo, G.; Mosca, H.; and Abel, P.: A New Approach for Atomistic Modeling of Pd/Cu(110) Surface Alloy Formation. (Submitted to Appl. Surf. Sci.)). Ternary alloy formation is a field yet to be fully explored experimentally. The computational tool, which is based on the BFS (Bozzolo, Ferrante, and Smith) method for the calculation of the energetics, consists of a small number of simple PCbased computer codes that deal with the different aspects of surface alloy formation. Two analysis modes are available within this package. The first mode provides an atom-by-atom description of real and virtual stages 1. during the process of surface alloying, based on the construction of catalogues of configurations where each configuration describes one possible atomic distribution. BFS analysis of this catalogue provides information on accessible states, possible ordering patterns, and details of island formation or film growth. More importantly, it provides insight into the evolution of the system. Software developed by the Computational Materials Group allows for the study of an arbitrary number of elements forming surface alloys, including an arbitrary number of surface atomic layers. The second mode involves large-scale temperature-dependent computer 2. simulations that use the BFS method for the energetics and provide information on the dynamic processes during surface alloying. These simulations require the implementation of Monte-Carlo-based codes with high efficiency within current workstation environments. This methodology capitalizes on the advantages of the BFS method: there are no restrictions on the number or type of elements or on the type of crystallographic structure considered. This removes any restrictions in the definition of the configuration catalogues used in the analytical calculations, thus allowing for the study of arbitrary ordering patterns, ultimately leading to the actual surface alloy structure. Moreover, the Monte Carlo numerical technique used for the large-scale simulations allows for a detailed visualization of the simulated process, the main advantage of this type of analysis being the ability to understand the underlying features that drive these processes. Because of the simplicity of the BFS method for e energetics used in these calculations, a detailed atom-by-atom analysis can be performed at any point in the simulation, providing necessary insight on the details of the process. The main objective of this research program is to develop a tool to guide experimenters in understanding and interpreting often unexpected results in alloy formation experiments. By reducing the computational effort without losing physical accuracy, we expect that powerful simulation tools will be developed in the immediate future, which will allow material scientists to easily visualize and analyze processes at a level not achievable experimentally.

Bozzolo, Guillermo H.↗
DOE OSTI2021

Quantifying Properties for a Mechanistic, Predictive Understanding of Aqueous Impact on Aging of Medium and Low Voltage AC and DC Cabling in Nuclear Power Plants

To address the gap in knowledge in understanding degradation in relevant service conditions, this project aimed to develop a mechanistic, predictive model of medium and low voltage cable failure based on the primary environmental degradation parameters of aqueous immersion time, temperature, and the oxidation extent. To do so, we took a two-fold approach toward evaluating degradation as related to the aqueous condition. First, we evaluated the chemical, mechanical, and electrical properties of polymers that comprise the cabling insulation under varied aqueous conditions at different temperatures. Secondly, we evaluated the performance of insulation materials under similar accelerated aging conditions in addition to the dielectric breakdown of the cabling under submersion conditions. Thermal oxidation and immersion of LDPE thin films were done using a Parr vessel and the extent of oxidation was monitored using the carbonyl index from ATR-FTIR spectrum, a ratio quantifying the carbonyl functionality produced on the LDPE film surface. Increasing the temperature, oxidation time, and oxygen pressure increased the water-vapor permeability and the carbonyl content. Thermogravimetric analysis of the films confirmed that there is indeed an initial weight loss of the smaller molecular weight molecules and volatiles. At higher temperatures, there is a secondary mass loss. Up to 70°C, the mass loss curves were similar until 80°C, when the aged films decreased at a much greater rate with respect to increasing temperature. This was seen in the 10, 50, 90% weight loss temperatures. The mechanical properties of sheet PP and HDPE as well as injection molded LDPE, HDPE, HDPP dog bones were plotted with respect to the time spent in the Parr vessel during thermal oxidation. For the sheet PP and HDPE, the UTS and the modulus of elasticity decreased while the elongation increased. The injection-molded LDPE, HDPE, and HDPP dog bones experienced a similar trend where the changes in the properties were within experimental uncertainty. This analysis of LDPE in dry and immersive oxygen-rich environments is an important baseline for future experiments looking at other degradation mechanisms. A predictive aging model for low-density polyethylene insulative cable housings was developed. This model revealed a combination of physical and chemical processes which lead to an increase in permeability and a decrease in strength of the insulator. In creating this predictive model, a gap of knowledge in the literature was revealed in two areas: understanding how the crystallinity of a polymer changes with time and visualizing the predicted pores formed through the insulator which leads to failure. Developing and using an adapted ATR-FTIR method, crystallinity was monitored and compared to bulk crystallinity found via DSC. Inhomogeneity in crystallinity changes were observed but the FTIR method but was found to be less reliable. Pore visualization was found utilizing electrical impedance spectroscopy saturated aged polyethylene films with synthesized citrate-capped gold nanoparticles, and chloroauric acid precursor. results indicate a decreasing impedance of the polyethylene films caused by the transport of ions through the film. SEM imaging was then performed for elemental analysis of the films and counter electrodes for the presence of ions. Chloride and gold ions were detected; however, no nanoparticles were found. This gives an estimated pore size of at least 0.3 nm through the aged film. Cyclic submergence of polyethylene and polypropylene in aqueous solutions of copper sulfate and Harrison’s solution were examined. It was discovered that aging in these mixed conditions and at 90°C did cause a small, yet significant increase to tensile strength for both unaged polyethylene and polypropylene. The PE and PP in cycled and submerged conditions did not have tensile strengths significantly different from dry-aged after 16 weeks. For both solutions, it was determined that capacitance increases with both water tree depth and cable temperature and is relatively independent of changing water tree AR. As for resistance, there is no apparent change between depth percentages of 10 and 80 for both solutions. This is because the water tree has not yet entered the conductor and thus there is no shorting current flowing through the tree yet. Between 80 and 100% water tree depth, it is observed that for both solutions there is a rapid decrease in resistance because the tree is now impacting the conductor shield and allowing shorting current to flow through it. The relationship between resistance and temperature was found to be inversed. As temperature increases, resistance decreases for both water tree solutions with distilled water having the most apparent difference across the range of temperatures simulated. With regards to water tree AR, resistance was found to be relatively independent at depth percentages of 10, 50, and 70, but there was some relationship at 100% depth for both solutions. At this depth, as water tree AR increased, resistance was found to also increase. Through the examination of voltage and electric field distribution plots, it was confirmed for EPR cables that increasing water tree depth leads to increasing distortion. The increase in shorting current at 100% depth for aqueous copper sulfate compared to distilled water was also visualized. Lastly, the relationship between shorting current and temperature for narrow water trees was visualized and validated using the distribution plots. This study has led to a better understanding of the effect of cable temperature, water tree solution, and geometry of the tree region on cable degradation. Furthermore, this study has documented simulation results of water tree degradation in EPR MV cables, which had been previously lacking in information. Future work in this area will add frequency into the mix and examine what effect it has on the rate of cable degradation as a result of water treeing. To compliment the work being carried by UMD on different polymer chemistries, harvested medium voltage (MV) cables were selected for accelerated aging study at ORNL. These medium voltage cables are representative of the vintages and materials that are currently in use at existing nuclear power plants (NPPs). These cables were immersed in water at 90°C and energized for a period of two years at elevated voltage. Partial discharge was measured during this period of time to track potential degradation. Unfortunately, due to the absence of failure, the extent of integration between the UMD techniques and the harvested MV insulation was limited. However, based on the findings and from the UMD research in the previous sections on polyethylene and polypropylene, follow-on characterization for harvested MV cable insulation should focus on sample preparation to take advantage of the UMD techniques to better understand mechanistic degradation in MV cable insulations. This would include: 1.) Utilization of solutions and impedance measurement insulation to track permeability in harvested and accelerated aging insulation samples, 2.) Adaptation of gold nanoparticle synthesis for electrical impedance spectroscopy with supporting SEM characterization to study pore formation in harvested and accelerated aged insulation samples, & 3.) ATR-FTIR crystallinity with multiple temperature dependent harvested MV cable insulation under different air and submerged accelerated aging. In addition, low voltage dielectric spectroscopy and high voltage dissipation factor, tan δ, could be implemented as a condition monitoring tool for harvested MV cable insulation in submerged environments.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2002

The footprints of visual attention in the Posner cueing paradigm revealed by classification images

In the Posner cueing paradigm, observers' performance in detecting a target is typically better in trials in which the target is present at the cued location than in trials in which the target appears at the uncued location. This effect can be explained in terms of a Bayesian observer where visual attention simply weights the information differently at the cued (attended) and uncued (unattended) locations without a change in the quality of processing at each location. Alternatively, it could also be explained in terms of visual attention changing the shape of the perceptual filter at the cued location. In this study, we use the classification image technique to compare the human perceptual filters at the cued and uncued locations in a contrast discrimination task. We did not find statistically significant differences between the shapes of the inferred perceptual filters across the two locations, nor did the observed differences account for the measured cueing effects in human observers. Instead, we found a difference in the magnitude of the classification images, supporting the idea that visual attention changes the weighting of information at the cued and uncued location, but does not change the quality of processing at each individual location.

Non-NASA Center↗
NASA NTRS2013

NASA Tech Briefs, September 2013

Topics include: ISS Ammonia Leak Detection Through X-Ray Fluorescence; A System for Measuring the Sway of the Vehicle Assembly Building; Fast, High-Precision Readout Circuit for Detector Arrays; Victim Simulator for Victim Detection Radar; Hydrometeor Size Distribution Measurements by Imaging the Attenuation of a Laser Spot; Quasi-Linear Circuit; High-Speed, High-Resolution Time-to-Digital Conversion; Li-Ion Battery and Supercapacitor Hybrid Design for Long Extravehicular Activities; Ultrasonic Low-Friction Containment Plate for Thermal and Ultrasonic Stir Weld Processes; High-Powered, Ultrasonically Assisted Thermal Stir Welding; Next-Generation MKIII Lightweight HUT/Hatch Assembly; Centrifugal Sieve for Gravity-Level-Independent Size; Segregation of Granular Materials; Ion Exchange Technology Development in Support of the Urine Processor Assembly; Nickel-Graphite Composite Compliant Interface and/or Hot Shoe Material; UltraSail CubeSat Solar Sail Flight Experiment; Mechanism for Deploying a Long, Thin-Film Antenna From a Rover; Counterflow Regolith Heat Exchanger; Acquisition and Retaining Granular Samples via a Rotating Coring Bit; Very-Low-Cost, Rugged Vacuum System; Medicine Delivery Device With Integrated Sterilization and Detection; FRET-Aptamer Assays for Bone Marker Assessment, C-Telopeptide, Creatinine, and Vitamin D; Multimode Directional Coupler for Utilization of Harmonic Frequencies from TWTAs; Dual-Polarization, Multi-Frequency Antenna Array for use with Hurricane Imaging Radiometer; Complementary Barrier Infrared Detector (CBIRD) Contact Methods; Autonomous Control of Space Nuclear Reactors; High-Power, High-Speed Electro-Optic Pockels Cell Modulator; Covariance Analysis Tool (G-CAT) for Computing Ascent, Descent, and Landing Errors; Enigma Version 12; Micrometeoroid and Orbital Debris (MMOD) Shield Ballistic Limit Analysis Program; Spitzer Telemetry Processing System; Planetary Protection Bioburden Analysis Program; Wing Leading Edge RCC Rapid Response Damage Prediction Tool (IMPACT2); ISSM: Ice Sheet System Model; Automated Loads Analysis System (ATLAS); Integrated Main Propulsion System Performance Reconstruction Process/Models. Phoenix Telemetry Processor; Contact Graph Routing Enhancements Developed in ION for DTN; GFEChutes Lo-Fi; Advanced Strategic and Tactical Relay Request Management for the Mars Relay Operations Service; Software for Generating Troposphere Corrections for InSAR Using GPS and Weather Model Data; Ionospheric Specifications for SAR Interferometry (ISSI); Implementation of a Wavefront-Sensing Algorithm; Sally Ride EarthKAM - Automated Image Geo-Referencing Using Google Earth Web Plug-In; Trade Space Specification Tool (TSST) for Rapid Mission Architecture (Version 1.2); Acoustic Emission Analysis Applet (AEAA) Software; Memory-Efficient Onboard Rock Segmentation; Advanced Multimission Operations System (ATMO); Robot Sequencing and Visualization Program (RSVP); Automating Hyperspectral Data for Rapid Response in Volcanic Emergencies; Raster-Based Approach to Solar Pressure Modeling; Space Images for NASA JPL Android Version; Kinect Engineering with Learning (KEWL); Spacecraft 3D Augmented Reality Mobile App; MPST Software: grl_pef_check; Real-Time Multimission Event Notification System for Mars Relay; SIM_EXPLORE: Software for Directed Exploration of Complex Systems; Mobile Timekeeping Application Built on Reverse-Engineered JPL Infrastructure; Advanced Query and Data Mining Capabilities for MaROS; Jettison Engineering Trajectory Tool; MPST Software: grl_suppdoc; PredGuid+A: Orion Entry Guidance Modified for Aerocapture; Planning Coverage Campaigns for Mission Design and Analysis: CLASP for DESDynl; and Space Place Prime.

Source record↗
DOE OSTI2020

Macro-physical Properties of Shallow Cumulus from Integrated ARM Observations (Final Report)

Fair-weather shallow cumuli (ShCu) play an important role in many climate-related processes. Irregular geometry of ShCu and their strong temporal and spatial variability make it challenging to observe ShCu holistically and to represent them correctly in climate models. To improve ShCu parameterizations, information on both vertically and horizontally resolved cloud properties is required. Commonly, the vertically resolved cloud properties are provided by zenith pointing lidar-radar observations with a very narrow field of view (FOV). Thus, these “pencil-beam” properties may not be representative of a larger surrounding area. Limited number of areal-averaged cloud properties, such as fractional sky cover (FSC), are offered typically by wide-FOV observations. The main goal of our project was to integrate advantages of the narrow-FOV (vertical structure of clouds) and wide-FOV (spatial arrangement of clouds) observations for an improved characterization of single-layer ShCu observed at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site for an 18-yr period (2000-2017). There are four major accomplishments of our project, First, an updated operational cloud classification for days with ShCu has been suggested and evaluated through a detailed comparison with the manually curated records. Our classification extends successfully the latest ARM cloud type Value Added Product (VAP) based on the Active Remote Sensing of Clouds (ARSCL) cloud product by incorporating both cloud fraction (CF) provided by narrow-FOV ceilometer data and FSC from wide-FOV images offered by a Total Sky Imager (TSI). Moreover, our classification allows one to identify impact of instrumentation changes at the SGP site, namely the transition to KAZRARSCL with the updated cloud radar, on the identification of periods with single-layer ShCu. Second, a new approach that resolves cloud area distributions for a given region (up to 4x4 km 2 ) has been suggested and cloud equivalent diameters (CEDs) have been estimated for the first time. These estimations have been performed over a wide range of cloud sizes (about 0.01–3.5 km) with high temporal resolution (30s) using wide-FOV TSI images and cloud base height (CBH) provided by complementary narrow-FOV lidar measurements. Our simple and computationally inexpensive approach offers a previously unavailable dataset for process studies in the convective boundary layer and evaluation of ShCu parameterizations in cloud-resolving models. Third, a long-term integrated record of ShCu macrophysical properties has been developed. The developed record represents the longest available compilation of events with ShCu and includes (i) a novel visualization of the spatial variability in cloud cover both along- and across-wind directions, (ii) updated estimates of narrow-FOV CF and wide-FOV FSC, (iii) updated narrow-FOV CBH, and (iv) complementary data, such as wind speed and direction from the 915-MHz Radar Wind Profiler (RWP) data. The developed record has been used successfully to assess conventional observational estimates of cloud cover and their sensitivity to the following two factors: (i) instrument-dependent cloud detection and data merging criteria and (ii) FOV configuration. Fourth, co-variability of the ShCu macrophysical properties and environmental parameters has been analyzed for a 3-yr period (2016-2018). Our initial analysis includes diurnal changes of FSCs obtained for clouds with small, moderate and large CEDs and several environmental parameters, such as lifted condensation level (LCL) and mixed layer height (zi). Preliminary results of our analysis suggest that the horizontal extent of ShCu is controlled substantially by the sign and magnitude of difference between these two parameters (zi-LCL): the CED tends to grow with increase of this difference (zi exceeds LCL). We have initiated relationships between the ShCu and key atmospheric parameters that control both the development and evolution of ShCu using our new data product, which combines effectively the advantages of narrow-FOV data offered by zenith pointing cloud radars and lidars and wide-FOV TSI images. While the latest instrumentation at the ARM sites may address these challenging relationships in the future, we believe that the historical ARM data at the SGP site has not yet been fully utilized. Overall, our data product can be used by researchers working on a wide range of climate-related projects. These projects may include (i) a comprehensive evaluation of outputs from the Large-Eddy Simulation (LES) and single-column models for their future improvement, (ii) the representativeness of “short-period” results obtained from the previous model and observational studies and (iii) the planning of future field campaigns with focus on improved understanding of the diurnal cycle of cumulus convection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

In situ study of nucleation and growth dynamics of Au nanoparticles on MoS 2 nanoflakes

Two-dimensional (2D) substrates decorated with metal nanoparticles offer new opportunities to achieve high-performance catalytic behavior. However, little is known on how the substrates control the nucleation and growth processes of the nanoparticles. This paper presents the visualization of dynamic nucleation and growth processes of gold nanoparticles on ultrathin MoS 2 nanoflakes by in situ liquid-cell transmission electron microscopy (TEM). The galvanic displacement resulting in Au nuclei formation on MoS 2 was observed in real time inside the liquid cell. We found that the growth mechanism of Au particles on pristine MoS 2 is in between diffusion-limited and reaction-limited, possibly due to the presence of electrochemical Ostwald ripening. A larger size distribution and more orientation variation is observed for the Au particles along the MoS 2 edge than on the interior. Differing from pristine MoS 2 , sulfur vacancies on MoS 2 induce Au particle diffusion and coalescence during the growth process. Density functional theory (DFT) calculations show that the size difference is because the exposed molybdenum atoms at the edge with dangling bonds can strongly interact with Au atoms, whereas sulfur atoms on the MoS 2 interior have no dangling bonds and weakly interact with gold atoms. In addition, S vacancies on MoS 2 generate strong nucleation centers that can promote diffusion and coalescence of Au nanoparticles. Furthermore, the present work provides key insights into the role of 2D materials in controlling the size and orientation of noble metal nanoparticles vital to the design of next generation catalysts.

2D materials↗
DOE OSTI2020

Evaluation of H-Canyon Ventilation Exhaust Tunnel Inspection at the Savannah River Site - 20052

The Department of Energy (DOE) H-Canyon Facility at the Savannah River Site (SRS) is operated by Savannah River Nuclear Solutions (SRNS). The Facility Structural Integrity Program [1] provides a process to inspect, evaluate, and document the conditions of passive safety related structural systems and components, their degradation mechanisms and their impact on the safety envelope. The Structural Integrity Program performs periodic visual inspections to confirm the systems and components can perform their safety function and recommend corrective actions before the function is compromised. The H-Canyon Exhaust Tunnel is part of the H-Canyon Ventilation System and is periodically inspected under the Structural Integrity Program. The Tunnel performs a Safety Class passive design feature function that is available 100% of the time. The underground reinforced concrete Tunnel directs Canyon process air from the Canyon to the Sand Filter System. The radiological and chemical airborne activity coupled with the physical configuration of the tunnel precludes a manned entry into the tunnel to perform the inspection. Unlike active systems that can be subjected to physical testing to confirm capability, qualification calculations, including modeling, is performed to assess passive design feature capabilities relative to operational, accident and natural phenomena hazards events. Input for these calculations come from knowledge of field conditions, design drawings and conservative assumptions. The knowledge of field conditions comes from concrete and soil testing activities as well as inspection data. The use of cameras on crawlers to obtain inspection data and how that data is used in qualification calculations is the focus of this discussion. Remote Tunnel inspections have been performed using cameras on a stick, referred to as pole cameras, prior to Calendar Year 2003 (CY2003). Since CY2003, six camera equipped crawlers have been used. Crawlers have evolved and improved by incorporating lessons learned from the previous inspections. The latest 2019 crawler successfully traveled the entire length of the tunnel between the Canyon and Sand Filter System and captured images that were previously unseen, thus providing key visual data for the Structural Integrity Program report [2]. The inspection data serves as 1) input to qualification calculations, 2) validates current conditions remain bounded by qualification calculations, and 3) documents a record of change over time. This data is vital for the effective evaluation of the tunnel structure and aids with predicting the life cycle service life safety envelope to protect the public, environment, and facility worker. The goal of the 2019 deployment was to use a more robust vehicle to ensure travel over obstacles, perform the inspection with higher resolution cameras, and have the ability to elevate the cameras to view surfaces previously obscured by abandon-in-place ducts in the tunnel. The inspection information will serve to strengthen the input basis for qualification calculations, versus using excessive conservative calculation assumptions. This 2019 inspection was compared with previous crawler visual evidence to identify change over time to project remaining service life and confirm field conditions remain bounded by qualification calculations. The results of the 2019 Exhaust Tunnel inspection continue to improve and provide a firm basis to make sound engineering judgements regarding tunnel capability and approaches in projecting remaining service life with an aging system and component. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Closed Loop Geothermal Working Group: GeoCLUSTER App, Subsurface Simulation Results, and Publications

To better understand the heat production, electricity generation performance, and economic viability of closed-loop geothermal systems in hot-dry rock, the Closed-Loop Geothermal Working Group -- a consortium of several national labs and academic institutions has tabulated time-dependent numerical solutions and levelized cost results of two popular closed-loop heat exchanger designs (u-tube and co-axial). The heat exchanger designs were evaluated for two working fluids (water and supercritical CO2) while varying seven continuous independent parameters of interest (mass flow rate, vertical depth, horizontal extent, borehole diameter, formation gradient, formation conductivity, and injection temperature). The corresponding numerical solutions (approximately 1.2 million per heat exchanger design) are stored as multi-dimensional HDF5 datasets and can be queried at off-grid points using multi-dimensional linear interpolation. A Python script was developed to query this database and estimate time-dependent electricity generation using an organic Rankine cycle (for water) or direct turbine expansion cycle (for CO2) and perform a cost assessment. This document aims to give an overview of the HDF5 database file and highlights how to read, visualize, and query quantities of interest (e.g., levelized cost of electricity, levelized cost of heat) using the accompanying Python scripts. Details regarding the capital, operation, and maintenance and levelized cost calculation using the techno-economic analysis script are provided. This data submission will contain results from the Closed Loop Geothermal Working Group study that are within the public domain, including publications, simulation results, databases, and computer codes. GeoCLUSTER is a Python-based web application created using Dash, an open-source framework built on top of Flask that streamlines the building of data dashboards. GeoCLUSTER provides users with a collection of interactive methods for streamlining the exploration and visualization of an HDF5 dataset. The GeoCluster app and database are contained in the compressed file geocluster_vx.zip, where the "x" refers to the version number. For example, geocluster_v1.zip is Version 1 of the app. This zip file also contains installation instructions. **To use the GeoCLUSTER app in the cloud, click the link to "GeoCLUSTER on AWS" in the Resources section below. To use the GeoCLUSTER app locally, download the geocluster_vx.zip to your computer and uncompress this file. When uncompressed this file comprises two directories and the geocluster_installation.pdf file. The geo-data app contains the HDF5 database in condensed format, and the GeoCLUSTER directory contains the GeoCLUSTER app in the subdirectory dash_app, as app.py. The geocluster_installation.pdf file provides instructions on installing Python, the needed Python modules, and then executing the app.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

PipeSight: A High-Performance Computing Platform for Pipeline Integrity Management

The Phase I feasibility study completed as part of this project has led to a number of innovative technologies being developed and has laid the foundation for a successful Phase II effort to commercialize a platform for managing the integrity of pipelines for the damage mechanisms of the new, hybrid-energy based economy. To ground the development efforts and direction of the project, an extensive market research and customer discovery effort was undertaken early in Phase I. Through this effort, a number of pipeline owners and operators were interviewed, and the following key findings were discovered about the pipeline industry: • Small pipeline operators do not have the central engineering groups necessary to perform their own independent analysis of inspection data, but instead rely on summarized tally sheets provided to them by inspection service providers. • The time it takes to go from an inspection to a completed engineering assessment, even for small segments of pipeline, can take anywhere from 30-120 days. During this delay, critical threats can (and have been known to) cause failures. • Uncertainty is often not accounted for in the assessment of pipeline integrity. The tally sheets provided by third-party service providers are almost always deterministic in nature, identifying threats that present a concern only to the current (not the future) integrity of the pipeline. • It is uncommon to apply the latest technologies to perform advanced assessments of damaged pipelines. There is a desire to use more advanced analysis capabilities to assess threats. Many pipeline operators indicated that they would often excavate a pipeline to perform an inspection and find that the damage was not as bad as they anticipated, thus using limited resources unnecessarily. Companies are not consistent in their use of inspection data to determine corrosion rates, and those that do only calculate deterministic corrosion rates. • The industry has prominently relied on time-based inspections but has recently started to transition to risk-based inspections. However, there appears to be no uniform guidance on how to do so while properly accounting for all sources of uncertainty. • Companies are not storing inspection data in a manner that allows for the ready determination of temporal trends. • Predictive maintenance principles and practices are beginning to be used by early adopters • Some pipelines are being re-purposed to transport different process fluids than they were designed for, e.g., H 2 and CO 2 rich process streams to serve the new hybrid-energy based economy, which are presenting new integrity concerns for the existing pipeline network that crisscrosses the United States. As a result of these discoveries, we were able to target the development efforts in Phase I to best serve the needs of the industry. In Phase I, we developed a way to correlate multiple large-scale scans of the pipeline to determine a probabilistic corrosion rate that accounts for all sources of error and uncertainty in the inspection process. This probabilistic corrosion rate can be used to predict the future thickness distribution of the pipe wall. We demonstrate how this analysis may be performed in an analytical fashion and has been implemented in such a manner that it can be readily distributed using GPU computing through integration of the Kokkos programming model. We also make a very novel extension of the analytical corrosion rate model to Bayesian Networks (an explainable AI technique) that can account for non-parametric distributions of corrosion rates. With the predictions made above for the probabilistic corrosion rate and corresponding future distribution of the pipe wall thickness, we can assess the integrity of the pipeline through the use of a probabilistic engineering assessment. We developed a novel screening data analysis approach that can rapidly identify ‘hotspots’ (local thin areas) where the integrity of the pipeline is a concern. Once more, we implemented this screening approach in C++ to leverage GPU computing via the Kokkos programming model. After the critical hotspots are identified, we developed a program that can automatically generate an advanced finite element model of the damaged regions. Since the number of damaged regions that require advanced analysis can number in the thousands, we integrated an open-source container-native workflow engine for orchestrating parallel jobs on the cloud. Initially, these advanced numerical models were only designed to account for loading due to internal pressure. However, in a slight pivot from the initial Phase I proposal, we developed a complete pipe stress analysis program (called Simflex) which can simulate the complete pipeline and its response to thermal expansion, pressure, thermal bowing, weight, wind, earthquake, support displacement, support friction and external forces. This pipe stress analysis program was written generically, to handle any piping system, but contains the features needed to model long pipelines (i.e., it incorporates a model for soil mechanics and can account for the nonlinear boundary conditions necessary to simulate long underground pipelines). This pipe stress analysis program can simulate any segment of the pipeline (simple or complex) under any set of conditions and loads, to determine the supplemental loads (axial forces and bending moments) at the location of damage. This enables the most accurate state of stress to be accounted for in the pipeline, which can prove critical when evaluating the integrity of a damaged region. In the process of developing the technologies to perform the integrity assessment of the pipeline, we also extended one of the industry standard approaches for performing the assessment of local thin areas that extend more in the circumferential direction than the longitudinal direction of the pipeline. This approach was presented to the API 579-1/AS ME FFS-1 steering committee in November 2021 for consideration in the next edition of the industry standard for Fitness-For-Service (expected to be released in 2023). To help pipeline operators make decisions with the results on any integrity assessment, we developed a new approach to the life-cycle management of pipelines which uses a Bayesian Decision Network. The network is designed to help pipeline operators plan and prioritize inspection activities and ultimately make smarter, more cost-effective decisions. The Bayesian approach accounts for all sources of uncertainty and carries them through to the final optimal decisions, providing a probabilistic framework for optimizing inspection intervals. The proof-of-concept networks developed in the feasibility study are complete, verified, and are focused on a subset of the pipeline. To expand this novel approach to the scale necessary for an entire network of pipelines in Phase II, we will leverage the DOE-funded Bengi solver for industrial-scale decision making with Bayesian Networks [22]. Once implemented, we will be able to provide the pipeline industry with a much-needed tool for optimal inspection planning using truly explainable artificial intelligence (XAI). To handle all of these advanced capabilities into a cloud-based platform, the architecture of the Equity Engineering Cloud (EEC) was extended to include Argo Workflows, a framework capable of distributing and managing a massive number of jobs that consume their own resources, such that thousands of serial finite element simulations can be run in parallel. As part of this substantial undertaking, we also integrated Argo Continuous Delivery (CD) into the EEC, to aid with the rapid prototyping and iterations that will be imperative to the success of the PipeSight platform’s Agile development process in Phase II. As part of the pipe stress analysis program, we also developed a custom visualizer that leverages the DOE-funded VTK visualization library. We added custom contouring capabilities and a means for interacting visually with both the inputs and outputs of the pipe stress analysis program. We also developed routines for automating the post-processing of the finite element simulations to determine if any failure criteria are met and to visualize the deformations, stresses and strains in ParaView using the exodus II file format (a subset of netCDF).

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

First Annual Report on Pulsed Thermal Tomography Nondestructive Evaluation of Additively Manufactured Reactor Materials and Components: Pulsed Thermal Tomography Nondestructive Examination of Additively Manufactured Reactor Materials and Components

Additive manufacturing (AM, or 3D printing) for commercial nuclear energy applications is an emerging method for cost-efficient manufacturing aimed at replacing aging nuclear reactor parts and reducing costs for new construction. Because of the geometry of metallic structures of interest for nuclear applications, which consist of planar primitives with no symmetry of revolution, limited options are available for nondestructive evaluation (NDE), either during or post manufacturing. Known material flaws in AM include low-density regions consisting of non-sintered powder, which have to be detected to ensure the safety of long-term performance nuclear reactor components. As a solution to NDE of AM, we are developing pulsed thermal tomography (PTT) models and depth inversion algorithms for 3D imaging. PTT has many advantages because the method is non-contact and allows for in-service, NDE of AM nuclear reactor parts. By analyzing transients of surface temperature response due to internal thermal resistances, one can obtain 3D reconstructions of material effusivity using a unique inversion algorithm developed at Argonne. This report provides results of preliminary performance evaluation of PTT capability in detection of flaws in metallic structures, as well imaging complex geometry AM structures. Performance of PTT was investigated through modeling and experiments, Computer models of PTT were developed using COMSOL Heat Transfer module. Defects were modeled as cylindrical flat bottom holes (FBH), which is a common model of calibrated material flaws in thermal tomography experiments. Materials considered in this study include stainless steel 316 (SS316), stainless steel 304 (SS304), and Inconel 718. Theoretical analyses were conducted to validate inversion of simulated PTT data with COMSOL for a plate. Subsequently, 3D reconstructions were performed on COMSOL simulations for FBH, revealing a decrease in spatial resolution over depth due to thermal diffusion. The results of this study show that the performance of the inversion algorithm for detecting smaller defects depends strongly on the depth of the defect as well as the incident heat flux. The size of detectable defect was estimated by fitting a Gaussian function to surface temperature profile. The criteria for detectability was taken as 20mK noise equivalent temperature difference (NETD), which is currently the sensitivity limit of high-performance infrared cameras. It was determined through computer simulations that the smallest detectable FBH in SS316 has a 50µm diameter and is located 0.5mm below the plate surface. Preliminary experimental studies were performed to evaluate PTT performance in detecting calibrated flaws using metallic specimens with AM defects simulated as drilled-in FBH’s. High strength Stainless Steel 316 and Inconel 718 alloys were considered, as well as lower grade Stainless Steel 304, Nickel 200, and Hastelloy C276. Specimens investigated in this report consisted of approximately 1/4in-thick plates made out of these alloys using conventional manufacturing methods. The diameters of FBH’s varied from 1mm to 8mm, and their depths below the plate flat surface varied between 1mm and 6mm. The size of the smallest FBH was limited to 1mm because conventional mechanical drills were used for creating the holes. PTT imaging results have shown that 1mm-diameter FBH located 1mm and 2mm below the surface were detectable. Larger size FBH were detectable at greater depth. For example, 6mm-diameter FBH could be detected at 8mm depth. Image contrast varied slightly between the specimens, with the best reconstructions obtained in SS316 and C276 plates. Finally, PTT capability in imaging several AM structures fabricated from Inconel 718 (IN718) powder feedstock was investigated. The AM structures have complex geometry, but do not have calibrated internal defects. The objective is to determine PTT settings, such as the total integration time, for imaging of representative AM structures. A new FLIR X8500sc IR camera with higher spatial resolution (1280x1024 at 181Hz frame rate) was integrated into the PTT laboratory setup. Software tools were developed to ensure compatibility of new camera files with legacy effusivity reconstruction codes. The new camera was used to obtain 3-D reconstructions of AM IN718 nozzle plate and a particle filter plate. It was shown that PTT can scan through 2/3in-thick Inconel 718 plates in approximately 15s. Several modes of 3-D data visualization were explored, including using ImageJ and MATLAB software packages. The results demonstrate the capability of PTT in 3-D imaging of complex geometry AM structures. Intermediate results demonstrated in this report will be developed further to investigate detection of flaws in AM specimens.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Leveraging 3D CAD and Virtual Reality in Design of the Calcine Disposition Project - 20512

Fluor Idaho's Calcine Disposition Project (CDP) is leveraging the power of three-dimensional (3D) visualization, light detection and ranging (Lidar), and virtual reality to improve designs and worker safety. With the use of computer-aided design (CAD) modeling software, the CDP has created 3D models of the facility and systems to aid in the design, development, and operation of its retrieval and transfer system. This by itself is nothing new and has been done in industry for years. The CDP project however is taking this a few steps farther with the use of Lidar and virtual reality software. The CRP is tasked with removing a radioactive granular material from stainless steel bins located in concrete vaults that were constructed in the late 1950's and early 1960's. While construction drawings are available, it is not certain that the drawings are as-built or how accurate they truly are. Part of the project requires precise placement of equipment on to bins around multiple pipe, electrical lines, and bin stiffeners. Relying solely on the existing drawings and models created from these drawings is highly risky and prone to failure. The first step to determine the precise location and orientation of all obstructions in the vault the project will use a Lidar system to scan the vault structure and precisely locate all vault and bin components. Once this information is obtained, the data will be integrated into the CAD model. The model will then be verified, and a true as-built model developed. Precise component location, size, and orientation will be used for final design and placement of retrieval components. The second step will be to take the verified 3D model and scan data and import it into virtual reality software creating a virtual world. While some work on the project requires all work to be completed remotely due to high radiation fields, some areas have low enough radiation levels that personnel can enter and perform work. With the use of a virtual reality headset, operators and technicians will be able to enter our virtual world and become familiar with the surroundings and perform work prior to ever entering the radiation environment. This training is expected to pay dividends by improving worker efficiency, reducing errors, and improving confidence that the work can be performed as expected. It will play a key aspect in as low as reasonably achievable (ALARA) principles. An additional aspect in using the Lidar scans and virtual environment, is that it will give engineers and managers an opportunity to revisit the facility, especially radiation areas, at any time to obtain design information and measurements that may be costly or impossible to obtain once the facility transitions into radiation operations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2015

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗
NASA NTRS2016

Forward Skirt Structural Testing on the Space Launch System (SLS) Program

Structural testing was performed to evaluate heritage forward skirts from the Space Shuttle program for use on the Space Launch System (SLS) program. One forward skirt is located in each solid rocket booster. Heritage forward skirts are aluminum 2219 welded structures. Loads are applied at the forward skirt thrust post and ball assembly. Testing was needed because SLS ascent loads are roughly 40% higher than Space Shuttle loads. Testing objectives were to determine margins of safety, demonstrate reliability, and validate analytical models. Two forward skirts were structurally tested using the test configuration. The test stand applied loads to the thrust post. Four hydraulic actuators were used to apply axial load and two hydraulic actuators were used to apply radial and tangential loads. The first test was referred to as FSTA-1 (Forward Skirt Structural Test Article) and was performed in April/May 2014. The purpose of FSTA-1 was to verify the ultimate capability of the forward skirt subjected to ascent ultimate loads. Testing consisted of two liftoff load cases taken to 100% limit load followed by an ascent load case taken to 110% limit load. The forward skirt was unloaded to no load after each test case. Lastly, the forward skirt was tested to 140% limit and then to failure using the ascent loads. The second test was referred to as FSTA-2 and performed in July/August of 2014. The purpose of FSTA-2 was to verify the ultimate capability of the forward skirt subjected to liftoff ultimate loads. Testing consisted of six liftoff load cases taken to 100% limit load followed by the six liftoff cases taken to 140% limit load. Two ascent load cases were then tested to 100% limit load. The forward skirt was unloaded to no load after each test case. Lastly, the forward skirt was tested to 140% limit and then to failure using the ascent loads. The forward skirts on FSTA-1 and FSTA-2 successfully carried all applied liftoff and ascent load cases. Both FSTA-1 and FSTA-2 were tested to failure by increasing the ascent loads. Failure occurred in the forward skirt thrust post radius. The forward skirts on FSTA-1 and FSTA-2 had nearly identical failure modes. FSTA-1 failed at 1.72 times limit load and FSTA-2 failed at 1.62 times limit load. This difference is primarily attributed to variation in material properties in the thrust post region. Test data were obtained from strain gages, deflection gages, ARAMIS digital strain measurement, acoustic emissions, and high-speed video. Strain gage data and ARAMIS strain were compared to finite element (FE) analysis predictions. Both the forward skirt and tooling were modeled. This allows the analysis to simulate the loading as close as possible to actual test configuration. FSTA-1 and FSTA-2 were instrumented with over 200 strain gages to ensure all possible failure modes could be captured. However, it turned out that three gages provided critical strain data. One was located in the post bore and two on the post radius. More gages were not specified due to space limitations and the desire to not interfere with the use of the ARAMIS system on the post radius. Measured strains were compared to analysis results for the load cycle to failure. Note that FSTA-1 gages were lost before failure was reached. FSTA-2 gages made it to the failure load but one of the radius gages was lost before testing began. This gage was not replaced because of the time and cost associated with disassembly of the test structure. Correlation to analysis was excellent for FSTA-1. FSTA-2 was not quite as good because there was more residual strain from previous load cycles. FSTA-2 was loaded and unloaded with 12 liftoff cases and two ascent cases before taking the skirt to failure. FSTA-1 only had two liftoff cases and one ascent case before taking the skirt to failure. The ARAMIS system was used to determine strain at the post radius by processing digital images of a speckled paint pattern. Digital cameras recorded images of the speckled paint pattern. ARAMIS strain results for FSTA-2 just prior to failure. Note a high strain location develops near the left side. This high strain compares well to analysis prediction for both FSTA-1 and FSTA-2. The strain at this location was also plotted versus limit load. Both FSTA-1 and FSTA-2 had excellent correlation between ARAMIS and analysis strains. Acoustic emission (AE) sensors were used to monitor for damage formation that may occur during testing (e.g., crack formation and growth or propagation). AE was very important because after disassembly of FSTA-1, a crack was observed in the ball fitting radius. The ball fitting did not crack on FSTA-2. AE data was used to reconstruct when the crack occurred. The AE energy versus time plot for FSTA. The energy increased considerably at 850 seconds (152% limit load), indicating a crack could have formed at this point. The only visual evidence found that could have corresponded to this was the crack that initiated in the ball fitting. The cracks in the forward skirt aluminum structures would likely have been lower energy due to a lower modulus and all that were found after failure correlated to occurring after the initial crack in the post radius. This was verified by high-speed cameras used to record the failure.

Lohrer, J. D.↗