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,099 records · Page 61

Detection of atmospheric pressure loading using very long baseline interferometry measurements

Loading of the Earth by the temporal redistribution of global atmospheric mass is likely to displace the positions of geodetic monuments by tens of millimeters both vertically and horizontally. Estimates of these displacements are determined by convolving National Meteorological Center (NMC) global values of atmospheric surface pressure with Farrell's elastic Green's functions. An analysis of the distances between radio telescopes determined by very long baseline interferometry (VLBI) between 1984 and 1992 reveals that in many of the cases studied there is a significant contribution to baseline length change due to atmospheric pressure loading. Our analysis covers intersite distances of between 1000 and 10,000 km and is restricted to those baselines measured more than 100 times. Accounting for the load effects (after first removing a best fit slope) reduces the weighted root-mean-square (WRMS) scatter of the baseline length residuals on 11 of the 22 baselines investigated. The slight degradation observed in the WRMS scatter on the remaining baselines is largely consistent with the expected statistical fluctuations when a small correction is applied to a data set having a much larger random noise. The results from all baselines are consistent with approximately 60% of the computed pressure contribution being present in the VLBI length determinations. Site dependent coefficients determined by fitting local pressure to the theoretical radial displacement are found to reproduce the deformation caused by the regional pressure to within 25% for most inland sites. The coefficients are less reliable at near coastal and island stations.

Vandam, T. M.↗

Microwave backscattering and emission model for grass canopies

Microwave radar and radiometer measurements of grasslands indicate a substantial reduction in sensor sensitivity to soil moisture in the presence of a thatch layer. When this layer is wet it masks changes in the underlying soil, making the canopy appear warm in the case of passive sensors (radiometer) and decreasing backscatter in the active case (scatterometer). A model for a grass canopy with thatch will be presented in this paper to explain this behavior and to compare with observations. The canopy model consists of three layers: grass, thatch, and the underlying soil. The grass blades are modeled by elongated elliptical discs and the thatch is modeled as a collection of disk shaped water droplets (i.e., the dry matter is neglected). The ground is homogeneous and flat. The distorted Born approximation is used to compute the radar cross section of this three layer canopy and the emissivity is computed from the radar cross section using the Peake formulation for the passive problem. Results are computed at L-band (1.4 GHz) and C-band (4.75 GHz) using canopy parameters (i.e., plant geometry, soil moisture, plant moisture, etc.) representative of Konza Prairie grasslands. The results are compared to C-band scatterometer measurements and L-band radiometer measurements at these grasslands.

TRT-THEORETICAL↗

Solution-Adaptive Cartesian Cell Approach for Viscous and Inviscid Flows

A Cartesian cell-based approach for adaptively refined solutions of the Euler and Navier-Stokes equations in two dimensions is presented. Grids about geometrically complicated bodies are generated automatically, by the recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, polygonal cut cells are created using modified polygon-clipping algorithms. The grid is stored in a binary tree data structure that provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite volume formulation. The convective terms are upwinded: A linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The results of a study comparing the accuracy and positivity of two classes of cell-centered, viscous gradient reconstruction procedures is briefly summarized. Adaptively refined solutions of the Navier-Stokes equations are shown using the more robust of these gradient reconstruction procedures, where the results computed by the Cartesian approach are compared to theory, experiment, and other accepted computational results for a series of low and moderate Reynolds number flows.

Coirier, William J.↗

Hybrid Analytical Technique for Nonlinear Vibration Analysis of Thin-Walled Beams

A two-step hybrid analytical technique is presented for the nonlinear vibration analysis of thin-walled beams. The first step involves the generation of various-order perturbation functions using the Linstedt-Poincare perturbation technique. The second step consists of using the perturbation functions as coordinate (or approximation) functions and then computing both the amplitudes of these functions and the nonlinear frequency of vibration via a direct variational procedure. The analytical formulation is based on a form of the geometrically nonlinear beam theory with the effects of in-plane inertia, rotatory inertia, and transverse shear deformation included. The effectiveness of the proposed technique is demonstrated by means of a numerical example of thin-walled beam with a doubly symmetric I-section. The solutions obtained using a single-spatial mode were compared with those obtained using multiple-spatial modes. The standard of comparison was taken to be the frequencies obtained by the direct integration/fast Fourier transform (FFT) technique. The nonlinear frequencies obtained by the hybrid technique were shown to converge to the corresponding ones obtained by the direct integration/fast Fourier transform (FFT) technique well beyond the range of applicability of the perturbation technique. The frequencies and total strain energy of the beam were overestimated by using a single-spatial mode.

Noor, Ahmed K.↗

Preliminary Work for Modeling the Propellers of an Aircraft as a Noise Source in an Acoustic Boundary Element Analysis

An algorithm for generating appropriate velocity boundary conditions for an acoustic boundary element analysis from the kinematics of an operating propeller is presented. It constitutes the initial phase of Integrating sophisticated rotorcraft models into a conventional boundary element analysis. Currently, the pressure field is computed by a linear approximation. An initial validation of the developed process was performed by comparing numerical results to test data for the external acoustic pressure on the surface of a tilt-rotor aircraft for one flight condition.

Vlahopoulos, Nickolas↗

DAO's Next Generation Physical-Space/Finite-Volume Data Assimilation System: Formulation and Initial Evaluation

The Physical-space/Finite-volume Data Assimilation System (fvDAS) is the next generation global atmospheric data assimilation system in development at the Data Assimilation Office at NASA's Goddard Space Flight Center. It is based on a new finite-volume general circulation model jointly developed by NASA and NCAR, and on the Physical-Space Statistical Analysis System (PSAS) developed at the DAO. In this talk we will describe the general system formulation, the adaptive quality control and general aspects of the error covariance modeling. The NASA-NCAR GCM is a completely new model which replaces the CEOs GCM used in the previous GEOS-1/2/3 Data Assimilation systems. A particular configuration of adaptive Statistical Quality Control and the Physical-space Statistical Analysis System (PSAS) are currently implemented in DAO's operational Data Assimilation System. However, the unique finite-volume formulation of the NASA-NCAR GCM, combined with the generality of the observation-space formulation of PSAS, provides for a very simple and accurate model-analysis interface. The system assimilates a variety of conventional and satellite observations. In particular, TOVS Level 1B radiances are assimilated using a 1-D variational scheme, both in clear sky and cloudy conditions. Computationally, the fvDAS runs approximately 10 times faster than the operational GEOS-Terra system. We will show that the next-generation fvDAS has much improved observation-minus-6hr forecast (O-F) statistics, as well as 5-day forecast skills. Top of the atmosphere radiation fields are in closer agreement with CERES measurements, with realistic precipitation and moisture fields. We will also show that the finite-volume formulation of the fvDAS produce assimilated fields which are more suitable for driving constituent transport models.

daSilva, A.↗

A Generalized Fluid Formulation for Turbomachinery Computations

A generalized formulation of the equations of motion of an arbitrary fluid are developed for the purpose of defining a common iterative algorithm for computational procedures. The method makes use of the equations of motion in conservation form with separate pseudo-time derivatives used for defining the numerical flux for a Riemann solver and the convergence algorithm. The partial differential equations are complemented by an thermodynamic and caloric equations of state of a complexity necessary for describing the fluid. Representative solutions with a new code based on this general equation formulation are provided for three turbomachinery problems. The first uses air as a working fluid while the second uses gaseous oxygen in a regime in which real gas effects are of little importance. These nearly perfect gas computations provide a basis for comparing with existing perfect gas code computations. The third case is for the flow of liquid oxygen through a turbine where real gas effects are significant. Vortex shedding predictions with the LOX formulations reduce the discrepancy between perfect gas computations and experiment by approximately an order of magnitude, thereby verifying the real gas formulation as well as providing an effective case where its capabilities are necessary.

Merkle, Charles L.↗

Development of a Nonequilibrium Radiative Heating Prediction Method for Coupled Flowfield Solutions

A method for predicting radiative heating and coupling effects in nonequilibrium flow-fields has been developed. The method resolves atomic lines with a minimum number of spectral points, and treats molecular radiation using the smeared band approximation. To further minimize computational time, the calculation is performed on an optimized spectrum, which is computed for each flow condition to enhance spectral resolution. Additional time savings are obtained by performing the radiation calculation on a subgrid optimally selected for accuracy. Representative results from the new method are compared to previous work to demonstrate that the speedup does not cause a loss of accuracy and is sufficient to make coupled solutions practical. The method is found to be a useful tool for studies of nonequilibrium flows.

Hartung, Lin C.↗

The Monoplane as a Lifting Vortex Surface

In Prandtl's airfoil theory the monoplane was replaced by a single lifting vortex line and yielded fairly practical results. However, the theory remained restricted to the straight wing. Yawed wings and those curved in flight direction could not be computed with this first approximation; for these the chordwise lift distribution must be taken into consideration. For the two-dimensional problem the transition from the lifting line to the lifting surface has been explained by Birnbaum. In the present report the transition to the three-dimensional problem is undertaken. The first fundamental problem involves the prediction of flow, profile, and drag for prescribed circulation distribution on the straight rectangular wing, the yawed wing for lateral boundaries parallel to the direction of flight, the swept-back wing, and the rectangular wing in slipping, with the necessary series developments for carrying through the calculations, the practical range of convergence of which does not comprise the wing tips or the break point of the swept-back wing. The second problem concerns the calculation of the circulation distribution with given profile for a slipping rectangular monoplane with flat profile and aspect ratio 6, and a rectangular wing with cambered profile and variable aspect ratio-the latter serving as check of the so-called conversion formulas of the airfoil theory.

Blenk, Hermann↗

A Reduced Dimension Static, Linearized Kalman Filter and Smoother

An approximate Kalman filter and smoother, based on approximations of the state estimation error covariance matrix, is described. Approximations include a reduction of the effective state dimension, use of a static asymptotic error limit, and a time-invariant linearization of the dynamic model for error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. Examples of use come from TOPEX/POSEIDON.

Kalman filter smoother modeling ocean simulation↗

Regression Verification Using Impact Summaries

Regression verification techniques are used to prove equivalence of syntactically similar programs. Checking equivalence of large programs, however, can be computationally expensive. Existing regression verification techniques rely on abstraction and decomposition techniques to reduce the computational effort of checking equivalence of the entire program. These techniques are sound but not complete. In this work, we propose a novel approach to improve scalability of regression verification by classifying the program behaviors generated during symbolic execution as either impacted or unimpacted. Our technique uses a combination of static analysis and symbolic execution to generate summaries of impacted program behaviors. The impact summaries are then checked for equivalence using an o-the-shelf decision procedure. We prove that our approach is both sound and complete for sequential programs, with respect to the depth bound of symbolic execution. Our evaluation on a set of sequential C artifacts shows that reducing the size of the summaries can help reduce the cost of software equivalence checking. Various reduction, abstraction, and compositional techniques have been developed to help scale software verification techniques to industrial-sized systems. Although such techniques have greatly increased the size and complexity of systems that can be checked, analysis of large software systems remains costly. Regression analysis techniques, e.g., regression testing [16], regression model checking [22], and regression verification [19], restrict the scope of the analysis by leveraging the differences between program versions. These techniques are based on the idea that if code is checked early in development, then subsequent versions can be checked against a prior (checked) version, leveraging the results of the previous analysis to reduce analysis cost of the current version. Regression verification addresses the problem of proving equivalence of closely related program versions [19]. These techniques compare two programs with a large degree of syntactic similarity to prove that portions of one program version are equivalent to the other. Regression verification can be used for guaranteeing backward compatibility, and for showing behavioral equivalence in programs with syntactic differences, e.g., when a program is refactored to improve its performance, maintainability, or readability. Existing regression verification techniques leverage similarities between program versions by using abstraction and decomposition techniques to improve scalability of the analysis [10, 12, 19]. The abstractions and decomposition in the these techniques, e.g., summaries of unchanged code [12] or semantically equivalent methods [19], compute an over-approximation of the program behaviors. The equivalence checking results of these techniques are sound but not complete-they may characterize programs as not functionally equivalent when, in fact, they are equivalent. In this work we describe a novel approach that leverages the impact of the differences between two programs for scaling regression verification. We partition program behaviors of each version into (a) behaviors impacted by the changes and (b) behaviors not impacted (unimpacted) by the changes. Only the impacted program behaviors are used during equivalence checking. We then prove that checking equivalence of the impacted program behaviors is equivalent to checking equivalence of all program behaviors for a given depth bound. In this work we use symbolic execution to generate the program behaviors and leverage control- and data-dependence information to facilitate the partitioning of program behaviors. The impacted program behaviors are termed as impact summaries. The dependence analyses that facilitate the generation of the impact summaries, we believe, could be used in conjunction with other abstraction and decomposition based approaches, [10, 12], as a complementary reduction technique. An evaluation of our regression verification technique shows that our approach is capable of leveraging similarities between program versions to reduce the size of the queries and the time required to check for logical equivalence. The main contributions of this work are: - A regression verification technique to generate impact summaries that can be checked for functional equivalence using an off-the-shelf decision procedure. - A proof that our approach is sound and complete with respect to the depth bound of symbolic execution. - An implementation of our technique using the LLVMcompiler infrastructure, the klee Symbolic Virtual Machine [4], and a variety of Satisfiability Modulo Theory (SMT) solvers, e.g., STP [7] and Z3 [6]. - An empirical evaluation on a set of C artifacts which shows that the use of impact summaries can reduce the cost of regression verification.

Backes, John↗

IKOS: Sound Static Program Analysis

IKOS (Inference Kernel for Open Static Analyzers) is a static analyzer for C/C++ based on the theory of Abstract Interpretation. It can detect or prove the absence of runtime errors (e.g, buffer overflows, integer overflows, null pointer dereferences, etc.) in the source code. IKOS uses Abstract Interpretation techniques to compute an over-approximation of all the reachable states of the program, thus it cannot miss a bug. In this talk, I will give an overview of the tool, then show how to apply it to a large software. I will present ikos-view, a web interface to examine the analysis results. I will discuss about methods to improve the analysis, such as adding code annotations, modeling library functions, and avoiding specific code patterns.

Arthaud, Maxime↗

Optimization of Depletion Modeling and Simulation for the High Flux Isotope Reactor

Monte Carlo based depletion tools used for the high-fidelity modeling and simulation of the High Flux Isotope Reactor (HFIR) come at a great computational cost; finding sufficient approximations is necessary to make the use of these tools feasible. The optimization of the neutronics and depletion model for the HFIR is based on two factors: (i) the explicit representation of the involute fuel plates with sets of polyhedra and (ii) the treatment of depletion mixtures and control element position during depletion calculations. A very fine representation (i.e., more polyhedra in the involute plate approximation) does not significantly improve simulation accuracy. The recommended representation closely represents the physical plates and ensures sufficient fidelity in regions with high flux gradients. Including the fissile targets in the central flux trap of the reactor as depletion mixtures has the greatest effect on the calculated cycle length, while localized effects (e.g., the burnup of specific isotopes or the power distribution evolution over the cycle) are more noticeable consequences of including a critical control element search or depleting burnable absorbers outside the fuel region.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Parameter distribution mapping in a gas turbine engine

A profile map of parameter values in a region is determined based on average values along linear paths through the region. In one example, a temperature map of a region of a gas turbine engine is created based on time-of-flight measurements from acoustic transceiver pairs arranged circumferentially around the region. A speed of sound for each transceiver pair is determined based on the time-of-flight measurements. An average temperature along each path is estimated from the time of flight, and the profile map is computed from the average temperatures. The profile map may be computed using a polynomial approximation technique, a grid optimization technique or a basis function technique.

42 ENGINEERING↗

Impact of Spatial Discretization on Reactivity Biases in Depleted TRIGA Fuel

Ongoing work at Kansas State University (KSU) and the University of Wisconsin involves the determination of fuel compositions in TRIGA reactors. Modeling of depletion requires a number of approximations to be made, which often includes selection of an appropriate spatial discretization for tracking material evolution in time. Incorporation of temperature feedback also requires such spatial zoning. Any approximation made leads to computational biases, i.e., differences between computed and accurate values of some parameter of interest. The objective of this study is to understand the potential magnitude of biases introduced by the selected spatial discretization. Specifically, biases in the reactivity of a square-pitch unit cell representative of a typical TRIGA fuel element were studied as a function of burnup for different radial and axial discretizations. Temperature effects were also studied by introducing a fixed, radial temperature profile. The results indicate a very fine radial discretization (> 10 radial regions) may be needed to eliminate biases (of statistical significance) from temperature feedback effects, while an axial discretization of {approx_equal}7 regions may be sufficient to eliminate biases due to evolution of the fuel composition. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Optimization of thermal protection systems for the space shuttle vehicle. Volume 1: Final report

A study performed to continue development of computational techniques for the Space Shuttle Thermal Protection System is reported. The resulting computer code was used to perform some additional optimization studies on several TPS configurations. The program was developed in Fortran 4 for the CDC 6400, and it was converted to Fortran 5 to be used for the Univac 1108. The computational methodology is developed in modular fashion to facilitate changes and updating of the techniques and to allow overlaying the computer code to fit into approximately 131,000 octal words of core storage. The program logic involves subroutines which handle input and output of information between computer and user, thermodynamic stress, dynamic, and weight/estimate analyses of a variety of panel configurations. These include metallic, ablative, RSI (with and without an underlying phase change material), and a thermodynamic analysis only of carbon-carbon systems applied to the leading edge and flat cover panels. Two different thermodynamic analyses are used. The first is a two-dimensional, explicit precedure with variable time steps which is used to describe the behavior of metallic and carbon-carbon leading edges. The second is a one-dimensional implicity technique used to predict temperature in the charring ablator and the noncharring RSI. The latter analysis is performed simply by suppressing the chemical reactions and pyrolysis of the TPS material.

Source record↗

Survey of the supporting research and technology for the thermal protection of the Galileo Probe

The Galileo Probe, which is scheduled to be launched in 1985 and to enter the hydrogen-helium atmosphere of Jupiter up to 1,475 days later, presents thermal protection problems that are far more difficult than those experienced in previous planetary entry missions. The high entry speed of the Probe will cause forebody heating rates orders of magnitude greater than those encountered in the Apollo and Pioneer Venus missions, severe afterbody heating from base-flow radiation, and thermochemical ablation rates for carbon phenolic that rival the free-stream mass flux. This paper presents a comprehensive survey of the experimental work and computational research that provide technological support for the Probe's heat-shield design effort. The survey includes atmospheric modeling; both approximate and first-principle computations of flow fields and heat-shield material response; base heating; turbulence modelling; new computational techniques; experimental heating and materials studies; code validation efforts; and a set of 'consensus' first-principle flow-field solutions through the entry maneuver, with predictions of the corresponding thermal protection requirements.

Howe, J. T.↗

Bias-Variance Trade-Off in Physics-Informed Neural Networks with Randomized Smoothing for High-Dimensional PDEs

Physics-Informed Neural Networks (PINNs) have triggered a paradigm shift in scientific computing, leveraging mesh-free properties and robust approximation capabilities. While proving effective for low-dimensional partial differential equations (PDEs), the computational cost of PINNs remains a hurdle in high-dimensional scenarios. This is particularly pronounced when computing high-order and high-dimensional derivatives in the physics-informed loss. Randomized Smoothing PINN (RS-PINN) introduces Gaussian noise for stochastic smoothing of the original neural net model, enabling the use of Monte Carlo methods for derivative approximation, which eliminates the need for costly automatic differentiation. Despite its computational efficiency, especially in the approximation of high-dimensional derivatives, RS-PINN introduces biases in both loss and gradients, negatively impacting convergence, especially when coupled with stochastic gradient descent (SGD) algorithms. We present a comprehensive analysis of biases in RS-PINN, attributing them to the nonlinearity of the Mean Squared Error (MSE) loss as well as the intrinsic nonlinearity of the PDE itself. We propose tailored bias correction techniques, delineating their application based on the order of PDE nonlinearity. The derivation of an unbiased RS-PINN allows for a detailed examination of its advantages and disadvantages compared to the biased version. Specifically, the biased version has a lower variance and runs faster than the unbiased version, but it is less accurate due to the bias. To optimize the bias-variance trade-off, we combine the two approaches in a hybrid method that balances the rapid convergence of the biased version with the high accuracy of the unbiased version. In addition to methodological contributions, we present an enhanced implementation of RS-PINN. Extensive experiments on diverse high-dimensional PDEs, including Fokker-Planck, Hamilton-Jacobi-Bellman (HJB), viscous Burgers’, Allen-Cahn, and Sine-Gordon equations, illustrate the bias-variance trade-off and highlight the effectiveness of the hybrid RS-PINN. Empirical guidelines are provided for selecting biased, unbiased, or hybrid versions, depending on the dimensionality and nonlinearity of the specific PDE problem.

97 MATHEMATICS AND COMPUTING↗