Search NASA⌕ Search

SEARCH · Search NASA

Results for “statistically designed experiments”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Dynamics and structure of turbulent premixed flames

In earlier work (Mantel & Bilger, 1994) the structure of the turbulent premixed flame was investigated using statistics based on conditional averaging with the reaction progress variable as the conditioning variable. The DNS data base of Trouve and Poinsot (1994) was used in this investigation. Attention was focused on the conditional dissipation and conditional axial velocity in the flame with a view to modeling these quantities for use in the conditional moment closure (CMC) approach to analysis of kinetics in premixed flames (Bilger, 1993). Two remarkable findings were made: there was almost no acceleration of the axial velocity in the flame front itself; and the conditional scalar dissipation remained as high, or higher, than that found in laminar premixed flames. The first finding was surprising since in laminar flames all the fluid acceleration occurs through the flame front, and this could be expected also for turbulent premixed flames at the flamelet limit. The finding gave hope of inventing a new approach to the dynamics of turbulent premixed flames through use of rapid distortion theory or an unsteady Bernoulli equation. This could lead to a new second order closure for turbulent premixed flames. The second finding was contrary to our measurements with laser diagnostics in lean hydrocarbon flames where it is found that conditional scalar dissipation drops dramatically below that for laminar flamelets when the turbulence intensity becomes high. Such behavior was not explainable with a one-step kinetic model, even at non-unity Lewis number. It could be due to depletion of H2 from the reaction zone by preferential diffusion. The capacity of the flame to generate radicals is critically dependent on the levels of H2 present (Bilger, et al., 1991). It seemed that a DNS computation with a multistep reduced mechanism would be worthwhile if a way could be found to make this feasible. Truly innovative approaches to complex problems often come only when there is the opportunity to work close at hand with the (in this case numerical) experimental data. Not only can one spot patterns and relationships in the data which could be important, but one can also get to know the limitations of the technique being used, so that when the next experiment is being designed it will address resolvable questions. A three-year grant from the Australian Research Council has enabled us to develop a small capability at the University of Sydney to work on DNS of turbulent reacting flow, and to analyze data bases generated at CTR. Collaboration between the University of Sydney and CTR is essential to this project and finding a workable modus operandum for this collaboration, given the constraints involved, has been a major objective of the past year's effort. The overall objectives of the project are: (1) to obtain a quantitative understanding of the dynamics of turbulent premixed flames at high turbulence levels with a view to developing improved second order closure models; and (2) to carry out new DNS experiments on turbulent premixed flames using a carefully chosen multistep reduced mechanism for the chemical kinetics, with a view to elucidating the laser diagnostic findings that are contrary to the findings for DNS using one-step kinetics. In this first year the objectives have been to make the existing CTR data base more accessible to coworkers at the University of Sydney, to make progress on understanding the dynamics of the flame in this existing CTR data base, and to carefully construct a suitable multistep reduced mechanism for use in a new set of DNS experiments on turbulent premixed flames.

Bilger, R. W.↗
DOE OSTI2025

Analysis and optimization of seismic monitoring networks with Bayesian optimal experimental design

SUMMARY Monitoring networks increasingly aim to assimilate data from a large number of diverse sensors covering many sensing modalities. Bayesian optimal experimental design (OED) seeks to identify data, sensor configurations or experiments which can optimally reduce uncertainty and hence increase the performance of a monitoring network. Information theory guides OED by formulating the choice of experiment or sensor placement as an optimization problem that maximizes the expected information gain (EIG) about quantities of interest given prior knowledge and models of expected observation data. Therefore, within the context of seismo-acoustic monitoring, we can use Bayesian OED to configure sensor networks by choosing sensor locations, types and fidelity in order to improve our ability to identify and locate seismic sources. In this work, we develop the framework necessary to use Bayesian OED to optimize a sensor network’s ability to locate seismic events from arrival time data of detected seismic phases at the regional-scale. This framework requires five elements: (i) A likelihood function that describes the distribution of detection and traveltime data from the sensor network, (ii) A prior distribution that describes a priori belief about seismic events, (iii) A Bayesian solver that uses a prior and likelihood to identify the posterior distribution of seismic events given the data, (iv) An algorithm to compute EIG about seismic events over a data set of hypothetical prior events, (v) An optimizer that finds a sensor network which maximizes EIG. Once we have developed this framework, we explore many relevant questions to monitoring such as: how to trade off sensor fidelity and earth model uncertainty; how sensor types, number and locations influence uncertainty; and how prior models and constraints influence sensor placement.

58 GEOSCIENCES↗
DOE OSTI2026

Ripening of Rh Nanoparticle Catalysts in Reverse Water–Gas Shift via a Data-Driven Model Combining Physics, Theory, and Experiment

Degradation via sintering is an ongoing challenge that impedes the broad commercial success of supported metallic nanoparticle catalysts. To mitigate degradation via informed catalyst design and process operations, here we aim to disambiguate the underlying mechanisms of sintering by combining theory and experiment in a quantitative framework. While mechanistic sintering models exist, they only model a single sintering pathway, even though multiple sintering mechanisms can occur simultaneously or dominate at different stages of the process. Data-driven machine learning models have emerged as a means to represent complex processes through data regression. However, machine learning models have very large data needs and lack mechanistic insights due to their black-box encoding. To develop an interpretive model of catalyst degradation via sintering, we constructed a hybrid model combining mechanistic “physics-based” models and data-driven methods to obtain both reliable predictions and mechanistic insights regarding experimentally observed sintering phenomena. Focusing on nanoparticle sintering in the Rh–TiO 2 catalyst for the reverse water–gas shift (RWGS) reaction, the hybrid model couples a mechanistic term for Ostwald ripening with energy values calculated via density functional theory (DFT) with a parametric, data-driven discrepancy function term for unmodeled mechanisms. The hybrid model is trained using Bayesian inference with data collected from small-angle X-ray scattering (SAXS) in situ experiments wherein average nanoparticle diameter versus time was measured at three relevant operating temperatures. The calibrated hybrid model results show that an Ostwald ripening-only model parameterized with fixed DFT energies does not fully capture the time and temperature dependence of the SAXS-observed sintering kinetics, and that an additional functional contribution, or DFT energy calibration, is required to reconcile simulation and experiment. Analysis of the hybrid-model error confirms that the hybrid model outperforms both the purely mechanistic and purely data-driven alternatives in terms of expected predictive accuracy for time-evolving average particle sizes. Furthermore, the results support the hypothesis that the Ostwald ripening mechanism is less important for explaining the sintering phenomena as operating temperature increases under an assumed fixed DFT parameterization. This could be explained in one of two ways: either latent, unmodeled sintering mechanisms dominate at higher temperatures, or the DFT uncertainty increases with temperature. The proposed modeling approach directly links theory to experiments and simulations via a statistical hybrid modeling framework and can be extended to other catalytic systems to improve predictive models and mechanistic understanding.

Bayesian hybrid modeling↗
NASA NTRS1998

The NASA Low-Pressure Turbine Flow Physics Program

An overview of the NASA Lewis Low-Pressure Turbine (LPT) Flow Physics Program will be presented. The program was established in response to the aero-engine industry's need for improved LPT efficiency and designs. Modern jet engines have four to seven LPT stages, significantly contributing to engine weight. In addition, there is a significant efficiency degradation between takeoff and cruise conditions, of up to 2 points. Reducing the weight and part count of the LPT and minimizing the efficiency degradation will translate into fuel savings. Accurate prediction methods of LPT flows and losses are needed to accomplish those improvements. The flow in LPT passages is at low Reynolds number, and is dominated by interplay of three basic mechanisms: transition, separation and wake interaction. The affecting parameters traditionally considered are Reynolds number, freestream turbulence intensity, wake frequency parameter, and the pressure distribution (loading). Three-dimensional effects and additional parameters, particularly turbulence characteristics like length scales, spectra and other statistics, as well as wake turbulence intensity and properties also play a role. The flow of most interest is on the suction surface, where large losses are generated as the flow tends to separate at the low Reynolds numbers. Ignoring wakes, a common flow scenario, there is laminar separation, followed by transition on the separation bubble and turbulent reattachment. If transition starts earlier the separation will be eliminated and the boundary layer will be attached leading to the well known bypass transition issues. In contrast, transition over a separation bubble is closer to free shear layer transition and was not investigated as well, particularly in the turbine environment. Unsteadiness created by wakes complicates the picture. Wakes induce earlier transition, and the calmed regions trailing the induced turbulent spots can delay or eliminate separation via shear stress modification. Three-dimensional flow physics and geometry will have strong effects. Altogether a very complex and challenging problem emerges. The objective of the program is to provide improved models and physical understanding of the complex flow, which are essential for accurate prediction of flow and losses in the LPT. Experimental, computational and analytical work as complementing and augmenting approaches are used. The program involves industry, universities and research institutes, and other government laboratories. It is characterized by strong interaction among participants, quick dissemination of results, and responsiveness to industry's needs. The presentation will describe the work elements. Highlighting some activities in progress are experiments on simulated blade suction surface in low-speed wind tunnels, on curved wall, and on a flat-plate, both with pressure gradient. In the area of computation, assessment of existing models is performed using RANS (Reynolds Averaged Navier Stokes) simulations. Laminar flow DNS was completed. Analytical studies of instability and receptivity in attached and separated flows were started. In the near future the program is moving to include wake effects and development of improved modeling. Experimental work in preparation stages are: (1) Addition of wakes to the curved tunnel experiment; (2) Low-speed rotating rig experiment on GE90 engine LPT; and (3) Transonic cascade. In the area of computation, it is expected to move from model assessment towards development of improved models. In addition, a new project of Large Eddy Simulation (LES) of LPT is to begin and will provide numerical data bases. It is planned to implement the emerging improved models in a multistage turbomachinery code and to validate against the GE90 engine LPT.

Ashpis, David E.↗
DOE OSTI2025

Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties

This project developed machine learning (ML) methods, lab data sets, and field data to advance geothermal exploration and geothermal energy production. The work had three focus areas. One involved the development of ML methods to use microearthquakes (MEQs) for imaging geothermal reservoir properties and improving subsurface characterization – most importantly the evolution of permeability within the evolving reservoir. This part of the work included development of ML approaches for automated MEQ location, focal mechanism determination and identification of earthquake precursors. The second area focused on using MEQ signals generated by geothermal exploration and production to predict the relationship between fluid injection and seismicity. Here, we extended to reservoir scale our success in using ML to predict laboratory earthquakes and fault zone stress state. The third focus area was on lab experiments. Here, we developed new ML models for lab earthquake prediction and identification of precursors to failure to improve earthquake forecasting and early warning in geothermal settings. Major outcomes of our work include ML models that learn from MEQ signals during geothermal exploration and production to predict induced seismicity. MEQs occur naturally in connection with drilling and energy production. We developed ML methods to use the seismic waves from these events to characterize the elastic, hydraulic and poromechanical properties of reservoirs. Our work illuminated fracture geometry and the evolution of fracture permeability by incorporating seismic coda wave analysis and ML methods to relate fluid injection and seismicity. We significantly expanded laboratory earthquake prediction to include methods that use both passive measurements of microearthquakes within the lab fault zones and also active source acoustic measurements of fault zone elastic properties. These methods can now predict fault zone stress state, time to failure and the magnitude of lab earthquakes. Our work showed that repetitive stick- slip failure events during frictional sliding (the lab equivalent of earthquakes) are preceded by a cascade of micro-failure events that radiate energy in a manner that foretells unstable failure – manifest as laboratory MEQs. We documented a mapping between fracture properties and statistical attributes of elastic radiation. We extended existing works to geothermal reservoir scale and developed ML methods to determine reservoir permeability, fracture properties, and their evolution during geothermal energy production. An attractive feature of ML algorithms is their ability to handle big datasets and reveal patterns and correlations that may remain invisible to conventional analyses. Our work connected data from field, laboratory and intermediate scales to study permeability, stress, strength, fracture stiffness and geometry. At the field scale we used data from the Newberry Volcano field site, UtahFORGE, EGS Collab, and also the Bedretto underground research lab in Switzerland. These data sets are bridging the gap between the lab scale, theory, and reservoir scale. Our work produced plain language summaries to improve public understanding of DOE research. We also developed openly distributed ML and seismicity datasets for use by all researchers and we published connections between induced seismicity in geothermal areas and reservoir properties including permeability, fracture properties, and stress state. Our models are designed for the large data sets of induced seismicity typically associated with geothermal sites. We produced labeled event catalogs and used them on geothermal data to assess how ML can facilitate geothermal production and exploration. All datasets are available on the GDR Productivity: The project produced 32 publications in peer reviewed journals (two are in review). It supported the work of 6 PhD students, 40 conference presentations, 6 keynote talks at national meetings, and mentoring and professional development for 4 postdoctoral fellows.

15 GEOTHERMAL ENERGY↗
NASA NTRS2004

Evaluation of a Multi-Parameter Sensor for Automated, Continuous Cell Culture Monitoring in Bioreactors

Compact and automated sensors are desired for assessing the health of cell cultures in biotechnology experiments in microgravity. Measurement of cell culture medium allows for the optirn.jzation of culture conditions on orbit to maximize cell growth and minimize unnecessary exchange of medium. While several discrete sensors exist to measure culture health, a multi-parameter sensor would simplify the experimental apparatus. One such sensor, the Paratrend 7, consists of three optical fibers for measuring pH, dissolved oxygen (p02), dissolved carbon dioxide (pC02) , and a thermocouple to measure temperature. The sensor bundle was designed for intra-arterial placement in clinical patients, and potentially can be used in NASA's Space Shuttle and International Space Station biotechnology program bioreactors. Methods: A Paratrend 7 sensor was placed at the outlet of a rotating-wall perfused vessel bioreactor system inoculated with BHK-21 (baby hamster kidney) cells. Cell culture medium (GTSF-2, composed of 40% minimum essential medium, 60% L-15 Leibovitz medium) was manually measured using a bench top blood gas analyzer (BGA, Ciba-Corning). Results: A Paratrend 7 sensor was used over a long-term (>120 day) cell culture experiment. The sensor was able to track changes in cell medium pH, p02, and pC02 due to the consumption of nutrients by the BHK-21. When compared to manually obtained BGA measurements, the sensor had good agreement for pH, p02, and pC02 with bias [and precision] of 0.02 [0.15], 1 mm Hg [18 mm Hg], and -4.0 mm Hg [8.0 mm Hg] respectively. The Paratrend oxygen sensor was recalibrated (offset) periodically due to drift. The bias for the raw (no offset or recalibration) oxygen measurements was 42 mm Hg [38 mm Hg]. The measured response (rise) time of the sensor was 20 +/- 4s for pH, 81 +/- 53s for pC02, 51 +/- 20s for p02. For long-term cell culture measurements, these response times are more than adequate. Based on these findings , the Paratrend sensor could offer automated, continuous monitoring of cell cultures with a temporal resolution of 1 minute, which is not attainable by sampling via handheld blood analyzer (i-STAT). Conclusion: The resulting bias and precision found in these cell culture-based studies is comparable to Paratrend sensor clinical results. Although the large error in p02 measurements (+/-18 mm Hg) may be acceptable for clinical applications, where Paratrend values are periodically adjusted to a BGA measurement, the O2 sensor in this bundle may not be reliable enough for the single-calibration requirement of sensors used in NASA's bioreactors. The pH and pC02 sensors in the bundle are reliable and stable over the measurement period, and can be used without recalibration to measure cell cultures in rn.jcrogravity biotechnology experiments. Future work will test additional Paratrend sensors to provide statistical assessment of sensor performance.

Pappas, D.↗
NASA NTRS2014

Focal Gray Matter Plasticity as a Function of Long Duration Head Down Tilted Bed Rest: Preliminary Results

Long duration spaceflight (i.e., 22 days or longer) has been associated with changes in sensorimotor systems, resulting in difficulties that astronauts experience with posture control, locomotion, and manual control. The microgravity environment is an important causal factor for spaceflight induced sensorimotor changes. Whether these sensorimotor changes are solely related to peripheral changes from reduced vestibular stimulation, body unloading, body fluid shifts or that they may be related to structural and functional brain changes is yet unknown. However, a recent study reported associations between microgravity and flattening of the posterior eye globe and protrusion of the optic nerve [1] possibly as the result of increased intracranial pressure due to microgravity induced bodily fluid shifts [3]. Moreover, elevated intracranial pressure has been related to white matter microstructural damage [2]. Thus, it is possible that spaceflight may affect brain structure and thereby cognitive functioning. Long duration head down tilt bed rest has been suggested as an exclusionary analog to study microgravity effects on the sensorimotor system [4]. Bed rest mimics microgravity in body unloading and bodily fluid shifts. In consideration of the health and performance of crewmembers both in- and post-flight, we are conducting a prospective longitudinal 70-day bed rest study as an analog to investigate the effects of microgravity on brain structure [5]. Here we present results of the first six subjects. Six subjects were assessed at 12 and 7 days before-, at 7, 30, and ~70 days in-, and at 8 and 12 days post 70 days of bed rest at the NASA bed rest facility in UTMB, Galveston, TX, USA. At each time point structural MRI scans (i.e., high resolution T1-weighted imaging and Diffusion Tensor Imaging (DTI)) were obtained using a 3T Siemens scanner. Focal changes over time in gray matter density were assessed using the voxel based morphometry 8 (VBM8) toolbox under SPM. Longitudinal processing in VBM8 includes linear registration of each scan to the mean of the subject and subsequently transforming all scans in to MNI space by applying the warp from the mean subject to MNI to the individual gray matter segmentations. Modulation was applied so that all images represented the volume of the original structure in native space. Voxel wise analysis was carried out on the gray matter images after smoothing, using a flexible factorial design with family wise error correction. Focal changes in white matter microstructural integrity were assessed using tract based spatial statistics (TBSS) as part of FMRIB software library (FSL). TBSS registers all DTI scans to standard space. It subsequently creates a study specific white matter skeleton of the major white matter tracts. For each subject, for each DTI metric (i.e. fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD)), the maximum value in a line perpendicular to the skeleton tract is projected to the skeleton. Non-parametric permutation based t-tests and ANOVA's were used for voxel-wise comparison of the skeletons. For both VBM and TBSS, comparison of pre bed rest measurements did not show significant differences. VBM analysis revealed decreased gray matter density in bilateral areas including the frontal medial cortex, the insular cortex and the caudate (see Figure) from 'pre to in bed rest'. Over the same time period, there was an increase in gray matter density in the cerebellum, occipital-, and parietal cortex, including the precuneus (see Figure). The majority of these changes did not recover from 'during to post bed rest'. TBSS analysis did not reveal significant changes in white matter microstructural integrity after correction for multiple comparisons. Uncorrected analyses (p<.015) revealed an increase in RD in the cerebellum and brainstem from pre bed rest to the first week in bed rest that did not recover post bed rest. Extended bed rest, which is an analog for microgravity, can result in gray matter changes and potentially in microstructural white matter changes in areas that are important for neuro motor behavior and cognition. These changes did not recover at two weeks post bed rest. Whether the effects of bed rest wear off at longer times post bed rest, and if they are associated with behavior are important questions that warrant further research.

Koppelmans, V.↗
DOE OSTI2022

In situ feature analysis for large-scale multiphase flow simulations

The study of multiphase flow is essential for designing chemical reactors such as fluidized bed reactors (FBR), as a detailed understanding of hydrodynamics is critical for optimizing reactor performance and stability. An FBR allows scientists to conduct different types of chemical reactions involving multiphase materials, especially interaction between gas and solids. During such complex chemical processes, the formation of void regions in the reactor, generally termed as bubbles, is an important phenomenon. The study of these bubbles has a deep implication in predicting the reactor’s overall efficiency. But physical experiments needed to understand bubble dynamics are costly and non-trivial due to the technical difficulties involved and harsh working conditions of the reactors. Therefore, to study such chemical processes and bubble dynamics, a state-of-the-art computational simulation MFIX-Exa is being developed. Despite the proven accuracy of MFIX-Exa in modeling bubbling phenomena, the large-scale output data prohibits the use of traditional post hoc analysis capabilities in both storage and I/O time. Herein, to address these issues and allow the application scientists to explore the bubble dynamics in an efficient and timely manner, we have developed an end-to-end analytics pipeline that enables in situ detection of bubbles, followed by a flexible post hoc visual exploration methodology of bubble dynamics. The proposed method enables interactive analysis of bubbles, along with quantification of several bubble characteristics, enabling experts to understand the bubble interactions in detail. Positive feedback from the experts has indicated the efficacy of the proposed approach for exploring bubble dynamics in very-large-scale multiphase flow simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - NLR Historical Wind

The U.S. Department of Energy and the National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis from variable sources, hydrogen compression and storage, and hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset represents part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure nationwide. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with various energy technologies. Future datasets will demonstrate how existing hydrogen fuel cell technologies can provide controllable, dispatchable, and variable power output for artificial intelligence (AI) data centers and other variable loads. This dataset entry describes hydrogen production by conducting a statistical analysis of historical wind data over a five-year period (2020-2025) from a single 1.5MW turbine manufactured by General Electric (GE) located at NLR’s Flatirons Campus, to generate an experimental test profile that was deployed on a 1.25-MW proton exchange membrane type MC250 electrolyzer system manufactured by Nel Hydrogen . [1] While the electrolyzer balance-of-plant supports up to 2.5 MW of electrolysis, NLR only has a single 1.25-MW electrolysis stack. The historical wind data provided several metrics, however, the analysis particularly focused on the measured power output by the wind turbine. The power output time series of data for each day was categorized by total energy generation and standard deviation, and the day that represented the highest combination of these two metrics was chosen – December 25th, 2022. This process was then repeated for a moving four-hour window within this day to identify the most statistically variable period. Finally, this four-hour period was scaled by 65% to match the 1.25 MW electrolyzer. The electrolysis system controls hydrogen production by varying DC current applied to the stack, from a maximum of 3000 A to a minimum safe operation of 300 A, or 10%. Because the current – voltage characteristic changes as the stack ages and efficiency degrades, the actual minimum safe operating power changes over time. The historical wind profiles were translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1 Hz frequency. For more details on the statistical analysis process, see the presentation labeled “ Public Reference Data for Megawatt-Scale Hydrogen Electrolysis” provided with each data entry. These datasets report relevant hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. Each .zip file represents a single wind turbine electrolysis experiment and is formatted as follows: {technology}_{scaling factor}-{electrolyzer ramp rate in amperes/second} For instance, “wind-GE1.5MW_0.65-400.zip” represents the hour-long experiment using historical data from the wind-GE1.5MW turbine, scaled to 65%, with the electrolyzer power supply set to a maximum ramp rate (gain and slew) of 400 A/s. Each .zip folder contains the following files: A .csv file containing raw data An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production, electrolysis power consumption, and wind power input. A PDF file detailing the historical wind data statistical analysis used to generate the wind profile. An experiment labeled “characterization_200.zip” demonstrates the MC250 electrolyzer steady-state response with 30-minute load steps for a total duration of 5 hours. Finally, a .csv file is provided with all simulated wind experiments combined into one dataset labeled "combined_historical_wind_experiments.csv". NLR also built an AI/machine-learning predictive model based on these datasets. The model ingests the electrolyzer current command in amperes, as well as various pressures and temperatures across the system, and predicts hydrogen output in kilograms per hour. The complete model can be found at https://huggingface.co/NatLabRockies/ptmelt-hydrogen-electrolysis [1] nelhydrogen.com/product/mc-series-electrolyser .

08 HYDROGEN↗
NASA NTRS2004

Gear Performance Improved by Coating

Gears, bearings, and other mechanical elements transmit loads through contacting surfaces. Even if properly designed, manufactured, installed, and maintained, gears and bearings will eventually fail because of the fatigue of the working surfaces. Economical means for extending the fatigue lives of gears and bearings are highly desired, and coatings offer the opportunity to engineer surfaces to extend the fatigue lives of mechanical components. A tungsten-containing diamondlike-carbon coating exhibiting high hardness, low friction, and good toughness was evaluated for application to spur gears. Fatigue testing was done at the NASA Glenn Research Center on both uncoated and coated spur gears. The results showed that the coating extended the surface fatigue lives of the gears by a factor of about 5 relative to the uncoated gears. For the experiments, a lot of spur test gears made from AISI 9310 gear steel were case-carburized and ground to aerospace specifications. The geometries of the 28-tooth, 8-pitch gears were verified as meeting American Gear Manufacturing Association (AGMA) quality class 12. One-half of the gears were randomly selected for coating. The method of coating was selected to achieve desired adherence, toughness, hardness, and low-friction characteristics. First the gears to be coated were prepared by blasting (vapor honing) with Al2O3 particles and cleaning. Then, the gears were provided with a thin adhesion layer of elemental chromium followed by magnetron sputtering of the outer coating consisting of carbon (70 at.%), hydrogen (15 at.%), tungsten (12 at.%), and nickel (3 at.%) (atomic percent at the surface). In total, the coating thickness was about 2.5 to 3 microns. As compared with the steel substrate, the coated surface was harder by a factor of about 2 and had a smaller elastic modulus. All gears were tested using a 5-centistoke synthetic oil, a 10,000-rpm rotation speed, and a hertzian contact stress of at least 1.7 GPa (250 ksi). Tests were run until either surface fatigue occurred or 300 million stress cycles were completed. Tests were run using either a pair of uncoated gears or a pair of coated gears (coated gears mated with uncoated gears were not evaluated). The fatigue test results, shown on Weibull coordinates in the graph, demonstrate that the coating provided substantially longer fatigue lives even though some of the coated gears endured larger stresses. The increase in fatigue life was a factor of about 5 and the statistical confidence for the improvement is high (greater than 99 percent). Examination of the tested gears revealed substantial reductions of total wear for coated gears in comparison to uncoated gears. The coated gear surface topography changed with running, with localized areas of the tooth surface becoming smoother with running. Theories explaining how coatings can extend gear fatigue lives are research topics for coating, tribology, and fatigue specialists. This work was done as a partnership between NASA, the U.S. Army Research Laboratory, United Technologies Research Corporation, and Sikorsky Aircraft.

Krantz, Timothy L.↗
DOE OSTI2023

A review of privacy in energy applications

As the distribution system continues to experience an increase in distributed energy resource (DER) and electric vehicle (EV) penetration, so does the need for new solutions that can help grid operators manage and leverage their capabilities. This will undoubtedly lead to new operational schemes and business opportunities that will transform the traditional consumer into a prosumer who will be more actively engaged in grid operations. Although the field is still under active development, many of the potential use cases presented in literature or industry are built upon edge computing, two-way communications, and other innovative computational constructs to attain their goals. However, at their core, many use cases assume a great level of data access to aid with the decision-making process, an assumption that may need to be revised to ensure fair and equitable operational processes are maintained. This may be particularly true as edge resources are predicted to participate in retail-side, many-to-many, or peer-to-peer markets and thus may lead to financial impacts if data access considerations are ignored. The need to revise data access mechanisms can be further justified by the introduction of new participants into the operational process, who do not have the same level of trust, nor the incentives to focus on energy delivery as their primary objective. At the same time, more consumers are becoming aware of their own data, and the potential impacts of its abuse. To help solution developers better understand these risks, this report has been developed to offer an initial introduction to the topic of privacy. This is achieved by 1) Highlighting the need for privacy-aware solutions; 2) Encouraging system designers to be inquisitive about the status quo; 3) Documenting the existing threat space; 4) Presenting and evaluating tools that may be helpful towards enabling better privacy postures; and 5) Making recommendations to encourage the adoption of better practices. From a technical perspective, the report focuses on evaluating two potential techniques by applying them to the Transactive Energy Space. Based on the obtained results, it can be established that differential privacy (DP) methods may have limited applicability when highly correlated, time-series data records need to be protected. However, DP may be a powerful tool when it is used to aggregate and analyze mid-size and large-size data sets in a more traditional statistical environment. The second tool under evaluation is threshold cryptography, which can guarantee complete secrecy (and thus privacy) but requires the establishment of key management procedures and dedicated communication channels for key coordination. Therefore, due to its increased computational overhead, the use of threshold cryptography must be weighted using a cost/benefit analysis on a per-application basis.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2018

A Scaling Study of Microbially-Enhanced Methane Production From Coal (Mecbm): Optimizing Nutrient Delivery For Maximized Methane Production

Microbially-enhanced coalbed methane (MECBM) recovery is a novel method to stimulate gas production in CBM reservoirs. We first develop suitable nutrient recipes for native coals from sites each in the Illinois Basin and in the San Juan Basin. The optimal concentrations for statistically significant parameters are identified and the methane yields at laboratory scale are obtained (Part I & II). We then demonstrate the importance of replicating truly in-situ conditions for bio- gasification studies by investigating methane productivity in a core holder simulating in-situ conditions. Based on the experimental results, the resulting gas production correlates positively with the permeability of the samples. We examine the deliverability of engineered microbial solutions to promote coal bioconversion in coals and demonstrate that greater than milli-darcy permeability is required under in-situ stress conditions (Part III). This constraint is supplemented by imaging and a fractal-based characterization approach to evaluate changes in the physical properties of coal due to bio-gasification. High resolution SEM imaging is used to observe swelling of the coal matrix with sub-micron wide continuous pores observed to persist even after microbial treatment (IV). We explore the potential of generating, propping-open and sustaining the functionality of hydraulic fractures in coal using different fracturing fluids (Part V). The alteration of the nano-pore structure due to MECBM treatment is examined by using small angle X-ray scattering (SAXS) to define the in-situ porosity evolution as a proxy to infer permeability evolution (Part VI). Coupled hydro-mechanical and bio-gasification modeling is applied to determine the effectiveness of nutrient delivery by hydraulic fracturing and its field implementation. Compared with both the unstimulated and hydro-shear-dilated cases of the natural fracture network, hydraulic fracturing demonstrates a significant improvement to enhance nutrient delivery. A preliminary analysis of field-scale methane yield is completed with consideration of the combined influence of permeability evolution, solute transport and microbial growth (Part VII). Developing nutrient solutions for Illinois Basin and San Juan Basin (Task 2): To enhance methane release from coal through bio-gasification, suitable nutrient solutions are needed to stimulate microbial activities toward coal depolymerization and conversion to biogas. Specific for bituminous coal in the Illinois basin, a nutrient recipe that can be used to enhance coal bio- gasification in situ is developed. Based on our previous experience working with Illinois coal and the characteristics of the formation water and the indigenous microbial community, four parameters, Fe powder, methanol, ethanol, and trace minerals are optimized through using a Box-Behnken design. The optimal condition predicted by the models is: Fe-powder at 74 mM, methanol at 97.91 mM, ethanol at 100 mM, and trace minerals at 100%. Under these conditions, the predicted methane yield and content is 1,417.35 ft 3 /ton and 80.7%, respectively. These results are then validated by experimental studies. In addition, each added component is evaluated in terms of its contribution to methane generated. Specifically, the role of coal in the bio-gasification process is studied against two other solid materials. Experimental data from different solid matrices prove that coal is indeed degraded by the studied microbial community.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRSDate not supplied

Aerothermal Analysis of the Dragonfly Titan Entry

The Dragonfly mission will send a rotorcraft lander to the surface of Saturn’s moon Titan as part of the New-Frontiers program. This will be the first spacecraft to land on Titan since the Huygens probe’s descent in January of 2005 and only the second spacecraft to enter Titan’s atmosphere. The Dragonfly entry capsule is significantly larger than the Huygens probe and will experience higher aerothermal environments. This poster will provide an overview of the aerothermodynamic models used in the Dragonfly aeroshell design process. Titan provides a unique entry environment that has several fundamental differences from the environments at the more common entry destinations of Earth and Mars. The lack of atmospheric oxygen significantly reduces heatshield recession (expected to be negligible), which simplifies material response analysis but also reduces the efficiency of the thermal protection system. The atmos-phere is composed of mostly nitrogen and trace amounts of methane. During hypersonic entry this methane dissociates and leads to the formation of the molecule CN, which is known to radiate strongly in the shock layer. Current Dragonfly heating predictions estimate that radiation contributes about 50% of the total heat flux along the forebody and up to 90% on the aftbody, making radiation a key element of aerothermal analysis. To properly address the importance of radiation, a unique radiative-heating correlation was developed that incorporates normal-shock equilibrium chemistry predictions based on work from the Mars 2020 program. These new heating correlations agree well with Computational Fluid Dynamics (CFD) predictions and enable large-scale trajectory analyses by providing rapid heating predictions along candidate trajectories. Dragonfly aerothermal environments are generated with the Data Parallel Line Relaxation (DPLR) finite-volume Navier-Stokes solver alongside the Nonequilibrium Radiative Transport and Spectra (NEQAIR) radiation transport code. DPLR produces shock-aligned flowfields that incorporate the effects of both chemical and thermal nonequilibrium while NEQAIR solves the radiation transport equations along lines-of-sight through the flowfield solutions to predict the radiative heating on the aeroshell. Titan’s atmosphere produces complex thermal and chemical nonequilibria and so classical radiation approximations that extrapolate from a single line-of-sight, such as the Tangent-Slab method, were found to overpredict peak radiative heat fluxes by 10-20%. Therefore, all radiative heating is directly calculated using 3-D transport for both the forebody and aftbody to reduce unnecessary conservatism. A 21-species finite-rate chemistry model, including electrons and ionized species, is used in the CFD to ensure accurate vehicle heating predictions (the level of ionization is particularly for accurate radiative flux predictions). This poster will cover: 1) The models and assumptions used to simulate entry into Titan’s atmosphere, 2) Dragonfly aerothermal engineering correlations alongside statistical heating values, and 3) Driving aerothermal considerations for a Titan entry such as CN formation and radiation-flowfield coupling.

EDL↗
NASA NTRSDate not supplied

Shoulder Postures in EVA Training in Reduced Gravity Analogues

Shoulder Postures in EVA Training in Reduced Gravity Analogues K. Guhl1, L. Vu2, H. Kim3, S. Rajulu4 1KBR Inc., Houston, TX, 2Aegis Aerospace Inc., Houston, TX, 3Leidos Innovations, Houston, TX, 4NASA Johnson Space Center, Houston, TX. During extravehicular activities (EVAs) and EVA training in both the Neutral Buoyancy Laboratory (NBL) and at the Active Response Gravity Offload System (ARGOS), crewmembers perform a variety of hand-intensive tasks with frequent arm/shoulder repositioning while wearing a pressurized spacesuit. As a result, crewmembers may experience ergonomic stressors such as awkward shoulder postures. The ergonomic shoulder risk is also compounded by limited or restricted shoulder mobility of the spacesuit, extreme work positions such as overhead tasks, and tasks with heavy tools and repetitive motions. Prolonged or frequent shoulder elevation and overhead work, in particular, can lead to excessive stresses and musculoskeletal injuries of the shoulder joints. Future EVA missions, specifically lunar surface EVAs, will also be longer in duration and thus increase the exposure to awkward shoulder postures. In this study, we aimed to assess the ergonomic risk of awkward shoulder postures in simulated lunar surface EVAs by quantifying when and how long the arms are raised above the chest level. We assessed video recordings of pilot lunar EVA simulations (3 lunar trials each in the NBL and at ARGOS and 2 microgravity EVA trials in the NBL). These runs consisted of both training and engineering test objectives. Actual demands for shoulder use varied for different EVA types and analogues, but many EVA runs have common tasks and require similar motion components. A video observation and event logging software was used to document the duration and occurrence of the subject’s arms being raised throughout the video recordings of each EVA run. An arm raised instance was classified as the arm being at 90 degrees or above with relation to gravity for lunar EVA training events and with relation to the body for microgravity EVA training events. Such events included EVA hardware maintenance or heavy geology sampling tool operations. Events where the arm load was partially supported by external objects, like climbing a ladder or leaning against a surface were separately identified and excluded. Statistical analysis was performed to summarize the timing, frequency, and durations of the arm raise events and compared across the different EVA tasks and analogue types. Preliminary observations indicated that the total duration and number of arm raised instances were surprisingly smaller for lunar surface EVA training as compared to microgravity EVA training. The observed difference may be attributed to differing task demands and unique environmental characteristics found in lunar EVAs in comparison to microgravity EVAs. A detailed statistical analysis will be performed between lunar surface EVA training events and microgravity EVA training events in the final submission. Overall, this analysis is expected to provide insight into how ergonomic recommendations can be refined for lunar EVA training with pressurized suits. It may also inform task design and influence suit padding design to better protect crewmembers during future lunar training.

Kaitlyn Lea Guhl↗
NASA NTRSDate not supplied

A Machine Learning Approach to Improve Air Traffic Management Initiatives

Collaborating closely with commercial air carriers and related organizations, the Federal Aviation Administration(FAA) regulates air traffic and ensures the safety and efficiency of air operations. Air traffic controllers make strategic decisions, such as delaying, rerouting, or canceling flights, partly based on guidance provided by the FAA’s Air TrafficControl System Command Center (ATCSCC). The guidance includes, among other things, control measures known asTraffic Management Initiatives (TMIs) designed to enhance safety and improve operational efficiency. TMIs play a crucial role in managing the demand and capacity within the U.S. National Airspace System (NAS). Two major TMIs that are routinely used (primarily to mitigate the adverse effects of bad weather) are Ground Delay Programs (GDPs) andGround Stops (GSs). In a GDP, flights destined for airports facing thunderstorm activity experience delays at their origin airports. This proactive approach minimizes the risk of routing aircraft through hazardous weather conditions and also replaces (fuel burning) airborne delays with ground delays. In a GS, a temporary restriction is imposed on the departure or arrival of aircraft at a specific airport or within a designated airspace. Although other TMIs (e.g., miles-in-trail) are also implemented as part of (air) traffic flow management in the NAS, the focus of this work is on GDPs and GSs. Since TMIs, by design, lead to flight delays or cancellations, it is crucial to put in place the right set of parameters(e.g., scope and duration of the GDP). For example, when the end time of a GDP extends beyond what is necessary, it imposes unnecessary delays on departing flights. This situation could occur as a result of inaccurate prediction of the(required) duration of the GDP based on the weather forecast. On the other hand, if a GDP ends prematurely before the underlying capacity constraints are resolved at the destination airport, it may result in airborne holding. The delicate balance lies in matching the termination of the GDP precisely with the resolution of capacity constraints, avoiding both the imposition of unnecessary ground delays and the need for airborne holding due to premature program termination.Failing to specify the right parameters for TMIs also leads to flight delays, creating a significant obstacle in managing the increasing traffic volumes causing increased work load for the controllers. To address this issue, we propose the integration of Machine Learning (ML) models in the traffic flow management(TFM) pipeline. In current operations, decisions are made by human experts based on extensive training, historical patterns, available traffic and weather data. Since we have an abundance of data from past events that tell us the likely impact of various TMIs, by ingesting historical data, properly trained ML models can offer valuable insights and aid human decision-making. With the FAA increasingly exploring advanced analytics, ML emerges as a focal point for enhancing TFM within the National Airspace System (NAS). As a first step, this study aims to provide traffic controllers with decision-making support for the issuance and adjustment of TMIs. Data analytics and machine learning have been previously employed to address some of the challenges associated with TMIs. Numerous studies have concentrated on various facets of TMI issuance, exploring factors influencing TMI parameters, including arrival rate, airport capacity, and delay prediction. For example, using weather forecasts, several statistical methods were used to produce probabilistic capacity profiles which in conjunction with deterministic models provided insights into the GDP planning process [1–4]. The downside of using deterministic models is that they rely on fixed inputs and predetermined rules, which lack the ability to account for the inherent uncertainty and variability present in real-world scenarios. In a separate series of studies, researchers aimed to predict the occurrences of GDPs and GSs. The majority of these studies utilized various supervised learning methods, including Decision Trees, Naive Bayes, Support VectorMachines, and Random Forests to analyze the influence of weather conditions and arrival demand on TMI incidents[5–8]. However, these studies primarily focused on predicting the incidence of TMIs without explicitly addressing the scope of TMIs, including their duration and their geographical coverage. Furthermore, the emphasis of these studies was largely on GDPs, given their higher frequency and longer duration when compared to GSs. A limited number of studies focused on predicting the parameters of TMIs, specifically addressing their duration and extent. In one such study focusing on optimizing the TMI parameters at San Francisco International Airport (SFO),the authors utilized a probabilistic forecast of fog [9]. They simulated various capacity scenarios based on the (fog)burn-off forecasts, selecting GDP parameters that minimized airborne and overall ground delays. However, this approach exclusively emphasizes stratus (fog) burn-off as the primary determinant of GDP and GS, neglecting other influential factors like severe weather events, runway closures, lower capacity than traffic demand, and other important variables. Given the complexity of predicting the TMI and determining its scope, we seek a more holistic approach. We aim to consider all significant factors that could impact TMIs and their parameters. What sets this research apart is the fusion of all data sources relevant to the issuance and adjustment of TMIs and it represents the first comprehensive attempt to optimize TMIs in this manner. Since this comprehensive solution involves various aspects, we break down the problem into smaller components and input all parameters into a unified model called the “TMI Adjuster”. Figure 1 shows the overall framework and the list of datasets used in each model. The objective of the TMI Adjuster module is to deliver reliable, consistent and expedited recommendations for the progression, adjustment, and termination of TMIs. The ML solution entails developing a pipeline capable of predicting the necessity of a TMI (e.g., GS or GDP) along with its various parameters. For example, in the case of a GS, this includes the scope of the GS either in terms of distance from the destination airport or based on pre-defined airspace sectors. Here, scope refers to those regions and departing airports that are subject to the GS. In this paper, we concentrate on the issuance of GSs in the three major airports in the New York area — LaGuardia(LGA), John F. Kennedy International (JFK), and Newark Liberty International (EWR). We fuse traffic, weather and other relevant aviation data from years 2017 to 2019 to train and validate the ML models. In particular, we use the following datasets: •Terminal Aerodrome Forecast (TAF): meteorological forecasts specific to each airport, issued four times a day, covering predefined time periods. •TMI data: includes all GSs and GDPs along with their respective parameters. •Aviation System Performance Metrics (ASPM): includes traffic related data such as aircraft delays, arrival, and departure rates. •Notices to Airmen (NOTAMs): utilized to extract runway closure data and manage interdependencies between terminals in close proximity. •Flight cancellation data •Airspace Flow Programs (AFP): includes information on flight airborne holdings caused by TMIs. The data preprocessing entails transforming ASPM, TMI, AFP, NOTAMs, and weather data into an hourly format and consolidating all datasets by merging them based on date and time as the primary key. The TMI Adjuster framework comprises two parallel models: one dedicated to GS and a second model focused on GDP. As previously mentioned, our specific focus is on the GS model as a multi-classification problem. In this framework, each data point of the GS model input summarizes ten hours of data. Specifically, the data loader for the GS model generates the input and output of the model as follows: at a given time step, the input includes the actual traffic, weather, and TMI data from the two-hour window before the time step, alongside the weather forecast and scheduled traffic for the next 8 hours starting from the time step. Based on this information, the output of the GS model for each time interval consists of three dimensions. The first dimension represents a binary decision on whether there should be a GS in place for the next hour or not. The second dimension is related to the scope of the GS in the United States, and the third dimension is related to the scope of the GS in Canada (i.e., to determine if the GS impacts airports in Canada).One of the challenges with TMI modeling is the sparsity of TMI events, particularly regarding its scope. To address this challenge in the scope of the GS model output, we implement grouping. The GS scope for the US region is defined based on a list of centers that should be included when the GS is in place. With 20 centers in the US, we utilized historical data to group them into 4 categories. In particular, we summarized our historical data in a graph format where nodes represent centers, and link weights are defined based on the co-occurrence of centers in the scope parameter ofTMIs. By identified strongly connected components in this graph, we were able to partition the centers into four groups. We consider two model structures for the GS Model. Firstly, a hierarchical classification model [10], where the human decision-making for a GS is of hierarchical nature. The decision-maker first decides whether there is a need fora GS, and if the answer is yes, determines the scope. A hierarchical classification model organizes the problem into a class hierarchy, typically a tree or a Directed Acyclic Graph (DAG) structure, and considers the dependency of the decision in the previous step to the next component [10]. Here, we employ the local classifier per level approach, which involves training one multi-class classifier for each level of the class hierarchy. The second structure is the independent structure. In this setting, as the name suggests, we do not consider the dependency of the decisions in the different dimensions of the output of the model. Instead, for each dimension, we train a multi-class classifier independently. Table 1 summarizes GS model statistics for training, validation and testing. The table documents the effect of limiting data to the time steps when there was actually a TMI in place or when a TMI had just terminated. This resulted in a more balanced distribution of the GS class(GS positive class)versus “No GS”(GS negative class), which might help the training process. While JFK and LGA follow very similar distributions, with 40% and 42% GS positive class respectively, EWR has proportionally fewer GS incidents at 28%. Our subsequent phase involves evaluating the performance of both hierarchical structure and independent structure using different state-of-the-art multi-class classifier models such as Random Forest, Decision Trees, K-nearest Neighbors, and Logistic Regression and forecast the duration and scope of the GSs.

Farzan Masrour Shalmani↗
DOE OSTI2018

The Relative Importance of Assumed Infrasound Source Terms and Effects of Atmospheric Models on the Linear Inversion of Infrasound Time Series at the Source Physics Experiment

We invert far-field infrasound data for the equivalent seismoacoustic time-domain moment tensor to assess the effects of variable atmospheric models and source phenomena. The infrasound data were produced by a series of underground chemical explosions that were conducted during the Source Physics Experiment (SPE), which was originally designed to study seismoacoustic signal phenomena. The first goal is to investigate the sensitivity of the inversion to the variability of the estimated atmospheric model. The second goal is to determine the relative contribution of two presumed source mechanisms to the observed infrasonic wavefield. Rather than using actual atmospheric observations to estimate the necessary atmospheric Green’s functions, we build a series of atmospheric models that rely on publicly available, regional-scale atmospheric observations. The atmospheric observations are summarized and interpolated onto a 3D grid to produce a model of sound speed at the time of the experiment. For each of four SPE acoustic datasets that we invert, we produced a suite of three atmospheric models for each chemical explosion event, based on 10 yrs of meteorological data: an average model, which averages the atmospheric conditions for 10 yrs prior to each SPE event, as well as two extrema models. To parameterize the inversion, we assume that the source of infrasonic energy results from the linear combination of explosion-induced surface spall and linear seismic-to-elastic mode conversion at the Earth’s free surface. We find that the inversion yields relatively repeatable results for the estimated spall source. Conversely, the estimated isotropic explosion source is highly variable. This suggests that 1) the majority of the observed acoustic energy is produced by the spall and/or 2) our modeling of the elastic energy, and the subsequent conversion to acoustic energy, is too simplistic.

58 GEOSCIENCES↗
NASA NTRS1998

Proposed Mars Surveyor Landing Sites in Northern Meridiani Sinus, Southern Elysium Planitia, and Argyre Planitia

Our objective is to propose two landing sites that the Mars Surveyor 2001 Lander and Athena Rover could go to on Mars that should meet the safety requirements of the spacecraft landing system and optimize surface operations (chiefly driven by power and communications requirements). An additional site within Argyre Planitia, initially proposed by Parker to the Mars Surveyor Landing Site program, is also proposed for potential consideration for post-2001 missions to Mars, as it is well outside the current latitude limits for the Athena Rover. All three sites are designed to be situated as close to a diversity of geologic units within a few kilometers of the landing site so that diversity can be placed in a geologic context. This objective is very different from the Mars Pathfinder requirement to land at a site with a maximum chance for containing a diversity of rocks within a few tens of meters of the lander. That requirement was driven by the Sojourner mobility limit of a few tens of meters. It can be argued that the Athena project, with its much larger mobility capability, might actually want to avoid such a site, because placing collected samples in geologic context would be difficult. While it has been argued, both before and after the Mars Pathfinder landing, that the provenance for local blocks may be determined by orbiter spectra, primarily from the MGS TES instrument, our ability to do so has yet to be demonstrated. Indeed, several months after conclusion of the Pathfinder mission, we have yet to reach a consensus on the composition of local materials. Our primary data set for selecting a landing site within the latitude and elevation constraints of the 2001 mission is the Viking Orbiter image archive. The site must be selected to place the landing ellipse so as to avoid obvious hazards, such as steep slopes, large or numerous craters, or abundant large knobs. For this purpose, we chose a resolution limit of better than 50 m/pixel. This necessarily excludes from the present study images from current and future orbiter spacecraft, until such data does become readily available. Within each proposed region, it may be possible to identify additional sites once these data become available. Second, the fine-component thermal inertia data, should be greater than about 5 or 6 cgs Units (10(exp -3) cal/sq cm s(exp -0.5)/K). Low thermal inertias imply dusty environments, which could pose a mobility hazard. Similarly, the albedo of the site should not be particularly high, which would also suggest dusty surfaces. Low albedos are preferred, as they often coincide with low Viking red:violet ratios and indicate less dusty surfaces. Next, the Modeled Block Abundance should also not be too high or too low. Based on the Viking Lander and Mars Pathfinder experiences, percentages of blocks should be on the order of 5-25%. Too many blocks could pose a hazard to the landing and mobility. Too few blocks could also indicate a dusty surface. Primary Landing Site: Northern Meridiani Sinus (Proposed by T. J. Parker and K., S. Edgett) Vital Statistics: (1) Latitude, Longitude: 0-3 N, 350-2 W. *Elevation (Viking): about0.5-1.5 Ian. (2) Viking Orbiter Image coverage: Excellent coverage by 15 - 25 m/pixel images (orbits 709A and 410B). Possible stereo coverage in region where two orbits overlap (probably small parallax angle, as these orbits are not listed in NASA Contractor Report 3501) (3) Albedo: about .18 -.26 (4) Block Abundance: 5-26% (5)Fine-Component Thermal Inertia: 5-9 cgs units This region consists of bright deposits similar to those described by Edgett et al, that also lie within a prominent dark albedo region. These deposits are flat-lying, to such a degree that they ramp against topography rather than draping over it. This led Edgett and Parker to suggest that they may be subaqueous sediments, possibly lacustrine or marine evaporites, laid down sometime from the late Noachian to middle Hesperian (age determination pending crater counts). A contact between this material and elevated, dissected highlands to the south was identified , and is described by Edgett et al. Our desire in proposing this landing site is to sample the edge of this deposit where it has been exposed through etching, presumably eolian deflation (the deposit, though in the highlands, is itself only lightly to moderately cratered). This should enable access to in situ stratigraphy. The actual landing site will be selected where slopes are not expected to be steep, such that the rover itself should be able to traverse them and sample layered materials on the way, either up or down the slope. Perhaps due to uncertainties at this time as to the friability or meter-scale roughness of the deposit, it might make sense to place the landing ellipse on the exhumed highland surface adjacent to the deflated margin of the deposit and plan on driving to the deposit rather than landing on it and driving downslope. This should also enable imaging the margin for evidence of layering should it prove too difficult to climb. A target ellipse on the highland surface should also allow Athena access to ancient Noachian highland materials, particularly if placed near crater ejecta or an inlier of knobby material. Secondary Landing Site: Southern Elysium Planitia (Proposed by T. J. Parker) Vital Statistics: (1) Latitude, Longitude: 1.5-3.5 S, 195-198 W. (2) Elevation (Viking): -1.0 km. (3) Viking Orbiter Image coverage: Excellent coverage by 15 - 25 m/pixel images (orbit 725). Possible stereo coverage between images from beginning and end of orbit that overlap (probably small parallax angle) (4) Albedo: about .27-.28 (5) Block Abundance: 4-7% (6) Fine-Component Thermal Inertia: about 3 cgs units This region consists of eroded knobby material, probably of Noachian age, though much of the crater population has been destroyed, that is onlapped at a sharp contact by an extensive plains unit in southern Elysium Planitia that is Amazonian in age. The plains materials have been attributed to unusually low-viscosity flood lavas from fissures south of the Elysium volcanic rise, or to lacustrine materials associated with a large, Amazonian lake at the source of Marte Vallis. Parker and Schenk presented evidence in support of the latter interpretation, though they attributed the putative shore morphology to an embayment of a northern plains ocean into the southern Elysium region. Detailed examination of the margin of the deposit, showing erosion, not simply burial, of small crater rims and fluidized ejecta blankets, also points to lacustrine or marine sedimentation rather than volcanic plains burial. The plains surface exhibits a "crusty" appearance that many researchers have attributed to pressure ridges in lava flows. In a lacustrine context, they also resemble pressure ridges in desiccated evaporite deposits and salt-rimmed pools (now dry) similar in scale and morphology to spectacular, hundred meter-scale pool rims in alkaline Lake Natron, East African Rift. The eroded highland margin surface adjacent to these plains appears to be fairly smooth, even at 15 m/pixel. Isolated knob inliers are scattered from a few kilometers to several tens of "kilometers apart. Heights of the knobs have not been measured yet but, based on experience with similar features in the Pathfinder landing ellipse, are probably typically on the order of several tens of meters high and smaller, though some of the largest knobs in the region are probably up to a few hundred meters high. Two craters larger than a kilometer in diameter, with fluidized deposits, lie nearby the proposed landing site. Very high-resolution images from MOC should help to determine whether a landing site navigable by the Athena rover could be placed in this region. The space between knobs and craters is large enough to enable placement of a target landing ellipse between them but still provide access to one or more of them and to the margin of the Elysium plains material. Post-2001 Mars Surveyor Landing Site: Argyre Planitia (Proposed by T. J. Parker) Vital Statistics: (1) Latitude, Longitude: 55-56 S, 41-43 W. (2) Elevation (Viking): 1.0 km. (3) Viking Orbiter Image coverage: Excellent coverage by 40 m/pixel images (orbits 567B, 568B, and 569B). Excellent stereo coverage with large parallax angles over the entire landing site region, and much of central and southern Argyre. (4) Albedo: about .23-.24 (5) Block Abundance: No data (6) Fine-Component Thermal Inertia: No data The floors of both the Argyre and Hellas basins contain etched layered materials that are probably thick accumulations of channel or lacustrine sediments. The deposits in Hellas are much more eroded than those in Argyre, and Hellas lacks a channel outlet. Argyre is unique in that Uzboi Vallis flowed out of the basin, requiring overflow of a standing body of water within Argyre. This makes it the largest impact basin on Mars with channels both draining into it and flowing out from it. Hellas' channels may be catastrophic flood channels, whereas Argyre was fed by modest-scale valley networks, though the outlet at Uzboi Vallis was a catastrophic flood Highland craters and basins of this kind should be high-priority landing targets for missions intended to focus on the search for either prebiotic organic materials or even simple fossil microorganisms. Basins with internally-draining valley networks should be preferred over flood channels, as they could have provided the long-term influx of water favorable to the origin of life. (Catastrophic floods are not conducive to fossil preservation, due to their very short durations and high transportation energies). They also afford an opportunity to study the evolution of the planet's climate and volatiles during the period of time between the late Noachian and early Hesperian, when a drastic change from a proposed early warm, wet climate to one more closely resembling the modern environment is thought to have occurred. Large basin

Parker, T. J.↗
NASA NTRS2014

Updates of the KArLE Experiment: New Libs Calibration Under High Vacuum for the Quantification of Potassium in Basalt for In Situ Geochronology

In planetary exploration, in situ absolute geochronology is one of the main important measurements that needs to be accomplished. Until now, on Mars, the age of the surface is only determined by crater density counting, which gives relative ages. These ages can have a lot of uncertainty as they depend on many parameters. More than that, the curves must be ties to absolute ages. Thus far, only the lost lander Beagle 2 was designed to conduct absolute geochronology measurements, though some recent attempts using MSL Curiosity show that this investigation is feasible and should be strongly encouraged for future flight. Experimental: The Potassium (K)-Argon Laser Experiment (KArLE) is being developed at MSFC through the NASA Planetary Instrument Definition and Development Program (PIDDP). The goal of this experiment is to provide in situ geochronology based on the K-Ar method. A laser ablates a rock under high vacuum, creating a plasma which is sensed by an optical spectrometer to do Laser Induced Breakdown Spectroscopy (LIBS). The ablated material frees gases, including radiogenic 40Ar,which is measured by a mass spectrometer (MS). As the potassium is a content and the 40Ar is a quantity, the ablated mass needed in order to relate them. The mass is given by the product of the ablated volume by the density of this material. So we determine the mineralogy of the ablated material with the LIBS spectra and images and calculate its density. The volume of the pit is measured by using microscopy. LIBS measurement of K under high vacuum: Three independant projects [1, 2, 3] including KArLE, are developing geochronological instruments based on this LA-LIBS-MS method. Despite several differences in their setup, all of them have validated the methods with analyses and ages. However, they all described difficulties with the LIBS measurements of K [3,4]. At ambient pressure, the quantification of K by LIBS on geological materials can be accurate [5]. However the protocol of the LA-LIBS-MS experiment required hundreds of shots under high vacuum in order to free enough 40Ar* to be measured by the QMS. This long duration of ablation may induces significant changes in the LIBS spectra. The pressure may increases by orders of magnitudewithin the chamber and the laser pit geometry can change the effectiveness of ablation and intensity of plasma light received. These effects introduce variation between the first and last spectra and so the quantification of K is more complex. The ablation of one crater can give, depending on the protocol of acquisition, from tens to hundreds of spectra. Protocol and results: We are in the process of further characterizing the variation introduced into LIBS spectra by the use of hundreds of laser shots, and definining a protocol that can be used to ensure accuracy and reporoducibility in the results.We are using natural rock powder standards fused in a furnace, as well as mars analog samples with known K content. We will show the result of the calibration and some new statistical approaches in order to apprehend the effects of the long time ablation on rocks under high vacuum.

Devismes, D.↗