Search NASA⌕ Search

SEARCH · Search NASA

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 1,441 records · Page 80

Improvements in Iterative Convergence of FUN3D Solutions

This paper presents a hierarchical adaptive nonlinear iteration method (HANIM) implemented in the NASA computational fluid dynamics code, FUN3D, to improve robustness and computational efficiency. In contrast to the legacy FUN3D iterative solver that relies on an approximate Jacobian, a simple multicolor Gauss-Seidel point-implicit iteration scheme, and linear Courant-Friedrichs-Lewy number (CFL) ramping, HANIM is based upon a hierarchy of modules including preconditioner, generalized conjugate residual, realizability check, nonlinear control, and CFL adaption modules. HANIM performance is systematically compared with the performance of the legacy solver of FUN3D and a baseline solver based on a preconditioner alone. Iterative solutions are compared for three benchmark cases: a subsonic separated flow around a hemisphere cylinder, a supersonic flow through a long duct, and a subsonic flow over the NASA wing-fuselage juncture model. Two Reynolds-averaged Navier-Stokes turbulence models are used in these computations, namely, the negative variant of the linear one-equation Spalart-Allmaras model and its nonlinear extension based on quadratic constitutive relations.

CFD↗

Determination of wind from NIMBUS 6 satellite sounding data

Objective methods of computing upper level and surface wind fields from NIMBUS 6 satellite sounding data are developed. These methods are evaluated by comparing satellite derived and rawinsonde wind fields on gridded constant pressure charts in four geographical regions. Satellite-derived and hourly observed surface wind fields are compared. Results indicate that the best satellite-derived wind on constant pressure charts is a geostrophic wind derived from highly smoothed fields of geopotential height. Satellite-derived winds computed in this manner and rawinsonde winds show similar circulation patterns except in areas of small height gradients. Magnitudes of the standard deviation of the differences between satellite derived and rawinsonde wind speeds range from approximately 3 to 12 m/sec on constant pressure charts and peak at the jet stream level. Fields of satellite-derived surface wind computed with the logarithmic wind law agree well with fields of observed surface wind in most regions. Magnitudes of the standard deviation of the differences in surface wind speed range from approximately 2 to 4 m/sec, and satellite derived surface winds are able to depict flow across a cold front and around a low pressure center.

Carle, W. E.↗

Exact and Approximate Probabilistic Symbolic Execution

Probabilistic software analysis seeks to quantify the likelihood of reaching a target event under uncertain environments. Recent approaches compute probabilities of execution paths using symbolic execution, but do not support nondeterminism. Nondeterminism arises naturally when no suitable probabilistic model can capture a program behavior, e.g., for multithreading or distributed systems. In this work, we propose a technique, based on symbolic execution, to synthesize schedulers that resolve nondeterminism to maximize the probability of reaching a target event. To scale to large systems, we also introduce approximate algorithms to search for good schedulers, speeding up established random sampling and reinforcement learning results through the quantification of path probabilities based on symbolic execution. We implemented the techniques in Symbolic PathFinder and evaluated them on nondeterministic Java programs. We show that our algorithms significantly improve upon a state-of- the-art statistical model checking algorithm, originally developed for Markov Decision Processes.

Symbolic Execution↗

On the use of the consistent mass matrix in flexible multibody systems

Two different methods are used in finite element structural analysis to define element mass distribution: consistent and lumped mass approximations. The lumped mass method is currently favored in multibody dynamics software packages because of its computational efficiency. This paper, however, presents a multibody dynamics formulation based on the consistent mass approach which should prove to be as computationally efficient as the lumped mass approach. It is shown that when component flexibility is modelled using modal coordinates, the equations of motion, as formulated by these two methods, are very similar in appearance, thereby facilitating their comparison. It is also shown, through simulation of a simple multibody system, that a system formulated with the lumped mass approximation will often require more finite elements and smaller integration step sizes than the corresponding consistent mass system.

Anthony, Tobin C.↗

The development and flight test of a deployable precision landing system for spacecraft recovery

A joint NASA Dryden Flight Research Facility and Johnson Space Center program was conducted to determine the feasibility of the autonomous recovery of a spacecraft using a ram-air parafoil system for the final stages of entry from space that included a precision landing. The feasibility of this system was studied using a flight model of a spacecraft in the generic shape of a flattened biconic which weighed approximately 150 lb and was flown under a commercially available, ram-air parachute. Key elements of the vehicle included the Global Positioning System guidance for navigation, flight control computer, ultrasonic sensing for terminal altitude, electronic compass, and onboard data recording. A flight test program was used to develop and refine the vehicle. This vehicle completed an autonomous flight from an altitude of 10,000 ft and a lateral offset of 1.7 miles which resulted in a precision flare and landing into the wind at a predetermined location. At times, the autonomous flight was conducted in the presence of winds approximately equal to vehicle airspeed. Several techniques for computing the winds postflight were evaluated. Future program objectives are also presented.

Sim, Alex G.↗

The Titan haze revisted: Magnetospheric energy sorces quantitative tholin yields

We present laboratory measurements of the radiation yields of complex organic solids produced from N2/CH4 gas mixtures containing 10 or 0.1% CH4. These tholins are thought to resemble organic aerosols produced in the atmospheres of Titan, Pluto, and Triton. The tholin yields are large compared to the total yield of gaseous products: nominally, 13 (C + N)/100 eV for Titan tholin and 2.1 (C + N)/100 eV for Triton tholin. High-energy magnetospheric electrons responsible for tholin production represents a class distinct from the plasma electrons considered in models of Titan's aiglow. Electrons with E greater than 20 keV provide an energy flux approximately 1 x 10(exp -2) erg/cm/sec, implying from our measured tholin yields a mass flux of 0.5 to 4.0 x 10(exp -14) g/sq cm/sec of tholin. (The corresponding thickness of the tholin sedimentary column accumulated over 4 Gyr on Titan's surface is 4 to 30 m). This figure is in agreement with required mass fluxes computed from recent radiative transfer and sedimentation models. If, however, theses results, derived from experiments at approximately 2 mb, are applied to lower pressure levels toward peak auroral electron energy deposition and scaled with pressure as the gas-phase organic yields, the derived tholin mass flux is at least an order of magnitude less. We attrribute this difference to the fact that tholin synthesis occurs well below the level of maximum electron energy depositon and to possible contributions to tholis from UV-derived C2-hydrocarbons. We conclude that Tita tholin, produced by magnetospheric electrons, is alone sufficient to supply at least a significant fraction of Titan's haze-a result consistent with the fact that the optical properties of Titan tholin, among all proposed material, are best at reproducing Titan's geometric albedo spectrum from near UV to mid-IR in light-scattering models.

Thompson, W. Reid↗

Design of the Cruise and Flap Airfoil for the X-57 Maxwell Distributed Electric Propulsion Aircraft

A computational and design study on an airfoil and high-lift flap for the X-57 Maxwell Distributed Electric Propulsion (DEP) testbed aircraft was conducted. The aircraft wing sizing study resulted in a wing area of 66.67 sq ft and aspect ratio of 15 with a design requirement of V(stall) = 58 KEAS, at a gross weight of 3,000 lb. To meet this goal an aircraft C(L,max) of 4.0 was required. The design cruise condition is 150 KTAS at 8,000 ft. This resulted in airfoil requirements of c(l) is approximately 0.90 for the cruise condition at Re = 2.35 x 10 (exp 6). A flapped airfoil with a c(l,max) of approximately 2.5 or greater, at Re = 1.0 x 10 (exp 6), was needed to have enough lift to meet the stall requirement with the DEP system. MSES computational analyses were conducted on the GAW-1, GAW-2, and the NACA 5415 airfoil sections, however they had limitations in either high drag or low c(l,max) on the cruise airfoil, which was the impetus for a new design. A design was conducted to develop a low drag airfoil for the X-57 cruise conditions with high c(l,max). The final design was the GNEW5BP93B airfoil with a minimum drag coefficient of c(d) = 0.0053 at c(l) = 0.90 and achieved laminar flow back to 69% chord on the upper surface and 62% chord on the lower surface. With fully turbulent flow, the drag increases to c(d) = 0.0120. The predicted maximum lift with turbulent flow is a c(l,max) of 1.95 at alpha = 19 deg. The airfoil is characterized by relatively flat pressure gradient regions on both surfaces at alpha = 0 deg, and aft camber to get extra lift out of the lower surface concave region. A 25% chord slotted flap was designed and analyzed with MSES for a 30 deg flap deflection. Additional 30 deg and 40 deg flap deflection analyses for two flap positions were conducted with USM3D using several turbulence models, for two angles of attack, to assess near c(l,max) with varied flap position. The maximum c(l) varied between 2.41 and 3.35. An infinite-span powered high-lift study was conducted on a GAW-1 constant chord 40 deg flapped airfoil section with FUN3D to quantify the airfoil lift increment that can be expected from a DEP system. The 16.7 hp/propeller blown wing increases the maximum C(L) from 3.45 to C(L) = 6.43, which is an effective q ratio of 1.86. This indicates that if the unblown high-lift flapped airfoil of the X-57 airplane achieves a c(l,max) of 2.78, then the high-lift augmentation blowing could yield a sectional lift coefficient of approximately 4.95 at c(l,max). Finally, a computational study was conducted with FUN3D on an infinite-span constant chord GAW-1 cruise airfoil to determine the impact of high-lift propeller diameter to wing chord ratio on the lift increment of the DEP system. A constant diameter propeller and nacelle size were used in the study. Three computational grids were made with airfoil chords of 0.5*chord, 1.0*chord, and 2.0*chord. Results of the propeller diameter to wing chord ratio study indicated that the blown to unblown C(L) ratio increased as the chord was decreased. However, because of the increase in relative size of the high-lift nacelle to the wing, which impacted wing lift performance, the study indicated that a propeller diameter to wing chord ratio of 1.0 gives the overall best maximum lift on the wing with the DEP system.

Viken, Jeffrey K.↗

Mir Cooperative Solar Array Flight Performance Data and Computational Analysis

The Mir Cooperative Solar Array (MCSA) was developed jointly by the United States (US) and Russia to provide approximately 6 kW of photovoltaic power to the Russian space station Mir. The MCSA was launched to Mir in November 1995 and installed on the Kvant-1 module in May 1996. Since the MCSA photovoltaic panel modules (PPMs) are nearly identical to those of the International Space Station (ISS) photovoltaic arrays, MCSA operation offered an opportunity to gather multi-year performance data on this technology prior to its implementation on ISS. Two specially designed test sequences were executed in June and December 1996 to measure MCSA performance. Each test period encompassed 3 orbital revolutions whereby the current produced by the MCSA channels was measured. The temperature of MCSA PPMs was also measured. To better interpret the MCSA flight data, a dedicated FORTRAN computer code was developed to predict the detailed thermal-electrical performance of the MCSA. Flight data compared very favorably with computational performance predictions. This indicated that the MCSA electrical performance was fully meeting pre-flight expectations. There were no measurable indications of unexpected or precipitous MCSA performance degradation due to contamination or other causes after 7 months of operation on orbit. Power delivered to the Mir bus was lower than desired as a consequence of the retrofitted power distribution cabling. The strong correlation of experimental and computational results further bolsters the confidence level of performance codes used in critical ISS electric power forecasting. In this paper, MCSA flight performance tests are described as well as the computational modeling behind the performance predictions.

Kerslake, Thomas W.↗

Computation of viscous flow in planar and axisymmetric ducts by an implicit marching procedure

A streamwise marching procedure, approximately 200 times faster than a full Navier-Stokes procedure with comparable accuracy, is presented for solving problems of compressible viscous subsonic flow. Results are presented and compared with experimental data for the cases of developing turbulent flow in a circular pipe; turbulent flow in a two-dimensional S-duct; and turbulent flow in a typical subsonic diffuser. Prior to each main marching step, a preliminary marching step is taken in which the integral continuity equation and an uncoupled form of the streamwise momentum equation are solved simultaneously to obtain the viscous pressure correction. During the main marching step the equations for continuity, streamwise momentum, cross-flow momentum, and energy are solved simultaneously as a coupled system using an implicit finite-difference method, with the viscous pressure correction treated as a source term. The analysis may be used for flows with both favorable and adverse pressure gradients and to predict the location of flow separation.

Towne, C. E.↗

Time-dependent Navier-Stokes computations of separated flows over airfoils

The application of both a central-difference and an upwind implicit approximate factorization Navier-Stokes algorithm to highly separated flow is described. Using the thin-layer approximation, both algorithms are employed to solve the low Reynolds number laminar flow around a circular cylinder with periodic shedding. Results agree with experimental Strouhal numbers to within about 10 percent. A grid refinement study is conducted using the upwind code for both the thin-layer and complete Navier-Stokes approximations. The circular cylinder solution is shown to be highly dependent upon both grid density outward from the body and grid extent. Finally, the thin-layer central-difference code is used to predict the time-accurate unsteady flow about two stalled airfoils. The NACA 0012 airfoil turbulent solutions are periodic with Strouhal numbers of about 0.1; a strong vortex is periodically shed from the leading edge. The LS(1)-0417 airfoil turbulent solution shows a high frequency oscillatory trailing edge separation. The 20 degree NACA 0012 airfoil solution is massively separated, but shows no evidence of periodicity.

Rumsey, C. L.↗

On recent advances and future research directions for computational fluid dynamics

This paper highlights some recent accomplishments regarding CFD numerical algorithm constructions for generation of discrete approximate solutions to classes of Reynolds-averaged Navier-Stokes equations. Following an overview of turbulent closure modeling, and development of appropriate conservation law systems, a Taylor weak-statement semi-discrete approximate solution algorithm is developed. Various forms for completion to the final linear algebra statement are cited, as are a range of candidate numerical linear algebra solution procedures. This development sequence emphasizes the key building blocks of a CFD RNS algorithm, including solution trial and test spaces, integration procedure and added numerical stability mechanisms. A range of numerical results are discussed focusing on key topics guiding future research directions.

Baker, A. J.↗

Numerical evaluation of the incomplete airy functions and their application to high frequency scattering and diffraction

The incomplete Airy integrals serve as canonical functions for the uniform ray optical solutions to several high frequency scattering and diffraction problems that involve a class of integrals characterized by two stationary points that are arbitrarily close to one another or to an integration endpoint. Integrals of such analytical properties describe transition region phenomena associated with composite shadow boundaries. An efficient and accurate method for computing the incomplete Airy functions would make the solutions to such problems useful for engineering purposes. Here, a convergent series solution form for the incomplete Airy functions is derived. Asymptotic expansions involving several terms were also developed and serve as large argument approximations. The combination of the series solution form with the asymptotic formulae provides for an efficient and accurate computation of the incomplete Airy functions. Validation of accuracy is accomplished using direct numerical integration data.

Constantinides, E. D.↗

An input adaptive, pursuit tracking model of the human opertor

Developed and evaluated is a simple model of the input adaptive behavior of the human operator (HO) in a pursuit tracking task in which the plant controlled consists of a pure gain. If it is assumed that the HO is approximately an optimal predictor using only position and velocity information, then there is a simple method of computing the values of the model parameters in terms of the autocorrelation function of the input signal. Experimental evidence indicates that the ability of the HO to use velocity information decreases with increasing signal velocity indicating that a biased estimator of the velocity weighting should be used. A suitable approximation is derived which has rapid convergence and low variance. The model thus derived is compared to actual subject transfer functions and is found to be in close agreement. In addition to tracking random processes the model can adapt to and track deterministic signals, such as sine waves, up to approximately the frequency at which human operators begin to track precognitively.

Ware, J. R.↗

Alternative approximation concepts for space frame synthesis

A structural synthesis methodology for the minimum mass design of 3-dimensionall frame-truss structures under multiple static loading conditions and subject to limits on displacements, rotations, stresses, local buckling, and element cross-sectional dimensions is presented. A variety of approximation concept options are employed to yield near optimum designs after no more than 10 structural analyses. Available options include: (A) formulation of the nonlinear mathematcal programming problem in either reciprocal section property (RSP) or cross-sectional dimension (CSD) space; (B) two alternative approximate problem structures in each design space; and (C) three distinct assumptions about element end-force variations. Fixed element, design element linking, and temporary constraint deletion features are also included. The solution of each approximate problem, in either its primal or dual form, is obtained using CONMIN, a feasible directions program. The frame-truss synthesis methodology is implemented in the COMPASS computer program and is used to solve a variety of problems. These problems were chosen so that, in addition to exercising the various approximation concepts options, the results could be compared with previously published work.

Lust, R. V.↗

Receptivity and Transition of Supersonic Boundary Layers Over Swept Wings

The receptivity, stability, and transition of three-dimensional supersonic boundary layers over (1) a swept cylinder, (2) a swept wing with a sharp leading edge, and (3) a swept wing with a blunt leading edge are numerically investigated for a free-stream Mach number of 3. These computations are compared to an earlier experimental and computational study performed by Archambaud et al.1 The steady flow fields with and without roughness elements are obtained by solving the full Navier-Stokes equations. The N-factors computed in this study at the transition onset locations reported in Ref. 1 for flow over the swept cylinder are approximately 16.5 for traveling crossflow disturbances and 9 for stationary disturbances. The N-factors for the traveling crossflow are high based on our past experiences. However, they are comparatively smaller than those reported by Archambaud et al., who found N-factor values in the range of 20 to 25 for traveling disturbances and 13 to 20 for stationary disturbances. Similarly, the N-factors computed in this study for the traveling and stationary disturbances for the flow over the sharp wing are approximately 7 and 2.5, respectively, and for the flow over the blunt wing are 6.5 and 4.8, respectively. Using the envelope method, Archambaud et al. obtained values of approximately 8.0 and 4.0 for the sharp wing case and 16.0 and 12.0 for the blunt wing case.

Balakumar, Ponnampalam↗

Long period nodal motion of sun synchronous orbits

An approximative model is formulated for assessing these perturbations that significantly affect long term modal motion of sun synchronous orbits. Computer simulations with several independent computer programs consider zonal and tesseral gravitational harmonics, third body gravitational disturbances induced by the sun and the moon, and atmospheric drag. A pendulum model consisting of evenzonal harmonics through order 4 and solar gravity dominated nodal motion approximation. This pendulum motion results from solar gravity inducing an inclination oscillation which couples into the nodal precession induced by the earth's oblateness. The pendulum model correlated well with simulations observed flight data.

Duck, K. I.↗

Observations on approximate integrations

A class of integration strategies is explored that fall in between the two extremes of symbolic integration and numerical quadrature, which are, respectively, aimed at the computer generation of answers in the form of exact expressions and numerical values. The theoretical advances in symbolic integration are discussed. Three major contexts of applications with attendant case studies, four possible types of strategies for approximate integration are examined. The feasibility and adequacy (or inadequancy) of MACSYMA for implementing these strategies are considered.

Ng, E. W.↗

Engineering flowfield method with angle-of-attack applications

An approximate inviscid flowfield method has been extended to include heat-transfer predictions using a technique to account for variable-entropy edge conditions. The engineering code computes the flowfield over hyperboloids, ellipsoids, paraboloids, and sphere cones at 0 deg angle of attack (AOA). For angle-of-attack applications, an approximation to sphere-cone streamline-spreading effects on the heat transfer along the windward and leeward rays and an empirical circumferential heating technique have been incorporated also in the method. The present engineering calculations yield good comparisons with existing pressure and heating data over sphere cones even at high incidence values with the restriction that the sonic-line location remain on the spherical cap.

Zoby, E. V.↗