Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear filtering”

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 631 records · Page 35

Nonlinear Coherent Optical Image Processing Using Logarithmic Transmittance of Bacteriorhodopsin Films

The transmission properties of some bacteriorhodopsin-film spatial light modulators are uniquely suited to allow nonlinear optical image-processing operations to be applied to images with multiplicative noise characteristics. A logarithmic amplitude-transmission characteristic of the film permits the conversion of multiplicative noise to additive noise, which may then be linearly filtered out in the Fourier plane of the transformed image. I present experimental results demonstrating the principle and the capability for several different image and noise situations, including deterministic noise and speckle. The bacteriorhodopsin film studied here displays the logarithmic transmission response for write intensities spanning a dynamic range greater than 2 orders of magnitude.

Downie, John D.↗

Development of Computational Aeroacoustics Code for Jet Noise and Flow Prediction

Accurate prediction of jet fan and exhaust plume flow and noise generation and propagation is very important in developing advanced aircraft engines that will pass current and future noise regulations. In jet fan flows as well as exhaust plumes, two major sources of noise are present: large-scale, coherent instabilities and small-scale turbulent eddies. In previous work for the NASA Glenn Research Center, three strategies have been explored in an effort to computationally predict the noise radiation from supersonic jet exhaust plumes. In order from the least expensive computationally to the most expensive computationally, these are: 1) Linearized Euler equations (LEE). 2) Very Large Eddy Simulations (VLES). 3) Large Eddy Simulations (LES). The first method solves the linearized Euler equations (LEE). These equations are obtained by linearizing about a given mean flow and the neglecting viscous effects. In this way, the noise from large-scale instabilities can be found for a given mean flow. The linearized Euler equations are computationally inexpensive, and have produced good noise results for supersonic jets where the large-scale instability noise dominates, as well as for the tone noise from a jet engine blade row. However, these linear equations do not predict the absolute magnitude of the noise; instead, only the relative magnitude is predicted. Also, the predicted disturbances do not modify the mean flow, removing a physical mechanism by which the amplitude of the disturbance may be controlled. Recent research for isolated airfoils' indicates that this may not affect the solution greatly at low frequencies. The second method addresses some of the concerns raised by the LEE method. In this approach, called Very Large Eddy Simulation (VLES), the unsteady Reynolds averaged Navier-Stokes equations are solved directly using a high-accuracy computational aeroacoustics numerical scheme. With the addition of a two-equation turbulence model and the use of a relatively coarse grid, the numerical solution is effectively filtered into a directly calculated mean flow with the small-scale turbulence being modeled, and an unsteady large-scale component that is also being directly calculated. In this way, the unsteady disturbances are calculated in a nonlinear way, with a direct effect on the mean flow. This method is not as fast as the LEE approach, but does have many advantages to recommend it; however, like the LEE approach, only the effect of the largest unsteady structures will be captured. An initial calculation was performed on a supersonic jet exhaust plume, with promising results, but the calculation was hampered by the explicit time marching scheme that was employed. This explicit scheme required a very small time step to resolve the nozzle boundary layer, which caused a long run time. Current work is focused on testing a lower-order implicit time marching method to combat this problem.

Keith, Theo G., Jr.↗

Adaptive Control for Uncertain Nonlinear Multi-Input Multi-Output Systems

Systems and methods of adaptive control for uncertain nonlinear multi-input multi-output systems in the presence of significant unmatched uncertainty with assured performance are provided. The need for gain-scheduling is eliminated through the use of bandwidth-limited (low-pass) filtering in the control channel, which appropriately attenuates the high frequencies typically appearing in fast adaptation situations and preserves the robustness margins in the presence of fast adaptation.

Cao, Chengyu↗

Bounded state space

This investigation is divided functionally into three different areas: (1) study of bounded state space, (2) nonlinear smoothing theory, and (3) system identification. (1) Study of bounded state space: necessary and sufficient conditions for an optimal control are obtained for a bounded state space optimal control problem. The difficulty of determining the so-called jump conditions is eliminated; however, the problem of determining the points where the response either enters or leaves the boundary still remains unsolved. (2) Nonlinear smoothing theory: nonlinear fixed-interval, fixed-point and fixed-lag smoothing of a random signal generated by a stochastic differential equation are investigated. Results on the asymptotic stability of a linear constant-parameter fixed-interval smoothing filter are obtained. (3) System identification: a particular stochastic modelling problem is solved. An Ito stochastic integral equation is used to mathematically model a black box having multiple inputs and multiple outputs. A new method for identifying system parameters is presented.

Eyman, E. D.↗

A class of all digital phase locked loops - Modeling and analysis

An all digital phase locked loop which tracks the phase of the incoming signal once per carrier cycle is proposed. The different elements and their functions, and the phase lock operation are explained in detail. The general digital loop operation is governed by a nonlinear difference equation from which a suitable model is developed. The lock range for the general model is derived. The performance of the digital loop for phase step and frequency step inputs for different levels of quantization without loop filter are studied. The analytical results are checked by simulating the actual system on the digital computer.

Reddy, C. P.↗

A Complementary Note to 'A Lag-1 Smoother Approach to System-Error Estimation': The Intrinsic Limitations of Residual Diagnostics

Recently, this author studied an approach to the estimation of system error based on combining observation residuals derived from a sequential filter and fixed lag-1 smoother. While extending the methodology to a variational formulation, experimenting with simple models and making sure consistency was found between the sequential and variational formulations, the limitations of the residual-based approach came clearly to the surface. This note uses the sequential assimilation application to simple nonlinear dynamics to highlight the issue. Only when some of the underlying error statistics are assumed known is it possible to estimate the unknown component. In general, when considerable uncertainties exist in the underlying statistics as a whole, attempts to obtain separate estimates of the various error covariances are bound to lead to misrepresentation of errors. The conclusions are particularly relevant to present-day attempts to estimate observation-error correlations from observation residual statistics. A brief illustration of the issue is also provided by comparing estimates of error correlations derived from a quasi-operational assimilation system and a corresponding Observing System Simulation Experiments framework.

Fixed-lag smoothing↗

Near optimum digital phase locked loops.

Near optimum digital phase locked loops are derived utilizing nonlinear estimation theory. Nonlinear approximations are employed to yield realizable loop structures. Baseband equivalent loop gains are derived which under high signal to noise ratio conditions may be calculated off-line. Additional simplifications are made which permit the application of the Kalman filter algorithms to determine the optimum loop filter. Performance is evaluated by a theoretical analysis and by simulation. Theoretical and simulated results are discussed and a comparison to analog results is made.

Polk, D. R.↗

A guidance and navigation system for continuous low-thrust vehicles

A midcourse guidance and navigation system for continuous low thrust vehicles was developed. The equinoctial elements are the state variables. Uncertainties are modelled statistically by random vector and stochastic processes. The motion of the vehicle and the measurements are described by nonlinear stochastic differential and difference equations respectively. A minimum time trajectory is defined; equations of motion and measurements are linearized about this trajectory. An exponential cost criterion is constructed and a linear feedback quidance law is derived. An extended Kalman filter is used for state estimation. A short mission using this system is simulated. It is indicated that this system is efficient for short missions, but longer missions require accurate trajectory and ground based measurements.

Jack-Chingtse, C.↗

Analysis of 12-700 micron reflectance spectra of three optical black samples

Normal incidence, specular reflectance spectra of three samples of optical black coatings have been measured in 10 filter passbands at effective wavelengths between 12 and 700 microns. Strong absorption, due to amorphous silicate material in the coatings, is found in the 15-50 micron region of each spectrum. Reflecting-layer model calculations are fitted by a nonlinear least-squares routine to the spectra at wavelengths beyond the absorption band to determine far-infrared optical parameters. It is found that, with a similar wavelength dependence, the absorption coefficient of a new black coating (ECP 2200) made by the Minnesota Mining and Manufacturing Company (3M) is less than 1/10 that of the discontinued coating, 3M Black Velvet Nextel. It is also shown that the extinction effect of a very rough substrate is negligible compared with the effects of increased upper surface roughness and coating thickness attendant upon fully covering a very rough substrate.

Smith, S. M.↗

Investigation of loss processes of Tm and Tm,Ho in YAG

The loss of excitation from various manifolds of Tm and Tm,Ho in YAG as a function of temperature and concentration is studied. Two probable loss mechanisms - a Tm up-conversion and a Ho up-conversion - are identified. A 785-nm CW diode laser with 400-nW peak power was focused to a small spot on the sample. The emission from the sample observed at 90 deg was monitored through a monochromator with slits open to 3 mm. Intensity of emission was measured by varying the power of the excitation source using a set of neutral density filters. Power is reported as the percentage of the peak power, and the intensity curves were normalized below 20 percent of transmission. The fact that there is emission above the pump energy indicates an up-conversion from excited manifolds. Nonlinear changes in the intensity of the emission from the Tm 3F4 manifold with the pump power reveals a loss of excitation from this manifold. The linear dependence of the 5I7 manifold emission with pump power at low Tm and high Ho concentrations and the gain of energy in the 5I6 manifold of Ho indicate that the 5I7 manifold loss is due to the coupling of Tm and Ho ions.

Armagan, G.↗

Ares I Flight Control System Design

The Ares I launch vehicle represents a challenging flex-body structural environment for flight control system design. This paper presents a design methodology for employing numerical optimization to develop the Ares I flight control system. The design objectives include attitude tracking accuracy and robust stability with respect to rigid body dynamics, propellant slosh, and flex. Under the assumption that the Ares I time-varying dynamics and control system can be frozen over a short period of time, the flight controllers are designed to stabilize all selected frozen-time launch control systems in the presence of parametric uncertainty. Flex filters in the flight control system are designed to minimize the flex components in the error signals before they are sent to the attitude controller. To ensure adequate response to guidance command, step response specifications are introduced as constraints in the optimization problem. Imposing these constraints minimizes performance degradation caused by the addition of the flex filters. The first stage bending filter design achieves stability by adding lag to the first structural frequency to phase stabilize the first flex mode while gain stabilizing the higher modes. The upper stage bending filter design gain stabilizes all the flex bending modes. The flight control system designs provided here have been demonstrated to provide stable first and second stage control systems in both Draper Ares Stability Analysis Tool (ASAT) and the MSFC 6DOF nonlinear time domain simulation.

Jang, Jiann-Woei↗

Anomaly Detection in Test Equipment via Sliding Mode Observers

Nonlinear observers were originally developed based on the ideas of variable structure control, and for the purpose of detecting disturbances in complex systems. In this anomaly detection application, these observers were designed for estimating the distributed state of fluid flow in a pipe described by a class of advection equations. The observer algorithm uses collected data in a piping system to estimate the distributed system state (pressure and velocity along a pipe containing liquid gas propellant flow) using only boundary measurements. These estimates are then used to further estimate and localize possible anomalies such as leaks or foreign objects, and instrumentation metering problems such as incorrect flow meter orifice plate size. The observer algorithm has the following parts: a mathematical model of the fluid flow, observer control algorithm, and an anomaly identification algorithm. The main functional operation of the algorithm is in creating the sliding mode in the observer system implemented as software. Once the sliding mode starts in the system, the equivalent value of the discontinuous function in sliding mode can be obtained by filtering out the high-frequency chattering component. In control theory, "observers" are dynamic algorithms for the online estimation of the current state of a dynamic system by measurements of an output of the system. Classical linear observers can provide optimal estimates of a system state in case of uncertainty modeled by white noise. For nonlinear cases, the theory of nonlinear observers has been developed and its success is mainly due to the sliding mode approach. Using the mathematical theory of variable structure systems with sliding modes, the observer algorithm is designed in such a way that it steers the output of the model to the output of the system obtained via a variety of sensors, in spite of possible mismatches between the assumed model and actual system. The unique properties of sliding mode control allow not only control of the model internal states to the states of the real-life system, but also identification of the disturbance or anomaly that may occur.

Solano, Wanda M.↗

The error associated with use of various forms of the divergence equation to diagnose geopotential and temperature

This study utilizes a model-generated data set to evaluate the errors associated with the diagnosis of geopotential height using three approximate forms of the divergence equation (DE). The DE was solved on the original sigma-coordinate surfaces of the model to exclude any error produced in the process of interpolating the data from observation points to the coordinate surface of the DE. An inviscid, nonlinear, divergent DE produced rms geopotential height errors as large as 3 m and rms temperature errors of over 26 C in the lower planetary boundary layer (PBL), with minimal errors above 800 mb. Elimination of the isobaric divergence in the velocity field caused additional height errors of 5-6 m and temperature errors of 0.4-1.2 C in the free troposphere. Use of a further-degraded DE without nonlinear terms caused height errors of 10-15 m in the upper troposphere with modestly increased temperature errors. Error fields at all levels had synoptic-scale features, with some mesoalpha structure. No spatial noise was apparent. Use of a temporal spectral filter to eliminate high-frequency modes from the model-generated data did not significantly influence the error associated with use of various forms of the DE.

Modica, George M.↗

Method and apparatus for predicting the direction of movement in machine vision

A computer-simulated cortical network is presented. The network is capable of computing the visibility of shifts in the direction of movement. Additionally, the network can compute the following: (1) the magnitude of the position difference between the test and background patterns; (2) localized contrast differences at different spatial scales analyzed by computing temporal gradients of the difference and sum of the outputs of paired even- and odd-symmetric bandpass filters convolved with the input pattern; and (3) the direction of a test pattern moved relative to a textured background. The direction of movement of an object in the field of view of a robotic vision system is detected in accordance with nonlinear Gabor function algorithms. The movement of objects relative to their background is used to infer the 3-dimensional structure and motion of object surfaces.

Lawton, Teri B.↗

Eigenfrequency constrained topology optimization of finite strain hyperelastic structures

This paper incorporates hyperelastic materials, nonlinear kinematics, and preloads in eigenfrequency constrained density–based topology optimization. The formulation allows for initial finite deformations and subsequent small harmonic oscillations. The optimization problem is solved by the method of moving asymptotes, and the gradients are calculated using the adjoint method. Both simple and degenerate eigenfrequencies are considered in the sensitivity analysis. A well-posed topology optimization problem is formulated by filtering the volume fraction field. Numerical issues associated with excessive distortion and spurious eigenmodes in void regions are reduced by removing low volume fraction elements. The optimization objective is to maximize stiffness subject to a lower bound on the fundamental eigenfrequency. Numerical examples show that the eigenfrequencies drastically change with the load magnitude, and that the optimization is able to produce designs with the desired fundamental eigenfrequency.

42 ENGINEERING↗

A Scale‐Dependent Analysis of the Barotropic Vorticity Budget in a Global Ocean Simulation

Abstract The climatological mean barotropic vorticity budget is analyzed to investigate the relative importance of surface wind stress, topography, planetary vorticity advection, and nonlinear advection in dynamical balances in a global ocean simulation. In addition to a pronounced regional variability in vorticity balances, the relative magnitudes of vorticity budget terms strongly depend on the length‐scale of interest. To carry out a length‐scale dependent vorticity analysis in different ocean basins, vorticity budget terms are spatially coarse‐grained. At length‐scales greater than 1,000 km, the dynamics closely follow the Topographic‐Sverdrup balance in which bottom pressure torque, surface wind stress curl and planetary vorticity advection terms are in balance. In contrast, when including all length‐scales resolved by the model, bottom pressure torque and nonlinear advection terms dominate the vorticity budget (Topographic‐Nonlinear balance), which suggests a prominent role of oceanic eddies, which are of km in size, and the associated bottom pressure anomalies in local vorticity balances at length‐scales smaller than 1,000 km. Overall, there is a transition from the Topographic‐Nonlinear regime at scales smaller than 1,000 km to the Topographic‐Sverdrup regime at length‐scales greater than 1,000 km. These dynamical balances hold across all ocean basins; however, interpretations of the dominant vorticity balances depend on the level of spatial filtering or the effective model resolution. On the other hand, the contribution of bottom and lateral friction terms in the barotropic vorticity budget remains small and is significant only near sea‐land boundaries, where bottom stress and horizontal viscous friction generally peak.

Khatri, Hemant↗

Measurement of the Relative and Absolute Microgravity Environment of Space Platforms, Recent Flight Results, and Determination of the Spatial Microgravity Field

Microgravity material processing in space, such as crystal growth, requires specific and accurate knowledge of the very low frequency, very low level microgravity. This means knowledge of absolute microgravity levels. Most acceleration measurement systems for space platforms only measure 'microvibrations,' the measurement of the absolute microgravity level being severely corrupted by bias errors, temperature hysteresis, nonlinearity by-products, and bias changes due to launch vibrations. An advanced system was developed incorporating a fixed 3-axis acceleration measurement system with improved signal conditioning and Invertible Accelerometers for absolute measurements. The higher sensitivity channels have lower frequency low pass filtering to meet crystal growers demands for recording of low level/low frequency disturbances. Also included is a pre-flight bias adjustment feature.

Bijvoet, J. A.↗

Atmospheric model development in support of SEASAT. Volume 3: Prediction models

As part of the SEASAT program of NASA, a set of four hemispheric, atmospheric models are developed. The four models use a polar stereographic grid in the horizontal and a sigma coordinate in the vertical. Conservation forms of the difference equations are integrated using either a fifteen or four minute time step on a 381 km or 127 km grid (at 60 deg N) for the 63 x 63 or 187 x 187 models respectively. A nonlinear pressure smoother is used along with momentum and temperature diffusion to help control computational noise. The horizontal boundary conditions are insulated slippery walls. Centered time differencing with time averaging of the pressure gradient force term in the momentum equations is used. Robert time filtering on of the temperature and moisture solutions is used for computational stability. The moisture and heat source/sink terms are modeled. Terms representing evaporation and large scale condensation, sensible heat exchange, parameterized cumulus convection and precipitation, and solar and terrestrial radiation are included. Initialization of the models is based on a pattern conservation technique to obtain objective analysis of the state parameter structure from the surface to 50 mb.

Wellck, R. E.↗