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 595 records · Page 33

Aliased tidal errors in TOPEX/POSEIDON sea surface height data

Alias periods and wavelengths for the M(sub 2, S(sub 2), N(sub 2), K(sub 1), O(sub 1), and P(sub 1) tidal constituents are calculated for TOPEX/POSEIDON. Alias wavelenghts calculated in previous studies are shown to be in error, and a correct method is presented. With the exception of the K(sub 1) constituent, all of these tidal aliases for TOPEX/POSEIDON have periods shorter than 90 days and are likely to be confounded with long-period sea surface height signals associated with real ocean processes. In particular, the correspondence between the periods and wavelengths of the M(sub 2) alias and annual baroclinic Rossby waves that plagued Geosat sea surface height data is avoided. The potential for aliasing residual tidal errors in smoothed estimates of sea surface height is calculated for the six tidal constituents. The potential for aliasing the lunar tidal constituents M(sub 2), N(sub 2) and O(sub 1) fluctuates with latitude and is different for estimates made at the crossovers of ascending and descending ground tracks than for estimates at points midway between crossovers. The potential for aliasing the solar tidal constituents S(sub 2), K(sub 1) and P(sub 1) varies smoothly with latitude. S(sub 2) is strongly aliased for latitudes within 50 degress of the equator, while K(sub 1) and P(sub 1) are only weakly aliased in that range. A weighted least squares method for estimating and removing residual tidal errors from TOPEX/POSEIDON sea surface height data is presented. A clear understanding of the nature of aliased tidal error in TOPEX/POSEIDON data aids the unambiguous identification of real propagating sea surface height signals. Unequivocal evidence of annual period, westward propagating waves in the North Atlantic is presented.

Schlax, Michael G.↗

Equivalent absorption coefficients generated from frequency probability distributions

A flexible and computationally accurate method of treating aerosol scattering in spectral regions in which gaseous absorption is important is described. In the method, line-by-line absorption coefficients are computed as a function of pressure, temperature, and absorber gas for the spectral region of interest. The coefficients are sorted into a probability distribution which is converted into a cumulative probability distribution, which in turn can be inverted due to its monotonic nature. The inverted distribution is a smooth curve giving the absorption coefficient as a function of an independent variable on the domain. The frequency integration of the radiative transfer equation can then be performed by a quadrature technique with values of the absorption coefficient determined from the inverted distribution curve. The method is illustrated by applying it to the 9.6 micron band of ozone.

Ackerman, T. P.↗

Analysis of high-incidence separated flow past airfoils

An unsteady Navier-Stokes (NS) analysis is developed and used to carefully examine high-incidence aerodynamic separated flows past airfoils. Clustered conformal C-grids are employed for the 12 percent thick symmetric Joukowski airfoil as well as for the NACA 0012 airfoil with a sharp trailing edge. The clustering is controlled by appropriate one-dimensional stretching transformations. An attempt is made to resolve many of the dominant scales of an unsteady flow with massive separation, while maintaining the transformation metrics to be smooth and continuous in the entire flow field. A fully implicit time-marching alternating-direction implicit-block Gaussian elimination (ADI-BGE) method is employed, in which no use is made of any explicit artificial dissipation. Detailed results are obtained for massively separated, unsteady flow past symmetric Joukowski and NACA 0012 airfoils.

Chia, K. N.↗

Applications of a Conservative Zonal Scheme to Transient and Geometrically Complex Problems

A conservative zoning technique, wherein the flow field for a finite-difference calculation is divided into several regions to simplify grid generation, is discussed and is applied in the solution of a two-dimensional problem of complex topology. Calculations are performed on two zonal, or patched, grid systems for the supersonic flow over a double-airfoil configuration. The solution is smooth and continuous across the zonal interfaces, and shock waves pass through the boundaries without distortion. In addition, the time accuracy of the zonal boundary method is verified by a two-zone cylinder calculation with a stationary inner and a rotating outer mesh. The feasibility of the zonal approach for use in the solution of geometrically complex and unsteady problems is thus demonstrated.

Hessenius, Kristin A.↗

Entropy Stable Method for the Euler Equations Revisited: Central Differencing via Entropy Splitting and SBP

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Ols- son & Oliger (1994), Gerritsen & Olsson (1996), and Yee et al. (2000) is revisited. The entropy splitting is a form of skew-symmetric splitting of the nonlinear Euler flux derivatives. Central differencing applied to the entropy splitting form of the Euler flux derivatives together with SBP difference operators will, hereafter, be referred to as entropy split schemes. This study is prompted by the recent growing interest in numerical methods for which a discrete entropy conservation law holds, a discrete global entropy conservation can be proved and/or the numerical method possesses a stable entropy in the framework of SBP difference operators and L2-energy norm estimate. The objective of this paper is to recast the entropy split scheme as the re- cent definition of an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for non- linear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grub- ber splitting.

Sjogreen, Bjorn↗

Entropy Stable Method for the Euler Equations Revisited: Central Differencing via Entropy Splitting and SBP

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Olsson & Oliger, Gerritsen & Olsson, and Yee et al. (15, 7, 37) is revisited. The objective of this paper is to prove for the first time that the entropy split scheme is an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for nonlinear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. The proof is to replace the spatial derivatives by summation-by-parts (SBP) difference operators in the entropy split form of the equations using the physical entropy of the Euler equations. The numerical boundary closure follows directly from the SBP operator. No additional numerical boundary procedure is required. In contrast, Tadmor-type entropy conserving schemes (31) using mathematical entropies and more recently in (35], do not naturally come with a numerical boundary closure and a generalized SBP operator has to be developed (18). Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grubber splitting.

Sjogreen, Bjorn↗

Development of Numerical Codes for Modeling Electromagnetic Behavior at High Frequencies Near Large Objects

A study into the problem of determining electromagnetic solutions at high frequencies for problems involving complex geometries, large sizes and multiple sources (e.g. antennas) has been initiated. Typical applications include the behavior of antennas (and radiators) installed on complex conducting structures (e.g. ships, aircrafts, etc..) with strong interactions between antennas, the radiation patterns, and electromagnetic signals is of great interest for electromagnetic compatibility control. This includes the overall performance evaluation and control of all on-board radiating systems, electromagnetic interference, and personnel radiation hazards. Electromagnetic computational capability exists at NASA LaRC, and many of the codes developed are based on the Moment Method (MM). However, the MM is computationally intensive, and this places a limit on the size of objects and structures that can be modeled. Here, two approaches are proposed: (i) a current-based hybrid scheme that combines the MM with Physical optics, and (ii) an Alternating Direction Implicit-Finite Difference Time Domain (ADI-FDTD) method. The essence of a hybrid technique is to split the overall scattering surface(s) into two regions: (a) a MM zone (MMZ) which can be used over any part of the given geometry, but is most essential over irregular and "non-smooth" geometries, and (b) a PO sub-region (POSR). Currents induced on the scattering and reflecting surfaces can then be computed in two ways depending on whether the region belonged to the MMZ or was part of the POSR. For the MMZ, the current calculations proceed in terms of basis functions with undetermined coefficients (as in the usual MM method), and the answer obtained by solving a system of linear equations. Over the POSR, conduction is obtained as a superposition of two contributions: (i) currents due to the incident magnetic field, and (ii) currents produced by the mutual induction from conduction within the MMZ. This effectively leads to a reduction in the size of linear equations from N to N - Npo with N being the total number of segments for the entire surface and Npo the number of segments over the POSR. The scheme would be appropriate for relatively large, flat surfaces, and at high frequencies. The ADI-FDTD scheme provides for both transient and steady state analyses. The restrictive Courant-Friedrich-Levy (CFL) condition on the time-step is removed, and so large time steps can be chosen even though the spatial grids are small. This report includes the problem definition, a detailed discussion of both the numerical techniques, and numerical implementations for simple surface geometries. Numerical solutions have been derived for a few simple situations.

Joshi, R. P.↗

Acoustic Liner Drag: Measurement Uncertainty Reduction and Application to Novel Perforate Geometries

In addition to developing acoustic liner concepts and characterizing their performance, the NASA Langley Liner Physics Team has investigated issues related to liner drag. A method to quantify relative drag of liner configurations was developed and employed to understand the effects of various liner features. It was observed that perforate shape could have a marked effect on the measured liner drag and led to the discovery of a low-drag geometry that cut the drag penalty between a perforate and a smooth wall by approximately 50%. Investigations of novel perforates to further reduce drag were stymied by measurement uncertainty that prevented further resolution of drag differences between configurations. The current study details efforts to understand the causes of this uncertainty and describe changes made to testing methods to reduce it. Previously evaluated perforates were retested with these improvements in the NASA Langley Grazing Flow Impedance Tube to determine their relative drag at a flow speed of Mach 0.5 without acoustic excitation. Their performance was compared to the previously identified low-drag geometry to see if further drag reductions could be realized.

drag↗

Acoustic Liner Drag: Measurement Uncertainty Reduction and Application to Novel Perforate Geometries

In addition to developing acoustic liner concepts and characterizing their performance, the NASA Langley Liner Physics Team has investigated issues related to liner drag. A method to quantify relative drag of liner configurations was developed and employed to understand the effects of various liner features. It was observed that perforate shape could have a marked effect on the measured liner drag and led to the discovery of a low-drag geometry that cut the drag penalty between a perforate and a smooth wall by approximately 50%. Investigations of novel perforates to further reduce drag were stymied by measurement uncertainty that prevented further resolution of drag differences between configurations. The current study details efforts to understand the causes of this uncertainty and describe changes made to testing methods to reduce it. Previously evaluated perforates were retested with these improvements in the NASA Langley Grazing Flow Impedance Tube to determine their relative drag at a flow speed of Mach 0.5 without acoustic excitation. Their performance was compared to the previously identified low-drag geometry to see if further drag reductions could be realized.

drag↗

Gaussian windows: A tool for exploring multivariate data

Presented here is a method for interactively exploring a large set of quantitative multivariate data, in order to estimate the shape of the underlying density function. It is assumed that the density function is more or less smooth, but no other specific assumptions are made concerning its structure. The local structure of the data in a given region may be examined by viewing the data through a Gaussian window, whose location and shape are chosen by the user. A Gaussian window is defined by giving each data point a weight based on a multivariate Gaussian function. The weighted sample mean and sample covariance matrix are then computed, using the weights attached to the data points. These quantities are used to compute an estimate of the shape of the density function in the window region. The local structure of the data is described by a method similar to the method of principal components. By taking many such local views of the data, we can form an idea of the structure of the data set. The method is applicable in any number of dimensions. The method can be used to find and describe simple structural features such as peaks, valleys, and saddle points in the density function, and also extended structures in higher dimensions. With some practice, we can apply our geometrical intuition to these structural features in any number of dimensions, so that we can think about and describe the structure of the data. Since the computations involved are relatively simple, the method can easily be implemented on a small computer.

Jaeckel, Louis A.↗

Accurate pointing of tungsten welding electrodes

Thoriated-tungsten is pointed accurately and quickly by using sodium nitrite. Point produced is smooth and no effort is necessary to hold the tungsten rod concentric. The chemically produced point can be used several times longer than ground points. This method reduces time and cost of preparing tungsten electrodes.

Ziegelmeier, P.↗

Implicit determination of multispectral scanner data variation over extended areas.

Description of a method by which variations in aerial multispectral scanner data can be accounted for when there is insufficient ground truth information. The data are suitably filtered and smoothed along the flight line, and the reduced results are presented in a form such that the investigator can readily determine the preprocessing that should be used and the spatial extent to which it should be applied. An algorithm is developed for the detection of cloud-shadowed areas.

Kriegler, F. J.↗

A simulator evaluation of the use of spoilers on a light aircraft

A fixed-base flight simulator was used to evaluate wing spoilers for longitudinal flight path control on a modified Cessna Cardinal aircraft. Spoilers which generated the proper pitching moment to maintain aircraft trim C sub L constant could be used as an effective descent rate control. More than 100 simulated ILS approaches were flown by evaluation pilots using both conventional methods and spoiler descent rate control. Three spoiler control schemes were evaluated during the ILS approaches. Using the spoilers for control, instrument approaches could be flown smoothly and precisely with constant airspeed and pitch attitude. While the spoilers could adequately control ILS approaches, a spoiler system with greater authority would be desirable for use in visual approaches.

Brainerd, C. H.↗

A radar image of Venus.

Radar scans of Venus have yielded a brightness map of a large portion of the surface. The bright area in the south (alpha) and the twin such areas in the north (beta and delta) were first discovered by spectral analysis of radar echos. When range-gating is also applied, their shapes are revealed, and they are seen to be roundish and about 1000 km across. Although radar brightness can be the result of either intrinsic reflectivity or surface roughness, polarization studies show these features to be rough (to the scale of the wavelength, 12.5 cm). Dark, circular areas can also be seen, many with bright central spots. The dark areas are probably smooth. The blurring of the equatorial strip is an artifact of the range-Doppler geometry; all resolution disappears at the equator. Another artifact of the method is the 'ghost', in the south, of the images of beta and delta. Such ghosts appear only at the eastern and western extremes of the map.

Goldstein, R. M.↗

Boundary-layer transition on a plate subjected to simultaneous spanwise and chordwise pressure gradients

The boundary-layer transition on a short plate was studied by means of the china-clay visual technique. The plate model was mounted in a wind tunnel so that it was subjected to small simultaneous spanwise and chordwise pressure gradients. Results of the experimental study, which was performed at three subsonic velocities, indicated that the transition pattern was appreciably curved in the spanwise direction but quite smooth and well behaved. Reasonable comparisons between predictions of transition and experiment were obtained from two finite-difference two-dimensional boundary-layer calculation methods which incorporated transition models based on the concept of a transition intermittency factor.

Boldman, D. R.↗

Multilevel techniques for nonelliptic problems

Multigrid and multilevel methods are extended to the solution of nonelliptic problems. A framework for analyzing these methods is established. A simple nonelliptic problem is given, and it is shown how a multilevel technique can be used for its solution. Emphasis is on smoothness properties of eigenvectors and attention is drawn to the possibility of conditioning the eigensystem so that eigenvectors have the desired smoothness properties.

Jespersen, D. C.↗

Thrust imbalance of solid rocket motor pairs on Space Shuttle flights

This analysis extends the investigation presented at the 17th Joint Propulsion Conference in 1981 to include fifteen sets of Space Shuttle flight data. The previous report dealt only with static test data and the first flight pair. The objective is to compare the authors' previous theoretical analysis of thrust imbalance with actual Space Shuttle performance. The theoretical prediction method, which involves a Monte Carlo technique, is reviewed briefly as are salient features of the flight instrumentation system and the statistical analysis. A scheme for smoothing flight data is discussed. The effects of changes in design parameters are discussed with special emphasis on the filament wound motor case being developed to replace the steel case. Good agreement between the predictions and the flight data is demonstrated.

Foster, W. A., Jr.↗

X-ray testing of the AXAF Technology Mirror Assembly (TMA) mirror

The X-ray performance of the Advanced X-ray Astrophysics Facility TMA in the 0.277-6.4-keV range. The precise figure and smooth surface originally specified (which were expected from the in-process optical and mechanical measurements) are verified. The paper primarily describes the test equipment and methods used and gives examples of the detailed on-axis and off-axis data.

Schwartz, D. A.↗