Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “planetary dynamics”

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,783 records · Page 99

Combined EDL-Mobility Planning for Planetary Missions

This paper presents an analysis framework for planetary missions that have coupled mobility and EDL (Entry-Descent-Landing) systems. Traditional systems engineering approaches to mobility missions such as MERs (Mars Exploration Rovers) and MSL (Mars Science Laboratory) independently study the EDL system and the mobility system, and does not perform explicit trade-off between them or risk minimization of the overall system. A major challenge is that EDL operation is inherently uncertain and its analysis results such as landing footprint are described using PDF (Probability Density Function). The proposed approach first builds a mobility cost-to-go map that encodes the driving cost of any point on the map to a science target location. The cost could include variety of metrics such as traverse distance, time, wheel rotation on soft soil, and closeness to hazards. It then convolves the mobility cost-to-go map with the landing PDF given by the EDL system, which provides a histogram of driving cost, which can be used to evaluate the overall risk of the mission. By capturing the coupling between EDL and mobility explicitly, this analysis framework enables quantitative tradeoff between EDL and mobility system performance, as well as the characterization of risks in a statistical way. The simulation results are presented with a realistic Mars terrain data

mission analysis↗

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

The environments experienced by vehicles upon entry into planetary atmospheres generally result in severe aerothermodynamic loading. These flows are characterized by the formation of strong shock waves, behind which high temperatures and non-equilibrium environments are generated, and the accurate prediction of internal energy relaxation and chemical kinetics becomes important. This manuscript focuses on the internal relaxation of the rotational energy of polyatomic molecules. Previous work has described the implementation of a discrete rotational energy model for diatomic molecules (linear rotators). Later, Gimelshein et. al extended this model to include polyatomic molecular internal energies. However, they did not include details on how to generate the particles at equilibrium conditions and centrifugal distortion was not included in their model. The purpose of the present study is to extend the afore-mentioned models to the treatment of discrete rotational energy in polyatomic molecules with the addition of centrifugal distortion. An overview of the background theory needed for the implementation of the model will be given followed by the model itself and verification of the model. For example, sampled versus Boltzmann distributions are compared in Figure 1 for CH3 (oblate rotator – two quantum numbers) and examples of thermal relaxation behavior for linear and non-linear rotators are presented in Figure 2. In addition, comparisons of thermophysical properties will be presented, as well as examples of adiabatic relaxation and application to hypersonic flow.

Rarefied Gas Dynamics↗

Spatiotemporal predictions of toxic urban plumes using deep learning

Industrial accidents, chemical spills, and structural fires can release large amounts of harmful materials that disperse into urban atmospheres and impact populated areas. Computer models are typically used to predict the transport of toxic plumes by solving fluid dynamical equations. However, these models can be computationally expensive due to the need for many grid cells to simulate turbulent flow and resolve individual buildings and streets. In emergency response situations, alternative methods are needed that can run quickly and adequately capture important spatiotemporal features. Here, we present a novel deep learning model called ST-GasNet inspired by the mathematical equations that govern the behavior of plumes as they disperse through the atmosphere. ST-GasNet learns the spatiotemporal dependencies from a limited set of temporal sequences of ground-level toxic urban plumes generated by a high-resolution large eddy simulation model. On independent sequences, ST-GasNet accurately predicts the late-time spatiotemporal evolution, given the early-time behavior as an input, even when a building splits a large plume into smaller plumes. By incorporating large-scale wind boundary condition information, ST-GasNet achieves a prediction accuracy of at least 90% on test data for the entire prediction period.

Civil and Environmental Engineering↗

Zero-Gravity Locomotion Simulators: New Ground-Based Analogs for Microgravity Exercise Simulation

Maintaining health and fitness in crewmembers during space missions is essential for preserving performance for mission-critical tasks. NASA's Exercise Countermeasures Project (ECP) provides space exploration exercise hardware and monitoring requirements that lead to devices that are reliable, meet medical, vehicle, and habitat constraints, and use minimal vehicle and crew resources. ECP will also develop and validate efficient exercise prescriptions that minimize daily time needed for completion of exercise yet maximize performance for mission activities. In meeting these mission goals, NASA Glenn Research Center (Cleveland, OH, USA), in collaboration with the Cleveland Clinic (Cleveland, Ohio, USA), has developed a suite of zero-gravity locomotion simulators and associated technologies to address the need for ground-based test analog capability for simulating in-flight (microgravity) and surface (partial-gravity) exercise to advance the health and safety of astronaut crews and the next generation of space explorers. Various research areas can be explored. These include improving crew comfort during exercise, and understanding joint kinematics and muscle activation pattern differences relative to external loading mechanisms. In addition, exercise protocol and hardware optimization can be investigated, along with characterizing system dynamic response and the physiological demand associated with advanced exercise device concepts and performance of critical mission tasks for Exploration class missions. Three zero-gravity locomotion simulators are currently in use and the research focus for each will be presented. All of the devices are based on a supine subject suspension system, which simulates a reduced gravity environment by completely or partially offloading the weight of the exercising test subject s body. A platform for mounting treadmill is positioned perpendicularly to the test subject. The Cleveland Clinic Zero-g Locomotion Simulator (ZLS) utilizes a pneumatic subject load device to apply a near constant gravity-replacement load to the test subject during exercise, and is currently used in conjunction with the General Clinical Research Center for evaluating exercise protocols using a bedrest analog. The enhanced ZLS (eZLS) at NASA Glenn Research Center features an offloaded treadmill that floats on a thin film of air and interfaces to a force reaction frame via variably-compliant isolators, or vibration isolation system. The isolators can be configured to simulate compliant interfaces to the vehicle, which affects mechanical loading to crewmembers during exercise, and has been used to validate system dynamic models for new countermeasures equipment designs, such as the second International Space Station treadmill slated for use in 2010. In the eZLS, the test subject and exercise device can be pitched at the appropriate angle for partial gravity simulations, such as lunar gravity (1/6th earth gravity). On both the eZLS and the NASA-Johnson Space Center standalone ZLS installed at the University of Texas Medical Branch in Galveston, Texas, USA, the subject's body weight relative to the treadmill is controlled via a linear motor subject load device (LM-SLD). The LM-SLD employs a force-feedback closed-loop control system to provide a relatively constant force to the test subject during locomotion, and is set and verified for subject safety prior to each session. Locomotion data were collected during parabolic flight and on the eZLS. The purpose was to determine the similarities and differences between locomotion in actual and simulated microgravity. Subjects attained greater amounts of hip flexion during walking and running during parabolic flight. During running, subjects had greater hip range of motion. Trunk motion was significantly less on the eZLS than during parabolic flight. Peak impact forces, loading rate, and impulse were greater on the eZLS than during parabolic while walking with a low external load (EL) and rning with a high EL. Activation timing differences existed between locations in all muscles except for the rectus femoris. The tibialis anterior and gluteus maximus were active for longer durations on the eZLS than in parabolic flight during walking. Ground reaction forces were greater with the LM-SLD than with bungees during eZLS locomotion. While the eZLS serves as a ground-based analog, researchers should be aware that subtle, but measurable, differences in kinematics and leg musculature activities exist between the environments. Aside from space applications, zero-gravity locomotion simulators may help medical researchers in the future with development of rehabilitative or therapeutic protocols for injured or ill patients. Zero-gravity locomotion simulators may be used as a ground-based test bed to support future missions for space exploration, and eventually may be used to simulate planetary locomotion in partial gravity environments, including the Moon and Mars. Figure: Zero-gravity Locomotion Simulator at the Cleveland Clinic, Cleveland, Ohio, USA

Perusek, Gail P.↗

Venus Global Reference Atmospheric Model (Venus-GRAM) Upgrades

Introduction: The Venus Global Reference Atmospheric Model (Venus-GRAM) is one of the most widely used engineering models of Venus’ atmosphere. The Venus-GRAM upgrades are being developed by NASA Marshall Space Flight Center (MSFC) and NASA Langley Research Center (LaRC). This presentation will provide details regarding the upgrades that have been made to Venus-GRAM and the ongoing objectives, tasks, and milestones related to the GRAM upgrades funded by the NASA Science Mission Directorate (SMD). Venus-GRAM: Venus-GRAM is an engineering-oriented atmospheric model that estimates mean values and statistical variations of the atmospheric properties of Venus. Venus-GRAM provides mean values and variability for any point in the atmosphere as well as seasonal, geographic, and altitude variations. Venus-GRAM outputs include atmospheric density, temperature, pressure, winds, and chemical composition along a user-defined path. It is extensively used by the engineering community because of its ability to create realistic dispersions. GRAMs have been integrated into high fidelity flight dynamic simulations of launch, entry, descent and landing (EDL), aerobraking and aerocapture. GRAMs are currently available for Earth, Mars, Venus, Neptune, Titan, and Uranus. The lower atmosphere model in Venus-GRAM (up to 250 km) is based on the Venus International Reference Atmosphere (VIRA) [1]. The Venus-GRAM thermosphere (250 to 1000 km) is based on a MSFC-developed model [2] which assumes an isothermal temperature profile initialized using VIRA conditions at 250 km [3]. The VIRA version included in Venus-GRAM includes Pioneer Venus Orbiter and Probe data as well as Venera probe data, but it does not include a solid planet model or a high-resolution gravity model [4]. Venus-GRAM Upgrade Status: Code Modernization. Venus-GRAM has been rearchitected from Fortran to a common object-oriented C++ framework called the GRAM Suite. This new architecture creates a common GRAM library of data models and utilities. The first C++ release of the rearchitected Venus-GRAM is a straight conversion from the latest Fortran version. Model Upgrades. The focus of the model upgrade task is to improve the atmosphere models in the existing GRAMs and to establish a foundation for developing GRAMs for additional destinations. The GRAM ephemeris has been upgraded to the NASA Navigation and Ancillary Information Facility (NAIF) SPICE toolkit (version N0066). The calculation of the speed of sound has also been improved in the GRAMs. In FY20, the GRAM project established contracts to improve the model data within Venus-GRAM. Hampton University is developing an empirical global model for Venus. The University of Wisconsin is reanalyzing the Venus Express radio occultation observations and analyzing the Akatsuki thermal imaging data. Upgraded Venus-GRAM Release. GRAM Suite Version 1.3 will be released in September 2021 and will contain the rearchitected Venus-GRAM, including the common GRAM framework and planet–specific code. A User Guide and Programmer’s Manual are released with all GRAMs. Conclusions: GRAMs are frequently used toolsets and vital in assessing effects of atmospheres on interplanetary spacecraft during the program life cycle process. Releases of the GRAM Suite, upgrades of the existing planetary GRAMs, and development of new planetary GRAMs are ongoing. Venus-GRAM atmosphere model upgrades will be included in the next phase of GRAM tasks. References: [1] Kliore, A. J. et al. (1985) ASR, 5, 11, 1-304. [2] Justh, H. L. et al. (2006) AIAA/AAS Astrodynamics Specialist Conference & Exhibit, Abstract AIAA-2006-6394. [3] Guide to Reference and Standard Atmosphere Models, BSR/AIAA G-003-2010. [4] Limaye, S. S. (2012), LPSC VEXAG Townhall Meeting. Acknowledgments: The authors gratefully acknowledge support from the NASA SMD.

atmospheric models↗

Lyα Observations of Comet C/2013 A1 (Siding Spring) Using MAVEN IUVS Echelle

The close approach of comet C/2013 A1 (Siding Spring) to Mars in 2014 October provided a unique opportunity to observe a dynamically new Oort cloud comet with potential for interaction with a planet’s atmosphere. The water-originating hydrogen coma of the comet extended to over 20 million km from the nucleus. Determining the properties of this coma contributes to characterizing the comet’s water content and production rate. The present study analyzes a unique data set of high spectral resolution UV observations of comet C/2013 A1 Siding Spring measured by the Mars Atmosphere and Volatile Evolution spacecraft. The Siding Spring observations capture Lyα emissions from the Martian corona, the interplanetary medium, as well as the cometary H and D reservoirs. The isolated cometary spectra are analyzed to reveal a velocity distribution of H atoms that are consistent with model estimates of H{sub 2}O photodissociated H emissions and of OH photodissociated H emissions, Doppler shifted from the main comet H emission line center by 18 km s{sup −1} and 8 km s{sup −1}, respectively. The variations in comet H brightness with distance from the nucleus are used to constrain cometary water production to a rate of 0.5 × 10{sup 28} molecules s{sup −1} at a time when Siding Spring was at 1.5 au, pre-perihelion.

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical modelling of impact seismic sources using the stress glut theory

SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantification of In-flight Physical Changes: Anthropometry and Neutral Body Posture

Currently, NASA does not have sufficient in-flight anthropometric data gathered to assess the impact of physical body shape and size changes on suit sizing. For developing future planetary and reduced gravity suits, NASA needs to quantify the impacts of microgravity on anthropometry, body posture, and neutral body postures (NBP) to ensure optimal crew performance, fit, and comfort. To obtain these impacts, anthropometric data, circumference, length, height, breadth, and depth for body segments (i.e. chest, waist, bicep, thigh, calf) from astronauts for pre, in-, and postflight conditions needs to be collected. Once this data has been collected, a comparison between pre, in-, and postflight anthropometric values will be analyzed, yielding microgravity factors. The NBP will be used to determined body posture (joint angle) changes between subjects throughout the duration of a mission. Data collection, starting with Increments 37/38, is still in progress with the completion of 3 out of 12 subjects. NASA suit engineers and NASA's Extravehicular Activity (EVA) Project Office have identified that suit fit in microgravity could become an issue. It has been noted that crewmembers often need to adjust their suit sizing once they are in orbit. This adjustment could be due to microgravity effects on anthropometry and postural changes, and is necessary to ensure optimal crew performance, fit, and comfort in space. To date, the only data collected to determine the effects of microgravity on physical human changes have been during Skylab, STS-57, and a recent HRP study on seated height changes due to spinal elongation (Spinal Elongation, Master Task List [MTL] #221). The Skylab and the STS-57 studies found that there is a distinct neutral body posture (NBP) based on photographs. The still photographs showed that there is a distinguishable posture with the arms raised and the shoulder abducted; and, in addition, the knees were flexed with noticeable hip flexion and the foot plantar flexed [1,2]. This is the one standard set of body joint angles for a NBP in microgravity. A recent simulated microgravity NBP study [3] has shown an individual variability and inconsistencies in defining NBP. This variation may be influenced by spinal growth, the type of suit fit, and other potential anthropometry factors such as spinal curvature, age, and gender. The variation aspect of this essential data is required for all kinds of space device designs (e.g. suits, habitat, mobility aids, etc.). The method proposed considers the dynamic nature of body movement and will use a measurement technique to continually monitor posture and develop a probability likelihood of the natural posture and how the NBP postures are affected by anthropometry. Additionally, Skylab studies found that crewmembers experienced a stature growth of up to 3%. The data included 3 crewmembers that showed that there is a bi-phasic stature growth once the crew enters into weightlessness. However, the Spinal Elongation study identified that the crewmembers could experience about a 6% growth in seated height and a 3% stature growth, when exposed to microgravity. The results prove that not all anthropometric measurements have the same microgravity percent growth factor. For EVA and suit engineers to properly update the sizing protocol for microgravity, they need additional anthropometric data from space missions. Hence, this study is aimed to gather additional in-flight anthropometric measurements, such as length, depth, breadth, and circumference, to determine the changes to body shape and size due to microgravity effects. It is anticipated that by recording the potential changes to body shape and size, a better suit sizing protocol will be developed for ISS and other space missions. In essence, this study will help NASA quantify the impacts of microgravity on anthropometry to ensure optimal crew performance, fit, and comfort. This study will use simplistic data collection techniques, 3D laser scanning, digital still, and video data, and perform photogrammetric analyses to determine the changes that occur to the body shape, size, and NBP when exposed to a microgravity environment.

Young, K. S.↗

Toward Improved Land Surface Initialization in Support of Regional WRF Forecasts at the Kenya Meteorological Service (KMS)

Flooding and drought are two key forecasting challenges for the Kenya Meteorological Service (KMS). Atmospheric processes leading to excessive precipitation and/or prolonged drought can be quite sensitive to the state of the land surface, which interacts with the planetary boundary layer (PBL) of the atmosphere providing a source of heat and moisture. The development and evolution of precipitation systems are affected by heat and moisture fluxes from the land surface, particularly within weakly-sheared environments such as in the tropics and sub-tropics. These heat and moisture fluxes during the day can be strongly influenced by land cover, vegetation, and soil moisture content. Therefore, it is important to represent the land surface state as accurately as possible in land surface and numerical weather prediction (NWP) models. Enhanced regional modeling capabilities have the potential to improve forecast guidance in support of daily operations and high-impact weather over eastern Africa. KMS currently runs a configuration of the Weather Research and Forecasting (WRF) NWP model in real time to support its daily forecasting operations, making use of the NOAA/National Weather Service (NWS) Science and Training Resource Center's Environmental Modeling System (EMS) to manage and produce the KMS-WRF runs on a regional grid over eastern Africa. Two organizations at the NASA Marshall Space Flight Center in Huntsville, AL, SERVIR and the Shortterm Prediction Research and Transition (SPoRT) Center, have established a working partnership with KMS for enhancing its regional modeling capabilities through new datasets and tools. To accomplish this goal, SPoRT and SERVIR is providing enhanced, experimental land surface initialization datasets and model verification capabilities to KMS as part of this collaboration. To produce a land-surface initialization more consistent with the resolution of the KMS-WRF runs, the NASA Land Information System (LIS) is run at a comparable resolution to provide real-time, daily soil initialization data in place of data interpolated from the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS) model soil moisture and temperature fields. Additionally, realtime green vegetation fraction (GVF) data from the Visible Infrared Imaging Radiometer Suite (VIIRS) on the Suomi National Polar-orbiting Partnership (Suomi- NPP) satellite will be incorporated into the KMS-WRF runs, once it becomes publicly available from the National Environmental Satellite Data and Information Service (NESDIS). Finally, model verification capabilities will be transitioned to KMS using the Model Evaluation Tools (MET; Brown et al. 2009) package in conjunction with a dynamic scripting package developed by SPoRT (Zavodsky et al. 2014), to help quantify possible improvements in simulated temperature, moisture and precipitation resulting from the experimental land surface initialization. Furthermore, the transition of these MET tools will enable KMS to monitor model forecast accuracy in near real time. This paper presents preliminary efforts to improve land surface model initialization over eastern Africa in support of operations at KMS. The remainder of this extended abstract is organized as follows: The collaborating organizations involved in the project are described in Section 2; background information on LIS and the configuration for eastern Africa is presented in Section 3; the WRF configuration used in this modeling experiment is described in Section 4; sample experimental WRF output with and without LIS initialization data are given in Section 5; a summary is given in Section 6 followed by acknowledgements and references.

Case, Johnathan L.↗

Large Vehicle Lunar Landing Surface Interaction and In-Situ Resource Based Risk Mitigation: Landing & Launch Pads

A key capability required for the exploration of planetary bodies is the ability to land on the surface. Previous work performed by NASA and other institutions has primarily focused on landing small spacecraft on planetary surfaces and the associated small-to-medium thrusters required for the soft landing. In the case of human exploration—particularly the establishment of long duration exploration and habitation outposts—the ability to land large landers, such as the SpaceX Starship, is necessary. These larger landing systems require the use of more powerful engines, with higher engine exhaust temperatures and higher landing loads. Understanding the excavation of material by the engines, as well as the potential for the landing legs to sink into the subsurface, is key in ensuring reliable and safe landings. A further improvement in landing reliability can be achieved by constructing landing / launch pads, especially with in-situ resources. Some material excavation by the plume is inevitable, leaving at least a portion of the surface scoured and uneven under the lander and ejecting regolith particles and rocks at very high velocities. One possible solution would be to robotically build landing / launch pads (ideally autonomously) at the destination using in-situ materials. In this case, the first one or few landers will need to land on unimproved surfaces at higher risk; however, they would bring the required equipment to build the landing pads with mostly local resources, thus increasing the reliability of safe landing for subsequent larger landers. A number of methods to build in-situ landing and launch pads have already been developed. These methods include no, or some, addition of required binder additives to the local regolith material, different processing approaches and result in varying landing pad strengths. A sub-scale rocket engine plume, was used to simulate some of the conditions of a landing on the Moon to assess the effectiveness of various materials for an in-situ built landing pad, The GO2/GCH4rocket engine fired on a 1m2area coupons of representative pad materials. The results will allow continued development towards materials that satisfy the landing pad properties required for the effective risk reduction and increased reliability for landing people and equipment on the lunar surface. This work contained two parts: (1) computer modeling of a large rocket engine plume interacting with regolith on the Moon, using the Granular Gas Flow Solver (GGFS) provided by CFD Research Corporation as well as other computational fluid dynamics codes (CFD) such as Loci/CHEM. (2) Developing landing/launch pad materials that could be used for in-situ construction on the lunar surface in the future, to mitigate the calculated effects of a large vehicle rocket engine landing and launching on the Moon.

Lunar↗

Large Vehicle Lunar Landing Surface Interaction and In-Situ Resource Based Risk Mitigation: Landing & Launch Pads

A key capability required for the exploration of planetary bodies is the ability to land on the surface. Previous work performed by NASA and other institutions has primarily focused on landing small spacecraft on planetary surfaces and the associated small-to-medium thrusters required for the soft landing. In the case of human exploration—particularly the establishment of long duration exploration and habitation outposts—the ability to land large landers, such as the SpaceX Starship, is necessary. These larger landing systems require the use of more powerful engines, with higher engine exhaust temperatures and higher landing loads. Understanding the excavation of material by the engines, as well as the potential for the landing legs to sink into the subsurface, is key in ensuring reliable and safe landings. A further improvement in landing reliability can be achieved by constructing landing / launch pads, especially with in-situ resources. Some material excavation by the plume is inevitable, leaving at least a portion of the surface scoured and uneven under the lander and ejecting regolith particles and rocks at very high velocities. One possible solution would be to robotically build landing / launch pads (ideally autonomously) at the destination using in-situ materials. In this case, the first one or few landers will need to land on unimproved surfaces at higher risk; however, they would bring the required equipment to build the landing pads with mostly local resources, thus increasing the reliability of safe landing for subsequent larger landers. A number of methods to build in-situ landing and launch pads have already been developed. These methods include no or some addition of required binder additives to the local regolith material, different processing approaches and result in varying landing pad strengths. A sub-scale rocket engine plume, was used to simulate some of the conditions of a landing on the Moon to assess the effectiveness of various materials for an in-situ built landing pad, The GO2/GCH4 rocket engine fired on 1m2 area coupons of representative pad materials. The results will allow continued development towards materials that satisfy the landing pad properties required for the effective risk reduction and increased reliability for landing people and equipment on the lunar surface. This work contained two parts: (1) computer modeling of a large rocket engine plume interacting with regolith on the Moon, using the Granular Gas Flow Solver (GGFS) provided by CFD Research Corporation as well as other computational fluid dynamics codes (CFD) such as Loci/CHEM. (2) Developing landing/launch pad materials that could be used for in-situ construction on the lunar surface in the future, to mitigate the calculated effects of a large vehicle rocket engine landing and launching on the Moon. The resulting computed values of plume impingement surface temperature, stagnation pressure, gas velocity, shear stress and heat flux were then matched as closely as possible in the Earth’s atmosphere in a sub-scale rocket engine GO2/CH4 test which was provided by Masten Space Systems in Mojave, California. The rocket engine was mounted on a test stand with vertical translation capabilities so that the landing operations of a lander could be simulated (Figure 1).The pad materials test coupons were placed at the surface of a large bin containing simulated lunar basalt regolith and subjected to a test firing as shown in Figure 2. The results of this testing will be presented with related findings, analysis and discussion.

Plume Surface Interaction↗

A Holistic DSMC Transport Database for Re-Entry and Ablation Modeling

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity solutions of environments containing combined continuum/rarified flow. The use of DSMC, a stochastic, particle-based method, is necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally very computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the used thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. While the transport models governing CFD and DSMC simulations are starkly different, they can effectively be linked via their use of reduced Chapman-Enskog collision integrals. In CFD, these integrals are typically stored as fitted polynomial expressions and used to directly compute gas transport properties via mixing rules or the full Chapman-Enskog formulation. In DSMC, they can be used to derive the collision parameters needed for the phenomenological collision cross-section models that govern particle interactions, via a Nelder-Mead optimization scheme. The goal of this work is to provide a unified DSMC transport database encompassing the vast majority of known gas species encountered during atmospheric entry, on Earth or any other Solar body. This goal is largely possible due to recently performed ab-initio quantum chemistry calculations. Combined with other high-fidelity literature sources, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 17000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute Variable Soft Sphere (VSS) collision model parameters for DSMC, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database will be able to be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below. Equilibrium mixture transport properties for a 19-species Titan atmospheric model [4] are computed using both fitted VSS parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches is less than 5% for the entire temperature range.

M R Gosma↗

Quantification of ln-Flight Physical Changes: Anthropometry and Neutral Body Posture

Currently, NASA does not have sufficient in-flight anthropometric data to assess the impact of changes in body shape and size. For developing future planetary and reduced-gravity suits, NASA needs to quantify the impacts of microgravity on anthropometry and body posture to ensure optimal crew performance, fit, and comfort. To obtain data on these changes, circumference, length, height, breadth, and depth for body segments (chest, waist, bicep, thigh, calf) from astronauts for preflight, in-flight, and post-flight conditions needs to be collected. Once these data have been collected, pre-flight, in-flight, and post-flight anthropometric values will be compared, yielding microgravity factors. The neutral body posture (NBP) will also be measured, to determine body posture (joint angle) changes between subjects throughout the duration of a mission. Data collection, starting with Increments 37/38, is still in progress but has been completed for 6 out of 9 subjects. NASA suit engineers and NASA's Extravehicular Activity (EVA) Project Office have identified that suit fit in microgravity could become an issue. It has been noted that crewmembers often need to adjust their suit sizing once they are in orbit. This adjustment could be due to microgravity effects on anthropometry and postural changes, and is necessary to ensure optimal crew performance, fit, and comfort in space. To date, the only data collected to determine the effects of microgravity on physical human changes were collected during Skylab 4, the Apollo-Soyuz Test Project (ASTP), Space Shuttle mission STS-57, and a recent HRP study on seated height changes due to spinal elongation (Spinal Elongation, Master Task List [MTL] #221). The Skylab 4, ASTP, and the STS-57 studies found that, according to photographs, a distinct NBP exists. The still photographs showed a distinguishable posture with the arms raised and the shoulders abducted; in addition, the knees are flexed, with noticeable hip flexion, and the foot is plantar flexed [1,2]. This combination is considered to be the standard set of body joint angles for an NBP in microgravity. A recent simulated microgravity NBP study [3] showed individual variability and inconsistencies in defining NBP. This variation may be influenced by spinal growth, and other potential anthropometry factors such as spinal curvature, age, and gender. Data on the variation of this posture data is required for all kinds of space device designs (such as suits, habitat, and mobility aids). The method proposed in this study considers the dynamic nature of body movement and will use a measurement technique to continually monitor posture and develop a probability likelihood of the neutral posture and how the NBP postures are affected by anthropometry. Additionally, Skylab studies found that crewmembers experienced a stature growth of up to 3 percent. The data included 3 crewmembers and showed that a biphasic stature growth occurs once the crewmember enters into weightlessness. However, the HRP Spinal Elongation study showed that crewmembers could experience about a 6 percent growth in seated height and a 3 percent stature growth, when exposed to microgravity. The results of that study prove that not all anthropometric measurements have the same microgravity percent growth factor. For EVA and suit engineers to properly update the sizing protocol for microgravity, they need additional anthropometric data from space missions. Hence, this study is aimed at gathering additional in-flight anthropometric measurements, such as length, depth, breadth, and circumference, to determine the changes to body shape and size caused by microgravity effects. It is anticipated that by recording the potential changes to body shape and size, NASA will develop a better suit sizing protocol for the International Space Station and other space missions. In essence, this study will help NASA quantify the impacts of microgravity on anthropometry to ensure optimal crew performance, fit, and comfort. This study will use standard anthropometry data collection techniques, 3D laser scanning, digital still photography, and video data, and perform photogrammetric analyses to determine the changes that occur to the body shape and size, and to NBP, when the human body is exposed to a microgravity environment.

Young, K. S.↗

Freeze-induced crystallization: An overlooked pathway for mineral genesis in natural waters

Natural ice plays salient roles in making the Earth habitable and sustainable. Previously overlooked, its role in chemical processes is now of emerging interest, particularly due to the freeze concentration effect, which can substantially promote chemical reactions during ice formation. We demonstrate here that ice formation can serve as a dynamic and unique pathway for mineral genesis. Freezing solutions containing dissolved manganese and carbonates produced rhodochrosite (Mn II CO 3 ) even under slightly undersaturated conditions. At room temperature, by contrast, this occurred when the solution saturation level was increased by ca. 30,000 times. The cryogenic rhodochrosite formed spherical aggregates of nano-polycrystallites, distinctly different from the cubic monocrystalline particles observed at room temperature. The distinct feature likely resulted from the combined effects of the intensified supersaturation induced by the freeze concentration effect and the low temperatures within liquid-like layers, conditions that make liquid-like layers an exceptional environment for mineral genesis, unlike typical natural water systems. The cryogenic rhodochrosite formation was successfully demonstrated using in situ, real-time X-ray absorption spectroscopy (XAS), enabling direct observation of freeze-induced solid formation. In conclusion, our findings reveal that freeze-induced crystallization may be an active mineralization pathway, potentially influencing elemental cycles within the cryosphere and contributing to minerals with distinguishing properties and reactivities in the environment.

Earth, Atmospheric, and Planetary Sciences↗

Aeolian Sand Transport in the Planetary Context: Respective Roles of Aerodynamic and Bed-Dilatancy Thresholds

The traditional view of aeolian sand transport generally estimates flux from the perspective of aerodynamic forces creating the airborne grain population, although it has been recognized that "reptation" causes a significant part of the total airborne flux; reptation involves both ballistic injection of grains into the air stream by the impact of saltating grains as well as the "nudging" of surface grains into a creeping motion. Whilst aerodynamic forces may initiate sand motion, it is proposed here that within a fully-matured grain cloud, flux is actually governed by two thresholds: an aerodynamic threshold, and a bed-dilatancy threshold. It is the latter which controls the reptation population, and its significance increases proportionally with transport energy. Because we only have experience with terrestrial sand transport, extrapolations of aeolian theory to Mars and Venus have adjusted only the aerodynamic factor, taking gravitational forces and atmospheric density as the prime variables in the aerodynamic equations, but neglecting reptation. The basis for our perspective on the importance of reptation and bed dilatancy is a set of experiments that were designed to simulate sand transport across the surface of a martian dune. Using a modified sporting crossbow in which a sand-impelling sabot replaced the bolt-firing mechanism, individual grains of sand were fired at loose sand targets with glancing angles typical of saltation impact; grains were projected at about 80 m/s to simulate velocities commensurate with those predicted for extreme martian aeolian conditions. The sabot impelling method permitted study of individual impacts without the masking effect of bed mobilization encountered in wind-tunnel studies. At these martian impact velocities, grains produced small craters formed by the ejection of several hundred grains from the bed. Unexpectedly, the craters were not elongated, despite glancing impact; the craters were very close to circular in planform. High-speed photography showed them to grow in both diameter and depth after the impactor had ricochetted from the crater site. The delayed response of the bed was "explosive" in nature, and created a miniature ejecta curtain spreading upward and outward for many centimeters for impact of 100-300 micron-diameter grains into similar material. Elastic energy deposited in the bed by the impacting grain creates a subsurface stress regime or "quasi-Boussinesq" compression field. Elastic recovery of the bed occurs by dilatancy; shear stresses suddenly convert the grains from closed to open packing, and grains are consequently able to eject themselves forcefully from the impact site. Random jostling of the grains causes radial homogenization of stress vectors and a resulting circular crater. There is a great temptation to draw parallels with cratering produced by meteorite impacts, but a rigorous search for common modelling ground between the two phenomena has not been conducted at this time. For every impact of an aerodynamically energized grain, there are several hundred grains ejected into the wind for the high-energy transport that might occur on Mars. Many of these grains will themselves become subject to the boundary layer's aerodynamic lift forces (their motion will not immediately die and add to the creep population), and these grains will become indistinguishable from those lifted entirely by aerodynamic forces. As each grain impacts the bed, it will eject even more grains into the flow. A cascading effect will take place, but because it must be finite in its growth, damping will occur as the number of grains set in motion causes mid-air collisions that prevent much of the impact energy from reaching the surface of the bed -thus creating a dynamic equilibrium in a high-density saltation cloud. It is apparent that for a given impact energy, the stress field permits a smaller volume of grains to convert to open packing as the size of the bed grains increases, or as the energy of the "percussive" grain decreases (by decrease in velocity or mass). Thus, the mass of the "repercussive" grain population that is ejected from the impact site becomes a function of the scale of the stress field in relation to the scale of the bed material (self-similarity being applicable if both bed size and energy are simultaneously adjusted). In other words, in a very high energy aeolian system where an aerodynamically raised grain can ballistically raise many more grains, the amount of material lifted into the wind becomes largely a function of a dilatancy threshold. If this threshold is exceeded, grains are repercussively injected into the saltation cloud. The "dilatancy threshold" may be defined in terms of the saltation percussive force required to convert the bed, through elastic response, from a closed to an open packing system. If open packing cannot be created, the grains cannot escape from the impact site, even though the elastic deformation and percussive force may be able to reorganize the grains with respect to one another. As the crossbow experiments showed, for an ever-increasing bed grain size, a point is reached when no material can be moved because the energy of the percussive grain is insufficient to dilate the relatively coarse bed. Although this seems to be stating the obvious -- that too little energy will not cause the bed to splash -- the consequences of exceeding the "splash threshold" by dilatancy are not so obvious for high-energy aeolian transport. It is noted that the force required to elastically dilate the bed has to overcome Coulombic grain attractions such as dipole-dipole coupling, dielectric, monopole, contact-induced dipole attractions, van der Waals forces, molecular monolayer capillary forces, as well as the mechanical interlocking frictional resistance of the grains. On Mars, it is predicted that the dilatancy threshold may be the prime control of grain flux. Additional information is contained in the original.

Marshall, J. R.↗

Effects of Lower Troposphere Vertical Mixing on Simulated Clouds and Precipitation Over the Amazon During the Wet Season

Planetary boundary layer (PBL) schemes parameterize unresolved turbulent mixing within the PBL and free troposphere (FT). Previous studies reported that precipitation simulation over the Amazon in South America is quite sensitive to PBL schemes and the exact relationship between the turbulent mixing and precipitation processes is, however, not disentangled. In this study, regional climate simulations over the Amazon in January–February 2019 are examined at process level to understand the precipitation sensitivity to PBL scheme. The focus is on two PBL schemes, the Yonsei University (YSU) scheme, and the asymmetric convective model v2 (ACM2) scheme, which show the largest difference in the simulated precipitation. During daytime, while the FT clouds simulated by YSU dissipate, clouds simulated by ACM2 maintain because of enhanced moisture supply due to the enhanced vertical moisture relay transport process: (a) vertical mixing within PBL transports surface moisture to the PBL top, and (b) FT mixing feeds the moisture into the FT cloud deck. Due to the thick cloud deck over Amazon simulated by ACM2, surface radiative heating is reduced and consequently the convective available potential energy is reduced. As a result, precipitation is weaker from ACM2. Two key parameters dictating the vertical mixing are identified, p , an exponent determining boundary layer mixing and λ , a scale dictating FT mixing. Sensitivity simulations with altered p , λ , and other treatments within YSU and ACM2 confirm the precipitation sensitivity. The FT mixing in the presence of clouds appears most critical to explain the sensitivity between YSU and ACM2.

54 ENVIRONMENTAL SCIENCES↗

Why Atens Enjoy Enhanced Accessibility For Human Space Flight

In the context of human space flight (HSF), the concept of near-Earth object (NEO) accessibility is highly subjective. Whether or not a particular NEO is accessible critically depends on mass, performance, and reliability of interplanetary HSF systems yet to be designed. Such systems would certainly include propulsion and crew life support with adequate shielding from both solar flares and galactic cosmic radiation. Equally critical architecture options are relevant to NEO accessibility. These options are also far from being determined and include the number of launches supporting an HSF mission, together with whether or not consumables are to be pre-emplaced at the destination. Until the unknowns of HSF to NEOs come into clearer focus, the notion of relative accessibility is of great utility. Imagine a group of NEOs, each with nearly equal HSF merit determined from their individual characteristics relating to crew safety, scientific return, resource utilization, and planetary defense. The more accessible members of this group are more likely to be explored first. A highly accessible NEO could conceivably be deferred in favor of a less accessible HSF destination because the latter is more accessible during a programmatically desirable launch compliant mission trajectory solutions detected in association with a specific NEO. The known NEO population is then surveyed to illustrate in which regions of heliocentric semi-major axis, eccentricity, and inclination (a, e, i) space NEOs with large n values are mapped. The (a, e, i) mapping is also formatted such that membership in each of four NEO orbit classifications, as defined below, is evident. Amors have orbits everywhere superior to (outside of) Earth's. An Amor is therefore defined to have perihelion between 1.017 astronomical units (AU) and the maximum NEO value of 1.3 AU. As of 0 hrs Universal Time on 1 January 2011 (UT epoch 2011.0), Amors numbered 2855 in the Jet Propulsion Laboratory (JPL) Small-Body Database (SBDB), comprising 37.7% of known NEOs. Apollos have orbits crossing Earth's with periods greater than Earth's. An Apollo is therefore defined to have perihelion less than 1.017 AU and a greater than 1.0 AU. As of 2011.0 UT, Apollos numbered 4080 in the SBDB, comprising 53.9% of known NEOs. Atens have orbits crossing Earth's with periods less than Earth's. An Aten is therefore defined to have aphelion greater than 0.983 AU and a less than 1.0 AU. As of 2011.0 UT, Atens numbered 618 in the SBDB, comprising 8.2% of known NEOs. Atiras have orbits everywhere inferior to (inside of) Earth's. An Atira is therefore defined to have aphelion less than 0.983 AU. As of 2011.0 UT, Atiras numbered 11 in the SBDB, comprising 0.1% of known NEOs. It is no surprise that the largest n values are chiefly associated with Apollos and Atens. Because these orbits cross Earth's, distance to be covered in a given round trip mission time delta-t can be far less than is possible for Amors or Atiras . This delta-t or the sum of mission propulsive impulse magnitudes delta-v can more frequently be minimized to enhance NHATS compliance for Apollos and Atens than is generally the case for Amors and Atiras. A less intuitive trend in NHATS results is that Atens nearly outnumber the more numerous Apollos among the most compliant NEOs as measured by n. This trend is completely out of proportion to the degree Atens are represented among the known NEO population. A theory based on geocentric NEO dynamics is presented by this paper to explain why Atens enjoy inherently greater accessibility than do Apollos. Another trend evident from mapping into (a, e, i) space is the dearth of known NEOs at low e when a < 1 AU. Underrepresentation of Atens and Atiras in the NEO catalog is at least in part attributable to observing exclusively from a perspective near Earth. Generally inferior Aten and Atira orbits are rarely, if ever, in Earth's night sky. Until a comprehensive NEO survey is conducted from an appropriateegion remote from Earth, the theory developed in this paper indicates a substantial fraction of the most accessible NEOs will remain unknown. season. Such a season is really yet another undetermined HSF architecture option. A launch season's duration will likely be measured in weeks, and it will be utilized at an indeterminate point almost certainly more than a decade in the future when HSF programmatic maturity is sufficient. Furthermore, current knowledge of the NEO population relevant to HSF is far from complete. In the 100-m-diameter class of greatest interest, only a few percent of the estimated NEO population is known [2, Figure 2.4]. Therefore, any known, lost, or fictitious NEO in a highly accessible orbit is a potential HSF destination of merit. Even if lost, fictitious, small, or hazardous, such a potential target (or another in a similar orbit) may ultimately prove to be an early HSF destination when the pertinent NEO population is more thoroughly catalogued and NEO orbits are more thoroughly maintained at high accuracy. This paper first reviews methodology and pertinent results from NASA-sponsored research performed in late 2010 and dubbed NEO HSF Accessible Targets Study (NHATS, pronounced as "gnats"). A useful accessibility metric developed during this study is n, the tally of NHATS-

Barbee, Brent↗

Feeding Ten Billion People Is Possible Within Four Terrestrial Planetary Boundaries

Global agriculture puts heavy pressure on planetary boundaries, posing the challenge to achieve future food security without compromising Earth system resilience. On the basis of process-detailed, spatially explicit representation of four interlinked planetary boundaries (biosphere integrity, land-system change, freshwater use, nitrogen flows) and agricultural systems in an internally consistent model framework, we here show that almost half of current global food production depends on planetary boundary transgressions. Hotspot regions, mainly in Asia, even face simultaneous transgression of multiple underlying local boundaries. If these boundaries were strictly respected, the present food system could provide a balanced diet (2,355 kcal per capita per day) for 3.4 billion people only. However, as we also demonstrate, transformation towards more sustainable production and consumption patterns could support 10.2 billion people within the planetary boundaries analysed. Key prerequisites are spatially redistributed cropland, improved water–nutrient management, food waste reduction and dietary changes. Adoption of the Sustainable Development Goals by all nations in 2015 is the first ever commitment to a world development path that safeguards the stability of the Earth system as a prerequisite for meeting universal human standards1. The longstanding challenge of achieving food security through sustainable agriculture is particularly acute in this context as world agriculture is a leading cause for the current transgressions of multiple planetary boundaries (PBs) globally and regionally2–5. The PB framework is a comprehensive scientific attempt to synoptically define our planet’s biogeophysical limits to anthropogenic interference. It suggests bounds to nine interacting processes that together delineate a Holocene-like Earth system state. The Holocene is chosen as the reference state as it is the only period known to provide a safe operating space for a world population of several billion people, and according to a precautionary principle, the PBs are set in sufficient distance from processes that may critically undermine Earth system resilience and global sustainability. A challenging question, thus, is whether human development goals such as food security can be met while maintaining multiple PBs along with their subglobal manifestations. Further PB transgressions could jeopardize the chances of providing sufficient food for a world population projected to be wealthier and reach >9 billion by 2050. This conundrum portrays a tradeoff between Earth’s biophysical carrying capacity and humankind’s rising food demand, calling in response for radical rethinking of food production and consumption patterns6–9. Yield gap closures, avoidance of excessive input use, shifts towards less resource-demanding diets, food waste reductions and efficient international trade are crucial options for sustainably increasing the food supply10–15. For example, enhancing water-use efficiency on irrigated and rain-fed farms can triple or quadruple crop yields in low-performing systems, suggesting possible global gains of >20% (ref. 16). Even higher gains appear feasible through globally optimized configurations of the land-use pattern17, and cutting food losses by half could generate food for another billion people18. Thus, collective large-scale implementation of such options could sustain food for a further growing world population19. Yet achieving this within a safe operating space as defined by PBs requires not only a halt to but actually a reversal of existing PB transgressions. Previous studies suggest that such a reconciliation might be possible, but these were based on aggregate representations of PBs (not accounting for the spatial patterns of limits, transgressions and interactions) or considered only one boundary in isolation17,20–23. Here, we systematically quantify to what extent current food production depends on local to global transgressions of the PBs for biosphere integrity, land-system change, freshwater use and nitrogen (N) flows, along with the potential of a range of solutions to avoid these transgressions and still increase food supply (Table 1). To this end, we configured an internally consistent process-based model of the terrestrial biosphere including agriculture (LPJmL) with multiple spatially distributed PBs and their interactions. LPJmL is among the longest-established and best-evaluated biosphere models, showing robust performance regarding simulation of, for example, carbon, water and crop yield dynamics (Supplementary Figs. 1 and 2 and Supplementary Table 1; see ref. 24 for a comprehensive benchmarking and Supplementary Methods for more detail on model evaluations). In principle following established definitions4, we refine the computation of some PBs with respect to their regional patterns and interactions (Methods), providing globally gridded precautionary limits to human interference with the Earth system at a level of great detail. In particular, we account for the evidence that many PBs need to be represented spatially explicitly4 to cover their

Gerten, Dieter↗