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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Photometric Calibration of Consumer Video Cameras

Equipment and techniques have been developed to implement a method of photometric calibration of consumer video cameras for imaging of objects that are sufficiently narrow or sufficiently distant to be optically equivalent to point or line sources. Heretofore, it has been difficult to calibrate consumer video cameras, especially in cases of image saturation, because they exhibit nonlinear responses with dynamic ranges much smaller than those of scientific-grade video cameras. The present method not only takes this difficulty in stride but also makes it possible to extend effective dynamic ranges to several powers of ten beyond saturation levels. The method will likely be primarily useful in astronomical photometry. There are also potential commercial applications in medical and industrial imaging of point or line sources in the presence of saturation.This development was prompted by the need to measure brightnesses of debris in amateur video images of the breakup of the Space Shuttle Columbia. The purpose of these measurements is to use the brightness values to estimate relative masses of debris objects. In most of the images, the brightness of the main body of Columbia was found to exceed the dynamic ranges of the cameras. A similar problem arose a few years ago in the analysis of video images of Leonid meteors. The present method is a refined version of the calibration method developed to solve the Leonid calibration problem. In this method, one performs an endto- end calibration of the entire imaging system, including not only the imaging optics and imaging photodetector array but also analog tape recording and playback equipment (if used) and any frame grabber or other analog-to-digital converter (if used). To automatically incorporate the effects of nonlinearity and any other distortions into the calibration, the calibration images are processed in precisely the same manner as are the images of meteors, space-shuttle debris, or other objects that one seeks to analyze. The light source used to generate the calibration images is an artificial variable star comprising a Newtonian collimator illuminated by a light source modulated by a rotating variable neutral- density filter. This source acts as a point source, the brightness of which varies at a known rate. A video camera to be calibrated is aimed at this source. Fixed neutral-density filters are inserted in or removed from the light path as needed to make the video image of the source appear to fluctuate between dark and saturated bright. The resulting video-image data are analyzed by use of custom software that determines the integrated signal in each video frame and determines the system response curve (measured output signal versus input brightness). These determinations constitute the calibration, which is thereafter used in automatic, frame-by-frame processing of the data from the video images to be analyzed.

Suggs, Robert↗
NASA NTRS1997

Large-Eddy Simulation of Transition to Turbulence in Boundary Layers

Large-eddy simulation results for laminar-to-turbulent transition in a spatially developing boundary layer are presented. The disturbances are ingested into a laminar flow through an unsteady suction-and-blowing strip. The filtered, three-dimensional time- dependent Navier-Stokes equations are integrated numerically using spectral, high-order finite-difference, and three-stage low-storage Runge-Kutta methods. The buffer-domain technique is used for the outflow boundary condition. The localized dynamic model used to parameterize the subgrid-scale stresses begins to have a significant impact at the beginning of the nonlinear transition (or intermittency) region. The flow structures commonly found in experiments are also observed in the present simulation; the computed linear instability modes and secondary instability lambda-vortex structures are in agreement with the experiments, and the streak-like-structures and turbulent statistics compare with both the experiments and the theory. The physics captured in the present LES are consistent with the experiments and the full Navier-Stokes simulation (DNS), at a significant fraction of the DNS cost. A comparison of the results obtained with several SGS models shows that the localized model gives accurate results both in a statistical sense and in terms of predicting the dynamics of the energy-carrying eddies, without ad hoc adjustments.

Huai, Xiao-Li↗
DOE OSTI2021

Deep Learning Segmentation of Complex Features in Atomic-Resolution Phase-Contrast Transmission Electron Microscopy Images

Phase-contrast transmission electron microscopy (TEM) is a powerful tool for imaging the local atomic structure of materials. TEM has been used heavily in studies of defect structures of two-dimensional materials such as monolayer graphene due to its high dose efficiency. However, phase-contrast imaging can produce complex nonlinear contrast, even for weakly scattering samples. It is, therefore, difficult to develop fully automated analysis routines for phase-contrast TEM studies using conventional image processing tools. For automated analysis of large sample regions of graphene, one of the key problems is segmentation between the structure of interest and unwanted structures such as surface contaminant layers. In this study, we compare the performance of a conventional Bragg filtering method with a deep learning routine based on the U-Net architecture. We show that the deep learning method is more general, simpler to apply in practice, and produces more accurate and robust results than the conventional algorithm. We provide easily adaptable source code for all results in this paper and discuss potential applications for deep learning in fully automated TEM image analysis.

47 OTHER INSTRUMENTATION↗
DOE OSTI2014

Assessment of volumetric noise and resolution performance for linear and nonlinear CT reconstruction methods

Purpose: For nonlinear iterative image reconstructions (IR), the computed tomography (CT) noise and resolution properties can depend on the specific imaging conditions, such as lesion contrast and image noise level. Therefore, it is imperative to develop a reliable method to measure the noise and resolution properties under clinically relevant conditions. This study aimed to develop a robust methodology to measure the three-dimensional CT noise and resolution properties under such conditions and to provide guidelines to achieve desirable levels of accuracy and precision. Methods: The methodology was developed based on a previously reported CT image quality phantom. In this methodology, CT noise properties are measured in the uniform region of the phantom in terms of a task-based 3D noise-power spectrum (NPS{sub task}). The in-plane resolution properties are measured in terms of the task transfer function (TTF) by applying a radial edge technique to the rod inserts in the phantom. The z-direction resolution properties are measured from a supplemental phantom, also in terms of the TTF. To account for the possible nonlinearity of IR, the NPS{sub task} is measured with respect to the noise magnitude, and the TTF with respect to noise magnitude and edge contrast. To determine the accuracy and precision of the methodology, images of known noise and resolution properties were simulated. The NPS{sub task} and TTF were measured on the simulated images and compared to the truth, with criteria established to achieve NPS{sub task} and TTF measurements with <10% error. To demonstrate the utility of this methodology, measurements were performed on a commercial CT system using five dose levels, two slice thicknesses, and three reconstruction algorithms (filtered backprojection, FBP; iterative reconstruction in imaging space, IRIS; and sinogram affirmed iterative reconstruction with strengths of 5, SAFIRE5). Results: To achieve NPS{sub task} measurements with <10% error, the number of regions of interest needed to be greater than 65. To achieve TTF measurements with <10% error, the contrast-to-noise ratio of the edge needed to be ≥15, achievable by averaging multiple slices across the same edge. The NPS{sub task} measured on a commercial CT system showed IR's reduced noise (IRIS, 30% and SAFIRE5, 55%) and “waxier” texture (peak frequencies: FBP, 0.25 mm{sup −1}; IRIS, 0.23 mm{sup −1}; and SAFIRE5, 0.16 mm{sup −1}). The TTF measured within the axial plane showed improved in-plane resolution with SAFIRE5 at the TTF 50% frequency, f{sub 50} (FBP, 0.36–0.41 mm{sup −1}; SAFIRE5, 0.37–0.46 mm{sup −1}). The TTF measured along the axial direction showed improved z-direction resolution with thinner slice thickness (f{sub 50}: 0.6 mm, 0.35–0.79 mm{sup −1}; 1.5 mm, 0.22–0.3 mm{sup −1}) and with SAFIRE5 (f{sub 50}: FBP, 0.35–0.52 mm{sup −1}; SAFIRE5, 0.42–0.79 mm{sup −1}). Both in-plane and z-direction resolution of SAFIRE5 showed strong dependency on contrast, reflecting SAFIRE5's nonlinearity. Conclusions: A methodology was developed to measure three-dimensional CT noise and resolution properties for iterative reconstruction, especially at challenging measurement conditions with low contrast and high image noise. The methodology also demonstrated its utility for evaluating commercial CT systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Semi-Implicit Computation of Fast Modes in a Scheme Integrating Slow Modes by a Leapfrog Method Based on a Selective Implicit Time Filter

Abstract A scheme for integration of atmospheric equations containing terms with differing time scales is developed. The method employs a filtered leapfrog scheme utilizing a fourth-order implicit time filter with one function evaluation per time step to compute slow-propagating phenomena such as advection and rotation. The terms involving fast-propagating modes are handled implicitly with an unconditionally stable method that permits application of larger time steps and faster computations compared to fully explicit treatment. Implementation using explicit and recurrent formulation is provided. Stability analysis demonstrates that the method is conditionally stable for any combination of frequencies involved in the slow and fast terms as they approach the origin. The implicit filter used in the method damps the computational modes without noticeably sacrificing the accuracy of the physical mode. TheO[(Δt 4 )] accuracy for amplitude errors achieved by the implicitly filtered leapfrog is preserved in applications where terms responsible for fast propagation are integrated with a semi-implicit method. Detailed formulation of the method for soundproof nonhydrostatic anelastic equations is provided. Procedures for implementation in global spectral shallow-water models are also given. Examples comparing numerical and analytical solutions for linear gravity waves demonstrate the accuracy of the scheme. The performance is also shown in more practical nonlinear applications, where numerical solutions accomplished by the method are evaluated against those computed from a scheme where the slow terms are handled by the third-order Runge–Kutta scheme. It demonstrates that the method is able to accurately resolve fine-scale dynamics of Kelvin–Helmholtz shear instabilities, the evolution of density current, and nonlinear drifts of twin tropical cyclones.

Meteorology & Atmospheric Sciences↗
NASA NTRS1973

A guidance and navigation system for continuous low thrust vehicles

A midcourse guidance and navigation system for continuous low thrust vehicles is described. A set of orbit elements, known as the equinoctial elements, are selected as the state variables. The 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 nominal trajectory is defined and the equation of motion and the measurement equation are linearized about this nominal trajectory. An exponential cost criterion is constructed and a linear feedback guidance law is derived to control the thrusting direction of the engine. Using this guidance law, the vehicle will fly in a trajectory neighboring the nominal trajectory. The extended Kalman filter is used for state estimation. Finally a short mission using this system is simulated. The results indicate that this system is very efficient for short missions.

Tse, C. J. C.↗
NASA NTRS2002

Observer-Based Magnetic Bearing Controller Developed for Aerospace Flywheels

A prototype of a versatile, observer-based magnetic bearing controller for aerospace flywheels was successfully developed and demonstrated on a magnetic bearing test rig (see the photograph) and an actual flywheel module. The objective of this development included a fast, yet low risk, control development process, and a robust, high-performance controller for a large variety of flywheels. This required a good system model, an efficient development procedure, and a model-based controller that addressed the key problems associated with flywheel and bearing imbalance, sensor error, and vibration. The model used in this control development and tuning procedure included the flexible rotor dynamics and motor-induced vibrations. Such a model was essential for low-risk scheduling of speed-dependent control parameters and for reliable evaluation of novel control strategies. The successfully tested control prototype utilized an extended Kalman filter to estimate the true rotor principal-axis motion from the raw sensor position feedback. For control refinement, the extended Kalman filter also estimated and eliminated the combined effects of mass-imbalance and sensor runouts from the input data. A key advantage of the design based on the extended Kalman filter is its ability to accurately estimate both the rotor's principal-axis position and gyroscopic rates with the least amount of phase lag. This is important for control parameter scheduling to dampen the gyroscopic motions. Because of large uncertainties in the magnetic bearing and imbalance characteristics, this state-estimation scheme alone is insufficient for containing the rotor motion within the desired 1-mil excursion radius. A nonlinear gain adjustment based on an estimation of the principal-axis orbit size was needed to provide a coarse (nonoptimal), but robust, control of the orbit growth. Control current minimization was achieved with a (steepest gradient) search of synchronous errors in the principal-axis position input data. Actual flywheel tests of this observer-based controller (developed entirely in-house) at the NASA Glenn Research Center showed that the model correctly predicted the rotor orbit growth as a function of rotational speed, and it demonstrated the capability of gain adjustments to arrest this growth. Data from these tests on an actual flywheel module spun to 26,000 rpm proved that the controller was able to contain the shaft motion to within much less than 0.5 mils of radial excursion with axis currents less than 300 mA in root-mean-square estimate. The test speed range was limited because of thermal expansion concerns for this particular flywheel unit, not because of any deficiency in the controller. Simulations for this unit indicated that the controller should be robust up to its top operating speed of 60,000 rpm. Aside from these important achievements, and most significantly, it took less than 1 week to adapt this controller from the simple test rig to the actual flywheel and to demonstrate full five-axis levitation and control. This demonstration showed that both the controller and the model-based development and tuning framework are easily adaptable to a wide range of rotors and bearing configurations and, hence, are capable of reducing design risks and costs for many future flywheel technology developments.

Le, Dzu K.↗
NASA NTRS2013

SDR Input Power Estimation Algorithms

The General Dynamics (GD) S-Band software defined radio (SDR) in the Space Communications and Navigation (SCAN) Testbed on the International Space Station (ISS) provides experimenters an opportunity to develop and demonstrate experimental waveforms in space. The SDR has an analog and a digital automatic gain control (AGC) and the response of the AGCs to changes in SDR input power and temperature was characterized prior to the launch and installation of the SCAN Testbed on the ISS. The AGCs were used to estimate the SDR input power and SNR of the received signal and the characterization results showed a nonlinear response to SDR input power and temperature. In order to estimate the SDR input from the AGCs, three algorithms were developed and implemented on the ground software of the SCAN Testbed. The algorithms include a linear straight line estimator, which used the digital AGC and the temperature to estimate the SDR input power over a narrower section of the SDR input power range. There is a linear adaptive filter algorithm that uses both AGCs and the temperature to estimate the SDR input power over a wide input power range. Finally, an algorithm that uses neural networks was designed to estimate the input power over a wide range. This paper describes the algorithms in detail and their associated performance in estimating the SDR input power.

Nappier, Jennifer M.↗
NASA NTRS1989

Singular spectrum analysis in nonlinear dynamics, with applications to paleoclimatic time series

Two dimensions of a dynamical system given by experimental time series are distinguished. Statistical dimension gives a theoretical upper bound for the minimal number of degrees of freedom required to describe the attractor up to the accuracy of the data, taking into account sampling and noise problems. The dynamical dimension is the intrinsic dimension of the attractor and does not depend on the quality of the data. Singular Spectrum Analysis (SSA) provides estimates of the statistical dimension. SSA also describes the main physical phenomena reflected by the data. It gives adaptive spectral filters associated with the dominant oscillations of the system and clarifies the noise characteristics of the data. SSA is applied to four paleoclimatic records. The principal climatic oscillations and the regime changes in their amplitude are detected. About 10 degrees of freedom are statistically significant in the data. Large noise and insufficient sample length do not allow reliable estimates of the dynamical dimension.

Vautard, R.↗
NASA NTRS2006

Pattern-Recognition Algorithm for Locking Laser Frequency

A computer program serves as part of a feedback control system that locks the frequency of a laser to one of the spectral peaks of cesium atoms in an optical absorption cell. The system analyzes a saturation absorption spectrum to find a target peak and commands a laser-frequency-control circuit to minimize an error signal representing the difference between the laser frequency and the target peak. The program implements an algorithm consisting of the following steps: Acquire a saturation absorption signal while scanning the laser through the frequency range of interest. Condition the signal by use of convolution filtering. Detect peaks. Match the peaks in the signal to a pattern of known spectral peaks by use of a pattern-recognition algorithm. Add missing peaks. Tune the laser to the desired peak and thereafter lock onto this peak. Finding and locking onto the desired peak is a challenging problem, given that the saturation absorption signal includes noise and other spurious signal components; the problem is further complicated by nonlinearity and shifting of the voltage-to-frequency correspondence. The pattern-recognition algorithm, which is based on Hausdorff distance, is what enables the program to meet these challenges.

Karayan, Vahag↗
NASA NTRS1999

Some Continuum Aspects of Data Assimilation

A long-sought goal in data assimilation is to build more fully the time dimension into the data assimilation process. Truly four-dimensional data assimilation requires evolving second-moment information, namely the estimation error covariance, along with the state estimate itself. The continuum evolution equations for the estimation error covariance constitute a system of partial differential equations in six space dimensions, forced by a model error covariance and by the observations themselves. The high dimensionality of this system poses a difficult computational problem. Numerous methods have been proposed for approximate, discrete solution of this system of equations. Four-dimensional variational schemes solve these equations indirectly, while schemes based on Kalman filtering solve the equations more directly, usually on a small subspace of the full six-dimensional space. In both cases, most proposed solution methods are not derived from the continuum covariance evolution equations themselves. Instead, most methods simply inherit the discretization used for evolving the state estimate, either in the form of a discrete tangent linear model or, in the case of ensemble schemes, in the form of a fully nonlinear discrete model. In this lecture we show that solutions of the continuum covariance evolution equations possess simple properties that are not readily satisfied upon discretization of the equations. we give numerical examples illustrating that "inherited" discretizations sometimes fail to produce meaningfully accurate solutions. Finally, we suggest alternative discretization methods that may yield more faithful approximate solutions of the continuum problem.

Cohn, Stephen E.↗
DOE OSTI2016

Stabilization of numerical interchange in spectral-element magnetohydrodynamics

In this study, auxiliary numerical projections of the divergence of flow velocity and vorticity parallel to magnetic field are developed and tested for the purpose of suppressing unphysical interchange instability in magnetohydrodynamic simulations. The numerical instability arises with equal-order C 0 finite- and spectral-element expansions of the flow velocity, magnetic field, and pressure and is sensitive to behavior at the limit of resolution. The auxiliary projections are motivated by physical field-line bending, and coercive responses to the projections are added to the flow-velocity equation. Their incomplete expansions are limited to the highest-order orthogonal polynomial in at least one coordinate of the spectral elements. Cylindrical eigenmode computations show that the projections induce convergence from the stable side with first-order ideal-MHD equations during h-refinement and p-refinement. Hyperbolic and parabolic projections and responses are compared, together with different methods for avoiding magnetic divergence error. Lastly, the projections are also shown to be effective in linear and nonlinear time-dependent computations with the NIMROD code [C. R. Sovinec, et al., J. Comput. Phys. 195 (2004) 355-386], provided that the projections introduce numerical dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Characterization of temporal variations and feedback timescales of exhaust gas recirculation gas properties using high-speed diode laser absorption spectroscopy for next-cycle control of cyclic variability

Dilute combustion offers efficiency gains in boosted gasoline direct injection engines both through knock-limit extension and thermodynamic advantages (i.e. the effect of γ on cycle efficiency), but is limited by cyclic variability at high dilution levels. Past studies have shown that the cycle-to-cycle dynamics are a combination of deterministic and stochastic effects. The deterministic causes of cyclic variations, which arise from feedback due to exhaust gas recirculation, imply the possibility of using active control strategies for dilution limit extension. While internal exhaust gas recirculation will largely provide a next-cycle effect (short-timescale feedback), the feedback of external exhaust gas recirculation will have an effect after a delay of several cycles (long timescale). Therefore, control strategies aiming to improve engine stability at dilution limit may have to account for both short- and long-timescale feedback pathways. This study shows the results of a study examining the extent to which variations in exhaust gas recirculation composition are preserved along the exhaust gas recirculation flow path and thus the relative importance and information content of the long-timescale feedback pathway. To characterize the filtering or retention of cycle-resolved feedback information, high-speed (1–5 kHz) CO 2 concentration measurements were performed simultaneously at three different locations along the low-pressure external exhaust gas recirculation loop of a four-cylinder General Motors gasoline direct injection engine using a multiplexed two-color diode laser absorption spectroscopy sensor system during steady-state and transient engine operation at various exhaust gas recirculation levels. It was determined that cycle-resolved feedback propagates through internal residual gases but is filtered out by the low-pressure exhaust gas recirculation flow system and do not reach the intake manifold. Intermediate variations driven by flow rate and compositional changes are also distinguished and identified.

Dilute combustion↗
NASA NTRS1978

Precision attitude determination for multimission spacecraft

Attitude determination algorithms for a multimission spacecraft are derived and their performance analyzed. The attitude determination system is composed of a strapdown Inertial Reference Unit (IRU), two fixed head star trackers, and an onboard computer. IRU data is processed to maintain real-time knowledge of spacecraft attitude relative to an inertial reference frame. Star tracker data is processed using Kalman filtering techniques to estimate and correct the attitude determination errors and the gyro drift compensation errors. The results of a star availability analysis for stellar, solar and earth pointing missions are presented. Linear covariance analysis techniques are used to evaluate nominal attitude determination performance, the effects of sensor measurement accuracy variations, the effects of errors in knowledge of sensor measurement accuracy, and the effects of star tracker misalignment errors. Results of a nonlinear simulation analysis of attitude determination performance are also presented. These analyses show that precision attitude determination for stellar, solar and earth pointing missions is achieved.

Murrell, J. W.↗
NASA NTRS1990

Estimating short-period dynamics using an extended Kalman filter

An extended Kalman filter is used to estimate the parameters of a low-order model from aircraft transient response data. The low-order model is a state-space model derived from the short-period approximation of the longitudinal aircraft dynamics. The model corresponds to the pitch rate to stick force transfer function currently used in flying qualities analysis. The parameters are estimated from flight data as well as from a 6-DOF nonlinear simulation of the aircraft. These two estimates are then compared, and the discrepancies noted. The low-order model is able to satisfactorily match both flight data and simulation data from a high-order computer simulation.

Bauer, Jeffrey E.↗
NASA NTRS2002

Improved equivalent circuit for twin slot terahertz receivers

Series-fed coplanar waveguide embedding circuits are being developed for terahertz mixers using, in particular, submicron-sized superconducting devices, such as hot electron bolometers as the nonlinear element. Although these mixers show promising performance, they usually also show a considerable downward shift in the center frequency, when compared with simulations obtained by using simplified models. This makes it very difficult to design low-noise mixers for a given THz frequency. This shiftis principally caused by parasitics due to the extremely small details (in terms of wavelength) of the device, and by the electrical properties of the RF choke filter in the DC/IF line. In this paper, we present an improved equivalent network model of such mixer circuits which agrees with measured results at THz frequencies and we propose a new set of THz bolometric mixers that have been fabricated and are currently being tested.

HEB slot antennas CPW↗
DOE OSTI2023

On unifying randomized methods for inverse problems

This work unifies the analysis of various randomized methods for solving linear and nonlinear inverse problems with Gaussian priors by framing the problem in a stochastic optimization setting. By doing so, we show that many randomized methods are variants of a sample average approximation (SAA). More importantly, we are able to prove a single theoretical result that guarantees the asymptotic convergence for a variety of randomized methods. Additionally, viewing randomized methods as an SAA enables us to prove, for the first time, a single non-asymptotic error result that holds for randomized methods under consideration. Another important consequence of our unified framework is that it allows us to discover new randomization methods. Here, we present various numerical results for linear, nonlinear, algebraic, and PDE-constrained inverse problems that verify the theoretical convergence results and provide a discussion on the apparently different convergence rates and the behavior for various randomized methods.

42 ENGINEERING↗
NASA NTRS1973

An all digital phase locked loop for synchronization of a sinusoidal signal embedded in white Gaussian noise

An all digital phase locked loop which tracks the phase of the incoming sinusoidal signal once per carrier cycle is proposed. The different elements and their functions and the phase lock operation are explained in detail. The nonlinear difference equations which govern the operation of the digital loop when the incoming signal is embedded in white Gaussian noise are derived, and a suitable model is specified. The performance of the digital loop is considered for the synchronization of a sinusoidal signal. For this, the noise term is suitably modelled which allows specification of the output probabilities for the two level quantizer in the loop at any given phase error. The loop filter considered increases the probability of proper phase correction. The phase error states in modulo two-pi forms a finite state Markov chain which enables the calculation of steady state probabilities, RMS phase error, transient response and mean time for cycle skipping.

Reddy, C. P.↗