Search NASASearch

SEARCH · Search NASA

Results for “incremental dynamic analysis”

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.

108 records · Page 6

Characterization of Low-Velocity Impact Damage in Toughened Non-Crimp Fabric Laminates

The damage response of resin-infused composites was investigated using quasi-static indentation (QSI) and instrumented drop-weight impact testing. The test specimens were composed of six layers of non-crimp fabrics where each layer contained +45°, 90°, −45°, and 0° plies. The goals of the testing were to identify damage responses that may differ from traditional thermoset materials and provide damage characterization data to evaluate current damage prediction analyses and aid improvements to current damage models. The damage was characterized using X-ray computed tomography (CT) and ultrasonic testing (UT). Incremental QSI loading was used on a single specimen to investigate damage growth in nine steps at progressively higher loads. Realtime acoustic emission data were used to determine the unloading points. Unlike typical thermoset composites, the force-displacement was smooth and showed no obvious indications of damage. The primary damage consisted of large lines of fiber fracture on the plies near the impact surface and delaminations through the thickness. Nine specimens were impacted at nominal impact energies of 25 ft-lbs, 30 ft-lbs, and 35 ft-lbs. The force histories at each energy followed nearly identical paths. Again, no obvious signs of damage were observed in any of the force histories. The general damage response in the dynamic tests was similar to the QSI-loaded specimens. However, repeat specimens impacted with a given impact energy showed variations in the damage. Some characteristics of the damage response were different than those typically observed in thermoset composites and will need to be evaluated using current analysis methods.

Composite Materials

Use Computer-Aided Tools to Parallelize Large CFD Applications

Porting applications to high performance parallel computers is always a challenging task. It is time consuming and costly. With rapid progressing in hardware architectures and increasing complexity of real applications in recent years, the problem becomes even more sever. Today, scalability and high performance are mostly involving handwritten parallel programs using message-passing libraries (e.g. MPI). However, this process is very difficult and often error-prone. The recent reemergence of shared memory parallel (SMP) architectures, such as the cache coherent Non-Uniform Memory Access (ccNUMA) architecture used in the SGI Origin 2000, show good prospects for scaling beyond hundreds of processors. Programming on an SMP is simplified by working in a globally accessible address space. The user can supply compiler directives, such as OpenMP, to parallelize the code. As an industry standard for portable implementation of parallel programs for SMPs, OpenMP is a set of compiler directives and callable runtime library routines that extend Fortran, C and C++ to express shared memory parallelism. It promises an incremental path for parallel conversion of existing software, as well as scalability and performance for a complete rewrite or an entirely new development. Perhaps the main disadvantage of programming with directives is that inserted directives may not necessarily enhance performance. In the worst cases, it can create erroneous results. While vendors have provided tools to perform error-checking and profiling, automation in directive insertion is very limited and often failed on large programs, primarily due to the lack of a thorough enough data dependence analysis. To overcome the deficiency, we have developed a toolkit, CAPO, to automatically insert OpenMP directives in Fortran programs and apply certain degrees of optimization. CAPO is aimed at taking advantage of detailed inter-procedural dependence analysis provided by CAPTools, developed by the University of Greenwich, to reduce potential errors made by users. Earlier tests on NAS Benchmarks and ARC3D have demonstrated good success of this tool. In this study, we have applied CAPO to parallelize three large applications in the area of computational fluid dynamics (CFD): OVERFLOW, TLNS3D and INS3D. These codes are widely used for solving Navier-Stokes equations with complicated boundary conditions and turbulence model in multiple zones. Each one comprises of from 50K to 1,00k lines of FORTRAN77. As an example, CAPO took 77 hours to complete the data dependence analysis of OVERFLOW on a workstation (SGI, 175MHz, R10K processor). A fair amount of effort was spent on correcting false dependencies due to lack of necessary knowledge during the analysis. Even so, CAPO provides an easy way for user to interact with the parallelization process. The OpenMP version was generated within a day after the analysis was completed. Due to sequential algorithms involved, code sections in TLNS3D and INS3D need to be restructured by hand to produce more efficient parallel codes. An included figure shows preliminary test results of the generated OVERFLOW with several test cases in single zone. The MPI data points for the small test case were taken from a handcoded MPI version. As we can see, CAPO's version has achieved 18 fold speed up on 32 nodes of the SGI O2K. For the small test case, it outperformed the MPI version. These results are very encouraging, but further work is needed. For example, although CAPO attempts to place directives on the outer- most parallel loops in an interprocedural framework, it does not insert directives based on the best manual strategy. In particular, it lacks the support of parallelization at the multi-zone level. Future work will emphasize on the development of methodology to work in a multi-zone level and with a hybrid approach. Development of tools to perform more complicated code transformation is also needed.

Jin, H.

Monitoring the Microgravity Environment Quality On-board the International Space Station Using Soft Computing Techniques: Preliminary System Performance Results - Part 2

This paper presents the preliminary performance results of the artificial intelligence monitoring system in full operational mode using near real time acceleration data downlinked from the International Space Station. Preliminary microgravity environment characterization analysis result for the International Space Station (Increment-2), using the monitoring system is presented. Also, comparison between the system predicted performance based on ground test data for the US laboratory "Destiny" module and actual on-orbit performance, using measured acceleration data from the U.S. laboratory module of the International Space Station is presented. Finally, preliminary on-orbit disturbance magnitude levels are presented for the Experiment of Physics of Colloids in Space, which are compared with on ground test data. The ground test data for the Experiment of Physics of Colloids in Space were acquired from the Microgravity Emission Laboratory, located at the NASA Glenn Research Center, Cleveland, Ohio. The artificial intelligence was developed by the NASA Glenn Principal Investigator Microgravity Services Project to help the principal investigator teams identify the primary vibratory disturbance sources that are active, at any moment of time, on-board the International Space Station, which might impact the microgravity environment their experiments are exposed to. From the Principal Investigator Microgravity Services' web site, the principal investigator teams can monitor via a dynamic graphical display, implemented in Java, in near real time, which event(s) is/are on, such as crew activities, pumps, fans, centrifuges, compressor, crew exercise, structural modes, etc., and decide whether or not to run their experiments, whenever that is an option, based on the acceleration magnitude and frequency sensitivity associated with that experiment. This monitoring system detects primarily the vibratory disturbance sources. The system has built-in capability to detect both known and unknown vibratory disturbance sources. Several soft computing techniques such as Kohonen's Self-Organizing Feature Map, Learning Vector Quantization, Back-Propagation Neural Networks, and Fuzzy Logic were used to design the system.

Jules, Kenol

Biomolecule Sequencer: Next-Generation DNA Sequencing Technology for In-Flight Environmental Monitoring, Research, and Beyond

On the International Space Station (ISS), technologies capable of rapid microbial identification and disease diagnostics are not currently available. NASA still relies upon sample return for comprehensive, molecular-based sample characterization. Next-generation DNA sequencing is a powerful approach for identifying microorganisms in air, water, and surfaces onboard spacecraft. The Biomolecule Sequencer payload, manifested to SpaceX-9 and scheduled on the Increment 4748 research plan (June 2016), will assess the functionality of a commercially-available next-generation DNA sequencer in the microgravity environment of ISS. The MinION device from Oxford Nanopore Technologies (Oxford, UK) measures picoamp changes in electrical current dependent on nucleotide sequences of the DNA strand migrating through nanopores in the system. The hardware is exceptionally small (9.5 x 3.2 x 1.6 cm), lightweight (120 grams), and powered only by a USB connection. For the ISS technology demonstration, the Biomolecule Sequencer will be powered by a Microsoft Surface Pro3. Ground-prepared samples containing lambda bacteriophage, Escherichia coli, and mouse genomic DNA, will be launched and stored frozen on the ISS until experiment initiation. Immediately prior to sequencing, a crew member will collect and thaw frozen DNA samples, connect the sequencer to the Surface Pro3, inject thawed samples into a MinION flow cell, and initiate sequencing. At the completion of the sequencing run, data will be downlinked for ground analysis. Identical, synchronous ground controls will be used for data comparisons to determine sequencer functionality, run-time sequence, current dynamics, and overall accuracy. We will present our latest results from the ISS flight experiment the first time DNA has ever been sequenced in space and discuss the many potential applications of the Biomolecule Sequencer for environmental monitoring, medical diagnostics, higher fidelity and more adaptable Space Biology Human Research Program investigations, and even life detection experiments for astrobiology missions.

Sequencer

Imaging spectrophotometry of the nuclear outflow of NGC 1068

This observational program (in conjunction with R. B. Tully (IfA, Honolulu), and J. Bland (Rice U., Houston)) aims to constrain the kinematic organization and dominant excitation mechanisms of ionized gas in active galaxies. More generally, researchers are interested in the dynamics of radiative, supersonic flows in the Interstellar Medium (ISM). Imaging Fabry-Perot interferometers and low-noise Charge Coupled Devices (CCDs) are used for complete spatial coverage of the complex gas distribution in circumnuclear narrow-line regions (NLRs). Extranuclear emission line widths in NLRs can exceed 3000 km s(-1), so to avoid inter-order confusion researchers use an etalon of 4000 km s(-1) free spectral range to map (N II) lambda lambda 6548, 6583 and H alpha. To maximize spatial resolution, researchers select nearby active systems independent of luminosity but known to possess interesting morphologies and/or high-velocity extranuclear ionized gas. Monochromatic images Full Width Half Maximum (FWHM) approx. 65 km s(-1) have thus far been obtained in 1 second or better seeing at the U. Hawaii 2.2m, CFH 3.6m, and CTIO 4.0m telescopes. These are stacked into grids of line profiles, of spectrophotometric quality, at sub-arcsecond increments across a 3 second field. To handle the approx. 20,000 to 300,000 useful spectra that arise from a typical night's work, researchers have developed a complete analysis and reduction package for VAX and Sun image workstations. Reduction involves parametrization of approx. 10 to the 8th raw data points to a few maps (e.g., velocities of each kinematic subsystem, continuum-free line fluxes) containing approx. 10 to the 5th pixels. Researchers identify kinematic and structural symmetries by examining these maps and the point to point variations of the synthesized line profiles. The combination of monochromatic images and full spatial sampling of line profiles has allowed them to isolate such symmetries and has led to reliable kinematic deprojections.

Cecil, Gerald

File Specification for GEOS-5 FP (Forward Processing)

The GEOS-5 FP Atmospheric Data Assimilation System (GEOS-5 ADAS) uses an analysis developed jointly with NOAA's National Centers for Environmental Prediction (NCEP), which allows the Global Modeling and Assimilation Office (GMAO) to take advantage of the developments at NCEP and the Joint Center for Satellite Data Assimilation (JCSDA). The GEOS-5 AGCM uses the finite-volume dynamics (Lin, 2004) integrated with various physics packages (e.g, Bacmeister et al., 2006), under the Earth System Modeling Framework (ESMF) including the Catchment Land Surface Model (CLSM) (e.g., Koster et al., 2000). The GSI analysis is a three-dimensional variational (3DVar) analysis applied in grid-point space to facilitate the implementation of anisotropic, inhomogeneous covariances (e.g., Wu et al., 2002; Derber et al., 2003). The GSI implementation for GEOS-5 FP incorporates a set of recursive filters that produce approximately Gaussian smoothing kernels and isotropic correlation functions. The GEOS-5 ADAS is documented in Rienecker et al. (2008). More recent updates to the model are presented in Molod et al. (2011). The GEOS-5 system actively assimilates roughly 2 × 10(exp 6) observations for each analysis, including about 7.5 × 10(exp 5) AIRS radiance data. The input stream is roughly twice this volume, but because of the large volume, the data are thinned commensurate with the analysis grid to reduce the computational burden. Data are also rejected from the analysis through quality control procedures designed to detect, for example, the presence of cloud. To minimize the spurious periodic perturbations of the analysis, GEOS-5 FP uses the Incremental Analysis Update (IAU) technique developed by Bloom et al. (1996). More details of this procedure are given in Appendix A. The assimilation is performed at a horizontal resolution of 0.3125-degree longitude by 0.25- degree latitude and at 72 levels, extending to 0.01 hPa. All products are generated at the native resolution of the horizontal grid. The majority of data products are time-averaged, but four instantaneous products are also available. Hourly data intervals are used for two-dimensional products, while 3-hourly intervals are used for three-dimensional products. These may be on the model's native 72-layer vertical grid or at 42 pressure surfaces extending to 0.1 hPa. This document describes the gridded output files produced by the GMAO near real-time operational FP, using the most recent version of the GEOS-5 assimilation system. Additional details about variables listed in this file specification can be found in a separate document, the GEOS-5 File Specification Variable Definition Glossary. Documentation about the current access methods for products described in this document can be found on the GMAO products page: http://gmao.gsfc.nasa.gov/products/.

GSI

Analysis of Piezoelectric Structural Sensors with Emergent Computing Techniques

The purpose of this project was to try to interpret the results of some tests that were performed earlier this year and to demonstrate a possible use of emergence in computing to solve IVHM problems. The test data used was collected with piezoelectric sensors to detect mechanical changes in structures. This project team was included of Dr. Doug Ramers and Dr. Abdul Jallob of the Summer Faculty Fellowship Program, Arnaldo Colon-Lopez - a student intern from the University of Puerto Rico of Turabo, and John Lassister and Bob Engberg of the Structural and Dynamics Test Group. The tests were performed by Bob Engberg to compare the performance two types of piezoelectric (piezo) sensors, Pb(Zr(sub 1-1)Ti(sub x))O3, which we will label PZT, and Pb(Zn(sub 1/3)Nb(sub 2/3))O3-PbTiO, which we will label SCP. The tests were conducted under varying temperature and pressure conditions. One set of tests was done by varying water pressure inside an aluminum liner covered with carbon-fiber composite layers (a cylindrical "bottle" with domed ends) and the other by varying temperatures down to cryogenic levels on some specially prepared composite panels. This report discusses the data from the pressure study. The study of the temperature results was not completed in time for this report. The particular sensing done with these piezo sensors is accomplished by the sensor generating an controlled vibration that is transmitted into the structure to which the sensor is attached, and the same sensor then responding to the induced vibration of the structure. There is a relationship between the mechanical impedance of the structure and the resulting electrical impedance produced in the in the piezo sensor. The impedance is also a function of the excitation frequency. Changes in the real part of impendance signature relative to an original reference signature indicate a change in the coupled structure that could be the results of damage or strain. The water pressure tests were conducted by pressurizing the bottle on a test stand, and running sweeps of excitations frequencies for each of the piezo sensors and recording the resulting impedance. The sweeps were limited to 401 points by the available analyzer, and it was decided to perform individual sweeps at five different excitation frequency ranges. The frequency ranges used for the PZTs were different in two of the five ranges from the ranges used for the SCP. The bottles were pressurized to empty (no water), 0psig, 77 psig, 155 psig, 227 psig in nearly uniform increments of about 77psi. One of each of the two types of piezo sensors was fastened on to the bottle surface at two locations: about midway between the ends on cylindrical portion of the bottle and at the very edge of one of the end domes. The data was collected in files by sensor type (2 cases), by location (2 cases), by frequency range (5 cases), and pressure (5cases) to produce 100 data sets of 401 impedances. After familiarization with the piezo sensing technology and obtaining the data, the team developed a set of questions to try to answer regarding the data and made assignments of responsibilities. The next section lists the questions, and the remainder of the report describes the data analysis work performed by Dr. Ramers. This includes a discussion of the data, the approach to answering the question using statistical techniques, the use of an emergent system to investigate the data where statistical techniques were not usable, conclusions regarding the data, and recommendations.

Ramers, Douglas L.

Structural Dynamics Observations in Space Launch System Green Run Hot Fire Testing

The Space Launch System (SLS) Core Stage (CS) Thrust Vector Control (TVC) system is comprised of eight mechanical feedback Shuttle heritage Type III TVC actuators and four RS-25 engines, each attached to a Shuttle heritage gimbal block/bearing. Two actuators are used to move each engine in two planes perpendicular to one another (i.e., pitch and yaw). The TVC system design leverages hardware from the Space Shuttle program as well as new hardware designed specifically for the Core Stage. The Green Run Hot Fire (GRHF) of the SLS Core Stage provided a flight-like ground test environment for verification of integrated vehicle TVC performance. A TVC model coupled to a vehicle structural dynamic model has been developed previously and incrementally validated in subsystem tests and simulations. Still, some aspects of TVC performance in GRHF were not anticipated. The ensuing investigation demonstrated the need for well-instrumented test environments, various levels of modeling fidelity, test-representative structural models, and caution in reuse of legacy components. This paper is the sixth installment in a seven-paper series surveying the design, engineering, test validation, and flight performance of the Core Stage Thrust Vector Control system. It introduces the salient structural dynamic phenomena uncovered in ambient and hot fire testing. During the Green Run test campaign, a comparison of ambient and hot fire step responses showed a significant change in apparent damping due to the presence of friction, challenging long standing assumptions that friction could be neglected. Additionally, the characteristic response of the engine and thrust structure during GRHF proved to be more complex than anticipated, as evidenced by the available actuator, thrust structure, and engine measurements. While the string-potentiometer based test instrumentation was intended to allow for reconstruction of the engine angles along the two control axes, the geometric placement, location uncertainty, and responses in overlapping frequency spectra revealed additional phenomena requiring further analysis and post-processing. The observations from both modal and frequency response testing during the Green Run ambient and hot fire configurations led to Engine and Core Stage FEM (finite element model) updates. When evidence of unexpected engine motion was found in engine section accelerometer data, the authors pursued additional structural analysis leading to FEM updates associated with the TVC gimbal and thrust structure. Through collaboration between structures, TVC, and flight control disciplines, the test-informed models and root-cause analysis led to confident flight rationale for the first flight of the SLS launch vehicle.

Richard K. Moore

Structural Dynamics Observations in Space Launch System Green Run Hot Fire Testing

The Space Launch System (SLS) Core Stage (CS) Thrust Vector Control (TVC) system is comprised of eight mechanical feedback Shuttle heritage Type III TVC actuators and four RS-25 engines, each attached to a Shuttle heritage gimbal block/bearing. Two actuators are used to move each engine in two planes perpendicular to one another (i.e., pitch and yaw). The TVC system design leverages hardware from the Space Shuttle program as well as new hardware designed specifically for the Core Stage. The Green Run Hot Fire (GRHF) of the SLS Core Stage provided a flight-like ground test environment for verification of integrated vehicle TVC performance. A TVC model coupled to a vehicle structural dynamic model has been developed previously and incrementally validated in subsystem tests and simulations. Still, some aspects of TVC performance in GRHF were not anticipated. The ensuing investigation demonstrated the need for well-instrumented test environments, various levels of modeling fidelity, test-representative structural models, and caution in reuse of legacy components. This paper is the sixth installment in a seven-paper series surveying the design, engineering, test validation, and flight performance of the Core Stage Thrust Vector Control system. It introduces the salient structural dynamic phenomena uncovered in ambient and hot fire testing. During the Green Run test campaign, a comparison of ambient and hot fire step responses showed a significant change in apparent damping due to the presence of friction, challenging long standing assumptions that friction could be neglected. Additionally, the characteristic response of the engine and thrust structure during GRHF proved to be more complex than anticipated, as evidenced by the available actuator, thrust structure, and engine measurements. While the string-potentiometer based test instrumentation was intended to allow for reconstruction of the engine angles along the two control axes, the geometric placement, location uncertainty, and responses in overlapping frequency spectra revealed additional phenomena requiring further analysis and post-processing. The observations from both modal and frequency response testing during the Green Run ambient and hot fire configurations led to Engine and Core Stage FEM (finite element model) updates. When evidence of unexpected engine motion was found in engine section accelerometer data, the authors pursued additional structural analysis leading to FEM updates associated with the TVC gimbal and thrust structure. Through collaboration between structures, TVC, and flight control disciplines, the test-informed models and root-cause analysis led to confident flight rationale for the first flight of the SLS launch vehicle.

Richard Moore

File Specification for GEOS-5 FP-IT (Forward Processing for Instrument Teams)

The GEOS-5 FP-IT Atmospheric Data Assimilation System (GEOS-5 ADAS) uses an analysis developed jointly with NOAA's National Centers for Environmental Prediction (NCEP), which allows the Global Modeling and Assimilation Office (GMAO) to take advantage of the developments at NCEP and the Joint Center for Satellite Data Assimilation (JCSDA). The GEOS-5 AGCM uses the finite-volume dynamics (Lin, 2004) integrated with various physics packages (e.g, Bacmeister et al., 2006), under the Earth System Modeling Framework (ESMF) including the Catchment Land Surface Model (CLSM) (e.g., Koster et al., 2000). The GSI analysis is a three-dimensional variational (3DVar) analysis applied in grid-point space to facilitate the implementation of anisotropic, inhomogeneous covariances (e.g., Wu et al., 2002; Derber et al., 2003). The GSI implementation for GEOS-5 FP-IT incorporates a set of recursive filters that produce approximately Gaussian smoothing kernels and isotropic correlation functions. The GEOS-5 ADAS is documented in Rienecker et al. (2008). More recent updates to the model are presented in Molod et al. (2011). The GEOS-5 system actively assimilates roughly 2 × 10(exp 6) observations for each analysis, including about 7.5 × 10(exp 5) AIRS radiance data. The input stream is roughly twice this volume, but because of the large volume, the data are thinned commensurate with the analysis grid to reduce the computational burden. Data are also rejected from the analysis through quality control procedures designed to detect, for example, the presence of cloud. To minimize the spurious periodic perturbations of the analysis, GEOS-5 FP-IT uses the Incremental Analysis Update (IAU) technique developed by Bloom et al. (1996). More details of this procedure are given in Appendix A. The analysis is performed at a horizontal resolution of 0.625-degree longitude by 0.5-degree latitude and at 72 levels, extending to 0.01 hPa. All products are generated at the native resolution of the horizontal grid. The majority of data products are time-averaged, but four instantaneous products are also available. Hourly data intervals are used for two-dimensional products, while 3-hourly intervals are used for three-dimensional products. These may be on the model's native 72-layer vertical grid or at 42 pressure surfaces extending to 0.1 hPa. This document describes the gridded output files produced by the GMAO near real-time operational GEOS-5 FP-IT processing in support of the EOS instrument teams. Additional details about variables listed in this file specification can be found in a separate document, the GEOS-5 File Specification Variable Definition Glossary.

3D Var

The Hazard Components of Representative Key Risks. The Physical Climate Perspective

The framework of Representative Key Risks (RKRs) has been adopted by the Intergovernmental Panel on Climate Change Working Group II (WGII) to categorize, assess and communicate a wide range of regional and sectoral key risks from climate change. These are risks expected to become severe due to the potentially detrimental convergence of changing climate conditions with the exposure and vulnerability of human and natural systems. Other papers in this special issue treat each of eight RKRs holistically by assessing their current status and future evolution as a result of this convergence. However, in these papers, such assessment cannot always be organized according to a systematic gradation of climatic changes. Often the big-picture evolution of risk has to be extrapolated from either qualitative effects of “low”, “medium” and “high” warming, or limited/focused analysis of the consequences of particular mitigation choices (e.g., benefits of limiting warming to 1.5 or 2C), together with consideration of the socio-economic context and possible adaptation choices. In this study we offer a representation – as systematic as possible given current literature and assessments – of the future evolution of the hazard components of RKRs. We identify the relevant hazards for each RKR, based upon the WGII authors’ assessment, and we report on their current state and expected future changes in magnitude, intensity and/or frequency, linking these changes to Global Warming Levels (GWLs) to the extent possible. We draw on the assessment of changes in climatic impact-drivers relevant to RKRs described in the 6th Assessment Report by Working Group I supplemented when needed by more recent literature. For some of these quantities - like regional trends in oceanic and atmospheric temperature and precipitation, some heat and precipitation extremes, permafrost thaw and Northern Hemisphere snow cover - a strong and quantitative relationship with increasing GWLs has been identified. For others - like frequency and intensity of tropical cyclones and extra-tropical storms, and fire weather - that link can only be described qualitatively. For some processes - like the behavior of ice sheets, or changes in circulation dynamics - large uncertainties about the effects of different GWLs remain, and for a few others - like ocean pH and air pollution - the composition of the scenario of anthropogenic emissions is most relevant, rather than the warming reached. In almost all cases, however, the basic message remains that every small increment in CO2 concentration in the atmosphere and associated warming will bring changes in climate phenomena that will contribute to increasing risk of impacts on human and natural systems, in the absence of compensating changes in these systems’ exposure and vulnerability, and in the absence of effective adaptation. Our picture of the evolution of RKR-relevant climatic impact-drivers complements and enriches the treatment of RKRs in the other papers in at least two ways: by filling in their often only cursory or limited representation of the physical climate aspects driving impacts, and by providing a fuller representation of their future potential evolution, an important component – if never the only one – of the future evolution of risk severity.

Representative Key Risks

Interannual Variability in Soil Trace Gas (CO2, N2O, NO) Fluxes and Analysis of Controllers

Interannual variability in flux rates of biogenic trace gases must be quantified in order to understand the differences between short-term trends and actual long-term change in biosphere-atmosphere interactions. We simulated interannual patterns (1983-1988) of global trace gas fluxes from soils using the NASA Ames model version of CASA (Carnegie-Ames-Stanford Approach) in a transient simulation mode. This ecosystem model has been recalibrated for simulations driven by satellite vegetation index data from the NOAA Advanced Very High Resolution Radiometer (AVHRR) over the mid-1980s. The predicted interannual pattern of soil heterotropic CO2 emissions indicates that relatively large increases in global carbon flux from soils occurred about three years following the strong El Nino Southern Oscillation (ENSO) event of 1983. Results for the years 1986 and 1987 showed an annual increment of +1 Pg (1015 g) C-CO2 emitted from soils, which tended to dampen the estimated global increase in net ecosystem production with about a two year lag period relative to plant carbon fixation. Zonal discrimination of model results implies that 80-90 percent of the yearly positive increments in soil CO2 emission during 1986-87 were attributable to soil organic matter decomposition in the low-latitudes (between 30 N and 30 S). Soils of the northern middle-latitude zone (between 30 N and 60 N) accounted for the residual of these annual increments. Total annual emissions of nitrogen trace gases (N2O and NO) from soils were estimated to vary from 2-4 percent over the time period modeled, a level of variability which is consistent with predicted interannual fluctuations in global soil CO2 fluxes. Interannual variability of precipitation in tropical and subtropical zones (30 N to 20 S appeared to drive the dynamic inverse relationship between higher annual emissions of NO versus emissions of N2O. Global mean emission rates from natural (heterotrophic) soil sources over the period modeled (1983-1988) were estimated at 57.1 Pg C-CO2yr-1, 9.8Tg (1012 g) N-NO yr-1, and 9.7 Tg N-N2O yr-1. Chemical fertilizer contributions to global soil N gas fluxes were estimated at between 1.3 to 7.3 Tg N-NO yr-1, and 1.2 to 4.0 Tg N-N2O yr-1.

Potter, C.

Integrated Simulation Design Challenges to Support TPS Repair Operations

During the Orbiter Repair Maneuver (ORM) operations planned for Return to Flight (RTF), the Shuttle Remote Manipulator System (SRMS) must grapple the International Space Station (ISS), undock the Orbiter, maneuver it through a long duration trajectory, and orient it to an EVA crewman poised at the end of the Space Station Remote Manipulator System (SSRMS) to facilitate the repair of the Thermal Protection System (TPS). Once repair has been completed and confirmed, then the SRMS proceeds back through the trajectory to dock the Orbiter to the Orbiter Docking System. In order to support analysis of the complex dynamic interactions of the integrated system formed by the Orbiter, ISS, SRMS, and SSMS during the ORM, simulation tools used for previous ‘nominal’ mission support required substantial enhancements. These upgrades were necessary to provide analysts with the capabilities needed to study integrated system performance. Prevalent throughout this ORM operation is a dynamically varying topology. In other words, the ORM starts with the SRMS grappled to the mated Shuttle/ISS stack (closed loop topology), moves to an open loop chain topology consisting of the Shuttle, SRMS, and ISS, and then, at the repair configuration, extends the chain topology to one consisting of the Shuttle, SMS, ISS, and SSRMS/EVA crewman. The resulting long dynamic chain of vehicles and manipulators may exhibit significant motion between the Shuttle worksite and the EVA crewman due to the system flexibility throughout the topology (particularly within the SRMS/SSRMS joints and links). Since the attachment points of both manipulators span the flexible structure of the ISS, simulation analysis may also need to take that into consideration. Moreover, due to the lengthy time duration associated with the maneuver and repair, orbital effects become a factor and require the ISS vehicle control system to maintain active attitude control. Several facets of the ORM operation make the associated analytical efforts different from previous mission support, including: (1) the magnitude of the SRMS handled payload (Le., Orbiter class), (2) the orbital effects induced on the integrated system consisting of the large Shuttle and ISS masses connected by a light flexible SRMS, (3) long duration environmental consequences due to the lengthy operational times associated with the maneuver and repair of the TPS, (4) active attitude control (as opposed to free drift) interacting with the SRMS and SSRMS manipulators (also due to the length of the maneuver and repair), (5) relative dynamics between the EVA crewman and thc worksite influenced by the extended flexible topology. In order to meet these analysis challenges, an ORM simulation architecture was developed leveraging upon numerous pre-existing simulation elements to analyze the various subsystems individually. For example, core manipulator subsystem simulations for both the SRMS and SSRMS were originally combined to provide the dual-arm dynamics topology simulation (in the absence of orbital dynamics and vehicle control). This capability was later merged with the simulation used to analyze SRMS loading with a heavy payload in the orbital environment with an active payload control system (in this case, the ISS Attitude Control System (ACS)), configured for the ORM. The resulting worksite dynamics simulation, based off of the modified ORM simulation, provided the extended topological chain of vehicles and manipulators, while taking into account the orbital effects of both the Shuttle and ISS (as well as its ACS). Verification and validation (V&V) of these integrated simulations became a challenge in itself. A systematic approach needed to be developed such that integration simulation results could be tested against previous constituent simulations upon which these simulations were built. General V&V categories included: (1) core orbital state propagation, (2), stand-alone SRMS, (3) stand-alone SSRMS, (4) stand-alone ISS ACS, (5) integrated Shuttle, SRMS, ISS (with active ACS) in the orbital environment, and (5) dual-arm SRMS/SSRMS dynamics topology. Integrated simulation V&V run suites were created and correlated to verification runs from subsystem simulations, in order to establish the validity of the results. This paper discusses the simulation design challenges encountered while developing simulation capabilities to mirror the ORM operations. The paper also describes the incremental build approach that was utilized, starting with the subsystem simulation elements and integration into increasing more complex simulations until the resulting ORM worksite dynamics simulation had been assembled. Furthermore, the paper presents an overall integrated simulation V&V methodology based upon a subsystem level testing, integrated comparisons, and phased checkout.

Thermal Protection System

Application of OpenFOAM to Plume Impingement in Space Environments

After 30 years of continuous human presence in low-earth orbit, NASA is returning to the moon and eventually will go to Mars. Travelling beyond low earth orbit requires NASA to learn how humans can live in Deep Space environments – beyond the protection of Earth’s magnetosphere and at distances from Earth that prevent a quick return in case of trouble. To this end, NASA is constructing the Lunar Gateway, an ISS-like space station to be put in orbit around the moon to act as a home base for Lunar exploration for NASA astronauts. The Gateway Lunar outpost will be built incrementally, via modules which will arrive at separate times and dock to the existing structure. The incremental addition of Gateway modules, and the docking of visiting vehicles, is achieved via a sequence of firings from the approaching body’s onboard reaction control system (RCS) thrusters to achieve the required approach trajectory. The typical hypergolic chemical RCS thrusters work by firing hot gases to produce adverse thrust and the needed change in velocity to safely finish the docking process. The exhaust gas from the RCS thrusters form plumes that expand into the vacuum of space and can impinge onto the outer surfaces of the Lunar Gateway, causing unwanted forces and moments, heat loads, sediment deposition, and in extreme cases, even surface erosion - all mechanisms that can damage the Lunar Gateway and must be minimized. Both permanent and visiting modules will have this RCS thruster exhaust impingement problem. This research aims to establish existing OpenFOAM solvers as a methodology for improving simulation techniques of rocket exhaust plume impingement in space environments. The flow structure of a plume in a space environment is complex; a plume that originates from a hypergolic chemical RCS thruster and expands into a vacuum will experience several regimes of rarefication. This range includes the continuum flow in the rocket nozzle through the fully rarefied free molecular flow further from the nozzle. The flow physics is different at these two extremes, and as such, the simulation approach for plumes is generally divided into a traditional computational fluid dynamics (CFD) simulation in and near the nozzle which is coupled to a subsequent direct simulation Monte Carlo (DSMC) simulation. At this time, the scope of this research is developing, verifying, and validating a method using existing solvers in the OpenFOAM framework for performing coupled CFD/DSMC calculations to determine the extent of plume impingement loading on generic space structures. This presentation will detail code-to-code comparisons between the hyStrath dsmcFoam+ solver, developed using OpenFOAM and available as open-source, and NASA’s in-house DSMC Analysis Code (DAC). Comparisons to several open-source publication findings using DAC [3,4] are presented, and advantages of using an OpenFOAM based solver are also discussed. The presentation concludes with a discussion of future work, and a plan for coupling the dsmcFoam+ solver with CFD simulations of chemical rocket engines for unified coupled plume simulation.

DSMC

Contracting Quality Early in the Lifecycle Using AS9145 Data Deliverables

Development schedules and a highly dynamic supply chain are a challenge to developers of complex systems produced at low volume. Flaws in designs, parts and materials availability problems, poor manufacturability, and a lack of knowledge about critical items and key process attributes can be realized well before traditional second-party quality assurance activities begin. Supplier audits and product inspections may have little mitigating effect once these foundational problems have been realized. Their impacts can be significant lifecycle disruption, cost overruns, inability to deliver to plan, and even project cancellation. AS9145, Requirements for Advanced Product Quality Planning and Production Part Approval Process, can be used to drive quality engineering practices into early development lifecycles to significantly reduce this late-cycle risk and to reduce the cost of quality overall. Since its initial publication in 2016, it has had very limited adoption by the DoD, no adoption by NASA, and sparse adoption in the aerospace and defense supply chain. A task group within the Aerospace Industries Association's (AIA) Joint Strategic Quality Council (JSQC) identified that both acquirers and suppliers see as AS9145 as a cost-adder and are hesitant to use it as an alternative to late-stage-heavy quality assurance approaches. A lack of prior use creates large capability gaps in request-for-proposal (RFP) teams, proposal teams, suppliers’ quality management systems (QMS), and in experienced personnel executing the early lifecycle approach. To create a more realizable on-ramp for using AS9145 in the space and defense sectors, the AIA JSQC task team created five deliverable requirements descriptions (DRDs) that can be used in a contract to begin to engage both parties in early lifecycle quality engineering and quality assurance activities, that reduce exposure to late-stage cost and schedule collapse due to unidentified risks in design, supply chain, and manufacturability. These DRDs drive the parties to engage in planning and analysis discussions early on to understand what production risks can be known and how they will focus resources based on safety criticality and the key elements of design and construction. The suppliers and acquirers who will produce the data and information required by the DRD will be able to incrementally evolve their QMS and the acquirer will incrementally be able to track and understand the benefits of cost shifting from late to early development phases. A white paper describing this approach and the five recommended DRDs will be published by the AIA in late 2024 or early 2025.

Jeannette Plante

Quantifying Leak Rates and Vacuum Impact on Granular Soil Column Pneumatic Seal

Introduction: All In-Situ Resource Utilization (ISRU) regolith reactor systems that aim to extract gases or vapors must maintain a seal along the flow path of the regolith from the inlet hopper to where extracted gases or vapors are generated. Using a column of regolith at the inlet has been proposed as an alternative to systems using series of mechanical valves to reduce mass and complexity. Preliminary investigations into the use of regolith itself as a barrier to escaping fluids have been conducted, which this research seeks to expand upon. Different vertical regolith simulant plug diameters and heights were tested under vacuum with dry nitrogen to observe leak rates and approximate steady state pressure gradient equilibria. The effect of absolute ambient pressure was also studied by repeating tests at incremental pressures ranging from 10e-5 Torr up to 1 atm to corroborate ambient bench tests from previous work. The impact of regolith dynamics was addressed to emulate continuous-feed regolith reaction processes by measuring leak rates through continuously draining simulant columns. OB-1 lunar highlands regolith simulant was used. Experimental Setup: Regolith simulant columns were constructed using a motorized butterfly valve separating upper acrylic tubes of various diameters from a larger, lower aluminum pipe plenum (See Figure 1). The butterfly valve was cracked open such that gas could freely pass through whereas regolith would bridge and stay suspended as a cohesive plug above in the acrylic portion. These assemblies were loaded into a vacuum chamber, and a compressed dry nitrogen line was routed into the lower plenum. Pressure transducers measured the absolute pressure of the lower plenum, which was compared to absolute chamber pressure to monitor the pressure gradient across the regolith simulant. A precision mass flow controller enabled fine control of plenum pressure. The plenum pressure was quasi-statically increased to set points informed by previous work [1] before shutting off the gas supply to allow the now pressurized gas volume to dissipate through the simulant into the chamber in what was referred to as a ‘leak test.’ In a separate set of tests, the pressure was allowed to quasi-statically build until visual disturbances in the soil were observed via remote camera feeds. Initial leak rate tests were repeated on the disturbed soil to observe changes in leak rates between initial as-poured and post-bubbling soil conditions. This test schedule was applied to simulant columns of three different diameters (approximately 2, 4, and 6 inches) as well as two different heights (around 8 and 16 inches) and was repeated at least three times for each configuration. Plenums were also pressurized while the butterfly valves were opened slightly further to facilitate slow, controlled, constant flowrate draining of the simulant to observe the impact that continuous-feed systems might have on leak rate through the regolith column. Finally, the regolith columns were reset and refilled and static leak rate tests were repeated approximately every decade in ambient pressure from 10e-5 Torr up to ambient conditions to observe the impact of transitioning out of the molecular flow regime on regolith column seal performance. Analysis: Trends with respect to each factor examined were constructed and reported alongside commentary on the most likely responsible physical phenomena. Confidence intervals were applied to the predictive findings of the data by quantifying instrumentation error and accounting for probabilistic impacts of curve fitting the generic exponential function onto pressure decay data. The leak rate is assumed directly proportional to the exponential term coefficient, while the steady-state static pressure gradient is assumed to be the asymptotic limit of each curve fit. Finally, a generic model developed by Ogino et al. [2] and later tuned for ambient testing was re-tuned using vacuum data to estimate leak rates on the lunar surface. Results: Pressure decay rates decreased up to 27-fold under vacuum when compared to ambient pressure, as shown in Figure 2. Simulant column diameter was found to exponentially correlate to both regolith column leak rate and steady-state equilibrium pressure delta. Height was found to be linearly related to leak rate –doubling the regolith column height consistently yielded around a 40% reduction in leak rate. Dynamic regolith draining was found to slightly increase leak rate, although to a lower degree than expected. Conclusion: Because of the transition into molecular gas flow through regolith, vacuum testing is required to properly estimate performance. Regolith columns should be designed as narrow as possible to avoid bridging and minimize leak rate, and as tall as is feasible to both further reduce leak rate and to help prevent bubbling, which has been shown to degrade the quality of the regolith seal. Disrupted soil was found to exhibit higher leak rates than pristine, as-poured soil, indicating the need to consider the impact of recovery from depressurization events. Regolith columns appear equally feasible for both dosed and continuous-feed reactors.

ISRU

Ventilatory responses to dynamic exercise elicited by intramuscular sensors

PURPOSE: Eight subjects, aged 27.0+/-1.6 yr, performed incremental workload cycling to investigate the contribution of skeletal muscle mechano- and metaboreceptors to ventilatory control during dynamic exercise. METHODS: Each subject performed four bouts of exercise: exercise with no intervention (CON); exercise with bilateral thigh cuffs inflated to 90 mm Hg (CUFF); exercise with application of lower-body positive pressure (LBPP) to 45 torr (PP); and exercise with 90 mm Hg thigh cuff inflation and 45 torr LBPP (CUFF+PP). Ventilatory responses and pulmonary gas exchange variables were collected breath-by-breath with concomitant measurement of leg intramuscular pressure. RESULTS: Ventilation (VE) was significantly elevated from CON during PP and CUFF+PP at workloads corresponding to > or = 60% CON peak oxygen uptake (VO2peak) and during CUFF at workloads > or = 80% CON VO2peak, P < 0.05. The VO2 at which ventilatory threshold occurred was significantly reduced from CON (2.17+/-0.28 L x min(-1)) to 1.60+/-0.19 L x min(-1), 1.45+/-0.15 L x min(-1), and 1.15+/-0.11 L x min(-1) during CUFF, PP, and CUFF+PP, respectively. The slope of the linear regression describing the VE/CO2 output relationship was increased from CON by approximately 22% during CUFF, 40% during PP, and 41% during CUFF+PP. CONCLUSIONS: As intramuscular pressure was significantly elevated immediately upon application of LBPP during PP and CUFF+PP without a concomitant increase in VE, it seems unlikely that LBPP-induced increases in VE can be attributed to activation of the mechanoreflex. These findings suggest that LBPP-induced reductions in perfusion pressure and decreases in venous outflow resulting from inflation of bilateral thigh cuffs may generate a metabolite sensitive intramuscular ventilatory stimulus.

Non-NASA Center

An Assessment of Molten Metal Detachment Hazards During Electron Beam Welding in Space

The safety issue has been raised with regards to potential molten metal detachments from the weld pool and cold filler wire during electron beam welding in space. This investigation was undertaken to evaluate if molten metal could detach and come in contact with astronauts and burn through the fabric of the astronauts' Extravehicular Mobility Unit (EMU) during electron beam welding in space. Molten metal detachments from either the weld/cut substrate or weld wire could present harm to a astronaut if the detachment was to burn through the fabric of the EMU. Theoretical models were developed to predict the possibility and size of the molten metal detachment hazards during the electron beam welding exercises at Low Earth Orbit (LEO). The primary molten metal detachment concerns were those cases of molten metal separation from the metal surface due to metal cutting, weld pool splashing, entrainment and release of molten metal due to filler wire snap-out from the weld puddle, and molten metal accumulation and release from the end of the weld wire. Some possible ways of obtaining molten metal drop detachments would include an impulse force, or bump, to the weld sample, cut surface, or filler wire. Theoretical models were developed for these detachment concerns from principles of impact and kinetic energies, surface tension, drop geometry, surface energies, and particle dynamics. The surface tension represents the force opposing the liquid metal drop from detaching whereas the weight of the liquid metal droplet represents a force that is tending to detach the molten metal drop. Theoretical calculations have indicated that only a small amount of energy is required to detach a liquid metal drop; however, much of the energy of an impact is absorbed in the sample or weld plate before it reaches the metal drop on the cut edge or surface. The tendency for detachment is directly proportional to the weld pool radius and metal density and inversely proportional to the surface tension of the liquid metal. For a detachment the initial kinetic energy of the weld pool with respect to the plate has to exceed the energy to form the extra surface required for the detachment of the pool. The difficulty is in transferring the energy from the point of impact through the plate and sample to the cut edge. It is likely that not all of the kinetic energy is available for detaching the pool; some may be sequestered in weld pool oscillations. The coefficient of restitution for the collision will be lower than one if irreversible deformation, for example plastic flow deformation, takes place during the collision. Thus determining the amount of energy from an impact that actually reaches the molten metal droplet is critical. Various molten metal detachment scenarios were tested experimentally in an enclosed vacuum chamber using the Ukrainian Universal Hand Tool, an electron beam welder designed for space welding. The experimental testing was performed in a 4 ft. X 4 ft. vacuum chamber at Marshall Space Flight Center, evacuated to vacuum levels of at least 50 microTorr, and also some welding garment material was utilized to observe the effect of the molten metal detachments on the material. A "carillon" apparatus consisting of four pendulum hammer strikers, each weighing approximately 3.65 lbs, raised to predetermined specific heights was used to apply an impact force to the weld sample/plate during electron beam welding and cutting exercises. The strikers were released by switching on an electric motor to rotate a pin holding wires retaining the strikers at desired heights. The specimens were suspended so as to be free to respond to the blows with a sudden velocity increment. The specimens were mounted on a hinged plate for minimizing effective mass with the option to fasten it down so as to raise its effective mass closer to that anticipated for an actual space welding scenario. Measurements were made of the impact energy and the horizontal fling distances of the detached metal drops. It was not particularly easy to generate the detachments for this experiment. This document presents the details of the theoretical modeling effort and a summary of the experimental effort to measure molten metal drop detachments from terrestrial electron beam welding in the enclosed vacuum chamber. The results of the experimental effort have shown that molten metal detachments can occur from the sample/weld plate only if a sufficiently large impact force is applied to the weld plate. A "weld pool detachment parameter" was determined to indicate whether detachment would occur. Detachment can be either full or partial (dripping), Partial detachment means that the weld pool detached from one side of the liquid-solid boundary so as to leave a hole at the puddle site but remained attached over part of the liquid-solid boundary and dripped down the plate with no fully detached material detected. Full detachment, however, does not necessarily mean that the whole pool fully detached; in some cases only a smaller portion of the pool detached, the remainder dripping down the plate. The weld pool detachment parameter according to theory and according to the empirical data allows a determination of whether full detachments might occur. Theoretical calculations indicated titanium alloy would be the most difficult from which to detach molten metal droplets followed by stainless steel and then by aluminum. The experimental results were for the most part consistent with the theoretical analysis and predictions. The above theory is applicable to other situations as desired for assessing the potential for molten metal detachments.

Fragomeni, James M.