Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical techniques”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,135 records · Page 63

Student Support for Research in Hierarchical Control and Trajectory Planning

Generally, classical polynomial splines tend to exhibit unwanted undulations. In this work, we discuss a technique, based on control principles, for eliminating these undulations and increasing the smoothness properties of the spline interpolants. We give a generalization of the classical polynomial splines and show that this generalization is, in fact, a family of splines that covers the broad spectrum of polynomial, trigonometric and exponential splines. A particular element in this family is determined by the appropriate control data. It is shown that this technique is easy to implement. Several numerical and curve-fitting examples are given to illustrate the advantages of this technique over the classical approach. Finally, we discuss the convergence properties of the interpolant.

Martin, Clyde F.↗

Numerical solution of equations governing longitudinal suspension line wave motion during the parachute unfurling process

Equations are presented which govern the dynamics of the lines-first parachute unfurling process, including wave motion in the parachute suspension lines. Techniques are developed for obtaining numerical solutions to the governing equations. Histories of tension at test data, and generally good agreement is observed. Errors in computed results are attributed to several areas of uncertainty, the most significant being a poorly defined boundary condition on the wave motion at the vehicle-suspension line boundary.

Poole, L. R.↗

Development of Optical Diagnostic Techniques for Microgravity Materials Processing

Materials processing including crystal growth, either under a gravity environment on ground or a microgravity environment in space, involves complicated phenomena of fluid motions in gas or liquid phases as well as interaction of various species. To obtain important physical insight, it is very necessary to provide gross-field optical diagnostics for monitoring various physical properties. Materials processing inhibits easy access by ordinary instruments and thus characterizing gross-field physical properties is very challenging. Typical properties of importance can be fluid velocity, temperature, and species concentration for fluids, and surface topology and defects for solids. Observing surface grow rate during crystal growth is also important. Material microstructures, i.e., integrity of crystal structures, is strongly influenced by the existence of thermally-induced flow as well as local nucleation of particles during solidification, which may act in many detrimental ways. In both ground-based and microgravity experiments, the nature of product property changes resulting from three-dimensional fluid or particle motions need be characterized. Gross-field diagnostics is thus required to identify their effects on product defects and process deficiencies. The quantitative visualization techniques can also be used for validation of numerical modeling. For optical nonintrusive gross-field diagnostic techniques, two approaches were developed as summer projects. One optical approach allows us to provide information of species concentration and temperature for monitoring in real time. The other approach, that is, the concept which is formulated for detection of surface topography measurement can provide unprecedented spatial resolution during crystal growth.

Cha, Soyoung Stephen↗

Microstress Analysis of Periodic Composites

Local elastic fields in the unit cell of a periodic composite are examined numerically with an integral equation approach. Techniques of Fourier series and Green's functions are used to construct the integral equations. Numerical solutions are obtained using the Fourier series approach with rectangular subvolume elements. Specific results are given for a tungsten/copper metal matrix composite.

Walker, Kevin P.↗

A Tower-based Prototype VHF/UHF Radar for Subsurface Sensing: System Description and Data Inversion Results

Knowledge of subsurface characteristics such as permittivity variations and layering structure could provide a breakthrough in many terrestrial and planetary science disciplines. For Earth science, knowledge of subsurface and subcanopy soil moisture layers can enable the estimation of vertical flow in the soil column linking surface hydrologic processes with that in the subsurface. For planetary science, determining the existence of subsurface water and ice is regarded as one of the most critical information needs for the study of the origins of the solar system. The subsurface in general can be described as several near-parallel layers with rough interfaces. Each homogenous rough layer can be defined by its average thickness, permittivity, and rms interface roughness assuming a known surface spectral distribution. As the number and depth of layers increase, the number of measurements needed to invert for the layer unknowns also increases, and deeper penetration capability would be required. To nondestructively calculate the characteristics of the rough layers, a multifrequency polarimetric radar backscattering approach can be used. One such system is that we have developed for data prototyping of the Microwave Observatory of Subcanopy and Subsurface (MOSS) mission concept. A tower-mounted radar makes backscattering measurements at VHF, UHF, and L-band frequencies. The radar is a pulsed CW system, which uses the same wideband antenna to transmit and receive the signals at all three frequencies. To focus the beam at various incidence angles within the beamwidth of the antenna, the tower is moved vertically and measurements made at each position. The signals are coherently summed to achieve focusing and image formation in the subsurface. This requires an estimate of wave velocity profiles. To solve the inverse scattering problem for subsurface velocity profile simultaneously with radar focusing, we use an iterative technique based on a forward numerical solution of the layered rough surface problem. The layers are each defined in terms of a small number of unknown distributions as given above. An a priori estimate of the solution is first assumed, based on which the forward problem is solved for the backscattered measurements. This is compared with the measured data and using iterative techniques an update to the solution for the unknowns is calculated. The process continues until convergence is achieved. Numerical results will be shown using actual radar data acquired with the MOSS tower radar system in Arizona in Fall 2003, and compared with in-situ measurements.

Moghaddam, Mahta↗

Numerical solutions of transonic flows by parametric differentiation and integral equation techniques

The paper discusses results of an exploratory study of the advantages obtained by combining integral equation methods with the method of parametric differentiation in the treatment of transonic flow problems. In the proposed method, the nonlinear unsteady transonic flow equation for small perturbations is transformed into a linear equation by the use of the method of parametric differentiation. The linear equation is split into a pair of weakly coupled partial differential equations by writing the transformed perturbation potential as the sum of a steady component and an unsteady component. The solution of the steady equation as an integral equation is based on Ogana's treatment (1978). As a test case, the formulation is applied to predict the steady transonic flow over a nonlifting parabolic-arc airfoil.

Sivaneri, N. T.↗

The 3D widgets for exploratory scientific visualization

Computational fluid dynamics (CFD) techniques are used to simulate flows of fluids like air or water around such objects as airplanes and automobiles. These techniques usually generate very large amounts of numerical data which are difficult to understand without using graphical scientific visualization techniques. There are a number of commercial scientific visualization applications available today which allow scientists to control visualization tools via textual and/or 2D user interfaces. However, these user interfaces are often difficult to use. We believe that 3D direct-manipulation techniques for interactively controlling visualization tools will provide opportunities for powerful and useful interfaces with which scientists can more effectively explore their datasets. A few systems have been developed which use these techniques. In this paper, we will present a variety of 3D interaction techniques for manipulating parameters of visualization tools used to explore CFD datasets, and discuss in detail various techniques for positioning tools in a 3D scene.

Herndon, Kenneth P.↗

Modeling of Convective-Stratiform Precipitation Processes: Sensitivity to Partitioning Methods and Numerical Advection Schemes

Six different convective-stratiform separation techniques, including a new technique that utilizes the ratio of vertical and terminal velocities, are compared and evaluated using two-dimensional numerical simulations of a tropical [Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE)] and midlatitude continental [Preliminary Regional Experiment for STORM-Central (PRESTORM)] squall line. The simulations are made using two different numerical advection schemes: 4th order and positive definite advection. Comparisons are made in terms of rainfall, cloud coverage, mass fluxes, apparent heating and moistening, mean hydrometeor profiles, CFADs (Contoured Frequency with Altitude Diagrams), microphysics, and latent heating retrieval. Overall, it was found that the different separation techniques produced results that qualitatively agreed. However, the quantitative differences were significant. Observational comparisons were unable to conclusively evaluate the performance of the techniques. Latent heating retrieval was shown to be sensitive to the use of separation technique mainly due to the stratiform region for methods that found very little stratiform rain. The midlatitude PRESTORM simulation was found to be nearly invariant with respect to advection type for most quantities while for TOGA COARE fourth order advection produced numerous shallow convective cores and positive definite advection fewer cells that were both broader and deeper penetrating above the freezing level.

Lang, Steve↗

Two alternative ways for solving the coordination problem in multilevel optimization

Two techniques for formulating the coupling between levels in multilevel optimization by linear decomposition, proposed as improvements over the original formulation, now several years old, that relied on explicit equality constraints which were shown by application experience as occasionally causing numerical difficulties. The two new techniques represent the coupling without using explicit equality constraints, thus avoiding the above diffuculties and also reducing computational cost of the procedure. The old and new formulations are presented in detail and illustrated by an example of a structural optimization. A generic version of the improved algorithm is also developed for applications to multidisciplinary systems not limited to structures.

Sobieszczanski-Sobieski, Jaroslaw↗

Cloud Filtering Using a Bi-Spectral Spatial Coherence Approach

The research in this paper focuses on describing a technique developed for cloud filtering using a bi-spectral approach on GOES-8/9 Imager data. The application was developed for use with infrared retrievals of geophysical parameters in mind, where cloud cover contaminates the derived product. However, numerous potential applications of the technique exist. The technique will be described and a preliminary validation of the algorithm will be presented. Although initially based on the spatial coherence approach from Coakley and Brethereton (1982), it has evolved to utilize a difference image of the I I and 3.9 micrometer channels on the GOES-8/9 Imager. This image is very similar to that produced by Nelson and Ellrod (1996). During the daytime the technique makes use of the varying solar reflectance in the 3.9 micrometer channel by clouds and land to identify cloudy pixels. While at night, the technique makes use of the varying emissivity of the clouds in the scene to discriminate between clear and cloudy pixels. The technique applies three basic threshold tests to produce the final cloud filtered image: 1) a standard deviation threshold to detect the spatial variance in the scene, 2) a difference threshold between adjacent pixels, and 3) a simple infrared temperature threshold. The first test is applied to the entire image at once, then in a second pass the next two tests are applied. The final infrared temperature threshold is only meant to identify the coldest clouds that might pass the previous tests. The technique performs well during the daytime, while nighttime performance is degraded but is promising. The technique has proven to be robust and shows great promise of meeting its original goal of cloud filtering for use in an infrared retrieval algorithm for use in climate studies.

Guillory, Anthony R.↗

Cloud Filtering Using a Bi-Spectral Spatial Coherence Approach

The research in this paper focuses on describing a technique developed for cloud filtering using a bi-spectral approach on GOES-8/9 Imager data. The application was developed for use with infrared retrievals of geophysical parameters in mind, where cloud cover contaminates the derived product. However, numerous potential applications of the technique exist. The technique will be described and a preliminary validation of the algorithm will be presented. Although initially based on the spatial coherence approach from Coakley and Brethereton (1982), it has evolved to utilize a difference image of the 11 and 3.9 micrometer channels on the GOES-8/9 Imager. This image is very similar to that produced by Nelson and Ellrod (1996). During the daytime the technique makes use of the varying solar reflectance in the 3.9 micrometer channel by clouds and land to identify cloudy pixels. While at night, the technique makes use of the varying emissivity of the clouds in the scene to discriminate between clear and cloudy pixels. The technique applies three basic threshold tests to produce the final cloud filtered image: 1) a standard deviation threshold to detect the spatial variance in the scene, 2) a difference threshold between adjacent pixels, and 3) a simple infrared temperature threshold. The first test is applied to the entire image at once, then in a second pass the next two tests are applied. The final infrared temperature threshold is only meant to identify the coldest clouds that might pass the previous tests. The technique performs well during the daytime, while nighttime performance is degraded but is promising. The technique has proven to be robust and shows great promise of meeting its original goal of cloud filtering for use in an infrared retrieval algorithm for use in climate studies.

Guillory, Anthony R.↗

An analysis of the accuracy of a parameter optimization

The numerical operations involved in a currently used optimization technique are discussed and analyzed with special attention to the numerical accuracy. Alternative methods for deriving linear system transfer functions, finding the relationships between the transfer function coefficients and the design parameters, and solving a matrix equation are presented for more accurate and cost effective solutions.

Baram, Y.↗

Computation of radiation view factors for surfaces with obstructed views of each other

Radiation calculations begin with a determination of view (or shape) factors describing the geometrical configuration of the exchanging surfaces. These viewfactor calculations can be very involved, time consuming and expensive. This paper describes the development of an algorithm for situations in which the view between two surfaces is obscured by one or more intervening surfaces. A technique based upon the combination of the Contour Integration and Double Area Integration methods and the acceleration techniques necessary for its efficient use are described. Numerical examples are given to demonstrate the use of the technique and to describe some of the problems inherent in treating radiation in enclosures.

Emery, A. F.↗

Effective Use of Cray Supercomputers

The effective use of a supercomputer depends on many aspects, including the ability of users to write efficient programs to use the resources of the system in an optimal manner. However, it is the responsibility of the system managers of these systems to insure the maximum effectiveness of the overall system is achieved. Many varying techniques have been developed at the Numerical Aerodynamic Simulation (NAS) Program to advance the management of these critical systems. Many of the issues and techniques used for managing supercomputers are common to multi-user UNIX systems, regardless of the version of UNIX or the power of the hardware. However, a UNICOS supercomputer presents some special challenges and requires additional features and tools to be developed to effectively manage the system. Only part of the challenge is related to performance monitoring and improvement. Much of the responsibility of the system manager is to provide fair and consistent access to the system resources, which is at times a difficult problem. After an introduction to the environment at the Numerical Aerodynamic Simulation Project is given as background this paper will first discuss the areas which are common to UNDC system management. Then the paper will discuss the specific areas of UNICOS which must be used to operate the system efficiently. The paper goes on to discuss the methods of supporting individual users in order to increase their effectiveness and the efficiency of their programs. This is accomplished through a professional support staff who interact on a daily basis to support the NAS scientific client community.

Kramer, William T. C.↗

Surrogate-equation technique for simulation of steady inviscid flow

A numerical procedure for the iterative solution of inviscid flow problems is described, and its utility for the calculation of steady subsonic and transonic flow fields is demonstrated. Application of the surrogate equation technique defined herein allows the formulation of stable, fully conservative, type dependent finite difference equations for use in obtaining numerical solutions to systems of first order partial differential equations, such as the steady state Euler equations. Steady, two dimensional solutions to the Euler equations for both subsonic, rotational flow and supersonic flow and to the small disturbance equations for transonic flow are presented.

Johnson, G. M.↗

Prioritization methodology for chemical replacement

Since United States of America federal legislation has required ozone depleting chemicals (class 1 & 2) to be banned from production, The National Aeronautics and Space Administration (NASA) and industry have been required to find other chemicals and methods to replace these target chemicals. This project was initiated as a development of a prioritization methodology suitable for assessing and ranking existing processes for replacement 'urgency.' The methodology was produced in the form of a workbook (NASA Technical Paper 3421). The final workbook contains two tools, one for evaluation and one for prioritization. The two tools are interconnected in that they were developed from one central theme - chemical replacement due to imposed laws and regulations. This workbook provides matrices, detailed explanations of how to use them, and a detailed methodology for prioritization of replacement technology. The main objective is to provide a GUIDELINE to help direct the research for replacement technology. The approach for prioritization called for a system which would result in a numerical rating for the chemicals and processes being assessed. A Quality Function Deployment (QFD) technique was used in order to determine numerical values which would correspond to the concerns raised and their respective importance to the process. This workbook defines the approach and the application of the QFD matrix. This technique: (1) provides a standard database for technology that can be easily reviewed, and (2) provides a standard format for information when requesting resources for further research for chemical replacement technology. Originally, this workbook was to be used for Class 1 and Class 2 chemicals, but it was specifically designed to be flexible enough to be used for any chemical used in a process (if the chemical and/or process needs to be replaced). The methodology consists of comparison matrices (and the smaller comparison components) which allow replacement technology to be quantitatively compared in several categories, and a QFD matrix which allows process/chemical pairs to be rated against one another for importance (using consistent categories). Depending on the need for application, one can choose the part(s) needed or have the methodology completed in its entirety. For example, if a program needs to show the risk of changing a process/chemical one may choose to use part of Matrix A and Matrix C. If a chemical is being used, and the process must be changed; one might use the Process Concerns part of Matrix D for the existing process and all possible replacement processes. If an overall analysis of a program is needed, one may request the QFD to be completed.

Cruit, Wendy↗

Finite-element numerical modeling of atmospheric turbulent boundary layer

A dynamic turbulent boundary-layer model in the neutral atmosphere is constructed, using a dynamic turbulent equation of the eddy viscosity coefficient for momentum derived from the relationship among the turbulent dissipation rate, the turbulent kinetic energy and the eddy viscosity coefficient, with aid of the turbulent second-order closure scheme. A finite-element technique was used for the numerical integration. In preliminary results, the behavior of the neutral planetary boundary layer agrees well with the available data and with the existing elaborate turbulent models, using a finite-difference scheme. The proposed dynamic formulation of the eddy viscosity coefficient for momentum is particularly attractive and can provide a viable alternative approach to study atmospheric turbulence, diffusion and air pollution.

Lee, H. N.↗

Nonlinear transient analysis of joint dominated structures

A residual force technique is presented that can perform the transient analyses of large, flexible, and joint dominated structures. The technique permits substantial size reduction in the number of degrees of freedom describing the nonlinear structural models and can account for such nonlinear joint phenomena as free-play and hysteresis. In general, joints can have arbitrary force-state map representations but these are used in the form of residual force maps. One essential feature of the technique is to replace the arbitrary force-state maps describing the nonlinear joints with residual force maps describing the truss links. The main advantage of this replacement is that the incrementally small relative displacements and velocities across a joint are not monitored directly thereby avoiding numerical difficulties. Instead, very small and 'soft' residual forces are defined giving a numerically attractive form for the equations of motion and thereby permitting numerically stable integration algorithms. The technique was successfully applied to the transient analyses of a large 58 bay, 60 meter truss having nonlinear joints. A method to perform link testing is also presented.

Chapman, J. M.↗