Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized algorithm”

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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Application of the Loci-Based CFD Code Chem at MSFC: Preliminary Results

Contents include the following: 1. Objectives. Concentrate on determining the qualitative accuracy, performance and robustness of the Chen code. 2. What is the Loci-Chem CFD code? Density-based, finite volume, generalized unstructured grid, Navier-Stokes solver. The algorithm was implemented using the Loci framework, which allows implementation issues such as parallel processing to be handled transparently to the coding of the CFD algorithm. 3. Application to Bifurcating Duct problem. Flow splits from single duct to two ducts. 4. Application to single element injector. 5. Application to PSU RBCC rig. 6. 90 degree elbow benchmark problem. 7. Future work.

West, Jeff S.↗
NASA NTRS2005

Quadratic Programming for Allocating Control Effort

A computer program calculates an optimal allocation of control effort in a system that includes redundant control actuators. The program implements an iterative (but otherwise single-stage) algorithm of the quadratic-programming type. In general, in the quadratic-programming problem, one seeks the values of a set of variables that minimize a quadratic cost function, subject to a set of linear equality and inequality constraints. In this program, the cost function combines control effort (typically quantified in terms of energy or fuel consumed) and control residuals (differences between commanded and sensed values of variables to be controlled). In comparison with prior control-allocation software, this program offers approximately equal accuracy but much greater computational efficiency. In addition, this program offers flexibility, robustness to actuation failures, and a capability for selective enforcement of control requirements. The computational efficiency of this program makes it suitable for such complex, real-time applications as controlling redundant aircraft actuators or redundant spacecraft thrusters. The program is written in the C language for execution in a UNIX operating system.

Singh, Gurkirpal↗
NASA NTRS2016

New Satellite Project Aerosol-UA: Remote Sensing of Aerosols in the Terrestrial Atmosphere

We discuss the development of the Ukrainian space project Aerosol-UA which has the following three main objectives: (1) to monitor the spatial distribution of key characteristics of terrestrial tropospheric and stratospheric aerosols; (2) to provide a comprehensive observational database enabling accurate quantitative estimates of the aerosol contribution to the energy budget of the climate system; and (3) quantify the contribution of anthropogenic aerosols to climate and ecological processes. The remote sensing concept of the project is based on precise orbital measurements of the intensity and polarization of sunlight scattered by the atmosphere and the surface with a scanning polarimeter accompanied by a wide-angle multispectral imager-polarimeter. Preparations have already been made for the development of the instrument suite for the Aerosol-UA project, in particular, of the multi-channel scanning polarimeter (ScanPol) designed for remote sensing studies of the global distribution of aerosol and cloud properties (such as particle size, morphology, and composition) in the terrestrial atmosphere by polarimetric and spectrophotometric measurements of the scattered sunlight in a wide range of wavelengths and viewing directions from which a scene location is observed. ScanPol is accompanied by multispectral wide-angle imager-polarimeter (MSIP) that serves to collect information on cloud conditions and Earths surface image. Various components of the polarimeter ScanPol have been prototyped, including the opto-mechanical and electronic assemblies and the scanning mirror controller. Preliminary synthetic data simulations for the retrieval of aerosol parameters over land surfaces have been performed using the Generalized Retrieval of Aerosol and Surface Properties (GRASP) algorithm. Methods for the validation of satellite data using ground-based observations of aerosol properties are also discussed. We assume that designing, building, and launching into orbit a multi-functional high-precision scanning polarimeter and an imager-polarimeter should make a significant contribution to the study of natural and anthropogenic aerosols and their climatic and ecological effects.

remote sensing↗
NASA NTRS2016

Atmospheric Correction for Satellite Ocean Color Radiometry

This tutorial is an introduction to atmospheric correction in general and also documentation of the atmospheric correction algorithms currently implemented by the NASA Ocean Biology Processing Group (OBPG) for processing ocean color data from satellite-borne sensors such as MODIS and VIIRS. The intended audience is graduate students or others who are encountering this topic for the first time. The tutorial is in two parts. Part I discusses the generic atmospheric correction problem. The magnitude and nature of the problem are first illustrated with numerical results generated by a coupled ocean-atmosphere radiative transfer model. That code allow the various contributions (Rayleigh and aerosol path radiance, surface reflectance, water-leaving radiance, etc.) to the topof- the-atmosphere (TOA) radiance to be separated out. Particular attention is then paid to the definition, calculation, and interpretation of the so-called "exact normalized water-leaving radiance" and its equivalent reflectance. Part I ends with chapters on the calculation of direct and diffuse atmospheric transmittances, and on how vicarious calibration is performed. Part II then describes one by one the particular algorithms currently used by the OBPG to effect the various steps of the atmospheric correction process, viz. the corrections for absorption and scattering by gases and aerosols, Sun and sky reflectance by the sea surface and whitecaps, and finally corrections for sensor out-of-band response and polarization effects. One goal of the tutorial-guided by teaching needs- is to distill the results of dozens of papers published over several decades of research in atmospheric correction for ocean color remote sensing.

MODIS↗
NASA NTRS2017

NASA Computational Case Study: Spectral Energy Distribution Fitting

The need for faster, more efficient algorithms is an important aspect of scientific computing. Generally, scientists are only exposed to computational issues that arise in their field. Thus, collaboration between a numerical analyst and a scientist is becoming necessary for scientific computing. The purpose of this case study is to expose computer scientists to processes that an astronomer would use to obtain useful results from raw data. For example, astronomers are interested in determining the properties of galaxies and measuring changes in those properties as a function of time throughout cosmic history. To do so, they use certain models that are designed and refined over time via observations at different wavelengths of the light spectrum. The process of matching these models with observed data from studying celestial bodies is referred to as Spectral Energy Distribution (SED) fitting. In this case study, we learn how to perform the SED fit. This process requires knowledge of both the astronomical and computational issues involved when fitting flux, the total energy from a source as seen from Earth, to a set of physical templates. Once a fit is complete, one can classify the source and estimate a number of physical parameters. The goal is to demonstrate how the computer science skill set can be used in the scientific community and to possibly improve one or more of the computational aspects of this problem. The following provides some background on the astronomy issues, including information on the spectral energy distribution and related physical parameters, and the computational issues, including the fitting procedure..

modeling↗
NASA NTRS2021

Validation of the aerosol optical property products derived by the GRASP/Component approach from multi-angular polarimetric observations

The multi-angle polarimetric satellite observations are helpful for improving the retrievals of aerosol parameters. However, practical applications of polarization technology are still limited because of complexity of measurement and interpretation of polarimetric observations. In this study, we analyze the performance of a new component approach developed in the frame of the Generalized Retrieval of Atmosphere and Surface Properties (GRASP) algorithm. In addition to aerosol optical properties including particle size distribution, non-sphericity and index of refraction that are commonly derived from multi-angle radiance and polarization measurements, the GRASP/Component approach also provides some information about aerosol composition. Specifically, in this approach aerosol is modeled as an internal mixture of several components with distinctly different chemical compositions and known refractive indices. The approach is intended not only to provide additional insight on aerosol composition but also to improve retrieval of basic aerosol optical properties. This study presents comprehensive validation and evaluation of Aerosol Optical Depth (AOD), Ångstrom exponent (AE), fine mode AOD (AODF), coarse mode AOD (AODC), and single scattering albedo (SSA) as retrieved by the GRASP/Component approach. The GRASP/Component products include aerosol retrievals using two different aerosol component mixing rules, i.e., Maxwell-Garnett (MG) effective medium approximation and a simple Volume- Weighted averaging (VW). The differences between the results obtained using these two assumptions are also discussed. The obtained results show that the aerosol optical property products of GRASP/Component approach have good agreement with the ground-based AERONET measurements, which is comparable to other PARASOL/ GRASP approaches. Specifically, the AOD retrieved by GRASP/Component approach show high correlation and nearly no bias both over land and ocean, as compared with AERONET. The more detailed aerosol properties such as AE, AODF, AODC and SSA also show one of the best comparisons with AERONET. These improvements can probably be attributed to the use of the additional physical constraints on spectral dependence of the complex refractive index and the reduction of total number of aerosol parameters directly retrieved in the GRASP/Component approach. In addition, the choice of mixing rules had no significant effect on optical retrievals. With the exception of SSA, the results obtained based on the MG mixing rule were found to be slightly better over those obtained using VW mixing rule, especially for bias.

Brent N Holben↗
NASA NTRS2022

Local Scale (3-M) Soil Moisture Mapping Using SMAP and Planet Superdove

A capability for mapping meter-level resolution soil moisture with frequent temporal sampling over large regions is essential for quantifying local-scale environmental heterogeneity and eco-hydrologic behavior. However, available surface soil moisture (SSM) products generally involve much coarser grain sizes ranging from 30 m to several 10s of kilometers. Hence a new method is proposed to estimate 3-m resolution SSM using a combination of multi-sensor fusion, machine- learning (ML) and Cumulative Distribution Function (CDF) matching approaches. This method established favorable SSM correspondence between 3-m pixels and overlying 9-km grid cells from overlapping Planet SuperDove (PSD) observations and NASA Soil Moisture Active-Passive (SMAP) mission products. The resulting 3-m SSM predictions showed improved accuracy by reducing ab- solute bias and RMSE by ~0.01 cm3/cm3 over the original SMAP data in relation to in-situ soil moisture measurements for the Australian Yanco region, while preserving the high sampling frequency (1-3 day global revisit) and sensitivity to surface wetness (R 0.865) from SMAP. Heterogeneous soil moisture distributions varying with vegetation biomass gradients and irrigation regimes were generally captured within a selected study area. Further algorithm refinement and implementation for regional applications will allow for improvement in water resources management, precision agriculture, and disaster forecasts and responses.

soil moisture↗
DOE OSTI2021

Real time evolution for ultracompact Hamiltonian eigenstates on quantum hardware

In this work we present a detailed analysis of variational quantum phase estimation (VQPE), a method based on real-time evolution for ground and excited state estimation on near-term hardware. We derive the theoretical ground on which the approach stands, and demonstrate that it provides one of the most compact variational expansions to date for solving strongly correlated Hamiltonians. At the center of VQPE lies a set of equations, with a simple geometrical interpretation, which provides conditions for the time evolution grid in order to decouple eigenstates out of the set of time evolved expansion states, and connects the method to the classical filter diagonalization algorithm. Further, we introduce what we call the unitary formulation of VQPE, in which the number of matrix elements that need to be measured scales linearly with the number of expansion states, and we provide an analysis of the effects of noise which substantially improves previous considerations. The unitary formulation allows for a direct comparison to iterative phase estimation. Our results mark VQPE as both a natural and highly efficient quantum algorithm for ground and excited state calculations of general many-body systems. We demonstrate a hardware implementation of VQPE for the transverse field Ising model. Further, we illustrate its power on a paradigmatic example of strong correlation (Cr2 in the SVP basis set), and show that it is possible to reach chemical accuracy with as few as ~50 timesteps.

Klymko, Katherine↗
NASA NTRS1992

Laboratory for Engineering Man/Machine Systems (LEMS): System identification, model reduction and deconvolution filtering using Fourier based modulating signals and high order statistics

Several important problems in the fields of signal processing and model identification, such as system structure identification, frequency response determination, high order model reduction, high resolution frequency analysis, deconvolution filtering, and etc. Each of these topics involves a wide range of applications and has received considerable attention. Using the Fourier based sinusoidal modulating signals, it is shown that a discrete autoregressive model can be constructed for the least squares identification of continuous systems. Some identification algorithms are presented for both SISO and MIMO systems frequency response determination using only transient data. Also, several new schemes for model reduction were developed. Based upon the complex sinusoidal modulating signals, a parametric least squares algorithm for high resolution frequency estimation is proposed. Numerical examples show that the proposed algorithm gives better performance than the usual. Also, the problem was studied of deconvolution and parameter identification of a general noncausal nonminimum phase ARMA system driven by non-Gaussian stationary random processes. Algorithms are introduced for inverse cumulant estimation, both in the frequency domain via the FFT algorithms and in the domain via the least squares algorithm.

Pan, Jianqiang↗
NASA NTRS2012

TRMM Version 7 Level 3 Gridded Monthly Accumulations of GPROF Precipitation Retrievals

In July 2011, improved versions of the retrieval algorithms were approved for TRMM. All data starting with June 2011 are produced only with the version 7 code. At the same time, version 7 reprocessing of all TRMM mission data was started. By the end of August 2011, the 14+ years of the reprocessed mission data became available online to users. This reprocessing provided the opportunity to redo and enhance upon an analysis of V7 impacts on L3 data accumulations that was presented at the 2010 EGU General Assembly. This paper will discuss the impact of algorithm changes made in th GPROF retrieval on the Level 2 swath products. Perhaps the most important change in that retrieval was to replacement of a model based a priori database with one created from Precipitation Radar (PR) and TMI brightness temperature (Tb) data. The radar pays a major role in the V7 GPROF (GPROF2010) in determining existence of rain. The level 2 retrieval algorithm also introduced a field providing the probability of rain. This combined use of the PR has some impact on the retrievals and created areas, particularly over ocean, where many areas of low-probability precipitation are retrieved whereas in version 6, these areas contained zero rain rates. This paper will discuss how these impacts get translated to the space/time averaged monthly products that use the GPROF retrievals. The level 3 products discussed are the gridded text product 3G68 and the standard 3A12 and 3B31 products. The paper provides an overview of the changes and explanation of how the level 3 products dealt with the change in the retrieval approach. Using the .25 deg x .25 degree grid, the paper will show that agreement between the swath product and the level 3 remains very high. It will also present comparisons of V6 and V7 GPROF retrievals as seen both at the swath level and the level 3 time/space gridded accumulations. It will show that the various L3 products based on GPROF level 2 retrievals are in close agreement. The paper concludes by outlining some of the challenges of the TRMM version 7 level 3 products.

Stocker, E. F.↗
NASA NTRS2006

Solver and Turbulence Model Upgrades to OVERFLOW 2 for Unsteady and High-Speed Applications

An implicit unfactored SSOR algorithm has been added to the overset Navier-Stokes CFD code OVERFLOW 2 for unsteady and moving body applications. The HLLEM and HLLC third-order spatial upwind convective flux models have been added for high-speed flow applications. A generalized upwind transport equation has been added for solution of the two-equation turbulence models and the species equations. The generalized transport equation is solved using an unfactored SSOR implicit algorithm. Three hybrid RANS/DES turbulence models have been added for unsteady flow applications. Wall function boundary conditions that include compressibility and heat transfer effects have been also been added to OVERFLOW 2.

Nichols, Robert H.↗
NASA NTRS1990

Partitioning sparse matrices with eigenvectors of graphs

The problem of computing a small vertex separator in a graph arises in the context of computing a good ordering for the parallel factorization of sparse, symmetric matrices. An algebraic approach for computing vertex separators is considered in this paper. It is shown that lower bounds on separator sizes can be obtained in terms of the eigenvalues of the Laplacian matrix associated with a graph. The Laplacian eigenvectors of grid graphs can be computed from Kronecker products involving the eigenvectors of path graphs, and these eigenvectors can be used to compute good separators in grid graphs. A heuristic algorithm is designed to compute a vertex separator in a general graph by first computing an edge separator in the graph from an eigenvector of the Laplacian matrix, and then using a maximum matching in a subgraph to compute the vertex separator. Results on the quality of the separators computed by the spectral algorithm are presented, and these are compared with separators obtained from other algorithms for computing separators. Finally, the time required to compute the Laplacian eigenvector is reported, and the accuracy with which the eigenvector must be computed to obtain good separators is considered. The spectral algorithm has the advantage that it can be implemented on a medium-size multiprocessor in a straightforward manner.

Pothen, Alex↗
DOE OSTI2022

Quantum advantage for noisy channel discrimination

Many quantum mechanical experiments can be viewed as multiround interactive protocols between known quantum circuits and an unknown quantum process. Fully quantum “coherent” access to the unknown process is known to provide an advantage in many discrimination tasks compared to when only incoherent access is permitted, but it is unclear if this advantage persists when the process is noisy. Here, we show that a quantum advantage can be maintained when distinguishing between two noisy single-qubit rotation channels. Numerical and analytical calculations reveal a distinct transition between the performance of fully coherent and fully incoherent protocols as a function of noise strength. Moreover, the size of the region of coherent quantum advantage shrinks inverse polynomially in the number of channel uses, and in an intermediate regime an improved strategy is a hybrid of fully coherent and fully incoherent subroutines. The fully coherent protocol is based on quantum signal processing, suggesting a generalizable algorithmic framework for the study of quantum advantage in the presence of realistic noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1992

An extension to Schneider's general paradigm for fault-tolerant clock synchronization

In 1987, Schneider presented a general paradigm that provides a single proof of a number of fault tolerant clock synchronization algorithms. His proof was subsequently subjected to the rigor of mechanical verification by Shankar. However, both Schneider and Shankar assumed a condition Shankar refers to as a bounded delay. This condition states that the elapsed time between synchronization events (i.e., the time that the local process applies an adjustment to its logical clock) is bounded. This property is really a result of the algorithm and should not be assumed in a proof of correctness. This paper remedies this by providing a proof of this property in the context of the general paradigm proposed by Schneider. The argument given is a generalization of Welch and Lynch's proof of a related property for their algorithm.

Miner, Paul S.↗
NASA NTRS1977

Evaluation of algorithms for geological thermal-inertia mapping

The errors incurred in producing a thermal inertia map are of three general types: measurement, analysis, and model simplification. To emphasize the geophysical relevance of these errors, they were expressed in terms of uncertainty in thermal inertia and compared with the thermal inertia values of geologic materials. Thus the applications and practical limitations of the technique were illustrated. All errors were calculated using the parameter values appropriate to a site at the Raft River, Id. Although these error values serve to illustrate the magnitudes that can be expected from the three general types of errors, extrapolation to other sites should be done using parameter values particular to the area. Three surface temperature algorithms were evaluated: linear Fourier series, finite difference, and Laplace transform. In terms of resulting errors in thermal inertia, the Laplace transform method is the most accurate (260 TIU), the forward finite difference method is intermediate (300 TIU), and the linear Fourier series method the least accurate (460 TIU).

Miller, S. H.↗
NASA NTRS1978

Development of a three-dimensional Navier-Stokes code on CDC star-100 computer

A three-dimensional code in body-fitted coordinates was developed using MacCormack's algorithm. The code is structured to be compatible with any general configuration, provided that the metric coefficients for the transformation are available. The governing equations are developed in primitive variables in order to facilitate the incorporation of physical boundary conditions and turbulence-closure models. MacCormack's two-step, unsplit, time-marching algorithm is used to solve the unsteady Navier-Stokes equations until steady-state solution is achieved. Cases discussed include (1) flat plate in supersonic free stream; (2) supersonic flow along an axial corner; (3) subsonic flow in an axial corner at M infinity = 0.95; and (4) supersonic flow in an axial corner at M infinity 1.5.

Vatsa, V. N.↗
NASA NTRS2015

An Efficient Solution Method for Multibody Systems with Loops Using Multiple Processors

This paper describes a multibody dynamics algorithm formulated for parallel implementation on multiprocessor computing platforms using the divide-and-conquer approach. The system of interest is a general topology of rigid and elastic articulated bodies with or without loops. The algorithm divides the multibody system into a number of smaller sets of bodies in chain or tree structures, called "branches" at convenient joints called "connection points", and uses an Order-N (O (N)) approach to formulate the dynamics of each branch in terms of the unknown spatial connection forces. The equations of motion for the branches, leaving the connection forces as unknowns, are implemented in separate processors in parallel for computational efficiency, and the equations for all the unknown connection forces are synthesized and solved in one or several processors. The performances of two implementations of this divide-and-conquer algorithm in multiple processors are compared with an existing method implemented on a single processor.

Multibody dynamics↗