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 1,117 records · Page 62

Controls on Martian Hydrothermal Systems: Application to Valley Network and Magnetic Anomaly Formation

Models of hydrothermal groundwater circulation can quantify limits to the role of hydrothermal activity in Martian crustal processes. We present here the results of numerical simulations of convection in a porous medium due to the presence of a hot intruded magma chamber. The parameter space includes magma chamber depth, volume, aspect ratio, and host rock permeability and porosity. A primary goal of the models is the computation of surface discharge. Discharge increases approximately linearly with chamber volume, decreases weakly with depth (at low geothermal gradients), and is maximized for equant-shaped chambers. Discharge increases linearly with permeability until limited by the energy available from the intrusion. Changes in the average porosity are balanced by changes in flow velocity and therefore have little effect. Water/rock ratios of approximately 0.1, obtained by other workers from models based on the mineralogy of the Shergotty meteorite, imply minimum permeabilities of 10(exp -16) sq m2 during hydrothermal alteration. If substantial vapor volumes are required for soil alteration, the permeability must exceed 10(exp -15) sq m. The principal application of our model is to test the viability of hydrothermal circulation as the primary process responsible for the broad spatial correlation of Martian valley networks with magnetic anomalies. For host rock permeabilities as low as 10(exp -17) sq m and intrusion volumes as low as 50 cu km, the total discharge due to intrusions building that part of the southern highlands crust associated with magnetic anomalies spans a comparable range as the inferred discharge from the overlying valley networks.

Harrison, Keith P.↗

Insulin Delivery System

When Programmable Implantable Medication System (PIMS) is implanted in human body, it delivers precise programmed amounts of insulin over long periods of time. Mini-Med Technologies has been refining the Technologies since initial development at APL. The size of a hockey puck, and encased in titanium shell, PIMS holds about 2 1/2 teaspoons of insulin at a programmed basal rate. If a change in measured blood sugar level dictates a different dose, the patient can vary the amount of insulin delivered by holding a small radio transceiver over the implanted system and dialing in a specific program held in the PIMS computer memory. Insulin refills are accomplished approximately 4 times a year by hypodermic needle.

Source record↗

Combining Particle Filters and Consistency-Based Approaches for Monitoring and Diagnosis of Stochastic Hybrid Systems

Fault detection and isolation are critical tasks to ensure correct operation of systems. When we consider stochastic hybrid systems, diagnosis algorithms need to track both the discrete mode and the continuous state of the system in the presence of noise. Deterministic techniques like Livingstone cannot deal with the stochasticity in the system and models. Conversely Bayesian belief update techniques such as particle filters may require many computational resources to get a good approximation of the true belief state. In this paper we propose a fault detection and isolation architecture for stochastic hybrid systems that combines look-ahead Rao-Blackwellized Particle Filters (RBPF) with the Livingstone 3 (L3) diagnosis engine. In this approach RBPF is used to track the nominal behavior, a novel n-step prediction scheme is used for fault detection and L3 is used to generate a set of candidates that are consistent with the discrepant observations which then continue to be tracked by the RBPF scheme.

Narasimhan, Sriram↗

An Engineering Aerodynamic Heating Method for Hypersonic Flow

A capability to calculate surface heating rates has been incorporated in an approximate three-dimensional inviscid technique. Surface streamlines are calculated from the inviscid solution, and the axisymmetric analog is then used along with a set of approximate convective-heating equations to compute the surface heat transfer. The method is applied to blunted axisymmetric and three-dimensional ellipsoidal cones at angle of attack for the laminar flow of a perfect gas. The method is also applicable to turbulent and equilibrium-air conditions. The present technique predicts surface heating rates that compare favorably with experimental (ground-test and flight) data and numerical solutions of the Navier-Stokes (NS) and viscous shock-layer (VSL) equations. The new technique represents a significant improvement over current engineering aerothermal methods with only a modest increase in computational effort.

Riley, Christopher J.↗

Adaptive Metropolis Sampling with Product Distributions

The Metropolis-Hastings (MH) algorithm is a way to sample a provided target distribution pi(z). It works by repeatedly sampling a separate proposal distribution T(x,x') to generate a random walk {x(t)}. We consider a modification of the MH algorithm in which T is dynamically updated during the walk. The update at time t uses the {x(t' less than t)} to estimate the product distribution that has the least Kullback-Leibler distance to pi. That estimate is the information-theoretically optimal mean-field approximation to pi. We demonstrate through computer experiments that our algorithm produces samples that are superior to those of the conventional MH algorithm.

Wolpert, David H.↗

HEC Applications on Columbia Project

NASA's Columbia system consists of a cluster of twenty 512 processor SGI Altix systems. Each of these systems is 3 TFLOP/s in peak performance - approximately the same as the entire compute capability at NAS just one year ago. Each 512p system is a single system image machine with one Linunx O5, one high performance file system, and one globally shared memory. The NAS Terascale Applications Group (TAG) is chartered to assist in scaling NASA's mission critical codes to at least 512p in order to significantly improve emergency response during flight operations, as well as provide significant improvements in the codes. and rate of scientific discovery across the scientifc disciplines within NASA's Missions. Recent accomplishments are 4x improvements to codes in the ocean modeling community, 10x performance improvements in a number of computational fluid dynamics codes used in aero-vehicle design, and 5x improvements in a number of space science codes dealing in extreme physics. The TAG group will continue its scaling work to 2048p and beyond (10240 cpus) as the Columbia system becomes fully operational and the upgrades to the SGI NUMAlink memory fabric are in place. The NUMlink uprades dramatically improve system scalability for a single application. These upgrades will allow a number of codes to execute faster at higher fidelity than ever before on any other system, thus increasing the rate of scientific discovery even further

Taft, Jim↗

Methodology for Flight Relevant Arc-Jet Testing of Flexible Thermal Protection Systems

A methodology to correlate flight aeroheating environments to the arc-jet environment is presented. For a desired hot-wall flight heating rate, the methodology provides the arcjet bulk enthalpy for the corresponding cold-wall heating rate. A series of analyses were conducted to examine the effects of the test sample model holder geometry to the overall performance of the test sample. The analyses were compared with arc-jet test samples and challenges and issues are presented. The transient flight environment was calculated for the Hypersonic Inflatable Aerodynamic Decelerator (HIAD) Earth Atmospheric Reentry Test (HEART) vehicle, which is a planned demonstration vehicle using a large inflatable, flexible thermal protection system to reenter the Earth's atmosphere from the International Space Station. A series of correlations were developed to define the relevant arc-jet test environment to properly approximate the HEART flight environment. The computed arcjet environments were compared with the measured arc-jet values to define the uncertainty of the correlated environment. The results show that for a given flight surface heat flux and a fully-catalytic TPS, the flight relevant arc-jet heat flux increases with the arc-jet bulk enthalpy while for a non-catalytic TPS the arc-jet heat flux decreases with the bulk enthalpy.

Mazaheri, Alireza↗

Spacecube V2.0 Micro Single Board Computer

A single board computer system radiation hardened for space flight includes a printed circuit board having a top side and bottom side; a reconfigurable field programmable gate array (FPGA) processor device disposed on the top side; a connector disposed on the top side; a plurality of peripheral components mounted on the bottom side; and wherein a size of the single board computer system is not greater than approximately 7 cm.times.7 cm.

Petrick, David J.↗

Using Adaptive Mesh Refinement to Study Grid Resolution Effects for Shock/Boundary-Layer Interactions

Adaptive Mesh Refinement (AMR) promises a much more computationally efficient meansto obtain a discrete approximation to a continuous boundary value problem of a specifiedaccuracy than classic isotropic grid refinement. The AMR capability of OVERFLOW is utilizedto provide estimates of the exact analytical solutions to problems of interest to turbulencemodeling. Predictions of surface pressure and skin friction, essentially the state of stress at thesurface, shows little difference with grids believed to be "grid resolved." Velocity profiles, on theother hand, show marked differences in flows with shocks. The AMR method, as implementedin OVERFLOW2.2k, appears to provide the ability to produce arbitrarily accurate solutionsat a predictable cost much smaller than classic uniform mesh refinement.

Adaptive Mesh Refinement↗

Progress on Transonic Flutter and Shock Buffet Computationsin Support of the Third Aeroelastic Prediction Workshop

This paper reports on the progress of the NASA Langley team contributions to the third Aeroelastic Pre-diction Workshop’s (AePW-3) High Angle Working Group (HAWG). The primary objectives of HAWG is to predict the fluter dynamic pressure of the NASA Benchmark Supercritical Wing (BSCW) configuration at Mach 0.8 and 5◦angle of attack. The secondary objective is to determine if a shock-buffet onset is present at or near that flow condition. The computational results are obtained using FUN3D, an unstructured grid Reynolds-averaged Navier-Stokes solver developed at the NASA Langley Research Center. The preliminary analysis results show a computationally-obtained flutter dynamic pressure of approximately 120 psf. Initial results describing unforced BSCW unsteady flow environment at flutter condition are also presented.

Pawel Chwalowski↗

Homography generation for image registration in inlier-poor domains

A method for efficient image registration between two images in the presence of inlier-poor domains includes receiving a set of candidate correspondences between the two images. An approximate homography between the two images is generated based upon a first correspondence in the correspondences. The set of candidate correspondences is filtered to identify inlier correspondences based upon the approximate homography. A candidate homography is computed based upon the inlier correspondences. The candidate homography can be selected as a final homography between the two images based upon a support of the candidate homography against the set of candidate correspondences. An image registration is performed between the two images based upon the candidate homography being selected as the final homography.

Gonzales, Antonio↗

A Bayesian Learning Approach to Wireless Outdoor Heatmap Construction using Deep Gaussian Process

We present a novel Bayesian learning approach to outdoor radio heatmap construction utilizing deep Gaussian process (GP). The proposed approach employs a two-layer hierarchy which consists of two cascaded Gaussian processes that are capable of modeling more complex input-output relations than standard single-layer Gaussian processes. Since deriving the exact model likelihood is challenging, a lower bound is optimized instead so that gradient descent-based methods can be performed to find out the optimal model parameters. Typically, inducing points are used in GPs to facilitate low-rank approximation of covariance (kernel) matrices for computation speedup. However, the inaccuracy induced by inducing points can accumulate when stacking multiple layers of GP which may hinder the performance of deep GP. Moreover, since inducing points need to be learned, having them at all layers of deep GP also incurs computational burden. To overcome the above challenges, in contrast to the canonical deep GP model, we use a modified architecture where a full standard GP resides in the first layer and inducing points are only introduced for the second layer. This modified architecture strikes a balance between model accuracy and training complexity. In the proposed model, the noise parameter of the first GP layer is also eliminated to improve the training efficiency as the noise parameter at the output of the second layer suffices to model the uncertainty in the output. The proposed approach is evaluated on real-world datasets, in the form of location-Received Signal Strength (RSS) pairs, collected from the Platform for Open Wireless Data-driven Experimental Research (POWDER) located at the campus of the University of Utah. Experiment results show that the proposed approach can achieve smaller prediction errors on various training and testing data configurations than DNN-based and GP-based methods.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Studies of electron collisions with polyatomic molecules using distributed-memory parallel computers

Elastic electron scattering cross sections from 5-30 eV are reported for the molecules C2H4, C2H6, C3H8, Si2H6, and GeH4, obtained using an implementation of the Schwinger multichannel method for distributed-memory parallel computer architectures. These results, obtained within the static-exchange approximation, are in generally good agreement with the available experimental data. These calculations demonstrate the potential of highly parallel computation in the study of collisions between low-energy electrons and polyatomic gases. The computational methodology discussed is also directly applicable to the calculation of elastic cross sections at higher levels of approximation (target polarization) and of electronic excitation cross sections.

Winstead, C.↗

Fast All-sky Radiation Model for Solar applications (FARMS) [SWR-16-18]

The Fast All-sky Radiation Model for Solar applications (FARMS) is used to compute cloudy irradiance. Radiative transfer (RT) models simulating broadband solar radiation have been widely used by atmospheric scientists to model solar resources for various energy applications such as operational forecasting. Due to the complexity of solving the RT equation, the computation under cloudy conditions can be extremely time consuming though many approximations (e.g. two-stream approach and delta-M truncation scheme) have been utilized. Thus, a more efficient RT model is crucial for model developers as a new option for approximating solar radiation at the land surface with minimal loss of accuracy. We have developed a fast all-sky radiation model for solar applications (FARMS) using the simplified clear-sky RT model, REST2, and simulated cloud transmittances and reflectances from the Rapid Radiation Transfer Model (RRTM) with a sixteen-stream Discrete Ordinates Radiative Transfer (DISORT). Simulated lookup tables (LUTs) of cloud transmittances and reflectances were created by varying cloud optical thicknesses, cloud particle sizes, and solar zenith angles. Equations with optimized parameters were fitted to the cloud transmittances and reflectances to develop the model. Using this model the all-sky solar irradiance at the land surface can be computed rapidly by combining REST2 with the cloud transmittances and reflectances. This new RT model is more than 1000 times faster than those currently utilized in solar resource assessment and forecasting since it does not explicitly solve the RT equation for each individual cloud condition. Our results indicate the accuracy of the fast radiative transfer model is comparable to or better than two-stream approximation in term of computing cloud transmittance and solar radiation.

Xie, Yu↗

Fast All-sky Radiation Model for Solar applications (FARMS) 2023 [SWR-16-18]

The Fast All-sky Radiation Model for Solar applications (FARMS) is used to compute cloudy irradiance. Radiative transfer (RT) models simulating broadband solar radiation have been widely used by atmospheric scientists to model solar resources for various energy applications such as operational forecasting. Due to the complexity of solving the RT equation, the computation under cloudy conditions can be extremely time consuming though many approximations (e.g. two-stream approach and delta-M truncation scheme) have been utilized. Thus, a more efficient RT model is crucial for model developers as a new option for approximating solar radiation at the land surface with minimal loss of accuracy. We have developed a fast all-sky radiation model for solar applications (FARMS) using the simplified clear-sky RT model, REST2, and simulated cloud transmittances and reflectances from the Rapid Radiation Transfer Model (RRTM) with a sixteen-stream Discrete Ordinates Radiative Transfer (DISORT). Simulated lookup tables (LUTs) of cloud transmittances and reflectances were created by varying cloud optical thicknesses, cloud particle sizes, and solar zenith angles. Equations with optimized parameters were fitted to the cloud transmittances and reflectances to develop the model. Using this model the all-sky solar irradiance at the land surface can be computed rapidly by combining REST2 with the cloud transmittances and reflectances. This new RT model is more than 1000 times faster than those currently utilized in solar resource assessment and forecasting since it does not explicitly solve the RT equation for each individual cloud condition. Our results indicate the accuracy of the fast radiative transfer model is comparable to or better than two-stream approximation in term of computing cloud transmittance and solar radiation.

Xie, Yu↗

Efficient method for approximating nonlinear dynamics: applications to uncertainty propagation and estimation

High-order Taylor series expansions can be used to model nonlinear dynamics at the cost of integrating a large set of variational equations to obtain high-order state-transition tensors (STTs). This paper presents an innovative technique for approximating the high-order STTs that reduces significantly the computational cost by retaining only the dominant secular terms. We propagate the low-order partial derivatives of Kepler’s equation, which only requires the integration of six additional equations to extend an n-th order approximation to order (n + 1). The approximation stems from the Lindstedt-Poincare procedure and exploits the stability properties of orbital motion. Since the method makes no dynamical assumptions, it can accommodate any source of orbital perturbations. We show how the approximation of the second-order STT significantly increases the accuracy of the linear method for uncertainty propagation with only a small computational overhead. Finally, we derive a high-order approximate extended Kalman filter that implements the proposed approximation of the STT and improves the performance of linear filters. Examples of application with different perturbation sources include the heliocentric orbit of an asteroid, an orbiter around Europa, and an Earth-orbiting satellite.

Park, Ryan S↗

Multi Model Monte Carlo with Python (MXMCPy)

Multi Model Monte Carlo with Python (\mxmc {}) is a software package developed as a general capability for computing the statistics of outputs from an expensive, high-fidelity model by leveraging faster, low-fidelity models for speedup. Motivated by uncertainty propagation problems where classical Monte Carlo (MC) simulation is computationally intractable, various multi-model MC approaches have recently emerged that yield unbiased estimators with significantly reduced variance relative to MC for the same cost. These existing methods include multi-level Monte Carlo (MLMC), multi-fidelity Monte Carlo (MFMC), and approximate control variates (ACV). Given a fixed computational budget and a collection of models with varying cost/accuracy, each method seeks a sample allocation strategy across the models that results in an estimator with optimal variance reduction. \mxmc {} is a versatile tool that enables convenient access to many existing multi-model MC approaches within one modular and extensible package. With \mxmc {}, users can easily compare existing methods to determine the best choice for their particular problem, while developers have a basis for implementing and sharing new variance reduction approaches. This report introduces the \mxmc {} software, providing a summary of the problem-solving workflow for users as well as a brief overview of the code layout for developers.

Geoffrey F Bomarito↗

Structure-preserving neural networks for the regularized entropy-based closure of a linear, kinetic, radiative transport equation

The main challenge of large-scale numerical simulation of radiation transport is the high memory and computation time requirements of discretization methods for kinetic equations. In this work, we derive and investigate a neural network-based approximation to the entropy-based closure method to accurately compute the solution of the multi-dimensional moment system with a low memory footprint and competitive computational time. We extend methods developed for the standard entropy-based closure to the regularized entropy-based closures. The main idea is to interpret structure-preserving neural network approximations of the regularized entropy-based closure as a two-stage approximation to the original entropy-based closure. We conduct a numerical analysis of this approximation and investigate optimal parameter choices. Our numerical experiments demonstrate that the method has a much lower memory footprint than traditional methods with competitive computation times and simulation accuracy. The code and all trained networks are provided on GitHub.

entropy closure↗