Search NASA⌕ Search

SEARCH · Search NASA

Results for “algorithm timings”

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 289 records · Page 16

Towards developing robust algorithms for solving partial differential equations on MIMD machines

Methods for efficient computation of numerical algorithms on a wide variety of MIMD machines are proposed. These techniques reorganize the data dependency patterns to improve the processor utilization. The model problem finds the time-accurate solution to a parabolic partial differential equation discretized in space and implicitly marched forward in time. The algorithms are extensions of Jacobi and SOR. The extensions consist of iterating over a window of several timesteps, allowing efficient overlap of computation with communication. The methods increase the degree to which work can be performed while data are communicated between processors. The effect of the window size and of domain partitioning on the system performance is examined both by implementing the algorithm on a simulated multiprocessor system.

Saltz, Joel H.↗

A streamwise upwind algorithm for the Euler and Navier-Stokes equations applied to transonic flows

A new algorithm was developed for the Euler and Navier-Stokes equations that uses upwind differencing based on the streamwise direction. This algorithm is time accurate and can be used in codes for calculating unsteady transonic flows over wings. Such codes can be used for the flutter analysis of wings. In this algorithm, the coordinate system is locally rotated to align with the streamwise direction. For differencing the convective terms in the streamwise direction, a new form of flux splitting is employed, in which the biasing depends on the local Mach number. In the plane perpendicular to the stream direction, the new flux splitting uses the condition of no flow in that local plane. By using a locally rotated coordinate system, the convective flux vector biasing depends on the total Mach number. Hence, the switching of the flux vector biasing occurs across shock waves and the proper domain of dependence is used in supersonic regions. For comparison, many other upwind methods switch differencing based on Mach number of shock wave in multidimensional flows. The formulas for the convective flux vector differencing do not contain any user specified parameters. So, the amount of numerical dissipation is automatically determined.

Goorjian, P. M.↗

Concurrent and vectorized mixed time, explicit nonlinear structural dynamics algorithms

A nonlinear structural dynamics program with an element library that exploits parallel processing is described. The aim is to exploit scheduling-allocation so that parallel processing and vectorization can effectively be treated in a general purpose program with explicit time integration and different time steps in different parts of the mesh. The program uses an element group scheme, which, as a by-product, also provides an automatic scheme for assigning different time steps to different parts of the mesh. The program has been tested on the Alliant FX/8; it shows a fivefold improvement in speed over compiler optimization.

Belytschko, Ted↗

Automatic voice recognition using traditional and artificial neural network approaches

The main objective of this research is to develop an algorithm for isolated-word recognition. This research is focused on digital signal analysis rather than linguistic analysis of speech. Features extraction is carried out by applying a Linear Predictive Coding (LPC) algorithm with order of 10. Continuous-word and speaker independent recognition will be considered in future study after accomplishing this isolated word research. To examine the similarity between the reference and the training sets, two approaches are explored. The first is implementing traditional pattern recognition techniques where a dynamic time warping algorithm is applied to align the two sets and calculate the probability of matching by measuring the Euclidean distance between the two sets. The second is implementing a backpropagation artificial neural net model with three layers as the pattern classifier. The adaptation rule implemented in this network is the generalized least mean square (LMS) rule. The first approach has been accomplished. A vocabulary of 50 words was selected and tested. The accuracy of the algorithm was found to be around 85 percent. The second approach is in progress at the present time.

Botros, Nazeih M.↗

A New Cell-Centered Implicit Numerical Scheme for Ions in the 2-D Axisymmetric Code Hall2de

We present a new algorithm in the Hall2De code to simulate the ion hydrodynamics in the acceleration channel and near plume regions of Hall-effect thrusters. This implementation constitutes an upgrade of the capabilities built in the Hall2De code. The equations of mass conservation and momentum for unmagnetized ions are solved using a conservative, finite-volume, cell-centered scheme on a magnetic-field-aligned grid. Major computational savings are achieved by making use of an implicit predictor/multi-corrector algorithm for time evolution. Inaccuracies in the prediction of the motion of low-energy ions in the near plume in hydrodynamics approaches are addressed by implementing a multi-fluid algorithm that tracks ions of different energies separately. A wide range of comparisons with measurements are performed to validate the new ion algorithms. Several numerical experiments with the location and value of the anomalous collision frequency are also presented. Differences in the plasma properties in the near-plume between the single fluid and multi-fluid approaches are discussed. We complete our validation by comparing predicted erosion rates at the channel walls of the thruster with measurements. Erosion rates predicted by the plasma properties obtained from simulations replicate accurately measured rates of erosion within the uncertainty range of the sputtering models employed.

electric propulsion↗

Attitude Determination by Using Horizon and Sun Sensors

The Pointing and Alignment Workstation (PAWS) developed by Teledyne Brown Engineering (TBE) has successfully supported the first and second Atmospheric Laboratory for Applications and Science (ATLAS 1, 2) spacelab missions for NASA. The primary PAWS objective was to provide realtime pointing information to instruments whose line of-sight is dependent on Shuttle attitude and to study/quantify the causes and effects of Shuttle and payload pointing errors. In addition to Shuttle IMU attitude information, PAWS used atmospheric science sensors data to determine the spacecraft attitude. PAWS successfully achieved these goals by acquiring and processing data during the ATLAS 1, 2 mission. This paper presents the attitude determination algorithm real time processing, and results of post mission analysis. The findings of this study include the quality of the horizon sensor and IMU measurements as well as accuracy of attitude processor algorithm.

Huang, Allen K. H.↗

Exponential Convergence of the Tracking Error in Adaptive Systems Without Persistent Excitation

Persistent excitation conditions which ensure parameter convergence in adaptive algorithms have been studied by many researchers. Here, conditions are investigated for exponential convergence of the tracking error in feedforward adaptive systems without persistent excitation. Particular attention is paid to the continuous-time LMS algorithm in the overparametrized case. Results are presented.

control↗

A generalized memory test algorithm

A general algorithm for testing digital computer memory is presented. The test checks that (1) every bit can be cleared and set in each memory work, and (2) bits are not erroneously cleared and/or set elsewhere in memory at the same time. The algorithm can be applied to any size memory block and any size memory word. It is concise and efficient, requiring the very few cycles through memory. For example, a test of 16-bit-word-size memory requries only 384 cycles through memory. Approximately 15 seconds were required to test a 32K block of such memory, using a microcomputer having a cycle time of 133 nanoseconds.

Milner, E. J.↗

Earth radiation budget measurements - Pre-ERBE, ERBE, and CERES

The paper reviews experiments dealing with the earth's radiation budget for the last two decades. Attention is focused on the major areas of concern in the satellite measurements: instrument calibration, angular dependence of the radiation leaving the atmosphere, and space and time dependence of the flux fields. It is pointed out that the current radiation-budget measurements come from the Earth Radiation Budget Experiment (ERBE) with improved calibration, estimates of the instantaneous flux, and algorithms for time interpolation. A new set of measurements will be produced from an investigation of the Clouds and the Earth's Radiant Energy System (CERES), a part of the Earth Observing System (EOS). CERES is expected to produce the radiation field from the top of the atmosphere to the earth's surface.

Barkstrom, Bruce R.↗

Phase Sensor for Aligning a Segmented Telescope Mirror

A phase sensor has been developed for use in aligning a segmented telescope mirror to within a fraction of a wavelength in piston. (As used here, piston signifies displacement of a mirror segment along the optical axis of the telescope.) Such precise alignment is necessary in order to realize the full benefit of the large aperture achievable through segmentation. This phase sensor is achromatic. It is based on two-wavelength shearing interferometry, and can be modified to utilize an extended or broad-band (e.g., white) light source. The sensor optics include a ruled diffraction grating and an imaging lens. The sensor can measure the piston shift between segments as well as aberrations of the segments. It can measure the surface error of an individual segment, making it possible to compensate for the error with optimal amount(s) of piston and/or tilt. The precise capture range of the sensor depends partly on the telescope design; the largest relative piston shifts measurable by use of this sensor are of the order of 100 m. The accuracy of the sensor also depends partly on the telescope design; in general, the accuracy is sufficient to enable alignment to within approximately half a wavelength. The interferometric image is digitized and processed by a simple algorithm in real time, and the output of the algorithm can be used to maintain alignment in real time, even in the presence of atmospheric turbulence. The sensor is robust. Through calibration, it can be made insensitive to (and, hence, tolerant of) misalignments and aberrations of its own optics, most aberrations of the telescope as a whole (in contradistinction to aberrations of individual segments), and most aberrations introduced by atmospheric turbulence

Stahl, Philip↗

Accelerating Time-Varying Hardware Volume Rendering Using TSP Trees and Color-Based Error Metrics

This paper describes a new hardware volume rendering algorithm for time-varying data. The algorithm uses the Time-Space Partitioning (TSP) tree data structure to identify regions within the data that have spatial or temporal coherence. By using this coherence, the rendering algorithm can improve performance when the volume data is larger than the texture memory capacity by decreasing the amount of textures required. This coherence can also allow improved speed by appropriately rendering flat-shaded polygons instead of textured polygons, and by not rendering transparent regions. To reduce the polygonization overhead caused by the use of the hierarchical data structure, we introduce an optimization method using polygon templates. The paper also introduces new color-based error metrics, which more accurately identify coherent regions compared to the earlier scalar-based metrics. By showing experimental results from runs using different data sets and error metrics, we demonstrate that the new methods give substantial improvements in volume rendering performance.

Ellsworth, David↗

Evaluating Retrieval Algorithm Climate Stability: Estimating 3D Optical Thickness Bias Distributions by Cloud Type

Detecting climate trends on large spatiotemporal scales requires accurate, stable measurements and stable retrieval algorithms. We strive to estimate how time-variant retrieval algorithm biases may impact trend detection. Here we focus on the 3D cloud optical thickness (τc) bias, which is among the largest in passive cloud retrieval algorithms. If this bias is time dependent, a possibility with potential decadal changes in cloud morphology, it may obscure genuine trends in τc. Although previous studies have evaluated the cloud- and sun-view geometry-dependent 3D τc bias on small spatial scales, before our current study none have evaluated the stability of this well-known bias on climate-relevant large spatiotemporal scales. These studies must estimate large scale distributions of the 3D τc bias by cloud type and estimate how cloud type amount may change between two climate states. We employ a novel approach to estimate large scale distributions of 3D τc using a proxy of the bias that quantifies the departure of clouds from satisfying the 1D radiative transfer assumption used in passive τc retrievals. This existing globally-distributed proxy is an angular consistency metric that was developed using fused Moderate-Resolution Imaging Spectroradiometer (MODIS) and Multi-angle Imaging Spectroradiometer (MISR) measurements. Calculating the 3D τc bias and the proxy, for known cloud fields enables us to establish statistical relationships between these two quantities, which can be used to calculate large-scale distributions of the 3D τc bias. This approach limits the number of 3D radiative transfer simulations required to only those needed to estimate a statistical relationship between the 3D τc bias for known cloud fields and a proxy of the bias. It is likely that future studies will be needed to evaluate retrieval algorithm bias stability for other geophysical variables as the community develops climate data records from satellite observations and their retrievals. This must be done in addition to monitoring and correcting measurement errors and uncertainties and understanding their impact on retrieved essential climate variables.

Yolanda Shea↗

Parallel Simulation of Unsteady Turbulent Flames

Time-accurate simulation of turbulent flames in high Reynolds number flows is a challenging task since both fluid dynamics and combustion must be modeled accurately. To numerically simulate this phenomenon, very large computer resources (both time and memory) are required. Although current vector supercomputers are capable of providing adequate resources for simulations of this nature, the high cost and their limited availability, makes practical use of such machines less than satisfactory. At the same time, the explicit time integration algorithms used in unsteady flow simulations often possess a very high degree of parallelism, making them very amenable to efficient implementation on large-scale parallel computers. Under these circumstances, distributed memory parallel computers offer an excellent near-term solution for greatly increased computational speed and memory, at a cost that may render the unsteady simulations of the type discussed above more feasible and affordable.This paper discusses the study of unsteady turbulent flames using a simulation algorithm that is capable of retaining high parallel efficiency on distributed memory parallel architectures. Numerical studies are carried out using large-eddy simulation (LES). In LES, the scales larger than the grid are computed using a time- and space-accurate scheme, while the unresolved small scales are modeled using eddy viscosity based subgrid models. This is acceptable for the moment/energy closure since the small scales primarily provide a dissipative mechanism for the energy transferred from the large scales. However, for combustion to occur, the species must first undergo mixing at the small scales and then come into molecular contact. Therefore, global models cannot be used. Recently, a new model for turbulent combustion was developed, in which the combustion is modeled, within the subgrid (small-scales) using a methodology that simulates the mixing and the molecular transport and the chemical kinetics within each LES grid cell. Finite-rate kinetics can be included without any closure and this approach actually provides a means to predict the turbulent rates and the turbulent flame speed. The subgrid combustion model requires resolution of the local time scales associated with small-scale mixing, molecular diffusion and chemical kinetics and, therefore, within each grid cell, a significant amount of computations must be carried out before the large-scale (LES resolved) effects are incorporated. Therefore, this approach is uniquely suited for parallel processing and has been implemented on various systems such as: Intel Paragon, IBM SP-2, Cray T3D and SGI Power Challenge (PC) using the system independent Message Passing Interface (MPI) compiler. In this paper, timing data on these machines is reported along with some characteristic results.

Menon, Suresh↗

Advanced MHD Algorithm for Solar and Space Science: lst Year Semi Annual Progress Report

We report progress for the development of MH4D for the first and second quarters of FY2004, December 29, 2002 - June 6, 2003. The present version of MH4D can now solve the full viscous and resistive MHD equations using either an explicit or a semi-implicit time advancement algorithm. In this report we describe progress in the following areas. During the two last quarters we have presented poster at the EGS-AGU-EUG Joint Assembly in Nice, France, April 6-11, 2003, and a poster at the 2003 International Sherwood Theory Conference in Corpus Christi, Texas, April 28-30 2003. In the area of code development, we have implemented the MHD equations and the semi-implicit algorithm. The new features have been tested.

Schnack, Dalton D.↗

Assessment of the Accuracy and Consistency of Cloud Property Retrievals from Passive Radiometers with Varying Spectral Information

Over the past few decades, space-based passive radiometers have been used to study the spatial distribution of clouds and some of their basic properties. Modern radiometers have significantly higher spatial resolution and collect more spectral information (i.e., have more spectral bands) than early instruments, and these advances have allowed researchers to characterize cloud physical and microphysical properties with greater detail. There is significant interest in utilizing as much of the existing satellite record as possible to study global and regional changes in cloud occurrence and distribution as well as the associated changes in reflected shortwave and outgoing longwave fluxes. Cloud property retrieval algorithms such as the NASA Satellite Cloud Optical and Radiative Property System (SatCORPS) have traditionally utilized as much spectral information as possible in order to maximize the accuracy of the retrievals. However, changes to the retrieval algorithm over time in order to take advantage of progressively more spectral information results in inconsistencies in cloud properties which in turn results in inconsistencies with products derived from these properties. It is not possible to implement modern state-of-the-art algorithms to earlier instruments which lack diverse spectral information, but many instruments since the early 2000s generally have approximately similar spectral bands in the visible (~0.65 microns), longwave infrared (~11 microns), and in many cases, the midwave infrared (~3.8 microns) and water vapor absorption (~6.7 microns), portions of the electromagnetic spectrum. Different versions of SatCORPS have been implemented with either the full set of spectral information collected by imagers such as GOES-16 and MODIS or with a reduced set of spectral bands in order to be more consistent with past imagers that lack such diverse spectral information. This study uses independent cloud property retrievals from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) to assess the accuracy of the retrievals (e.g., cloud fraction, thermodynamic phase, altitude) from the different retrieval methods and to demonstrate the advantages and disadvantages of using a reduced set of common spectral bands in order to achieve temporal and spatial consistency across different instruments.

Christopher R. Yost↗

Studies in astronomical time series analysis. I - Modeling random processes in the time domain

Several random process models in the time domain are defined and discussed. Attention is given to the moving average model, the autoregressive model, and relationships between and combinations of these models. Consideration is then given to methods for investigating pulse structure, procedures of model construction, computational methods, and numerical experiments. A FORTRAN algorithm of time series analysis has been developed which is relatively stable numerically. Results of test cases are given to study the effect of adding noise and of different distributions for the pulse amplitudes. A preliminary analysis of the light curve of the quasar 3C 272 is considered as an example.

Scargle, J. D.↗

Numerical solution of the incompressible Navier-Stokes equations for steady-state and time-dependent problems

An algorithm for the solution of the incompressible Navier-Stokes equations in three-dimensional generalized curvilinear coordinates is presented. The algorithm can be used to compute both steady-state and time-dependent flow problems. The algorithm is based on the method of artificial compressibility and uses a higher-order flux-difference splitting technique for the convective terms and a second-order central difference for the viscous terms. The steady-state solution of flow through a square duct with a 90 deg bend is computed and the results are compared with experimental data. Good agreement is observed. A comparison with an analytically known exact solution is then performed to verify the time accuracy of the algorithm. Finally, the flow through an artificial heart configuration with moving boundaries is calculated and presented.

Rogers, Stuart E.↗