Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “successive linear approximation”

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 55 records · Page 3

Implicit solvers for unstructured meshes

Implicit methods were developed and tested for unstructured mesh computations. The approximate system which arises from the Newton linearization of the nonlinear evolution operator is solved by using the preconditioned GMRES (Generalized Minimum Residual) technique. Three different preconditioners were studied, namely, the incomplete LU factorization (ILU), block diagonal factorization, and the symmetric successive over relaxation (SSOR). The preconditioners were optimized to have good vectorization properties. SSOR and ILU were also studied as iterative schemes. The various methods are compared over a wide range of problems. Ordering of the unknowns, which affects the convergence of these sparse matrix iterative methods, is also studied. Results are presented for inviscid and turbulent viscous calculations on single and multielement airfoil configurations using globally and adaptively generated meshes.

Venkatakrishnan, V.↗

Relaxation methods in fluid mechanics

The present work considers the iterative solution of a coupled set of difference equations and examines methods that carry successive approximates to a state that is invariant with further iteration and independent of the initial guess. Methods are studied with regard to their efficiency and economy of computer resources. The basic principles of classical relaxation are set forth, with attention confined to linear elliptic equations. This discussion involves the evaluation of the spectral radius that is the magnitude of the eigenvalue with largest modulus. The subject of relaxation is then related to the study of ordinary differential equations and hyperbolic partial differential equations. Problems that occur when linearly dependent eigenvectors appear in the relaxation matrix are discussed, leading to multiply connected eigenvalues in the Jordan canonical form. Finally, a brief survey of relaxation methods used in aerodynamics is given.

Lomax, H.↗

Changes in the performance characteristics of a GaAs near infrared light emitting diode when exposed to various current and thermal stresses

The changes that occurred in the optical and electrical characteristics of a near infrared, GaAs light emitting diode, when operated under various levels and combinations of current and thermal stresses are discussed. A total of forty parts were operated for two thousand hours under eight different sets of dc current and ambient temperature conditions. Degradation in the radiant optical power of these devices was thirty-four percent when operated at their rated current and an ambient temperature of 298K (25 C). Derating the current and/or the thermal stress reduced the degradation of this parameter in approximately a linear manner. All degraded devices behaved similarly, exhibiting rapid nonlinear degradation followed by a gradual linear degradation and finally a period of stable operation. An attempt was made to correlate initial device condition to degradation during stress testing, but met with little success.

Thomas, E. F., Jr.↗

Sensitivity derivatives for three dimensional supersonic Euler code using incremental iterative strategy

In a recent work, an incremental strategy was proposed to iteratively solve the very large systems of linear equations that are required to obtain quasianalytical sensitivity derivatives from advanced computational fluid dynamics (CFD) codes. The technique was sucessfully demonstrated for two large two-dimensional problems: a subsonic and a transonic airfoil. The principal feature of this incremental iterative stategy is that it allows the use of the identical approximate coefficient matrix operator and algorithm to solve the nonlinear flow and the linear sensitivity equations; at convergence, the accuracy of the sensitivity derivatives is not compromised. This feature allows a comparatively straightforward extension of the methodology to three-dimensional problems; this extension is successfully demonstrated in the present study for a space-marching solution of the three-dimensional Euler equations over a Mach 2.4 blended wing-body configuration.

Korivi, Vamshi Mohan↗

Projection methods for the numerical solution of Markov chain models

Projection methods for computing stationary probability distributions for Markov chain models are presented. A general projection method is a method which seeks an approximation from a subspace of small dimension to the original problem. Thus, the original matrix problem of size N is approximated by one of dimension m, typically much smaller than N. A particularly successful class of methods based on this principle is that of Krylov subspace methods which utilize subspaces of the form span(v,av,...,A(exp m-1)v). These methods are effective in solving linear systems and eigenvalue problems (Lanczos, Arnoldi,...) as well as nonlinear equations. They can be combined with more traditional iterative methods such as successive overrelaxation, symmetric successive overrelaxation, or with incomplete factorization methods to enhance convergence.

Saad, Youcef↗

Cosmic ray propagation in interplanetary space

The validity of the test-particle picture, the approximation of static fields, and the spatial-diffusion approximation are discussed in a general way before specific technical assumptions are introduced. It is argued that the spatial-diffusion equation for the intensity per unit energy has a much wider range of applicability than the kinetic (Fokker-Planck) equation it is derived from. This gives strong weight to the phenomenological propagation theory. The general success (and possible failure at small energies) of the phenomenological theory for the modulation of galactic cosmic rays and solar events is described. Apparent effects such as the 'free boundary' are given disproportionate weight since they establish the connection with the detailed plasma physics of the solar wind. Greatest attention is paid to the pitch-angle diffusion theory. A general theory is presented which removes the well-known secularities of the quasi-linear approximation. The possible breakdown of any pitch-angle diffusion theory at very small energies is perhaps connected with the observed 'turn up' of the spectrum at low energies. A first attempt to derive the spatial dependence of the diffusion coefficient in the solar cavity, using such a divergence free scattering theory, is described and compared with recent observations out to 5 AU.

Voelk, H. J.↗

COLUTA: Custom 8-Channel 15-bit 40-MSPS ADC for the ATLAS Liquid Argon Calorimeter Readout

The stringent radiation-tolerance, quantization requirements, and the need for seamless integration of the on-detector readout electronics chain require the design of a full-custom analog-to-digital converter (ADC) for the 182,468 channels in the ATLAS Liquid Argon (LAr) calorimeter at the Large Hadron Collider (LHC) at CERN, as part of the High-Luminosity LHC (HL-LHC) upgrade. Each of the 8 channels in the prototype 65 nm CMOS ADC with 15-bit resolution and >68 dB signal-to-noise-and-distortion ratio (SNDR), or equivalently >11 effective number of bits (ENOB), consists of a Multiplying-DAC (MDAC) and a successive-approximation (SAR) ADC. A digital data processing unit (DDPU) calibrates and serially transmits the quantized data. Initial lab characterization of 18 chips shows a minimum SNDR of 69.5 dB at full-scale at about 5 MHz and a maximum power consumption of 1.17 W; differential non-linearity (DNL) measurements show no missing codes. Additional validation of the ADC is ongoing before approximately 70k chips can be mass-produced for the upgrade.

Rui Xu↗

Recursive algorithms for two-dimensional smoothing using bicubic hermite polynomial

It is noted that in the past, smoothing splines originated from approximation theory have been successfully applied to data filtering and image smoothing problems. Even though the nonrecursive technique of smoothing splines gives an optimal solution, the amount of computation increases rapidly with the size of the two-dimensional data. A derivation is presented here of quarter-plane filtering algorithms that provide smoothed estimates of function values and their derivatives by fitting two-dimensional smoothing splines in a recursive manner. The derivation procedure sheds light on specific problems encountered in two-dimensional filtering problems. What is more, the amount of computation for this recursive processor increases only linearly with the size of the two-dimensional data. Because of certain approximations introduced in its derivation, this recursive processor becomes suboptimal.

Kim, C. S.↗

Development and testing of stable, invariant, isoparametric curvilinear 2- and 3-D hybrid-stress elements

Linear and quadratic Serendipity hybrid-stress elements are examined in respect of stability, coordinate invariance, and optimality. A formulation based upon symmetry group theory successfully addresses these issues in undistorted geometries and is fully detailed for plane elements. The resulting least-order stable invariant stress polynomials can be applied as astute approximations in distorted cases through a variety of tensor components and variational principles. A distortion sensitivity study for two- and three-dimensional elements provides favorable numerical comparisons with the assumed displacement method.

Punch, E. F.↗

Integrated Stirling Convertor and Hall Thruster Test Conducted

An important aspect of implementing Stirling Radioisotope Generators on future NASA missions is the integration of the generator and controller with potential spacecraft loads. Some recent studies have indicated that the combination of Stirling Radioisotope Generators and electric propulsion devices offer significant trip time and payload fraction benefits for deep space missions. A test was devised to begin to understand the interactions between Stirling generators and electric thrusters. An electrically heated RG- 350 (350-W output) Stirling convertor, designed and built by Stirling Technology Company of Kennewick, Washington, under a NASA Small Business Innovation Research agreement, was coupled to a 300-W SPT-50 Hall-effect thruster built for NASA by the Moscow Aviation Institute (RIAME). The RG-350 and the SPT-50 shown, were installed in adjacent vacuum chamber ports at NASA Glenn Research Center's Electric Propulsion Laboratory, Vacuum Facility 8. The Stirling electrical controller interfaced directly with the Hall thruster power-processing unit, both of which were located outside of the vacuum chamber. The power-processing unit accepted the 48 Vdc output from the Stirling controller and distributed the power to all the loads of the SPT-50, including the magnets, keeper, heater, and discharge. On February 28, 2001, the Glenn test team successfully operated the Hall-effect thruster with the Stirling convertor. This is the world's first known test of a dynamic power source with electric propulsion. The RG-350 successfully managed the transition from the purely resistive load bank within the Stirling controller to the highly capacitive power-processing unit load. At the time of the demonstration, the Stirling convertor was operating at a hot temperature of 530 C and a cold temperature of -6 C. The linear alternator was producing approximately 250 W at 109 Vac, while the power-processing unit was drawing 175 W at 48 Vdc. The majority of power was delivered to the Hall thruster discharge circuit operating at 115 Vdc and 0.9 A. Testing planned for late 2001 will examine the possibility of directly driving the Hall thruster discharge circuit using rectified and filtered output from the Stirling alternator.

Mason, Lee S.↗

Delay Line Detectors for the UVCS and Sumer Instruments on the SOHO Satellite

Microchannel plate based detectors with cross delay line image readout have been rapidly implemented for the SUMER and UVCS instruments aboard the Solar Orbiting Heliospheric Observatory (SOHO) mission to be launched in July 1995. In October 1993 a fast track program to build and characterize detectors and detector control electronics was initiated. We present the detector system design for the SOHO UVCS and SUMER detector programs, and results from the detector test program. Two deliverable detectors have been built at this point, a demonstration model for UVCS, and the flight Ly alpha detector for UVCS, both of which are to be delivered in the next few weeks. Test results have also been obtained with one other demonstration detector system. The detector format is 26mm x 9mm, with 1024 x 360 digitized pixels, using a low resistance Z stack of microchannel plates (MCP's) and a multilayer cross delay line anode (XDL). This configuration provides gains of approximately 2 x 10(exp 7) with good pulse height distributions (less than 50% FWHM) under uniform flood illumination, and background levels typical for this configuration (approximately 0.6 event cm (exp -2)sec(exp -1)). Local counting rates up to about 400 events/pixel/sec have been achieved with no degradation of the MCP gain. The detector and event encoding electronics achieves about 25 millimeter FVHM with good linearity (plus or minus approximately 1 pixel) and is stable to high global counting rates (greater than 4 x 10(exp 5) events sec(exp -1)). Flat field images are dominated by MCP fixed pattern noise and are stable, but the MCP multifiber modulation usually expected is uncharacteristically absent. The detector and electronics have also successfully passed both thermal vacuum and vibration tests.

Seigmund, O. H. W.↗

Nondestructive Determination of Bond Strength

Although many nondestructive techniques have been applied to detect disbonds in adhesive joints, no absolutely reliable nondestructive method has been developed to detect poor adhesion and evaluate the strength of bonded joints prior to the present work which used nonlinear ultrasonic methods to investigate adhesive bond cure conditions. Previously, a variety of linear and nonlinear ultrasonic methods with water coupling had been used to study aluminum-adhesive-aluminum laminates, prepared under different adhesive curing conditions, for possible bond strength determination. Therefore, in the course of this research effort, a variety of finite-amplitude experimental methods which could possibly differentiate various cure conditions were investigated, including normal and oblique incidence approaches based on nonlinear harmonic generation as well as several non-collinear two-wave interaction approaches. Test samples were mechanically scanned in various ways with respect to the focus of a transmitting transducer operated at several variable excitation frequencies and excitation levels. Even when powerful sample-related resonances were exploited by means of a frequency scanning approach, it was very difficult to isolate the nonlinear characteristics of adhesive bonds. However, a multi-frequency multi-power approach was quite successful and reliable. Ultrasonic tone burst signals at increasing power levels, over a wide frequency range, were transmitted through each bond specimen to determine its excitation dependent nonlinear harmonic resonance behavior. Relative amplitude changes were observed particularly in the higher harmonic spectral data and analyzed using a local displacement and strain analysis in the linear approximation. Two analysis approaches of the excitation-dependent data at specific resonances were found to be quite promising. One of these approaches may represent a very robust algorithm for classifying an adhesive bond as being properly cured or not. Another approach, in addition to differentiation between various cure conditions, may even provide information with respect to the bond strength. Several technical papers were published during the course of this research and a summary is presented in the Ph.D. dissertation of Tobias P. Berndt, a graduate student financially supported by this NASA Grant.

Source record↗

Optimizing the Zeldovich approximation

We have recently learned that the Zeldovich approximation can be successfully used for a far wider range of gravitational instability scenarios than formerly proposed; we study here how to extend this range. In previous work (Coles, Melott and Shandarin 1993, hereafter CMS) we studied the accuracy of several analytic approximations to gravitational clustering in the mildly nonlinear regime. We found that what we called the 'truncated Zeldovich approximation' (TZA) was better than any other (except in one case the ordinary Zeldovich approximation) over a wide range from linear to mildly nonlinear (sigma approximately 3) regimes. TZA was specified by setting Fourier amplitudes equal to zero for all wavenumbers greater than k(sub nl), where k(sub nl) marks the transition to the nonlinear regime. Here, we study the cross correlation of generalized TZA with a group of n-body simulations for three shapes of window function: sharp k-truncation (as in CMS), a tophat in coordinate space, or a Gaussian. We also study the variation in the crosscorrelation as a function of initial truncation scale within each type. We find that k-truncation, which was so much better than other things tried in CMS, is the worst of these three window shapes. We find that a Gaussian window e(exp(-k(exp 2)/2k(exp 2, sub G))) applied to the initial Fourier amplitudes is the best choice. It produces a greatly improved crosscorrelation in those cases which most needed improvement, e.g. those with more small-scale power in the initial conditions. The optimum choice of kG for the Gaussian window is (a somewhat spectrum-dependent) 1 to 1.5 times k(sub nl). Although all three windows produce similar power spectra and density distribution functions after application of the Zeldovich approximation, the agreement of the phases of the Fourier components with the n-body simulation is better for the Gaussian window. We therefore ascribe the success of the best-choice Gaussian window to its superior treatment of phases in the nonlinear regime. We also report on the accuracy of particle positions and velocities produced by TZA.

Melott, Adrian L.↗

Air fluorescence detection of large air showers below the horizon

In the interest of exploring the cosmic ray spectrum at energies greater than 10 to the 18th power eV, where flux rates at the Earth's surface drop below 100 yr(-1) km(-2) sr(-1), cosmic ray physicists have been forced to construct ever larger detectors in order to collect useful amounts of data in reasonable lengths of time. At present, the ultimate example of this trend is the Fly's Eye system in Utah, which uses the atmosphere around an array of skyward-looking photomultiplier tubes. The air acts as a scintillator to give detecting areas as large as 5000 square kilometers sr (for highest energy events). This experiment has revealed structure (and a possible cutoff) in the ultra-high energy region above 10 o the 19th power eV. The success of the Fly's Eye experiment provides impetus for continuing the development of larger detectors to make accessible even higher energies. However, due to the rapidly falling flux, a tenfold increase in observable energy would call for a hundredfold increase in the detecting area. But, the cost of expanding the Fly's Eye detecting area will approximately scale linearly with area. It is for these reasons that the authors have proposed a new approach to using the atmosphere as a scintillator; one which will require fewer photomultipliers, less hardware (thus being less extensive), yet will provide position and shower size information.

Halverson, P.↗

Numerical simulation of boundary-layer transition

The transition to turbulence in boundary layers was investigated by direct numerical solution of the nonlinear, three-dimensional, incompressible Navier-Stokes equations in the half-infinite domain over a flat plate. Periodicity was imposed in the streamwise and spanwise directions. A body force was applied to approximate the effect of a nonparallel mean flow. The numerical method was spectra, based on Fourier series and Jacobi polynomials, and used divergence-free basis functions. Extremely rapid convergence was obtained when solving the linear Orr-Sommerfeld equation. The early nonlinear and three-dimensional stages of transition, in a boundary layer disturbed by a vibrating ribbon, were successfully simulated. Excellent qualitative agreement was observed with either experiments or weakly nonlinear theories. In particular, the breakdown pattern was staggered or nonstaggered depending on the disturbance amplitude.

Spalart, P. R.↗

Generative Representations for Automated Design of Robots

A method of automated design of complex, modular robots involves an evolutionary process in which generative representations of designs are used. The term generative representations as used here signifies, loosely, representations that consist of or include algorithms, computer programs, and the like, wherein encoded designs can reuse elements of their encoding and thereby evolve toward greater complexity. Automated design of robots through synthetic evolutionary processes has already been demonstrated, but it is not clear whether genetically inspired search algorithms can yield designs that are sufficiently complex for practical engineering. The ultimate success of such algorithms as tools for automation of design depends on the scaling properties of representations of designs. A nongenerative representation (one in which each element of the encoded design is used at most once in translating to the design) scales linearly with the number of elements. Search algorithms that use nongenerative representations quickly become intractable (search times vary approximately exponentially with numbers of design elements), and thus are not amenable to scaling to complex designs. Generative representations are compact representations and were devised as means to circumvent the above-mentioned fundamental restriction on scalability. In the present method, a robot is defined by a compact programmatic form (its generative representation) and the evolutionary variation takes place on this form. The evolutionary process is an iterative one, wherein each cycle consists of the following steps: 1. Generative representations are generated in an evolutionary subprocess. 2. Each generative representation is a program that, when compiled, produces an assembly procedure. 3. In a computational simulation, a constructor executes an assembly procedure to generate a robot. 4. A physical-simulation program tests the performance of a simulated constructed robot, evaluating the performance according to a fitness criterion to yield a figure of merit that is fed back into the evolutionary subprocess of the next iteration. In comparison with prior approaches to automated evolutionary design of robots, the use of generative representations offers two advantages: First, a generative representation enables the reuse of components in regular and hierarchical ways and thereby serves a systematic means of creating more complex modules out of simpler ones. Second, the evolved generative representation may capture intrinsic properties of the design problem, so that variations in the representations move through the design space more effectively than do equivalent variations in a nongenerative representation. This method has been demonstrated by using it to design some robots that move, variously, by walking, rolling, or sliding. Some of the robots were built (see figure). Although these robots are very simple, in comparison with robots designed by humans, their structures are more regular, modular, hierarchical, and complex than are those of evolved designs of comparable functionality synthesized by use of nongenerative representations.

Homby, Gregory S.↗

Suboptimal schemes for atmospheric data assimilation based on the Kalman filter

This work is directed toward approximating the evolution of forecast error covariances for data assimilation. The performance of different algorithms based on simplification of the standard Kalman filter (KF) is studied. These are suboptimal schemes (SOSs) when compared to the KF, which is optimal for linear problems with known statistics. The SOSs considered here are several versions of optimal interpolation (OI), a scheme for height error variance advection, and a simplified KF in which the full height error covariance is advected. To employ a methodology for exact comparison among these schemes, a linear environment is maintained, in which a beta-plane shallow-water model linearized about a constant zonal flow is chosen for the test-bed dynamics. The results show that constructing dynamically balanced forecast error covariances rather than using conventional geostrophically balanced ones is essential for successful performance of any SOS. A posteriori initialization of SOSs to compensate for model - data imbalance sometimes results in poor performance. Instead, properly constructed dynamically balanced forecast error covariances eliminate the need for initialization. When the SOSs studied here make use of dynamically balanced forecast error covariances, the difference among their performances progresses naturally from conventional OI to the KF. In fact, the results suggest that even modest enhancements of OI, such as including an approximate dynamical equation for height error variances while leaving height error correlation structure homogeneous, go a long way toward achieving the performance of the KF, provided that dynamically balanced cross-covariances are constructed and that model errors are accounted for properly. The results indicate that such enhancements are necessary if unconventional data are to have a positive impact.

Todling, Ricardo↗

Structural Modeling of a Five-Meter Thin Film Inflatable Antenna/Concentrator

Inflatable structures have been the subject of renewed interest in recent years for space applications such as communications antennas, solar thermal propulsion, and space solar power. A major advantage of using inflatable structures in space is their extremely light weight. An obvious second advantage is on-orbit deployability and related space savings in the launch configuration. A recent technology demonstrator flight for inflatable structures was the Inflatable Antenna Experiment (IAE) that was deployed on orbit from the Shuttle Orbiter. Although difficulty was encountered in the inflation/deployment phase, the flight was successful overall and provided valuable experience in the use of such structures. Several papers on static structural analysis of inflated cylinders have been written, describing different techniques such as linear shell theory, and nonlinear and variational methods, but very little work had been done in dynamics of inflatable structures until recent years. In 1988 Leonard indicated that elastic beam bending modes could be utilized in approximating lower-order frequencies of inflatable beams. Main, et al. wrote a very significant 1995 paper describing results of modal tests of inflated cantilever beams and the determination of effective material properties. Changes in material properties for different pressures were also discussed, and the beam model was used in a more complex structure. The paper demonstrated that conventional finite element analysis packages could be very useful in the analysis of complex inflatable structures. The purposes of this paper are to discuss the methodology for dynamically characterizing a large 5-meter thin film inflatable reflector, and to discuss the test arrangement and results. Nonlinear finite element modal results are compared to modal test data. The work is significant and of considerable interest to researchers because of 1) the large size of the structure, making it useful for scaling studies, and 2) application of commercially available finite element software for modeling pressurized thin-film structures.

Smalley, Kurt B.↗