Search NASA⌕ Search

SEARCH · Search NASA

Results for “methods: data analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,675 records · Page 93

Single-Mode Projection Filters for Modal Parameter Identification for Flexible Structures

Single-mode projection filters are developed for eigensystem parameter identification from both analytical results and test data. Explicit formulations of these projection filters are derived using the orthogonal matrices of the controllability and observability matrices in the general sense. A global minimum optimization algorithm is applied to update the filter parameters by using the interval analysis method. The updated modal parameters represent the characteristics of the test data. For illustration of this new approach, a numerical simulation for the MAST beam structure is shown by using a one-dimensional global optimization algorithm to identify modal frequencies and damping. Another numerical simulation of a ten-mode structure is also presented by using a two-dimensional global optimization algorithm to illustrate the feasibility of the new method. The projection filters are practical for parallel processing implementation.

Huang, Jen-Kuang↗

Single-mode projection filters for modal parameter identification for flexible structures

Single-mode projection filters are developed for eigensystem parameter identification from both analytical results and test data. Explicit formulations of these projection filters are derived using the orthogonal matrices of the controllability and observability matrices in the general sense. A global minimum optimization algorithm is applied to update the filter parameters by using the interval analysis method. The updated modal parameters represent the characteristics of the test data. For illustration of this new approach, a numerical simulation for the MAST beam structure is shown by using a one-dimensional global optimization algorithm to identify modal frequencies and damping. The projection filters are practical for parallel processing implementation.

Huang, Jen-Kuang↗

A Comparison of Methods for Assessing Space Suit Joint Ranges of Motion

Through the Advanced Exploration Systems (AES) Program, NASA is attempting to use the vast collection of space suit mobility data from 50 years worth of space suit testing to build predictive analysis tools to aid in early architecture decisions for future missions and exploration programs. However, the design engineers must first understand if and how data generated by different methodologies can be compared directly and used in an essentially interchangeable manner. To address this question, the isolated joint range of motion data from two different test series were compared. Both data sets were generated from participants wearing the Mark III Space Suit Technology Demonstrator (MK-III), Waist Entry I-suit (WEI), and minimal clothing. Additionally the two tests shared a common test subject that allowed for within subject comparisons of the methods that greatly reduced the number of variables in play. The tests varied in their methodologies: the Space Suit Comparative Technologies Evaluation used 2-D photogrammetry to analyze isolated ranges of motion while the Constellation space suit benchmarking and requirements development used 3-D motion capture to evaluate both isolated and functional joint ranges of motion. The isolated data from both test series were compared graphically, as percent differences, and by simple statistical analysis. The results indicated that while the methods generate results that are statistically the same (significance level p= 0.01), the differences are significant enough in the practical sense to make direct comparisons ill advised. The concluding recommendations propose direction for how to bridge the data gaps and address future mobility data collection to allow for backward compatibility.

Aitchison, Lindsay T.↗

Tip vortex geometry of a hovering helicopter rotor in ground effect

The wide-field shadowgraph method has been used to photograph the tip vortices of a hovering helicopter rotor in ground effect. The shadowgraphs were used to obtain quantitative measurements of the rotor tip vortex geometry both in and out of ground effect. Many important phenomena are visible in the rotor wake using this method. These include the variation in descent and contraction rates of the tip vortices in ground effect, and the interaction between tip vortices in the far wake. The tip vortex geometry from the shadowgraphs is compared with the tip vortex geometry predicted using a free wake hover performance analysis. The free wake analysis accurately predicts the tip vortex geometry both in and out of ground effect. Performance data from the test is compared with the performance predicted using several methods, including the free wake analysis. All methods provided reasonable predictions for the helicopter performance in ground effect.

Light, Jeffrey S.↗

Structural and Stratigraphic Controls on Methane Hydrate occurrence and distribution: Gulf of Mexico, Walker Ridge 313 and Green Canyon 955: Final Report

The goal of this project was to determine structural and stratigraphic controls on hydrate occurrence and distribution in Green Canyon (GC) 955 and Walker Ridge (WR) 313 blocks using seismic and well data. Gas hydrate was discovered in these blocks in coarse- and fine-grained sediments during the 2009 Joint Industrial project (JIP) Leg 11 drilling expedition. Although the immediate interest of the exploration community is exclusively hydrate which is present in coarse–grained sediments, factors that control hydrate and free gas distribution in the two blocks and whether coarse and fine-grained hydrate-bearing units are related in any manner, formed the core of this research. The project spanned from 10/01/2012 to 07/31/2016. In the project, in both the leased blocks, the interval spanning the gas hydrate stability zone (GHSZ) was characterized using a joint analysis of sparse Ocean Bottom Seismic (OBS) and dense, surface–towed multichannel seismic (MCS) data. The project team had the luxury of calibrating their results with two well logs. Advance processing methods such as depth migration and full-waveform inversion (FWI) were used for seismic data analysis. Hydrate quantification was achieved through interpretation of the FWI velocity field using appropriate rock physics models at both blocks. The seismic modeling/inversion methodology (common to both GC955 and WR313 blocks) was as follows. First, the MCS data were depth migrated using a P-wave velocity (VP) model constructed using inversion of reflection arrival times of a few (four in both cases) key horizons carefully picked in the OBS data to farthest possible offsets. Then, the resolution of the traveltime VP model was improved to wavelength scale by inverting OBS gathers up to the highest frequency possible (21.75 Hz for GC955 and 17.5 for WR313) using FWI. Finally, the hydrate saturation (or the volume fraction) was estimated at the well location assuming one of the other hydrate morphology (filling the primary or the secondary porosity) was extrapolated out from the wells using the FWI VP as a guide. General outcomes were as follows. First and foremost, an imaging methodology using sparse seismic data, which is easily replicable at other sites with similar datasets, has been demonstrated. The end product of this methodology at both the leased blocks is quantitative estimates of hydrate distribution. Second, at both locations there is strong evidence that the base of the GHSZ, which does not appear as a clear Bottom Simulating Reflection (BSR), manifests in the VP perturbations created by FWI, suggesting that FWI is sensitive to subtle compositional changes in shallow sediments and establishes it as a valuable tool for investigations of hydrate-bearing basins. Third, through joint interpretation of the depth migrated image and the FWI VP model, how structure and stratigraphy jointly determine hydrate and free gas distribution in both blocks could be clearly visualized. The joint interpretation also suggests that the coarse and fine grained hydrate-bearing sediments at both leased are connected. Site specific results, in addition to general results, are as follows. At GC955 the overlying fine-grained hydrate-bearing unit could have been sourced from the underlying hydrate coarse-grained channel-levee complex through a chimney feature. The channel-levee system at GC955 is compartmentalized by faults, of which only a few may be impermeable. Although compartmentalized, the channel-levee system in the GC955 as a whole might be in communication except selected zones. At WR313 the overlying fine-grained fracture-filled hydrate unit appears to be sourced from below the GHSZ. The reason that only a particular fine-grained unit has hydrate, despite having lower porosity that the bounding units, could be the presence of secondary porosity (such as those formed from clay dewatering under compaction). In conclusion, the project was a pioneering effort in in joint analysis of OBS and MCS datasets for advancing the knowledge about a hydrate and free–gas system dynamics using advanced processing methods such as FWI and depth migration. Results obtained in this project can greatly advance the tools and techniques used for delineating specific hydrate prospects. Results obtained in this project can also be seamlessly incorporated into other DOE funded project on modeling the potential productivity and commercial viability of hydrate from sand-dominated reservoirs. The OBS and MCS data in this project were acquired in 2012 (after the JIP II drilling) by the USGS and therefore the results are a posteriori. Nonetheless, the seismic inversion workflow established through this project can be used to generate various what-if quantification scenarios even in absence of logs and serve as a valuable tool for guiding drilling operations. Results from this project can augment other DOE sponsored projects on determining the commercial viability of methane production from the Gulf of Mexico.

03 NATURAL GAS↗

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗

Reduction and Analysis of Phosphor Thermography Data With the IHEAT Software Package

Detailed aeroheating information is critical to the successful design of a thermal protection system (TPS) for an aerospace vehicle. This report describes NASA Langley Research Center's (LaRC) two-color relative-intensity phosphor thermography method and the IHEAT software package which is used for the efficient data reduction and analysis of the phosphor image data. Development of theory is provided for a new weighted two-color relative-intensity fluorescence theory for quantitatively determining surface temperatures on hypersonic wind tunnel models; an improved application of the one-dimensional conduction theory for use in determining global heating mappings; and extrapolation of wind tunnel data to flight surface temperatures. The phosphor methodology at LaRC is presented including descriptions of phosphor model fabrication, test facilities and phosphor video acquisition systems. A discussion of the calibration procedures, data reduction and data analysis is given. Estimates of the total uncertainties (with a 95% confidence level) associated with the phosphor technique are shown to be approximately 8 to 10 percent in the Langley's 31-Inch Mach 10 Tunnel and 7 to 10 percent in the 20-Inch Mach 6 Tunnel. A comparison with thin-film measurements using two-inch radius hemispheres shows the phosphor data to be within 7 percent of thin-film measurements and to agree even better with predictions via a LATCH computational fluid dynamics solution (CFD). Good agreement between phosphor data and LAURA CFD computations on the forebody of a vertical takeoff/vertical lander configuration at four angles of attack is also shown. In addition, a comparison is given between Mach 6 phosphor data and laminar and turbulent solutions generated using the LAURA, GASP and LATCH CFD codes. Finally, the extrapolation method developed in this report is applied to the X-34 configuration with good agreement between the phosphor extrapolation and LAURA flight surface temperature predictions. The phosphor process outlined in the paper is believed to provide the aerothermodynamic community with a valuable capability for rapidly obtaining (4 to 5 weeks) detailed heating information needed in TPS design.

Merski, N. Ronald↗

THE INFLUENCE OF CHANGES IN NUCLEAR COVARIANCE DATA ON THE CALCULATION OF Ck FOR HIGHLY ENRICHED URANIUM SOLUTION SYSTEMS

Since their initial development nearly 20 years ago, sensitivity/uncertainty (S/U) techniques have been increasingly applied to criticality safety validation. These techniques can generally be applied quickly and for a range of purposes after the initial investment of calculating sensitivities for a safety analysis model and a range of potential benchmark experiments. This paper discusses these applications and how they have been affected by changes in covariance data libraries.Critical experiment selection is the most prominent area in which S/U methods impact validation efforts. These methods allow for rigorous calculation of similarity between safety analysis models and potential benchmarks. The most widely used of these methods is the integral parameter ck, although other integral parameters such as E and some non-integral parameters have also been explored. The ck parameter calculates the fraction of nuclear data–induced uncertainty shared between an application and an experiment; this metric is predicated on the assessment that nuclear data errors are the most likely source of computational bias. Therefore, values of ck are inherently dependent on the nuclear covariance library used to calculate them and inevitably change when these data change.Past studies have shown how much the ck values can vary, especially in difficult validation cases like burnup credit. These significant changes are the result of differences in the relative uncertainties of nuclides with large sensitivities, and they imply differences in benchmark applicability to the safety analysis model. This paper demonstrates that more common criticality safety models involving single fissile species have significantly less variability of ck between covariance libraries, thus enhancing confidence in the application of S/U methods for experiment selection in these scenarios. This paper presents case studies examining the impact of these changes on some hypothetical safety analysis systems.

Clarity, Justin↗

Hydrogen effects on thermal diffusivity and electrical resistivity of zircaloy cladding

Thermal diffusivity measurements on Zirconium-based cladding materials have historically been a challenge due to the difficulty to measure on specimens with curved geometries, including nuclear grade Zircaloy cladding materials. Here, in this work, we first used laser flash analysis method and four-probe configuration method to measure the thermal diffusivity and electrical resistivity of Zircaloy tubes respectively, which show good agreement with Zircaloy plates in this work, as well as previously published data. The consistent results proved the applicability of the laser flash analysis setup and four-probe configuration method for investigating thermal diffusivity and electrical resistivity of Zircaloy tubes. We further investigated the hydrogen effect on thermal diffusivity and electrical resistivity of Zircaloy. Hydrogen plays significant roles in the thermal diffusivity of Zircaloy, which depends on the hydrogen concentration. For higher hydrogen contents (1130 and 1820 wppm in this work), where phase transformation (α-Zr + δ-hydride → α-Zr + β-Zr) occurs at 567 °C, thermal diffusivity decreases as a function of temperature at the α-Zr + δ-hydride phase regime, while increase at higher temperature at the α-Zr + β-Zr phase regime. For low hydrogen concentration, where hydride dissolved into α-Zr matrix phase at higher temperature, the thermal diffusivity is lower than non-hydrided Zircaloy-4 with similar temperature-dependent trends. Such observation demonstrates the hydrogen effects on reducing thermal diffusivity of α-Zr phase. Hydrogen increases the electrical resistivity of Zircaloy. Similar to thermal diffusivity results, the phase transformation causes a reversal of temperature-dependent trends in electrical resistivity results in Zircaloy-4 with higher hydrogen contents.

36 MATERIALS SCIENCE↗

Updated ASME design correlations and qualification plan for powder bed fusion 316H stainless steel

This report provides an update on the Advanced Materials and Manufacturing Technologies (AMMT) program effort to qualify Laser-Powder Bed Fusion (L-PBF) 316H stainless steel for use with the ASME Boiler & Pressure Vessel Code Section III, Division 5 rules. The report summarizes progress in testing and characterizing L-PBF material at elevated temperatures by providing preliminary design data for L-PBF 316H and by comparing the elevated temperature performance of the L-PBF material to wrought and conventional fusion welded 316H. The report then updates the initial AMMT qualification plan for L-PBF 316H, originally developed in 2023, to update the accelerated qualification strategy adopted in that plan to account for the new high temperature test data. The report also explores a few methods for further accelerating the qualification process using machine learning techniques to supplement the more conventional, empirical analysis methods typically used by ASME to correlate and extrapolate time-dependent material test data.

36 MATERIALS SCIENCE↗

Changes in Oculomotor Behavior and PVT Reaction Time During One Night of Sleep Deprivation

INTRODUCTION: The Psychomotor Vigilance Test (PVT) is a widely used objective measure of sustained attention, alertness, and fatigue. Previous studies consistently show that prolonged wakefulness and disrupted sleep patterns lead to increased lapses of attention and slower reaction times. These findings have critical implications for various professions, including healthcare, aviation, transportation safety, and spaceflight. In this study, we employed linear mixed models (LMM) to investigate hourly-measured 5-minute PVT reaction times during one night of total sleep deprivation. Our goal was to characterize the contribution of homeostatic sleep pressure (time awake) and circadian phase (salivary melatonin levels) on the dynamics of PVT reaction time. MATERIALS & METHODS: Data from twelve human participants were used in the analysis. Participants were healthy non-smokers, aged 18 and 40, with normal sleep habits defined as Pittsburg Sleep Quality Index scores < 5, and Morningness – Eveningness Questionnaire scores > 42 and < 58, respectively. Participants followed a constant-routine protocol and arrived at the NASA sleep laboratory approximately 1 to 2 hours after their habitual waking time. Tympanic temperature was measured every 30 min to provide a real-time estimate of circadian phase. Participants were required to stay awake throughout the laboratory experiment, which continued until their tympanic temperature had returned to baseline levels. This approach ensured that we could capture each participant's circadian trough and recovery, typically occurring between 24 and 26 hours after waking. The PVT and saliva melatonin were measured every hour, starting approximately from 3 hours after awakening. More information about the measurement protocol can be found in the publication by Stone et al. We used PVT reaction times between 100 and 10000 ms to calculate Reciprocal Reaction Times (RRTs). LMM models were employed to assess the contribution of time awake and saliva melatonin levels on changes in RRTs. Two models were created: the first model treated time awake as a fixed effect and the second model considered both the salivary melatonin level and time awake as fixed effects. Participants were included as a random effect in both models. We compared the models using an Analysis of Variance (ANOVA). The models were computed using the 'lmer' function from the 'lme4' package in R, and LMM p-values were determined using Satterthwaite's method. RESULTS: We found that time awake had a statistically significant effect on PVT RRTs during one night of sleep deprivation [χ2 = 2132.4, df = 1, p < 0.001]. In addition, the model including both saliva melatonin levels and time awake was statistically significantly better than one that uses just time awake [χ2 = 17.4, df = 1, p < 0.001]. DISCUSSION & CONCLUSION: As expected, these preliminary results suggest that the time awake and circadian phase each contribute to performance reductions as measured by the PVT. Future analyses will explore the contributions of other factors to the model.

sleep deprivation↗

Machine Learning for Automated Extraction of Building Geometry

As data science comes to buildings, the promise of using machine learning and novel sources of data has received much attention. Advances in machine learning and computer vision algorithms, combined with increased access to unstructured data (e.g., images and text), have created an opportunity for automated extraction of building characteristics – cost-effectively, and at scale. Acquisition of features such as footprint are time consuming and costly to acquire with today’s manual methods, but can be streamlined through intelligent software-based solutions applied to satellite images. When combined with aerial RGB and thermal images, full 3D geometries and thermal maps can be constructed to determine additional characteristics such as window to wall ratio, height, number of stories and envelope thermal characteristics. In this paper we present three contributions to accelerate these high potential opportunities: (1) a methodical analysis of how these features can be integrated into today’s simulation and data driven software tools to enhance efficiency measure identification and owner/operator decision making; (2) development and accuracy testing of open source deep neural network methods to extract building footprints from satellite imagery, including the curation and application of openly available GIS datasets for training and continued development by others; and (3) an open framework for drone-based image capture and creation of 3D building geometries. This work represents an important bridge between high-level studies that span diverse application areas and those that detail point solutions yet cannot be easily replicated or extended.

Touzani, Samir↗

Voltron Compatible Whole Building Root-Fault Detection and Diagnosis

Buildings consume more than 40% of primary energy in the U.S. Malfunctioning sensors, components, and control systems, as well as degrading components in Heating, Ventilating and Air-conditioning (HVAC) and lighting systems are main reasons for energy waste and unsatisfactory indoor environment. Field studies have shown that an energy saving of 5-30% of total building energy consumption and improved indoor air quality can be achieved by simply applying automated fault detection and diagnosis (AFDD), followed by corrections even if these are not done in real time. Extensive research has been made on the development of component level AFDD tools in the past two decades. However, in recent years, the development of method for a whole building level AFDD has received much attention. Here, a whole building fault refers to a fault happening in one subsystem but impact on more than one subsystems, or has a significant impact on building performances, such as building’s energy consumption. Existing component level AFDD tools could fail to detect or give false alarm for faults that trigger abnormalities in multiple subsystems. Isolating such whole building faults, which have impacts on multiple subsystems, is also challenging by simply using component based AFDD solutions. On the other hand, it might be more cost-effective for medium- or small-sized commercial buildings that have subsystems (e.g. primary cooling subsystem, air loop subsystem, heating subsystem, lighting subsystem, etc.) to install one whole building level AFDD tool, rather than many component level AFDD tools, to focus on significant faults. The focus of this project is firstly to provide training opportunities in building AFDD research area for undergraduate students and to develop a workforce pipeline in building sciences. The objective of the proposed research activities is to develop cost-effective and scalable AFDD solutions for whole building level faults. The developed AFDD solution is VOLTTRONTM compatible to further increase its plug-n-play capacity and the market penetration. VOLTTRONTM is an open-source innovative distributed control and sensing software platform developed by the Pacific Northwest National Laboratory. A VOLTTRONTM compatible AFDD tool brings the benefits such as low-cost deployment, and better building to grid integration, etc. The performance of the developed whole building AFDD tools is evaluated using data collected from a mixed-use medium-sized university campus building. Another significant goal of this project is to engage undergraduate students in research activities and expose them to building control and AFDD fields. A demonstration building – Nesbitt Hall at Drexel University is identified for this study. This demonstration building has a typical HVAC system configuration, which includes a water-cooled chiller system, three variable air volume (VAV) air handling unit (AHU) systems, and a hydronic heating system. Faults, that are expected to have a whole building level impact, are artificially implemented in this demonstration building for three different seasons. Building data have been collected in this project, which include fault free building automation system (BAS) data and BAS data that contain system behaviors caused by artificially implemented and naturally occurred faults. These collected data are used to evaluate the developed whole building AFDD tools. A data-driven method, which includes a weather and schedule-based pattern matching (WPM) method and Feature based Principal Component Analysis (FPCA) method, is developed for whole building level fault detection. Parameter sensitivity tests, including snapshot window size, data sample searching pool size, etc., are implemented to evaluate the impacts of various WPM-FPCA parameters. A data-driven and expert knowledge/rule-based method using Bayesian Network (BN) and WPM is developed for whole building level fault diagnosis. The developed WPM-BN method, which includes a two-layer BN structure model and BN parameters, are either learned from baseline data or developed from expert knowledge. The developed whole building AFDD tools are evaluated using BAS data collected from the demonstration building. Using the evaluation data set, the WPM-FPCA fault detection method is evaluated under two different Principle Component (PC) retention rates, i.e., 0.65 retention rate and 0.95 retention rate. It is observed that when using 0.65 PC retention rate, the WPM-FPCA fault detection method achieves a 79% fault detection rate (this rate reaches to 85% when including the six make-up fault test cases discussed in the next bullet), while when using 0.95 PC retention rate, the WPM-FPCA fault detection method achieves 100% fault detection rate (this rate remains to be 100% when including the six make-up fault test cases). Under both PC retention rates, the WPM-FPCA fault detection method demonstrates 0% false alarm rates. Using the evaluation data set, the developed WPM-BN fault diagnosis method has successfully diagnosed thirteen (13) faults out of the fourteen (14) fault test cases. The one mis-diagnosed case exhibits fewer fault symptoms in different subsystems, when compared with other fault test cases. Overall speaking, the developed whole building AFDD solutions demonstrate satisfactory fault detection/diagnosis accuracy with very low false alarm rate. Lastly, a survey and market study has been performed to understand the existing AFDD market, including its gaps and needs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Real-Time Parameter Estimation for Flexible Aircraft

A method for estimating aeroelastic stability and control derivatives for flexible aircraft is developed and demonstrated using flight test data for the X-56A subscale demonstrator. The method uses the equation-error approach with frequency-domain data, and can be applied post-flight or in real time during flight. The non-dimensional aeroelastic forces and moments and the explanatory variables (including generalized displacement, rate, and acceleration states for the vibration modes) are estimated using a finite element model and onboard sensor measurements in both a least squares and Kalman filtering framework. The data are then transformed into the frequency domain for parameter estimation using equation error. This method can result in a more efficient analysis than with other iterative methods, and can leverage existing statistical tools for model structure determination, data collinearity detection, combining multiple maneuvers or prior information, and others to improve model quality.

Grauer, Jared A.↗

Generic hypersonic inlet module analysis

A computational study associated with an internal inlet drag analysis was performed for a generic hypersonic inlet module. The purpose of this study was to determine the feasibility of computing the internal drag force for a generic scramjet engine module using computational methods. The computational study consisted of obtaining two-dimensional (2D) and three-dimensional (3D) computational fluid dynamics (CFD) solutions using the Euler and parabolized Navier-Stokes (PNS) equations. The solution accuracy was assessed by comparisons with experimenal pitot pressure data. The CFD analysis indicates that the 3D PNS solutions show the best agreement with experimental pitot pressure data. The internal inlet drag analysis consisted of obtaining drag force predictions based on experimental data and 3D CFD solutions. A comparative assessment of each of the drag prediction methods is made and the sensitivity of CFD drag values to computational procedures is documented. The analysis indicates that the CFD drag predictions are highly sensitive to the computational procedure used.

Cockrell, Charles E., Jr.↗

Generic Hypersonic Inlet Module Analysis

A computational study associated with an internal inlet drag analysis was performed for a generic hypersonic inlet module. The purpose of this study was to determine the feasibility of computing the internal drag force for a generic scramjet engine module using computational methods. The computational study consisted of obtaining two-dimensional (2D) and three-dimensional (3D) computational fluid dynamics (CFD) solutions using the Euler and parabolized Navier-Stokes (PNS) equations. The solution accuracy was assessed by comparisons with experimental pitot pressure data. The CFD analysis indicates that the 3D PNS solutions show the best agreement with experimental pitot pressure data. The internal inlet drag analysis consisted of obtaining drag force predictions based on experimental data and 3D CFD solutions. A comparative assessment of each of the drag prediction methods is made and the sensitivity of CFD drag values to computational procedures is documented. The analysis indicates that the CFD drag predictions are highly sensitive to the computational procedure used.

Cockrell, Chares E., Jr.↗

Measurements of mantle wave velocities and inversion for lateral heterogeneity and anisotropy. II - Analysis by the single-station method

In the present investigation, the single-station method reported by Brune et al. (1960) is utilized for an analysis of long-period Love(G) and Rayleigh(R) waves recorded on digital seismic networks. The analysis was conducted to study the lateral heterogeneity of surface wave velocities. The data set is examined, and a description is presented of the single-station method. Attention is given to an error analysis for velocity measurements, the estimation of geographical distribution of surface wave velocities, the global distribution of surface wave velocities, and the correlation between the surface wave velocities and the heat flow on the geoid. The conducted measurements and inversions of surface wave velocities are used as a basis to derive certain conclusions. It is found that the application of the single-station method to long-period surface waves recorded on digital networks makes it possible to reach an accuracy level comparable to great circle velocity measurements.

Nakanishi, I.↗

Evaluation of Piloted Inputs for Onboard Frequency Response Estimation

Frequency response estimation results are presented using piloted inputs and a real-time estimation method recently developed for multisine inputs. A nonlinear simulation of the F-16 and a Piper Saratoga research aircraft were subjected to different piloted test inputs while the short period stabilator/elevator to pitch rate frequency response was estimated. Results show that the method can produce accurate results using wide-band piloted inputs instead of multisines. A new metric is introduced for evaluating which data points to include in the analysis and recommendations are provided for applying this method with piloted inputs.

Grauer, Jared A.↗