Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smoothing 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 991 records · Page 55

Partial and Interaction Spline Models for the Semiparametric Estimation of Functions of Several Variables

A partial spline model is a model for a response as a function of several variables, which is the sum of a smooth function of several variables and a parametric function of the same plus possibly some other variables. Partial spline models in one and several variables, with direct and indirect data, with Gaussian errors and as an extension of GLIM to partially penalized GLIM models are described. Application to the modeling of change of regime in several variables is described. Interaction splines are introduced and described and their potential use for modeling non-linear interactions between variables by semiparametric methods is noted. Reference is made to recent work in efficient computational methods.

Wahba, Grace↗

MEDLI2: MISP Inferred Aerothermal Environment and Flow Transition Assessment

The Mars Entry, Descent, and Landing Instrumentation 2 (MEDLI2) sensor suite on the Mars2020 mission contained multiple sensors on the aeroshell to measure the aerothermal environment during entry into the Martian atmosphere. These sensors performed superbly and successfully returned forebody and aftbody heating measurements. Analysis of MEDLI2 data indicated flow transitioning from a laminar to turbulent state on the heatshield. No evidence of flow transition was observed on the backshell. Two methods were used to estimate flow transition times on the heatshield: (1) temperature gradient of near-surface thermocouple data and (2) heat flux gradient from an inverse reconstruction approach using thermocouple data and material response modeling. Both methods produced similar transition times with an estimated accuracy of ±1 s. To assess various transition criteria, transition parameters were evaluated at each sensor location using flow field solutions from computational fluid dynamics (CFD) simulations. The idea was to use conservative values inferred from MEDLI2 data as transition criteria for other Mars missions. To test this hypothesis, MEDLI data from the Mars Science Laboratory (MSL) mission was used to compare predicted vs. actual flow transition times. The comparisons suggest smooth wall transition criteria are not well-suited in modeling the rapid progression of a turbulent transition front. Transition criteria containing a roughness element parameter agreed better with the flight data. In summary, critical transition values derived from MEDLI2 data may be used as a starting point in constructing a flow transition model for future Mars missions.

Chun Y Tang↗

Development of an adaptive hp-version finite element method for computational optimal control

In this research effort, the usefulness of hp-version finite elements and adaptive solution-refinement techniques in generating numerical solutions to optimal control problems has been investigated. Under NAG-939, a general FORTRAN code was developed which approximated solutions to optimal control problems with control constraints and state constraints. Within that methodology, to get high-order accuracy in solutions, the finite element mesh would have to be refined repeatedly through bisection of the entire mesh in a given phase. In the current research effort, the order of the shape functions in each element has been made a variable, giving more flexibility in error reduction and smoothing. Similarly, individual elements can each be subdivided into many pieces, depending on the local error indicator, while other parts of the mesh remain coarsely discretized. The problem remains to reduce and smooth the error while still keeping computational effort reasonable enough to calculate time histories in a short enough time for on-board applications.

Hodges, Dewey H.↗

River and lake ice conditions as determined from AIRSAR imagery

Synthetic aperture radar (SAR) imagery data can provide information on types and distribution of river and lake ice needed for studying river ice processes and dynamics, monitoring ice during winter navigation, and formulating ice control strategies. Visible and IR remote sensing systems cannot provide such data and present field methods are inadequate for characterizing ice conditions over long river reaches. Our ongoing analysis of JPL's AIRSAR imagery data and concurrent ground truth of ice conditions on the Tanana River and surrounding lakes near Fairbanks, Alaska, in March 1988, has resulted in several findings: hummocked ice covers and zones of variable ice surface roughness within them can be differentiated; C- and L-band data are more sensitive than P-band to the range of surface roughnesses encountered; smooth, level ice that is clear or contains small bubbles produces little backscatter; snow-covered river ice, whether rough or smooth, is distinguishable from snow-covered river sediments on exposed river beds and unvegetated bars; and open water leads are readily distinguished.

Melloh, Rae A.↗

Ultrasonic angle beam standard reflector

A method that provides an impression profile in a reference standard material utilized in inspecting critically stressed components with pulsed ultrasound is described. A die stamp having an I letter is used to impress the surface of a reference material. The die stamp is placed against the surface and struck with an inertia imparting member to impress the I in the reference standard material. Upset may appear on the surface as a result of the impression and is removed to form a smooth surface. The stamping and upset removal is repeated until the entire surface area of a depth control platform on the die stamp uniformly contacts the material surface. The I impression profile in the reference standard material is utilized for reflecting pulsed ultrasonic beams for inspection purposes.

Berry, R. F., Jr.↗

Mixed Linear/Square-Root Encoded Single-Slope Ramp Provides Low-Noise ADC with High Linearity for Focal Plane Arrays

Single-slope analog-to-digital converters (ADCs) are particularly useful for onchip digitization in focal plane arrays (FPAs) because of their inherent monotonicity, relative simplicity, and efficiency for column-parallel applications, but they are comparatively slow. Squareroot encoding can allow the number of code values to be reduced without loss of signal-to-noise ratio (SNR) by keeping the quantization noise just below the signal shot noise. This encoding can be implemented directly by using a quadratic ramp. The reduction in the number of code values can substantially increase the quantization speed. However, in an FPA, the fixed pattern noise (FPN) limits the use of small quantization steps at low signal levels. If the zero-point is adjusted so that the lowest column is onscale, the other columns, including those at the center of the distribution, will be pushed up the ramp where the quantization noise is higher. Additionally, the finite frequency response of the ramp buffer amplifier and the comparator distort the shape of the ramp, so that the effective ramp value at the time the comparator trips differs from the intended value, resulting in errors. Allowing increased settling time decreases the quantization speed, while increasing the bandwidth increases the noise. The FPN problem is solved by breaking the ramp into two portions, with some fraction of the available code values allocated to a linear ramp and the remainder to a quadratic ramp. To avoid large transients, both the value and the slope of the linear and quadratic portions should be equal where they join. The span of the linear portion must cover the minimum offset, but not necessarily the maximum, since the fraction of the pixels above the upper limit will still be correctly quantized, albeit with increased quantization noise. The required linear span, maximum signal and ratio of quantization noise to shot noise at high signal, along with the continuity requirement, determines the number of code values that must be allocated to each portion. The distortion problem is solved by using a lookup table to convert captured code values back to signal levels. The values in this table will be similar to the intended ramp value, but with a correction for the finite bandwidth effects. Continuous-time comparators are used, and their bandwidth is set below the step rate, which smoothes the ramp and reduces the noise. No settling time is needed, as would be the case for clocked comparators, but the low bandwidth enhances the distortion of the non-linear portion. This is corrected by use of a return lookup table, which differs from the one used to generate the ramp. The return lookup table is obtained by calibrating against a stepped precision DC reference. This results in a residual non-linearity well below the quantization noise. This method can also compensate for differential non-linearity (DNL) in the DAC used to generate the ramp. The use of a ramp with a combination of linear and quadratic portions for a single-slope ADC is novel. The number of steps is minimized by keeping the step size just below the photon shot noise. This in turn maximizes the speed of the conversion. High resolution is maintained by keeping small quantization steps at low signals, and noise is minimized by allowing the lowest analog bandwidth, all without increasing the quantization noise. A calibrated return lookup table allows the system to maintain excellent linearity.

Wrigley, Chris J.↗

Optimized Umkehr Profile Algorithm for Ozone Trend Analyses

The long-term record of Umkehr measurements from four NOAA Dobson spectrophotometers was reprocessed after updates to the instrument calibration procedures. In addition, a new data quality-control tool was developed for the Dobson automation software (WinDobson). This paper presents a comparison of Dobson Umkehr ozone profiles from NOAA ozone network stations (Boulder, OHP, MLO, Lauder) against several satellite records, including Aura Microwave Limb Sounder (MLS; ver. 4.2), and combined SBUV and OMPS records (NASA AGG and NOAA COH). A subset of satellite data is selected to match Dobson Umkehr observations at each station spatially (distance less than 200 km) and temporally (within 24 hours). Umkehr Averaging Kernels (AKs) are applied to vertically smooth all overpass satellite profiles prior to comparisons. The station Umkehr record consists of several instrumental records, which have different optical characterizations, and thus instrument-specific stray light contributes to the data processing errors and creates step changes in the record. This work evaluates the overall quality of Umkehr long-term measurements at NOAA ground-based stations and assesses the impact of the instrumental changes on the stability of the Umkehr ozone profile record. This paper describes a method designed to correct biases and discontinuities in the retrieved Umkehr profile that originate from the Dobson calibration process, repair, or optical realignment of the instrument. The M2GMI and GMI CTM ozone profile model output matched to station location and date of observation is used to evaluate instrumental step changes in the Umkehr record. Homogenization of the Umkehr record and discussion of the apparent stray light error in retrieved ozone profiles are the focus of this paper. Homogenization of ground-based records is of great importance for studies of long-term ozone trends and climate change.

Umkehr↗

Supercritical tests of a self-optimizing, variable-Camber wind tunnel model

A testing procedure was used in a 16-foot Transonic Propulsion Wind Tunnel which leads to optimum wing airfoil sections without stopping the tunnel for model changes. Being experimental, the optimum shapes obtained incorporate various three-dimensional and nonlinear viscous and transonic effects not included in analytical optimization methods. The method is a closed-loop, computer-controlled, interactive procedure and employs a Self-Optimizing Flexible Technology wing semispan model that conformally adapts the airfoil section at two spanwise control stations to maximize or minimize various prescribed merit functions subject to both equality and inequality constraints. The model, which employed twelve independent hydraulic actuator systems and flexible skins, was also used for conventional testing. Although six of seven optimizations attempted were at least partially convergent, further improvements in model skin smoothness and hydraulic reliability are required to make the technique fully operational.

Levinsky, E. S.↗

Wind tunnel measurements of three-dimensional wakes of buildings

Measurements relevant to the effect of buildings on the low level atmospheric boundary layer are presented. A wind tunnel experiment was undertaken to determine the nature of the flow downstream from a gap between two transversely aligned, equal sized models of rectangular cross section. These building models were immersed in an equilibrium turbulent boundary layer which was developed on a smooth floor in a zero longitudinal pressure gradient. Measurements with an inclined (45 degree) hot-wire were made at key positions downstream of models arranged with a large, small, and no gap between them. Hot-wire theory is presented which enables computation of the three mean velocity components, U, V and W, as well as Reynolds stresses. These measurements permit understanding of the character of the wake downstream of laterally spaced buildings. Surface streamline patterns obtained by the oil film method were used to delineate the separation region to the rear of the buildings for a variety of spacings.

Logan, E., Jr.↗

A global approach for using kinematic redundancy to minimize base reactions of manipulators

An important consideration in the use of manipulators in microgravity environments is the minimization of the base reactions, i.e. the magnitude of the force and the moment exerted by the manipulator on its base as it performs its tasks. One approach which was proposed and implemented is to use the redundant degree of freedom in a kinematically redundant manipulator to plan manipulator trajectories to minimize base reactions. A global approach was developed for minimizing the magnitude of the base reactions for kinematically redundant manipulators which integrates the Partitioned Jacobian method of redundancy resolution, a 4-3-4 joint-trajectory representation and the minimization of a cost function which is the time-integral of the magnitude of the base reactions. The global approach was also compared with a local approach developed earlier for the case of point-to-point motion of a three degree-of-freedom planar manipulator with one redundant degree-of-freedom. The results show that the global approach is more effective in reducing and smoothing the base force while the local approach is superior in reducing the base moment.

Chung, C. L.↗

A solution adaptive grid procedure for an upwind parabolized flow solver

A solution adaptive grid procedure based on an error equi-distribution scheme is developed and applied to a Parabolized Navier-Stokes solver. An improved method for selecting weighting functions is introduced which involves normalizing a combination of flowfield gradients and curvature of a number of dependent variables and then selecting the largest at each point. The scheme re-distributes grid points line-by-line, with grid point motion controlled by forces analogous to tensional and torsional spring forces with the spring constants set equal to the weighting functions. Torsional terms are functions of the grid point positions along neighboring grid lines and provide grid smoothness and stability. A grid-fitting scheme is introduced for external flows in which the number of grid points in the freestream are reduced to a minimum. Results for several problems are presented to demonstrate the improvements obtainable with the solution adaptive grid procedure.

Harvey, Albert D.↗

AFTI/F-111 airplane mission adaptive wing operational flight evaluation technique using uplinked pilot command cues

NASA and the USAF have conducted a program to investigate aircraft performance improvements utilizing a mission adaptive wing (MAW). The MAW was designed and developed for the AFTI/F-111 variable-sweep aircraft to provide a hydraulically driven, smooth, and continuous variable camber of the trailing and leading edges as a function of maneuvering requirements or of flight conditions. The remotely augmented vehicle facility (RAV) at the NASA DFRF, as utilized in the MAW investigations, is described. The RAV was a dedicated, ground based, general purpose facility capable of receiving a data stream downlinked from a test vehicle, processing this data stream in a digital computer, and transmitting processed data back to the test vehicle. It is shown that this method of flight testing provides a technique that can evaluate highly dynamic maneuvers.

Kempel, Robert W.↗

Methods for extracting aerodynamic accelerations from Orbiter High Resolution Accelerometer Package flight data

This report presents a process for extracting the aerodynamic accelerations of the Shuttle Orbiter Vehicle from the High Resolution Accelerometer Package (HiRAP) flight data during reentry. The methods for obtaining low-level aerodynamic accelerations, principally in the rarefied flow regime, are applied to 10 Orbiter flights. The extraction process is presented using data obtained from Space Transportation System Flight 32 (Mission 61-C) as a typical example. This process involves correcting the HiRAP measurements for the effects of temperature bias and instrument offset from the Orbiter center of gravity, and removing acceleration data during times they are affected by thruster firings. The corrected data are then made continuous and smooth and are further enhanced by refining the temperature bias correction and removing effects of the auxiliary power unit actuation. The resulting data are the current best estimate of the Orbiter aerodynamic accelerations during reentry and will be used for further analyses of the Orbiter aerodynamics and the upper atmosphere characteristics.

Thompson, J. M.↗

An early solar dynamo prediction: Cycle 23 is approximately cycle 22

In this paper, we briefly review the 'dynamo' and 'geomagnetic precursor' methods of long-term solar activity forecasting. These methods depend upon the most basic aspect of dynamo theory to predict future activity, future magnetic field arises directly from the magnification of pre-existing magnetic field. We then generalize the dynamo technique, allowing the method to be used at any phase of the solar cycle, through the development of the 'Solar Dynamo Amplitude' (SODA) index. This index is sensitive to the magnetic flux trapped within the Sun's convection zone but insensitive to the phase of the solar cycle. Since magnetic fields inside the Sun can become buoyant, one may think of the acronym SODA as describing the amount of buoyant flux. Using the present value of the SODA index, we estimate that the next cycle's smoothed peak activity will be about 210 +/- 30 solar flux units for the 10.7 cm radio flux and a sunspot number of 170 +/- 25. This suggests that solar cycle #23 will be large, comparable to cycle #22. The estimated peak is expected to occur near 1999.7 +/- 1 year. Since the current approach is novel (using data prior to solar minimum), these estimates may improve when the upcoming solar minimum is reached.

Schatten, Kenneth H.↗

The dependence of wind stress on wave height and wind speed

Three near-neutral boundary layer data sets were investigated with the aim of finding a dependence of wind stress on both wind speed and significant wave height. The data set most representative of open-ocean wave height, wind speed, and momentum flux conditions, was selected and analyzed by means of the least-squares method to produce a new parameterization for the wind stress as a function of both wind speed and significant wave height. This study shows that the wind stress, and consequently the drag coefficient, decreases with increasing wave height for a fixed wind speed. The study also shows that the curvature of the wind profile decreases with increasing wave height and that the C(DN) = A + BU-bar(10) form for the drag coefficient parameterization is inadequate. A drag coefficient that applies to both smooth and rough flows is proposed. These results are more applicable for open-ocean deep-water conditions and less applicable for sheltered, closed, shallow water sites.

Blake, Reginald A.↗

Theory for the effects of turbulence in a planetary atmosphere on radio occultation

Rytov's method is used to formulate the correlation functions for log-amplitude and phase fluctuations for both spherical and plane wave propagation in a turbulent medium whose correlation function for refractive index fluctuations is described by the product of a function of the average coordinate and a function of the difference coordinate. The results are applied to the study of radio occultation effects due to the atmosphere of Venus in the case of a flyby space probe. It is assumed that turbulence of isotropic and smoothly varying characteristics occurs in the Venusian atmosphere.

Woo, R.↗

Gravitational spectra from direct measurements

A simple rapid method is described for determining the spectrum of a surface field from harmonic analysis of direct measurements along great circle arcs. The method is shown to give excellent overall trends to very high degree from even a few short arcs of satellite data. Three examples are taken with perfect measurements of satellite tracking over a planet made up of hundreds of point-masses using (1) altimetric heights from a low orbiting spacecraft, (2) velocity residuals between a low and a high satellite in circular orbits, and (3) range-rate data between a station at infinity and a satellite in highly eccentric orbit. In particular, the smoothed spectrum of the Earth's gravitational field is determined to about degree 400(50 km half wavelength) from 1 D x 1 D gravimetry and the equivalent of 11 revolutions of Geos 3 and Skylab altimetry. This measurement shows there is about 46 cm of geoid height remaining in the field beyond degree 180.

Wagner, C. A.↗

Increasing Durability of Flame-Sprayed Strain Gauges

Thermally sprayed dielectric ceramic coatings are the primary means of attaching strain and temperature gauges to hot-section rotating parts of turbine engines. As hot-section temperatures increase, lifetimes of installed gauges decrease, and seldom exceed one hour above 2,000 F (approx.1,100 C). Advanced engine components are expected to operate at temperatures approaching 2,200 F (approx.1,200 C), and the required high-temperature lifetime is 10 hours minimum. Typically, to enable a ceramic coating to adhere to the smooth surface of an engine component, a thermally sprayed NiCrAlY or NiCoCrAlY bond coat is applied to the smooth surface, thereby providing a textured surface to which the ceramic coat can adhere. The main failure mechanism of this system is decohesion and/or delamination at the interface between the ceramic top coat and the bond coat, caused by oxidation of the bond coat and stresses from the mismatch between the coefficients of thermal expansion of the ceramic top coat and the metallic bond coat. The approach taken to increase the high-temperature lifetime of a gauge attached to an engine component by the method described above involves (1) selective oxidation of the bond coat by means of a heat treatment in reduced oxygen partial pressure followed by (2) the application of a noble-metal diffusion barrier. In experiments to test this approach, heat treatments of NiCoCrAlY bond coats were carried out in a tube furnace in which, in each case, the temperature was alternately (1) increased at a rate of 3 C per minute and (2) held steady for one hour until the desired temperature was reached. The tube furnace was continuously purged with dry nitrogen gas. A final heat-treatment temperature range of 1,600 to 1,800 F (871 to 982 C) proved most beneficial.

Gregory, Otto J.↗