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 739 records · Page 41

Aircraft automatic-flight-control system with inversion of the model in the feed-forward path using a Newton-Raphson technique for the inversion

A new automatic flight control system concept suitable for aircraft with highly nonlinear aerodynamic and propulsion characteristics and which must operate over a wide flight envelope was investigated. This exact model follower inverts a complete nonlinear model of the aircraft as part of the feed-forward path. The inversion is accomplished by a Newton-Raphson trim of the model at each digital computer cycle time of 0.05 seconds. The combination of the inverse model and the actual aircraft in the feed-forward path alloys the translational and rotational regulators in the feedback path to be easily designed by linear methods. An explanation of the model inversion procedure is presented. An extensive set of simulation data for essentially the full flight envelope for a vertical attitude takeoff and landing aircraft (VATOL) is presented. These data demonstrate the successful, smooth, and precise control that can be achieved with this concept. The trajectory includes conventional flight from 200 to 900 ft/sec with path accelerations and decelerations, altitude changes of over 6000 ft and 2g and 3g turns. Vertical attitude maneuvering as a tail sitter along all axes is demonstrated. A transition trajectory from 200 ft/sec in conventional flight to stationary hover in the vertical attitude includes satisfactory operation through lift-cure slope reversal as attitude goes from horizontal to vertical at constant altitude. A vertical attitude takeoff from stationary hover to conventional flight is also demonstrated.

Smith, G. A.↗

High resolution upwind schemes for the three-dimensional incompressible Navier-Stokes equations

Based on flux-difference splitting, implicit high resolution schemes are constructed for efficient computations of steady-state solutions to the three-dimensional, incompressible Navier-Stokes equations in curvilinear coordinates. These schemes use first-order accurate Euler backward-time differencing and second-order central differencing for the viscous shear fluxes. Up to third-order accurate upwind differencing is achieved through a reconstruction of the solution from its cell averages. The reconstruction is accomplished by linear interpolation, where the node stencils are selected such that in regions of smooth solution the flow is highly resolved while spurious oscillations in regions of rapid changes in gradient are still suppressed. Fairly rapid convergence to steady-state solutions is attained with a completely vectorizable hybrid time-marching method. Flows around a sharp-edged delta wing are computed with the maximum accuracy of the upwind-differencing restricted to first-, second-, and third-order, to illustrate the effect of accuracy on the global and on the local vortical flow fields. The results are validated with experimental data.

Hartwich, PETER-M.↗

Photometric geodesy of main-belt asteroids. II - Analysis of lightcurves for poles, periods, and shapes

The assumption that asteroids can be modeled as smooth, featureless, triaxial ellipsoids that rotate about their shortest axes is presently used to study all but one of the 26 asteroids treated in the Weidenschilling et al. (1987) 'photometric geodesy' program. Rotational poles derived from three independent methods are used to determine each asteroid's sidereal period and triaxial ellipsoid axial ratios, together with their associated photometric parameters. The asteroids appear to have rotational poles that do not lie along their orbital planes.

Drummond, J. D.↗

On constructing three-dimensional overlapping grids with CMPGRD

Techniques for the construction of three-dimensional composite overlapping grids, using the grid construction program CMPGRD, are described. The overlapping approach can be used to generate grids for regions of complicated geometry. The grids can be constructed to be smooth and free from coordinate singularities. The ability to create smooth grids for complicated regions is an important first step towards the accurate numerical solution of partial differential equations. The creation of grids for surfaces defined by cross-sections such as an airplane wing is described. A method for integrating the patched surfaces generated by a computer aided design (CAD) package with the CMPGRD program is described. Additionally, the creation of grids in regions where surfaces intersect is described.

Henshaw, William D.↗

A Localized Tau Method PDE Solver

In this paper we present a new form of the collocation method that allows one to find very accurate solutions to time marching problems without the unwelcome appearance of Gibb's phenomenon oscillations. The basic method is applicable to any partial differential equation whose solution is a continuous, albeit possibly rapidly varying function. Discontinuous functions are dealt with by replacing the function in a small neighborhood of the discontinuity with a spline that smoothly connects the function segments on either side of the discontinuity. This will be demonstrated when the solution to the inviscid Burgers equation is discussed.

Cottam, Russell↗

Orion MPCV Touchdown Detection Threshold Development and Testing

A robust method of detecting Orion Multi ]Purpose Crew Vehicle (MPCV) splashdown is necessary to ensure crew and hardware safety during descent and after touchdown. The proposed method uses a triple redundant system to inhibit Reaction Control System (RCS) thruster firings, detach parachute risers from the vehicle, and transition to the post ]landing segment of the Flight Software (FSW). The vehicle crew is the prime input for touchdown detection, followed by an autonomous FSW algorithm, and finally a strictly time based backup timer. RCS thrusters must be inhibited before submersion in water to protect against possible damage due to firing these jets under water. In addition, neglecting to declare touchdown will not allow the vehicle to transition to post ]landing activities such as activating the Crew Module Up ]righting System (CMUS), resulting in possible loss of communication and difficult recovery. A previous AIAA paper gAssessment of an Automated Touchdown Detection Algorithm for the Orion Crew Module h concluded that a strictly Inertial Measurement Unit (IMU) based detection method using an acceleration spike algorithm had the highest safety margins and shortest detection times of other methods considered. That study utilized finite element simulations of vehicle splashdown, generated by LS ]DYNA, which were expanded to a larger set of results using a Kriging surface fit. The study also used the Decelerator Systems Simulation (DSS) to generate flight dynamics during vehicle descent under parachutes. Proto ]type IMU and FSW MATLAB models provided the basis for initial algorithm development and testing. This paper documents an in ]depth trade study, using the same dynamics data and MATLAB simulations as the earlier work, to further develop the acceleration detection method. By studying the combined effects of data rate, filtering on the rotational acceleration correction, data persistence limits and values of acceleration thresholds, an optimal configuration was determined. The lever arm calculation, which removes the centripetal acceleration caused by vehicle rotation, requires that the vehicle angular acceleration be derived from vehicle body rates, necessitating the addition of a 2nd order filter to smooth the data. It was determined that using 200 Hz data directly from the vehicle IMU outperforms the 40 Hz FSW data rate. Data persistence counter values and acceleration thresholds were balanced in order to meet desired safety and performance. The algorithm proved to exhibit ample safety margin against early detection while under parachutes, and adequate performance upon vehicle splashdown. Fall times from algorithm initiation were also studied, and a backup timer length was chosen to provide a large safety margin, yet still trigger detection before CMUS inflation. This timer serves as a backup to the primary acceleration detection method. Additionally, these parameters were tested for safety on actual flight test data, demonstrating expected safety margins.

Daum, Jared↗

Determination of discontinuities in the derivative of a piecewise-smooth function on the basis of the results of measurements

A piecewise-smooth function with discontinuity in the first derivative on a given interval is considered. The values of the function are measured at a sequence of points in the interval and a random error is included in the measurements. A method is proposed to estimate the position of the discontinuity in the derivative. Regression lines are associated with each measurement point and account for k - 1 points preceding or following the point. The estimate for the position of the discontinuity is the measurement point with the largest angle between the regression lines. The error in the estimate is analyzed and the results are verified.

Gurin, L. S.↗

Carbide coated fibers in graphite-aluminum composites

The study of protective-coupling layers of refractory metal carbides on the graphite fibers prior to their incorporation into composites is presented. Such layers should be directly wettable by liquid aluminum and should act as diffusion barriers to prevent the formation of aluminum carbide. Chemical vapor deposition was used to uniformly deposit thin, smooth, continuous coats of ZrC on the carbon fibers of tows derived from both rayon and polyacrylonitrile. A wet chemical coating of the fibers, followed by high-temperature treatment, was used, and showed promise as an alternative coating method. Experiments were performed to demonstrate the ability of aluminum alloys to wet carbide surfaces. Titanium carbide, zirconium carbide and carbide-coated graphite surfaces were successfully wetted. Results indicate that initial attempts to wet surfaces of ZrC-coated carbon fibers appear successful.

Imprescia, R. J.↗

NASA research related to sailplane airfoils

The theoretical methods and experimental facilities at the NASA Langley Research Center have been employed to conduct investigations of sailplane airfoils. The unique and powerful capabilities of the Eppler Program have been used to design and analyze many airfoils and to smooth several Wortmann airfoils. Wind-tunnel investigations of two sailplane airfoils have been conducted in the Langley low-turbulence pressure tunnel. A procedure for sailplane performance improvement has been outlined.

Somers, D. M.↗

Updating Landsat-derived land-cover maps using change detection and masking techniques

The California Integrated Remote Sensing System's San Bernardino County Project was devised to study the utilization of a data base at a number of jurisdictional levels. The present paper discusses the implementation of change-detection and masking techniques in the updating of Landsat-derived land-cover maps. A baseline landcover classification was first created from a 1976 image, then the adjusted 1976 image was compared with a 1979 scene by the techniques of (1) multidate image classification, (2) difference image-distribution tails thresholding, (3) difference image classification, and (4) multi-dimensional chi-square analysis of a difference image. The union of the results of methods 1, 3 and 4 was used to create a mask of possible change areas between 1976 and 1979, which served to limit analysis of the update image and reduce comparison errors in unchanged areas. The techniques of spatial smoothing of change-detection products, and of combining results of difference change-detection algorithms are also shown to improve Landsat change-detection accuracies.

Likens, W.↗

High-order ENO schemes applied to two- and three-dimensional compressible flow

High order essentially non-oscillatory (ENO) finite difference schemes are applied to the 2-D and 3-D compressible Euler and Navier-Stokes equations. Practical issues, such as vectorization, efficiency of coding, cost comparison with other numerical methods, and accuracy degeneracy effects, are discussed. Numerical examples are provided which are representative of computational problems of current interest in transition and turbulence physics. These require both nonoscillatory shock capturing and high resolution for detailed structures in the smooth regions and demonstrate the advantage of ENO schemes.

Shu, Chi-Wang↗

Total ozone and aerosol optical depths inferred from radiometric measurements in the Chappuis absorption band

A second-derivative smoothing technique, commonly used in inversion work, is applied to the problem of inferring total columnar ozone amounts and aerosol optical depths. The application is unique in that the unknowns may be solved for directly without employing standard inversion methods. It is shown, however, that by employing inversion constraints, better solutions are normally obtained. The method is quite versatile and able to deal with varying total ozone and various aerosol size distributions. The technique is applied first in simulation, then to 119 days of measurements taken in Tucson, Arizona, that are compared to TOMS values for the same dates. The technique is also applied to two measurements taken at Mauna Loa, Hawaii, for which Dobson ozone values are available in addition to the TOMS values, and the results agree to within 15 percent. It is also shown through simulations that additional information can be obtained from measurements outside the Chappuis band. This approach reduces the bias and spread of the estimated total ozone and is unique in that it uses measurements from both the Chappuis and Huggins absorption bands.

Flittner, D. E.↗

Phased-Array Measurements of Single Flow Hot Jets

A 16 microphone phased-array system has been successfully applied to measure jet noise source distributions. In this study, a round convergent nozzle was tested at various hot and cold flow conditions: acoustic Mach numbers are between 0.35 and 1.6 and static temperature ratios are varied from cold to 2.7. The classical beamforming method was applied on narrowband frequencies. From the measured source distributions locations of peak strength were tracked and found to be very consistent between adjacent narrowband frequencies. In low speed heated and unheated jets, the peak source locations vary smoothly from the nozzle exit to downstream as the frequency is decreased. When the static temperature ratio was kept constant, the peak source position moved downstream with increasing acoustic Mach number for the Strouhal numbers smaller than about 1.5. It was also noted that the peak source locations of low frequencies occur farther downstream than the end of potential core.

Bridges, James↗

Algorithm for NO2 Vertical Column Retrieval from the Ozone Monitoring Instrument

We describe the operational algorithm for the retrieval of stratospheric, tropospheric, and total column densities of nitrogen dioxide NO2 from earthshine radiances measured by the Ozone Monitoring Instrument (OMI), aboard the EOS-Aura satellite. The algorithm uses the DOAS method for the retrieval of slant column NO densities. Air mass factors (AMFs) calculated from a stratospheric NO2 profile are used to make initial estimates of the vertical column density. Using data collected over a 24-h period, a smooth estimate of the global stratospheric field is constructed. Where the initial vertical column densities exceed the estimated stratospheric field, we infer the presence of tropospheric NO2, and recalculate the vertical column density (VCD) using an AMF calculated from an assumed tropospheric NO2 profile. The parameters that control the operational algorithm were selected with the aid of a set of data assembled from stratospheric and tropospheric chemical transport models. We apply the optimized algorithm to OMI data and present global maps of NO2 VCDs for the first time.

Bucsela, Eric J.↗

Incremental States for Precise On-Orbit Relative Knowledge in Formation Flight

High precision and close proximity formation flight is an enabling technology for future space missions and requires an on-board relative navigation capability that is accurate to the mm-level and robust to formation parameters. Common estimation techniques linearize the entire formation about one spacecraft’s position, resulting in degraded accuracy due to linearization errors when separation distances become large. Additionally, methods which decouple absolute and relative state estimates usually require ad-hoc methods to incorporate the estimates together. This work discusses an alternative ”incremental” formulation of the relative navigation problem which is in variant to formation size, robust to coupling between absolute and relative dynamics, and can undergo similarity transformations to smoothly incorporate either absolute or relative information without numerical issues. A specific example of this architecture is presented in the context of formation navigation using Carrier-Differential Global Positioning System measurements, and is compared to a traditional leader-linearized filter. In the presence of accurate measurements, the incremental architecture is shown to reduce linearization errors and mean estimate errors by up to three orders of magnitude without the incorporation of new sensor information.

Seubert, Carl↗

Recursive inverse kinematics for robot arms via Kalman filtering and Bryson-Frazier smoothing

This paper applies linear filtering and smoothing theory to solve recursively the inverse kinematics problem for serial multilink manipulators. This problem is to find a set of joint angles that achieve a prescribed tip position and/or orientation. A widely applicable numerical search solution is presented. The approach finds the minimum of a generalized distance between the desired and the actual manipulator tip position and/or orientation. Both a first-order steepest-descent gradient search and a second-order Newton-Raphson search are developed. The optimal relaxation factor required for the steepest descent method is computed recursively using an outward/inward procedure similar to those used typically for recursive inverse dynamics calculations. The second-order search requires evaluation of a gradient and an approximate Hessian. A Gauss-Markov approach is used to approximate the Hessian matrix in terms of products of first-order derivatives. This matrix is inverted recursively using a two-stage process of inward Kalman filtering followed by outward smoothing. This two-stage process is analogous to that recently developed by the author to solve by means of spatial filtering and smoothing the forward dynamics problem for serial manipulators.

Rodriguez, G.↗

An effective multigrid method for high-speed flows

The use is considered of a multigrid method with central differencing to solve the Navier-Stokes equations for high speed flows. The time dependent form of the equations is integrated with a Runge-Kutta scheme accelerated by local time stepping and variable coefficient implicit residual smoothing. Of particular importance are the details of the numerical dissipation formulation, especially the switch between the second and fourth difference terms. Solutions are given for 2-D laminar flow over a circular cylinder and a 15 deg compression ramp.

Swanson, R. C.↗