Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “filter methods”

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,459 records · Page 81

Numerical Simulation of Turbulent Jets with Rectangular Cross-Section

Three-dimensional turbulent jets with rectangular cross-section are simulated with a finite-difference numerical method. The full Navier-Stokes equations are solved at low Reynolds numbers, whereas at the high Reynolds numbers filtered forms of the equations are solved along with a sub-grid scale model to approximate effects of the unresolved scales. A 2-N storage, third-order Runge-Kutta scheme is used for temporal discretization and a fourth-order compact scheme is used for spatial discretization. Computations are performed for different inlet conditions which represent different types of jet forcing. The phenomenon of axis-switching is observed, and it is confirmed that this is based on self-induction of the vorticity field. Budgets of the mean streamwise velocity show that convection is balanced by gradients of the Reynolds stresses and the pressure.

Wilson, Robert V.↗

Subauditory Speech Recognition based on EMG/EPG Signals

Sub-vocal electromyogram/electro palatogram (EMG/EPG) signal classification is demonstrated as a method for silent speech recognition. Recorded electrode signals from the larynx and sublingual areas below the jaw are noise filtered and transformed into features using complex dual quad tree wavelet transforms. Feature sets for six sub-vocally pronounced words are trained using a trust region scaled conjugate gradient neural network. Real time signals for previously unseen patterns are classified into categories suitable for primitive control of graphic objects. Feature construction, recognition accuracy and an approach for extension of the technique to a variety of real world application areas are presented.

Jorgensen, Charles↗

Measurement and Modeling of Vertically Resolved Aerosol Optical Properties and Radiative Fluxes Over the ARM SGP Site

In order to meet one of its goals - to relate observations of radiative fluxes and radiances to the atmospheric composition - the Department of Energy's Atmospheric Radiation Measurement (ARM) program has pursued measurements and modeling activities that attempt to determine how aerosols impact atmospheric radiative transfer, both directly and indirectly. However, significant discrepancies between aerosol properties measured in situ or remotely remain. One of the objectives of the Aerosol Intensive Operational Period (TOP) conducted by ARM in May 2003 at the ARM Southern Great Plains (SGP) site in north central Oklahoma was to examine and hopefully reduce these differences. The IOP involved airborne measurements from two airplanes over the heavily instrumented SGP site. We give an overview of airborne results obtained aboard the Center for Interdisciplinary Remotely-Piloted Aircraft Studies (CIRPAS) Twin Otter aircraft. The Twin Otter performed 16 research flights over the SGP site. The aircraft carried instrumentation to perform in-situ measurements of aerosol absorption, scattering, extinction and particle size. This included such novel techniques as the photoacoustic and cavity ring-down methods for in-situ absorption (675 nm) and extinction (675 and 1550 nm) and a new multiwavelength, filter-based absorption photometer (467, 530, 660 nm). A newly developed instrument measured cloud condensation nucleus concentration (CCN) concentrations at two supersaturation levels. Aerosol optical depth and extinction (354-2139 nm) were measured with the NASA Ames Airborne Tracking 14-channel sunphotometer. Furthermore, up-and downwelling solar (broadband and spectral) and infrared radiation were measured using seven individual radiometers. Three up-looking radiometers werer mounted on a newly developed stabilized platform, keeping the instruments level up to aircraft pitch and roll angles of approximately 10(exp 0). This resulted in unprecedented continuous vertical profiles of radiative fluxes, which we will compare to modeled fluxes using the aforementioned data as input.

Schmid, B.↗

Methods and Systems for Characterization of an Anomaly Using Infrared Flash Thermography

A method for characterizing an anomaly in a material comprises (a) extracting contrast data; (b) measuring a contrast evolution; (c) filtering the contrast evolution; (d) measuring a peak amplitude of the contrast evolution; (d) determining a diameter and a depth of the anomaly, and (e) repeating the step of determining the diameter and the depth of the anomaly until a change in the estimate of the depth is less than a set value. The step of determining the diameter and the depth of the anomaly comprises estimating the depth using a diameter constant C.sub.D equal to one for the first iteration of determining the diameter and the depth; estimating the diameter; and comparing the estimate of the depth of the anomaly after each iteration of estimating to the prior estimate of the depth to calculate the change in the estimate of the depth of the anomaly.

Koshti, Ajay M.↗

GNSS Ephemeris with Graceful Degradation and Measurement Fusion

A method for providing an extended propagation ephemeris model for a satellite in Earth orbit, the method includes obtaining a satellite's orbital position over a first period of time, applying a least square estimation filter to determine coefficients defining osculating Keplarian orbital elements and harmonic perturbation parameters associated with a coordinate system defining an extended propagation ephemeris model that can be used to estimate the satellite's position during the first period, wherein the osculating Keplarian orbital elements include semi-major axis of the satellite (a), eccentricity of the satellite (e), inclination of the satellite (i), right ascension of ascending node of the satellite (.OMEGA.), true anomaly (.theta.*), and argument of periapsis (.omega.), applying the least square estimation filter to determine a dominant frequency of the true anomaly, and applying a Fourier transform to determine dominant frequencies of the harmonic perturbation parameters.

Garrison, James Levi↗

Model-Based Fatigue Prognosis of Fiber-Reinforced Laminates Exhibiting Concurrent Damage Mechanisms

Prognostics of large composite structures is a topic of increasing interest in the field of structural health monitoring for aerospace, civil, and mechanical systems. Along with recent advancements in real-time structural health data acquisition and processing for damage detection and characterization, model-based stochastic methods for life prediction are showing promising results in the literature. Among various model-based approaches, particle-filtering algorithms are particularly capable in coping with uncertainties associated with the process. These include uncertainties about information on the damage extent and the inherent uncertainties of the damage propagation process. Some efforts have shown successful applications of particle filtering-based frameworks for predicting the matrix crack evolution and structural stiffness degradation caused by repetitive fatigue loads. Effects of other damage modes such as delamination, however, are not incorporated in these works. It is well established that delamination and matrix cracks not only co-exist in most laminate structures during the fatigue degradation process but also affect each other's progression. Furthermore, delamination significantly alters the stress-state in the laminates and accelerates the material degradation leading to catastrophic failure. Therefore, the work presented herein proposes a particle filtering-based framework for predicting a structure's remaining useful life with consideration of multiple co-existing damage-mechanisms. The framework uses an energy-based model from the composite modeling literature. The multiple damage-mode model has been shown to suitably estimate the energy release rate of cross-ply laminates as affected by matrix cracks and delamination modes. The model is also able to estimate the reduction in stiffness of the damaged laminate. This information is then used in the algorithms for life prediction capabilities. First, a brief summary of the energy-based damage model is provided. Then, the paper describes how the model is embedded within the prognostic framework and how the prognostics performance is assessed using observations from run-to-failure experiments

Prognostics↗

Robotic Specialization in Autonomous Robotic Structural Assembly

Robotic in-space assembly of large space structures is a long-term NASA goal to reduce launch costs and enable larger scale missions. Recently, researchers have proposed using discrete lattice building blocks and co-designed robots to build high-performance, scalable primary structure for various on-orbit and surface applications. These robots would locomote on the lattice and work in teams to build and reconfigure building-blocks into functional structure. However, the most reliable and efficient robotic system architecture, characterized by the number of different robotic 'species' and the allocation of functionality between species, is an open question. To address this problem, we decompose the robotic building-block assembly task into functional primitives and, in simulation, study the performance of the the variety of possible resulting architectures. For a set consisting of five process types (move self, move block, move friend, align bock, fasten block), we describe a method of feature space exploration and ranking based on energy and reliability cost functions. The solution space is enumerated, filtered for unique solutions, and evaluated against energy and reliability cost functions for various simulated build sizes. We find that a 2 species system, dividing the five mentioned process types between one unit cell transport robot and one fastening robot, results in the lowest energy cost system, at some cost to reliability. This system enables fastening functionality to occupy the build front while reducing the need for that functional mass to travel back and forth from a feed station. Because the details of a robot design affect the weighting and final allocation of functionality, a sensitivity analysis was conducted to evaluate the effect of changing mass allocations on architecture performance. Future systems with additional functionalities such as repair, inspection, and others may use this process to analyze and determine alternative robot architectures.

Bernus, Borbala↗

Applying Time-Parallelization to Turbulent Flows

Parallelization of the temporal domain is explored for the solution of turbulent flows. Multigrid reduction-in-time (MGRIT) is used to advance the large-scale fluid dynamics in time sequentially on the coarsest space-time grid but propagate the information in time parallel on all other levels. The goal of this process is to accurately and efficiently resolve the coarse-scale turbulence structure and use that to drive the fine-scales of the turbulent flow. The extra forcing from nonlinear multigrid facilitates the coupling and interaction between fine and coarse scales, through which the multiscale nonlinear physics is properly captured. Adaptive mesh refinement is employed to finely resolve only the regions with strong gradients, which provides further computational efficiency. The underlying computational fluid dynamics solver is a fourth-order finite-volume scheme with the standard 4-stage Runge-Kutta method. An advanced approach is devised and implemented to enable MGRIT to solve highly turbulent flows successfully. Furthermore, the method is applied to solve a Taylor-Green vortex problem and a doubleshear-layer turbulent mixing flow. Results are promising, validating that MGRIT with the filtering approach has the potential to efficiently solve general turbulent flows.

Computational Fluid Dynamics↗

Orbit Determination for the Lunar Reconnaissance Orbiter Using an Extended Kalman Filter

Orbit determination (OD) analysis results are presented for the Lunar Reconnaissance Orbiter (LRO) using a commercially available Extended Kalman Filter, Analytical Graphics' Orbit Determination Tool Kit (ODTK). Process noise models for lunar gravity and solar radiation pressure (SRP) are described and OD results employing the models are presented. Definitive accuracy using ODTK meets mission requirements and is better than that achieved using the operational LRO OD tool, the Goddard Trajectory Determination System (GTDS). Results demonstrate that a Vasicek stochastic model produces better estimates of the coefficient of solar radiation pressure than a Gauss-Markov model, and prediction accuracy using a Vasicek model meets mission requirements over the analysis span. Modeling the effect of antenna motion on range-rate tracking considerably improves residuals and filter-smoother consistency. Inclusion of off-axis SRP process noise and generalized process noise improves filter performance for both definitive and predicted accuracy. Definitive accuracy from the smoother is better than achieved using GTDS and is close to that achieved by precision OD methods used to generate definitive science orbits. Use of a multi-plate dynamic spacecraft area model with ODTK's force model plugin capability provides additional improvements in predicted accuracy.

Orbit Determination↗

Micropulsed Plasma Thrusters for Attitude Control of a Low-Earth-Orbiting CubeSat

This study presents a 3-Unit CubeSat design with commercial-off-the-shelf hardware, Teflon-fueled micropulsed plasma thrusters, and an attitude determination and control approach. The micropulsed plasma thruster is sized by the impulse bit and pulse frequency required for continuous compensation of expected maximum disturbance torques at altitudes between 400 and 1000 km, as well as to perform stabilization of up to 20 deg /s and slew maneuvers of up to 180 deg. The study involves realistic power constraints anticipated on the 3-Unit CubeSat. Attitude estimation is implemented using the q method for static attitude determination of the quaternion using pairs of the spacecraft-sun and magnetic-field vectors. The quaternion estimate and the gyroscope measurements are used with an extended Kalman filter to obtain the attitude estimates. Proportional-derivative control algorithms use the static attitude estimates in order to calculate the torque required to compensate for the disturbance torques and to achieve specified stabilization and slewing maneuvers or combinations. The controller includes a thruster-allocation method, which determines the optimal utilization of the available thrusters and introduces redundancy in case of failure. Simulation results are presented for a 3-Unit CubeSat under detumbling, pointing, and pointing and spinning scenarios, as well as comparisons between the thruster-allocation and the paired-firing methods under thruster failure.

Gatsonis, Nikolaos A.↗

Quantitative Gas Property Measurements By Filtered Rayleigh Scattering: A Review

Filtered Rayleigh scattering (FRS) is a laser-based diagnostic technique used to nonintrusively quantify various thermodynamic properties of a light-scattering gas. The backbone of FRS is the molecular filtering of Rayleigh scattered light. This concept was initially introduced by the atmospheric LIDAR community before being adopted within the aerospace research field in the early 1990s. Since then, FRS has matured into a versatile quantitative diagnostic tool and has found use in a variety of flow regimes ranging from sub- to supersonic speeds in both reacting and nonreacting environments. This adoption can be attributed to the wealth of information that can be obtained via FRS, including the gas density, pressure, temperature, velocity, species composition, or, in some cases, several of these properties at once. This article reviews the current state of FRS methodology in recovering such gas properties. As knowledge of the fundamentals of Rayleigh scattering and spectral light filtering is crucial to the design of an FRS experiment, we begin by briefly reviewing these areas. Subsequently, we conduct a survey of experimental design strategies, assumptions, and data reduction methods used to measure different gas properties using FRS. We conclude the review with a short discussion on quantification of experimental uncertainty and future trends in FRS.

filtered Rayleigh scattering↗

A dynamic likelihood approach to filtering transport processes: advection-diffusion dynamics

A Bayesian data assimilation scheme is formulated for advection-dominated advective and diffusive evolutionary problems, based upon the Dynamic Likelihood (DLF) approach to filtering. The DLF was developed specifically for hyperbolic problems –waves–, and in this paper, it is extended via a split step formulation, to handle advection-diffusion problems. In the dynamic likelihood approach, observations and their statistics are used to propagate probabilities along characteristics, evolving the likelihood in time. The estimate posterior thus inherits phase information. For advection-diffusion the advective part of the time evolution is handled on the basis of observations alone, while the diffusive part is informed through the model as well as observations. We expect, and indeed show here, that in advection-dominated problems, the DLF approach produces better estimates than other assimilation approaches, particularly when the observations are sparse and have low uncertainty. The added computational expense of the method is cubic in the total number of observations over time, which is on the same order of magnitude as a standard Kalman filter and can be mitigated by bounding the number of forward propagated observations, discarding the least informative data.

97 MATHEMATICS AND COMPUTING↗

Data Fusion of In-Flight Aerothermodynamic Heating Measurements Using Kalman Filtering

On February 18th, 2021, the Mars 2020 entry system successfully delivered the Perserverance rover to the surface of Mars at Jezero Crater. The entry capsule carried a set of instrumentation installed on the heat shield and backshell, named the Mars 2020 Entry, Descent, and Landing Instrumentation. The instruments include pressure transducers, thermocouples, heat flux gauges, and radiometers to measure the aerodynamic and aerothermodynamic performance of the entry vehicle. Three of these sensors, a thermocouple plug, heat flux gauge, and a radiometer, are located in close proximity on the backshell. Each sensor is exposed to roughly the same environment, but measured these environments in different ways, each with its own set of modeling and measurement error complications. This paper develops a method for blending each of these measurements together in a single algorithm to produce estimates of the aerothermodynamic environments at that location on the backshell. The approach makes use of the Kalman filter methodology for solving state estimation problems. The filter has a predictor/corrector structure in which simplified process models are used to propagate the aeroheating states forward through time, which are then updated based on the measurement data. The method has been tested on simulated datasets.

Christopher D. Karlgaard↗

The Behavior of Filters and Smoothers for Strongly Nonlinear Dynamics

The Kalman filter is the optimal filter in the presence of known Gaussian error statistics and linear dynamics. Filter extension to nonlinear dynamics is non trivial in the sense of appropriately representing high order moments of the statistics. Monte Carlo, ensemble-based, methods have been advocated as the methodology for representing high order moments without any questionable closure assumptions (e.g., Miller 1994). Investigation along these lines has been conducted for highly idealized dynamics such as the strongly nonlinear Lorenz (1963) model as well as more realistic models of the oceans (Evensen and van Leeuwen 1996) and atmosphere (Houtekamer and Mitchell 1998). A few relevant issues in this context are related to the necessary number of ensemble members to properly represent the error statistics and, the necessary modifications in the usual filter equations to allow for correct update of the ensemble members (Burgers 1998). The ensemble technique has also been applied to the problem of smoothing for which similar questions apply. Ensemble smoother examples, however, seem to quite puzzling in that results of state estimate are worse than for their filter analogue (Evensen 1997). In this study, we use concepts in probability theory to revisit the ensemble methodology for filtering and smoothing in data assimilation. We use Lorenz (1963) model to test and compare the behavior of a variety implementations of ensemble filters. We also implement ensemble smoothers that are able to perform better than their filter counterparts. A discussion of feasibility of these techniques to large data assimilation problems will be given at the time of the conference.

Zhu, Yanqiu↗

Methodology to Define Delivery Accuracy Under Current Day ATC Operations

In order to enable arrival management concepts and solutions in a NextGen environment, ground- based sequencing and scheduling functions have been developed to support metering operations in the National Airspace System. These sequencing and scheduling algorithms as well as tools are designed to aid air traffic controllers in developing an overall arrival strategy. The ground systems being developed will support the management of aircraft to their Scheduled Times of Arrival (STAs) at flow-constrained meter points. This paper presents a methodology for determining the undelayed delivery accuracy for current day air traffic control operations. This new method analyzes the undelayed delivery accuracy at meter points in order to understand changes of desired flow rates as well as enabling definition of metrics that will allow near-future ground automation tools to successfully achieve desired separation at the meter points. This enables aircraft to meet their STAs while performing high precision arrivals. The research presents a possible implementation that would allow delivery performance of current tools to be estimated and delivery accuracy requirements for future tools to be defined, which allows analysis of Estimated Time of Arrival (ETA) accuracy for Time-Based Flow Management (TBFM) and the FAA's Traffic Management Advisor (TMA). TMA is a deployed system that generates scheduled time-of-arrival constraints for en- route air traffic controllers in the US. This new method of automated analysis provides a repeatable evaluation of the delay metrics for current day traffic, new releases of TMA, implementation of different tools, and across different airspace environments. This method utilizes a wide set of data from the Operational TMA-TBFM Repository (OTTR) system, which processes raw data collected by the FAA from operational TMA systems at all ARTCCs in the nation. The OTTR system generates daily reports concerning ATC status, intent and actions. Due to its availability, ease of use, and vast collection of data across several airspaces it was determined that the OTTR data set would be the best method to utilize moving forward with this analysis. The particular variables needed for further analysis were determined along with the necessary OTTR reports, by working closely with the repository team additional analysis reports were developed that provided key ETA and STA information at the freeze horizon. One major benefit of the OTTR data is that using the correct reports the data across several airports could be analyzed over large periods of time. The OTTR data processes the TBFM data daily and is stored in various formats across several airspaces. This allowed us to develop our own parsing methods and raw data processing that would not rely on other computationally expensive tools that perform more in depth analysis of similar sets of data. The majority of this work consisted of the development of the ability to filter flights to create a subset of flights that could be considered undelayed, which is defined as a flight at the freeze horizon with an ETA and STA difference that was minimal or close to zero. This was a broad method that allowed the consideration of a large data set which consisted of all the traffic across a two month period in 2013, the hottest and coldest months, arriving into four airports: George Bush Intercontinental, Denver International, Los Angeles International, and Phoenix Sky Harbor.

delivery accuracy↗

Polarimetric Imaging using Two Photoelastic Modulators

A method of polarimetric imaging, now undergoing development, involves the use of two photoelastic modulators in series, driven at equal amplitude but at different frequencies. The net effect on a beam of light is to cause (1) the direction of its polarization to rotate at the average of two excitation frequencies and (2) the amplitude of its polarization to be modulated at the beat frequency (the difference between the two excitation frequencies). The resulting modulated optical light beam is made to pass through a polarizing filter and is detected at the beat frequency, which can be chosen to equal the frame rate of an electronic camera or the rate of sampling the outputs of photodetectors in an array. The method was conceived to satisfy a need to perform highly accurate polarimetric imaging, without cross-talk between polarization channels, at frame rates of the order of tens of hertz. The use of electro-optical modulators is necessitated by a need to obtain accuracy greater than that attainable by use of static polarizing filters over separate fixed detectors. For imaging, photoelastic modulators are preferable to such other electrio-optical modulators as Kerr cells and Pockels cells in that photoelastic modulators operate at lower voltages, have greater angular acceptances, and are easier to use. Prior to the conception of the present method, polarimetric imaging at frame rates of tens of hertz using photoelastic modulators was not possible because the resonance frequencies of photoelastic modulators usually lie in the range from about 20 to about 100 kHz.

Wang, Yu↗

Baseband pulse shaping techniques for nonlinearly amplified pi/4-QPSK and QAM systems

A new generation of multi-stage pi/4-shifted QPSK and of superposed quadrature-amplitude-modulated (SQAM) modulators-coherent demodulators (modems) and of continuous phase modulated (CPM)-gaussian premodulation filtered minimum-shift-keying (MGMSK) systems is proposed and studied. These modems will lead to bandwidth and power efficient satellite communications systems designs. As an illustrative application, a baseband processing technique pi/4-controlled transition PSK (pi/4-CTPSK) is described. To develop a cost and power efficient design strategy, we assume that nonlinear, fully saturated high power amplifiers (HPA) are utilized in the satellite earth station transmitter and in the satellite transponder. Modem structures which could lead to application specific integrated circuit (ASIC) satellite on-board processing universal modem applications are also considered. Multistate GMSK (i.e., MGMSK) signal generation methods by means of two or more RF combined nonlinearly amplified SQAM modems and by one multistate (in-phase and quadrature-baseband premodulation filtered-superposed) SQAM architecture and one RF nonlinear amplifier are studied. During the SQAM modem development phase we investigate the potential system advantages of the pi/4-shifted logic. The bandwidth efficiency of the proposed multistate GMSK and baseband filtered PAM-FM modulator (a new class in the CPM family) will be significantly higher than that of conventional G-MSK systems. To optimize the practical P(sub e) = f((E sub b)/(N sub o)) performance we consider improved coherent demodulation MGMSK structures such as deviated-frequency locking coherent demodulators. For relative low bit rate SATCOM applications, e.g., bit rates less than 300 kb/s, phase noise tracking cancellation (for fixed site earth station) and phase noise cancellation as well as Doppler compensation (for satellite to mobile earth station) applications may be required. We study digital channel sounding methods which could cancel the phase noise-caused degradations of CPM and GMSK modems.

Feher, Kamilo↗

Underlying Fundamentals of Kalman Filtering for River Network Modeling

The grand challenge of producing hydrometeorological estimates every time and everywhere has motivated the fusion of sparse observations with dense numerical models, with a particular interest on discharge in river modeling. Ensemble methods are largely preferred as they enable the estimation of error properties, but at the expense of computational load and generally with underestimations. These imperfect stochastic estimates motivate the use of correction methods, that is, error localization and inflation, although the physical justifications for their optimality are limited. The purpose of this study is to use one of the simplest forms of data assimilation when applied to river modeling and reveal the underlying mechanisms impacting its performance. Our framework based on assimilating daily averaged in situ discharge measurements to correct daily averaged runoff was tested over a 4-yr case study of two rivers in Texas. Results show that under optimal conditions of inflation and localization, discharge simulations are consistently improved such that the mean values of Nash–Sutcliffe efficiency are enhanced from211.32 to 0.55 at observed gauges and from212.24 to21.10 at validation gauges. Yet, parameters controlling the inflation and the localization have a large impact on the performance. Further investigations of these sensitivities showed that optimal inflation occurs when compensating exactly for discrepancies in the magnitude of errors while optimal localization matches the distance traveled during one assimilation window. These results may be applicable to more advanced data assimilation methods as well as for larger applications motivated by upcoming river-observing satellite missions, such as NASA’s Surface Water and Ocean Topography mission.

Streamflow↗