Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sorting”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,423 records · Page 79

Generation-based Evolutionary Tool for the Optimization of Constellations (GenETOC)

With the rapid growth in the capabilities of smaller satellites, satellite architectures that replace a single, extremely capable spacecraft with multiple, cheaper ones are gaining in popularity. Unfortunately, the orbit design process for constellations can be significantly more involved, especiallywhen the relative placement of the individual spacecraft within the constellation is not constrained by mission and/or science objectives. Optimizing a satellite constellation in the presence of multiple, competing objectives is a highly complex problem to which many traditional mathematical optimization methods cannot be applied and few tools exist to help mission designers search for promising candidate mission designs. The Generation-based Evolutionary Tool for the Optimization of Constellations (GenETOC) has been created to search for near-optimal constellation design options. GenETOC combines a modified version of the Non-dominated Sorting Genetic Algorithm II (NSGA II) with STK Components libraries (a 3rdparty .NET package created by Analytical Graphics Inc.) to create a framework that enables a mission designer to generate a simulation that models the design problem and obtain a family of potential, near-optimal solutions that can be investigated more in detail. GenETOC was developed in C# using the .NET framework with Windows Presentation Foundation (WPF) serving as the framework from which to create the graphical user interface (GUI). GenETOC user inputs can be categorized into three major data components: definition of the problem (areas of interest, satellite decision parameters, and sensor configurations), definition of performance objectives, and specification of the genetic algorithm (GA) parameters. In the problem definition component, the user is prompted to define the areas of interest against which the performance metrics will be computed, define the sensor parameters and attach them to specific spacecraft, select which satellite orbital parameters will be added to the decision space of the GA, and specify the range of desired values for each optimization parameter. For performance objectives, the user is presented with a list of available coverage and revisit performance based calculation options from which two metrics are chosen to serve as the objective functions that the GA will use to evaluate solutions during the optimization process. Finally, the definition of the GA parameters provides user control over the number of generations (number of optimization iterations), the population size (number of candidate constellations created in each generation), and the adaptive mutation and crossover threshold values (control parameters for how frequently each process occurs during the optimization). GenETOC has been extensively tested to verify the individual components of the optimization process. The GA has been tested against a suite of GA test problems to confirm convergence to the known two and three-dimensional Pareto fronts. The coverage and revisit performance metrics obtained in GenETOC are compared with STK desktop scenarios, confirming the constellations are being appropriately modeled within GenETOC simulations. A walkthrough of a simple, example problem is provided to illustrate the workings of GenETOC and to demonstrate the output available to the mission designer.

mission design↗

An Updated Secondary Lunar Meteoroid Ejecta Model for Engineering Design

Introduction: The surface of the Moon is constantly being bombarded by a flux of meteoroids of various sizes. Impacts due to these meteoroids produce secondary ejecta material at much lower speeds but with a total mass larger than the original impactor. Details about the secondary ejecta are important for planning missions on the lunar surface. In this work, an updated ejecta model is presented called the Meteoroid Model of Secondary Ejecta (MeMoSeE), to re-place the Apollo-era ejecta model, NASA SP-8013 [1], in the SLS-SPEC 159 Design Specification for Natural Environments (DSNE) [2]. The model pro-duces secondary ejecta flux environments for a user-specified location on the lunar surface, and sorts the incoming secondary flux by angular direction and speed. Methods: MeMoSeE is separated into three parts: the inputs, the conversion step, and the integration of fluxes. Inputs to the model include the primary meteoroid fluxes and the primary near-Earth object (NEO) fluxes. Meteoroid fluxes, both asteroidal and cometary, are calculated using the Meteoroid Engineering Model (MEM3) [3] for different locations on the Moon. For each surface location, an ephemeris is generated using the JPL HORIZONS System [4] that feeds into MEM3. The NEO fluxes are approximated by the high-density population of MEM3 (i.e., only the directionality), where the speed distribution of the NEO fluxes is renormalized to match observations [5]. The regolith properties are used as defined in the DSNE [2]. The conversion step utilizes scaling laws given by Housen & Holsapple 2011 [6] to convert the primary impactor flux to the total mass of secondary ejecta. The ejecta distribution, at the point-of-impact (POI), is separated into a zenith angle and azimuthal angle distribution. The zenith angle distribution follows a beta distribution where the peak depends on the impact altitude angle and the impact azimuth [7]. We employ an ejecta azimuth distribution that is based on Rival & Mandeville 1999 [8] which focuses ejecta in the downstream direction for more oblique impacts. Finally, during the integration step, we sum secondary ejecta number fluxes at a particular region-of-interest (ROI) that originated from many POI locations over the entire surface of the Moon. We keep track of both altitude and azimuth angle bins as well as a range of speed bins, following the igloo gridding as done in MEM3 [9]. The ejecta particle size distribution and density is assumed to be the same as the lunar regolith [10, 2]. Results: The primary fluxes are computed for one Metonic cycle (19 years) for various locations over the lunar surface with a fixed orientation. Both the angular and speed distributions of the primary fluxes are dependent on the latitude and longitude. The over-all primary fluxes show a roughly 13% increase from the eastern limb to the western limb. In general, the speed distribution of the secondary fluxes span from a user-defined minimum speed to the escape speed of the Moon (2.38 km/s), roughly following a power-law relation [6]. Different parts of the speed distribution come from different primary impact locations on the Moon. The secondary ejecta is dominated by the slowest speeds, where these particles originate nearby the ROI. For speeds around 71% of the escape speed, the secondary ejecta originates from locations near the antipodal point. On the other hand, for secondary ejecta speeds that exceed roughly 90% the escape speed, the ejecta particles come from all over the lunar surface to the ROI. Comparing the secondary ejecta fluxes from MeMoSeE with NASA SP-8013 [1], there is a reduction by about 2-3 orders of magnitude for secondary ejecta particles greater than 1 µg. These estimates agree with recent findings from Bjorkman & Christiansen 2019 [11]. The secondary ejecta fluxes are also compared with the primary fluxes, where the ejecta fluxes are roughly an order of magnitude greater than the primary fluxes. References: [1] Cour-Palais, B. G., (1969) NASA SP-8013. [2] NASA SLS-SPEC-159 Rev. H (2020). [3] Moorhead, A. V., et al. (2019) JS&R, 1-17. [4] Giorgini, J. D., (2015) IAUGA, 29, 2256293. [5] Moorhead, A. V., (2020) Memo OSMA/MEO/Lunar-001. [6] Housen, K. R., and Holsapple, K. A., (2011) Icarus, 211(1), 856-875. [7] Gault, D. E., and Wedekind, J. A., (1978) L&PSCP, 9, 3843-3875. [8] Rival, M., and Mandeville, J., (1999) Space Debris, 1(1), 45-57. [9] Moorhead, A. V., (2019) MEM3 User Guide. [10] Carrier, W. D., (1973) The Moon, 6(3-4), 250-263. [11] Bjorkman, M. D. & Christiansen, E. L., (2019) ODC

impacts↗

A Systematic Approach to Hardware Qualification

A systematic approach for the development of a hardware qualification approach is described. This approach stems from the need to break the "Catch-22" of not being able to fly new technology because it hasn't flown. A physics of failure approach is used to identify failure modes and the impact and likelihood of these failures on the mission requirements is plotted in a Requirements Matrix. These same failure modes are plotted against the effectiveness of the available Preventions, Analyses, Control and Tests (PACTs) at screening for, or eliminating, these failure modes in a Test Effectiveness Matrix. Matrix multiplication results in a ranked set of PACTs which can be sorted according to cost and redundancy with other PACTs. This, and other information which results from the process, will allow project managers to make more informed decisions regarding the cost and risk tradeoffs inherent in any qualification program.

hardware↗

Radar Observations of Asteroid 1620 Geographos

Goldstone radar observations of Geographos from Aug. 28 through Sep. 2, 1994, yield over 400 delay-Doppler images whose linear spatial resolutions range from ~75 m to ~151 m, and 138 pairs of dual-polarization (OC, SC) spectra with one-dimensional resolution of 103 m. Each data type provides thorough rotational coverage. The images contain an intrinsic north/south ambiguity, but the equatorial view allows accurate determination of the shape of the radar-facing part of the asteroid's pole-on silhouette at any rotation phase. Sums of co-registered images that cover nearly a full rotation have defined the extremely elongated shape of that silhouette (Ostro et al. 1995, Nature 375, 474-477). Here we present individual images and co-registered sums over ~30û of rotation phase that show the silhouette's structural characteristics in finer detail and also reveal numerous contrast features "inside" the silhouette. Those features include several candidate craters as well as indications of other sorts of large-scale topographic relief, including a prominent central indentation. Protuberances at the asteroid's ends may be related to the pattern of ejecta removal and deposition caused by the asteroid's gravity field. The asteroid's surface is homogeneous and displays only modest roughness at centimeter-to-meter scales. Our estimates of radar cross section and the currently available constraints on the asteroid's dimensions are consistent with a near-surface bulk density between 2 and 3 g cm-3. The delay-Doppler trajectory of Geographos' center of mass has been determined to about 200 m on Aug. 28 and to about 100 m on Aug. 31, an improvement of two orders of magnitude over pre-observation ephemerides.

Asteroid↗

Applications of the Dynamic N-Dimensional K-Vector

The n-dimensional k-vector (NDKV) is an appealing alternative to binary tress for resolving complex queries in large relational databases. The method has excelled in several applications involving static databases. The present paper extends the theory supporting the NDKV to handle dynamic databases, where the data is updated frequently. This includes deleting records, adding new entries, or editing existing elements. The merit of this new version of the NDKV, the dynamic n-dimensional k-vector (DNDKV), is that it is no longer necessary to recompute the entire k-vector (the main structure that indexes the data) every time a record changes. The algorithm updates the four constituents of the standard NDKV on the fly: the database, sorted database, index, and k-vector tables. As a result, the DNDKV becomes comparable in terms of capabilities and flexibility to stateof-the-art storage engines relying on structured query languages (SQL). The performance of the DNDKV is assessed by running typical read/write operations on a database that contains millions of pre-computed missions to celestial bodies. This database requires frequent updates whenever an orbit solution is refined or new bodies are discovered. The DNDKV is faster than rebuilding the k-vector tables completely, provided that the number of elements being added or removed is not excessively large. Direct runtime comparisons with MySQL suggest that the DNDKV is several times faster for reading but might be slower for writing and updating the database. One limit of the technique is the elements being added must be within the range of the current k-vector tables. If this is not the case, the technique cannot be used and the k-vector tables must be rebuilt from scratch.

Mortari, Daniele↗

Modeling the Radiative Effect on Microphysics in Cirrus Clouds against Satellite Observations

The radiative effect on microphysics 23 (REM) plays an important role in the dew/frost formation near the surface. How REM impacts cirrus clouds is investigated in this paper, using bin microphysical model simulations and coincident data of the CloudSat and Global Precipitation Measurement (GPM) satellites. REM affects ice crystal spectrum with two types: radiative cooling and warming. Radiative cooling, as predicted by the bin model simulations, favors the formation of horizontally oriented ice crystals (HOICs), but radiative warming does not. Hence, a test of REM can be transformed to a test of HOICs, because HOICs can be measured by the microwave polarization observations of the GPM Microwave Imager (GMI) at 166 GHz. To analyze the GMI data for their HOIC distribution, clouds are sorted into four groups with different optical depth and altitude, based on the radiative cooling/warming ratio (or eta) computed with satellite-retrieved ice water content. Their HOIC distributions (e.g., the mid-level thick clouds have more HOICs than the high-level ones) agree well with those predicted by the bin model simulations. The general agreement between the GMI observations and bin model simulations suggests that REM is common in cirrus clouds and impacts cirrus clouds significantly.

microphysics↗

An Updated Secondary Lunar Meteoroid Ejecta Model for Engineering Design

The surface of the Moon is constantly being bombarded by a flux of meteoroids of various sizes. Impacts due to these meteoroids produce secondary ejecta material at much lower speeds but with a total mass larger than the original impactor. Details about the secondary ejecta are important for planning missions on the lunar surface. In this work, an updated ejecta model is presented called the Meteoroid Model of Secondary Ejecta (MeMoSeE), to replace the Apollo-era ejecta model, NASA SP-8013 [1], in the SLS-SPEC 159 Design Specification for Natural Environments (DSNE) [2]. The model produces secondary ejecta flux environments for a user-specified location on the lunar surface, and sorts the incoming secondary flux by angular direction and speed.

Secondary Ejecta↗

Radio Frequency Plasma Synthesis of Boron Nitride Nanotubes (BNNTs) for Structural Applications: Part III

The emergence of commercial quantities of BNNTs facilitates the exploration of structural material applications. BNNT-reinforced composites must exhibit extraordinary mechanical properties in order to justify the cost premium. Prevailing theory indicates that suitable BNNTs must possess an aspect ratio (AR= length/diameter) approaching 10,000 to provide effective discontinuous reinforcement under tensile loading. A survey of commercial products shows that although the quantities are sufficient, lots of uniform morphology are unavailable. Therefore, the opportunity exists to supplement procedures either during synthesis operations or post-synthesis processing in order to deliver homogeneous batches of high AR BNNTs. An assessment of relevant technologies reveals established methods compatible with sorting of BNNTs by physical dimensions, i.e., size and/or shape. Electrostatic fractionation emerges as the leading candidate, with the potential to isolate high-purity, low-defect, high AR product. The adaptation of existing equipment to the handling of fibrous particulate on the nanoscale will be paramount to success. Although currently focused on BNNTs, the concepts outlined may be applied to nanotubes synthesized from any organic or inorganic compounds. This survey sets the stage for delivering reproducible lots of BNNTs that satisfy the dimensional requirements for use in acreage structural materials.

induction plasma, BNNT, aspect ratio, fractionatio↗

Particle Interaction Physics Model Formulation for Plume-Surface Interaction Erosion and Cratering

As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The Predictive Simulation Capability (PSC) Element is focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the predictive simulation capability, the GCD Program also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper will present the status of implementing and maturing particle-particle interaction constituent physics models essential in simulating the landing surface granular material flow under PSI effects. This gas-particle multi-phase interaction modeling of plume impingement flow on the extra-terrestrial soil material is performed with the Gas-Granular Flow Solver (GGFS) addressed in a companion paper. The response of regolith particle flow induced by lander PSI requires accurate representation of the regolith granular material fluidic behavior and gas-granular interactions. The lunar regolith, as the extreme example, is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking particle shapes for the very jagged particles. The combination of particle shape and size distribution has been identified as major drivers in the complex particle flow response and resulting crater shape characteristics of extraterrestrial granular material. Constituent models for spherical particles can be formulated directly from particle kinetics theory. Complex particle shapes can be modeled by gluing together elemental spherical shapes into composite particles, requiring a Discrete Element Model (DEM) particle kinetics modeling approach to extract data and formulate constituent models. Mixture constituent models for poly-disperse mixtures (i.e, containing distribution of particle sizes) have recently been developed. The required non-spherical particle mixture granular material response closure models are then obtained through small-scale unit physics DEM simulations for the range of particle shapes, mixtures and packing densities. The granular material response closure models are then implemented in the Eulerian granular flow formulation. This DEM-based constituent model extraction process and formulation of poly-disperse particle mixtures has been successfully developed by small business and academic partners in the development of the Gas-Granular Flow Solver (GGFS) simulation program simulation framework. The currently implemented capabilities have reached the capability level of modeling bi-disperse, non-spherical particle mixtures is being continuously extended towards computational modeling of full range irregular particle mixtures. Under the GCD project, this technology is being further developed, transferred to NASA analysts, and matured towards application readiness. The predictive simulation capability team under the GCD project has acquired the modeling tools and processes of the DEM based constituent model formulation from the GGFS development team and is developing the capability to replicate the existing process. This is the first important step towards the ability of the NASA team to independently perform such model development in a production setting. Further efforts are underway to migrate the DEM based model simulation process performed with the academic based tools to more capable Open Source, highly parallelized simulation tools for efficient operation on NASA HPC assets. Evaluation of the currently implemented (such as mono-disperse and bi-disperse spherical and irregular shape particle constituent model applications) and continuously evolving full-range particle physics models in the GGFS tool is performed by the NASA team to advance application readiness of the simulations. Application testing for complex PSI erosions and cratering scenarios such as the Apollo LM is performed for axi-symmetric and full 3D simulations to aid the tool developers in achieving practical application readiness for NASA projects. Important validation and application testing will further be performed against experimental data generated under the GCD PSI project experimental component.

Peter A Liever↗

Modeling the Radiative Effect on Microphysics in Cirrus Clouds Against Satellite Observations

The radiative effect on microphysics (REM) plays an important role in the dew/frost formation near the surface. How REM impacts cirrus clouds is investigated in this study, using bin microphysical model simulations and coincident data of the CloudSat and Global Precipitation Measurement (GPM) satellites. REM affects ice crystal spectrum with two types: radiative cooling and warming. Radiative cooling, as predicted by the bin-model simulations, favors the formation of horizontally oriented ice crystals (HOICs), but radiative warming does not. Hence, a test of REM can be transformed to a test of HOICs, because HOICs can be measured by the microwave polarization observations of the GPM Microwave Imager (GMI) at 166 GHz. To analyze the GMI data for their HOIC distribution, clouds are sorted into four groups with different optical depth and altitude, based on the radiative cooling/warming ratio (or eta) computed with satellite-retrieved ice water content. Their HOIC distributions (e.g., the midlevel thick clouds have more HOICs than the high-level ones) agree well with those predicted by the bin-model simulations. The general agreement between the GMI observations and bin-model simulations suggests that REM is common in cirrus clouds and impacts cirrus clouds significantly.

clouds↗

Enabling Single Cell Research in Hollow Hydrogel Microparticles

Studying biology at the single cell level enables novel discoveries through ultra-high throughput identification and isolation of rare cells and interesting subpopulations. Encapsulation of single cells in hollow hydrogel microparticles, or PicoShells, enables easy visualization, manipulation, and assaying of single cells. The PicoShell’s porous hydrogel shell allows diffusion of nutrients and wastes into and out of the hollow liquid core, allowing uninhibited cell growth. The porous shell also allows diffusion of assay reagents, enabling solution changes, and multi-step assays. Finally, single cells and their progeny can be tracked through time, for example to measure single cell growth rates. However, PicoShells are fabricated using microfluidic techniques, which require specialized expertise, equipment, and facilities. This limits PicoShell’s widespread utility for single cell biology. For Saccharomyces cerevisiae and other desiccation-tolerant microbes, desiccating PicoShells containing cells of interest offers a solution, enabling those without microfluidic capabilities to obtain dried, shelf-stable PicoShells containing their cells of interest from a specialized microfluidics facility. Researchers could then rehydrate the PicoShells and perform their single cell assays, bypassing the microfluidic steps. Desiccated PicoShells could also be flown in space, where desiccation is already a common practice for microbes (i.e. BioSentinel and BioNutrients), bringing single cell analysis capabilities to space. In this work, we develop a method to desiccate PicoShells containing S. cerevisiae while maintaining characteristics similar to fresh, non-desiccated PicoShells. We show good microparticle morphology and hydrogel pore size as well as high yeast viability in line with previous studies. We also demonstrate sorting of a fast-growing yeast population using size-based filtration, as S. cerevisiae in PicoShells are able to physically stretch the microparticle as they grow, increasing the PicoShell diameter. These efforts contribute to the feasibility of leveraging microparticle PicoShells for single cell space biology research on Earth and in space.

Simon Ng↗

Constraints on the Origin of Mercury’s Large Core from Core-Mantle Differentiation Models

Mercury’s core is notoriously large when compared to other planets of our solar system. The origin of this large core is still uncertain. Available data on the surface composition and internal structure of Mercury from the past MESSENGER mission and future data collected by BepiColombo will continue to provide clues to Mercury’s formation. Here, we present results combining experimental data on elemental distribution between core and mantle with spacecraft data, to estimate bulk Mercury composition. We applied this strategy to major elements (Fe, Si, Mg, Al and O), as well as minor elements (Cr and Ti), and compared derived compositions to chondritic data and the chemical compositions of the other terrestrial planets. Our results show that Mercury has a chemical composition significantly different from all known materials of the solar system. In addition, numerous scenarios were proposed to explain Mercury’s structure, including “chaotic models” such as giant impacts and “orderly models” such as aerodynamic sorting. Here, we tested whether Mercury’s composition can be explained by mantle stripping by impacts. We will show how such a scenario reconciles several features of Mercury’s geochemistry with chondritic data. We will also discuss the outlook of additional constraints from supplementary data potentially collected by BepiColombo.

Mercury↗

Habitable Working and Living Spaces in the Deep Space Science Vessel

The Deep Space Science Vessel (DSSV) is a conceptual design of a very large spacecraft intended as a mobile microgravity science platform. It represents an effort to conceptualize the systems and internal architectures needed to support a forty-eight-person crew for multi-year missions throughout the inner solar system and potentially beyond. The general arrangement of the DSSV is a modular spacecraft. Considering only the pressurized, habitable elements, the DSSV includes two large pressurized modules (Hab and Lab) docked together along with eight smaller node modules completing the habitable configuration. The Hab and Lab modules are docked side by side, with the longitudinal axes parallel to each other and a single docking port connecting them. A node module is docked to each dome on the Hab and Lab modules, such that the two small node modules on each end of the large Hab and Lab domes are docked both to each other and to the Hab and Lab. This creates a “racetrack” translation path on the interior. The other four node modules are docked at the center of the Hab and Lab modules, two on each side, creating another “racetrack” translation path perpendicular to the first one. The DSSV has a generally vertical orientation. Within the Hab and Lab modules the decks are perpendicular to the longitudinal axis. The node modules are oriented with their longitudinal axes perpendicular to those of the Hab and Lab modules, so the nodes have a horizontal internal orientation. Decks are numbered one through nine, with decks two and eight representing the “top” and “bottom” domes of the Hab and Lab modules. The Hab Module decks encompass the ship’s galley, crew quarters, waste, hygiene, portions of crew exercise, and group recreation. The Lab Module contains food production facilities, the life science lab, the ship’s infirmary, the physical science lab, and the maintenance and fabrication workshop. The Node Modules serve a number of different functions. The two Exercise Nodes are docked to deck five of both the Hab and Lab modules and contain most of the aerobic and resistive exercise devices. Also docked to deck five but on the opposite side of the Hab and Lab modules are the Observation Deck and Space Café. These two node modules provide social gathering space for small numbers of crew. The Mission Operations Node is on deck nine directly above the Galley and contains spacecraft monitoring and commanding capabilities. Docked to it is the EVA Operations Node, which contains suit maintenance and storage. (The airlock is a separate, external element docked to the EVA Operations Node.) The Subsystems Node is on deck one beneath the maintenance and fabrication workshop. Docked to it is the Stowage Node. This node does not house DSSV primary stowage, but is a staging point where stowage brought in from logistics modules can be sorted, unpacked, or repacked as needed prior to distribution to the appropriate sections of the spacecraft.

Deep Space Habitat↗

Soyuz Landing Risk Characterization

INTRODUCTION United States On-orbit Segment (USOS) astronauts have been returning to Earth from the International Space Station (ISS) aboard the Russian Soyuz vehicle since 2003. The energetics associated with Soyuz landings are comparable to our next generation of US space vehicles (Orion, Crew Dragon, and Starliner). Soyuz also lands crew in a similar state of spaceflight deconditioning that is anticipated using the future vehicles. A number of injuries have now been observed during these Soyuz landings. Currently it is unknown how Soyuz landing accelerations and the numbers and types of associated injuries relate to the new occupant protection requirements levied upon future space vehicles. Understanding this relationship will allow better quantification of the risk of injury for future spacecraft designs. METHODS An estimate of the occurrences of injury during Soyuz landings has been determined using data contained in NASA flight medicine databases and supplemented with data collected from crewmembers and flight surgeons. At the time of this abstract, 96 USOS astronauts have returned to Earth in the Soyuz (59 US, 31 International Partner, and 6 spaceflight participants). Of that total, data from 80 individual crew landings has been collected through this study. Injuries are classified as either contusions/abrasions, Class I (minor, no medical follow-up), Class II (moderate, medical follow-up necessary), Class III (severe), and Class IV (life threatening). Current NASA occupant protection standards are primarily based on the Multi-axial Dynamic Response Index (MDRI), which is a simple lumped-parameter spring, mass, damper model tuned to human impact responses in three orthogonal axes. RESULTS Using what is known about Soyuz landing dynamics from airborne and drop test data, the MDRI predicts essentially no Class IV or III injuries, < 1% Class II injuries, and < 5% Class I injuries. Approximately one-third of all USOS crewmembers for whom data has been gathered have suffered some sort of injury related to landing in the Soyuz vehicle. Most of these injuries fall into the contusion/abrasion category (21.3 %). Class I injuries have been sustained by 5.0% of all crew, and Class II injuries have been observed in 3.8%. No Class III or IV injuries have been sustained to date. Forty-three injuries have been documented, with some crewmembers sustaining more than one injury in a single landing. The knee was the most likely area of injury (n=9), followed by the back (n=8), and then the arm (n=7). CONCLUSION Soyuz landing dynamics are the most comparable analog environment for the anticipated loads associated with Orion and commercial crew vehicle landings. The MDRI under predicts the percentage of Class II injuries experienced during Soyuz landings. The MDRI is better at predicting Class I injury occurrences. However, if contusions/abrasions are included in the minor injury Class I category, then the MDRI also vastly under predicts those injuries as well. Given this data, it will be important to re-evaluate the expected injury rates for future vehicles, and work toward improving the predictive tools for spaceflight use. Spaceflight deconditioning may account for some of the under prediction, and work is ongoing to assess the relationship between landing injuries and impact tolerance changes due to spaceflight deconditioning.

N Newby↗

Voxel Based Three-Dimensional Topology Optimization of Heat Exchanger Fins

Increasing interest in novel aircraft propulsion systems and potential for unwanted heat generation, or capture and reuse of waste heat, may require increasingly lightweight and high performing heat exchangers. Advances in manufacturing technologies have shown potential to create complex designs, but design tools need more flexibility. This study utilizes genetic algorithm-driven topology optimization to develop high performance heat exchanger fins for critical applications such as aerospace. The solid domain is generated using voxel representation where a voxel value of 1 indicates the solid domain and a voxel value of 0 refers to the fluid domain. The use of voxel representation somewhat matches the digitization of a model that is required to fabricate using additive manufacturing, and also allows for a highly unconstrained geometry. To test the topology optimization approach, a three-dimensional(3D)baseline offset strip fin geometry in steady laminar flow(Reynolds number = 215)with conjugate heat transfer(simultaneous solution of solid and fluid temperature fields)is optimized. New designs are generated using the genetic algorithm (GA) and sent to evaluation by the CFD software OpenFOAM; then the GA sorts and selects the reproduction pool for the following generation. This process is repeated for 60 generations. The study also investigates the effect off in material on the performance of the GA and the resulting designs. The results show that the optimal designs have overall performance improvement of 18% relative to the baseline. Additionally, a fin constructed of a lower conductivity material (such as an Inconel super alloy that might be necessary for waste heat recovery applications)results in lower overall performance improvement (11%)and optimal designs with higher pressure drop relative to their baseline, and relative to optimal designs produced using higher conductivity materials.

3D Optimization↗

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A Liever↗

Machine Learning Based Path Planning for Improved Rover Navigation

Enhanced AutoNav (ENav), the baseline surface navigation software for NASA’s Perseverance rover, sorts a list of candidate paths for the rover to traverse, then uses the Approximate Clearance Evaluation (ACE) algorithm to evaluate whether the most highly ranked paths are safe. ACE is crucial for maintaining the safety of the rover, but is computationally expensive. If the most promising candidates in the list of paths are all found to be infeasible, ENav must continue to search the list and run time-consuming ACE evaluations until a feasible path is found. In this paper, we present two heuristics that, given a terrain heightmap around the rover, produce cost estimates that more effectively rank the candidate paths before ACE evaluation. The first heuristic uses Sobel operators and convolution to incorporate the cost of traversing high-gradient terrain. The second heuristic uses a machine learning (ML) model to predict areas that will be deemed untraversable by ACE. We used physics simulations to collect training data for the ML model and to run Monte Carlo trials to quantify navigation performance across a variety of terrains with various slopes and rock distributions. Compared to ENav's baseline performance, integrating the heuristics can lead to a significant reduction in ACE evaluations and average computation time per planning cycle, increase path efficiency, and maintain or improve the rate of successful traverses. This strategy of targeting specific bottlenecks with ML while maintaining the original ACE safety checks provides an example of how ML can be infused into planetary science missions and other safety-critical software.

Yue, Yisong↗

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A. Liever↗