PuMA and Multiscale Modeling
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Camera-based wildlife monitoring is often overwhelmed by non-target triggers and slowed by manual review or cloud-dependent inference, which can prevent timely intervention for high stakes human–wildlife conflicts. Our key contribution is a deployable, fully offline edge vision sensor that achieves near-real-time, highly accurate wildlife event classification by combining detector-based empty-image suppression with a lightweight classifier trained with a staged transfer-learning curriculum. Specifically, Stage 1 uses a pretrained You Only Look Once (YOLO)-family detector for permissive animal localization and empty-trigger suppression, and Stage 2 uses a lightweight EfficientNet-based binary classifier to confirm puma on detector crops and gate downstream actions. Our design is robust to low-quality nighttime monochrome imagery (motion blur, low contrast, illumination artifacts, and partial-body captures) and operates using commercially available components in connectivity-limited settings. In field deployments running since May 2025, end-to-end latency from camera trigger to action command is approximately 4 s. Ablation studies using a dataset of labeled wildlife images (pumas, not pumas) show that the two-stage approach substantially reduces false alarms in identifying pumas relative to a full-image classifier while maintaining high recall. On the held-out test set (N = 1434 events), the proposed two-stage cascade achieves precision 0.983, recall 0.975, F1 0.979, accuracy 0.986, and balanced accuracy 0.983, with only 8 false positives and 12 false negatives. The system can be easily adapted for other species, as demonstrated by rapid retraining of the second stage to classify ringtails. Downstream responses (e.g., notifications and optional audio/light outputs) provide flexible actuation capabilities that can be configured to support intervention.
Cometary grains containing large amounts of carbon and/or organic matter were discovered by in situ measurements of cometary dust composition during VEGA and GIOTTO fly-by missions. In accordance with the classification for the data of PUMA-1 and PUMA-2 mass-spectrometers on board the VEGA spacecraft, particles with a ratio of C to any rock-forming element (Mg, Si, Fe, Ca etc.) greater than 10, were categorized as CHON. There are 464 such particles in PUMA-1 data and 51 in PUMA-2 data. Application of cluster analysis to these grains revealed several distinct compositional classes, namely: (H,C,N,O), (H,C,N), (H,C), (H,C,O), (C,N), (C,O), (C,N,O), and (C). Similar classes were identified among particles analyzed by PIA. Also, about a third of all particles fell into groups (H) and (O) characterized by abundances of these elements beyond chemically reasonable limits.
Blade section normal force and pitching moment were investigated for six rotors operating at transition and high speeds: H-34 in flight and wind tunnel, SA 330 (research Puma), SA 349/2, UH-60A full-scale and BO-105 model (HART-I). The measured data from flight and wind tunnel tests were compared with calculations obtained using the comprehensive analysis CAMRAD II. The calculations were made using two free wake models: rolled-up and multiple-trailer with consolidation models. At transition speed, there is fair to good agreement for the blade section normal force between the test data and analysis for the H-34, research Puma, and SA 349/2 with the rolled-up wake. The calculated airloads differ significantly from the measurements for the UH-60A and BO-105. Better correlation is obtained for the UH-60A and BO-105 by using the multiple-trailer with consolidation wake model. In the high speed condition, the analysis shows generally good agreement with the research Puma flight data in both magnitude and phase. However, poor agreement is obtained for the other rotors examined. The analysis shows that the aerodynamic tip design (chord length and quarter chord location) of the Puma has an important influence on the phase correlation.
The goal of this presentation is to provide a general overview of the multi-scale modeling formulation to determine if there is additional surface recession in Thermal Protection Systems (TPS) materials as a result of mechanical erosion due to high shear conditions during atmospheric entry. This modeling process is performed at different scales by leveraging two computational frameworks developed at NASA: the Porous Microstructure Analysis (PuMA) software, and the Porous material Analysis Toolbox based on OpenFOAM (PATO). PuMA specializes in computing effective macro-scale material properties by performing material response simulations on 3D digital micro-scale representations of porous micro-structures. The modeling of TPS materials at the micro-scale is essential to understand how they behave as part of a heat shield assembly. The first part of the presentation will detail the implementation of PuMA’s cell-centered finite volume elasticity solver, which allows the computation of macro-scale effective mechanical properties of heterogeneous and anisotropic materials such as fibrous and woven TPS composites. These homogenized mechanical properties are used by PATO’s mechanical erosion model to predict the TPS material’s recession at a larger scale. This work will also provide some examples of multi-scale analysis from the fiber level up to the unit cell. The second part of the presentation will focus on the macro-scale approach to determine if erosion at the heat shield’s surface occurs due to mechanical and thermal loads experienced during atmospheric entry. To accomplish this, a solid mechanics module was integrated within PATO enabling it to model the potential mechanical erosion in three steps: first, after obtaining the effective mechanical properties with PuMA, the implemented stress analysis solver computes the stress and the displacement fields for the TPS material using the wall shear stress tensor, computed using a CFD solver, as boundary conditions; then, regions on the surface where the stress meets the failure criteria are identified; finally, the failed material is removed and the mesh is redistributed accordingly. The outcome is a model capable of predicting the total recession in the material due to surface chemistry and mechanical erosion.
The development of a microjoule-class pulsed Doppler lidar and deployment of this compact system on mobile platforms such as aircraft, ships, or trucks have opened a new opportunity to characterize the dynamics of complex mesoscale wind flows. The PickUp-based Mobile Atmospheric Sounder (PUMAS) truck-based lidar system was recently used during the American Wake Experiment (AWAKEN) to assess the general structure of boundary layer (BL) wind and turbulence around wind turbines in central Oklahoma. Wind speed profiles averaged over PUMAS transects influenced by the operating turbines (waked flow) show a 1–2 m s −1 reduction compared to mean undisturbed (free flow) wind speed profiles. Spatial variability in wind speed was observed in time–height cross-sections at different distances from turbines. The wind speeds were about 9–12 m s −1 at 6 km distance compared to 5–7 m s −1 at the transects near the turbines. The PUMAS dataset from AWAKEN demonstrated the capability of the mobile Doppler lidar system to document spatial variability in wind flows at different distances from wind turbines and obtain quantitative estimates of wind speed reduction in the waked flow. The high-frequency, simultaneous measurements of the horizontal and vertical winds provide a new approach for characterizing dynamic processes critical for wind farm wake analyses.
In order to understand the influence of etch pits on the material properties of carbon fiber microstructures, we use the Porous Microstructure Analysis (PuMA) software to create realistic pitted structures from virgin fiber structures and then analyze their degradation as the degree of pitting increases. For that purpose, we developed a pitting module within PuMA that creates pits with user-defined geometry (shape, size) and distribution (defect density) on the surface of arbitrary microstructures. In this work, we consider FiberForm microstructures, which are commonly used as the base material for ablative thermal protection systems (TPS). Using the new pitting module, we mimic the evolution of etch pits on FiberForm as it occurs from oxidation in ablative heating environments, by first importing the virgin FiberForm structure from a micrograph scan and then imposing different pitting conditions by varying pit radii and defect densities. We then leverage PuMA’s capabilities to calculate material properties for each structure. The results clearly show the degradation of thermal conductivity and structural integrity as etch pits grow, and furthermore highlight the complex evolution of the surface topology, which results in significant changes of the flowgeometry and chemical gas-surface interactions. Furthermore, a comparison of these changes with the case of shrinking fibers, which has heretofore been the standard for modelling ablation of carbon structures, yields significant differences. For example, we observe the formation of small chunks of fibers as pits grow in size, which could be released into the surrounding flow environment as spallation particles. This motivates us to continue our efforts of modelling pitting as an integral part of the structural degradation of carbon fiber structures in oxidation environments, which will ultimately improve our understanding of failure mechanisms in those materials.
In order to understand the influence of etch pits on the material properties of carbon fiber microstructures, we use the Porous Microstructure Analysis (PuMA) software to create realistic pitted structures from virgin fiber structures and then analyze their degradation as the degree of pitting increases. For that purpose, we developed a pitting module within PuMA that creates pits with user-defined geometry (shape, size) and distribution (defect density) on the surface of arbitrary microstructures. In this work, we consider FiberForm microstructures, which are commonly used as the base material for ablative thermal protection systems (TPS). Using the new pitting module, we mimic the evolution of etch pits on FiberForm as it occurs from oxidation in ablative heating environments, by first importing the virgin FiberForm structure from a micrograph scan and then imposing different pitting conditions by varying pit radii and defect densities. We then leverage PuMA’s capabilities to calculate material properties for each structure. The results clearly show the degradation of thermal conductivity and structural integrity as etch pits grow, and furthermore highlight the complex evolution of the surface topology, which results in significant changes of the flowgeometry and chemical gas-surface interactions. Furthermore, a comparison of these changes with the case of shrinking fibers, which has heretofore been the standard for modelling ablation of carbon structures, yields significant differences. For example, we observe the formation of small chunks of fibers as pits grow in size, which could be released into the surrounding flow environment as spallation particles. This motivates us to continue our efforts of modelling pitting as an integral part of the structural degradation of carbon fiber structures in oxidation environments, which will ultimately improve our understanding of failure mechanisms in those materials.
The Integrated Robotics System Simulation (ROBSIM) was used to evaluate the performance of the PUMA 560 arm as applied to testing of inertial sensors. Results of this effort were used in the design and development of a feasibility test environment using a PUMA 560 arm. The implemented facility demonstrated the ability to perform conventional static inertial instrument tests (rotation and tumble). The facility included an efficient data acquisitions capability along with a precision test servomechanism function resulting in various data presentations which are included in the paper. Analysis of inertial instrument testing accuracy, repeatability and noise characteristics are provided for the PUMA 560 as well as for other possible commercial arm configurations. Another integral aspect of the effort was an in-depth economic analysis and comparison of robot arm testing versus use of contemporary precision test equipment.
Gross motion control of robotic manipulators typically requires significant on-line computations to compensate for nonlinear dynamics due to gravity, Coriolis, centripetal, and friction nonlinearities. One controller proposed by Luo and Saridis avoids these computations by feeding back joint acceleration and torque. This study implements the controller on a Puma 600 robotic manipulator. Joint acceleration measurement is obtained by measuring linear accelerations of each joint, and deriving a computationally efficient transformation from the linear measurements to the angular accelerations. Torque feedback is obtained by using the previous torque sent to the joints. The implementation has stability problems on the Puma 600 due to the extremely high gains inherent in the feedback structure. Since these high gains excite frequency modes in the Puma 600, the algorithm is modified to decrease the gain inherent in the feedback structure. The resulting compensator is stable and insensitive to high frequency unmodeled dynamics. Moreover, a second compensator is proposed which uses acceleration and torque feedback, but still allows nonlinear terms to be fed forward. Thus, by feeding the increment in the easily calculated gravity terms forward, improved responses are obtained. Both proposed compensators are implemented, and the real time results are compared to those obtained with the computed torque algorithm.
Blade flap bending moments are investigated for six rotors operating at transition and high speeds: H-34 in flight and wind tunnel, SA 330 (research Puma), SA 349/2, UH-60A full-scale, and BO-105 model (HART-I). The measured data from flight and wind tunnel tests are compared with calculations obtained using the comprehensive analysis CAMRAD II. The calculations mere made using two free wake models: rolled-up and multiple-trailer with consolidation models. At transition speed, there is fair to good agreement for the flap bending moment between the test data and analysis for the H-34, research Puma, and SA 349/2 with the rolled-up wake. The calculated flap bending moments differ significantly from measurements for the UH-60A and BO-105. Better correlation is obtained for the UH-60A by using the multiple-trailer with consolidation wake model. Although the multiple-trailer with consolidation wake model shows good correlation on the normal force for the BO-105, the same analysis shows poor correlation on the flap bending moment. In the high speed condition, the analysis shows generally good agreement with the research Puma flight data in both magnitude and phase. However, poor agreement is obtained for the other rotors examined. Although the analysis significantly underpredicts the vibratory normal force on the advancing side for the H-34, the vibratory bending moment correlation is fair to good on both magnitude and phase.
Blade flap and chord bending and torsion moments are investigated for six rotors operating at transition and high speed: H-34 in flight and wind tunnel, SA 330 (research Puma), SA 349/2, UH-60A full-scale, and BO- 105 model (HART-I). The measured data from flight and wind tunnel tests are compared with calculations obtained using the comprehensive analysis CAMRAD II. The calculations were made using two free wake models: rolled-up and multiple-trailer with consolidation models. At transition speed, there is fair to good agreement for the flap and chord bending moments between the test data and analysis for the H-34, research Puma, and SA 349/2. Torsion moment correlation, in general, is fair to good for all the rotors investigated. Better flap bending and torsion moment correlation is obtained for the UH-60A and BO-105 rotors by using the multiple-trailer with consolidation wake model. In the high speed condition, the analysis shows generally better correlation in magnitude than in phase for the flap bending and torsion moments. However, a significant underprediction of chord bending moment is observed for the research Puma and UH-60A. The poor chord bending moment correlation appears to be caused by the airloads model, not the structural dynamics.
The Portable Unit for Metabolic Analysis (PUMA) is an instrument that measures several quantities indicative of human metabolic function. Specifically, this instrument makes time-resolved measurements of temperature, pressure, flow, and the partial pressures of oxygen and carbon dioxide in breath during both inhalation and exhalation. Portable instruments for measuring these quantities have been commercially available, but the response times of those instruments are too long to enable temporal resolution of phenomena on the time scales of human respiration cycles. In contrast, the response time of the PUMA is significantly shorter than characteristic times of human respiration phenomena, making it possible to analyze varying metabolic parameters, not only on sequential breath cycles but also at successive phases of inhalation and exhalation within the same breath cycle. In operation, the PUMA is positioned to sample breath near the subject s mouth. Commercial off-the-shelf sensors are used for three of the measurements: a miniature pressure transducer for pressure, a thermistor for temperature, and an ultrasonic sensor for flow. Sensors developed at Glenn Research Center are used for measuring the partial pressures of oxygen and carbon dioxide: The carbon dioxide sensor exploits the relatively strong absorption of infrared light by carbon dioxide. Light from an infrared source passes through the stream of inhaled or exhaled gas and is focused on an infrared- sensitive photodetector. The oxygen sensor exploits the effect of oxygen in quenching the fluorescence of ruthenium-doped organic molecules in a dye on the tip of an optical fiber. A blue laser diode is used to excite the fluorescence, and the optical fiber carries the fluorescent light to a photodiode, the temporal variation of the output of which bears a known relationship with the rate of quenching of fluorescence and, hence, with the partial pressure of oxygen. The outputs of the sensors are digitized, preprocessed by a small onboard computer, and then sent wirelessly to a desktop computer, where the collected data are analyzed and displayed. In addition to the raw data on temperature, pressure, flow, and mole fractions of oxygen and carbon dioxide, the display can include volumetric oxygen consumption, volumetric carbon dioxide production, respiratory equivalent ratio, and volumetric flow rate of exhaled gas.
Updates on NASA‘s efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. The PMM effort is part of the Entry Systems Modeling (ESM) project under NASA’s Game Changing Development (GCD) program. To reduce the need for extensive testing and accelerate the design cycle process, ESM is developing simulation and modeling tools that enable the characterization of the properties of thermal protection materials and their response to extremely hot plasma. The Porous Microstructure Analysis (PuMA) software has been developed to compute effective material properties and perform material response simulations on digitized microstructures of porous media. PuMA is able to import three-dimensional digital images obtained from X-ray microtomography or to generate artificial microstructures that mimic real materials. PuMA also provides a module for interactive 3D visualizations. Version 3, which was recently released as open-source, includes modules to compute simple morphological properties such as porosity, volume fractions, pore diameter, and specific surface area. Additional capabilities include the determination of effective thermal and electrical conductivity (both radiative and solid conduction - including the ability to simulate local anisotropy for the latter); effective diffusivity and tortuosity from the continuum to the rarefied regime; techniques to determine the local material orientation, as well as mechanical properties (elasticity coefficients), and permeability. Computed properties are then used to inform a macro-scale material response model, such as those implemented in the Porous material Analysis Toolbox based on OpenFOAM (PATO) software developed within ESM. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale, assuming Local Thermal Equilibrium. Recent efforts include the development of a mechanical erosion model as well as a unified model allowing an intrinsic coupling between fluid and material. Comparison to flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI] and Mars 2020 MEDLI2) is critical in order to validate these computational tools. Examples of ablative material response using the code will be presented, including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrated the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains through the use of massively parallel computations.
Tracking morphological changes of materials during heating is crucial to understand its response in fire protection, biofuel production, thermal protection systems (TPS) for hypersonic flight. As materials are heated, they undergo physical and chemical changes due to water loss, stretching or shrinking, pyrolysis and chemical reactions in the ambient environment. The effects of these changes can have a profound impact on the material’s performance, indicated by changes in on the porosity and volume. While materials such as wood shrink as they pyrolyze and lose mass, others swell due to their inherent characteristics when exposed to heat [1]. This study focuses on experiments conducted at the Advanced Light Source (ALS) beamline 8.3.2, where in situ micro-computed tomography (µ-CT) is performed on materials as they are being pyrolyzed. Through in situ µ-CT, the change in total volume and porosity can be obtained in real-time, allowing for better understanding of the underlying thermophysical and chemical processes as a function of temperature. This study also focuses on the implementation of the Porous Microstructure Analysis software (PuMA) [2] to obtain thermal conductivity, permeability, and other properties of the material from the 3D tomographies. The information gained from these tomographies will supplement microscale model development of material morphological change and will aid macroscale modeling for high-temperature applications. For this study, Room Temperature Vulcanizing silicone (RTV) [3-5] is heated from room temperature to 1000°C using an infrared lamp heating system, and tomographies are continuously collected as the sample is heated. The tomographies are then segmented to obtain solid and void phases, from which estimates of pore size, porosity and total volume are extracted as a function of temperature. PuMA is deployed on the segmented tomographies to obtain thermal conductivity, permeability, and other properties as a function of temperature. Preliminary results show that RTV first intumesces (swells) as pyrolysis begins, due to build-up of pyrolysis gases in closed pores, and then shrinks significantly as more open pores are formed and the pyrolysis gases outgas. Pore network visualization of the tomographies using OpenPNM [6] showed the increase in pore connectivity with increase in temperature. Future work will focus on using PuMA to obtain macroscopic properties of RTV as a function of temperature.