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,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

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↗
NASA NTRS1998

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.↗
NASA NTRS1998

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.↗
DOE OSTI2024

Lattice Green’s Functions for High-Order Finite Difference Stencils

Lattice Green's Functions (LGFs) are fundamental solutions to discretized linear operators, and as such they are a useful tool for solving discretized elliptic PDEs on domains that are unbounded in one or more directions. The majority of existing numerical solvers that make use of LGFs rely on a second-order discretization and operate on domains with free-space boundary conditions in all directions. Under these conditions, fast expansion methods are available that enable precomputation of 2D or 3D LGFs in linear time, avoiding the need for brute-force multi-dimensional quadrature of numerically unstable integrals. Here we focus on higher-order discretizations of the Laplace operator on domains with more general boundary conditions, by (1) providing an algorithm for fast and accurate evaluation of the LGFs associated with high-order dimension-split centered finite differences on unbounded domains, and (2) deriving closed-form expressions for the LGFs associated with both dimension-split and Mehrstellen discretizations on domains with one unbounded dimension. Through numerical experiments we demonstrate that these techniques provide LGF evaluations with near machine-precision accuracy, and that the resulting LGFs allow for numerically consistent solutions to high-order discretizations of the Poisson's equation on fully or partially unbounded 3D domains.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1974

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.↗
NASA NTRS1981

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.↗
NASA NTRS1989

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.↗
NASA NTRS1981

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.↗
NASA NTRS1995

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↗
NASA NTRS1979

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.↗
NASA NTRS1987

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.↗
NASA NTRS2017

Automatic Estimation of Verified Floating-Point Round-Off Errors via Static Analysis

This paper introduces a static analysis technique for computing formally verified round-off error bounds of floating-point functional expressions. The technique is based on a denotational semantics that computes a symbolic estimation of floating-point round-o errors along with a proof certificate that ensures its correctness. The symbolic estimation can be evaluated on concrete inputs using rigorous enclosure methods to produce formally verified numerical error bounds. The proposed technique is implemented in the prototype research tool PRECiSA (Program Round-o Error Certifier via Static Analysis) and used in the verification of floating-point programs of interest to NASA.

Moscato, Mariano↗
NASA NTRS1985

Analysis of transitional separation bubbles on infinite swept wings

A previously developed two-dimensional local inviscid-viscous interaction technique for the analysis of airfoil transitional separation bubbles, ALESEP (Airfoil Leading Edge Separation), has been extended for the calculation of transitional separation bubbles over infinite swept wings. As part of this effort, Roberts' empirical correlation, which is interpreted as a separated flow empirical extension of Mack's stability theory for attached flows, has been incorporated into the ALESEP procedure for the prediction of the transition location within the separation bubble. A series of two-dimensional calculations are presented as a verification of the prediction capability of the interaction technique with this transition model. Numerical tests have shown that this two-dimensional natural transition correlation may also be applied to transitional separation bubbles over infinite swept wings. Results of the interaction procedure are compared with Horton's detailed experimental data for separated flow over a swept plate which demonstrates the accuracy of the present technique. The principal conclusion of this paper is that the prediction of transitional separation bubbles over two-dimensional or infinite swept geometries is now possible using the present interacting boundary layer approach.

Davis, R. L.↗
NASA NTRS1990

Theory of biaxial graded-index optical fiber

A biaxial graded-index fiber with a homogeneous cladding is studied. Two methods, wave equation and matrix differential equation, of formulating the problem and their respective solutions are discussed. For the wave equation formulation of the problem it is shown that for the case of a diagonal permittivity tensor the longitudinal electric and magnetic fields satisfy a pair of coupled second-order differential equations. Also, a generalized dispersion relation is derived in terms of the solutions for the longitudinal electric and magnetic fields. For the case of a step-index fiber, either isotropic or uniaxial, these differential equations can be solved exactly in terms of Bessel functions. For the cases of an istropic graded-index and a uniaxial graded-index fiber, a solution using the Wentzel, Krammers and Brillouin (WKB) approximation technique is shown. Results for some particular permittivity profiles are presented. Also the WKB solutions is compared with the vector solution found by Kurtz and Streifer. For the matrix formulation it is shown that the tangential components of the electric and magnetic fields satisfy a system of four first-order differential equations which can be conveniently written in matrix form. For the special case of meridional modes, the system of equations splits into two systems of two equations. A general iterative technique, asymptotic partitioning of systems of equations, for solving systems of differential equations is presented. As a simple example, Bessel's differential equation is written in matrix form and is solved using this asymptotic technique. Low order solutions for particular examples of a biaxial and uniaxial graded-index fiber are presented. Finally numerical results obtained using the asymptotic technique are presented for particular examples of isotropic and uniaxial step-index fibers and isotropic, uniaxial and biaxial graded-index fibers.

Kawalko, Stephen F.↗
NASA NTRS1984

A spectral element method for fluid dynamics - Laminar flow in a channel expansion

A spectral element method that combines the generality of the finite element method with the accuracy of spectral techniques is proposed for the numerical solution of the incompressible Navier-Stokes equations. In the spectral element discretization, the computational domain is broken into a series of elements, and the velocity in each element is represented as a high-order Lagrangian interpolant through Chebyshev collocation points. The hyperbolic piece of the governing equations is then treated with an explicit collocation scheme, while the pressure and viscous contributions are treated implicitly with a projection operator derived from a variational principle. The implementation of the technique is demonstrated on a one-dimensional inflow-outflow advection-diffusion equation, and the method is then applied to laminar two-dimensional (separated) flow in a channel expansion. Comparisons are made with experiment and previous numerical work.

Patera, A. T.↗
NASA NTRS1990

A probabilistic approach to the dynamic analysis of ducts subjected to multibase harmonic and random excitation

The dynamic behavior of the discharge duct of the high-pressure oxidizer turbopump of a cryogenic rocket motor is investigated analytically. The probabilistic analysis program NESSUS (Numerical Evaluation of Stochastic Structures Under Stress; Cruse et al., 1988) is used to treat the uncertainties due to random and harmonic excitation (e.g., pump noise, pump-induced harmonics, and combustion noise), variations in engine inlet pressure, and changes in system damping. The load modeling procedure, the variation in power-spectral density in different zones of the engine structure, and the dynamic structural-analysis technique are described, and the numerical results of the NESSUS analysis are presented in extensive tables and graphs and discussed in detail.

Debchaudhury, Amit↗
NASA NTRS2001

Prediction of Business Jet Airloads Using The Overflow Navier-Stokes Code

The objective of this work is to evaluate the application of Navier-Stokes computational fluid dynamics technology, for the purpose of predicting off-design condition airloads on a business jet configuration in the transonic regime. The NASA Navier-Stokes flow solver OVERFLOW with Chimera overset grid capability, availability of several numerical schemes and convergence acceleration techniques was selected for this work. A set of scripts which have been compiled to reduce the time required for the grid generation process are described. Several turbulence models are evaluated in the presence of separated flow regions on the wing. Computed results are compared to available wind tunnel data for two Mach numbers and a range of angles-of-attack. Comparisons of wing surface pressure from numerical simulation and wind tunnel measurements show good agreement up to fairly high angles-of-attack.

Bounajem, Elias↗
NASA NTRS2000

Evaluation of Bogus Vortex Techniques with Four-Dimensional Variational Data Assimilation

The effectiveness of techniques for creating "bogus" vortices in numerical simulations of hurricanes is examined by using the Penn State/NCAR nonhydrostatic mesoscale model (MM5) and its adjoint system. A series of four-dimensional variational data assimilation (4-D VAR) experiments is conducted to generate an initial vortex for Hurricane Georges (1998) in the Atlantic Ocean by assimilating bogus sea-level pressure and surface wind information into the mesoscale numerical model. Several different strategies are tested for improving the vortex representation. The initial vortices produced by the 4-D VAR technique are able to reproduce many of the structural features of mature hurricanes. The vortices also result in significant improvements to the hurricane forecasts in terms of both intensity and track. In particular, with assimilation of only bogus sea-level pressure information, the response in the wind field is contained largely within the divergent component, with strong convergence leading to strong upward motion near the center. Although the intensity of the initial vortex seems to be well represented, a dramatic spin down of the storm occurs within the first 6 h of the forecast. With assimilation of bogus surface wind data only, an expected dominance of the rotational component of the wind field is generated, but the minimum pressure is adjusted inadequately compared to the actual hurricane minimum pressure. Only when both the bogus surface pressure and wind information are assimilated together does the model produce a vortex that represents the actual intensity of the hurricane and results in significant improvements to forecasts of both hurricane intensity and track.

Pu, Zhao-Xia↗