Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Incremental Computing”

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.

504 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Computational Investigation of Nominally-Orthogonal Pneumatic Active Flow Control for High-Lift Systems

We explore the feasibility of using nominally-orthogonal jets as active aerodynamic load control for multi-element high-lift systems, and whether the nominally-orthogonal jets can offer a variety of performance improvements. These nominally-orthogonal jets inject momentum normal to the airfoil surface near the flap trailing edge, where they create a vortex that entrains flow from the opposing side and change the airfoil circulation. Lift-enhancement opportunities of trailing edge nominally-orthogonal jets have previously been studied by Malavard et al. and Blaylock et al. on single-element airfoils; however, their effect on drag was not thoroughly investigated. In this study, we investigate two- dimensional nominally-orthogonal jet effects on both lift and drag on a two-element airfoil, NLR7301. We utilize Chimera Grid Tools to generate structured curvilinear overset grids, and the Reynolds-averaged Navier-Stokes solver OVERFLOW-2 to solve for the flow field around the airfoil. We perform various computational sensitivity studies on the baseline airfoil without a jet to validate computational results against benchmark experimental data. Using a Chimera overset grid topology, we demonstrate a similar lift-enhancement effect between a nominally-orthogonal jet and a nominally-orthogonal physical tab employed at the same location on the studied airfoil. After we introduce the nominally-orthogonal jet concept, we investigate nominally-orthogonal jets with various momentum coefficient settings, C(sub µ) = 0.00 − 0.04 and present a lift-enhancement relationship ∆C(sub l) ≃ 3.59 (C(sub µ)) for this airfoil. We discuss that utilizing a nominally-orthogonal jet with C(sub µ) = 0.01 can shift the linear region of the lift curve by a ∆C(sub l) = 0.36 for the pressure side jet and by a ∆C(sub l) = −0.27 for the suction side jet. Employing a nominally-orthogonal jet is also shown effective in altering the drag. To study the impact, we carry out a drag decomposition study in the form of drag polars. We show for a given C(sub l) = 2.50, a nominally-orthogonal jet with C(sub µ) = 0.01 on the pressure and suction side of the airfoil results in 113 drag count decrements and 41 drag count increments, respectively, compared to the baseline airfoil with no jet. These results show that large and controllable changes in aerodynamic performance can be achieved by relatively small active flow control inputs using the nominally-orthogonal jets presented in this study.

Nominally-Orthogonal↗
NASA NTRS2012

Personal Computer Transport Analysis Program

The Personal Computer Transport Analysis Program (PCTAP) is C++ software used for analysis of thermal fluid systems. The program predicts thermal fluid system and component transients. The output consists of temperatures, flow rates, pressures, delta pressures, tank quantities, and gas quantities in the air, along with air scrubbing component performance. PCTAP s solution process assumes that the tubes in the system are well insulated so that only the heat transfer between fluid and tube wall and between adjacent tubes is modeled. The system described in the model file is broken down into its individual components; i.e., tubes, cold plates, heat exchangers, etc. A solution vector is built from the components and a flow is then simulated with fluid being transferred from one component to the next. The solution vector of components in the model file is built at the initiation of the run. This solution vector is simply a list of components in the order of their inlet dependency on other components. The component parameters are updated in the order in which they appear in the list at every time step. Once the solution vectors have been determined, PCTAP cycles through the components in the solution vector, executing their outlet function for each time-step increment.

DiStefano, Frank, III↗
NASA NTRS2014

Advanced Modeling and Uncertainty Quantification for Flight Dynamics; Interim Results and Challenges

As part of the NASA Vehicle Systems Safety Technologies (VSST), Assuring Safe and Effective Aircraft Control Under Hazardous Conditions (Technical Challenge #3), an effort is underway within Boeing Research and Technology (BR&T) to address Advanced Modeling and Uncertainty Quantification for Flight Dynamics (VSST1-7). The scope of the effort is to develop and evaluate advanced multidisciplinary flight dynamics modeling techniques, including integrated uncertainties, to facilitate higher fidelity response characterization of current and future aircraft configurations approaching and during loss-of-control conditions. This approach is to incorporate multiple flight dynamics modeling methods for aerodynamics, structures, and propulsion, including experimental, computational, and analytical. Also to be included are techniques for data integration and uncertainty characterization and quantification. This research shall introduce new and updated multidisciplinary modeling and simulation technologies designed to improve the ability to characterize airplane response in off-nominal flight conditions. The research shall also introduce new techniques for uncertainty modeling that will provide a unified database model comprised of multiple sources, as well as an uncertainty bounds database for each data source such that a full vehicle uncertainty analysis is possible even when approaching or beyond Loss of Control boundaries. Methodologies developed as part of this research shall be instrumental in predicting and mitigating loss of control precursors and events directly linked to causal and contributing factors, such as stall, failures, damage, or icing. The tasks will include utilizing the BR&T Water Tunnel to collect static and dynamic data to be compared to the GTM extended WT database, characterizing flight dynamics in off-nominal conditions, developing tools for structural load estimation under dynamic conditions, devising methods for integrating various modeling elements into a real-time simulation capability, generating techniques for uncertainty modeling that draw data from multiple modeling sources, and providing a unified database model that includes nominal plus increments for each flight condition. This paper presents status of testing in the BR&T water tunnel and analysis of the resulting data and efforts to characterize these data using alternative modeling methods. Program challenges and issues are also presented.

Hyde, David C.↗
NASA NTRS2015

A Module for Assimilating Hyperspectral Infrared Retrieved Profiles into the Gridpoint Statistical Interpolation System for Unique Forecasting Applications

Hyperspectral infrared sounder radiance data are assimilated into operational modeling systems however the process is computationally expensive and only approximately 1% of available data are assimilated due to data thinning as well as the fact that radiances are restricted to cloud-free fields of view. In contrast, the number of hyperspectral infrared profiles assimilated is much higher since the retrieved profiles can be assimilated in some partly cloudy scenes due to profile coupling other data, such as microwave or neural networks, as first guesses to the retrieval process. As the operational data assimilation community attempts to assimilate cloud-affected radiances, it is possible that the use of retrieved profiles might offer an alternative methodology that is less complex and more computationally efficient to solve this problem. The NASA Short-term Prediction Research and Transition (SPoRT) Center has assimilated hyperspectral infrared retrieved profiles into Weather Research and Forecasting Model (WRF) simulations using the Gridpoint Statistical Interpolation (GSI) System. Early research at SPoRT demonstrated improved initial conditions when assimilating Atmospheric Infrared Sounder (AIRS) thermodynamic profiles into WRF (using WRF-Var and assigning more appropriate error weighting to the profiles) to improve regional analysis and heavy precipitation forecasts. Successful early work has led to more recent research utilizing WRF and GSI for applications including the assimilation of AIRS profiles to improve WRF forecasts of atmospheric rivers and assimilation of AIRS, Cross-track Infrared and Microwave Sounding Suite (CrIMSS), and Infrared Atmospheric Sounding Interferometer (IASI) profiles to improve model representation of tropopause folds and associated non-convective wind events. Although more hyperspectral infrared retrieved profiles can be assimilated into model forecasts, one disadvantage is the retrieved profiles have traditionally been assigned the same error values as the rawinsonde observations when assimilated with GSI. Typically, satellitederived profile errors are larger and more difficult to quantify than traditional rawinsonde observations (especially in the boundary layer), so it is important to appropriately assign observation errors within GSI to eliminate potential spurious innovations and analysis increments that can sometimes arise when using retrieved profiles. The goal of this study is to describe modifications to the GSI source code to more appropriately assimilate hyperspectral infrared retrieved profiles and outline preliminary results that show the differences between a model simulation that assimilated the profiles as rawinsonde observations and one that assimilated the profiles in a module with the appropriate error values.

Berndt, Emily↗
NASA NTRS2017

Applying Graph Theory to Problems in Air Traffic Management

Graph theory is used to investigate three different problems arising in air traffic management. First, using a polynomial reduction from a graph partitioning problem, it is shown that both the airspace sectorization problem and its incremental counterpart, the sector combination problem are NP-hard, in general, under several simple workload models. Second, using a polynomial time reduction from maximum independent set in graphs, it is shown that for any fixed e, the problem of finding a solution to the minimum delay scheduling problem in traffic flow management that is guaranteed to be within n1-e of the optimal, where n is the number of aircraft in the problem instance, is NP-hard. Finally, a problem arising in precision arrival scheduling is formulated and solved using graph reachability. These results demonstrate that graph theory provides a powerful framework for modeling, reasoning about, and devising algorithmic solutions to diverse problems arising in air traffic management.

precision arrival scheduling↗
NASA NTRS2000

Comparison of Analytical and CAA Solution(s): Unsteady Response of a Rectilinear Swept Cascade to an Incident Gust

This problem is aimed at testing the ability of computational aeroacoustic methods to calculate the acoustic pressure field generated as a result of the interaction of a convected harmonic gust with a rectilinear cascade of swept flat plates. All relevant geometric information is the same as the unsteady response of an isolated finite span swept airfoil to an incident gust problem with the exception of the sweep angle alpha, which is now taken to be a variable. Assume a cascade stagger angle of zero (with the x-axis aligned with the chord) and a gap-to-chord ratio of h/c = 1. The mean flow Mach number M(sub 0), gust frequency omega, gust amplitude A, and chordwise wavenumber k(sub x) are the same. But for the cascade problem take k(sub y) = pi and k(sub z) = 0. The appropriate physical scales are also the same, as is the governing equation and boundary conditions. Naturally, the impermeability condition now applies to the entire cascade. For this problem, determine the amplitude of upstream-radiated rms acoustic pressure as a function of the sweep angle at the specified frequency. Specifically, show the variations of rms acoustic pressure amplitude at the upstream location (-5c, 0, l/2) for sweep angle alpha in the range (0.0 deg, 30.0 deg). Use sweep angle increments no larger than 2.5 deg. Express the results in dB using the rms pressure value for alpha = 0.0 deg as the reference level. The appropriate FORTRAN output statement should read "WRITE(IUNIT,*) alpha, 20-log10 (p(sub rms)(alpha)/p(sub rms)(0))".

Envia, Edmane↗
NASA NTRS2019

An Evaluation and Recommendations for Further CFD Research Based on the NASA Common Research Model (CRM) Analysis from the AIAA Drag Prediction Workshop (DPW) Series

The six Drag Prediction Workshops (DPW) have drawn a total of 125 participants submitting 225 data entries. Over half of these occurred in the last three workshops dealing with the NASA Common Research Model (CRM). Participants have ranged from North America, Europe, Asia, and South America, representing Government, Industry, Academia, and Commercial/Vendors. The DPW Series has generated a tremendous amount of CFD data, freely available to the public domain. A continuous improvement of the results as measured by: the reduction of the spread of drag predictions for the specified “nominal” cruise point, accuracy of drag increment, and the angle-of-attack sweep behavior, has been noted. Although it has been difficult separating out the parts and pieces that contribute to errors in CFD drag predictions, part of the improvements over the years can be attributed to grid resolution. The scatter in predicted drag in the continuum has been reduced dramatically from plus or minus 50 counts in DPW-I to around plus or minus 5 counts for DPW-V and - VI. During the six workshops, the grid size has grown dramatically. The average size of the medium wing-body meshes in DPW-I through DPW-IV have been 3.2, 5.4, 7.8 and 10.9 million, respectively. This represents a growth rate of approximately 17 percent per year during the eight years between DPW-I and DPW-IV. For DPW-VI, this number has grown from 25 to 50 million points for the various families of grids available. The finest level grids have increased steadily, from just over 3 million grid points in DPW-I to 225 million for the WBNP (Wing-Body-Nacelle-Pylon) in DPW-VI. While progress has been made, persistent problems with grid generation and separated flow prediction continue. The experiences and recommendations presented here are consistent with those espoused in the “CFD Vision 2030 Study: A Path to Revolutionary Computational Aerosciences,” NASA/CR-2014-218178.

Tinoco, Edward N.↗
NASA NTRSDate not supplied

Application of OpenFOAM to Plume Impingement in Space Environments

After 30 years of continuous human presence in low-earth orbit, NASA is returning to the moon and eventually will go to Mars. Travelling beyond low earth orbit requires NASA to learn how humans can live in Deep Space environments – beyond the protection of Earth’s magnetosphere and at distances from Earth that prevent a quick return in case of trouble. To this end, NASA is constructing the Lunar Gateway, an ISS-like space station to be put in orbit around the moon to act as a home base for Lunar exploration for NASA astronauts. The Gateway Lunar outpost will be built incrementally, via modules which will arrive at separate times and dock to the existing structure. The incremental addition of Gateway modules, and the docking of visiting vehicles, is achieved via a sequence of firings from the approaching body’s onboard reaction control system (RCS) thrusters to achieve the required approach trajectory. The typical hypergolic chemical RCS thrusters work by firing hot gases to produce adverse thrust and the needed change in velocity to safely finish the docking process. The exhaust gas from the RCS thrusters form plumes that expand into the vacuum of space and can impinge onto the outer surfaces of the Lunar Gateway, causing unwanted forces and moments, heat loads, sediment deposition, and in extreme cases, even surface erosion - all mechanisms that can damage the Lunar Gateway and must be minimized. Both permanent and visiting modules will have this RCS thruster exhaust impingement problem. This research aims to establish existing OpenFOAM solvers as a methodology for improving simulation techniques of rocket exhaust plume impingement in space environments. The flow structure of a plume in a space environment is complex; a plume that originates from a hypergolic chemical RCS thruster and expands into a vacuum will experience several regimes of rarefication. This range includes the continuum flow in the rocket nozzle through the fully rarefied free molecular flow further from the nozzle. The flow physics is different at these two extremes, and as such, the simulation approach for plumes is generally divided into a traditional computational fluid dynamics (CFD) simulation in and near the nozzle which is coupled to a subsequent direct simulation Monte Carlo (DSMC) simulation. At this time, the scope of this research is developing, verifying, and validating a method using existing solvers in the OpenFOAM framework for performing coupled CFD/DSMC calculations to determine the extent of plume impingement loading on generic space structures. This presentation will detail code-to-code comparisons between the hyStrath dsmcFoam+ solver, developed using OpenFOAM and available as open-source, and NASA’s in-house DSMC Analysis Code (DAC). Comparisons to several open-source publication findings using DAC [3,4] are presented, and advantages of using an OpenFOAM based solver are also discussed. The presentation concludes with a discussion of future work, and a plan for coupling the dsmcFoam+ solver with CFD simulations of chemical rocket engines for unified coupled plume simulation.

DSMC↗
NASA NTRS1993

Development and use of computational techniques in Army Aviation research and development programs for crash resistant helicopter technology

During the 1960's over 30 full-scale aircraft crash tests were conducted by the Flight Safety Foundation under contract to the Aviation Applied Technology Directorate (AATD) of the U.S. Army Aviation Systems Command (AVSCOM). The purpose of these tests were to conduct crash injury investigations that would provide a basis for the formulation of sound crash resistance design criteria for light fixed-wing and rotary wing aircraft. This resulted in the Crash Survival Design Criteria Designer's Guide which was first published in 1967 and has been revised numerous times, the last being in 1989. Full-scale aircraft crash testing is an expensive way to investigate structural deformations of occupied spaces and to determine the decelerative loadings experienced by occupants in a crash. This gave initial impetus to the U.S. Army to develop analytical methods to predict the dynamic response of aircraft structures in a crash. It was believed that such analytical tools could be very useful in the preliminary design stage of a new helicopter system which is required to demonstrate a level of crash resistance and had to be more cost effective than full-scale crash tests or numerous component design support tests. From an economic point of view, it is more efficient to optimize for the incorporation of crash resistance features early in the design stage. However, during preliminary design it is doubtful if sufficient design details, which influence the exact plastic deformation shape of structural elements, will be available. The availability of simple procedures to predict energy absorption and load-deformation characteristics will allow the designer to initiate valuable cost, weight, and geometry tradeoff studies. The development of these procedures will require some testing of typical specimens. This testing should, as a minimum, verify the validity of proposed procedures for providing pertinent nonlinear load-deformation data. It was hoped that through the use of these analytical models, the designer could optimize aircraft design for crash resistance from both a weight and cost increment standpoint, thus enhancing the acceptance of the design criteria for crash resistance.

Burrows, Leroy T.↗
NASA NTRS2002

Analysis of High Order Difference Methods for Multiscale Complex Compressible Flows

Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes with incremental studies was initiated. Here we further refine the analysis on, and improve the understanding of the adaptive numerical dissipation control strategy. Basically, the development of these schemes focuses on high order nondissipative schemes and takes advantage of the progress that has been made for the last 30 years in numerical methods for conservation laws, such as techniques for imposing boundary conditions, techniques for stability at shock waves, and techniques for stable and accurate long-time integration. We concentrate on high order centered spatial discretizations and a fourth-order Runge-Kutta temporal discretizations as the base scheme. Near the bound-aries, the base scheme has stable boundary difference operators. To further enhance stability, the split form of the inviscid flux derivatives is frequently used for smooth flow problems. To enhance nonlinear stability, linear high order numerical dissipations are employed away from discontinuities, and nonlinear filters are employed after each time step in order to suppress spurious oscillations near discontinuities to minimize the smearing of turbulent fluctuations. Although these schemes are built from many components, each of which is well-known, it is not entirely obvious how the different components be best connected. For example, the nonlinear filter could instead have been built into the spatial discretization, so that it would have been activated at each stage in the Runge-Kutta time stepping. We could think of a mechanism that activates the split form of the equations only at some parts of the domain. Another issue is how to define good sensors for determining in which parts of the computational domain a certain feature should be filtered by the appropriate numerical dissipation. For the present study we employ a wavelet technique introduced in as sensors. Here, the method is briefly described with selected numerical experiments.

Sjoegreen, Bjoern↗
NASA NTRSDate not supplied

Understanding and Utilizing PBL Height Data from Multiple Observing Systems in the GEOS System

The accuracy of PBL height simulation is a key issue in many applications including forecasting near surface meteorology and air quality, however, it is a very challenging problem due to the lack of not only comprehensive, global Planetary Boundary Layer (PBL) observations but also a strategy and infrastructure to utilize PBL height data from a variety of sensors. Following the designation of PBL as an incubation class observable in the 2017 Decadal Survey, the PBL Incubation Study Team Report [14] made clear that “a future global PBL observing system requires modeling and data assimilation as essential components.” There is an urgent need for global modeling development in order to utilize Program of Record (POR) observations, assess their impacts, and identify gaps to be filled by future PBL missions. Our overall objective is to develop PBL data assimilation capabilities in the NASA Global Earth Observing System (GEOS), focusing on PBL height from multiple observing systems, to support the assessment and use of future PBL observations. The NASA GEOS system is composed of the GEOS global atmospheric general circulation model (AGCM) and the atmospheric data assimilation system (ADAS). The PBL parameterizations include the “Lock” K-profile scheme driven by surface and cloud-top buoyancy fluxes ([4]), and the “Louis” local scheme for stable conditions based on the Richardson number ([5]). Above the mixed layer defined by the Lock surface plume, shallow cumulus convection is represented by the mass flux scheme of [9]. Additional parameterizations are summarized in [1]. The ADAS employs the hybrid 4D Ensemble- Variational (EnVar) configuration ([15]), with the ensemble providing flow-dependent background error covariance information. The resultant analysis increments are fed back to the forecast model through the 4D incremental analysis update (IAU) approach ([11]). In this study, PBL height data are being or have been generated from radiosondes, GNSS RO, satellite (CATS, CALIPSO and ICESat-2) and ground-based (MPLNET) lidars, and wind profiler. Investigations have been conducted to specify quality marks for PBL height retrievals for the data assimilation purpose. These PBL height data have different strengths and weaknesses ([2], [3], [6], [7], [8], [10]), and the satellite PBL height data provide better global coverage and complement in-situ PBL height data. Radiosondes offer high accuracy and in situ measurement of temperature and humidity profiles, but with poor spatio-temporal sampling. The in-situ observing systems like MPLNET and wind profiler provide long history of PBL height records at each station. The GNSS RO based PBL height is retrieved based on the sharp gradients in refractivity profile that represent the fine vertical structure of temperature and moisture changes above the PBL. However, not all RO refractivity profiles reach the surface depending on location and regime, and RO refractivity retrievals can be negatively biased below 2km. The PBL height data from satellite lidars provide high resolution along track PBL height retrievals, but over land they are affected by previous day convective PBL aerosol and strongly associated with mixing layer and retrievals cannot be made below thick, attenuating clouds. A successful assimilation of PBL height data requires a thorough understanding of the observing method and the retrieval algorithm for each observing system in order to use the PBL height data from multiple observing systems properly. Due to the sensitivity of PBL height data to the observing method and choice of algorithm, it is important to use a model definition appropriate for each observation type to compute differences between PBL height data and model PBL height (OmFs). The GEOS model currently includes two PBL height definitions suitable for direct comparison with observed PBL height, and additional definitions are being added in this study. Evaluation of different model PBL height definitions is underway. Meanwhile, efforts have been made in the GEOS data assimilation system to develop PBL height data assimilation capability. PBL height data can be assimilated using two different approaches. The traditional approach is to construct an observation operator and its tangent linear and adjoint, which link control variables to PBL height data from each observing system. This observation operator can be very complicated, e.g., the lidar-based PBL height observation operator includes the backscatter lidar forward observation operator, the algorithm to derive PBL height from attenuated total backscatter, interpolation, and calculations handling the mismatch between observed and model scales. The other approach is to augment PBL height to the control variable vector, and it is adopted in this study. The latter approach was also used in previous studies, e.g., the assimilation of PBL height data from radiosonde and aircraft in the Real Time Mesoscale Analysis (RTMA) system for a dispersion modelling study ([13]); the PBL height assimilation study using lidar PBL height data at Greensburg, Kansas for a field campaign ([12]). The PBL height assimilation from multiple observing systems in this study allows us to take advantage of the diverse PBL height data that provide much better global coverage collectively under different meteorological conditions and with different temporal and spatial scales. As all the PBL heights are tightly coupled with the PBL thermodynamic variables, the strong correlations, which are provided by the 4D ensemble forecast, enable PBL height data from various sources to interact and combine coherently and provide additional information for PBL temperature and moisture fields. The results of comparisons among PBL height data from different sources and the evaluation of the model PBL height definitions with the PBL height data will be presented, and the PBL height data synergy strategies and preliminary results will also be discussed at the conference.

Y. Zhu↗
NASA NTRS2002

Using Clustering to Establish Climate Regimes from PCM Output

A multivariate statistical clustering technique--based on the k-means algorithm of Hartigan has been used to extract patterns of climatological significance from 200 years of general circulation model (GCM) output. Originally developed and implemented on a Beowulf-style parallel computer constructed by Hoffman and Hargrove from surplus commodity desktop PCs, the high performance parallel clustering algorithm was previously applied to the derivation of ecoregions from map stacks of 9 and 25 geophysical conditions or variables for the conterminous U.S. at a resolution of 1 sq km. Now applied both across space and through time, the clustering technique yields temporally-varying climate regimes predicted by transient runs of the Parallel Climate Model (PCM). Using a business-as-usual (BAU) scenario and clustering four fields of significance to the global water cycle (surface temperature, precipitation, soil moisture, and snow depth) from 1871 through 2098, the authors' analysis shows an increase in spatial area occupied by the cluster or climate regime which typifies desert regions (i.e., an increase in desertification) and a decrease in the spatial area occupied by the climate regime typifying winter-time high latitude perma-frost regions. The patterns of cluster changes have been analyzed to understand the predicted variability in the water cycle on global and continental scales. In addition, representative climate regimes were determined by taking three 10-year averages of the fields 100 years apart for northern hemisphere winter (December, January, and February) and summer (June, July, and August). The result is global maps of typical seasonal climate regimes for 100 years in the past, for the present, and for 100 years into the future. Using three-dimensional data or phase space representations of these climate regimes (i.e., the cluster centroids), the authors demonstrate the portion of this phase space occupied by the land surface at all points in space and time. Any single spot on the globe will exist in one of these climate regimes at any single point in time. By incrementing time, that same spot will trace out a trajectory or orbit between and among these climate regimes (or atmospheric states) in phase (or state) space. When a geographic region enters a state it never previously visited, a climatic change is said to have occurred. Tracing out the entire trajectory of a single spot on the globe yields a 'manifold' in state space representing the shape of its predicted climate occupancy. This sort of analysis enables a researcher to more easily grasp the multivariate behavior of the climate system.

Oglesby, Robert↗
NASA NTRS2003

Control of Interacting Vortex Flows at Subsonic and Transonic Speeds Using Passive Porosity

A wind tunnel experiment was conducted in the NASA Langley Research Center (LaRC) 8-foot Transonic Pressure Tunnel (TPT) to determine the effects of passive surface porosity on vortex flow interactions about a general research fighter configuration at subsonic and transonic speeds. Flow- through porosity was applied to a wind leading-edge extension (LEX) mounted to a 65 deg cropped delta wind model to promote large nose-down pitching moment increments at high angles of attack. Porosity decreased the vorticity shed from the LEX, which weakened the LEX vortex and altered the global interactions of the LEX and wing vortices at high angles of attack. Six-component forces and moments and wing upper surface static pressure distributions were obtained at free- stream Mach numbers of 0.50, 0.85, and 1.20, Reynolds number of 2.5(10(exp-6) per foot, angles of attack up to 30 deg and angles of sideslip to plus or minus 8 deg. The off-surface flow field was visualized in selected cross-planes using a laser vapor screen flow visualization technique. Test data were obtained with a centerline vertical tail and with alternate twin, wing-mounted vertical fins having 0 deg and 30 deg cant angles. In addition, the porosity of the LEX was compartmentalized to determine the sensitivity of the vortex- dominated aerodynamics to the location and level of porosity applied to the LEX.

NATO FURNISHED↗
NASA NTRS1999

Optical Flow-Field Techniques Used for Measurements in High-Speed Centrifugal Compressors

The overall performance of a centrifugal compressor depends on the performance of the impeller and diffuser as well as on the interactions occurring between these components. Accurate measurements of the flow fields in each component are needed to develop computational models that can be used in compressor design codes. These measurements must be made simultaneously over an area that covers both components so that researchers can understand the interactions occurring between the two components. Optical measurement techniques are being used at the NASA Lewis Research Center to measure the velocity fields present in both the impeller and diffuser of a 4:1 pressure ratio centrifugal compressor operating at several conditions ranging from design flow to surge. Laser Doppler Velocimetry (LDV) was used to measure the intrablade flows present in the impeller, and the results were compared with analyses obtained from two three-dimensional viscous codes. The development of a region of low throughflow velocity fluid within this high-speed impeller was examined and compared with a similar region first observed in a large low-speed centrifugal impeller at Lewis. Particle Image Velocimetry (PIV) is a relatively new technique that has been applied to measuring the diffuser flow fields. PIV can collect data rapidly in the diffuser while avoiding the light-reflection problems that are often encountered when LDV is used. The Particle Image Velocimeter employs a sheet of pulsed laser light that is introduced into the diffuser in a quasi-radial direction through an optical probe inserted near the diffuser discharge. The light sheet is positioned such that its centerline is parallel to the hub and shroud surfaces and such that it is parallel to the diffuser vane, thereby avoiding reflections from the solid surfaces. Seed particles small enough to follow the diffuser flow are introduced into the compressor at an upstream location. A high-speed charge-coupled discharge (CCD) camera is synchronized to the laser pulse rate; this allows it to capture images of seed particle position that are separated by a small increment in time. A crosscorrelation of a particle's position in two consecutive images provides an estimate of flow velocity and direction. Multiple image pairs obtained in rapid succession at a particular flow condition provide enough measurements for statistical significance. PIV provides simultaneous velocity measurements over the entire plane that is illuminated by the light sheet instead of at a single point, as is the case when LDV is used. PIV has a further advantage in that the laser light pulse can be triggered by an external source such as a high-response pressure transducer. This feature will allow PIV to synchronize flow imaging to physical phenomena such as rotating stall or stall precursor waves. We hope that this technique can be used to obtain images of the flow field during and just prior to stall.

Skoch, Gary J.↗
NASA NTRS2016

Serious Gaming for Test & Evaluation of Clean-Slate (Ab Initio) National Airspace System (NAS) Designs

Incremental approaches to air transportation system development inherit current architectural constraints, which, in turn, place hard bounds on system capacity, efficiency of performance, and complexity. To enable airspace operations of the future, a clean-slate (ab initio) airspace design(s) must be considered. This ab initio National Airspace System (NAS) must be capable of accommodating increased traffic density, a broader diversity of aircraft, and on-demand mobility. System and subsystem designs should scale to accommodate the inevitable demand for airspace services that include large numbers of autonomous Unmanned Aerial Vehicles and a paradigm shift in general aviation (e.g., personal air vehicles) in addition to more traditional aerial vehicles such as commercial jetliners and weather balloons. The complex and adaptive nature of ab initio designs for the future NAS requires new approaches to validation, adding a significant physical experimentation component to analytical and simulation tools. In addition to software modeling and simulation, the ability to exercise system solutions in a flight environment will be an essential aspect of validation. The NASA Langley Research Center (LaRC) Autonomy Incubator seeks to develop a flight simulation infrastructure for ab initio modeling and simulation that assumes no specific NAS architecture and models vehicle-to-vehicle behavior to examine interactions and emergent behaviors among hundreds of intelligent aerial agents exhibiting collaborative, cooperative, coordinative, selfish, and malicious behaviors. The air transportation system of the future will be a complex adaptive system (CAS) characterized by complex and sometimes unpredictable (or unpredicted) behaviors that result from temporal and spatial interactions among large numbers of participants. A CAS not only evolves with a changing environment and adapts to it, it is closely coupled to all systems that constitute the environment. Thus, the ecosystem that contains the system and other systems evolves with the CAS as well. The effects of the emerging adaptation and co-evolution are difficult to capture with only combined mathematical and computational experimentation. Therefore, an ab initio flight simulation environment must accommodate individual vehicles, groups of self-organizing vehicles, and large-scale infrastructure behavior. Inspired by Massively Multiplayer Online Role Playing Games (MMORPG) and Serious Gaming, the proposed ab initio simulation environment is similar to online gaming environments in which player participants interact with each other, affect their environment, and expect the simulation to persist and change regardless of any individual player's active participation.

Allen, B. Danette↗
NASA NTRS2022

22” ADP Fan Rig Liners Design Report

A liner design study was conducted as part of a cooperative effort between five (5) government/industry teams that together seek to demonstrate the technology for designing and manufacturing acoustic liners that are twenty five percent (25%) more efficient than 1992 technology liners. The study emphasized teaming collaboration. The improved liners were designed by Pratt & Whitney and Boeing Airplane Company and built by Rohr Inc, and will be tested in the NASA/P& W 22-inch ADP fan rig at NASA Lewis Research Center's 9' x 15' wind tunnel. Design guidelines and decisions were made collectively during monthly design review telecons and at the formal final design review. The tools that were used to design the new improved liners were not new, but the process that was developed as part of this study that led into the evaluation and selection of the final and best liner designs is new. Until now, the liner design process as practiced by different industry teams varied greatly. Some procedures were based on empirical liner attenuation databases which could not adequately account for engine-to-engine hardwall fan noise spectral differences, while others were entirely theoretical. And regardless of which method one chose, the major difficulty was still the lack of understanding and knowledge of the actual hardwall fan source noise modal structure and farfield SPL spectra. The NASA-led effort to attempt actual measurements of the fan source noise modal structure by means of the rotating microphone array is expected to contribute significantly to the understanding of the nature of the fan tone noise modes, but the application to broadband noise is still a long way off. The new design process included consideration of the measured fan tone modes, but for the majority of the spectra a separate, systematic process was used. It is a common practice in the engine/nacelle industry that acoustic liners for new products are designed before actual measured hardwall engine far-field noise spectra are available. Target noise spectra are normally derived from existing engine noise databases with some adjustments to absolute SPLs and frequencies. This practice has been acceptable as long as the new engine was a derivative of the base engine. However, in the case of the ADP, the transition from a current engine base is too - great, and the adjustments could not adequately account for the quantum changes in SPL spectral differences. Examples were the 1992 single degree of freedom (SDOF) inlet and aft liners (designed by P&W) that would be used as the baseline liners against which the new improved liners' efficien¬cies would -be measured. As will be shown, these baseline liners have been found to be deeper than the desire optimum depths. The new liner design process begins with the selection of the target hard wall fan noise spectra. These spectra were obtained by scaling up (5.91 scale factor) the measured hardwall fan spectral data from the 22” ADP fan rig. Next, the modal energy contents for each 1/3-octave band center frequency of these hardwall fan noise spectra are estimated. (Within each 1/3-octave band center frequency, the model energy distribution approximation is for both tone and broadband noise). For the inlet noise, P&W uses a derivative of the NASA Lewis (Ed Rice's) method of classifying; and grouping propagating modes by their cutoff ratios. The next step is to assign energy level to each group of modes having the same cutoff ratio. In Ed Rice's model, the modes were grouped according to ten (10) cutoff ratio intervals with center values located at 1.026, 1.085, 1.155, 1.24, 1.35, 1.49, 1.69, 2.0, and 4.47 (with equal number of modes in each interval). The modes were assumed to have equal energy. P&W’s model expands the cutoff ratio-mode grouping into two hundred (200) smaller cutoff ratio intervals with center values located f1".'m 1.003 to 11.5 in 199 increasing incremental intervals. Next, P&W's model uses a "2-parameter" normal distribution as a template to assign energy levels to these 200 pre-determined cutoff ratio values (for each 1/3-octave band center frequency). In the past, P&W had conducted an extensive study to determine what "2-parameter" values are appropriate for fullscale inlet liners, and had developed a set of twenty-four (24) "2-pararneter" values (i.e. one for each 1/3-octave frequency band) that when used in Ed Rice's Inlet Attenuation Prediction Method produced predicted liner attenuations that closely matched measured liner attenuations from several P&W's engines. In the absence of actual measured tone and broadband modal data from the 22" ADP rig, the process will use P&W's proprietary set of "2-parameter" values for this liner design study. For the aft noise, Boeing uses a modal energy approximation that the "transport energy" of each propagating mode is equal. This approximation is almost the same as the equal energy per mode approximation, except for modes that are near cutoff. Boeing's model forces these modes to have lower energy levels. Both P&W and Boeing agreed that the "transport energy" approximation should work well for the ADP aft fan noise which appears to be dominated by broadband noise. The design process then proceeds to calculate the optimum liner impedances for each frequency in both the inlet and the aft. These optimum impedances represent the target impedances that the designed liners should have. Liners with impedances matching the optimum impedances at all frequencies are "ideal" liners. These ideal liners are theoretically the best liners. Unfortunately, it has been showed that it is impossible to design and build such ideal liners. The next best liners are ones that have impedances matching the optimum impedances over some frequency range (not all frequencies as for the ideal). This is accomplished by the use of "frequency weightings". Several of P&W's and Boeing's existing computer decks were used for optimizing and matching the designed liner impedances to the target optimum values (with the various frequency weightings specified). The optimization produces liner candidates with predicted liner impedances and descriptions of their physical liner characteristics. These candidate liner impedances are then used to predict their spectral attenuation characteristics which are then used together with the target hardwall fan noise spectra to determine the resulting treated noise spectra and PNLT values. Further optimization around the selected candidate designs yield final designs that are best in PNLT attenuations. Use of the optimization decks allowed a large number of liner candidates to be screened in a relatively short period of time. This design process is systematic and is efficient. The new inlet SDOF liner design obtained from the improved process was predicted to be 34% more effective (per unit area) than the 1992 baseline inlet liner. This inlet liner design is a 112 rayl wovenwiremesh facesheet over a 0.312-inch deep honeycomb core. The new aft SDOF liner design was predicted to be 52% more efficient than the aft baseline liner. The new aft SDOF liner is "segmented" with a shallower liner on the core cowl (inner duct wall) and a deeper liner on the fan cowl (outer duct wall). The shallower liner is a 70.6 rayl woven-wiremesh facesheet over a 0.141-inch deep honeycomb core. The deeper liner is a 68 rayl woven-wiremesh facesheet over a 0.309-inch deep honeycomb core. All liner dimensions are for the 22-inch model-scale ADP fan rig liners. The selected advanced liners are "segmented" double-layer (DDOF) liners for the inlet and aft locations. Also, for the inlet, a bulk liner with ceramic foam for wider broadband noise absorption was also selected. Triple-layer liner designs were not considered since the model scaled liners ( 1/5. 91) were dimensionally too small to be built correctly, and irrin earlier concept study, Boeing found a triple-layer to have only very small benefits over a double-layer. The inlet DDOF liner was predicted to be 83% more effective than the baseline. The inlet DDOF design is a 78 rayl facesheet over a 0.080 top cavity depth, a 68 rayl septum and a 0.227-inch bottom cavity depth. The inlet bulk liner is a 60 rayl facesheet over a 0.33-inch deep honeycomb filled with high temperature (HTP) ceramic foam with a density of 4.8 lb/cu.ft and a flow resistivity of 167 rayl/cm. The inlet bulk liner was predicted to be 83% more effective than the baseline liner. The aft DDOF liners are segmented with a shallower liner on the core cowl and a deeper liner on the fan cowl. The shallow DDOF liner is a 49.8 rayl facesheet over a 0.093-inch top cavity depth, a septum of 88.2 rayls over a 0.181-inch bottom cavity depth. The deep DDOF liner is a 12.9 rayl facesheet over a 0.140-inch top cavity depth, a septum of 53.1 ray ls over a 0.258-inch botom cavity depth. The segmented aft DDOF liners were predicted to be 86% more efficient than the baseline liner.

turbofan acoustic treatment↗
NASA NTRS2021

NASA Electric Vertical Takeoff and Landing (eVTOL) Aircraft Technology for Public Services – A White Paper

History has shown that our personal life is highly dependent on the technology that people have developed. A strategic scan of the aerospace environment at the beginning of the 21st century strongly suggests that the world might be approaching a new age of airpower—the era of electrified/hybrid aircraft propulsion. Undeniably, starting from the Montgolfier Brothers balloon flight in 1783, to the Wright Brothers piston engine flight in 1903, and the jet engine of the 1960s, or the space age of today, one can say that leaps in propulsion technology have marked the different ages of human flight. The technological advancements, brought at the beginning of 21st century by the revolution in data exchange, computational power, sensors, wireless communication, internet, and autonomy, contributed to the vision of this new age of propulsion we are approaching. Historically, conventional vertical takeoff and landing (VTOL) aircraft have been equipped with propulsion units relying on complex internal combustion machines (turbines, piston engines, for example), and complex mechanical arrangements (gearboxes, shafts, variable pitch propeller). By contrast, electric VTOL aircraft (eVTOL)1 rely on simpler propulsion units (electric motors and in some cases fixed-pitch propellers). This promotes redundancy and improves tolerance to failures, in turn improving safety. The use of simpler electric propulsion units should also allow significant acquisition and operating cost reductions. Whether full-electric (relying solely on batteries) or hybrid-electric (relying on a combination of batteries, fuel-powered engines, and generators.), eVTOLs are also expected to generate less noise and air pollution than conventional aircraft with similar payloads. According to the 2019 Annual Review of IATA (International Air Transport Association) [ref.1], due to an expected increase in air transport traffic by 5% every year and a doubling of air transport passenger numbers to 8.2 billion by 2037 significant challenges are posed to the aviation industry. Furthermore, this report does not factor in the expected demand for short-range (intra-city) air transportation, which is in development and yet to be operational. The increased demand to fly creates a responsibility to expand in a sustainable manner and an endeavor to develop more environmentally-friendly aircraft. eVTOL aircraft, either piloted or autonomous, is gathering considerable interest worldwide. Modern and novel full-electric or hybrid-electric eVTOL configurations enable a new paradigm shift in air transportation as the aviation industry remains committed to its goals of carbon-neutral growth from 2020 onwards and cutting CO2 emissions to half 2005 levels by 2050. While electric power has been used for decades, recent developments in mobile electric/hybrid propulsion coupled with advanced materials and autonomous systems may create the possibility to transition into the next age of air mobility propelled by electric/hybrid VTOL aircraft technology. Although eVTOL aircraft might seem like an incremental improvement or even a counterintuitive regression with regard to past VTOL development, it has in fact the potential to transform air mobility across a wide range of government applications. Previous transformations in aviation generated dramatic leaps in performance, but the cost was commensurate with performance, limiting quantity produced. This next age appears to take a different approach. Performance may not increase, but at this moment technology is poised for future urban mobility that will spawn commercial passenger drone services, that is, autonomous (pilotless) air taxis and thereby add a new dimension to the urban transportation mix of the future [ref. 2]. Advances in electric propulsion, autonomous flight technology, and 5G communication networks will enable this fast new-growing market to become a reality. It is now time to envision the introduction of electric/hybrid eVTOL aircraft for Public Services2. We believe that in the next decades eVTOL aircraft will have the potential to become an essential tool to Public Service agencies around the world in applications such as firefighting, public safety, search and rescue, disaster relief and law enforcement. This is due to several major factors. • First, with the increasing popularity of small, unmanned aircraft vehicles (UAVs) or drones, many companies today are focusing on the development of passenger UAVs designed to accommodate up to five passengers or equivalent cargo payload. Many such configurations are electric or hybrid-electric designs with VTOL capabilities. Several of these projects have started a flight test program and many more are expected to be in the experimental and development phase in 2020. Such revolutionary vehicles could be in commercial operations by 2030. These eVTOL systems could be ready for selected Public Services missions even sooner. • Second, although these advanced eVTOL vehicles under development still need access to fuel (hybrid) and/or electric charging capability, they can take off and land from almost anywhere. Therefore, such vehicles, both manned and unmanned can be successfully integrated for the critical missions of the Public Services with extra deployment flexibilities. • Third, advancement in electric propulsion systems in the automotive industry together with NASA’s leading efforts in electrification of aircraft propulsion systems, FAA’s ongoing active eVTOL certification programs, and EASA’s proposed framework for the certification of electric/hybrid small category VTOL aircraft in Europe [ref. 3] will help accelerate industry electric propulsion system development and integration. • Finally, eVTOL vehicles could be deployed for Public Services sooner than air taxi or other commercial applications, since Public Services missions may be more easily approved based on specific mission criteria, localized airworthiness authority for public-use aircraft3, and are normally operating under centralized airspace management and control by the theater command. Moreover, public perception and acceptance are generally less of a concern when operations save lives and benefit the wider community. The prioritized introduction of eVTOL aircraft in Public Services is ambitious, but we believe it is achievable in the coming decades if fundamental enablers (people and technologies) are engaged in defining the objectives and needs of these missions. The revolution that is currently taking place in eVTOL aircraft represents an unprecedented opportunity to develop a safer, more affordable, more available and more environmentally friendly future of vertical flight. To ensure that these novel aircraft meet the future expectations of Public Services, it is essential to take a collaborative and multi-disciplinary approach to their development, across engineering disciplines, policy-making, program management, business case development, manufacturing, and flight demonstrations. It should be noted that the term eVTOL (in the near term) used throughout this publication implies aircraft capable of transporting up to 5 persons which may or may not include a pilot if operated fully autonomously, assuming an average of 200 pounds (91 kg) per person or equivalent payload and a range up to 60 miles plus a suitable reserve. Hybrid or hydrogen powered eVTOLs would have greater range. For example, a “3-seat” eVTOL aircraft may only be able to carry two fully equipped firemen, and payload capacity is more relevant when used for the supply mission. Moreover, this paper concentrates on the “last-mile” solutions with a deployment time of no more than 6 hours. Although not specifically discussed in this document, it is understood that the future of Transformative Vertical Flight in general and Public Services, in particular, will also involve smaller UAVs that will undoubtedly play a crucial role in future aerial operations. For example, smaller unmanned aircraft may be used to dispatch medical supplies, portable filtration systems or perform the Search task of future Search and Rescue (SAR) operations. Close collaboration between the aircraft industry, the Civil Aviation Authorities (CAA), e.g., Federal Aviation Administration (FAA), European Aviation Safety Agency (EASA), Transport Canada Civil Aviation (TCCA) and the Department of Defense (DoD) certifiers, will help identify Public Services requirements, define expectations and limit development cost and timescales. Take the US Air Force Agility Prime as an example, the majority of the eVTOL application opportunities and mission elements identified are in line with the NASA TVF WG-4 objectives and use cases. Together, it forms a strong partnership to accelerate the development, certification, and practical deployment for public service missions. The US Air Force Agility Prime has been a collaboration partner on this white paper, and provided valuable input and recommendations. Most of the eVTOL public service mission elements discussed in this paper and additional use cases envisioned by the NASA TVF WG-4 team are shared by the Agility Prime program. The focus and efforts of the Agility Prime in product and system development, industry and government partnership, accelerated certifications as well as early test and deployment are totally in sync with the path forward recommended by this white paper. This kind of collaboration and partnership will help enable the practical use of the eVTOL for public service missions, benefit the eVTOL public acceptance, and accelerate the eVTOL industry revolution.

Johnny T. Doo↗
NASA NTRS2005

Analysis of Piezoelectric Structural Sensors with Emergent Computing Techniques

The purpose of this project was to try to interpret the results of some tests that were performed earlier this year and to demonstrate a possible use of emergence in computing to solve IVHM problems. The test data used was collected with piezoelectric sensors to detect mechanical changes in structures. This project team was included of Dr. Doug Ramers and Dr. Abdul Jallob of the Summer Faculty Fellowship Program, Arnaldo Colon-Lopez - a student intern from the University of Puerto Rico of Turabo, and John Lassister and Bob Engberg of the Structural and Dynamics Test Group. The tests were performed by Bob Engberg to compare the performance two types of piezoelectric (piezo) sensors, Pb(Zr(sub 1-1)Ti(sub x))O3, which we will label PZT, and Pb(Zn(sub 1/3)Nb(sub 2/3))O3-PbTiO, which we will label SCP. The tests were conducted under varying temperature and pressure conditions. One set of tests was done by varying water pressure inside an aluminum liner covered with carbon-fiber composite layers (a cylindrical "bottle" with domed ends) and the other by varying temperatures down to cryogenic levels on some specially prepared composite panels. This report discusses the data from the pressure study. The study of the temperature results was not completed in time for this report. The particular sensing done with these piezo sensors is accomplished by the sensor generating an controlled vibration that is transmitted into the structure to which the sensor is attached, and the same sensor then responding to the induced vibration of the structure. There is a relationship between the mechanical impedance of the structure and the resulting electrical impedance produced in the in the piezo sensor. The impedance is also a function of the excitation frequency. Changes in the real part of impendance signature relative to an original reference signature indicate a change in the coupled structure that could be the results of damage or strain. The water pressure tests were conducted by pressurizing the bottle on a test stand, and running sweeps of excitations frequencies for each of the piezo sensors and recording the resulting impedance. The sweeps were limited to 401 points by the available analyzer, and it was decided to perform individual sweeps at five different excitation frequency ranges. The frequency ranges used for the PZTs were different in two of the five ranges from the ranges used for the SCP. The bottles were pressurized to empty (no water), 0psig, 77 psig, 155 psig, 227 psig in nearly uniform increments of about 77psi. One of each of the two types of piezo sensors was fastened on to the bottle surface at two locations: about midway between the ends on cylindrical portion of the bottle and at the very edge of one of the end domes. The data was collected in files by sensor type (2 cases), by location (2 cases), by frequency range (5 cases), and pressure (5cases) to produce 100 data sets of 401 impedances. After familiarization with the piezo sensing technology and obtaining the data, the team developed a set of questions to try to answer regarding the data and made assignments of responsibilities. The next section lists the questions, and the remainder of the report describes the data analysis work performed by Dr. Ramers. This includes a discussion of the data, the approach to answering the question using statistical techniques, the use of an emergent system to investigate the data where statistical techniques were not usable, conclusions regarding the data, and recommendations.

Ramers, Douglas L.↗