Search NASA⌕ Search

SEARCH · Search NASA

Results for “iterative regularization”

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 73 records · Page 4

Time-resolved diffusion tomographic imaging in highly scattering turbid media

A method for imaging objects in highly scattering turbid media. According to one embodiment of the invention, the method involves using a plurality of intersecting source/detectors sets and time-resolving equipment to generate a plurality of time-resolved intensity curves for the diffusive component of light emergent from the medium. For each of the curves, the intensities at a plurality of times are then inputted into the following inverse reconstruction algorithm to form an image of the medium: X.sup.(k+1).spsp.T =?Y.sup.T W+X.sup.(k).spsp.T .LAMBDA.!?W.sup.T W+.LAMBDA.!.sup.-1 wherein W is a matrix relating output at detector position r.sub.d, at time t, to source at position r.sub.s, .LAMBDA. is a regularization matrix, chosen for convenience to be diagonal, but selected in a way related to the ratio of the noise, to fluctuations in the absorption (or diffusion) X.sub.j that we are trying to determine: .LAMBDA..sub.ij =.lambda..sub.j .delta..sub.ij with .lambda..sub.j = /<.DELTA.Xj.DELTA.Xj> Here Y is the data collected at the detectors, and X.sup.k is the kth iterate toward the desired absoption information.

Alfano, Robert R.↗

Applications of partial orderings to the study of positive definiteness, monotonicity, and convergence of iterative methods for linear systems.

Most of the known results concerning convergence of iterative methods for solving linear systems involve either positive definiteness or monotonicity. In this paper a new concept, called K-semipositivity, is introduced, which provides a link between convergence theory, monotonicity, and positive definiteness. By using this concept, together with partial orderings in Euclidean n-space, several new convergence theorems are proved. Application to Jacobi's methods and the theory of regular splittings shows the usefulness of these new results.

Vandergraft, J. S.↗

Reducing V&V Cost of Flight Critical Systems: Myth or Reality?

This paper presents an overview of NASA research program on the V&V of flight critical systems. Five years ago, NASA started an effort to reduce the cost and possibly increase the effectiveness of V&V for flight critical systems. It is the right time to take a look back and realize what progress has been made. This paper describes our overall approach and the tools introduced to address different phases of the software lifecycle. For example, we have improved testing by developing a statistical learning approach tor defining test cases. The tool automatically identifies possible unsafe conditions by analyzing outliers in output data; using an iterative learning process, it can then generate more test cases that represent potentially unsafe regions of operation. At the code level, we have developed and made available as open source a static analyzer for C and C++ programs called IKOS. We have shown that IKOS is very precise in the analysis of embedded C programs (very few false positives) and a bit less for regular C and C++ code. At the design level, in collaboration with our NRA partners, we have developed a suite of analysis tools for Simulink models. The analysis is done in a compositional framework for scalability.

Brat, Guillaume P.↗

Inverse solutions for laminar boundary-layer flows with separation and reattachment

Numerical solutions of the laminar, incompressible boundary layer equations are presented for flows involving separation and reattachment. Regular solutions are obtained with an inverse approach in which either the displacement thickness or the skin friction is specified; the pressure is deduced from the solution. A vorticity-stream-function formulation of the boundary layer equations is used to eliminate the unknown pressure. Solutions of the resulting finite difference equations, in which the flow direction is taken into account, are obtained by several global iteration schemes which are stable and have unconditional diagonal dominance. Results are compared with Klineberg and Steger's separated boundary layer calculations, and with Briley's solution of Navier-Stokes equations for a separated region. In addition, an approximate technique is presented in which the streamwise convection of vorticity is set equal to zero in the reversed flow region; such a technique results in a quick forward marching procedure for separated flows.

Carter, J. E.↗

User's manual for tooth contact analysis of face-milled spiral bevel gears with given machine-tool settings

Research was performed to develop a computer program that will: (1) simulate the meshing and bearing contact for face milled spiral beval gears with given machine tool settings; and (2) to obtain the output, some of the data is required for hydrodynamic analysis. It is assumed that the machine tool settings and the blank data will be taken from the Gleason summaries. The theoretical aspects of the program are based on 'Local Synthesis and Tooth Contact Analysis of Face Mill Milled Spiral Bevel Gears'. The difference between the computer programs developed herein and the other one is as follows: (1) the mean contact point of tooth surfaces for gears with given machine tool settings must be determined iteratively, while parameters (H and V) are changed (H represents displacement along the pinion axis, V represents the gear displacement that is perpendicular to the plane drawn through the axes of the pinion and the gear of their initial positions), this means that when V differs from zero, the axis of the pionion and the gear are crossed but not intersected; (2) in addition to the regular output data (transmission errors and bearing contact), the new computer program provides information about the contacting force for each contact point and the sliding and the so-called rolling velocity. The following topics are covered: (1) instructions for the users as to how to insert the input data; (2) explanations regarding the output data; (3) numerical example; and (4) listing of the program.

Litvin, Faydor L.↗

Building Automation System Replacement and Improvements during Concurrent Integration and Testing for the Spacecraft Systems Development and Integration Facility

The Spacecraft Systems Development and Integration Facility (SSDIF), the premier cleanroom for NASA Goddard Space Flight Center (GSFC), has gone through multiple iterations of changes to its building automation system (BAS). NASA flight projects require more stringent environmental controls for temperature, humidity, air exchange rates, and pressurization. This led to increasingly difficult circumstances with replacing and improving the BAS during Integration and Testing (I&T) operations. The SSDIF was partially transitioned to a new BAS during November and December of 2023. Transition occurred in a live environment, while regular I&T activities resumed for the Roman Space Telescope. This resulted in the development and execution of a comprehensive method of procedure. Strategies such as sequencing equipment transition, and using temporary controls for critical equipment, avoided detriments to cleanroom environment. While the new BAS was being built, data points on the previous system would become obsolete, leading to data gaps and uncertainties. Critical data points were distinguished and trended on the new system, providing key insights on equipment functionality, while the BAS was completed. Control hunting and environmental deviations were persistent cases with the SSDIF. In some instances, proportional-integral-derivative (PID) tuning was implemented on current logic to bring smooth PID loop hunting. In other situations, controls logic was completely altered by changing setpoints or sequences. This provided greater laminar flow, more stable pressurization, and improved temperature and humidity controls. This paper documents the transitionary procedure to a new BAS, the improvements made, and challenges faced when altering a BAS during I&T operations.

Cleanroom↗

Solutions for laminar boundary layers with separation and reattachment

Numerical solutions of the laminar, incompressible boundary-layer equations are presented for flows involving separation and reattachment. Regular solutions are obtained with an inverse approach in which either the displacement thickness or the skin friction is specified, and the pressure is deduced from the solution. A stream function/vorticity formulation of the boundary-layer equations is used to eliminate the unknown pressure. Solutions of the resulting finite-difference equations, in which the flow direction is taken into account, are obtained by a successive column iteration scheme. Results are compared with Klineberg and Steger's (1974) separated boundary-layer calculations, and with Briley's (1971) solution of the Navier-Stokes equations for a separated region.

Carter, J. E.↗

Development of a Daily GRACE Mascon Solution for Terrestrial Water Storage

The Gravity Recovery and Climate Experiment (GRACE) and GRACE-Follow On missions have provided a global history of terrestrial water storage changes since 2002. Traditional GRACE products resolve monthly time-variable gravity at spatial resolutions of 300-500 km, with many recent efforts focusing on regularized mass concentration (mascon) solutions to better resolve signals spatially. However, monthly resolution inhibits the applicability of GRACE to investigations of sub-monthly signals. This study presents a new daily mass change solution estimated as deviations from the NASA Goddard Space Flight Center (GSFC) converged monthly mascon product and quantifies the fundamental trade-off between temporal and spatial resolution in GRACE-only solutions. As an iteration of the GSFC product, this daily solution represents a mixing of monthly information at higher relative spatial resolution with daily information at lower relative spatial resolution, resolving each temporal scale at the best resolution achievable without requiring hydrological model dependencies in the estimation. The resolution of the daily mascons is 300-400 km at high latitudes and 600-1,000 km lower, depending on each daily orbital track and the proximity of each mascon to constraint region boundaries (e.g., coastlines). Through simulations and model comparisons, we demonstrate that daily signals over areas larger than 400,000 km(2) are recoverable, with basins larger than 800,000 km(2) exhibiting strong signal recovery relative to leakage. This analysis establishes baselines for daily signal recovery from GRACE in the context of longer time scale solutions, characterizes leakage inherent in daily GRACE information, and creates new opportunities for applying GRACE to investigations of sub-monthly signals.

M. J. Croteau↗

The GEOS Retrospective Data Assimilation System: The 6-hour lag case

The fixed-lag Kalman smoother (FLKS) has been proposed as a framework to construct data assimilation procedures capable of producing high-quality climate research datasets. Fixed-lag Kalman smoother-based systems, referred to as retrospective data assimilation systems, are an extension to three-dimensional filtering procedures with the added capability of incorporating observations not only in the past and present time of the estimate, but also at future times. A variety of simplifications are necessary to render retrospective assimilation procedures practical. In this article, we present an FLKS-based retrospective data assimilation system implementation for the Goddard Earth Observing System (GOES) Data Assimilation System (DAS). The practicality of this implementation comes from the practicality of its underlying (filter) analysis system, i.e., the physical-space statistical analysis system (PSAS). The behavior of two schemes is studied here. The first retrospective analysis (RA) scheme is designed simply to update the regular PSAS analyses with observations available at times ahead of the regular analysis times. Although our GEOS DAS implementation is general, results are only presented for when observations 6-hours ahead of the analysis time are used to update the PSAS analyses and thereby to calculate the so-called lag-1 retrospective analyses. Consistency tests for this RA scheme show that the lag-1 retrospective analyses indeed have better 6-hour predictive skills than the predictions from the regular analyses. This motivates the introduction of the second retrospective analysis scheme which, at each analysis time, uses the 6-hour retrospective analysis to replace the first-guess normally used in the PSAS analysis, and therefore allows the calculation of a revised (filter) PSAS analysis. Since in this scheme the lag-1 retrospective analyses influence the filter results, this procedure is referred to as the retrospective-based iterative analysis (RIA) scheme. Results from the RIA scheme indicate its potential for improving the overall quality of the assimilation.

Zhu, Yan-Qiu↗

Parallel architectures for iterative methods on adaptive, block structured grids

A parallel computer architecture well suited to the solution of partial differential equations in complicated geometries is proposed. Algorithms for partial differential equations contain a great deal of parallelism. But this parallelism can be difficult to exploit, particularly on complex problems. One approach to extraction of this parallelism is the use of special purpose architectures tuned to a given problem class. The architecture proposed here is tuned to boundary value problems on complex domains. An adaptive elliptic algorithm which maps effectively onto the proposed architecture is considered in detail. Two levels of parallelism are exploited by the proposed architecture. First, by making use of the freedom one has in grid generation, one can construct grids which are locally regular, permitting a one to one mapping of grids to systolic style processor arrays, at least over small regions. All local parallelism can be extracted by this approach. Second, though there may be a regular global structure to the grids constructed, there will be parallelism at this level. One approach to finding and exploiting this parallelism is to use an architecture having a number of processor clusters connected by a switching network. The use of such a network creates a highly flexible architecture which automatically configures to the problem being solved.

Gannon, D.↗

Software for Refining or Coarsening Computational Grids

A computer program performs calculations for refinement or coarsening of computational grids of the type called 'structured' (signifying that they are geometrically regular and/or are specified by relatively simple algebraic expressions). This program is designed to facilitate analysis of the numerical effects of changing structured grids utilized in computational fluid dynamics (CFD) software. Unlike prior grid-refinement and -coarsening programs, this program is not limited to doubling or halving: the user can specify any refinement or coarsening ratio, which can have a noninteger value. In addition to this ratio, the program accepts, as input, a grid file and the associated restart file, which is basically a file containing the most recent iteration of flow-field variables computed on the grid. The program then refines or coarsens the grid as specified, while maintaining the geometry and the stretching characteristics of the original grid. The program can interpolate from the input restart file to create a restart file for the refined or coarsened grid. The program provides a graphical user interface that facilitates the entry of input data for the grid-generation and restart-interpolation routines.

Daines, Russell↗

Software for Refining or Coarsening Computational Grids

A computer program performs calculations for refinement or coarsening of computational grids of the type called "structured" (signifying that they are geometrically regular and/or are specified by relatively simple algebraic expressions). This program is designed to facilitate analysis of the numerical effects of changing structured grids utilized in computational fluid dynamics (CFD) software. Unlike prior grid-refinement and -coarsening programs, this program is not limited to doubling or halving: the user can specify any refinement or coarsening ratio, which can have a noninteger value. In addition to this ratio, the program accepts, as input, a grid file and the associated restart file, which is basically a file containing the most recent iteration of flow-field variables computed on the grid. The program then refines or coarsens the grid as specified, while maintaining the geometry and the stretching characteristics of the original grid. The program can interpolate from the input restart file to create a restart file for the refined or coarsened grid. The program provides a graphical user interface that facilitates the entry of input data for the grid-generation and restart-interpolation routines.

Daines, Russell↗

Software for Refining or Coarsening Computational Grids

A computer program performs calculations for refinement or coarsening of computational grids of the type called structured (signifying that they are geometrically regular and/or are specified by relatively simple algebraic expressions). This program is designed to facilitate analysis of the numerical effects of changing structured grids utilized in computational fluid dynamics (CFD) software. Unlike prior grid-refinement and -coarsening programs, this program is not limited to doubling or halving: the user can specify any refinement or coarsening ratio, which can have a noninteger value. In addition to this ratio, the program accepts, as input, a grid file and the associated restart file, which is basically a file containing the most recent iteration of flow-field variables computed on the grid. The program then refines or coarsens the grid as specified, while maintaining the geometry and the stretching characteristics of the original grid. The program can interpolate from the input restart file to create a restart file for the refined or coarsened grid. The program provides a graphical user interface that facilitates the entry of input data for the grid-generation and restart-interpolation routines.

Daines, Russell↗

A Highly Accurate Technique for the Treatment of Flow Equations at the Polar Axis in Cylindrical Coordinates using Series Expansions

Numerical methods for solving the flow equations in cylindrical or spherical coordinates should be able to capture the behavior of the exact solution near the regions where the particular form of the governing equations is singular. In this work we focus on the treatment of these numerical singularities for finite-differences methods by reinterpreting the regularity conditions developed in the context of pseudo-spectral methods. A generally applicable numerical method for treating the singularities present at the polar axis, when nonaxisymmetric flows are solved in cylindrical, coordinates using highly accurate finite differences schemes (e.g., Pade schemes) on non-staggered grids, is presented. Governing equations for the flow at the polar axis are derived using series expansions near r=0. The only information needed to calculate the coefficients in these equations are the values of the flow variables and their radial derivatives at the previous iteration (or time) level. These derivatives, which are multi-valued at the polar axis, are calculated without dropping the accuracy of the numerical method using a mapping of the flow domain from (0,R)*(0,2pi) to (-R,R)*(0,pi), where R is the radius of the computational domain. This allows the radial derivatives to be evaluated using high-order differencing schemes (e.g., compact schemes) at points located on the polar axis. The proposed technique is illustrated by results from simulations of laminar-forced jets and turbulent compressible jets using large eddy simulation (LES) methods. In term of the general robustness of the numerical method and smoothness of the solution close to the polar axis, the present results compare very favorably to similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r=0, following the method proposed recently by Mohseni and Colonius (1). Extension of the method described here for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical or spherical coordinates with finite-differences schemes of various level of accuracy is immediate.

Constantinescu, George S.↗

Efficacy of Code Optimization on Cache-based Processors

The current common wisdom in the U.S. is that the powerful, cost-effective supercomputers of tomorrow will be based on commodity (RISC) micro-processors with cache memories. Already, most distributed systems in the world use such hardware as building blocks. This shift away from vector supercomputers and towards cache-based systems has brought about a change in programming paradigm, even when ignoring issues of parallelism. Vector machines require inner-loop independence and regular, non-pathological memory strides (usually this means: non-power-of-two strides) to allow efficient vectorization of array operations. Cache-based systems require spatial and temporal locality of data, so that data once read from main memory and stored in high-speed cache memory is used optimally before being written back to main memory. This means that the most cache-friendly array operations are those that feature zero or unit stride, so that each unit of data read from main memory (a cache line) contains information for the next iteration in the loop. Moreover, loops ought to be 'fat', meaning that as many operations as possible are performed on cache data-provided instruction caches do not overflow and enough registers are available. If unit stride is not possible, for example because of some data dependency, then care must be taken to avoid pathological strides, just ads on vector computers. For cache-based systems the issues are more complex, due to the effects of associativity and of non-unit block (cache line) size. But there is more to the story. Most modern micro-processors are superscalar, which means that they can issue several (arithmetic) instructions per clock cycle, provided that there are enough independent instructions in the loop body. This is another argument for providing fat loop bodies. With these restrictions, it appears fairly straightforward to produce code that will run efficiently on any cache-based system. It can be argued that although some of the important computational algorithms employed at NASA Ames require different programming styles on vector machines and cache-based machines, respectively, neither architecture class appeared to be favored by particular algorithms in principle. Practice tells us that the situation is more complicated. This report presents observations and some analysis of performance tuning for cache-based systems. We point out several counterintuitive results that serve as a cautionary reminder that memory accesses are not the only factors that determine performance, and that within the class of cache-based systems, significant differences exist.

VanderWijngaart, Rob F.↗

Deconvolution of SNPP VIIRS Solar Diffuser Bidirectional Reflectance Distribution Function On-orbit Change Factor

The Earth-observing Visible Infrared Imaging Radiometer Suite (VIIRS) on the Suomi National Polar-orbiting Partnership satellite regularly calibrates its reflective solar bands (RSBs), primarily through observing an onboard sunlit solar diffuser (SD). The on-orbit change of the value of the SD bidirectional reflectance distribution function (BRDF) is quantified by a numerical factor, called the H-factor, and is determined by the onboard SD stability monitor (SDSM). Because the spectral response function of an SDSM detector spreads in wavelength, the directly measured H-factor is the true H-factor convolved with the spectral response function. To find the true H-factor, we use the traditional direct method and an innovative iterative approach to separately deconvolve the measured H-factor. Our iterative approach relies on two properties of the SDSM detector spectral response function: the central peak width is narrow enough so that the H-factor does not change much over the peak width, and the dominance of the spectral response function’s integral with respect to the wavelength over the width. The iterative approach is more accurate, of a smaller noise impact, much more flexible in terms of interpolation and extrapolation of function values, and faster. We have used deconvolved H-factors to calibrate the NASA SNPP VIIRS RSB Collections 1 and 2 Level-1B products.

SNPP VIIRS↗

Cyclic Oxidation Modeling and Life Prediction

The cyclic oxidation process can be described as an iterative scale growth and spallation sequence by a number of similar models. Model input variable include oxide scale type and growth parameters, spalling geometry, spall constant, and cycle duration. Outputs include net weight change, the amounts of retained and spalled oxide, the total oxygen and metal consumed, and the terminal rates of weight loss and metal consumption. All models and their variations produce a number of similar characteristic features. In general, spalling and material consumption increase to a steady state rate, at which point the retained scale approaches a constant and the rate of weight loss becomes linear. For one model, this regularity was demonstrated as dimensionless, universal expressions, obtained by normalizing the variables by critical performance factors. These insights were enabled through the use of the COSP for Windows cyclic oxidation spalling program.

Smialek, James L.↗

Multilevel Algorithm for Atmospheric Data Assimilation

A multiscale algorithm for the problem of optimal statistical interpolation of observed data has been developed. This problem includes the calculation of the vector of the 'analyzed' (best estimated) atmosphere flow field w(sup a) by the formula: w(sup a) = w(sup f) + P(sup f) H(sup T) y, where the quantity y is defined by the equation (H P(sup f) H(sup T) + R)y = w(sup o) - H w(sup f), using the given model forecast first guess w(sup f) and the vector of observations w(sup o); H is an interpolation operator from the regular grid to the observation network, P(sup f) is the forecast error covariance matrix, and R is the observation error covariance matrix. At this initial stage the case of univariate analysis of single level radiosonde height data is considered. The matrix R is assumed to be diagonal, and the matrix P(sup f) is assumed to be given by the formula P(sub ij)(sup f) = sigma(sub i)(sup f) mu(sub ij) sigma(sub j)(sub f), where mu(sub ij) is a smooth, decreasing function of the distance between the i-th and the j-th points. In this paper we describe a multiscale iterative process based on a multiresolution, simultaneous displacement technique and a localized variational calculation of iteration parameters.

Brandt, Achi↗