Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation test beds”

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 793 records · Page 44

Core-Shell Oxidative Aromatization Catalysts for Single Step Liquefaction of Distributed Shale Gas (Final Technical Report)

The objective of this project was to design and demonstrate a core-shell structured multifunctional catalyst to convert the light (dry) components of shale gas into liquid aromatic compounds (primarily benzene and toluene) in a single step. Operated in a modular oxidative aromatization system (OAS) under a cyclic redox scheme, the novel catalyst and process can significantly improve the value and transportability of distributed shale gas. Since the project started, each quarter addressed a different set of tasks related to the completion of the milestone detailed in the project award. The yearly summaries of these tasks are summarized below: Q1-Q4: • Conducted project planning and literature search. • Investigated a number of SHC redox catalysts using thermogravimetric analysis and fixed-bed reactor experiments. • Initiated process modeling towards generating two process models for the methane DHA base case and OAS process. • Developed DHA catalysts capable of producing >500 g/kg-cat-hr aromatics at 80% or greater aromatics selectivity at 700°C. Q5-Q8: • Developed alternative approaches with sequential bed configurations to enhance the aromatic yields based on OCM+DHA • Improved the zeolite synthesis efficiency by using the microwave-assisted technique and investigated the synthesis conditions on the zeolite yield, crystalline structure and morphology • Constructed a set of Aspen Plus process models with significant energy savings for OAS as compared to the base case non-oxidative DHA. • Adapted conventional hydrothermal method to be applicable to the microwave synthesizer unit for more efficient catalyst synthesis. • Studied the structure of the OCM catalyst and the dispersion of the carbonate in the redox reactions and in methane flow with Raman Spectroscopy. Q9-Q12: • Scaled up the catalyst synthesis with the microwave synthesis method. Based on its performance, procedural characterizations and catalytic performance testing were further conducted for the new microwave synthesized catalysts with the newly-developed product analysis procedure. • Developed the reaction system setup for the C2-DHA or OCM+DHA reaction product and achieved a better product collection-analysis method for the aromatic products with an improved carbon balance. The product from the OCM reaction exhibited complicated effects on the DHA catalyst. • Conducted additional OCM catalyst characterization using Near Ambient Pressure X-ray Photoelectron Spectroscopy and in situ Raman characterization • Validated the significant energy savings for OAS as compared to the base case non-oxidative DHA. Successfully set up the simulation model for the OCM+DHA+SHC reaction system based on the updated experimental results from NCSU. Q13-End of project: • Synthesized new zeolite catalysts by the microwave method, conducted characterizations (XRD, SEM, and TEM) and catalytic behavior testing. • Explored the “wet” C 2 H 6 and C 2 H 4 DHA reactions with using steam co-feed. A subsequent reduction as the regeneration step can regenerate the DHA catalyst and recover 99% activity of the fresh performance. • Achieved a 15.3% single-pass aromatic yield from methane by rationally combining the OCM and DHA at different temperatures. • Conducted a 105-hour stability test with an improved regeneration procedure, with an average aromatic yield of 13.8%. • Developed new catalyst and achieved a record-high 23.2% yield.

03 NATURAL GAS↗

Development of a Forced Oscillation System for Measuring Dynamic Derivatives of Fluidic Vehicles

A new Forced Oscillation System (FOS) has been designed and built at NASA Langley Research Center that provides new capabilities for aerodynamic researchers to investigate the dynamic derivatives of vehicle configurations. Test vehicles may include high performance and general aviation aircraft, re-entry spacecraft, submarines and other fluidic vehicles. The measured data from forced oscillation testing is used in damping characteristic studies and in simulation databases for control algorithm development and performance analyses. The newly developed FOS hardware provides new flexibility for conducting dynamic derivative studies. The design is based on a tracking principle where a desired motion profile is achieved via a fast closed-loop positional controller. The motion profile for the tracking system is numerically generated and thus not limited to sinusoidal motion. This approach permits non-traditional profiles such as constant velocity and Schroeder sweeps. Also, the new system permits changes in profile parameters including nominal offset angle, waveform, and associated parameters such as amplitude and frequency. Most importantly, the changes may be made remotely without halting the FOS and the tunnel. System requirements, system analysis, and the resulting design are addressed for a new FOS in the 12-Foot Low-Speed Wind Tunnel (LSWT). The overall system including mechanical, electrical, and control subsystems is described. The design is complete, and the FOS has been built and installed in the 12-Foot LSWT. System integration and testing have verified design intent and safe operation. Currently it is being validated for wind-tunnel operations and aerodynamic testing. The system is a potential major enhancement to forced oscillation studies. The productivity gain from the motion profile automation will shorten the testing cycles needed for control surface and aircraft control algorithm development. The new motion capabilities also will serve as a test bed for researchers to study and to improve and/or alter future forced oscillation testing techniques.

Trieu, B. C.↗

High-Temperature Gas-Cooled Reactors Multiphysics Simulation Demonstration and Code Validation

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor applications. The first part of this effort focuses on the running-in verification of Griffin’s multiphysics capabilities, specifically for simulating the evolution of pebble-bed reactor cores from startup to equilibrium. Since Fiscal Year 2024, improvements and enhancements have been implemented in Griffin, including simplifying the process to specify streamlines and developing the online cross-section generation capability. In the absence of validation data, code-to-code comparisons are conducted with kugelpy, showing good agreement for integral quantities like k-eff predictions and predictions for maximum power density. However, accuracy issues are noted for more detailed quantities like the spatial distribution of fission rate densities which will require further work to address. The second part of this report presents an improved System Analysis Module (SAM) core channel model where the effects of cross flow are considered during the pressurized loss of forced cooling transient, resulting in an improved agreement of the predicted pebble temperature with respect to the predictions from the SAM 2D porous media model. Additionally, the wall channeling effect due to variable porosity at the near wall region of the core is also investigated. Furthermore, to demonstrate Griffin’s online cross-section generation capability, a Multiphysics simulation is performed by coupling Griffin to the SAM core channel model. In the third part of the report, as a part of the Organisation for Economic Co-operation and Development/Nuclear Energy Agency (OECD/NEA) thermal-hydraulic code validation benchmark activity for a high-temperature gas-cooled reactor, the High Temperature Test Facility (HTTF) is investigated first using the NekRS computational fluid dynamics (CFD) code to study the flow mixing phenomenon in the lower plenum of the facility. Then, code-to-code and code-to-data comparisons are performed for Test PG27, which is a pressurized conduction cooldown (PCC) test, using five different codes by six organizations from five countries. The different simulations show good agreements in terms of the general trend but there are differences in some results such as the peak temperatures of different regions and heat removal rate.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SAM Code Development for Source Term Modeling in Fluoride-salt-cooled High-temperature Reactors

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. Many concepts leverage higher operating temperatures for improved efficiency with a variety of coolants and structural materials to support those needs. Such design choices may also pose unique radiological source term risks. There are continuous efforts to incorporate new physics and phenomena relevant to advanced reactor concepts, including capabilities to enable system-level source term modeling. The Pebble-Bed Fluoride-salt-cooled High-temperature Reactor (PB-FHR) is a promising candidate among advanced nuclear reactor concepts with its improved passive safety characteristics and high thermal efficiency. In addition to past efforts to support the development and utilization of SAM for PB-FHR safety analysis, the species transport modeling capabilities in SAM have been extended to simulate source term phenomena in concepts like the PB-FHR. The PB-FHR concept utilizes pebble-form TRISO fuel and FLiBe salt coolant to provide robust barriers to the release of almost all fission products and radiological source terms. However, tritium poses a unique radiological risk due to its significant production from neutron interactions with 6 Li and 9 Be in the FLiBe as well as its high mobility at elevated temperatures where it can permeate through structural metals. Secondly, the use of graphite at elevated temperatures poses the risk of oxidation damage to any structural components or fuel pebbles that are exposed to accidental air ingress. This report summarizes progress made in modeling these phenomena in SAM, which leverages the system-level multi-physics thermal hydraulic simulation to support an effective engineering-scale source term modeling capability. A tritium transport model is developed to simulate the various transport pathways in FHR and MSR concepts and is presented through verification and validation examples as well as demonstrations simulating experimental test loops and reference plant FHR models. A graphite oxidation model is also introduced to calculate local oxidation rates and is presented with initial validation results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

High-fidelity multiphysics load following and accidental transient modeling of microreactors using NEAMS tools: Application of NEAMS codes to perform multiphysics modeling analyses of micro-reactor concepts

The feasibility of modeling microreactors using high-fidelity models with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) tools is investigated in this report. Three overarching questions guided this research: can NEAMS tools readily be applied for high-fidelity multiphysics modeling of different types of transients in microreactor designs; how accurate are the results obtained; and are improvements needed in accuracy or user experience of NEAMS tools, especially considering newly developed capabilities? This work builds upon FY-2022 work, and two microreactor concepts considering heat pipe (HP-MR) and gas-cooled (GC-MR) technologies were further analyzed using high-fidelity multiphysics simulations. The NEAMS tools considered and coupled within the MultiApp environment are Griffin for neutronics, BISON for thermo-mechanics, Sockeye for heat pipe modeling (in HP-MR), SAM for 1D Fluid – 3D solid modeling of coolant channels and system modeling of balance of plant components (in GC-MR), and the SWIFT code for hydrogen redistribution in hydride moderator. The Heat Pipe MicroReactor (HP-MR) concept was further analyzed in FY-2023 to demonstrate the stochastic TRISO failure modeling capability in BISON to check operational limits of the TRISO fuel. A new full-core Gas-Cooled MicroReactor (GC-MR) model was developed based on the initial assembly-model used in Y-2022 and used for steady-state and accidental depressurization transient simulations. Accuracy of the simulations performed was assessed through 1) verification analyses completed on the different physics with code-to-code comparison, and 2) validation of the multiphysics simulations based on modeling of the Kilopower Reactor Using Stirling Technology (KRUSTY) experiment. In FY-2023, the mesh and model of KRUSTY was updated to closely match publicly available data, and the neutronic model was verified and validated against experimental control rod worth measurements. The multiphysics model of KRUSTY was developed and used for steady-state analysis and for modeling reactivity insertion transient. The calculated power increase and stabilization agrees well with experimental data following adjustment in fuel thermal expansion coefficient. As an important component of this project, the ANL team gathered experience with a wide range of NEAMS tools: the MOOSE Mesh System, Griffin, BISON, SWIFT, Sockeye, SAM, Workbench, and the MOOSE MultiApp System, and provided assessment of new capabilities. Noteworthy are the user assessment of the “vapor-only” flow model in Sockeye and development of a multiphysics startup transient in HP-MR unit cell for use as tutorial in Sockeye. The full-core GC-MR model was used for assessment of SAM for balance of plant modeling and for demonstrating the SWIFT code capability for hydrogen redistribution modeling in multiphysics transient analyses. In this process, several bugs/issues were identified and reported to developers. Finally, the assembly GC-MR model developed in FY-2022 coupling Griffin, BISON and SAM through flow blockage and rod ejection transients was published to the National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB). The Heat Pipe MicroReactor (HP-MR) concept high-fidelity multiphysics coupling of Griffin/BISON/Sockeye in load-following and heat pipe failure transients was also published on the VTB. Those submissions are enabling thorough review of these models as well as wide distribution to industry, regulator, and university users. In this analysis, several new research questions were uncovered, and follow-up analyses are recommended to further improve some models, consider additional transients, and continue development of VTB models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Process-dependent anisotropic thermal conductivity of laser powder bed fusion AlSi10Mg: impact of microstructure and aluminum-silicon interfaces

Purpose AlSi10Mg alloy is commonly used in laser powder bed fusion due to its printability, relatively high thermal conductivity, low density and good mechanical properties. However, the thermal conductivity of as-built materials as a function of processing (energy density, laser power, laser scanning speed, support structure) and build orientation, are not well explored in the literature. This study aims to elucidate the relationship between processing, microstructure, and thermal conductivity. Design/methodology/approach The thermal conductivity of laser powder bed fusion (L-PBF) AlSi10Mg samples are investigated by the flash diffusivity and frequency domain thermoreflectance (FDTR) techniques. Thermal conductivities are linked to the microstructure of L-PBF AlSi10Mg, which changes with processing conditions. The through-plane exceeded the in-plane thermal conductivity for all energy densities. A co-located thermal conductivity map by frequency domain thermoreflectance (FDTR) and crystallographic grain orientation map by electron backscattered diffraction (EBSD) was used to investigate the effect of microstructure on thermal conductivity. Findings The highest through-plane thermal conductivity (136 ± 2 W/m-K) was achieved at 59 J/mm 3 and exceeded the values reported previously. The in-plane thermal conductivity peaked at 117 ± 2 W/m-K at 50 J/mm 3 . The trend of thermal conductivity reducing with energy density at similar porosity was primarily due to the reduced grain size producing more Al-Si interfaces that pose thermal resistance. At these interfaces, thermal energy must convert from electrons in the aluminum to phonons in the silicon. The co-located thermal conductivity and crystallographic grain orientation maps confirmed that larger colonies of columnar grains have higher thermal conductivity compared to smaller columnar grains. Practical implications The thermal properties of AlSi10Mg are crucial to heat transfer applications including additively manufactured heatsinks, cold plates, vapor chambers, heat pipes, enclosures and heat exchangers. Additionally, thermal-based nondestructive testing methods require these properties for applications such as defect detection and simulation of L-PBF processes. Industrial standards for L-PBF processes and components can use the data for thermal applications. Originality/value To the best of the authors’ knowledge, this paper is the first to make coupled thermal conductivity maps that were matched to microstructure for L-PBF AlSi10Mg aluminum alloy. This was achieved by a unique in-house thermal conductivity mapping setup and relating the data to local SEM EBSD maps. This provides the first conclusive proof that larger grain sizes can achieve higher thermal conductivity for this processing method and material system. This study also shows that control of the solidification can result in higher thermal conductivity. It was also the first to find that the build substrate (with or without support) has a large effect on thermal conductivity.

Engineering↗

Technology test bed engine real-time failure control

The Real-Time Failure Control (RTFC) program involves development of a failure detection algorithm, for the Space Shuttle Main Engine (SSME). This failure detection approach is signal-based and entails monitoring SSME measurement signals based on predetermined as well as on-line computed mean and standard deviation values. Twenty-four engine measurements are monitored in the algorithm and provisions are made to add more parameters if needed. Each of the first values of every measurement signal at the algorithm start is checked against safety limits placed around a pre-computed engine-to-engine mean value (MV) with a bandwidth equal to a given multiple of the pre-computed standard deviation (SD). If several parameters are out of the bounds of these limits a failure is signaled. During the first two seconds (after algorithm start) a moving average (MA) and a SD is computed on-line in real-time. The moving average of each parameter is computed by averaging the incoming signal measurement with the four most recent previous signal measurements. The moving average is updated at every sampling interval (40 msec) and is checked against a similar safety band around the initial signal value for each parameter. If several anomalies are registered, a failure is signaled by the algorithm. At the end of the two-second interval the MA is fixed as the mean value for the rest of the algorithm operation and a safety band is placed above and below this value equal to a multiple of the computed SD. However, the safety band is adjusted by adjusting the mean value when propellant tank repressurization and venting take place. 'Influence Coefficients' are used to make the necessary adjustments to the safety limits of those parameters that are affected by repressurization and venting or valve closure and opening. The MA is, in both cases, continuously updated and checked against the safety band. Once more, if several parameters exceed the limits a failure is signaled. At the start of every scheduled power transient the algorithm is stopped. It is re-initiated after two seconds from the termination of the power transient and the process is repeated. The final report is divided into four major sections. The most encompassing of all is the discussion section that has sub-sections on: (1) RTFC algorithm development, (2) RTFC simulations; (3) RTFC current limitations; and (4) enhancements planned for.

Panossian, Hagop V.↗

Tone Noise and Nearfield Pressure Produced by Jet-Cavity Interaction

Cavity flow resonance can cause numerous problems in aerospace applications. While our long-term goal is to understand cavity flows well enough to devise effective cavity resonance suppression techniques, this paper describes a fundamental study of resonant tones produced by jet-cavity interaction at subsonic and supersonic speeds. Our specific jet-cavity configuration can also be used as a test bed for evaluating active and passive flow resonance control concepts. Two significant findings emerge from this study. 1) Originally, we expected that tones produced by jet-cavity interaction would resemble cavity tones or jet tones or would involve some simple combinations of each. The experimental data do not support these expectations: instead, the jet cavity interaction produce a unique set of tones. We propose simple yet and physically insightful correlations for these tones. Although the pressure patterns on the cavity floor display very complex variations with the Mach number for a length/depth = 8 cavity, the tones correspond to the acoustic modes of the cavity-independent of flow. For a length/ depth = 3 cavity, however, a surprise emerges: the pressure patterns on the cavity floor are not so complex but the tones depend significantly on the flow. Additionally, we examine the role of external feedback unique to jet-cavity interaction. 2) Previous research led us to expect that traditional classifications (open, transitional, or closed) for cavities in an infinite flight stream would be insensitive to small changes in Mach number and would depend primarily on cavity length/depth ratios. Use of the novel high resolution photoluminescent pressure sensitive paint shows that the classifications are actually quite sensitive to jet Mach number for a length/depth = 8 cavity. However, these classifications provide no guidance whatsoever for tone amplitude or frequency. Detailed experimental data and insights presented here will assist researchers who are performing numerical simulations of jet-cavity flows as a first step toward devising resonance suppression methods.

Raman, Ganesh↗

Chapter 2: Evaluating a Concentrating Solar Power Plantas an Extended-Duration Peaking Resource

We explore the ability of a concentrating solar power (CSP) plant with thermal energy storage (TES) to provide peaking capacity. We focus on future power systems, wherein net load patterns may be significantly different than they are today (e.g., due to higher renewable-energy penetrations). We examine 28 locations in the southwestern United States over an 18-year period. The hourly operation of the CSP plants are simulated to determine their potential to provide energy during an eight-hour peak-load window for each day up to 365 days per year. Our result shows that for the large majority of locations and years, CSP plants with certain configurations (i.e., in terms of solar field and TES sizes) can provide nearly 100% peak-load capacity. We examine also the amount of supplemental energy (e.g.,by using natural gas as a supplemental thermal-energy source) that would be required to ensure that a CSP plant could serve the eight highest-load hours of every day of the year. We find that in most cases, a CSP plant supplemented with natural gas would require less than 5% of the fuel that is used by a natural-gas fired power plant providing the same level of reliable capacity. A series of sensitivity analyses show that these results are robust to the number of peak-load hours and days that are considered and the configuration of the CSP plant.

capacity value↗

Isotherm Modeling and Techno-Economic Analysis of Contactor Technologies for New Tetraamine-Appended MOF for NGCC Applications

Recently, a family of tetraamine-functionalized metal-organic frameworks (MOF) has been reported as promising sorbent materials for capturing CO2 from flue gas conditions relevant to natural gas combined cycle (NGCC) applications. The main advantages of these materials are their two-step cooperative CO2 adsorption, which gives rise to unusual two step-shaped CO2 adsorption profiles and their high thermal stability. This work presents the modelling of the two-transition isotherm of the tetraamine-appended MOF, N,N'-bis(3-aminopropyl)-1,4-diaminobutane (3-4-3)-appended Mg2(dobpdc), and the techno-economic analysis (TEA) of carbon capture processes utilizing this sorbent. Due to the unusual isotherm shapes of the experimental CO2 adsorption data for tetraamine-appended Mg2(dobpdc) and the strong nonlinearity of CO2 loading with respect to temperature and pressure, we tested two different models which use logistic functions for representing the different isotherm behaviors in the different pressure ranges. The first model uses the quadratic isotherm model in the low-pressure region, the Langmuir isotherm model in the middle pressure range, and the dual site Langmuir isotherm model in the high-pressure range. To model the transition between regions we used the arctangent functions independent of temperature and the thermal effect was accounted by using the Clausius-Clapeyron relation]. The second model is an extension of the weighted dual-site Langmuir isotherm model presented by Hughes et al.. In this extended model, the dual-site Langmuir isotherm is employed in the three transition regions, using temperature-dependent logistic functions to activate or deactivate the isotherm model in the low, middle, and high-pressure ranges. Both models fit the experimental data quite well with root mean squared errors (RMSE’s) of 0.41 and 0.17 for model 1 and model 2, respectively. Since model 2 resulted in a lower RMSE, it was leveraged for the development of the gas/solid contactor models used by the TEA. Specifically, two different contactor models, an axial-flow fixed bed and moving bed contactor, were developed as part of this work. These models are dynamic, pressure-driven, and consist of mass, energy, and momentum conservation equations. A kinetic model was also developed by performing parameter estimation using experimental fixed bed breakthrough data. These models are then used to simulate CO2 capture processes from the flue gas generated from a ~600 gross MW NGCC power plant. A cost model was developed which considers the capital cost of the reactors and the significant operating costs such as steam and electricity. Using NETL’s Framework for the Optimization and Quantification of Uncertainty of Uncertainty and Surrogates tool (FOQUS), which has the capability of linking models built using numerous modelling platforms with derivative-free optimization solvers, a techno-economic optimization of the carbon capture processes was performed which minimizes the cost of capture.

Caballero, Daison↗

Online Monitoring of Medium Voltage Cable Systems with Spread Spectrum Time Domain and Frequency Domain Reflectometry

In-service failures of wave energy convertor (WEC) cable systems can have a significant cost and power availability impact. Close parallel research 2019 data showed > 1B£ and 9 Terra-Watt-Hours associated with global off-shore wind (OSW) cable failures (Strang-Moran 2020). OSW is a closely related technology but currently is significantly cheaper than WEC technology. For wave energy to compete, the problem of reliable cable transmission must be mitigated. This project develops isolation technology to allow online high frequency reflectometry testing of medium voltage cables (1 to 10 kV and higher) without arcing or damage to the test instrument. Online spread spectrum time domain reflectometry (SSTDR) testing has been established for low voltage cable systems in the aircraft and rail industry and the ability to detect and locate cable flaws of interest is well understood. Extending reflectometry testing to medium voltage systems could enable detection of cable damage before failures occur thereby allowing repair and replacement of damaged cable segments to be scheduled and managed. The seedling project succeeded to pass and receive high frequency SSTDR signals onto a cable up to 1 kV using a parallel trace isolation circuit board that can be connected onto the test cable. The approach used a novel circuit design for which an invention disclosure has been filed. A proposed sapling project would extend the technology toward the higher operating voltages used by WEC systems, thereby enabling online SSTDR cable monitoring. The goal of the seedling project was to extend the capability of the ARENA cable/motor test bed to address medium voltages and to develop a high pass filter isolation architecture to protect the reflectometry instrument from the low frequency (DC – 60 Hz) line voltage while allowing the high frequency diagnostic signal to pass to and from the test instrument to the live line. Initial efforts focused on passive LCR filter circuits to reduce 60 Hz levels below 10 volts from a 10 kV line while allowing the MHz high frequency chirps to pass onto the cables and for mV signals to be detected. We discovered that the parasitic loss behavior of real high voltage components precluded this approach from working. An alternate approach was adapted for the electric field to couple between two parallel traces on a printed circuit board much like a radio-frequency coupler. The challenge here was and is to have the parallel traces close enough to each other to effectively pass the high frequency chirp onto the live line and receive any reflected signal from any encountered impedance change along the cable. This reflected signal will be in the mV range. The traces however must be far enough apart to not allow arcing on the board. A design with 3 mm spacing was determined to allow the high frequency signal to pass onto the live line and receive the mV signal back into the instrument while reducing the 60 Hz voltage amplitude by >80 dB (more than a factor of 10,000) without allowing arcing from across the parallel traces. This was confirmed by simulation and test.

16 TIDAL AND WAVE POWER↗

NASA Stennis Space Center Integrated System Health Management Test Bed and Development Capabilities

Integrated System Health Management (ISHM) is a capability that focuses on determining the condition (health) of every element in a complex System (detect anomalies, diagnose causes, prognosis of future anomalies), and provide data, information, and knowledge (DIaK)-not just data-to control systems for safe and effective operation. This capability is currently done by large teams of people, primarily from ground, but needs to be embedded on-board systems to a higher degree to enable NASA's new Exploration Mission (long term travel and stay in space), while increasing safety and decreasing life cycle costs of spacecraft (vehicles; platforms; bases or outposts; and ground test, launch, and processing operations). The topics related to this capability include: 1) ISHM Related News Articles; 2) ISHM Vision For Exploration; 3) Layers Representing How ISHM is Currently Performed; 4) ISHM Testbeds & Prototypes at NASA SSC; 5) ISHM Functional Capability Level (FCL); 6) ISHM Functional Capability Level (FCL) and Technology Readiness Level (TRL); 7) Core Elements: Capabilities Needed; 8) Core Elements; 9) Open Systems Architecture for Condition-Based Maintenance (OSA-CBM); 10) Core Elements: Architecture, taxonomy, and ontology (ATO) for DIaK management; 11) Core Elements: ATO for DIaK Management; 12) ISHM Architecture Physical Implementation; 13) Core Elements: Standards; 14) Systematic Implementation; 15) Sketch of Work Phasing; 16) Interrelationship Between Traditional Avionics Systems, Time Critical ISHM and Advanced ISHM; 17) Testbeds and On-Board ISHM; 18) Testbed Requirements: RETS AND ISS; 19) Sustainable Development and Validation Process; 20) Development of on-board ISHM; 21) Taxonomy/Ontology of Object Oriented Implementation; 22) ISHM Capability on the E1 Test Stand Hydraulic System; 23) Define Relationships to Embed Intelligence; 24) Intelligent Elements Physical and Virtual; 25) ISHM Testbeds and Prototypes at SSC Current Implementations; 26) Trailer-Mounted RETS; 27) Modeling and Simulation; 28) Summary ISHM Testbed Environments; 29) Data Mining - ARC; 30) Transitioning ISHM to Support NASA Missions; 31) Feature Detection Routines; 32) Sample Features Detected in SSC Test Stand Data; and 33) Health Assessment Database (DIaK Repository).

Figueroa, Fernando↗

Large Scale Testing of a Foam/Multilayer Insulation Thermal Control System (TCS) for Cryogenic Upper Stages

The development of high energy cryogenic upper stages is essential for the efficient delivery of large payloads to various destinations envisioned in future programs. A key element in such upper stages is cryogenic fluid management (CFM) advanced development/technology. Due to the cost of and limited opportunities for orbital experiments, ground testing must be employed to the fullest extent possible. Therefore, a system level test bed termed the Multipurpose Hydrogen Test Bed (MHTB), which is representative in size and shape (3 meter diameter by 3 meter long with a volume of 18 cubic meters) of a fully integrated space transportation vehicle liquid hydrogen propellant tank has been established. To date, upper stage studies have often baselined the foam/multilayer insulation (FMLI) combination concept; however, hardware experience with the concept is minimal and was therefore selected for the MHTB. The foam element (isofoam SS-1 171 with an average thickness of 3.5 centimeters) is designed to protect against ground hold/ascent flight environments, and allows for the use of a dry nitrogen purge as opposed to the more complex/heavy helium purge subsystem normally required with MLI in cryogenic applications. The MLI (45 layers of Double Aluminized Mylar with Dacron spacers) provides protection in the vacuum environment of space and is designed for an on-orbit storage period of 45 days. Several unique features were incorporated in the MLI concept and included: variable density MLI (reduces weight and radiation losses by changing the layer density), larger but fewer DAM perforations for venting during ascent to orbit (reduces radiation losses), and roll wrap installation of the MLI with a commercially established process to lower assembly man-hours and reduce seam heat leak. Thermal performance testing of the MHTB TCS was conducted during three test series conducted between September 1995 and May 1996. Results for the ground hold portion of the tests were as expected producing an average heat leak of 63 WattS/M2 at an average foam surface temperature of 170 K. The results of the simulated orbit hold test interval produced heat leaks ranging from 0.085 to 0.22 Watts/squareM at warm boundary temperatures of 164K and 305K, respectively. When compared to the performance for a traditional MLI system, a 60% reduction in orbital heat leak or boiloff was measured. Overall, the MHTB TCS demonstrated satisfactory performance for all mission phases required of a cryogenic upper stage.

Hastings, Leon↗

Optics Program Modified for Multithreaded Parallel Computing

A powerful high-performance computer program for simulating and analyzing adaptive and controlled optical systems has been developed by modifying the serial version of the Modeling and Analysis for Controlled Optical Systems (MACOS) program to impart capabilities for multithreaded parallel processing on computing systems ranging from supercomputers down to Symmetric Multiprocessing (SMP) personal computers. The modifications included the incorporation of OpenMP, a portable and widely supported application interface software, that can be used to explicitly add multithreaded parallelism to an application program under a shared-memory programming model. OpenMP was applied to parallelize ray-tracing calculations, one of the major computing components in MACOS. Multithreading is also used in the diffraction propagation of light in MACOS based on pthreads [POSIX Thread, (where "POSIX" signifies a portable operating system for UNIX)]. In tests of the parallelized version of MACOS, the speedup in ray-tracing calculations was found to be linear, or proportional to the number of processors, while the speedup in diffraction calculations ranged from 50 to 60 percent, depending on the type and number of processors. The parallelized version of MACOS is portable, and, to the user, its interface is basically the same as that of the original serial version of MACOS.

Lou, John↗

Performance of Pulsed Thermal Tomography Imaging with Machine Learning-Based Classification of Defects in Additively Manufactured Structures

Additive manufacturing (AM) is an emerging method for cost-efficient fabrication of complex topology nuclear reactor parts from high-strength corrosion resistance alloys, such as stainless steel and Inconel. AM of metallic structures for nuclear energy applications is currently based on laser powder bed fusion (LPBF) process. Some of the challenges with using LPBF method for nuclear manufacturing include the possibility of introducing pores into metallic structures. Integrity of AM structures needs to be evaluated nondestructively because material flaws could lead to premature failures in high temperature nuclear reactor environment. Currently, there exist limited capabilities to evaluate actual AM structures non-destructively. Pulsed Thermal Tomography Imaging (PTT) provides a capability for non-destructive evaluation (NDE) of subsurface defects in arbitrary size structures. The PTT method is based on recording material surface temperature transients with infrared (IR) camera following thermal pulse delivered on material surface with flash light. The PTT method has advantages for NDE of actual AM structures because the method involves one-sided non-contact measurements and fast processing of large sample areas captured in one image. Following initial qualification of an AM component for deployment in a nuclear reactor, a PTT system can also be used for in-service nondestructive evaluation (NDE) applications. In this report, we describe recent progress in enhancing PTT capabilities in detecting and visualizing microscopic defects in metallic specimens. The thermal tomography (TT) algorithm obtains depth reconstructions of spatial effusivity from the data cube of sequentially recorded surface temperatures. However, interpretation of TT images is non-trivial because of blurring of images with increasing depth. To address this challenge, we have developed a deep learning convolutional neural network (CNN) to classify size and orientation subsurface defects in simulated TT images. CNN is trained on a database of TT images created for a set of simulated metallic structures with elliptical subsurface voids. Test of CNN performance demonstrate the ability to classify radii and angular orientation of subsurface defects in TT images. In addition, we have shown that CNN trained on elliptical defects is capable of classifying irregular-shaped defects obtained from scanning electron microscopy (SEM) of stainless steel sections printed with LPBF.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Streamlined Approach to Spatial Mapping of Complex 3D Thermal Boundary Condition Data

Spatial mapping of 1D, 2D, and particularly 3D thermal boundary condition data is essential for accurate thermal simulation of spacecraft and launch vehicles. This information typically is available in simple ASCII file format but can be problematic to use for thermal simulations. Difficulties that can be encountered include the ability to directly use the available file format, transformation of the data into the desired reference coordinate system in 3D space, and accurate interpolation of the data onto the thermal calculation domain. This paper presents a streamlined approach to spatial mapping of complex unsymmetrical 3D cold wall heat flux data using the fields’ capability in Simcenter3D. A recent example of 3D thermal plume simulation at MSFC will be employed to illustrate the use of text, csv, and Excel files, specifying coordinate systems for transformation of data in 3D space, and use of the various interpolation schemes available for spatial mapping. Spatial boundary condition verification is also very important, and the spacecraft thermal analyst needs effective visualization tools to develop confidence in the boundary condition definitions. Visualization of the spacecraft plume raw field data, the calculated field data, and the interpolated field data onto the thermal calculation domain will be discussed, with a particular focus on real-time viewing of the interpolated data for available interpolation schemes. The objective of this paper is to outline the modeling techniques used with the integrated thermal analysis of a terrestrial lunar lander demonstrator. Specifically, the techniques outlined herein focus on spatial mapping of convective and radiative exhaust plume results from a computational fluid dynamics (CFD) analysis. Two cases were investigated, however, the modeling techniques that were utilized translate to both cases seamlessly. Those two cases were as follows: • Lunar lander demonstrator firing on-pad (worst case). • Lunar lander demonstrator hovering at an elevated position above the launch pad. For the purposes of information control, this paper will focus on the modeling techniques only, using altered or dimensionless data where necessary. BACKGROUND The thermal model examined in this paper is of the lunar lander demonstrator known as XL-1T (terrestrial), born from a collaborative effort between Masten Space Systems (MSS) and NASA, under NASA’s Lunar CATALYST (also known as Lunar Cargo Transportation and Landing by Soft Touchdown) initiative. The lander is a reusable vertical takeoff/vertical landing (VTVL) test bed which is controlled by four throttleable main engines utilizing green hypergolic propellants.

Hawkins, Robert F.↗

Slip Tendency Analysis of Fracture Networks to Determine Suitability of Candidate Testbeds for the EGS Collab Hydroshear Experiment

Experiment 2 of the EGS Collab project is aimed at testing stimulation by hydro-shearing of existing natural fractures, versus Experiment 1, which is focused on hydraulic fracturing a rock mass to enhance permeability. The main criterion for the testbed selection in Experiment 2 is the presence of an interconnected network of fractures, at least 10 meters in extent, and with orientations such that shear slip can be induced at injection pressures less than the minimum horizontal stress (Shmin). The feasibility analysis for this experiment requires a well-constrained stress state along with well-characterized fracture networks. The fracture systems at two candidate locations in the Sanford Underground Research Facility (SURF), the 4850 Level (number refers to depth below ground surface in feet) and the 4100 Level, have been characterized to different extents. The 4850 Level has a well-characterized Discrete Fracture Network (DFN) from borehole (drilled for an experimental test bed) and drift observations. In contrast, the 4100 level has a fracture network characterized only by observations from the drift wall, which provides little constraint on fracture extents. This paper will present assessments for the interpreted slip potential for the natural fractures at the two locations. Data uncertainties are addressed by performing a probabilistic analysis that takes into account the uncertainty in the stress state, uncertainty in the fracture properties, and preliminary borehole locations. There is strong evidence of shear stimulation of some natural fractures in Experiment 1 on the 4850 level. We use examples of natural fractures strongly linked to shear stimulation in Experiment 1 from multiple monitoring indicators to test the consistency of the stress model and guide the uncertainty interpretation. This analysis is intended to guide the site selection process for Experiment 2 by highlighting the fracture orientations that are likely to be shear-stimulated in a majority of modeled realizations.

40 EE - Geothermal Technologies Office (EE-4G)↗

An Orthogonal Recursive Bisection (ORB) Based Time Advancement Algorithm for CFD-DEM Solvers

The time integration of the granular phase in coupled computational fluid dynamics (CFD) – discrete element method (DEM) simulations presents a unique computational challenge brought about by the large variations in particle collisional time scales. Particles in the dilute regions of the computational domain can be advanced with large time steps while dense regions require much smaller time increments. However, the time step size in most solvers is globally set as the limit for accuracy and stability imposed by the collisions and is typically orders of magnitude less than that required away from collisions. This work addresses this precise issue and provides a strategy to avoid the use of a global conservative small time step size for the entire set of particles.A novel time stepping algorithm for CFD-DEM solvers using a partitioning approach using orthogonal recursive bisection (ORB) that allows for variable time steps among particles is described and its computational performance is compared against baseline explicit methods, typically used in several CFD-DEM solvers. ORB has advantages of being relatively quick and easy to update incrementally and has the required heuristic behavior (i.e., it will split the region in half with a cluster on each side) when groups of particles are well separated (clustered). The algorithm presented in this work uses a local time stepping approach to resolve collisional time scales for subsets of particles that are present at the leaves of the ORB, thereby resulting in substantial reduction of computational cost. The parallel implementation of this method where a ``knapsack” algorithm is used in tandem with ORB for effective load-balancing is also presented, where a best possible partitioning is obtained based on number of particles and local time-stepping costs. The algorithm is tested against benchmark problems with varying particle distributions that include fluidized bed and riser flow scenarios. Preliminary results indicate that the approach is 2-3X faster than traditional explicit methods for problems that involve both dense and dilute regions, while maintaining the same level of accuracy.

adaptive timestepping↗