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 325 records · Page 18

Low-energy electron-atom bremsstrahlung

This paper extends recent work providing an elementary calculation of bremsstrahlung and opacity associated with the scattering of low-energy electrons by neutral atoms and molecules. The method applies when the scattering potential is 'hard' or when the collision time is short, applying the classical soft-photon emission probability formula for arbitrary bremsstrahlung photon energy. However, now, in addition to correcting the probability factor for finite bremsstrahlung photon energy, another factor corrects for the reduced phase space available to the outgoing electron. The bremsstrahlung cross section and opacity are then computed directly from the elastic scattering cross section, determined experimentally or calculated; a small (approximately 10 percent) correction is computed from the anisotropic term in the elastic scattering cross section. The opacity is evaluated for electron scattering by H, He, and H2 using experimentally determined values for the elastic scattering cross section, and is compared with more elaborate calculations. The agreement is good (within 10 percent), indicating an accuracy for the general method comparable to variations among the results of different elaborate theoretical computations. The agreement seems to validate the basic approximation of short collision time even at large bremsstrahlung photon energy for electron energies and temperatures up to a few eV.

Gould, R. J.↗

Phase equilibrium of methane and nitrogen at low temperatures - Application to Titan

Since the vapor phase composition of Titan's methane-nitrogen lower atmosphere is uniquely determined as a function of the Gibbs phase rule, these data are presently computed via integration of the Gibbs-Duhem equation. The thermodynamic consistency of published measurements and calculations of the vapor phase composition is then examined, and the saturated mole fraction of gaseous methane is computed as a function of altitude up to the 700-mbar level. The mole fraction is found to lie approximately halfway between that computed from Raoult's law, for a gas in equilibrium with an ideal solution of liquid nitrogen and methane, and that for a gas in equilibrium with pure liquid methane.

Kouvaris, Louis C.↗

Practical Calculation of Second-order Supersonic Flow past Nonlifting Bodies of Revolution

Calculation of second-order supersonic flow past bodies of revolution at zero angle of attack is described in detail, and reduced to routine computation. Use of an approximate tangency condition is shown to increase the accuracy for bodies with corners. Tables of basic functions and standard computing forms are presented. The procedure is summarized so that one can apply it without necessarily understanding the details of the theory. A sample calculation is given, and several examples are compared with solutions calculated by the method of characteristics.

Van Dyke, Milton D↗

Sensitivity analysis and approximation methods for general eigenvalue problems

Optimization of dynamic systems involving complex non-hermitian matrices is often computationally expensive. Major contributors to the computational expense are the sensitivity analysis and reanalysis of a modified design. The present work seeks to alleviate this computational burden by identifying efficient sensitivity analysis and approximate reanalysis methods. For the algebraic eigenvalue problem involving non-hermitian matrices, algorithms for sensitivity analysis and approximate reanalysis are classified, compared and evaluated for efficiency and accuracy. Proper eigenvector normalization is discussed. An improved method for calculating derivatives of eigenvectors is proposed based on a more rational normalization condition and taking advantage of matrix sparsity. Important numerical aspects of this method are also discussed. To alleviate the problem of reanalysis, various approximation methods for eigenvalues are proposed and evaluated. Linear and quadratic approximations are based directly on the Taylor series. Several approximation methods are developed based on the generalized Rayleigh quotient for the eigenvalue problem. Approximation methods based on trace theorem give high accuracy without needing any derivatives. Operation counts for the computation of the approximations are given. General recommendations are made for the selection of appropriate approximation technique as a function of the matrix size, number of design variables, number of eigenvalues of interest and the number of design points at which approximation is sought.

Murthy, D. V.↗

An Abstract Interpretation Framework for the Round-Off Error Analysis of Floating-Point Programs

This paper presents an abstract interpretation framework for the round-off error analysis of floating-point programs. This framework defines a parametric abstract analysis that computes, for each combination of ideal and floating-point execution path of the program, a sound over-approximation of the accumulated floating-point round-off error that may occur. In addition, a Boolean expression that characterizes the input values leading to the computed error approximation is also computed. An abstraction on the control flow of the program is proposed to mitigate the explosion of the number of elements generated by the analysis. Additionally, a widening operator is defined to ensure the convergence of recursive functions and loops. An instantiation of this framework is implemented in the prototype tool PRECiSA that generates formal proof certificates stating the correctness of the computed round-off errors.

Titolo, Laura↗

Use of a computer model in the understanding of erythropoietic control mechanisms

During an eight-week visit approximately 200 simulations using the computer model for the regulation of erythopoiesis were carries out in four general areas: with the human model simulating hypoxia and dehydration, evaluation of the simulation of dehydration using the mouse model. The experiments led to two considerations for the models. Firstly, a direct relationship between erythropoietin concentration and bone marrow sensitivity to the hormone and, secondly, a partial correction of tissue hypoxia prior to compensation by an increased hematocrit. This latter change in particular produced a better simuation of the effects of hypoxia on plasma erythropoietin concentrations.

Dunn, C. D. R.↗

A vectorization of the Hess McDonnell Douglas potential flow program NUED for the STAR-100 computer

The computer program NUED for analyzing potential flow about arbitrary three dimensional lifting bodies using the panel method was modified to use vector operations and run on the STAR-100 computer. A high speed of computation and ability to approximate the body surface with a large number of panels are characteristics of NUEDV. The new program shows that vector operations can be readily implemented in programs of this type to increase the computational speed on the STAR-100 computer. The virtual memory architecture of the STAR-100 facilitates the use of large numbers of panels to approximate the body surface.

Boney, L. R.↗

Three-dimensional inelastic approximate analysis code (MOMM)

The Mechanics of Materials Model (MOMM) is one of a series of new stand-alone three dimensional nonlinear structural analysis codes. Incorporation of a general purpose finite element computer code into the hot section design process was severely limited by the high costs involved. MOMM is a stiffness method finite element code that uses an internally generated network of beams to characterize hot section component behavior. The method was proposed as a fast, easy to use, computationally efficient tool for approximate analyses. MOMM incorporates a wide variety of analysis capabilities, material models, and load type specifiers instrumental for the analysis of hot section components.

Meister, Jeffrey P.↗

HAL/SM language specification

A programming language is presented for the flight software of the NASA Space Shuttle program. It is intended to satisfy virtually all of the flight software requirements of the space shuttle. To achieve this, it incorporates a wide range of features, including applications-oriented data types and organizations, real time control mechanisms, and constructs for systems programming tasks. It is a higher order language designed to allow programmers, analysts, and engineers to communicate with the computer in a form approximating natural mathematical expression. Parts of the English language are combined with standard notation to provide a tool that readily encourages programming without demanding computer hardware expertise. Block diagrams and flow charts are included. The semantics of the language is discussed.

Williams, G. P. W., Jr.↗

Exploring Quantum State Preparation Using Tensor Networks and Sparse Wavefunction Simulations

The variational quantum eigenvalue solver is a powerful hybrid quantum-classical approach that has been suggested as a candidate method to run on near-term quantum hardware for computing ground state electronic energies of molecular systems. However, even for small molecules, the number of variational parameters and qubits required to minimize the electronic energy is beyond the reach of current quantum computers except for small basis sets. We explore a new paradigm for state preparation where we test how much of the optimization can be approximately prepared with classical computers to reduce the number of optimization steps performed using a quantum device. By adapting a recent algorithm for the factorized form of the UCC ansatz, we can study molecular electronic structure problems with up to 64 qubits. In addition, we also test a related approach of using tensor networks to optimize quantum circuits in order to benchmark various lattice models. We present results using these approaches and discuss strategies for incorporating these ideas into variational algorithms involving near-term quantum computers. Our results help demonstrate the strength of the UCC ansatz and address pressing questions about optimal initial parameterizations and circuit construction.

quantum computing↗

An Investigation of Multipath Effects on the GPS System During Auto-Rendezvous and Capture

The proposed use of a Cargo Transport Vehicle (CTV) to carry hardware to the Space Station Freedom (SSF) during the construction phase of the SSF project requires remote maneuvering of the CTV. The CTV is not a manned vehicle. Obtaining the relative positions of the CTV and SSF for remote auto-rendezvous and capture (AR&C) scenarios will rely heavily on the Global Positioning System (GPS). The GPS system is expected to guide the CTV up to a distance of 100 to 300 meters from the SSF. At some point within this range, an optical docking system will take over the remote guidance for capture. During any remote guidance by GPS it is possible that significant multipath signals may be caused by large objects in the vicinity of the module being remotely guided. This could alter the position obtained by the GPS system from the actual position. Due to the nature of the GPS signals, it has been estimated that if the difference in distance between the Line of Sight (LOS) path and the multipath is greater than 300 meters, the GPS system is capable of discriminating between the direct signal and the reflected (or multipath) signal. However, if the path difference is less than 300 meters, one must be concerned. This report details the work accomplished by the Electromagnetic Simulations Laboratory at Marquette University over the period December 1993 to May 1995. This work is an investigation of the strength and phase of a multipath signal arriving at the CTV relative to the direct or line of sight (LOS) signal. The signal originates at a GPS satellite in half geo-stationary orbit and takes two paths to the CTV: (1) the direct or LOS path from the GPS satellite to the CTV; and (2) a scattered path from the GPS satellite to the SSF module and then to the CTV. The scattering from a cylinder has been computed using the physical optics approximation for the current. No other approximations or assumptions have been made including no assumptions regarding the far field or Fresnel field approximations. The integrations required to obtain the scattered field have been computed numerically using an N dimensional Romberg integration. The total scattered electric field is then projected onto the RCP component in the direction of propagation only. The direct or line of sight signal is then used to compute the relative strength and phase of the scattered field. The trajectory of the CTV has been parameterized into 4,214 points that are calculated for each of the geometries investigated. The motion of the CTV between points is small enough for the magnitude data (dB down from direct signal) to appear very smooth; however, because of the distances and wavelengths involved, the phase of the scattered field relative to the direct signal varies very rapidly.

Richie, James E.↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper presents the results of an ongoing study into the linear stability characteristics of the boundary layer flow over the common research model with natural laminar flow (CRMNLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work involving parallel stability computations of a boundary layer flow based on the conical flow approximation has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien- Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. Here, we examine the effects of the simplifying approximations in both basic state computation and the stability analysis, with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, the basic states are computed by using full Navier-Stokes equations and the stability analysis is performed by using a nonorthogonal coordinate system that allows a clear distinction between planar TS and CF instabilities. Furthermore, the effects of curvature and nonparallel mean flow have been included in the stability computations based on the parabolized stability equations (PSE). The fully turbulent Reynolds-Averaged-Navier-Stokes (RANS) mean flow solutions show good agreement with the measured wall pressure distribution. Viscous-inviscid interactive effects are observed to be important because the shock fronts along the suction surface are influenced by the imposed transition front. The stability results confirm the previous findings related to TS amplification within the inboard region of the wing and the dominance of stationary CF modes in the outboard region. However, given the close proximity of the measured transition front and the dual shock system within the outer part of the wing, the onset of transition may well be shock limited within the outboard region. In general, the transition criterion based on the dual N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at Re MAC = 15 million and AoA = 1.5, 2 degrees; however, the low values of the correlating N-factors at Re MAC = 17.5 million support the hypothesis that the measured transition at the higher Reynolds number may have been strongly influenced by the merging of turbulent wedges that originate from surface imperfections near the leading edge.

Boundary layer transition↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper presents the results of an ongoing study into the linear stability characteristics of the boundary layer flow over the common research model with natural laminar flow (CRMNLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work involving parallel stability computations of a boundary layer flow based on the conical flow approximation has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien- Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. Here, we examine the effects of the simplifying approximations in both basic state computation and the stability analysis, with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, the basic states are computed by using full Navier-Stokes equations and the stability analysis is performed by using a nonorthogonal coordinate system that allows a clear distinction between planar TS and CF instabilities. Furthermore, the effects of curvature and nonparallel mean flow have been included in the stability computations based on the parabolized stability equations (PSE). The fully turbulent Reynolds-Averaged-Navier-Stokes (RANS) mean flow solutions show good agreement with the measured wall pressure distribution. Viscous-inviscid interactive effects are observed to be important because the shock fronts along the suction surface are influenced by the imposed transition front. The stability results confirm the previous findings related to TS amplification within the inboard region of the wing and the dominance of stationary CF modes in the outboard region. However, given the close proximity of the measured transition front and the dual shock system within the outer part of the wing, the onset of transition may well be shock limited within the outboard region. In general, the transition criterion based on the dual N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at Re MAC = 15 million and AoA = 1.5, 2 degrees; however, the low values of the correlating N-factors at Re MAC = 17.5 million support the hypothesis that the measured transition at the higher Reynolds number may have been strongly influenced by the merging of turbulent wedges that originate from surface imperfections near the leading edge.

Boundary layer transition↗

Hybrid transfinite element modeling/analysis of nonlinear heat conduction problems involving phase change

The present paper describes the applicability of hybrid transfinite element modeling/analysis formulations for nonlinear heat conduction problems involving phase change. The methodology is based on application of transform approaches and classical Galerkin schemes with finite element formulations to maintain the modeling versatility and numerical features for computational analysis. In addition, in conjunction with the above, the effects due to latent heat are modeled using enthalpy formulations to enable a physically realistic approximation to be dealt computationally for materials exhibiting phase change within a narrow band of temperatures. Pertinent details of the approach and computational scheme adapted are described in technical detail. Numerical test cases of comparative nature are presented to demonstrate the applicability of the proposed formulations for numerical modeling/analysis of nonlinear heat conduction problems involving phase change.

Tamma, Kumar K.↗

Computational fluid dynamics and aerothermodynamics

Approximations applicable to the radiating, reacting, and conducting stagnation region of a hypervelocity vehicle were incorporated into a method for rapidly obtaining approximate solutions. This solution utilizes a coordinate system based upon the origin of the radiative losses and includes in a phenomenologically correct manner the effects of chemical and thermal nonequilibrium, and nonequilibrium, nongray radiative transfer. Results were presented which demonstrate the usefulness of the method and indicate which radiation parameters require further study and definition. Excellent comparisons were obtained with published results for the Fire2 data. An axisymmetric nonequilibrium inverse method was modified and extended and used to investigate and compare various vibration dissociation chemistry coupling models and radiative heat transfer approximations. The similarities, differences, and consequences of using these models in the Aero-assist Orbital Transfer Vehicles flight regime will be discussed.

Carlson, Leland A.↗

Nongray radiative gas analyses using the S-N discrete ordinates method

The S-N discrete ordinates method is applied to analyze radiative heat transfer in nongray gases. Spectral correlation between the terms in the equation of transfer is considered for black or nearly nonreflecting walls. Formulations to apply the S-N method using a narrow-band or the exponential wide-band model are presented. The net radiative wall heat fluxes and the radiative source distributions are obtained for uniform, parabolic, and boundary layer-type temperature profiles, as well as for a parabolic concentration profile. The narrow- and wide-band nongray solutions are compared with gray-band approximations using the same band models. The computational speed of the gray-band approximation is obtained at the expense of accuracy in the internal fluxes and radiative source distributions. The wall radiative flux predictions by the gray-band approximation are satisfactory.

Kim, T. K.↗