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.

SourceResearch connectionFrom yearThrough year

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

At least 595 records · Page 33

All sources NASA NTRS DOE OSTI

NASA NTRS2006

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↗

NASA NTRSDate not supplied

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↗

DOE OSTI2024

Harnessing complexity: Nonlinear optical phenomena in L-shapes, nanocrescents, and split-ring resonators

Here, we conduct systematic studies of the optical characteristics of plasmonic nanoparticles that exhibit C 2v symmetry. In particular, we analyze three distinct geometric configurations: an L-type shape, a crescent, and a split-ring resonator shaped like the Greek letter π. Optical properties are examined using the finite-difference time-domain method. It is demonstrated that all three shapes exhibit two prominent plasmon modes associated with the two axes of symmetry. This is in addition to a wide range of resonances observed at high frequencies corresponding to quadrupole modes and peaks due to sharp corners. Next, to facilitate nonlinear analysis, we employ a semiclassical hydrodynamic model, where the electron pressure term is explicitly accounted for. This model goes beyond the standard Drude description and enables capturing nonlocal and nonlinear effects. Employing this model enables us to rigorously examine the second-order angular resolved nonlinear optical response of these nanoparticles in each of the three configurations. Two pumping regimes are considered, namely, continuous wave (CW) and pulsed excitations. For CW pumping, we explore the properties of the second harmonic generation (SHG). Polarization and angle-resolved SHG spectra are obtained, revealing strong dependence on the nanoparticle geometry and incident wave polarization. The C 2v symmetry is shown to play a key role in determining the polarization states and selection rules of the SHG signal. For pulsed excitations, we discuss the phenomenon of broadband terahertz (THz) generation induced by the difference-frequency generation . It is shown that the THz emission spectra exhibit unique features attributed to the plasmonic resonances and symmetry of the nanoparticles. The polarization of the generated THz waves is also examined, revealing interesting patterns tied to the nanoparticle geometry. To gain deeper insight, we propose an analytical theory that agrees very well with the numerical experiments. The theory shows that the physical origin of the THz radiation is the mixing of various frequency components of the fundamental pulse by the second-order nonlinear susceptibility. An expression for the far-field THz intensity is derived in terms of the incident pulse parameters and the nonlinear response tensor of the nanoparticle. The results presented in this work offer new insights into the linear and nonlinear optical properties of nanoparticles with C 2v symmetry. The demonstrated strong SHG response and efficient broadband THz generation hold great promise for applications in nonlinear spectroscopy, nanophotonics, and optoelectronics. The proposed theoretical framework also provides a valuable tool for understanding and predicting the nonlinear behavior of other related nanostructures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2016

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↗

NASA NTRS2001

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.↗

NASA NTRS2017

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.↗

NASA NTRS2018

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↗

NASA NTRS1995

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.↗

NASA NTRS2008

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↗

DOE OSTI2024

Solid State Quantum Refrigeration Superconducting, Absorption and Measurement Based (Final Technical Report)

During this DOE grant, DE-SC0017890, in place for the past six years, all proposed research was carried out and published in peer-reviewed papers, as well as other projects that emerged during the research. In that effort the research team accomplished all proposed research, as well as many closely related research projects discovered and conceived of during the grant. These works included “Efficient Quantum Measurement Engines”, a work published in Physical Review Letters, giving a theory of quantum measurement-based engines, which uses quantum measurement as a resource. These engines are designed to efficiently convert energy from the stochastic quantum measurement process into useful work. Further publications include “Experimental Realization of a Quantum Dot Energy Harvester”, a joint theory and experimental work in collaboration with the group of Charles Smith in Cambridge, UK, as well as long time theoretical collaborators, Rafael Sánchez and Björn Sothmann. This work, featured as an Editor’s Suggestion in Physical Review Letters, realized an earlier theoretical proposal of ours, whereby two resonant tunneling quantum dots are connected to a central electronic cavity that is heated by a hot energy source. We also published “Superconducting Quantum Refrigerator: Breaking and Rejoining Cooper Pairs with Magnetic Field Cycles” a work done in collaboration with experimentalist Francesco Giazotto from ENS Pisa, Italy, which also resulted in a patent. This paper, published in Phys. Rev. Applied, advanced the concept of a cyclic fridge based on the normal/superconducting phase transition together with layered materials separated by tunnel junctions. We also completed the proposed research on a heat transistor, publishing “Thermal transistor and thermometer based on Coulomb-coupled conductors”, carried out as a collaboration between my group and theorists Splettstoesser (Lund U., Sweden), Sothmann (U. Duisburg-Essen, Germany), and Sánchez (U. Autónoma de Madrid, Spain). We carried out an analysis of a quantum coupled to a quantum point contact as a sensitive thermometer and heat transistor. We found the optimal statistical estimator for the temperature and compared it with experiments on the same type of devices. We also investigated autonomous quantum absorption refrigerators using quantum dots to cool by using a very hot thermal reservoir to drive heat between two other reservoirs. In the article “Quantifying the quantum heat contribution from a driven superconducting circuit”, we demonstrated that for a driven superconducting circuit, we showed heat flow provided by a hot source to the qubit can be switched on and off by varying external parameters, the frequency and the intensity of the driving. In the work “Stochastic thermodynamic cycles of a mesoscopic thermoelectric engine”, we reconsidered the autonomous thermoelectric heat engine in terms of underlying cycles. Rather than periodic behavior, the cycles were stochastic in nature. Nevertheless, by undertaking a graph theoretical analysis of the elementary transport processed, great quantitative and qualitative insight could be found. We also considered the quantum measurement process and showed that a quantum version of Maxwell’s demon could be related to the work extraction of a quantum system, closely related to arrow-of-time measures for quantum measurement, as described in our article “Thermodynamics of quantum measurement and Maxwell's demon's arrow of time”. This work was selected in Phys. Rev. A as an Editor’s Suggestion. A recent preprint titled “Cyclic Superconducting Quantum Refrigerators Using Guided Fluxon Propagation” accomplished an important piece of this grant: to propose a new kind of quantum refrigerator using the dynamics of fluxons in a type II superconductor. This invention envisioned a race-track type geometry where fluxons are confined. By applying a gradient of magnetic field together with electric current in a Corbino geometry, the circulating fluxons can actively cool a cold reservoir, realizing a new type of cyclic superconducting refrigerator. We also investigated the possibility of thermal control from different points of view. The application of quantum measurement to the system gives a new kind of control on the system of interest – we have pioneered this approach and shown that measurement can boost the thermal power of quantum engines as described in “Continuous measurement boosted adiabatic quantum thermal machines”. The ability to have heat flows on demand is an outstanding challenge, and we have provided new solutions to this problem in Thermal control across a chain of electronic nanocavities” for a chain of electron cavities using gating voltage control. The control methods using qubit/qubit coupling to create absorption fridges at their most fundamental level have also been developed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Structural and mechanical properties of monolayer amorphous carbon and boron nitride

Amorphous materials exhibit various characteristics that are not featured by crystals and can sometimes be tuned by their degree of disorder (DOD). Here, we report results on the mechanical properties of monolayer amorphous carbon (MAC) and monolayer amorphous boron nitride (maBN) with different DOD. The pertinent structures are obtained by kinetic-Monte-Carlo (kMC) simulations using machine-learning potentials (MLP) with density-functional-theory (DFT)-level accuracy. An intuitive order parameter, namely the areal fraction F x occupied by crystallites within the continuous random network, is proposed to describe the DOD. We find that F x captures the essence of the DOD: Samples with the same F x but different sizes and arrangements of crystallites, obtained using two distinct kMC procedures, have virtually identical radial distributions functions as well as bond-length and bond-angle distributions. Furthermore, by simulating the fracture process with molecular dynamics, we found that the mechanical responses of MAC and maBN before fracture are mainly determined by F x and are insensitive to the sizes and specific arrangements and to some extent the numbers and area distributions of the crystallites. The behavior of cracks in the two materials is analyzed and found to mainly propagate in meandering paths in the CRN region and to be influenced by crystallites in distinct ways that toughen the material. Furthermore, the present results reveal the relation between structure and mechanical properties in amorphous monolayers and may provide a universal toughening strategy for 2D materials.

2-dimensional systems↗

NASA NTRSDate not supplied

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↗

NASA NTRS2018

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 ↗

NASA NTRS2011

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↗

NASA NTRS2001

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↗

NASA NTRS2022

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↗

NASA NTRS2007

Development of a Tunable Electromechanical Acoustic Liner for Engine Nacelles

This report describes the development of a tunable electromechanical Helmholtz resonator (EMHR) for engine nacelles using smart materials technology. This effort addresses both near-term and long-term goals for tunable electromechanical acoustic liner technology for the Quiet Aircraft Technology (QAT) Program. Analytical models, i.e. lumped element model (LEM) and transfer matrix (TM) representation of the EMHR, have been developed to predict the acoustic behavior of the EMHR. The models have been implemented in a MATLAB program and used to compare with measurement results. Moreover, the prediction performance of models is further improved with the aid of parameter extraction of the piezoelectric backplate. The EMHR has been experimentally investigated using standard two-microphone method (TMM). The measurement results validated both the LEM and TM models of the EMHR. Good agreement between predicted and measured impedance is obtained. Short- and open circuit loads define the limits of the tuning range using resistive and capacitive loads. There is approximately a 9% tuning limit under these conditions for the non-optimized resonator configuration studied. Inductive shunt loads result in a 3 degree-of-freedom DOF) system and an enhanced tuning range of over 20% that is not restricted by the short- and open-circuit limits. Damping coefficient ' measurements for piezoelectric backplates in a vacuum chamber are also performed and indicate that the damping is dominated by the structural damping losses, such as compliant boundaries, and other intrinsic loss mechanisms. Based on models of the EMHR, a Pareto optimization design of the EMHR has been performed for the EMHR with non-inductive loads. The EMHR with non-inductive loads is a 2DOF system with two resonant fiequencies. The tuning ranges of the two resonant frequencies of the EMHR with non-inductive loads cannot be optimized simultaneously; a trade-off (i.e., a Pareto solution) must be reached. The Pareto solution provides the information for a designer that shows how design trade-offs can be used to satisfy specific design requirements. The optimization design of the EMHR with inductive loads aims at optimal tuning of these three resonant fiequencies. The results indicate that it is possible to keep the acoustic reactance of the resonator close to a constant over a given frequency range. An effort to mimic the second layer of the NASA 2DOF liner using a piezoelectric composite diaphragm has been made. The optimal acoustic reactance of the second layer of the NASA 2DOF liner is achieved using a thin PVDF composite diaphragm, but matching the acoustic resistance requires further investigation. Acoustic energy harvesting is achieved by connecting the EMHR to an energy reclamation circuit that converts the ac voltage signal across the piezoceramic to a conditioned dc signal. Energy harvesting experiment yields 16 m W continuous power for an incident SPL of 153 dB. Such a level is sufficient to power a variety of low power electronic devices. Finally, technology transfer has been achieved by converting the original NASA ZKTL FORTRAN code to a MATLAB code while incorporating the models of the EMHR. Initial studies indicate that the EMHR is a promising technology that may enable lowpower, light weight, tunable engine nacelle liners. This technology, however, is very immature, and additional developments are required. Recommendations for future work include testing of sample EMHR liner designs in NASA Langley s normal incidence dual-waveguide and the grazing-incidence flow facility to evaluating both the impedance characteristics as well as the energy reclamation abilities. Additional design work is required for more complex tuning circuits with greater performance. Poor electromechanical coupling limited the electromechanical tuning capabilities of the proof of concept EMHR. Different materials than those studies and perhaps novel composite material systems may dramatically improvehe electromechanical coupling. Such improvements are essential to improved mimicking of existing double layer liners.

Liu, Fei↗

NASA NTRS2001

Active Control of Rotating Stall Demonstrated for a Multistage Compressor With Inlet Distortion

Aircraft compressors can suffer debilitating consequences as a result of rotating stall and surge events caused by inlet distortions. This is particularly true of aircraft during takeoff, when the compressor is operating at peak performance close to the surge line. Significant research has been conducted by the NASA Glenn Research Center in the area of compressor stability enhancement through active and passive control methods. Most recently, an experiment was conducted at the Wright Patterson Air Force Base Research Laboratory on a two-stage fan with inlet guide vanes and inlet distortion. In this joint Small Business Innovation Research effort between Scientific Systems and Glenn, control of rotating stall was demonstrated in a multistage transonic fan. This twostage fan with inlet guide vanes was tested under clean and distorted inlet conditions. The compressor was also configured with a circumferential distortion screen capable of 180 of distortion and with 14 high-velocity injectors upstream of the first rotor. Twelve of these injectors could oscillate up to frequencies of 450 Hz. The additional two injectors were located next to each other and were used in concert with each other as a single, on/off, high-authority actuator. In a first test of injection in this multistage environment, 12 of the valves were opened 50 percent of their full stroke to assess steady injection through the compressor. This baseline injection is shown in the compressor characteristic of the following figure, and stall margin improvements are tracked from this baseline condition. The compressor was then tested with clean inlet conditions using 12 injectors and active control. Pressure disturbances were tracked before rotating stall, and a constant gain control scheme reduced the stalling mass flow by 10.8 percent over the baseline. With the distortion screen present in the inlet, a pole-zero cancellation control scheme was used to achieve a 6.4-percent decrease in stalling mass flow. These improvements also are shown in the figure. In a final experiment, actively controlled, high-frequency injection from the 12 valves was used in conjunction with the high-authority actuators. In this test, the stalling mass flow of the compressor was reduced by 27 percent as indicated in the graph. These results were obtained by injecting less than 2 percent of the total compressor throughflow into the rotor tip region via 14 injection ports. These results mark the first successful demonstration of actively controlled air injection as a stall-control strategy for multistage compressors operating at speeds typical of an actual gas turbine engine. A goal of continuing research is to determine the combination of air-injection parameters and control strategies that are most effective in providing stall control for both clean and distorted inlet flow conditions for multistage environments. Other goals include the demonstration of stall control at many locations along the core compressor and development and application of active stall control strategies that will be integral flightworthy components of onboard engine hardware.

VanSchalkwyk, Christian↗