Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “smoothing techniques”

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 307 records · Page 17

Spatial operator approach to flexible multibody system dynamics and control

The inverse and forward dynamics problems for flexible multibody systems were solved using the techniques of spatially recursive Kalman filtering and smoothing. These algorithms are easily developed using a set of identities associated with mass matrix factorization and inversion. These identities are easily derived using the spatial operator algebra developed by the author. Current work is aimed at computational experiments with the described algorithms and at modelling for control design of limber manipulator systems. It is also aimed at handling and manipulation of flexible objects.

Rodriguez, G.↗

Navier-Stokes turbine heat transfer predictions using two-equation turbulence closures

Navier-Stokes calculations were carried out in order to predict the heat-transfer rates on turbine blades. The calculations were performed using TRAF2D which is a k-epsilon, explicit, finite volume mass-averaged Navier-Stokes solver. Turbulence was modeled using Coakley's q-omega and Chien's k-epsilon two-equation models and the Baldwin-Lomax algebraic model. The model equations along with the flow equations were solved explicitly on a nonperiodic C grid. Implicit residual smoothing (IRS) or a combination of multigrid technique and IRS was applied to enhance convergence rates. Calculations were performed to predict the Stanton number distributions on the first stage vane and blade row as well as the second stage vane row of the SSME high-pressure fuel turbine. The comparison serves to highlight the weaknesses of the turbulence models for use in turbomachinery heat-transfer calculations.

Ameri, Ali A.↗

Navier-Stokes turbine heat transfer predictions using two-equation turbulence

Navier-Stokes calculations were carried out in order to predict the heat transfer rates on turbine blades. The calculations were performed using TRAF2D which is a two-dimensional, explicit, finite volume mass-averaged Navier-Stokes solver. Turbulence was modeled using q-omega and k-epsilon two-equation models and the Baldwin-Lomax algebraic model. The model equations along with the flow equations were solved explicitly on a non-periodic C grid. Implicit residual smoothing (IRS) or a combination of multigrid technique and IRS was applied to enhance convergence rates. Calculations were performed to predict the Stanton number distributions on the first stage vane and blade row as well as the second stage vane row of the Rocketdyne Space Shuttle Main Engine (SSME) high pressure fuel turbine. The comparison with the experimental results, although generally favorable, serves to highlight the weaknesses of the turbulence models and the possible areas of improving these models for use in turbomachinery heat transfer calculations.

Ameri, Ali A.↗

Aerodynamic design optimization via reduced Hessian SQP with solution refining

An all-at-once reduced Hessian Successive Quadratic Programming (SQP) scheme has been shown to be efficient for solving aerodynamic design optimization problems with a moderate number of design variables. This paper extends this scheme to allow solution refining. In particular, we introduce a reduced Hessian refining technique that is critical for making a smooth transition of the Hessian information from coarse grids to fine grids. Test results on a nozzle design using quasi-one-dimensional Euler equations show that through solution refining the efficiency and the robustness of the all-at-once reduced Hessian SQP scheme are significantly improved.

Feng, Dan↗

Aerodynamic design optimization via reduced Hessian SQP with solution refining

An all-at-once reduced Hessian Successive Quadratic Programming (SQP) scheme has been shown to be efficient for solving aerodynamic design optimization problems with a moderate number of design variables. This paper extends this scheme to allow solution refining. In particular, we introduce a reduced Hessian refining technique that is critical for making a smooth transition of the Hessian information from coarse grids to fine grids. Test results on a nozzle design using quasi-one-dimensional Euler equations show that through solution refining the efficiency and the robustness of the all-at-once reduced Hessian SQP scheme are significantly improved.

Feng, Dan↗

Low Temperature Growth of Nanostructured Diamond Films on Metals

The field of nanocrystalline diamond and tetrahedral amorphous carbon films has been the focus of intense experimental activity in the last few years for applications in field emission display devices, optical windows, and tribological coatings, The choice of substrate used in most studies has typically been silicon. For metals, however, the thermal expansion mismatch between the diamond film and substrate gives rise to thermal stress that often results in delamination of the film. To avoid this problem in conventional CVD deposition low substrate temperatures (less than 700 C) have been used, often with the incorporation of oxygen or carbon monoxide to the feedgas mixture. Conventionally grown CVD diamond films are also rough and would require post-deposition polishing for most applications. Therefore, there is an obvious need to develop techniques for deposition of well-adhered, smooth nano-structured diamond films on metals for various tribological applications. In our work, nanostructured diamond films are grown on a titanium alloy substrate using a two-step deposition process. The first step is performed at elevated temperature (820 C) for 30 minutes using a H2/CH4/N2 gas mixture in order to grow a thin (approx. 600 nm) nanostructured diamond layer and improve film adhesion. The remainder of the deposition involves growth at low temperature (less than 600 C) in a H2/CH4/O2 gas mixture. Laser reflectance Interferometry (LRI) pattern during growth of a nanostructured diamond film on Ti-6Al-4V alloy. The first 30 minutes are at a high temperature of 820 C and the rest of the film is grown at a low temperature of 580 T. The fringe pattern is observed till the very end due to extremely low surface roughness of 40 nm. The continuation of the smooth nanostructured diamond film growth during low temperature deposition is confirmed by in-situ laser reflectance interferometry and by post-deposition micro-Raman spectroscopy and surface profilometry. Similar experiments performed without the starting nanostructured diamond layer resulted in poorly adhered films with a more crystalline appearance and a higher surface roughness. This low temperature deposition of nanostructured diamond films on metals offers advantages in cases where high residual thermal stress leads to delamination at high temperatures.

Baker, Paul A.↗

Combinations of Earth Orientation Measurements: COMB94 and POLE94

A Kalman filter has been used to combine Earth orientation measurements taken by optical astrometry with SPACE94, a previously determined combination of Earth orientation measurements taken by space-geodetic techniques. This created two series of combined, smoothed, and interpolated polar motion: COMB94 covering 1962-1995 at 5-day intervals, and POLE94 covering 1900 to 1995 at 30.4375-day intervals.

Earth↗

Flex Fuel Optimized SI and HCCI Engine

The central objective of the proposed work is to demonstrate an HCCI (homogeneous charge compression ignition) capable SI (spark ignited) engine that is capable of fast and smooth mode transition between SI and HCCI combustion modes. The model-based control technique was used to develop and validate the proposed control strategy for the fast and smooth combustion mode transition based upon the developed control-oriented engine; and an HCCI capable SI engine was designed and constructed using production ready two-step valve-train with electrical variable valve timing actuating system. Finally, smooth combustion mode transition was demonstrated on a metal engine within eight engine cycles. The Chrysler turbocharged 2.0L I4 direct injection engine was selected as the base engine for the project and the engine was modified to fit the two-step valve with electrical variable valve timing actuating system. To develop the model-based control strategy for stable HCCI combustion and smooth combustion mode transition between SI and HCCI combustion, a control-oriented real-time engine model was developed and implemented into the MSU HIL (hardware-in-the-loop) simulation environment. The developed model was used to study the engine actuating system requirement for the smooth and fast combustion mode transition and to develop the proposed mode transition control strategy. Finally, a single cylinder optical engine was designed and fabricated for studying the HCCI combustion characteristics. Optical engine combustion tests were conducted in both SI and HCCI combustion modes and the test results were used to calibrate the developed control-oriented engine model. Intensive GT-Power simulations were conducted to determine the optimal valve lift (high and low) and the cam phasing range. Delphi was selected to be the supplier for the two-step valve-train and Denso to be the electrical variable valve timing system supplier. A test bench was constructed to develop control strategies for the electrical variable valve timing (VVT) actuating system and satisfactory electrical VVT responses were obtained. Target engine control system was designed and fabricated at MSU for both single-cylinder optical and multi-cylinder metal engines. Finally, the developed control-oriented engine model was successfully implemented into the HIL simulation environment. The Chrysler 2.0L I4 DI engine was modified to fit the two-step vale with electrical variable valve timing actuating system. A used prototype engine was used as the base engine and the cylinder head was modified for the two-step valve with electrical VVT actuating system. Engine validation tests indicated that cylinder #3 has very high blow-by and it cannot be reduced with new pistons and rings. Due to the time constraint, it was decided to convert the four-cylinder engine into a single cylinder engine by blocking both intake and exhaust ports of the unused cylinders. The model-based combustion mode transition control algorithm was developed in the MSU HIL simulation environment and the Simulink based control strategy was implemented into the target engine controller. With both single-cylinder metal engine and control strategy ready, stable HCCI combustion was achived with COV of 2.1% Motoring tests were conducted to validate the actuator transient operations including valve lift, electrical variable valve timing, electronic throttle, multiple spark and injection controls. After the actuator operations were confirmed, 15-cycle smooth combustion mode transition from SI to HCCI combustion was achieved; and fast 8-cycle smooth combustion mode transition followed. With a fast electrical variable valve timing actuator, the number of engine cycles required for mode transition can be reduced down to five. It was also found that the combustion mode transition is sensitive to the charge air and engine coolant temperatures and regulating the corresponding temperatures to the target levels during the combustion mode transition is the key for a smooth combustion mode transition. As a summary, the proposed combustion mode transition strategy using the hybrid combustion mode that starts with the SI combustion and ends with the HCCI combustion was experimentally validated on a metal engine. The proposed model-based control approach made it possible to complete the SI-HCCI combustion mode transition within eight engine cycles utilizing the well controlled hybrid combustion mode. Without intensive control-oriented engine modeling and HIL simulation study of using the hybrid combustion mode during the mode transition, it would be impossible to validate the proposed combustion mode transition strategy in a very short period.

33 ADVANCED PROPULSION SYSTEMS↗

A New and Fast Method for Smoothing Spectral Imaging Data

The Airborne Visible Infrared Imaging Spectrometer (AVIRIS) acquires spectral imaging data covering the 0.4 - 2.5 micron wavelength range in 224 10-nm-wide channels from a NASA ER-2 aircraft at 20 km. More than half of the spectral region is affected by atmospheric gaseous absorption. Over the past decade, several techniques have been used to remove atmospheric effects from AVIRIS data for the derivation of surface reflectance spectra. An operational atmosphere removal algorithm (ATREM), which is based on theoretical modeling of atmospheric absorption and scattering effects, has been developed and updated for deriving surface reflectance spectra from AVIRIS data. Due to small errors in assumed wavelengths and errors in line parameters compiled on the HITRAN database, small spikes (particularly near the centers of the 0.94- and 1.14-micron water vapor bands) are present in this spectrum. Similar small spikes are systematically present in entire ATREM output cubes. These spikes have distracted geologists who are interested in studying surface mineral features. A method based on the "global" fitting of spectra with low order polynomials or other functions for removing these weak spikes has recently been developed by Boardman (this volume). In this paper, we describe another technique, which fits spectra "locally" based on cubic spline smoothing, for quick post processing of ATREM apparent reflectance spectra derived from AVIRIS data. Results from our analysis of AVIRIS data acquired over Cuprite mining district in Nevada in June of 1995 are given. Comparisons between our smoothed spectra and those derived with the empirical line method are presented.

Gao, Bo-Cai↗

ISEE-C attitude determination using fine sun sensor data only

Techniques developed to determine the spin axis attitude using Fine Sun Sensor (FSS) data only are described. At any given time, the Sun angle specifies the orientation of the spin axis relative to the sunline. The instantaneous time rate of change of the sun angle is directly proportional to the orientation of the spin axis relative to a reference plane that is normal to the ecliptic. Thus, the spin axis attitude can be determined when sufficient data has been collected to accurately measure the rate of change of the sun angle. The uncertainties can be computed directly from the uncertainties in the coefficients of the smoothed sun angle curve. The FSS-only technique is unique in that ephemeris vectors are required only to transform the attitude results to more conventional coordinate frames. The combination of the mission geometry and the FSS accuracy make ISEE-C an ideal mission for applying this method. However, the technique can be used on other missions, such as spin stabilized geosynchronous missions.

Gunshol, L. P.↗

Discrete Surface Evolution and Mesh Deformation for Aircraft Icing Applications

Robust, automated mesh generation for problems with deforming geometries, such as ice accreting on aerodynamic surfaces, remains a challenging problem. Here we describe a technique to deform a discrete surface as it evolves due to the accretion of ice. The surface evolution algorithm is based on a smoothed, face-offsetting approach. We also describe a fast algebraic technique to propagate the computed surface deformations into the surrounding volume mesh while maintaining geometric mesh quality. Preliminary results presented here demonstrate the ecacy of the approach for a sphere with a prescribed accretion rate, a rime ice accretion, and a more complex glaze ice accretion.

ice accretion↗

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↗

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. Investigation along these lines has been conducted for highly idealized dynamics such as the strongly nonlinear Lorenz model as well as more realistic models of the means and atmosphere. 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 situations to allow for correct update of the ensemble members. The ensemble technique has also been applied to the problem of smoothing for which similar questions apply. Ensemble smoother examples, however, seem to be quite puzzling in that results state estimates are worse than for their filter analogue. In this study, we use concepts in probability theory to revisit the ensemble methodology for filtering and smoothing in data assimilation. We use the Lorenz model to test and compare the behavior of a variety of 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, Yanqui↗

Some adaptive filtering techniques applied to the passive remote sensing problem

A review is made of recursive statistical regression techniques incorporating past or past and future observations through smoothing and Kalman filtering, respectively; with results for the cases of the Tiros-N/MSU and Nimbus-6/Scams remote sensing satellite experiments. In response to the lack of a satisfactory model for the medium sounded, which is presently a major limitation on retrieval technique performance, a novel, global approach is proposed which casts the retrieval problem into the framework of adaptive filtering. A numerical implementation of such an adaptive system is presented, with a multilayer, semi-spectral general circulation model for the atmosphere being used to fine-tune the sensor as well as the dynamical equations of a Kalman filter. It is shown that the assimilation of radiometric data becomes a straightforward subproblem.

Toldalagi, P. M.↗

Development and evaluation of an inverse solution technique for studying helicopter maneuverability and agility

An inverse solution technique for determining the maximum maneuvering performance of a helicopter using smooth, pilotlike control inputs is presented. Also described is a pilot simulation experiment performed to investigate the accuracy of the solution resulting from this technique. The maneuverability and agility capability of the helicopter math model was varied by varying the pitch and roll damping, the maximum pitch and roll rate, and the maximum load-factor capability. Three maneuvers were investigated: a 180-deg turn, a longitudinal pop-up, and a lateral jink. The inverse solution technique yielded accurate predictions of pilot-in-the-loop maneuvering performance for two of the three maneuvers.

Whalley, Matthew S.↗

An improved multilevel Monte Carlo method for estimating probability distribution functions in stochastic oil reservoir simulations

In this paper, we develop an improved multilevel Monte Carlo (MLMC) method for estimating cumulative distribution functions (CDFs) of a quantity of interest, coming from numerical approximation of large-scale stochastic subsurface simulations. Compared with Monte Carlo (MC) methods, that require a significantly large number of high-fidelity model executions to achieve a prescribed accuracy when computing statistical expectations, MLMC methods were originally proposed to significantly reduce the computational cost with the use of multifidelity approximations. The improved performance of the MLMC methods depends strongly on the decay of the variance of the integrand as the level increases. However, the main challenge in estimating CDFs is that the integrand is a discontinuous indicator function whose variance decays slowly. To address this difficult task, we approximate the integrand using a smoothing function that accelerates the decay of the variance. In addition, we design a novel a posteriori optimization strategy to calibrate the smoothing function, so as to balance the computational gain and the approximation error. The combined proposed techniques are integrated into a very general and practical algorithm that can be applied to a wide range of subsurface problems for high-dimensional uncertainty quantification, such as a fine-grid oil reservoir model considered in this effort. The numerical results reveal that with the use of the calibrated smoothing function, the improved MLMC technique significantly reduces the computational complexity compared to the standard MC approach. Finally, we discuss several factors that affect the performance of the MLMC method and provide guidance for effective and efficient usage in practice.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A two-dimensional Euler solution for an unbladed jet engine configuration

A two dimensional nonaxisymmetric Euler solution in a geometric representation of a jet engine configuration without blades is presented. The domain, including internal and external flow, is covered with a multiblock grid. To construct the grid, a domain decomposition technique is used to subdivide the domain and smooth grids are dimensioned and placed in each block. The grid contains 44 blocks which cover the external field, the inlet, bypass duct, core duct and nozzle of the nonaxisymmetric engine configuration. The geometry is symmetric about the meanline of the hub, but the grid is not since there is no symmetry condition applied to the grid between the two halves. With a symmetric grid at zero angle of attack, the measures of the solution would cancel exactly. With an asymmetric grid, the solution will not necessarily be symmetric and the lift coefficient will not necessarily be zero. Thus, grid asymmetry can be exploited to verify the resolution of the solution. The solution may be verified on the basis of five theoretical quantities: conservation of mass and energy, deviation of the lift coefficient from zero, deviation of the drag coefficient from zero, deviations from constant entropy, and deviations in the pressure distributions over the symmetric surfaces of the components. This technique is suitable for obtaining numerical solutions in complex geometries and provides a foundation for complete engine throughflow calculations.

Stewart, Mark E. M.↗

Post-instability in continuous systems. I - Failure of differentiability of solutions in continuum mechanics

It is pointed out that mathematical models of continua are based on certain assumptions regarding functions which must be at least piece-wise differentiable. The assumption about smoothness of the functions makes it possible to use the mathematical technique of differentiable equations. However, this artificial mathematical limitation follows neither from the principles of mechanics nor from the definition of a continuum. The price paid for such a mathematical convenience is instability (in the class of smooth functions) of the solutions to the corresponding governing equations in some regions of the parameters. A new mathematical technique should, therefore, be developed to describe the solutions which are not necessarily differentiable. The present investigation is concerned with the criteria of applicability of the classical models of continua from the view point of stability of the corresponding solutions, postinstability models derived by reformulation of the original models, and postinstability models in enlarged classes of functions.

Zak, M.↗