Search NASA⌕ Search

SEARCH · Search NASA

Results for “visualizing data variation”

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.

393 records · Page 22

Integrative Modeling and Analysis of Fungal Central Carbon Metabolism

Over a thousand fungal genomes have been sequenced, yet manually curated genome-scale metabolic models (GEMs) are available for only a limited number of species. Moreover, these models have often been developed independently, leading to inconsistencies in namespaces, compartment definitions, and pathway representations that hinder comparative analysis, the systematic reuse of prior curation efforts, and the integration of consolidated metabolic knowledge. Here, we present the Consolidated Fungal Core Metabolism Model (CFCMM), constructed by integrating thirteen published fungal models spanning Ascomycota, Mucoromycota, and both Crabtree-positive and Crabtree-negative yeasts. We harmonized metabolites and reactions into a non-redundant shared ModelSEED ontological space, standardized compartmentalization, and refined gene–protein–reaction (GPR) rules. Using pathway-level visualization and systematic gap detection, we further improved the integrated network through literature-guided curation to correct stoichiometry, stereospecificity, and pathway architecture. Orthologous protein family reconstruction and functional annotation workflows were used to validate and inform GPR associations, with particular emphasis on ambiguous enzyme superfamilies and membrane-associated components. Using the resulting CFCMM, we built high-quality central carbon core models for each fungus and performed flux balance analysis to quantify ATP-yield variation under aerobic and anaerobic conditions, explicitly evaluating scenarios driven by differences in electron transport chain (ETC) composition. Simulations reproduced the expected fermentative yield of approximately 2 mmol ATP per mmol glucose under anaerobic conditions and separated the thirteen fungi into two bioenergetic groups under aerobic respiration based on Complex I status, with predicted yields of approximately 30 versus 22 mmol ATP per mmol glucose. Forcing flux through the alternative oxidase bypass further reduced ATP yields to approximately 12 and 4 mmol ATP per mmol glucose in Complex I-containing and Complex I-lacking fungi, respectively. Collectively, this work provides a manually curated, ModelSEED-consistent, and extensible fungal core metabolic template, deployed in DOE KBase as a resource for automated reconstruction of central carbon core models from any sequenced fungal genome. In addition, the CFCMM provides modular components for developing GEMs with more accurate energy predictions and enables robust comparative analyses of fungal bioenergetics and core metabolic diversity

59 BASIC BIOLOGICAL SCIENCES↗

VAIM-CFF: a variational autoencoder inverse mapper solution to Compton form factor extraction from deeply virtual exclusive reactions

We develop a new methodology for extracting Compton form factors (CFFs) from deeply virtual exclusive reactions such as the unpolarized DVCS cross section using a specialized inverse problem solver, a variational autoencoder inverse mapper (VAIM). The VAIM-CFF framework not only allows us access to a fitted solution set possibly containing multiple solutions in the extraction of all 8 CFFs from a single cross section measurement, but also accesses the lost information contained in the forward mapping from CFFs to cross section. We investigate various assumptions and their effects on the predicted CFFs such as cross section organization, number of extracted CFFs, use of uncertainty quantification technique, and inclusion of prior physics information. We then use dimensionality reduction techniques such as principal component analysis to visualize the missing physics information tracked in the latent space of the VAIM framework. Through re-framing the extraction of CFFs as an inverse problem, we gain access to fundamental properties of the problem not comprehensible in standard fitting methodologies: exploring the limits of the information encoded in deeply virtual exclusive experiments.

Accelerator Physics↗

SAR backscatter from coniferous forest gaps

A study is in progress comparing Airborne Synthetic Aperture Radar (AIRSAR) backscatter from coniferous forest plots containing gaps to backscatter from adjacent gap-free plots. Issues discussed are how do gaps in the range of 400 to 1600 sq m (approximately 4-14 pixels at intermediate incidence angles) affect forest backscatter statistics and what incidence angles, wavelengths, and polarizations are most sensitive to forest gaps. In order to visualize the slant-range imaging of forest and gaps, a simple conceptual model is used. This strictly qualitative model has led us to hypothesize that forest radar returns at short wavelengths (eg., C-band) and large incidence angles (e.g., 50 deg) should be most affected by the presence of gaps, whereas returns at long wavelengths and small angles should be least affected. Preliminary analysis of 1989 AIRSAR data from forest near Mt. Shasta supports the hypothesis. Current forest backscatter models such as MIMICS and Santa Barbara Discontinuous Canopy Backscatter Model have in several cases correctly predicted backscatter from forest stands based on inputs of measured or estimated forest parameters. These models do not, however, predict within-stand SAR scene texture, or 'intrinsic scene variability' as Ulaby et al. has referred to it. For instance, the Santa Barbara model, which may be the most spatially coupled of the existing models, is not truly spatial. Tree locations within a simulated pixel are distributed according to a Poisson process, as they are in many natural forests, but tree size is unrelated to location, which is not the case in nature. Furthermore, since pixels of a simulated stand are generated independently in the Santa Barbara model, spatial processes larger than one pixel are not modeled. Using a different approach, Oliver modeled scene texture based on an hypothetical forest geometry. His simulated scenes do not agree well with SAR data, perhaps due to the simple geometric model used. Insofar as texture is the expression of biological forest processes, such as succession and disease, and physical ones, such as fire and wind-throw, it contains useful information about the forest, and has value in image interpretation and classification. Forest gaps are undoubtedly important contributors to scene variance. By studying the localized effects of gaps on forest backscatter, guided by our qualitative model, we hope to understand more clearly the manner in which spatial heterogeneities in forests produce variations in backscatter, which collectively give rise to scene texture.

Day, John L.↗

Model, data, and code for paper "Modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x"

The data package includes data, model, and code that support the analyses and conclusions in the paper titled “modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x”. The primary goal of this paper is to demonstrate that key streamflow properties such as water depth, flow velocity, and dynamic pressure in a natural river at 30-kilometer scale over 5 years can be reliably and efficiently modeled using the computational framework presented in this paper. To support the paper, various data types from remote sensing, field observations, and computational models are used. Specific details are described as follows. Firstly, the river bathymetry data was obtained from a Light Detection and Ranging (LiDAR) survey. This data is then converted to a triangulated surface format, STL, for mesh generation in OpenFOAM. The STL data can be found in Model_Setups/BaseCase_2013To2015/constant/triSurface. The OpenFOAM mesh generated using this STL file can be found in constant/polyMesh. Other model setups, boundary and initial conditions can be found in /system and /0.org under folder BaseCase_2013To2015. A similar data structure can also be found in BaseCase_2018To2019 for the simulations during 2018 and 2019. Secondly, the OpenFOAM simulations need the upstream discharge and water depth information at the upstream boundary to drive the model. These data are generated from a one-dimensional hydraulic model and the data can be found under the folder Model_Setups /1D model Mass1 data. The mass1_65.csv and mass1_191.csv files include the results of the 1D model at the model inlet and outlet, respectively. The Matlab source code Mass1ToOFBC20182019.m is used to convert these data into OpenFOAM boundary condition setups.With the above OpenFOAM model, it can generate data for water surface elevation, flow velocity, and dynamic pressure. In this paper, the water surface elevation was measured at 7 locations during different periods between 2011 and 2019. The exact survey locations (see Fig1_SurveyLocations.txt) can be found in folder Fig_1. The variation of water stage over time at the 7 locations can be found in folder /Observation_WSE. The data type include .txt, .csv, .xlsx, and .mat. The .mat data can be loaded by Matlab.We also measured the flow velocities at 12 cross-sections along the river. At each cross-section, we recorded the x, y locations, depth, three velocity components u,v,w. These data are saved to a Matlab format which can be found under folder /Observation_Velocity and /Fig_1. The relative locations of velocity survey locations to the river bathymetry can be found in Figure 1c.The water stage data at the 7 locations from OpenFOAM, 1D, and 2D hydraulic models are also provided to evaluate the long-term performance of 3D models vs 1D/2D models. The water stage data for the 7 locations from OpenFOAM have been saved to .mat format and can be found in /OpenFOAM_WSE. The water stage data from the 1D model are saved in .csv format and can be found in /Mass1_WSE. The water stage from the 2D model is saved as .mat format and can be found in / Mass2_WSEIn addition, the OpenFOAM model outputs the information of hydrostatic and hydrodynamic pressure. They are saved as .mat format under folder /Fig_11/2013_1. As the files are too large, we only uploaded the data for January 2013. The area of different ratio of dynamic pressure to static pressure for all simulation range, i.e., 2013-2015, are saved to .mat format. They can be found in /Fig_11/PA. Further, the data of wall clock time versus the solution time of the OpenFOAM modeling are also saved to .mat format under folder /Fig_13/LogsMat. In summary, the data package contains seven data types, including .txt, .csv, .xlsx, .dat, .stl, .m, and .mat. The former 4 types can be directly open using a text editor or Microsoft Office. The .mat format needs to be read by Matlab. The Matlab source code .m files need to be run with Matlab. The OpenFOAM setups can be visualized in ParaView. The .stl file can be opened in ParaView or Blender. The data in subfolders Fig_1 to Fig_10 and Fig_12 are copied from the aforementioned data folders to generate specific figures for the paper. A readME.txt file is included in each subfolder to further describe how the data in each folder are generated and used to support the paper.Please use the data package's DOI to cite the data package. Please contact yunxiang.chen@pnnl.gov if you need more data related to the paper.

54 ENVIRONMENTAL SCIENCES↗

A finite element/level set model of polyurethane foam expansion and polymerization

Polyurethane foams are used widely for encapsulation and structural purposes because they are inexpensive, straightforward to process, amenable to a wide range of density variations (1 lb/ft3 - 50 lb/ft3), and able to fill complex molds quickly and effectively. Computational model of the filling and curing process are needed to reduce defects such as voids, out-of-specification density, density gradients, foam decomposition from high temperatures due to exotherms, and incomplete filling. This paper details the development of a computational fluid dynamics model of a moderate density PMDI structural foam, PMDI-10. PMDI is an isocyanate-based polyurethane foam, which is chemically blown with water. The polyol reacts with isocyanate to produces the polymer. PMDI- 10 is catalyzed giving it a short pot life: it foams and polymerizes to a solid within 5 minutes during normal processing. To achieve a higher density, the foam is over-packed to twice or more of its free rise density of 10 lb/ft3. The goal for modeling is to represent the expansion, filling of molds, and the polymerization of the foam. This will be used to reduce defects, optimize the mold design, troubleshoot the processed, and predict the final foam properties. A homogenized continuum model foaming and curing was developed based on reaction kinetics, documented in a recent paper; it uses a simplified mathematical formalism that decouples these two reactions. The chemo-rheology of PMDI is measured experimentally and fit to a generalized- Newtonian viscosity model that is dependent on the extent of cure, gas fraction, and temperature. The conservation equations, including the equations of motion, an energy balance, and three rate equations are solved via a stabilized finite element method. The equations are combined with a level set method to determine the location of the foam-gas interface as it evolves to fill the mold. Understanding the thermal history and loads on the foam due to exothermicity and oven curing is very important to the results, since the kinetics, viscosity, and other material properties are all sensitive to temperature. Results from the model are compared to experimental flow visualization data and post-test X-ray computed tomography (CT) data for the density. Several geometries are investigated including two configurations of a mock structural part and a bar geometry to specifically test the density model. We have found that the model predicts both average density and filling profiles well. However, it under predicts density gradients, especially in the gravity direction. Further model improvements are also discussed for future work.

36 MATERIALS SCIENCE↗

Optic Nerve Sheath Mechanics in VIIP Syndrome

Visual Impairment Intracranial Pressure (VIIP) syndrome is a major concern in current space medicine research. While the exact pathology of VIIP is not yet known, it is hypothesized that the microgravity-induced cephalad fluid shift increases intracranial pressure (ICP) and drives remodeling of the optic nerve sheath. To investigate this possibility, we are culturing optic nerve sheath dura mater samples under different pressures and investigating changes in tissue composition. To interpret results from this work, it is essential to first understand the biomechanical response of the optic nerve sheath dura mater to loading. Here, we investigated the effects of mechanical loading on the porcine optic nerve sheath.Porcine optic nerves (number: 6) were obtained immediately after death from a local abattoir. The optic nerve sheath (dura mater) was isolated from the optic nerve proper, leaving a hollow cylinder of connective tissue that was used for biomechanical characterization. We developed a custom mechanical testing system that allowed for unconfined lengthening, twisting, and circumferential distension of the dura mater during inflation and under fixed axial loading. To determine the effects of variations in ICP, the sample was inflated (0-60 millimeters Hg) and circumferential distension was simultaneously recorded. These tests were performed under variable axial loads (0.6 grams - 5.6 grams at increments of 1 gram) by attaching different weights to one end of the dura mater. Results and Conclusions: The samples demonstrated nonlinear behavior, similar to other soft connective tissue (Figure 1). Large increases in diameter were observed at lower transmural pressures (approximately 0 to 5 millimeters Hg), whereas only small diameter changes were observed at higher pressures. Particularly interesting was the existence of a cross-over point at a pressure of approximately 11 millimeters Hg. At this pressure, the same diameter is obtained for all axial loads applied to the tissue; i.e., as the axial load is varied, the diameter of the dura mater remains constant. This cross-over in the pressure-diameter curves occurred in all optic nerve sheaths that were tested, and may correspond with in vivo ICP levels for pigs. These data suggest that diameter of the dura mater of the optic nerve remains nearly constant in vivo despite being stretched axially. This may be a homeostatic mechanism aimed at maintaining target stresses/strains on the cells in the dura mater, and deviations from these stresses may play an important role in optic nerve sheath remodeling. Future studies will involve subjecting the dura mater to varying pressures and axial tensions for extended periods of time, while monitoring changes in the biomechanical properties. The data can then be used to study the effects of changes in ICP on the remodeling of the dura mater.

biomechanics↗

GEO Population Estimates using Optical Survey Data

Optical survey data taken using the NASA Michigan Orbital Debris Survey Telescope (MODEST) gives us an opportunity to statistically sample faint object population in the Geosynchronous (GEO) and near-GEO environment. This paper will summarize the MODEST survey work that has been conducted by NASA since 2002, and will outline the techniques employed to arrive at the current population estimates in the GEO environment for dim objects difficult to detect and track using current systems in the Space Surveillance Network (SSN). Some types of orbits have a higher detection rate based on what parts of the GEO belt is being observed, a straightforward statistical technique is used to debias these observations to arrive at an estimate of the total population potentially visible to the telescope. The size and magnitude distributions of these fainter debris objects are markedly different from the catalogued population. GEO debris consists of at least two different populations, one which follows the standard breakup power law and one which has anomalously high Area-to-Mass Ratios (1 to approx. 30 square meters per kilogram; a sheet of paper = approx. 13 square meters per kilogram). The Inter-Agency Space Debris Coordination Committee (IADC) is investigating objects in GEO orbits with anomalously high Area-to-Mass Ratios (AMRs). The ESA Space Debris Telescope discovered this population and has and its properties of inclinations (0 to 30 degrees), changing eccentricities (0 and 0.6), and mean motions (approx. 1 rev), will be presented. The accepted interpretation of this orbital behavior is that solar radiation pressure drives the perturbations causing time varying inclinations and eccentricities. The orbital parameters are unstable for this population and thus difficult to predict. Their dim visual magnitudes and photometric variability make observations a challenge. The IADC has enlisted a series of observatories (participating institutions: University of Michigan/CTIO, Astronomical Institute University of Bern, Boeing LTS / AMOS, Keldysh Institute of Applied Mathematics) at different longitudes. Complete observational coverage over periods of days to months will provide a better understanding of the properties, such as solar radiation pressure effects on orbital elements, size, shape, attitude, color variations, and spectral characteristics. Results from recent observational programs will be summarized, and includes a description of the orbit elements prediction processes, a summary of the metric tracking performance, and some photometric characteristics of this class of debris.

Barker, Edwin S.↗

AGR-1 Compact 4-1-1 Post-Irradiation Examination Results

Destructive post-irradiation examination was performed on AGR-1 fuel Compact 4-1-1, which was irradiated to a final compact-average burnup of 19.4% FIMA (fissions per initial metal atom) and a time-average, volume-average temperature of 1072°C. The analysis of this compact focused on characterizing the extent of fission product release from the particles and examining particles to determine the condition of the kernels and coating layers. The work included deconsolidation of the compact and leach-burn-leach analysis, visual inspection and gamma counting of individual particles, metallurgical preparation of selected particles, and examination of particle cross-sections with optical microscopy, electron microscopy, and elemental analysis. Deconsolidation-leach-burn-leach (DLBL) analysis revealed no particles with failed TRISO or failed SiC layers (as indicated by very low uranium inventory in all of the leach solutions). The total fractions of the predicted compact inventories of fission products Ce-144, Cs-134, Cs-137, and Sr-90 that were present in the compact outside of the SiC layers were <2×10 -6 , based on DLBL data. The Ag-110m fraction in the compact outside the SiC layers was 3.3×10 -2 , indicating appreciable release of silver through the intact coatings and subsequent retention in the OPyC layers or matrix. The Eu-154 fraction was 2.4×10 -4 , which is equivalent to the inventory in one average particle, and indicates a small but measurable level of release from the intact coatings. Gamma counting of 61 individual particles indicated no particles with anomalously low fission product retention. The average ratio of measured inventory to calculated inventory was close to a value of 1.0 for several fission product isotopes (Ce-144, Cs-134, and Cs-137), indicating good retention and reasonably good agreement with the predicted inventories. Measured-to-calculated (M/C) activity ratios for fission products Eu-154, Eu-155, Ru-106, Sb-125, and Zr-95 were significantly less than 1.0. However, as no significant release of these fission products from compacts was noted during previous analysis of the AGR-1 capsule components, the low M/C ratios are most likely an indication of a bias in the inventories predicted by physics simulations of the AGR-1 experiment. The distribution of Ag-110m M/C ratios was centered on a value of 1.02 and was fairly broad (standard deviation of 0.18, with values as high as 1.42 and as low as 0.68). Based on all data gathered to date, it is believed that silver retention in the particles was on average relatively high, but that the broad distribution in values among the particles represents significant variation in the inventory of Ag-110m generated in the particles. Ceramographic analysis of particle cross-sections revealed many of the characteristic microstructures often observed in irradiated AGR-1 particles from other fuel compacts. Palladium-rich fission product clusters were observed in the IPyC and SiC layers near the IPyC-SiC interface of three Compact 4-1-1 particle cross-sections. In spite of the presence of fission product clusters in the SiC layer, no significant corrosion or degradation of the layer was observed in any of the particles examined.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Reliability Assessment of Remote Visual Examination

Remote visual testing (RVT) is a commonly used nondestructive examination method for inservice inspection (ISI) of reactor internals to detect cracking and gross component failures. Despite widespread use, the detection reliability of RVT and the factors that impact overall RVT performance have been unresolved issues. This report describes the results from an assessment sponsored by the U.S. Nuclear Regulatory Commission and conducted by the Pacific Northwest National Laboratory, in cooperation with the Electric Power Research Institute, for evaluating the reliability of RVT methods currently being used for [reactor] in-vessel visual inspection. The goal was to assess the performance of commercially applied RVT examination procedures implemented by qualified personnel, as well as to identify and qualitatively assess enhancements to examination procedures for detecting flaws in test specimens. The assessment was performed over three phases of research, with the results of the first phase used to identify key controllable parameters impacting RVT performance and to design the next two phases (round-robin testing). Participants in the round-robin tests included teams of inspectors from commercial nuclear power ISI vendors. Participants were asked to determine if a specimen contained a crack and its approximate location, orientation, and length. The detection and location information were used to compute estimates of the probability of detection and false call probability for various scenarios such as the placement and orientation of the flaw on the specimen, team-to-team variations in detection performance, and effects of secondary review on RVT detection. Results showed that crack opening displacement (COD) is the dominant factor in the reliability of crack detection using commercially applied RVT procedures, with crack length being weakly correlated with detection probability. The implication is that RVT detection is likely heavily dependent on the contrast produced by the crack opening, with crack detection becoming less reliable as the COD decreases. Further, the results indicated RVT will be challenged when cracks are located in the vicinity of surface features such as scratches or weld ripples, or close to the edge of welds where shadowing and/or the presence of weld undercuts may complicate the ability to detect the crack. The assessment also reinforced earlier findings on the importance of lighting in flaw detection and appeared to align with other studies that find improved reliability when using multiple inspectors or analysts. Finally, the results pointed to the importance of practice with specimens that mimic conditions likely to be found in the field. Based on the findings and limitations of the assessment, a number of recommendations are made regarding best practices for improving the reliability of RVT in a field setting and for addressing remaining informational gaps. This report describes the design and noted limitations of the three phases of research, the analysis methodology for each phase, and the results of the research. The results of this assessment provide a benchmark set of data on the reliability of RVT for detecting cracking, assuming the implementation of field-like inspection procedures. The likely impact of several uncontrolled factors on RVT detection performance are discussed, and recommendations regarding the use of these results to assess field performance are provided. Finally, recent advances in RVT technology are briefly discussed and point to the potential need for continued research to evaluate the capability and effectiveness of the technique as improvements are implemented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DEVELOPMENT OF AN INFLIGHT COUNTERMEASURE TO MITIGATE POSTFLIGHT GAIT DYSFUNCTION

Following spaceflight crewmembers experience gait and postural instabilities due to inflight adaptive alterations in sensorimotor function. These changes can pose a risk to crew safety if nominal or emergency vehicle egress is required immediately following long-duration spaceflight. At present, no operational countermeasure is available to mitigate postflight locomotor disturbances. Therefore, the goal of this study is to develop an inflight training regimen that facilitates the recovery of locomotor function after long-duration spaceflight. The countermeasure we are developing is based on the concept of variable practice. During this type of training the subject gains experience producing the appropriate adaptive motor behavior under a variety of sensory conditions and response constraints. This countermeasure is built around current ISS treadmill exercise activities. Crewmembers will conduct their nominal inflight treadmill exercise while being exposed to variations in visual flow patterns. These variations will challenge the postural and locomotor systems repeatedly, thereby promoting adaptive reorganization in locomotor behavior. As a result of this training a subject learns to solve a class of motor problems, rather than a specific motor solution to one problem, Le., the subject learns response generalizability or the ability to "learn to learn" under a variety of environmental constraints. We anticipate that this training will accelerate recovery of postural and locomotor function during readaptation to gravitational environments following spaceflight facilitating neural adaptation to unit (Earth) and partial (Mars) gravity after long-duration spaceflight. The study calls for one group of subjects to perform the inflight treadmill training regimen while a control group of subjects performs only the nominal exercise procedures. Locomotor function in both groups is assessed before and after spaceflight using two tests of gait function: The Integrated Treadmill Locomotion Test (ITLT) and the Functional Mobility Test (FMT). The ITLT characterizes alterations in the integrated function of multiple sensorimotor subsystems responsible for the control of locomotion. This test calls for subjects to walk on a motorized treadmill while we assess changes in dynamic postural stability, head-trunk coordination, short-latency head stabilization responses, dynamic visual acuity, lower limb coordination strategies and gait cycle timing. To make these assessments we measure the following parameters while subjects walk on the treadmill: 1) full body 3-dimensional kinematics using a motion capture system (Motion Analysis Corp., Santa Rosa, CA); 2) the shock-wave transmitted from heel-strike to the head using triaxial accelerometers placed on the tibia and head (Entran, Fairfield, NJ); 3) vertical forces using an instrumented treadmill (Kistler Instrument Corp., Amherst, NY); 4) Dynamic visual acuity using Landolt Cs presented on a laptop computer located 4m from the eyes and 5) Gait cycle timing using foot-switches (Motion Lab Systems, Inc., Baton Rouge, LA) attached to the plantar surface of each shoe at the heel and toe. The FMT evaluates a subject's ability to perform challenging locomotor maneuvers similar to those encountered during an egress from a space vehicle. Subjects step over and duck under obstacles along with negotiating a series of pylons set up on a base of 10 cm thick medium density foam. The dependent measures for the FMT are time to complete the course and the number of obstacles touched. To date, we have collected pre and postflight locomotion data from Expeditions 5-9 who will serve as part of the control group for this study. Preliminary results comparing the recovery rates in gait control sub-systems obtained from the ITLT and FMT performance showed two recovery patterns: 1) a concordant recovery trend between gait control parameters and FMT performance indicating a restitution pattern of recovery and 2) gait controecovery that lagged recovery in FMT performance suggesting that improvement in locomotor function was attained through a pattern of substitution. These data suggest that recovery of postflight locomotor function may occur through adaptive mechanisms that lead to either restitution or substitution of function. Understanding the modes of postflight readaptation has implications for countermeasure development and testing and in astronaut postflight rehabilitation.

Bloomberg, J. J.↗

Development of an Inflight Countermeasure to Mitigate Postflight Gait Dysfunction

Following spaceflight crewmembers experience gait and postural instabilities due to inflight adaptive alterations in sensorimotor function. These changes can pose a risk to crew safety if nominal or emergency vehicle egress is required immediately following long-duration spaceflight. At present, no operational countermeasure is available to mitigate postflight locomotor disturbances. Therefore, the goal of this study is to develop an inflight training regimen that facilitates the recovery of locomotor function after long-duration spaceflight. The countermeasure we are developing is based on the concept of variable practice. During this type of training the subject gains experience producing the appropriate adaptive motor behavior under a variety of sensory conditions and response constraints. This countermeasure is built around current ISS treadmill exercise activities. Crewmembers will conduct their nominal inflight treadmill exercise while being exposed to variations in visual flow patterns. These variations will challenge the postural and locomotor systems repeatedly, thereby promoting adaptive reorganization in locomotor behavior. As a result of this training a subject learns to solve a class of motor problems, rather than a specific motor solution to one problem, Le., the subject learns response generalizability or the ability to "learn to learn" under a variety of environmental constraints. We anticipate that this training will accelerate recovery of postural and locomotor function during readaptation to gravitational environments following spaceflight facilitating neural adaptation to unit (Earth) and partial (Mars) gravity after long-duration spaceflight. The study calls for one group of subjects to perform the inflight treadmill training regimen while a control group of subjects performs only the nominal exercise procedures. Locomotor function in both groups is assessed before and after spaceflight using two tests of gait function: The Integrated Treadmill Locomotion Test (ITLT) and the Functional Mobility Test (FMT). The ITLT characterizes alterations in the integrated function of multiple sensorimotor subsystems responsible for the control of locomotion. This test calls for subjects to walk on a motorized treadmill while we assess changes in dynamic postural stability, head-trunk coordination, short-latency head stabilization responses, dynamic visual acuity, lower limb coordination strategies and gait cycle timing. To make these assessments we measure the following parameters while subjects walk on the treadmill: 1) full body 3-dimensional kinematics using a motion capture system (Motion Analysis Corp., Santa Rosa, CA); 2) the shock-wave transmitted from heel-strike to the head using triaxial accelerometers placed on the tibia and head (Entran, Fairfield, NJ); 3) vertical forces using an instumented treadmill (Kistler Instrument Corp., Amherst, NY); 4) Dynamic visual acuity using Landolt Cs presented on a laptop computer located 4m from the eyes and 5) Gait cycle timing using foot-switches (Motion Lab Systems, Inc., Baton Rouge, LA) attached to the plantar surface of each shoe at the heel and toe. The FMT evaluates s. subject's ability to perform challenging locomotor maneuvers similar to those encountered during an egress from a space vehicle. Subjects step over and duck under obstacles along with negotiating a series of pylons set up on a base of 10 cm thick medium density foam. The dependent measures for the FMT are time to complete the course and the number of obstacles touched. To date, we have collected pre and postflight locomotion data from Expeditions 5-9 who will serve as part of the control group for this study. Preliminary results comparing the recovery rates in gait control sub-systems obtained from the ITLT and FMT performance showed two recovery patterns: 1) a concordant recovery trend between gait control parameters and FMT performance indicating a restitution pattern of recovery and 2) gait controecovery that lagged recovery in FMT performance suggesting that improvement in locomotor function was attained through a pattern of substitution. These data suggest that recovery of postflight locomotor function may occur through adaptive mechanisms that lead to either restitution or substitution of function. Understanding the modes of postflight readaptation has implications for countermeasure development and testing and in astronaut postflight rehabilitation.

Bloomberg, J. J.↗

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

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

Lohrer, J. D.↗

EXCHANGE Campaign 1: A Community-Driven Baseline Characterization of Soils, Sediments, and Water Across Coastal Gradients

The EXploration of Coastal Hydrobiogeochemistry Across a Network of Gradients and Experiments (EXCHANGE) program is a consortium of scientists working together to improve our understanding of how the two-way exchange of water between estuaries or large lake lacustuaries and the terrestrial landscape influence the state and function of ecosystems across the coastal interface. EXCHANGE Campaign 1 (EC1) focuses on the spatial variation in biogeochemical structure and function at the coastal terrestrial-aquatic interface (TAI). In the Fall of 2021, the EXCHANGE Consortium gathered samples from 52 TAIs. Samples collected from EC1 were analyzed for bulk geochemical parameters, bulk physicochemical parameters, organic matter characteristics, and redox-sensitive elements.Please download ec1_README.pdf for a complete list of available data in each .zip folder, package version history, and detailed information about the project. This README will serve as the central place for EC1 Data Package updates. Experimental setup and v1 methods are documented in Myers-Pigg and Pennington et al., 2023 (https://doi.org/10.1038/s41597-023-02548-7).EC1 Data Package Structure:ec1_README.pdfec1_methods.pdfec1_metadata_v3.zip...ec1_dd.csv...ec1_flmd.csv...ec1_sample_catalog.csv...ec1_metadata_kitlevel.csv...ec1_metadata_collectionlevel.csv...ec1_data_collectionlevel.csv...ec1_igsn_metadata.csvec1_soil_v3.zipec1_sediment_v3.zipec1_water_v3.zipec1_processingscripts_v3.zipThis data package is on v3 and was originally published May 2023 (v1). Subsequent updates will be published here with new version numbers. Please see the Change History section in ec1_README.pdf for detailed changes.---Acknowledging EXCHANGE: General Support and Data Product UseWe ask that users of EXCHANGE data add the following acknowledgement when publishing data in scholarly articles and data repositories:"This research is based on work supported by COMPASS-FME, a multi-institutional project supported by the U.S. Department of Energy, Office of Science, Biological and Environmental Research as part of the Environmental System Science Program."

54 ENVIRONMENTAL SCIENCES↗

Particle-image Velocimetry Sensitivity through Automatic Differentiation

A commonly encountered difficulty in modern experimental fluid dynamics experiments is the need for full-flow visualization and measurement; traditional techniques such as hot-wire probes and laser Doppler velocimetry cannot provide such information without difficult scanning procedures. Particle-image velocimetry (PIV) is a measurement technique which can overcome this spatial limitation to provide real-time three-dimensional fluid velocity vectors for an entire plane within the flow, as opposed to a single point that is provided by other measurement techniques. To accomplish this, a PIV system seeds the fluid flow with tracer particles; these must have nearly the same density as that of the fluid and a high drag coefficient to ensure that they follow the fluid flow without deviation. A laser is sent through optics to expand its beam into a sheet which illuminates the plane of interest in the fluid. The particles reflect more light than the fluid, and their image is recorded using a high speed camera. Due to limitations in bandwidth, video cannot be taken fast enough to be useful, and instead pairs of images are captured with a known separation in time. The illuminated particles will move between the two images taken by the camera, and their relative motion can be used to compute the fluid velocity at that point. When a single camera is used, only two-dimensional velocity vectors can be computed, but with the addition of another camera this can be improved to full three-dimensional vectors. The flow's time evolution can be measured by taking series of image pairs and processing each pair to get the velocity field at the time between each image pair. Only the two-dimensional version will be considered in this work. To compute the velocity vectors of the flow, the images are broken up into interrogation regions, usually in a uniform grid. In the most basic PIV algorithm, corresponding regions in the image pairs are cross-correlated to compute the most likely displacement from one image to the next. This is used since the cross correlation will achieve a maximum value when the two interrogation regions line up the best. The location of the peak in the correlation map and scaling information for the optical setup is used to compute the particle displacement from one image to the next. This process is similar to holding two semi-transparent images and moving one until they match as well as possible. Sub-pixel displacements between images can be obtained by fitting a Gaussian function to the correlation peak. Knowing the displacement vector for each set of interrogation regions in the image pair, it is a simple matter to divide by the time gap between the two images to compute the fluid velocity at the point in time and space. Since there are many interrogation regions in each image pair, the velocity field for a sheet is sampled according to the arrangement of the interrogation regions. The sensitivity of the final computed vector to the input images is difficult to compute, and has led most considerations of uncertainty in the PIV field to apply Monte Carlo methods and correlations to quantify the systematic and random errors which result from application of various PIV algorithms. In this work, automatic differentiation (AD) has been used to trace the dependence of velocity vectors computed using the correlation-peak PIV algorithm back to the input values given to an artificial image generator. These initial results form the foundation for continued work in the application of AD to algorithms in PIV. They demonstrate that AD methods can be successfully applied to PIV algorithms without undue complexity. The derivatives of each individual vector with regard to the coefficients of the velocity field used to create the artificial images can be computed and the statistics of them considered collectively, or correlated between each other to quantify the effectiveness of a PIV algorithm. With the initial implementation complete, adding extra algorithms is relatively simple and considerable data can be generated for analysis. It is hoped that application of these methods can lead to improved PIV algorithms which demonstrate better sensitivity to small variations in particle displacement. While these methods can be used as an analysis tool for refinement of new algorithms, computation of actual resultant vector sensitivities to input images' individual pixels is the next logical step of this work. Such information will be analyzed with experimental data to produce a better understanding of how output vectors' uncertainty relates to the actual image pairs, and lead to automatic bad vector removal. (authors)

42 ENGINEERING↗

A Study of Nucleate Boiling with Forced Convection in Microgravity

The ultimate objective of basic studies of flow boiling in microgravity is to improve the understanding of the processes involved, as manifested by the ability to predict its behavior. This is not yet the case for boiling heat transfer even in earth gravity, despite the considerable research activity over the past 30 years. The elements that constitute the nucleate boiling process - nucleation, growth, motion, and collapse of the vapor bubbles (if the bulk liquid is subcooled) - are common to both pool and flow boiling. It is well known that the imposition of bulk liquid motion affects the vapor bubble behavior relative to pool boiling, but does not appear to significantly influence the heat transfer. Indeed, it has been recommended in the past that empirical correlations or experimental data of pool boiling be used for design purposes with forced convection nucleate boiling. It is anticipated that such will most certainly not be possible for boiling in microgravity, based on observations made with pool boiling in microgravity. In earth gravity buoyancy will act to remove the vapor bubbles from the vicinity of the heater surface regardless of how much the imposed bulk velocity is reduced, depending, of course, on the geometry of the system. Vapor bubbles have been observed to dramatically increase in size in pool boiling in microgravity, and the heat flux at which dryout took place was reduced considerably below what is generally termed the critical heat flux (CHF) in earth gravity, depending on the bulk liquid subcooling. However, at heat flux levels below dryout, the nucleate pool boiling process was enhanced considerably over that in earth gravity, in spite of the large vapor bubbles formed in microgravity and perhaps as a consequence. These large vapor bubbles tended to remain in the vicinity of the heater surface, and the enhanced heat transfer appeared to be associated with the presence of what variously has been referred to as a liquid microlayer between the bubble and the heater surface. The enhancement of the boiling process with low velocities in earth gravity for those orientations producing the formation of a liquid macrolayer described above, accompanied by "sliding" vapor bubbles, has been demonstrated. The enhancement was presented as a function of orientation, subcooling, and heated length, while a criterion for the heat transfer for mixed natural/forced convection nucleate boiling was given previously. A major unknown in the prediction and application of flow boiling heat transfer in microgravity is the upper limit of the heat flux for the onset of dryout (or critical heat flux - CHF), for given conditions of fluid-heater surfaces, including geometry, system pressure and bulk liquid subcooling. It is clearly understood that the behavior in microgravity will be no different than on earth with sufficiently high flow velocities, and would require no space experimentation. However, the boundary at which this takes place is still an unknown. Previous results of CHF measurements were presented for low velocity flow boiling at various orientations in earth gravity as a function of flow velocity and bulk liquid subcooling, along with preliminary measurements of bubble residence times on a flat heater surface. This showed promise as a parameter to be used in modeling the CHF, both in earth gravity and in microgravity. The objective of the work here is to draw attention to and show results of current modeling efforts for the CHF, with low velocities in earth gravity at different orientations and subcoolings. Many geometrical possibilities for a heater surface exist in flowing boiling, with boiling on the inner and outer surfaces of tubes perhaps being the most common. If the vapor bubble residence time on and departure size from the heater surface bear a relationship to the CHF, as results to be given indicate, it is important that visualization of and access to vapor bubble growth be conveniently available for research purposes. In addition, it is desirable to reduce the number of variables as much as possible in a fundamental study. These considerations dictated the use of a flat heater surface, which is rectangular in shape, 1.91 cm by 3.81 cm (0.75 x 1.5 inches), consisting either of a 400 Angstrom thick semi-transparent gold film sputtered on a quartz substrate which serves simultaneously as a heater and a resistance thermometer, or a copper substrate of the same size. The heater substrate is a disc which can be rotated so that the heated length in the flow direction can be changed from 1.91 to 3.81 cm (0.75 to 1.5 inches). The fluid is R-113, and the velocities can be varied between 0.5 cm/s and 60 cm/s. For a sufficiently low velocity the CHF can be modeled reasonably well at various orientations by the correlation for pool boiling corrected for the influence of bulk liquid subcooling, multiplied by the square root of q, the angle relative to horizontal. This arises from equating buoyancy and drag forces in the inverted positions where the vapor bubbles are held against the heater surface as they slide. A distortion of the measurements relative to pool boiling occurs as the flow velocity increases. In modeling this effect at different levels of subcooling it appeared appropriate to estimate the volumetric rate of vapor generation, using measurements of bubble frequency (or residence time), void fraction and average bubble boundary layer thickness. These were determined with the use of a platinum hot wire probe 0.025 mm in diameter by 1.3 mm long, applying a constant current to distinguish between contact with liquid or vapor. Two-dimensional spatial variations are obtained with a special mechanism to resolve displacements in increments of 0.025 mm. From such measurements it was determined that the fraction of the surface heat transfer resulting in evaporation varies inversely with the subcooling correction factor for the CHF. The measured inverse bubble residence time is normalized relative to that predicted for an infinite horizontal flat plate at the CHF, and is correlated well with the CHF normalized relative to that for pool boiling, for various orientation angles and subcooling levels. This correspondence is then combined with a normalizing factor for the energy flux leaving the heater surface at the CHF and the computed bubble radius at departure, determined from the balance between the outward velocity of the interface due to evaporation and the buoyance induced velocity of the center of mass of the bubble. The product of the CHF and the corresponding residence time was determined to be a constant for all orientations at a given bulk flow velocity and liquid subcooling, and must be determined empirically for each velocity and subcooling at present. It then becomes possible to predict the CHF for the different orientations, velocities, and subcoolings. These are compared with normalized measurements of the CHF for velocities ranging from 4 cm/s to 55 cm/s, subcoolings from 2.8 to 22.2 K, over orientations angles of 360 degrees.

Merte, Herman, Jr.↗