Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “approximate computing”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

CFD Simulations in Support of Shuttle Orbiter Contingency Abort Aerodynamic Database Enhancement

Modern Computational Fluid Dynamics (CFD) techniques were used to compute aerodynamic forces and moments of the Space Shuttle Orbiter in specific portions of contingency abort trajectory space. The trajectory space covers a Mach number range of 3.5-15, an angle-of-attack range of 20deg-60deg, an altitude range of 100-190 kft, and several different settings of the control surfaces (elevons, body flap, and speed brake). Presented here are details of the methodology and comparisons of computed aerodynamic coefficients against the values in the current Orbiter Operational Aerodynamic Data Book (OADB). While approximately 40 cases have been computed, only a sampling of the results is provided here. The computed results, in general, are in good agreement with the OADB data (i.e., within the uncertainty bands) for almost all the cases. However, in a limited number of high angle-of-attack cases (at Mach 15), there are significant differences between the computed results, especially the vehicle pitching moment, and the OADB data. A preliminary analysis of the data from the CFD simulations at Mach 15 shows that these differences can be attributed to real-gas/Mach number effects. The aerodynamic coefficients and detailed surface pressure distributions of the present simulations are being used by the Shuttle Program in the evaluation of the capabilities of the Orbiter in contingency abort scenarios.

Papadopoulos, Periklis E.↗
NASA NTRS2012

Electron Stark Broadening Database for Atomic N, O, and C Lines

A database for efficiently computing the electron Stark broadening line widths for atomic N, O, and C lines is constructed. The line width is expressed in terms of the electron number density and electronatom scattering cross sections based on the Baranger impact theory. The state-to-state cross sections are computed using the semiclassical approximation, in which the atom is treated quantum mechanically whereas the motion of the free electron follows a classical trajectory. These state-to-state cross sections are calculated based on newly compiled line lists. Each atomic line list consists of a careful merger of NIST, Vanderbilt, and TOPbase line datasets from wavelength 50 nm to 50 micrometers covering the VUV to IR spectral regions. There are over 10,000 lines in each atomic line list. The widths for each line are computed at 13 electron temperatures between 1,000 K 50,000 K. A linear least squares method using a four-term fractional power series is then employed to obtain an analytical fit for each line-width variation as a function of the electron temperature. The maximum L2 error of the analytic fits for all lines in our line lists is about 5%.

Liu, Yen↗
NASA NTRS1979

Simulation of three-dimensional compressible viscous flow on the Illiac IV computer

Complicated three-dimensional viscous transonic flows about bodies at high angles of attack are solved on the Illiac IV computer. It is shown that certain approximate forms of the compressible Reynolds-averaged Navier-Stokes equations can be computed about realistic three-dimensional geometries with relative ease on the Illiac IV. The ease and efficiency with which this can be done depend on the approximations made in the basic equations, the choice of the numerical algorithm used for the solution, and the data-base system that controls the data management and identifies and manipulates the vectors. A pencil data-base system is found to be particularly suitable for the approximations and numerical method chosen to produce the results presented. In addition, some comparisons are made of computer predictions with experimental results for various lows about hemisphere-cylinders in both subsonic and supersonic free streams. The same viscous model and numerical model are used, showing good qualitative agreement in the location of separation lines and pressure distributions.

Pulliam, T. H.↗
NASA NTRS2001

Multidisciplinary Design Optimization of a Full Vehicle with High Performance Computing

Multidisciplinary design optimization (MDO) of a full vehicle under the constraints of crashworthiness, NVH (Noise, Vibration and Harshness), durability, and other performance attributes is one of the imperative goals for automotive industry. However, it is often infeasible due to the lack of computational resources, robust simulation capabilities, and efficient optimization methodologies. This paper intends to move closer towards that goal by using parallel computers for the intensive computation and combining different approximations for dissimilar analyses in the MDO process. The MDO process presented in this paper is an extension of the previous work reported by Sobieski et al. In addition to the roof crush, two full vehicle crash modes are added: full frontal impact and 50% frontal offset crash. Instead of using an adaptive polynomial response surface method, this paper employs a DOE/RSM method for exploring the design space and constructing highly nonlinear crash functions. Two NMO strategies are used and results are compared. This paper demonstrates that with high performance computing, a conventionally intractable real world full vehicle multidisciplinary optimization problem considering all performance attributes with large number of design variables become feasible.

Yang, R. J.↗
NASA NTRS1974

Efficient stream distributions in radiative transfer theory

This paper discusses a new, computationally-efficient method for approximating the integro-differential equation of radiative transfer with a finite set of coupled differential equations for discrete streams. The method uses recommended spatial distributions of streams that are quite different from those typically used in that they are based on the symmetry of several regular Platonic solids. To facilitate the use of such distributions, an explicit, one-parameter relationship between the physical radiance and the abstract stream is formulated. The parameter is used to determine the minimum number of streams required in the radiative transfer model. Accuracy and computational efficiency are shown to be served best by choosing a stream distribution that is invariant to a large number of three space rotations. For various values of the above-mentioned parameter, the resulting recommended stream distribution is shown to be more computationally efficient than more conventional stream distributions. Finally, the incorporation of polarization in the stream definition is described.

Whitney, C.↗
NASA NTRS2016

The Transition from Thick to Thin Plate Wake Physics: Whither Vortex Shedding?

The near and very near wake of a flat plate with a circular trailing edge is investigated with data from direct numerical simulations. Computations were performed for six different combinations of the Reynolds numbers based on plate thickness (D) and boundary layer momentum thickness upstream of the trailing edge (theta). Unlike the case of the cylinder, these Reynolds numbers are independent parameters for the flat plate. The separating boundary layers are turbulent in all the cases investigated. One objective of the study is to understand the changes in the wake vortex shedding process as the plate thickness is reduced (increasing theta/D). The value of D varies by a factor of 16 and that of theta by approximately 5 in the computations. Vortex shedding is vigorous in the low theta/D cases with a substantial decrease in shedding intensity in the large theta/D cases. Other shedding characteristics are also significantly altered with increasing theta/D. A visualization of the shedding process in the different cases is provided and discussed. The basic shedding mechanism is explored in depth. The effect of changing theta/D on the time-averaged, near-wake velocity statistics is also discussed. A functional relationship between the shedding frequency and the Reynolds numbers mentioned above is obtained.

CFD↗
NASA NTRS1987

Numerical simulation of three dimensional transonic flows

The three-dimensional flow over a projectile has been computed using an implicit, approximately factored, partially flux-split algorithm. A simple composite grid scheme has been developed in which a single grid is partitioned into a series of smaller grids for applications which require an external large memory device such as the SSD of the CRAY X-MP/48, or multitasking. The accuracy and stability of the composite grid scheme has been tested by numerically simulating the flow over an ellipsoid at angle of attack and comparing the solution with a single grid solution. The flowfield over a projectile at M = 0.96 and 4 deg angle-of-attack has been computed using a fine grid, and compared with experiment.

Sahu, Jubaraj↗
NASA NTRS2015

A Module for Assimilating Hyperspectral Infrared Retrieved Profiles into the Gridpoint Statistical Interpolation System for Unique Forecasting Applications

Hyperspectral infrared sounder radiance data are assimilated into operational modeling systems however the process is computationally expensive and only approximately 1% of available data are assimilated due to data thinning as well as the fact that radiances are restricted to cloud-free fields of view. In contrast, the number of hyperspectral infrared profiles assimilated is much higher since the retrieved profiles can be assimilated in some partly cloudy scenes due to profile coupling other data, such as microwave or neural networks, as first guesses to the retrieval process. As the operational data assimilation community attempts to assimilate cloud-affected radiances, it is possible that the use of retrieved profiles might offer an alternative methodology that is less complex and more computationally efficient to solve this problem. The NASA Short-term Prediction Research and Transition (SPoRT) Center has assimilated hyperspectral infrared retrieved profiles into Weather Research and Forecasting Model (WRF) simulations using the Gridpoint Statistical Interpolation (GSI) System. Early research at SPoRT demonstrated improved initial conditions when assimilating Atmospheric Infrared Sounder (AIRS) thermodynamic profiles into WRF (using WRF-Var and assigning more appropriate error weighting to the profiles) to improve regional analysis and heavy precipitation forecasts. Successful early work has led to more recent research utilizing WRF and GSI for applications including the assimilation of AIRS profiles to improve WRF forecasts of atmospheric rivers and assimilation of AIRS, Cross-track Infrared and Microwave Sounding Suite (CrIMSS), and Infrared Atmospheric Sounding Interferometer (IASI) profiles to improve model representation of tropopause folds and associated non-convective wind events. Although more hyperspectral infrared retrieved profiles can be assimilated into model forecasts, one disadvantage is the retrieved profiles have traditionally been assigned the same error values as the rawinsonde observations when assimilated with GSI. Typically, satellitederived profile errors are larger and more difficult to quantify than traditional rawinsonde observations (especially in the boundary layer), so it is important to appropriately assign observation errors within GSI to eliminate potential spurious innovations and analysis increments that can sometimes arise when using retrieved profiles. The goal of this study is to describe modifications to the GSI source code to more appropriately assimilate hyperspectral infrared retrieved profiles and outline preliminary results that show the differences between a model simulation that assimilated the profiles as rawinsonde observations and one that assimilated the profiles in a module with the appropriate error values.

Berndt, Emily↗
NASA NTRS1990

Numerical simulation of three-dimensional transonic flows

The three-dimensional flow over a projectile has been computed using an implicit, approximately factored, partially flux-split algorithm. A simple composite grid scheme has been developed in which a single grid is partitioned into a series of smaller grids for applications which require an external large memory device such as the SSD of the CRAY X-MP/48 or multi-tasking. The accuracy and stability of the composite grid scheme have been tested by numerically simulating the flow over an ellipsoid at an angle of attack and comparing the solution with a single-grid solution. The flow field over a projectile at M = 0.96 and 1.1, and 4-deg angle of attack has been computed using a fine grid and compared with experiment.

Sahu, Jubaraj↗
NASA NTRS2020

Reaching for 20 Years with the IMERG Multi-Satellite Products

The latest releases of Global Precipitation Measurement (GPM) mission products cap five years of vigorous development cycle since the launch of the GPM Core Observatory, and these now provide datasets that are relatively homogeneous across the joint Tropical Rainfall Measuring Mission (TRMM) and GPM eras. Version 06 of the U.S. GPM team's Integrated Multi-satellitE Retrievals for GPM (IMERG) merged precipitation product enforces a consistent intercalibration for all precipitation products computed from individual satellites with the TRMM and GPM Core Observatory sensors as the TRMM- and GPM era calibrators, respectively, and incorporates monthly surface gauge data. The basic IMERG algorithm now features precipitation motion vectors (used to drive the Lagrangian interpolation, or "morphing") that arecomputed by tracking vertically integrated vapor fields analyzed in MERRA2 and GEOS5. This innovation provides globally complete coverage, expanding IMERG's coverage beyond the 60°N-S latitude band provided by IR-based vectors, although we continue to mask out precipitation over snowy/icy surfaces as unreliable. A second innovation is the Quality Index (QI) data field. The half-hourly QI is taken as the approximate Kalman Filter correlation computed in the morphing calculation.

Huffman, George J.↗
NASA NTRS1983

Electromagnetic backscattering from a layer of vegetation - A discrete approach

Electromagnetic backscattering from a layer of vegetation over a flat lossy ground is studied. The vegetated region is composed of leaves which are modeled by lossy circular dielectric discs with prescribed orientation statistics. The mean field in the layer is computed by using the Foldy approximation and a two-variable technique which is valid for small fractional volume. Simple expressions are given for the skin depth of the mean wave. Following this, the distorted Born approximation is used to calculate the co-polarized and cross polarized backscattering coefficients. These are then computed as a function of angle of incidence for different layer depths and leaf configurations.

Lang, R. H.↗
NASA NTRS1988

Mars rover technology development requirements

The technology development requirements for various Mars rover range capabilities are discussed, focusing on local navigation of the rover. The capabilities of two methods are compared. In one method, operators on the earth view stereo pictures sent by the rover and determine short traverse paths which the rover follows. The other method achieves more autonomous capability by using computer vision from orbital imagery with approximate long routes commanded from earth. The locomotion, navigation, ground operations, computation, power, thermal control, communications, sample acquisition, and analysis and preservation requirements are examined.

Bedard, Roger↗
NASA NTRS1988

Interaction and merging of vortex filaments

The asymptotic solutions of Navier-Stokes equations for vortex filaments of finite strength with small effective vortical cores are summarized with special emphasis placed on the physical meaning and the practical limit to the applicability of the asymptotic solution. Finite-difference solutions of Navier-Stokes equations for the marging of the filament(s) are described with a focus on the development of the approximate boundary conditions for the computational domain. An efficiency study employing a model problem is used to assess the advantages of the present approximate boundary condition method over previously used techniques. Applications of the present method are presented for the motion and decay of a 3:1 elliptic vortex ring, and for the merging process of a pair of coaxial vortex rings. A numerical procedure for the problem of local merging of vortex filaments, which requires the asymptotic analysis as well as the numerical Navier-Stokes solver, is also presented.

Liu, C. H.↗
NASA NTRS1989

Shape sensitivity analysis of flutter response of a laminated wing

A method is presented for calculating the shape sensitivity of a wing aeroelastic response with respect to changes in geometric shape. Yates' modified strip method is used in conjunction with Giles' equivalent plate analysis to predict the flutter speed, frequency, and reduced frequency of the wing. Three methods are used to calculate the sensitivity of the eigenvalue. The first method is purely a finite difference calculation of the eigenvalue derivative directly from the solution of the flutter problem corresponding to the two different values of the shape parameters. The second method uses an analytic expression for the eigenvalue sensitivities of a general complex matrix, where the derivatives of the aerodynamic, mass, and stiffness matrices are computed using a finite difference approximation. The third method also uses an analytic expression for the eigenvalue sensitivities, but the aerodynamic matrix is computed analytically. All three methods are found to be in good agreement with each other.

Kapania, Rakesh K.↗
NASA NTRS1994

A Full Navier-Stokes Analysis of Subsonic Diffuser of a Bifurcated 70/30 Supersonic Inlet for High Speed Civil Transport Application

A full Navier-Stokes analysis was performed to evaluate the performance of the subsonic diffuser of a NASA Lewis Research Center 70/30 mixed-compression bifurcated supersonic inlet for high speed civil transport application. The PARC3D code was used in the present study. The computations were also performed when approximately 2.5 percent of the engine mass flow was allowed to bypass through the engine bypass doors. The computational results were compared with the available experimental data which consisted of detailed Mach number and total pressure distribution along the entire length of the subsonic diffuser. The total pressure recovery, flow distortion, and crossflow velocity at the engine face were also calculated. The computed surface ramp and cowl pressure distributions were compared with experiments. Overall, the computational results compared well with experimental data. The present CFD analysis demonstrated that the bypass flow improves the total pressure recovery and lessens flow distortions at the engine face.

Kapoor, Kamlesh↗
NASA NTRS1972

Optimizing a liquid propellant rocket engine with an automated combustor design code (AUTOCOM)

A procedure for automatically designing a liquid propellant rocket engine combustion chamber in an optimal fashion is outlined. The procedure is contained in a digital computer code, AUTOCOM. The code is applied to an existing engine, and design modifications are generated which provide a substantial potential payload improvement over the existing design. Computer time requirements for this payload improvement were small, approximately four minutes in the CDC 6600 computer.

Hague, D. S.↗
NASA NTRS1973

On singlet s-wave electron-hydrogen scattering.

Discussion of various zeroth-order approximations to s-wave scattering of electrons by hydrogen atoms below the first excitation threshold. The formalism previously developed by the author (1967, 1968) is applied to Feshbach operators to derive integro-differential equations, with the optical-potential set equal to zero, for the singlet and triplet cases. Phase shifts of s-wave scattering are computed in the zeroth-order approximation of the Feshbach operator method and in the static-exchange approximation. It is found that the convergence of numerical computations is faster in the former approximation than in the latter.

Madan, R. N.↗
NASA NTRS1990

A review of reduced Navier-Stokes computations for compressible viscous flows

A reduced form of the Navier-Stokes equations, defined by a single composite of the Euler, boundary layer, and triple deck approximations, is considered for the computation of viscous interacting flows. Global pressure or pseudopotential relaxation methods are integrated with coupled sparse matrix direct solvers or coupled strongly implicit ILU algorithms to efficiently capture sharp shocks and regions of recirculation. Solutions are obtained for a variety of 2D and 3D geometries and for Mach numbers (M) spanning the range from incompressible (M = 0) to supersonic (M = 6) flow.

Rubin, S. G.↗