Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative limiters”

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 91 records · Page 5

Space Station Freedom electric power system photovoltaic power module integrated launch package

The launch of the Space Station Freedom solar power module requires a weight efficient structure that will include large components within the limited load capacity of the Space Shuttle cargo bay. The design iterations to meet these requirements have evolved from a proposal concept featuring a separate cradle and integrated equipment assembly (IEA), to a package that interfaces directly with the Shuttle. Size, weight, and cost have been reduced as a result.

Nathanson, Theodore H.↗

Interferometric reconstruction of three-dimensional high-speed aerodynamic flows

Holographic interferometry can be a very useful diagnostic tool in high-speed aerodynamic testing. During this summer research period, various possible approaches for accurately reconstructing three-dimensional flows from limited data were examined. The approach based on the combination of the following three techniques appears to be promising: (1) Continuous Local Basis Function Method - this computational tomographic method has a power to accurately reconstruct continuous regions and is appropriate from well-conditioned to moderately limited data; (2) Variable Basis Method - this computational tomographic method provides accuracy near discontinuities, i.e., shock regions, and is appropriate from moderately-limited to severely-limited data; and (3) Complementary Field Method - this is a general iterative reconstructor that can be coupled with any computational tomographic techniques. Mathematically, it can be shown that this method can provide better accuracy than the direct reconstruction as in a conventional approach. Our numerical simulation of experiments demonstrated improved reconstruction results even when these techniques were individually tested.

Cha, Soyoung Stephen↗

Tractable Pareto Optimization of Temporal Preferences

This paper focuses on temporal constraint problems where the objective is to optimize a set of local preferences for when events occur. In previous work, a subclass of these problems has been formalized as a generalization of Temporal CSPs, and a tractable strategy for optimization has been proposed, where global optimality is defined as maximizing the minimum of the component preference values. This criterion for optimality, which we call 'Weakest Link Optimization' (WLO), is known to have limited practical usefulness because solutions are compared only on the basis of their worst value; thus, there is no requirement to improve the other values. To address this limitation, we introduce a new algorithm that re-applies WLO iteratively in a way that leads to improvement of all the values. We show the value of this strategy by proving that, with suitable preference functions, the resulting solutions are Pareto Optimal.

Morris, Robert↗

Efficiency and Accuracy of Time-Accurate Turbulent Navier-Stokes Computations

The accuracy and efficiency of two types of subiterations in both explicit and implicit Navier-Stokes codes are explored for unsteady laminar circular-cylinder flow and unsteady turbulent flow over an 18-percent-thick circular-arc (biconvex) airfoil. Grid and time-step studies are used to assess the numerical accuracy of the methods. Nonsubiterative time-stepping schemes and schemes with physical time subiterations are subject to time-step limitations in practice that are removed by pseudo time sub-iterations. Computations for the circular-arc airfoil indicate that a one-equation turbulence model predicts the unsteady separated flow better than an algebraic turbulence model; also, the hysteresis with Mach number of the self-excited unsteadiness due to shock and boundary-layer separation is well predicted.

Rumsey, Christopher L.↗

What Information Theory Says About Best Response and About Binding Contracts

Product Distribution (PD) theory is the information-theoretic extension of conventional full- rationality game theory to bounded rational games. Here PD theory is used to investigate games in which the players use bounded rational best-response strategies. This investigation illuminates how to determine the optimal organization chart for a corporation, or more generally how to order the sequence of moves of the players / employees so as to optimize an overall objective function. It is then shown that in the continuum-time limit, bounded rational best response games result in a variant of the replicator dynamics of evolutionary game theory. This variant is then investigated for team games, in which the players share the same utility function, by showing that such continuum- limit bounded rational best response is identical to Newton-Raphson iterative optimization of the shared utility function. Next PD theory is used to investigate changing the coordinate system of the game, i.e., changing the mapping from the joint move of the players to the arguments in the utility functions. Such a change couples those arguments, essentially by making each players move be an offered binding contract.

Wolpert, David H.↗

Development and Implementation of a Generic Analysis Template for Structural-Thermal-Optical-Performance Modeling

Performance-related effects of system level temperature changes can be a key consideration in the design of many types of optical instruments. This is especially true for space-based imagers, which may require complex thermal control systems to maintain alignment of the optical components. Structural-Thermal-Optical-Performance (STOP) analysis is a multi-disciplinary process that can be used to assess the performance of these optical systems when subjected to the expected design environment. This type of analysis can be very time consuming, which makes it difficult to use as a trade study tool early in the project life cycle. In many cases, only one or two iterations can be performed over the course of a project. This limits the design space to best practices since it may be too difficult, or take too long, to test new concepts analytically. In order to overcome this challenge, automation, and a standard procedure for performing these studies is essential. A methodology was developed within the framework of the Comet software tool that captures the basic inputs, outputs, and processes used in most STOP analyses. This resulted in a generic, reusable analysis template that can be used for design trades for a variety of optical systems. The template captures much of the upfront setup such as meshing, boundary conditions, data transfer, naming conventions, and post-processing, and therefore saves time for each subsequent project. A description of the methodology and the analysis template is presented, and results are described for a simple telescope optical system.

Scola, Salvatore↗

Solution of the symmetric eigenproblem AX=lambda BX by delayed division

Delayed division is an iterative method for solving the linear eigenvalue problem AX = lambda BX for a limited number of small eigenvalues and their corresponding eigenvectors. The distinctive feature of the method is the reduction of the problem to an approximate triangular form by systematically dropping quadratic terms in the eigenvalue lambda. The report describes the pivoting strategy in the reduction and the method for preserving symmetry in submatrices at each reduction step. Along with the approximate triangular reduction, the report extends some techniques used in the method of inverse subspace iteration. Examples are included for problems of varying complexity.

Thurston, G. A.↗

Specifying spacecraft flexible appendage rigidity

As a method for specifying the required degree of rigidity of spacecraft flexible appendages, an analytical technique is proposed for establishing values for the frequency, damping ratio, and modal gain (deflection) of the first several bending modes. The shortcomings of the technique result from the limitations associated with the order of the equations that can be handled practically. An iterative method is prescribed for handling a system whose structural flexibility is described by more than one normal mode. The analytical technique is applied to specifying solar panel rigidity constraints for the NASA Space Telescope. The traditional nonanalytic procedure for specifying the required degree of rigidity of spacecraft flexible appendages has been to set a lower limit below which bending mode frequencies may not lie.

Seltzer, S. M.↗

Inlet Spillage Drag Predictions Using the AIRPLANE Code

AIRPLANE (Jameson/Baker) is a steady inviscid unstructured Euler flow solver. It has been validated on many HSR geometries. It is implemented as MESHPLANE, an unstructured mesh generator, and FLOPLANE, an iterative flow solver. The surface description from an Intergraph CAD system goes into MESHPLANE as collections of polygonal curves to generate the 3D mesh. The flow solver uses a multistage time stepping scheme with residual averaging to approach steady state, but R is not time accurate. The flow solver was ported from Cray to IBM SP2 by Wu-Sun Cheng (IBM); it could only be run on 4 CPUs at a time because of memory limitations. Meshes for the four cases had about 655,000 points in the flow field, about 3.9 million tetrahedra, about 77,500 points on the surface. The flow solver took about 23 wall seconds per iteration when using 4 CPUs. It took about eight and a half wall hours to run 1,300 iterations at a time (the queue limit is 10 hours). A revised version of FLOPLANE (Thomas) was used on up to 64 CPUs to finish up some calculations at the end. We had to turn on more communication when using more processors to eliminate noise that was contaminating the flow field; this added about 50% to the elapsed wall time per iteration when using 64 CPUs. This study involved computing lift and drag for a wing/body/nacelle configuration at Mach 0.9 and 4 degrees pitch. Four cases were considered, corresponding to four nacelle mass flow conditions.

Thomas, Scott D.↗

Numerical simulation of rarefied flow through a slit. I - Direct simulation Monte Carlo results

The pressure-driven flow of a rarefied monatomic gas through a two-dimensional slit is simulated using the direct simulation Monte Carlo technique. Of particular interest is the change in flow field structure as pressure ratio and Knudsen number are varied. Comparisons are made to quantify the limits of validity of free-molecular theory and approximate, nearly free-molecular iterative methods. Also addressed is the sensitivity of the numerical solutions to grid structure and boundary conditions. The free-molecular theory is found to predict quantitative flow field properties (e.g., centerline velocities or downstream flux profiles) reasonably well for large finite Knudsen number with the error dependent on the pressure ratio. The nearly free-molecular corrections are shown to have limited range of applicability. A previously derived parameter is found to correlate total mass flux well as a function of pressure ratio and Knudsen number over a large portion of the transitional regime.

Wadsworth, D. C.↗

A method of limit point calculation in finite element structural analysis

An approach is presented for the calculation of limit points for structures described by discrete coordinates, and whose governing equations derive from finite element concepts. The nonlinear load-displacement path of the imperfect structure is first traced by use of a direct iteration scheme and the determinant of the governing algebraic equations is calculated at each solution point. The limit point is then established by extrapolation and imposition of the condition of zero slope of the plot of load vs. determinant. Three problems are solved in illustration of the approach and in comparison with alternative procedures and test data.

Gallagher, R. H.↗

Low Density Parity Check Codes Based on Finite Geometries: A Rediscovery and More

Low density parity check (LDPC) codes with iterative decoding based on belief propagation achieve astonishing error performance close to Shannon limit. No algebraic or geometric method for constructing these codes has been reported and they are largely generated by computer search. As a result, encoding of long LDPC codes is in general very complex. This paper presents two classes of high rate LDPC codes whose constructions are based on finite Euclidean and projective geometries, respectively. These classes of codes a.re cyclic and have good constraint parameters and minimum distances. Cyclic structure adows the use of linear feedback shift registers for encoding. These finite geometry LDPC codes achieve very good error performance with either soft-decision iterative decoding based on belief propagation or Gallager's hard-decision bit flipping algorithm. These codes can be punctured or extended to obtain other good LDPC codes. A generalization of these codes is also presented.

Kou, Yu↗

A semi-direct procedure using a local relaxation factor and its application to an internal flow problem

Generally, fast direct solvers are not directly applicable to a nonseparable elliptic partial differential equation. This limitation, however, is circumvented by a semi-direct procedure, i.e., an iterative procedure using fast direct solvers. An efficient semi-direct procedure which is easy to implement and applicable to a variety of boundary conditions is presented. The current procedure also possesses other highly desirable properties, i.e.: (1) the convergence rate does not decrease with an increase of grid cell aspect ratio, and (2) the convergence rate is estimated using the coefficients of the partial differential equation being solved.

Chang, S. C.↗

Transient temperature distributions in simple conducting bodies steadily heated through a laminar boundary layer

An analysis is made of the transient heat-conduction effects in three simple semi-infinite bodies: the flat insulated plate, the conical shell, and the slender solid cone. The bodies are assumed to have constant initial temperatures and, at zero time, to begin to move at a constant speed and zero angle of attack through a homogeneous atmosphere. The heat input is taken as that through a laminar boundary layer. Radiation heat transfer and transverse temperature gradients are assumed to be zero. The appropriate heat-conduction equations are solved by an iteration method, the zeroeth-order terms describing the situation in the limit of small time. The method is presented and the solutions are calculated to three orders which are sufficient to give reasonably accurate results when the forward edge has attained one-half the total temperature rise (nose half-rise time). Flight Mach number and air properties occur as parameters in the result. Approximate expressions for the extent of the conduction region and nose half-rise times as functions of the parameters of the problem are presented. (author)

Parker, Hermon M↗

Recovering the fine structures in solar images

Several examples of the capability of the blind iterative deconvolution (BID) technique to recover the real point spread function, when limited a priori information is available about its characteristics. To demonstrate the potential of image post-processing for probing the fine scale and temporal variability of the solar atmosphere, the BID technique is applied to different samples of solar observations from space. The BID technique was originally proposed for correction of the effects of atmospheric turbulence on optical images. The processed images provide a detailed view of the spatial structure of the solar atmosphere at different heights in regions with different large-scale magnetic field structures.

Karovska, Margarita↗

Combined fast multipole-QR compression technique for solving electrically small to large structures for broadband applications

An approach that efficiently solves for a desired parameter of a system or device that can include both electrically large fast multipole method (FMM) elements, and electrically small QR elements. The system or device is setup as an oct-tree structure that can include regions of both the FMM type and the QR type. An iterative solver is then used to determine a first matrix vector product for any electrically large elements, and a second matrix vector product for any electrically small elements that are included in the structure. These matrix vector products for the electrically large elements and the electrically small elements are combined, and a net delta for a combination of the matrix vector products is determined. The iteration continues until a net delta is obtained that is within predefined limits. The matrix vector products that were last obtained are used to solve for the desired parameter.

Jandhyala, Vikram↗

An Evaluation of Size-Resolved Cloud Microphysics Scheme Numerics for use with Radar Observations Part II: Condensation and Evaporation

Accurate numerical modeling of clouds and precipitation is essential for weather forecasting and climate change research. While size-resolved (bin) cloud microphysics models predict particle size distributions without imposing shapes, results are subject to artificial size distribution broadening owing to numerical diffusion associated with various processes. Whereas Part 1 addressed collision-coalescence, here we investigate numerical diffusion that occurs in solving condensation and evaporation. In a parcel model framework, all of the numerical schemes examined converge to one solution of condensation and evaporation as the mass grid is refined, and the advection-based schemes are recommended over the reassigning schemes. Including Eulerian vertical advection in a column limits the convergence to some extent, but that limitation occurs at a sufficiently fine mass grid, and the number of iterations in solving vertical advection should be minimized to reduce numerical diffusion. Insubstantial numerical diffusion in solving condensation can be amplified if collision-coalescence is also active, which in turn can be substantially diminished if turbulence effects on collision are incorporated. Large-eddy simulations of a drizzling stratocumulus field reveal that changes in moments of Doppler spectra obtained using different mass grids are consistent with those obtained from the simplified framework, and that spectral moments obtained using a mass grid designed to effectively reduce numerical diffusion are generally closer to observations. Notable differences between the simulations and observations still exist, and our results suggest a need to consider whether factors other than numerical diffusion in the fundamental process schemes employed can cause such differences.

Cloud microphysics↗

Data System Implications Derived from User Application Requirements for Satellite Data

An investigation of the data system needs as driven by users of space acquired Earth observation data is documented. Two major categories of users, operational and research, are identified. Limiting data acquisition alleviates some of the delays in processing thus improving timeliness of the delivered product. Trade offs occur between timeliness and data distribution costs, and between data storage and reprocessing. The complexity of the data system requirements to apply space data to users' needs is such that no single analysis suffices to design and implement the optimum system. A series of iterations is required with analyses of the salient problems in a general way, followed by a limited implementation of benefit to some users with a continual upgrade in system capacity, functions, and applications served. The resulting most important requirement for the data system is flexibility to accommodate changing requirements as the system is implemented.

Neiers, J.↗