Search NASA⌕ Search

SEARCH · Search NASA

Results for “parameter continuation method”

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 469 records · Page 26

Generalized Software Architecture Applied to the Continuous Lunar Water Separation Process and the Lunar Greenhouse Amplifier

This innovation provides the user with autonomous on-screen monitoring, embedded computations, and tabulated output for two new processes. The software was originally written for the Continuous Lunar Water Separation Process (CLWSP), but was found to be general enough to be applicable to the Lunar Greenhouse Amplifier (LGA) as well, with minor alterations. The resultant program should have general applicability to many laboratory processes (see figure). The objective for these programs was to create a software application that would provide both autonomous monitoring and data storage, along with manual manipulation. The software also allows operators the ability to input experimental changes and comments in real time without modifying the code itself. Common process elements, such as thermocouples, pressure transducers, and relative humidity sensors, are easily incorporated into the program in various configurations, along with specialized devices such as photodiode sensors. The goal of the CLWSP research project is to design, build, and test a new method to continuously separate, capture, and quantify water from a gas stream. The application is any In-Situ Resource Utilization (ISRU) process that desires to extract or produce water from lunar or planetary regolith. The present work is aimed at circumventing current problems and ultimately producing a system capable of continuous operation at moderate temperatures that can be scaled over a large capacity range depending on the ISRU process. The goal of the LGA research project is to design, build, and test a new type of greenhouse that could be used on the moon or Mars. The LGA uses super greenhouse gases (SGGs) to absorb long-wavelength radiation, thus creating a highly efficient greenhouse at a future lunar or Mars outpost. Silica-based glass, although highly efficient at trapping heat, is heavy, fragile, and not suitable for space greenhouse applications. Plastics are much lighter and resilient, but are not efficient for absorbing longwavelength infrared radiation and therefore will lose more heat to the environment compared to glass. The LGA unit uses a transparent polymer antechamber that surrounds part of the greenhouse and encases the SGGs, thereby minimizing infrared losses through the plastic windows. With ambient temperatures at the lunar poles at 50 C, the LGA should provide a substantial enhancement to currently conceived lunar greenhouses. Positive results obtained from this project could lead to a future large-scale system capable of running autonomously on the Moon, Mars, and beyond. The software for both applications needs to run the entire units and all subprocesses; however, throughout testing, many variables and parameters need to be changed as more is learned about the system operation. The software provides the versatility to permit the software operation to change as the user requirements evolve.

Perusich, Stephen↗

Ice-Ridge Pile Up and the Genesis of Martian "Shorelines"

Unique geomorphologic features such as basin terraces exhibiting topographic continuity have been found within several Martian craters as shown in Viking, MOC, and THEMIS images. These features, showing similarity to terrestrial shorelines, have been mapped and cataloged with significant effort [1]. Currently, open wave action on the surface of paleolakes has been hypothesized as the geomorphologic agent responsible for the generation of these features [2]. As consequence, feature interpretations, including shorelines, wave-cut benches, and bars are, befittingly, lacustrine. Because such interpretations and their formation mechanisms have profound implications for the climate and potential biological history of Mars, confidence is crucial. The insight acquired through linked quantitative modeling of geomorphologic agents and processes is key to accurately interpreting these features. In this vein, recent studies [3,4] involving the water wave energy in theoretical open water basins on Mars show minimal erosional effects due to water waves under Martian conditions. Consequently, sub-glacial lake flattens the surface, produces a local velocity increase over the lake, and creates a deviation of the ice flow from the main flow direction [11]. These consequences of ice flow are observed at Lake Vostok, Antarctica an excellent Martian analogue [11]. Martian observations include reticulate terrain exhibiting sharp inter-connected ridges speculated to reflect the deposition and reworking of ice blocks at the periphery of ice-covered lakes throughout Hellas [12]. Our model determines to what extent ice, a terrestrial geomorphologic agent, can alter the Martian landscape. Method: We study the evolution of crater ice plugs as the formation mechanism of surface features frequently identified as shorelines. In particular, we perform model integrations involving parameters such as ice slope and purity, atmospheric pressure and temperature, crater shape and composition, and an energy balance between solar flux, geothermal flux, latent heat, and ablation. Our ultimate goal is to understand how an intracrater ice plug could create the observed shoreline features and how these

Barnhart, C. J.↗

Assessment of Bulk Absorber Properties for Multi-Layer Perforates in Porous Honeycomb Liners

CONTINUING progress in materials technology provides potential for improved acoustic liners for attenuating broadband fan noise emissions from aircraft engine nacelles. Conventional liners (local-reacting perforate-over-honeycomb structures) provide significant narrow-band attenuation, but limited attenuation over wide bandwidths. Two approaches for increasing attenuation bandwidth are to (1) replace the honeycomb structure with bulk material, or (2) cascade multiple layers of perforate/honeycomb structures. Usage of the first approach is limited because of mechanical and maintenance reasons, while multi-layer liners are limited to about three layers because of their additional mechanical complexity, depth and weight. The current research concerns a novel approach reported by the University of Cincinnati, in which a single-layer conventional liner is converted into an extended-reaction, broadband absorber by making the honeycomb core structure porous. This modified single-layer liner requires no increase in depth and weight, and minimal increase in mechanical complexity. Langley has initiated research to identify potential benefits of liner structures with porous cell walls. This research has two complementary goals: (1) develop and validate experimental techniques for treating multi-layer perforates (representative of the internal cells of a liner with porous cell walls) as 1-D bulk materials, and (2) develop analytical approaches to validate this bulk material assumption. If successful, the resultant model can then be used to design optimized porous honeycomb liners. The feasibility of treating an N-layer perforate system (N porous plates separated by uniform air gaps) as a one-dimensional bulk absorber is assessed using the Two-Thickness Method (TTM), which is commonly used to educe bulk material intrinsic acoustic parameters. Tests are conducted with discrete tone and random noise sources, over an SPL range sufficient to determine the nonlinearity of the test specimens, for test specimens consisting of 5, 10 and 15% porous plates. Measured impedances for two liner thicknesses (e.g., 12 and 24 layers) are used as input to the TTM to determine the characteristic impedance and propagation constant that characterize these liners as bulk absorbers. These parameters are then used to calculate the predicted impedance of liners with different thicknesses (e.g., 36 layers), and a comparison of predicted and measured impedances for these other thicknesses is used to determine the efficacy of this approach. Finally, an independent method is used to educe the propagation constant for a single representative sample, and excellent comparison between the results for this method and those for the TTM provides increased confidence in the results achieved with the TTM. In general, the results demonstrate these multi-layer perforates can be acceptably treated as bulk absorbers.

Jones, Michael G.↗

Advances in Space Radiation Physics and Transport

The space radiation environment is a complex mixture of particle types and energies originating from sources inside and outside of the galaxy. These environments may be modified by the heliospheric and geomagnetic conditions as well as planetary bodies and vehicle or habitat mass shielding. In low Earth orbit (LEO), the geomagnetic field deflects a portion of the galactic cosmic rays (GCR) and all but the most intense solar particle events (SPE). There are also dynamic belts of trapped electrons and protons with low to medium energy and intense particle count rates. In deep space, the GCR exposure is more severe than in LEO and varies inversely with solar activity. Unpredictable solar storms also present an acute risk to astronauts if adequate shielding is not provided. Near planetary surfaces such as the Earth, moon or Mars, secondary particles are produced when the ambient deep space radiation environment interacts with these surfaces and/or atmospheres. These secondary particles further complicate the local radiation environment and modify the associated health risks. Characterizing the radiation fields in this vast array of scenarios and environments is a challenging task and is currently accomplished with a combination of computational models and dosimetry. The computational tools include models for the ambient space radiation environment, mass shielding geometry, and atomic and nuclear interaction parameters. These models are then coupled to a radiation transport code to describe the radiation field at the location of interest within a vehicle or habitat. Many new advances in these models have been made in the last decade, and the present review article focuses on the progress and contributions made by workers and collaborators at NASA in the same time frame. Although great progress has been made, and models continue to improve, significant gaps remain and are discussed in the context of planned future missions. Of particular interest is the juxtaposition of various review committee findings regarding the accuracy and gaps of combined space radiation environment, physics, and transport models with the progress achieved over the past decade. While current models are now fully capable of characterizing radiation environments in the broad range of forecasted mission scenarios, it should be remembered that uncertainties still remain and need to be addressed.

Space radiation↗

Metallic Environmentally Resistant Coating Rapid Innovation Initiative

Lightweight alloys such as aluminum (Al) and titanium (Ti) are often specified for space systems to minimize mass while maintaining structural integrity [1,2]. Such alloys however, have poor tribological response (high friction and wear), especially in extreme space environments, which becomes worse with the additional presence of lunar regolith. This leads to short lifetimes and premature failures that will ultimately limit long term operations on the lunar surface [2]. This project is addressing this technology gap by developing advanced wear- and radiation-resistant coatings for lightweight parts to extend the life-time and sustainability of both lunar and Martian assets. Ceramics were considered for their high wear resistance, but were rejected because of their low fracture toughness, which would be especially problematic for structural components. The novel and existing coating technologies and deposition methods are being tested in this project. The coating materials are Boron Nitride-Aluminum (BN-Al), Nickle Titanium (NiTi), Aluminum Oxide (AlO3), Ti64 with hBN at 2 and 10 vol percent (Ti-2vol%hBN and Ti-10vol%hBN), and the deposition techniques are high pressure cold-spray (CS) and ambient and vacuum plasma-spray (APS and VPS) [3,4,5]. BN-Al, NiTi, Ti-2vol%hBN, and Ti-10vol%hBN were applied with all three deposition techniques, and AlO3 was applied only using the APS deposition technique. A tungstenite (WS2) film was applied to the NiTi VPS coating. The coating and deposition technique configurations are being tested against several key end-use performance parameters. The parameters include the capabilities of the coatings under wear environments such as regolith simulant, thermal cycling from high (120°C) to cryogenic (-173°C), high vacuum (~10-7 torr), and pre- and post-exposure to ionizing particle radiation. Wear tests include pin on disk, three-body abrasion, and surface erosion by high velocity regolith impacts. The initial down selection is being performed using on virgin and environmental exposure samples using pin on disk and three body abrasion wear tests. From preliminary assessment of this testing, a few configurations were eliminated. Ti-10vol%hBN could not be applied using CS application and NiTi could not be applied using APS or CS application so these configurations did not get tested. Early pin on disk wear tests showed poor wear performance of BN-Al so this coating was eliminated early in the testing process. AlO3 on Al substrate did not survive thermal cycling, but wear testing is continuing for AlO3 on Ti substrate. A more detailed analysis is being con-ducted to further reduce the number of configuration for phase II and III. Phase II testing will include conventionally and additively manufactured substrates with surface erosion testing. In phase III, the coatings will be applied to three mechanism types: channel and slot, ball and socket, and a hinge joint. Each will demonstrate a different type of wear incidence. The mechanisms and their base materials are of direct interest to the end users and infusion points: the Human Landing System (HLS) and the Lunar Surface Innovation Initiative (sustained lunar surface operations). The technology development project is based out of Marshall Space Flight Center, has partnered with Florida International University (Miami, FL) and Plasma Processes (Huntsville, Al), and is supported by a group of NASA mentors from different centers.

dust mitigation↗

Comparison of Interstitial Fluid pH, PCO2, PO2 with Venous Blood Values During Repetitive Handgrip Exercise

We evaluated the use of a small, fiber optic sensor to measure pH, PCO2 and PO2 from forearm muscle interstitial fluid (IF) during handgrip dynamometry. PURPOSE: Compare pH, PCO2 and PO2 values obtained from venous blood with those from the IF of the flexor digitorum superficialis (FDS) during three levels of exercise intensity. METHODS: Six subjects (5M/1F), average age 29+/-5 yrs, participated in the study. A venous catheter was placed in the retrograde direction in the antecubital space and a fiber optic sensor (Paratrend, Diametrics Medical, Inc.) was placed through a 22 G catheter into the FDS muscle under ultrasound guidance. After a 45 min rest period, subjects performed three 5-min bouts of repetitive handgrip exercise (2s contraction/1 s relaxation) at attempted levels of 15%, 30% and 45% of maximal voluntary contraction. The order of the exercise bouts was random with the second and third bouts started after blood lactate had returned to baseline. Venous blood was sampled every minute during exercise and analyzed with an I-Stat CG-4+ cartridge, while IF fiber optic sensor measurements were obtained every 2 s. Change from pre-exercise baseline to end of exercise was computed for pH, PCO2 and PO2. Blood and IF values were compared with a paired t-test. RESULTS: Baseline values for pH, PCO2 and PO2 were 7.37+/-0.02, 46+/-4 mm Hg, and 36+/-6 mm Hg respectively in blood and 7.39+/-0.02, 44+/-6 mm Hg, and 35+/-14 mm Hg in IF. Average changes over all exercise levels are noted in the Table below. For each parameter the exercise-induced change was at least twice as great in IF as in blood. In blood and IF, pH and PCO2 increases were directly related to exercise intensity. Change in venous PO2 was unrelated to exercise intensity, while IF PO2 decreased with increases in exercise intensity. CONCLUSIONS: Measurement of IF pH, PCO2 and PO2 is more sensitive to exercise intensity than measurement of the same parameters in venous blood and provides continuous assessment during and after exercise.

Hagan, Ronald Donald↗

Probabilistic Modeling of a Three-Stage Human Landing System Architecture

Space Policy Directive-1 has led to NASA partnerships with commercial entities on procurement which includes the development of the Human Landing System (HLS) [1]. With the goal of delivering human crew to the lunar surface by 2024, system uncertainties become an important obstacle to the maturation of multiple new, driving technologies and mission concepts of the HLS program. As unmitigated uncertainties have previously led to failed development programs, these risks and their impacts must be understood and handled to ensure program success [2]. Sources of uncertainty include novel engine designs and configurations, increased reliance on cryogenic fluid management(CFM), and refueling technologies—which propagate as high-level performance metrics such as overall propellant mass and engine performance. Also, the occurrence of operational uncertainties—e.g. launch conditions or need to abort during the mission—can cause cascading effects on the rest of the mission that are difficult to definitively quantify, and are outside the scope of control. These concrete examples and other occurrences can be categorized as either epistemic or aleatory uncertainties.Epistemic uncertainty arises due to a lack of knowledge and can be alleviated with design and program maturation. Aleatory uncertainty is due to the inherent randomness of the system and cannot be directly reduced, unlike epistemic uncertainty. Robust design and probabilistic methods can compensate for aleatory effects. A taxonomy of uncertainty is referred to for this work [3]. In this paper, a probabilistic methodology to handle uncertainties has been demonstrated on a three-element HLS concept [1, 4], which allows tracking of current best estimates of the concept and assessment of concept design robustness against uncertainties. A sample case has been completed for this abstract, and an expansion on the methodology will be included in the final paper. This methodology has two key parts: first, the creation of a dynamic architecture model of a three-element HLS concept; and second, its use with surrogate modeling and range estimating techniques to capture and propagate uncertainties. This abstract will cover the basics of the approach used, and further details and justifications will be in the final paper.The mission profile associated with this three-element concept (Fig 1) was modeled as a set of mission events that facilitated mass changes, idles, or spacecraft maneuvers. The mission profile scope starts with each element’s NRHO orbit insertion and aggregation and ends at post-sortie rendezvous with Orion. More detail on the mission profile will be in the final paper. The DYnamic Rocket EQuation Tool (DYREQT), a space systems synthesis and sizing framework used by NASA, was used as the physics framework to model the HLS architecture for applying the probabilistic methodology [5, 6]. Specifically, a parametric representation of the lander, ascent, and transfer elements and the mission profile of each element was established, with vehicle and mission parameters available as inputs to allow for a dynamic model. Each vehicle stage was modeled with high-level performance metrics, using Isp and propellant mass fraction (PMF) to remain parametric. For the probabilistic analysis, uncertainties of interest within the HLS concept were enumerated and represented as parameters within the DYREQT model as inputs for vehicle stages or mission profile events. These parameters were frozen at their nominal values for the purposes of baselining architecture performance and sizing the vehicle appropriately based on reference documentation [1]. Range estimating—a probabilistic method that combines Monte Carlo sampling, focus on critical parameters, and heuristics to assess risk and opportunities—is traditionally used with Mass Equipment Lists (MELs), but has been adapted with operational parameters as well as vehicle parameters in theDYREQT model to capture mission uncertainty alongside vehicle uncertainty [7, 3]. This method was selected due to its application and insight on a system from a bottom-up perspective, independence from historical rules of thumb, and ability to generate sensitivities based on design decisions and uncertainties. As a sample case for the abstract, the boiloff rates of the vehicle elements and the loiter times during the mission (simulating launch time variations and changing window of opportunities) were used with range estimating to provide preliminary results. To perform the range estimation portion of this methodology (depicted in Fig. 3, further details in final paper), the DYREQT model was sampled using a Design of Experiments (DoE) to efficiently explore the architecture design space with respect to the sample set of uncertainty parameters; 5,000 cases via Latin Hypercube Sampling were computed on the DYREQT architecture model. Then, the results were used to create surrogate models, multivariate regressions that can visualize hypercube trends in the design space, of the architecture with respect to the uncertainty parameters. Range estimating was applied to the surrogates instead of the actual models, which saves computational expense due to the bulk of cases needed for the Monte Carlo simulation as part of range estimating. Uncertainty parameters were sampled independently from triangular distributions using the DoE ranges as ‘min’ and ‘max’, and the nominal value as ‘most likely’. Based engineering intuition, some uncertainty parameters are correlated—e.g. if the main propellant has a high boil-off rate, the oxidizer should follow suit as both are related to CFM technology.While a Monte Carlo simulation samples all inputs as independent, the results would show model correlations; thus, it is efficient to sample the inputs as correlated. Using a correlation matrix constructed for the uncertainty parameters, previously independent samples were transformed to perform a Correlated Monte Carlo. A table for the DoE ranges and probability distribution parameters is shown in Table 1, and more details on Correlated Monte Carlo Simulations will be discussed in the final paper. The model’s resulting DoE showed that multivariate polynomial equations fit via least squares method captured its behavior accurately for the sample case. For the Correlated Monte Carlo Simulation, a positive correlation between fuel and oxidizer boiloff rates was used as a demonstration. 10,000 cases were computed with the surrogates and the launched masses for each vehicle element was collated. The results can be displayed in a probability density function (PDF), showing the impact of the uncertainty parameters chosen. Integrating the PDFs will yield a cumulative distribution function (CDF) that shows the cumulative probability of a given value on the x-axis. For the sample case, the elements’ launch mass margin was calculated and represented in as CDFs, as a demonstrated representation of figures of merit for the HLS concept. For the lander and ascent elements, the NRHO mass insertion limit is 16t; the transfer element has a limit of 30t [1]. It can be seen with Figure 2 that this probabilistic methodology can provide insight into mass margin with respect to the uncertainties being modeled. Currently, the results show that the lander (descent) vehicle element has the most restrictive design space; it is the only element to show a 10% probability of negative margin. Further analysis on the Monte Carlo results will show sensitivities for driving constraints and parameters for architecture feasibility, which can lead to establishing potential mission rules.The combination of range estimating with a parametric architecture model for HLS demonstrated the capability of this probabilistic methodology in a sample case. As the HLS development progresses, this methodology has the potential for keeping current best estimates of architecture performance for awarded concepts due to the flexibility in DYREQT’s modeling framework and its parametric nature. Concept maturation and increased epistemic knowledge can be injected into the model probabilistic modeling, and thus continue to track probability of mission success.

Stephanie Y Zhu↗

Local starburst galaxies and their descendants. Statistics from the Sloan Digital Sky Survey

Despite strong interest in the starburst phenomenon in extragalactic astronomy, the concept remains ill-defined. Here we use a strict definition of starburst to examine the statistical properties of starburst galaxies in the local universe. We also seek to establish links between starburst galaxies, post-starburst (hereafter postburst) galaxies, and active galaxies. Data were selected from the Sloan Digital Sky Survey DR7. We applied a novel method of treating dust attenuation and derive star formation rates, ages, and stellar masses assuming a two-component stellar population model. Dynamical masses are calculated from the width of the H-alpha line. These masses agree excellently with the photometric masses. The mass (gas+stars) range is approximately 10( exp 9) - 10(exp 11.5) solar mass. As a selection criterion for starburst galaxies, we use, the birthrate parameter, b = SFR/SFR, requiring that b is greater than 3. For postburst galaxies, we use, the equivalent width of Hdelta in absorption with the criterion EW (sub Hdelta_abs) is greater than 6 A. Results. We find that only 1% of star-forming galaxies are starburst galaxies. They contribute 3-6% to the stellar production and are therefore unimportant for the local star formation activity. The median starburst age is 70 Myr roughly independent of mass, indicating that star formation is mainly regulated by local feedback processes. The b-parameter strongly depends on burst age. Values close to b = 60 are found at ages approximately 10 Myr, while almost no starbursts are found at ages greater than 1 Gyr. The median baryonic burst mass fraction of sub-L galaxies is 5% and decreases slowly towards high masses. The median mass fraction of the recent burst in the postburst sample is 5-10%. A smaller fraction of the postburst galaxies, however, originates in non-bursting galaxies. The age-mass distribution of the postburst progenitors (with mass fractions is greater than 3%) is bimodal with a break at logM(solar mass ) 10.6, above which the ages are doubled. The starburst and postburst luminosity functions (LFs) follow each other closely until M(sub r ) (is) approximately -21, when active galactic nuclei (AGNs) begin to dominate. The postburst LF continues to follow the AGN LF, while starbursts become less significant. This suggests that the number of luminous starbursts is underestimated by about one dex at high luminosities, because of having large amounts of dust and/or being outshone by an AGN. It also indicates that the starburst phase preceded the AGN phase. Finally, we look at the conditions for global gas outflow caused by stellar feedback and find that massive starburst galaxies are susceptible to such outflows.

Statistical distributions↗

Association of American Geographers, Remote Sensing Specialty Group Special Issue of Geocarto International

This special issue continues the precedence of the Association of American Geographers (AAG), Remote Sensing Specialty Group (RSSG) for publishing selected articles in Geocarto International as a by-product from the AAG annual meeting. As editors, we issued earlier this year, a solicitation for papers to be published in a special issue of Geocarto International that were presented in RSSG-sponsored sessions at the 2001 AAG annual meeting held in New York City on February 27-March 3. Although not an absolute requisite for publication, the vast majority of the papers in this special issue were presented at this year's AAG meeting in New York. Other articles in this issue that were not part of a paper or poster session at the 2001 AAG meeting are authored by RSSG members. Under the auspices of the RSSG, this special Geocarto International issue provides even more compelling evidence of the inextricable linkage between remote sensing and geography. The papers in this special issue fall into four general themes: 1) Urban Analysis and Techniques for Urban Analysis; 2) Land Use/Land Cover Analysis; 3) Fire Modeling Assessment; and 4) Techniques. The first four papers herein are concerned with the use of remote sensing for analysis of urban areas, and with use or development of techniques to better characterize urban areas using remote sensing data. As the lead paper in this grouping, Rashed et al., examine the usage of spectral mixture analysis (SMA) for analyzing satellite imagery of urban areas as opposed to more 'standard' methods of classification. Here SMA has been applied to IRS-1C satellite multispectral imagery to extract measures that better describe the 'anatomy' of the greater Cairo, Egypt region. Following this paper, Weng and Lo describe how Landsat TM data have been used to monitor land cover types and to estimate biomass parameters within an urban environment. The research reported in this paper applies an integrated GIS (Geographic Information System) approach for detecting urban growth and assessing its impact on biomass in the Zhujiang Delta, China. The remaining two papers in this first grouping deal with improved techniques for characterizing and analyzing urban areas using remote sensing data. Myint examines the use of texture analysis to better classify urban features. Here wavelet analysis has been employed to assist in deriving a more robust classification of the urban environment from high spatial resolution, multispectral aircraft data. Mesev provides insight on how through the modification of the standard maximum likelihood image analysis technique, population census data can be used enhance the overall robustness of urban image classification through the modification of the standard maximum likelihood image analysis technique.

Allen, Thomas R.↗

Shape Effect Analysis of Aluminum Projectile Impact on Whipple Shields

The informed design with respect to hypervelocity collisions involving micrometeoroid and orbital debris (MMOD) is influential to the success of space missions. For an orbit comparable to that of the International Space Station, velocities for MMOD can range from 1 to 15 km/s, with an average velocity around 10 km/cu s. The high energy released during collisions at these speeds can result in damage to a spacecraft, or worst-case, loss of the spacecraft, thus outlining the importance of methods to predict the likelihood and extent of damage due to an impact. Through experimental testing and numerical simulations, substantial work has been conducted to better understand the effects of hypervelocity impacts (HVI) on spacecraft systems and shields; however, much of the work has been focused on spherical impacting particles. To improve environment models for the analysis of MMOD, a large-scale satellite break-up test was performed at the Arnold Engineering and Development Complex to better understand the varied impactor geometries that could be generated from a large impact. As a part of the post-experiment analysis, an undertaking to characterize the irregular fragments generated is currently being performed by the University of Florida under the management of NASA's Orbital Debris Program Office at Johnson Space Center (JSC). DebriSat was a representative, modern LEO satellite that was catastrophically broken up in a HVI test. The test chamber was lined with a soft-catch system of foam panels that captured the fragments after impact. Initial predictions put the number of fragments larger than 2mm generated from the HVI at roughly 85,000. The number of fragments thus far extracted from the foam panels has exceeded 100,000, with that number continuously increasing. The shapes of the fragments vary dependent upon the material. Carbon-fiber reinforced polymer pieces, for instance, are abundantly found as thin, flat slivers. The characterization of these fragments with respect to their mass, size, and material composition needs to be summarized in a form that can be used in MMOD analysis. The mechanism that brings these fragment traits into MMOD analysis is through ballistic limit equations (BLE) that have been developed largely for a few types of materials1. As a BLE provides the failure threshold for a shield or spacecraft component based on parameters such as the projectile impact velocity and size, and the target's materials, thickness, and configuration, it is used to design protective shields for spacecraft such as Whipple shields (WS) to an acceptable risk level. The majority of experiments and simulations to test shields and validate BLEs have, heretofore, largely used spheres as the impactor, not properly reflecting the irregular shapes of MMOD. This shortfall has motivated a numerical impact analysis study of HVI involving non-spherical geometries to identify key parameters that environment models should provide.

Carrasquilla, Maria J.↗

Real-Time Adaptive Drag Minimization Wind Tunnel Investigation of a Flexible Wing with Variable Camber Continuous Trailing Edge Flap System

This paper reports the results of a recently completed real-time adaptive drag minimization wind tunnel investigation of a highly flexible wing wind tunnel model equipped with the Variable Camber Continuous Trailing Flap (VCCTEF) technology at the University of Washington Aeronautical Laboratory (UWAL). The wind tunnel investigation is funded by NASA SBIR Phase II contract with Scientific Systems Company, Inc. (SSCI) and University of Washington (UW) as a subcontractor. The wind tunnel model is a sub-scale Common Research Model (CRM) wing constructed of foam core and fiberglass skin and is aeroelastically scaled to achieve a wing tip deflection of 10% of the wing semi-span which represents a typical wing tip deflection for a modern transport such as Boeing 787. The jig-shape twist of the CRM wing is optimized using a CART3D aero-structural model to achieve the minimum induced drag for the design cruise lift coefficient of 0.5. The wing is equipped with two chord wise cambered segments for each of the six span wise flap sections for a total of 12 individual flap segments that comprise the VCCTEF system. Each of the 12 flap segments is actively controlled by an electric servo-actuator. The real-time adaptive drag optimization strategy includes an on-board aerodynamic model identification, a model excitation, and a real-time drag optimization. The on-board aerodynamic model is constructed parametrically as a function of the angle of attack and flap positions to model the lift and drag coefficients of the wing. The lift coefficient models include a linear model and a second-order model. The drag coefficient models include a quadratic model and a higher-order up to 6th-order model to accurately model the drag coefficient at high angles of attack. The onboard aerodynamic model identification includes a recursive least-squares (RLS) algorithm and a batch least-squares (BLS) algorithm designed to estimate the model parameters. The model excitation method is designed to sample the input set that comprises the angle of attack and the flap positions. Three model excitation methods are developed: random excitation method, sweep method, and iterative angle-of-attack seeking method. The real-time drag optimization includes a generic algorithm developed by SSCI and several optimization methods developed by NASA which include a second-order gradient Newton-Raphson optimization method, an iterative gradient optimization method, a pseudo-inverse optimization method, an analytical optimization method, and an iterative refinement optimization method. The first wind tunnel test entry took place in September 2017. This test revealed major hardware issues and required further redesign of the flap servo mechanisms. The second test entry took place in April 2018. However, the test was not successful due to the issues with the onboard aerodynamic model identification RLS algorithm which incorrectly identified model parameters. This test also provides an experimental comparison study between the VCCTEF and a variable camber discrete trailing edge flap (VCDTEF) without the elastomer transition mechanisms. The experimental result confirms the benefit of the VCCTEF which produces lower drag by 5% than the VCDTEF. The third and final test entry took place in June 2018 after the issues with the RLS algorithm have been identified and corrected. Additional improvements were implemented. These include the BLS algorithm, the iterative angle-of-attack seeking method, the iterative gradient optimization method, and the pseudo-inverse optimization method. The test objectives were successfully demonstrated as the real-time drag optimization identifies several optimal solutions at off-design lift coefficients. The iterative gradient optimization method is found to achieve up to 4.7% drag reduction for the off-design lift coefficient of 0.7. The pseudo-inverse optimization method which does not require the drag coefficient model is found to be quite effective in reducing drag. Up to 9.4% drag reduction for the off-design lift coefficient of 0.7 is achieved with the pseudo-inverse optimization method. The wind tunnel investigation demonstrates the potential of real-time drag optimization technology. Several new capabilities are developed that could enable future adaptive wing technologies for flexible wings equipped with drag control devices such as the VCCTEF.

Drag Optimization↗

Modeling and simulation of a turbulent far wake

Work continued on two projects which had been started during previous years. Both projects involve calculations of the subsonic, turbulent far wake of a two-dimensional object at a Reynolds number of 1000 (based on wake momentum thickness). This flow was used as a test case for direct comparison of various turbulence models and a direct numerical simulations (DNS) of this flow were undertaken. In the turbulence model comparison studies, for any particular model tested, a unique self-similar solution was obtained far enough downstream, regardless of inlet conditions. Furthermore, different turbulence models led to different far-wake equilibrium solutions. No turbulence model could correctly predict all features of the turbulent far wake. For example, the spreading rate and turbulent shear stresses were underpredicted by all the standard models (both two-equation and full Reynolds stress models). In cases where a more correct spreading rate was achieved, it was at the expense of the turbulent kinetic energy, which was overpredicted. In general, the Algebraic Dissipation Rate Model of Gatski and Speziale, 1992, when added to any of the standard models, improved the results dramatically. Also, full Reynolds stress closure models did a much better job at predicting the shapes of both the mean and turbulence profiles, but the spreading rate was not significantly improved over that predicted by the simpler two-equation models. There are two main conclusions from these studies: First, in a comparison such as this, it is not enough to compare just one parameter, like the spreading rate. A good prediction for one parameter does not necessarily imply good predictions for all parameters in a flow. Second, since no turbulence model could correctly predict the turbulent far wake, much of the important physics of turbulent free shear flows is apparently lost by the assumptions inherent in today's methods of turbulence modeling; turbulence models must be improved. Direct simulations of this flow were begun last year in order to provide a data base through which some of the deficiencies of the existing turbulence models could be identified. Quantities such as the pressure-strain correlation, turbulent diffusion, and the dissipation rate tensor can be easily calculated from the DNS results, whereas these quantities are nearly impossible to measure experimentally. Improvements to existing turbulence models (and development of new models) require knowledge about flow quantities such as these. During this summer, diagnostics codes were written which will calculate the parameters mentioned above, along with other single-point and multi-point statistics. The DNS calculations are still in progress at the time of this writing. When these calculations are complete, the diagnostics codes will be applied so that the results can aid turbulence modelers. In addition, the results will show whether or not there exists a universal equilibrium turbulent far wake, independent of initial conditions.

Cimbala, John M.↗

Potential Use of a Bayesian Network for Discriminating Flash Type from Future GOES-R Geostationary Lightning Mapper (GLM) data

Continuous monitoring of the ratio of cloud flashes to ground flashes may provide a better understanding of thunderstorm dynamics, intensification, and evolution, and it may be useful in severe weather warning. The National Lighting Detection Network TM (NLDN) senses ground flashes with exceptional detection efficiency and accuracy over most of the continental United States. A proposed Geostationary Lightning Mapper (GLM) aboard the Geostationary Operational Environmental Satellite (GOES-R) will look at the western hemisphere, and among the lightning data products to be made available will be the fundamental optical flash parameters for both cloud and ground flashes: radiance, area, duration, number of optical groups, and number of optical events. Previous studies have demonstrated that the optical flash parameter statistics of ground and cloud lightning, which are observable from space, are significantly different. This study investigates a Bayesian network methodology for discriminating lightning flash type (ground or cloud) using the lightning optical data and ancillary GOES-R data. A Directed Acyclic Graph (DAG) is set up with lightning as a "root" and data observed by GLM as the "leaves." This allows for a direct calculation of the joint probability distribution function for the lighting type and radiance, area, etc. Initially, the conditional probabilities that will be required can be estimated from the Lightning Imaging Sensor (LIS) and the Optical Transient Detector (OTD) together with NLDN data. Directly manipulating the joint distribution will yield the conditional probability that a lightning flash is a ground flash given the evidence, which consists of the observed lightning optical data [and possibly cloud data retrieved from the GOES-R Advanced Baseline Imager (ABI) in a more mature Bayesian network configuration]. Later, actual GLM and NLDN data can be used to refine the estimates of the conditional probabilities used in the model; i.e., the Bayesian network is a learning network. Methods for efficient calculation of the conditional probabilities (e.g., an algorithm using junction trees), finding data conflicts, goodness of fit, and dealing with missing data will also be addressed.

Solakiewiz, Richard↗

Intracranial Effects of Artificial Gravity: A 3T MRI Study

INTRODUCTION Spaceflight associated neuro-ocular syndrome (SANS) is characterized by the development of optic disc edema, posterior globe flattening, choroidal/retinal folds and hyperopic refractive errors1. SANS is hypothesized to be a result of headward fluid shifts that invariably occurs in the microgravity environment. As a countermeasure, artificial gravity (AG) through centrifugation has been proposed to reduce this headward fluid shift, however there is no current proof of benefit. The goal of this study was to determine if the application of AG can prevent or reduce known changes in brain volumetry, internal carotid artery (ICA) stroke volume and cerebral spinal fluid (CSF) flow velocity that occurs during simulated chronic headward fluid shift using head down tilt bed rest (HDTBR) methodology2 as an indicator of countermeasure efficacy. METHODS Healthy volunteers were recruited for an IRB approved HDTBR study performed at the German Aerospace Center in Cologne, Germany. Strict six-degree HDTBR was used as a spaceflight analog to induce a continuous headward fluid shift. HDTBR was carried out for 60 days for all subjects. Short-arm centrifugation was utilized to generate AG equating to ~0.3g of acceleration at the level of the eye. The subjects were divided equally into three groups: NoAG (control; n=8), daily intermittent AG (6 x 5 min iAG; n=8), and daily continuous 30 min (cAG; n=8). All studies were performed on a single dedicated 3T MRI Scanner. Pulse-gated MRI phase-contrast flow imaging was used to quantify ICA stroke volume and peak-to-peak CSF flow velocity in the mid cerebral aqueduct. 3D-SPGR was acquired for volumetric segmentation of the brain and CSF spaces. MRI acquisitions were obtained at baseline (BDC), 14 days into HDTBR (HDTBR14), 52 days into HDTBR (HDTBR52) and 3-5 days after HDTBR (recovery, R+3/5).The data were analyzed by the mixed model, which included intervention and time (BDC, HDTBR 14, HDTBR 52, R+3/5) as the fixed effects and included subject as the random effect.RESULTS24 healthy subject volunteers (16 men, 8 women, mean age = 33 years ± 9 [standard deviation] and mean BMI = 24.3 kg/m2 ± 2.0) successfully completed all phases of the study. Strict six-degree HDTBR was characterized by progressive and statistically significant (p<.01) increases in mean combined brain and CSF volumes and mean aqueductal CSF peak-to-peak flow velocity, as well as statistically significant (p<.01) progressive decrease in mean ICA stroke volume from baseline to 52 days post intervention (Figs. 1-3). Compared to baseline, only combined brain and CSF volumes did not return to baseline values in the recovery period (p=NS). Neither iAG nor cAG exerted any significant effects on the measured MRI brain parameters as compared to HDTBR alone (p=NS). CONCLUSION Our results indicate that HDTBR at 6-degrees was effective in producing alterations in ICA stroke volume, aqueductal CSF flow velocity, and combined brain and CSF volumetric change that is associated with chronic headward fluid shift. Short duration, 30-min daily exposure to either iAG or cAG appears to be insufficient in preventing or reducing the effects of chronic HDTBR and thus may not be a suitable countermeasure as currently deployed. AG protocol modifications, including increased duration and magnitude of exposure, should be considered for future research.

L A Kramer↗

Optimal Estimation-Based Algorithm to Retrieve Aerosol Optical Properties for GEMS Measurements Over Asia

The Geostationary Environment Monitoring Spectrometer (GEMS) is scheduled to be in orbit in 2019 onboard the GEO-KOMPSAT 2B satellite and will continuously monitor air quality over Asia. The GEMS will make measurements in the UV spectrum (300-500 nm) with 0.6 nm resolution. In this study, an algorithm is developed to retrieve aerosol optical properties from UV-visible measurements for the future satellite instrument and is tested using 3 years of existing OMI L1B data. This algorithm provides aerosol optical depth (AOD), single scattering albedo (SSA) and aerosol layer height (ALH) using an optimized estimation method. The retrieved AOD shows good correlation with Aerosol Robotic Network (AERONET) AOD with correlation coefficients of 0.83, 0.73 and 0.80 for heavy-absorbing fine (HAF) particles, dust and non-absorbing (NA) particles, respectively. However, regression tests indicate underestimation and overestimation of HAF and NA AOD, respectively. In comparison with AOD from the OMI/Aura Near-UV Aerosol Optical Depth and Single Scattering Albedo 1-orbit L2 Swath 13 km x 24 km V003 (OMAERUV) algorithm, the retrieved AOD has a correlation coefficient of 0.86 and linear regression equation, AOD(sub GEMS) = 1.18AOD(sub OMAERUV) + 0.09. An uncertainty test based on a reference method, which estimates retrieval error by applying the algorithm to simulated radiance data, revealed that assumptions in the spectral dependency of aerosol absorptivity in the UV cause significant errors in aerosol property retrieval, particularly the SSA retrieval. Consequently, retrieved SSAs did not show good correlation with AERONET values. The ALH results were qualitatively compared with the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) products and were found to be well correlated for highly absorbing aerosols. The difference between the attenuated-backscatter-weighted height from CALIOP and retrieved ALH were mostly closed to zero when the retrieved AOD is higher than 0.8 and SSA is lower than 0.93. Although retrieval accuracy was not significantly improved, the simultaneous consistent retrieval of AOD, SSA and ALH alone demonstrates the value of this stand-alone algorithm, given their nature for error using other methods. The use of these properties as input parameters for the air mass factor calculation is expected to improve the retrieval of other trace gases over Asia.

Geostationary Environment Monitoring Spectrometer ↗

An Experimental Study of Upward Burning Over Long Solid Fuels: Facility Development and Comparison

As NASA's mission evolves, new spacecraft and habitat environments necessitate expanded study of materials flammability. Most of the upward burning tests to date, including the NASA standard material screening method NASA-STD-6001, have been conducted in small chambers where the flame often terminates before a steady state flame is established. In real environments, the same limitations may not be present. The use of long fuel samples would allow the flames to proceed in an unhindered manner. In order to explore sample size and chamber size effects, two large chambers were developed at NASA GRC under the Flame Prevention, Detection and Suppression (FPDS) project. The first was an existing vacuum facility, VF-13, located at NASA John Glenn Research Center. This 6350 liter chamber could accommodate fuels sample lengths up to 2 m. However, operational costs and restricted accessibility limited the test program, so a second laboratory scale facility was developed in parallel. By stacking additional two chambers on top of an existing combustion chamber facility, this 81 liter Stacked-chamber facility could accommodate a 1.5 m sample length. The larger volume, more ideal environment of VF-13 was used to obtain baseline data for comparison with the stacked chamber facility. In this way, the stacked chamber facility was intended for long term testing, with VF-13 as the proving ground. Four different solid fuels (adding machine paper, poster paper, PMMA plates, and Nomex fabric) were tested with fuel sample lengths up to 2 m. For thin samples (papers) with widths up to 5 cm, the flame reached a steady state length, which demonstrates that flame length may be stabilized even when the edge effects are reduced. For the thick PMMA plates, flames reached lengths up to 70 cm but were highly energetic and restricted by oxygen depletion. Tests with the Nomex fabric confirmed that the cyclic flame phenomena, observed in small facility tests, continued over longer sample. New features were also observed at the higher oxygen/pressure conditions available in the large chamber. Comparison of flame behavior between the two facilities under identical conditions revealed disparities, both qualitative and quantitative. This suggests that, in certain ranges of controlling parameters, chamber size and shape could be one of the parameters that affect the material flammability. If this proves to be true, it may limit the applicability of existing flammability data.

Kleinhenz, Julie↗

Numerical Simulation of One- and Two-Phase Flows in Propulsion Systems

In this report, we present some results of problems investigated during joint research between the Hampton University Fluid Mechanics and Acoustics Laboratory (HU/FM&AL), NASA GRC, and the LaRC Hyper-X Program. This work is supported by joint research between the NASA GRC and the Institute of Mechanics at Moscow State University (IM/MSU) in Russia under a CRDF grant. The main areas of current scientific interest of the HU/FM&AL include an investigation of the proposed and patented advanced methods for aircraft engine thrust and noise benefits. These methods are based on nontraditional 3D corrugated and composite nozzle, inlet, propeller and screw designs such as a Bluebell and Telescope nozzle, Mobius-shaped screw, etc. This is the main subject of our other projects, of which one is presented at the current conference. Here we analyze additional methods for exhaust jet noise reduction without essential thrust loss and even with thrust augmentation. Such additional approaches are: (1) to add some solid, fluid, or gas mass at discrete locations to the main supersonic gas stream to minimize the negative influence of strong shock waves formed in propulsion systems. This mass addition may be accompanied by heat addition to the main stream as a result of the fuel combustion or by cooling of this stream as a result of the liquid mass evaporation and boiling; (2) Use of porous or permeable nozzles and additional shells at the nozzle exit for preliminary cooling of the hot jet exhaust and pressure compensation for off-design conditions (so-called continuous ejector with small mass flow rate); and (3) to propose and analyze new effective methods of fuel injection into the flow stream in air-breathing engines. The research is focused on a wide regime of problems in the propulsion field as well as in experimental testing and theoretical and numerical simulation analyses for advanced aircraft and rocket engines. The FM&AL Team uses analytical methods, numerical simulations, and experimental tests at the Hampton University campus, NASA, and IM/MSU. The main outcomes during this reporting period are: (1) Publications: The AIAA Paper #01-1893 has been accepted for the AIAA/NAL-NASDA-ISAS 10th International Space Planes and Hypersonic Systems and Technologies Conference, 24-27 April 2001, Kyoto, Japan. The AIAA Paper #01-3204 has been accepted for presentation at the 37th AIAA/ASME/SAE/ASEE Joint Propulsion Conference being held on 08-11 July, in Salt Lake City, UT; (2) Grants and proposals: The HU/FM&AL was awarded the NASA grant NAG3-2495 in October 2000 and the CRDF award was granted to the NASA GRC-HU FM&AL and IM/MSU (Russia) in July 2000. A solicited proposal was submitted for the NASA NRA-01GRC-02 competition and two unsolicited proposals to NASA are in preparation; (3) Theory and numerical simulations: Analytical theory, numerical simulation, comparison of theoretical with experimental results, and modification of theoretical approaches, models, grids, etc. Such investigations have been conducted for three main problems: (a) Combustion efficiency optimization in the half-duct combustor system; (b) Drag reduction effects for blunt bodies with solid needles; and (c) Solid particle injection from the butt-end against a supersonic flow. The NASA CFL3D, HU/FM&AL, and IM/MSU GODUNOV codes were used and modified for solution of these problems. The codes are based on full Euler and Navier-Stokes solvers with and without nonequilibrium oxygen-nitrogen and air-hydrogen chemical reactions in laminar and turbulent gas flow regimes; (4) Experimental tests: In the IM/MSU supersonic wind tunnel, experimental tests were conducted with different number, location, and geometric parameters of solid needles mounted at the front of the butt-end in supersonic flow. Optimal parameters were determined that provide minimal butt-end drag in the stationary flow regime. Experimental and numerical simulation results are in good agreement; (5) Student Research Activity: Involvement of one graduate and two undergraduate students as research assistants in the current project.

Gilinsky, Mikhail↗

Investigating Stable Oxygen and Carbon Isotopic Variability in Speleothem Records Over the Last Millennium Using Multiple Isotope-Enabled Climate Models

The incorporation of water isotopologues into the hydrology of general circulation models (GCMs) facilitates the comparison between modeled and measured proxy data in paleoclimate archives. However, the variability and drivers of measured and modeled water isotopologues, as well as the diversity of their representation in different models, are not well constrained. Improving our understanding of this variability in past and present climates will help to better constrain future climate change projections and decrease their range of uncertainty. Speleothems are a precisely datable terrestrial paleoclimate archives and provide well-preserved (semi-)continuous multivariate isotope time series in the lower latitudes and mid-latitudes and are therefore well suited to assess climate and isotope variability on decadal and longer timescales. However, the relationships of speleothem oxygen and carbon isotopes to climate variables are influenced by site-specific parameters, and their comparison to GCMs is not always straightforward. Here we compare speleothem oxygen and carbon isotopic signatures from the Speleothem Isotopes Synthesis and Analysis database version 2 (SISALv2) to the output of five different water-isotope-enabled GCMs (ECHAM5-wiso, GISS-E2-R, iCESM, iHadCM3, and isoGSM) over the last millennium (850–1850 CE). We systematically evaluate differences and commonalities between the standardized model simulation outputs. The goal is to distinguish climatic drivers of variability for modeled isotopes and compare them to those of measured isotopes. We find strong regional differences in the oxygen isotope signatures between models that can partly be attributed to differences in modeled surface temperature. At low latitudes, precipitation amount is the dominant driver for stable water isotope variability; however, at cave locations the agreement between modeled temperature variability is higher than for precipitation variability. While modeled isotopic signatures at cave locations exhibited extreme events coinciding with changes in volcanic and solar forcing, such fingerprints are not apparent in the speleothem isotopes. This may be attributed to the lower temporal resolution of speleothem records compared to the events that are to be detected. Using spectral analysis, we can show that all models underestimate decadal and longer variability compared to speleothems (albeit to varying extents). We found that no model excels in all analyzed comparisons, although some perform better than the others in either mean or variability. Therefore, we advise a multi-model approach whenever comparing proxy data to modeled data. Considering karst and cave internal processes, e.g., through isotope-enabled karst models, may alter the variability in speleothem isotopes and play an important role in determining the most appropriate model. By exploring new ways of analyzing the relationship between the oxygen and carbon isotopes, their variability, and co-variability across timescales, we provide methods that may serve as a baseline for future studies with different models using, e.g., different isotopes, different climate archives, or different time periods.

isotope-enabled simulations↗