Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “smoothing technique”

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 181 records · Page 10

Uniformly high order accurate essentially non-oscillatory schemes 3

In this paper (a third in a series) the construction and the analysis of essentially non-oscillatory shock capturing methods for the approximation of hyperbolic conservation laws are presented. Also presented is a hierarchy of high order accurate schemes which generalizes Godunov's scheme and its second order accurate MUSCL extension to arbitrary order of accuracy. The design involves an essentially non-oscillatory piecewise polynomial reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell. The reconstruction algorithm is derived from a new interpolation technique that when applied to piecewise smooth data gives high-order accuracy whenever the function is smooth but avoids a Gibbs phenomenon at discontinuities. Unlike standard finite difference methods this procedure uses an adaptive stencil of grid points and consequently the resulting schemes are highly nonlinear.

Harten, A.↗

Enhanced image products by limb-correction and smoothing of VAS radiance measurements

Limb darkening and brightening affect the GOES satellite's VISSR Atmospheric Sounder instrument's imagery with varying satellite viewing angles and instrument noise; limb effects are also associated with a vertical shift of a wavelength's characteristic weighting function. The radiances in an image can be normalized to a reference zenith angle by limb-correcting the image. An additional undesirable effect is that of instrument noise, which varies as a function of channel characteristics and is greatest for channels sensitive to the upper atmosphere. The statistical structure function analysis method of Hillger et al. (1987) is used to obtain noise levels and optimal spatial averaging distances; the averaging domain size from the analysis is used to produce smoothed images, via a running average technique.

Jones, Andrew S.↗

A generalized scheme for three-dimensional hyperbolic grid generation

A hyperbolic grid generation scheme formulated from grid orthogonality and cell volume specification is generalized such that high quality three-dimensional grids can be obtained for a wide variety of geometries. The speed of the scheme is one to two orders of magnitude faster than typical elliptic grid generation methods. The robustness of the scheme is significantly enhanced by several new techniques. By using a spatially variable smoothing coefficient, orthogonality and smoothness of the grid are maintained around complex concave and convex surface topologies. A metric correction procedure is employed to guarantee that the grid marches out of a corner by bisecting the angles subtended by the neighboring points. Extra robustness at severe convex corners is achieved by special local treatments. Different extrapolation methods are used to provide smoothness at floating edges and axis regions. The versatility of the new hyperbolic grid generation scheme is demonstrated by three-dimensional grids generated for external components of the intergrated Space Shuttle vehicle and the SOFIA telescope.

Chan, William M.↗

Reduction of Marine Magnetic Data for Modeling the Main Field of the Earth

The marine data set archived at the National Geophysical Data Center (NGDC) consists of shipborne surveys conducted by various institutes worldwide. This data set spans four decades (1953, 1958, 1960-1987), and contains almost 13 million total intensity observations. These are often less than 1 km apart. These typically measure seafloor spreading anomalies with amplitudes of several hundred nanotesla (nT) which, since they originate in the crust, interfere with main field modeling. The source for these short wavelength features are confined within the magnetic crust (i.e., sources above the Curie isotherm). The main field, on the other hand, is of much longer wavelengths and originates within the earth's core. It is desirable to extract the long wavelength information from the marine data set for use in modeling the main field. This can be accomplished by averaging the data along the track. In addition, those data which are measured during periods of magnetic disturbance can be identified and eliminated. Thus, it should be possible to create a data set which has worldwide data distribution, spans several decades, is not contaminated with short wavelengths of the crustal field or with magnetic storm noise, and which is limited enough in size to be manageable for the main field modeling. The along track filtering described above has proved to be an effective means of condensing large numbers of shipborne magnetic data into a manageable and meaningful data set for main field modeling. Its simplicity and ability to adequately handle varying spatial and sampling constraints has outweighed consideration of more sophisticated approaches. This filtering technique also provides the benefits of smoothing out short wavelength crustal anomalies, discarding data recorded during magnetically noisy periods, and assigning reasonable error estimates to be used in the least square modeling. A useful data set now exists which spans 1953-1987.

Baldwin, R. T.↗

Large area soft x-ray collimator to facilitate x-ray optics testing

The first objective of this program is to design a nested conical foil x-ray optic which will collimate x-rays diverging from a point source. The collimator could then be employed in a small, inexpensive x-ray test stand which would be used to test various x-ray optics and detector systems. The second objective is to demonstrate the fabrication of the x-ray reflectors for this optic using lacquer-smoothing and zero-stress electroforming techniques.

Espy, Samuel L.↗

Detailed measurements in the SSME high pressure fuel turbine with smooth rotor blades

Several tests of the Rocketdyne configuration of the Space Shuttle Main Engine (SSME) High Pressure Fuel Turbopump (HPFTP) Turbine have been completed in the Turbine Test Equipment (TTE) at Marshall Space Flight Center. The tests involved using scaled performance parameters and model measurements to predict the performance of the turbine. The overall performance has been the primary objective of the tests to date, but more detailed measurements are also of interest. During the most recent test of the Rocketdyne configuration of the HPFTP turbine with smooth rotor blades, several different measurement techniques were used to study the turbine inlet and exit velocity profiles, boundary layer thicknesses, turbulence intensities, etc. Data has been obtained using various hot film probes and three-hole cobra probes. Laser Velocimeter measurements were also made. The test plan and test data will be presented and discussed as well as lessons learned on how to obtain the various types of data.

Hudson, Susan T.↗

Multistage Schemes with Multigrid for Euler and Navier-Strokes Equations: Components and Analysis

A class of explicit multistage time-stepping schemes with centered spatial differencing and multigrids are considered for the compressible Euler and Navier-Stokes equations. These schemes are the basis for a family of computer programs (flow codes with multigrid (FLOMG) series) currently used to solve a wide range of fluid dynamics problems, including internal and external flows. In this paper, the components of these multistage time-stepping schemes are defined, discussed, and in many cases analyzed to provide additional insight into their behavior. Special emphasis is given to numerical dissipation, stability of Runge-Kutta schemes, and the convergence acceleration techniques of multigrid and implicit residual smoothing. Both the Baldwin and Lomax algebraic equilibrium model and the Johnson and King one-half equation nonequilibrium model are used to establish turbulence closure. Implementation of these models is described.

Swanson, R. C.↗

Experimental Aerodynamic Simulation of Scallop Ice Accretion on a Swept Wing

Understanding the aerodynamic impact of swept-wing ice accretions is a crucial component of the design of modern aircraft. Computer-simulation tools are commonly used to approximate ice shapes, so the necessary level of detail or fidelity of those simulated ice shapes must be understood relative to high-fidelity representations of the ice. Previous tests were performed in the NASA Icing Research Tunnel to acquire high-fidelity ice shapes. From this database, full-span artificial ice shapes were designed and manufactured for both an 8.9%-scale and 13.3%-scale semispan wing model of the CRM65 which has been established as the full-scale baseline for this swept-wing project. These models were tested in the Walter H. Beech wind tunnel at Wichita State University and at the ONERA F1 facility, respectively. The data collected in the Wichita St. University wind tunnel provided a low-Reynolds number baseline study while the pressurized F1 facility produced data over a wide range of Reynolds and Mach numbers with the highest Reynolds number studied being approximately Re=11.9x106. Past work focused on only three different fidelity variations for ice shapes based on multiple icing conditions. This work presents a more detailed investigation into several fidelity representations of a single highly three-dimensional scallop ice accretion. Sensitivity to roughness size and application technique on a low-fidelity smooth ice shape is described. The data indicate that the aerodynamic performance is not especially sensitive to the grit variations. An ice accretion code was also used to generate ice shapes for aerodynamic testing and comparisons. These ice shapes have a general appearance like the low-fidelity smooth ice shapes, but in this case, the computer-generated ice shape is significantly smaller. As such, the impact of that ice shape on the aerodynamic performance of the wing is reduced compared to the smooth ice shape based on the icing experiment for those same conditions. Spanwise discontinuities were also introduced to a low-fidelity ice shape in an attempt to quantify the impact of those variation in the high-fidelity ice shape. While the lift data indicate good agreement between the high-fidelity ice shapes and the low-fidelity ice shapes and the low-fidelity ice shapes with spanwise discontinuities, a closer investigation of the data suggests potential, significant differences in the flowfield. These results were similar at both facilities over the wide range of test conditions utilized.

Woodard, Brian↗

Digital enhancement of computerized axial tomograms

A systematic evaluation has been conducted of certain digital image enhancement techniques performed in image space. Three types of images have been used, computer generated phantoms, tomograms of a synthetic phantom, and axial tomograms of human anatomy containing images of lesions, artificially introduced into the tomograms. Several types of smoothing, sharpening, and histogram modification have been explored. It has been concluded that the most useful enhancement techniques are a selective smoothing of singular picture elements, combined with contrast manipulation. The most useful tool in applying these techniques is the gray-scale histogram.

Roberts, E., Jr.↗

An interactive grid generation technique for turbomachinery

A combination algebraic/elliptic technique is presented for the generation of 3-D grids about turbomachinery blade rows for both axial and radial flow machinery. The technique is build around use of an advanced engineering workstation to construct several 2-D grids interactively on predetermined blade-to-blade surfaces. A 3-D grid is generated by interpolating these surface grids onto an axisymmetric grid. On each blade to blade surface, a grid is created using algebraic techniques near the blade to control orthogonality within the boundary layer region and elliptic techniques in the mid-passage to achieve smoothness. The interactive definition of bezier curves as internal boundaries is the key to simple construction. The approach is adapted for use with the average passage solution technique, although this is not a limitation for most other uses. A variety of examples are presented.

Beach, Tim↗

TOWARD A METHOD FOR SCALING HUMAN BODY MODELS IN AN IMU-BASED WORKFLOW

BACKGROUND Scaled biomechanical models can more accurately inform crew health decisions when tailored to the wide range of astronaut sizes. One component to improve scaling of existing models to better represent each unique astronaut’s size is the individual length scaling of limbs. Traditionally, these lengths are determined by motion capture or manual measurement. A new method is herein proposed for length scaling which can be done by measuring linear and angular accelerations at a desired point during isolated motion around a point of rotation, then calculating the distance between the desired point and point of rotation. When an Inertial Measurement Unit (IMU) device is placed at the distal point of a limb, the isolated motion is about that limb’s proximal joint. For example, to measure forearm length, an IMU is placed at the wrist, the point of rotation is at the elbow, and the isolated motion is forearm flexion and extension. These calculated lengths are then used to scale models to each unique astronaut’s size, thereby improving the applicability of the model. This method of scaling limb segments can be used for any limb that has an easily defined proximal joint for the limb to rotate around including hands, arms, legs, feet. Utilizing IMUs for data collection also provides the synergistic ability to record data without a dedicated space in a room with many cameras, therefore reducing the data collection footprint, or record data where optical motion capture is not possible, such as inside a spacesuit. METHODS AND RESULTS To test this method, upper body data collection was performed with 5 Xsens DOT IMUs on a single subject. IMUs consist of an accelerometer, a gyroscope, and a magnetometer which collect linear acceleration, angular velocity, and magnetic fluctuations, respectively. Before any ground-based laboratory collection, the magnetic fluctuations are used to correct the heading of the IMU in space relative to the Earth’s magnetic field. The direct measurement of angular velocity is integrated to calculate angular acceleration. Then the linear acceleration ( a ) and angular acceleration (α) are solved using r = at/α to calculate the radius, which in this case is the distance between the IMU and the point of rotation (i.e., segment length). The IMU must be placed at the most distal point of the limb being measured (i.e., ankle if measuring lower leg length) and the test plan must consist of an isolated motion about that limb’s proximal joint (i.e., knee flexion and extension if measuring lower leg length). The distances (radii) calculated at every time interval were filtered (bandpass filter keeping 5th-90th percentile data) to eliminate outliers and spurious data that occur when the isolated motion was stopped or nearly stopped. The remaining distances were averaged, resulting in the calculated limb length. Scaling factors were then computed by dividing the calculated limb length by the unscaled model’s length. These scale factors are plugged into the Scale Tool in OpenSim [1,2] to apply the scaling to the OpenSim Full Body Rajagopal Model [3,4]. Manual measurements of limb lengths were taken before data collection started and used for comparing against the calculated lengths. The Anthropometric Survey of US Army Personnel (ANSUR II) [5] was also used as a third source of reference for limb length measurements. The following measurements were retrieved from the subject before data collection: 34.5 cm from L1 to C7 (thorax), 25.7 cm from C7 Joint Center (JC) to head vertex (neck and head), 36.3 cm from shoulder JC to elbow JC (humerus), 29.5 cm from elbow JC to wrist JC (forearm), and 16.2 cm from clavicle to acromion (clavicle). Of those five, forearm and humerus lengths were calculated using this proposed method to obtain preliminary results. The forearm length after filtering and averaging was calculated to be 37.6 cm. This is a 28% overestimation from the measured forearm length (29.5 cm). The humerus length after filtering and averaging was calculated to be 47.8 cm. This is a 31% difference from the measured subject length (36.3 cm). Sources of error include imperfect isolated motion (method currently expects that motion should be perfectly circular in a 2D plane, include no rotation of the IMU, and be relatively smooth; a more secure IMU attachment method will help), unrefined filter techniques (removed highest and lowest values with 20% high and low pass filters and no smoothing filters), arbitrary removal of stopped or near stopped data (kept data for only a short range before and after the angular velocity peaking), and a more representative method for removing gravitational acceleration is needed (current method is to zero all accelerations against a baseline taken just before the isolated motion which does not account for the gravitational acceleration changed due to IMU rotation during movement). Addressing these error sources will improve the accuracy of the limb length calculation. Next steps include creating a method for whole-body scaling estimation using individual limb scale factors. Continued pursuit of these techniques is expected to enable acquiring anthropometric information using only IMUs in real-time.

E. K. Marecki↗

Numerical comparison of convective heat transfer augmentation devices used in cooling channels of hypersonic vehicles

Hypersonic vehicles are exposed to extreme thermal conditions compared to subsonic aircraft; therefore, some level of thermal management is required to protect the materials used. Normally, hypersonic vehicles experience the highest temperatures in the nozzle throat, and aircraft and propulsion system leading edges. Convective heat transfer augmentation techniques can be used in the thermal management system to increase heat transfer of the cooling channels in those areas. The techniques studied in this report are pin-fin, offset-fin, ribbed and straight roughened channel. A smooth straight channel is used as the baseline for comparing the techniques. SINDA '85, a lumped parameter finite difference thermal analyzer, is used to model the channels. Subroutines are added to model the fluid flow assuming steady one dimensional compressible flow with heat addition and friction. Correlations for convective heat transfer and friction are used in conjunction with the fluid flow analysis mentioned. As expected, the pin-fin arrangement has the highest heat transfer coefficient and the largest pressure drop. All the other devices fall in between the pin-fin and smooth straight channel. The selection of the best heat augmentation method depends on the design requirements. A good approach may be a channel using a combination of the techniques. For instance, several rows of pin-fins may be located at the region of highest heat flux, surrounded by some of the other techniques. Thus, the heat transfer coefficient is maximized at the region of highest heat flux while the pressure drop is not excessive.

Maldonado, Jaime J.↗

Influence of gaseous hydrogen on the mechanical properties of high temperature alloys

Tensile tests of six nickel-base and one cobalt-base alloy were conducted in 34.5 MN/sq m helium and hydrogen environments at temperatures from 297 K to 1,088 K. Mechanical properties tests of the nickel-base alloy MAR M-246 (Hf modified), in two cast conditions, were conducted in gaseous environments at temperatures from 297 K to 1,144 K and pressures from one atmosphere to 34.5 MN/sq m. The objective of this program was to obtain the mechanical properties of the various alloys proposed for use in space propulsion systems in a pure hydrogen environment at different temperatures and to compare with the mechanical properties in helium at the same conditions. All testing was conducted on solid specimens exposed to external gaseous pressure. Smooth and notched tensile properties were determined using ASTM tensile testing techniques, and creep-rupture life was determined using ASTM creep-rupture techniques. Low-cycle fatigue life was established by constant total strain and constant stress testing using smooth specimens and a closed-loop test machine.

Source record↗

Wavelet Sparse Approximate Inverse Preconditioners

There is an increasing interest in using sparse approximate inverses as preconditioners for Krylov subspace iterative methods. Recent studies of Grote and Huckle and Chow and Saad also show that sparse approximate inverse preconditioner can be effective for a variety of matrices, e.g. Harwell-Boeing collections. Nonetheless a drawback is that it requires rapid decay of the inverse entries so that sparse approximate inverse is possible. However, for the class of matrices that, come from elliptic PDE problems, this assumption may not necessarily hold. Our main idea is to look for a basis, other than the standard one, such that a sparse representation of the inverse is feasible. A crucial observation is that the kind of matrices we are interested in typically have a piecewise smooth inverse. We exploit this fact, by applying wavelet techniques to construct a better sparse approximate inverse in the wavelet basis. We shall justify theoretically and numerically that our approach is effective for matrices with smooth inverse. We emphasize that in this paper we have only presented the idea of wavelet approximate inverses and demonstrated its potential but have not yet developed a highly refined and efficient algorithm.

Chan, Tony F.↗

Method of measuring elastic strain distribution in specimens used for high frequency fatigue testing

A technique using a point contact piezoelectric detector has been developed which permits the displacement amplitude distributions over the entire length of high frequency fatigue specimens to be measured. Elastic strain distributions can be derived with a computer analysis that uses curve smoothing and numerical derivative methods. The equipment and techniques for measuring and processing the data are described, and displacement amplitudes, strain distributions, and magnification factors for straight rods and dumbbell-shaped specimens are obtained and shown to agree with theoretically predicted values. It is shown that the technique presented here can be used to design specimens with unusual displacement amplitudes and strain distributions.

Sirian, C. R.↗

The expanded LaGrangian system for constrained optimization problems

Smooth penalty functions can be combined with numerical continuation/bifurcation techniques to produce a class of robust and fast algorithms for constrainted optimization problems. The key to the development of these algorithms is the Expanded Lagrangian System which is derived and analyzed in this work. This parameterized system of nonlinear equations contains the penalty path as a solution, provides a smooth homotopy into the first-order necessary conditions, and yields a global optimization technique. Furthermore, the inevitable ill-conditioning present in a sequential optimization algorithm is removed for three penalty methods: the quadratic penalty function for equality constraints, and the logarithmic barrier function (an interior method) and the quadratic loss function (an interior method) for inequality constraints. Although these techniques apply to optimization in general and to linear and nonlinear programming, calculus of variations, optimal control and parameter identification in particular, the development is primarily within the context of nonlinear programming.

Poore, A. B.↗

The expanded Lagrangian system for constrained optimization problems

Smooth penalty functions can be combined with numerical continuation/bifurcation techniques to produce a class of robust and fast algorithms for constrained optimization problems. The key to the development of these algorithms is the Expanded Lagrangian System which is derived and analyzed in this work. This parameterized system of nonlinear equations contains the penalty path as a solution, provides a smooth homotopy into the first-order necessary conditions, and yields a global optimization technique. Furthermore, the inevitable ill-conditioning present in a sequential optimization algorithm is removed for three penalty methods: the quadratic penalty function for equality constraints, and the logarithmic barrier function (an interior method) and the quadratic loss function (an interior method) for inequality constraints. Although these techniques apply to optimization in general and to linear and nonlinear programming, calculus of variations, optimal control and parameter identification in particular, the development is primarily within the context of nonlinear programming.

Poore, A. B.↗

Scale effects: HCMM data simulation. Usage of filtering techniques for scaling-up simulations

Image reduction used to simulate increase in altitude of an acquisition platform is equivalent to data smoothing, and can be achieved either by neighborhood averaging or by filtering techniques. The averaging approach is limited for accurate simulation. A filtering method is described which was based on the hypothesis that all changes due to altitude increase can be represented by a point spread function. Determination of the scale function and factor are discussed as well as implementation of the filtering. Filtering can be either in the spatial or frequency domain. In the spatial domain, filtering consists of the convolution of the image with the weights mask, and then of the declination of the points according to the appropriates scale factor. A simulation of an aircraft day image in the infrared channel is examined.

Digennaro, V.↗