Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence analysis”

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 487 records · Page 27

Vertical profiles of heating derived from IR-based precipitation estimates during FGGE SOP-1

This paper examines a technique for retrieving from geostationary IR data the vertical profiles of heating and cooling due to moist diabatic processes. First, GOES IR imagery is used to estimate precipitation fields which are independent of fields inferred from residuals in heat budget analysis based on the FGGE level III-b data. Vertical distributions of the associated heating are then obtained using thermodynamic data from the level III-b analysis, one-dimensional cloud models, and the satellite-estimated precipitation. The technique was applied to infer heating in the South Pacific convergence zone during a portion of FGEE SOP-1, and the results were compared with heat-budget calculations made using the ECMWF analyses.

Robertson, Franklin R.↗

Stability analysis of multigrid acceleration methods for the solution of partial differential equations

A calculation is made of the stability of various relaxation schemes for the numerical solution of partial differential equations. A multigrid acceleration method is introduced, and its effects on stability are explored. A detailed stability analysis of a simple case is carried out and verified by numerical experiment. It is shown that the use of multigrids can speed convergence by several orders of magnitude without adversely affecting stability.

Fay, John F.↗

Flutter and Divergence Analysis using the Generalized Aeroelastic Analysis Method

The Generalized Aeroelastic Analysis Method (GAAM) is applied to the analysis of three well-studied checkcases: restrained and unrestrained airfoil models, and a wing model. An eigenvalue iteration procedure is used for converging upon roots of the complex stability matrix. For the airfoil models, exact root loci are given which clearly illustrate the nature of the flutter and divergence instabilities. The singularities involved are enumerated, including an additional pole at the origin for the unrestrained airfoil case and the emergence of an additional pole on the positive real axis at the divergence speed for the restrained airfoil case. Inconsistencies and differences among published aeroelastic root loci and the new, exact results are discussed and resolved. The generalization of a Doublet Lattice Method computer code is described and the code is applied to the calculation of root loci for the wing model for incompressible and for subsonic flow conditions. The error introduced in the reduction of the singular integral equation underlying the unsteady lifting surface theory to a linear algebraic equation is discussed. Acknowledging this inherent error, the solutions of the algebraic equation by GAAM are termed 'exact.' The singularities of the problem are discussed and exponential series approximations used in the evaluation of the kernel function shown to introduce a dense collection of poles and zeroes on the negative real axis. Again, inconsistencies and differences among published aeroelastic root loci and the new 'exact' results are discussed and resolved. In all cases, aeroelastic flutter and divergence speeds and frequencies are in good agreement with published results. The GAAM solution procedure allows complete control over Mach number, velocity, density, and complex frequency. Thus all points on the computed root loci can be matched-point, consistent solutions without recourse to complex mode tracking logic or dataset interpolation, as in the k and p-k solution methods.

Edwards, John W.↗

Real time on-chip sequential adaptive principal component analysis for data feature extraction and image compression

In this paper, we present a new, simple, and optimized hardware architecture sequential learning technique for adaptive Principle Component Analysis (PCA) which will help optimize the hardware implementation in VLSI and to overcome the difficulties of the traditional gradient descent in learning convergence and hardware implementation.

PCA feature extraction image compression real time↗

Finite element analysis of transonic flows in highly staggered cascades

Finite element method is employed in the solution of full-potential equation for analyzing steady, transonic flows through a series of highly staggered cascades. Convergence characteristics and efficiency of the developed numerical procedure are discussed. Accuracy of the mathematical model is investigated in terms of the employed computational grid. In particular, the importance of accurate representation of the trailing edge together with the shock on the obtained flow results are studied. Numerical results include choked flows, for a cascade of NACA 0012 airfoils with pitch/chord length = 1, and stagger angle of 45 deg.

Akay, H. U.↗

Flow in a porous nozzle with massive wall injection

An analytical and experimental investigation has been conducted to determine the effect of massive wall injection on the flow characteristics in a nozzle. The experiments were performed on a water table with a porous-nozzle test section. This had 45 deg and 15 deg half angles of convergence and divergence, respectively, throat radius of 2.5 inches, and throat width of 3 inches. The hydraulic analogy was employed to qualitatively extend the results to a compressible gas flow through the nozzle. An analysis of the water table flow was made using a one-dimensional flow assumption in the continuity and momentum equations. An analysis of a compressible flow in a nozzle was made in a manner analogous to that for the water flow. It is shown that the effect of blowing is to move the sonic position downstream of the geometric throat. Similar results were determined for the incompressible water table flow. Limited photographic results are presented for an injection of air, CO2, and Freon-12 into a main-stream air flow in a convergent-divergent nozzle. Schlieren photographs were used to visualize the flow.

Kinney, R. B.↗

Analysis of a Shock-Associated Noise Prediction Model Using Measured Jet Far-Field Noise Data

A code for predicting supersonic jet broadband shock-associated noise was assessed us- ing a database containing noise measurements of a jet issuing from a convergent nozzle. The jet was operated at 24 conditions covering six fully expanded Mach numbers with four total temperature ratios. To enable comparisons of the predicted shock-associated noise component spectra with data, the measured total jet noise spectra were separated into mixing noise and shock-associated noise component spectra. Comparisons between predicted and measured shock-associated noise component spectra were used to identify de ciencies in the prediction model. Proposed revisions to the model, based on a study of the overall sound pressure levels for the shock-associated noise component of the mea- sured data, a sensitivity analysis of the model parameters with emphasis on the de nition of the convection velocity parameter, and a least-squares t of the predicted to the mea- sured shock-associated noise component spectra, resulted in a new de nition for the source strength spectrum in the model. An error analysis showed that the average error in the predicted spectra was reduced by as much as 3.5 dB for the revised model relative to the average error for the original model.

Dahl, Milo D.↗

A Downscaling Analysis of the Urban Influence on Rainfall: TRMM Satellite Component AMS Conference on Satellite Meteorology and Oceanography

A recent publication by Shepherd et al. (2002) demonstrated the feasibility of using TRMM precipitation radar (PR) estimates to identify precipitation anomalies caused by urbanization. The approach is particularly useful for investigating this global process because TRMM data span large portions of the globe and comprise an extended temporal dataset. Recent literature suggests that urbanized regions of Houston, Texas may be influencing lightning and precipitation formation over and downwind of the city. Possible mechanisms include: (1) enhanced convergence through interactions between the sea breeze, Galveston bay breeze, and urban heat island circulations, (2) enhanced convergence due to increased surface roughness over the city and/or destabilization of the boundary layer by the UHI, or (3) enhanced cloud condensation nuclei due to urban and industrial aerosol sources. In this study, a downscaling analysis of spatial and temporal trends in rainfall around the Houston Area is being conducted. The downscaling analysis concept involves identifying and quantifying urban rainfall anomalies at progressively smaller spatial and temporal scales using the TRMM satellite, ground-based radar, and a dense network of rain gauges. The goal is to test the hypothesis that the Houston urban district and regions in the climatological downwind region of the city exhibit enhanced rainfall amounts relative to the climatological upwind regions. TRMM was launched in 1997 and currently operates in a low-inclination (35 deg), non-sun-synchronous orbit at an altitude of 402 km (350 km prior to August 2001). The satellite analysis follows the methodologies described in Shepherd et al. (2002). Nearly five years of TRMM PR-derived mean monthly rainfall estimates are utilized to produce annual and warm season isohyetal analyses around Houston. Early results indicate that rainfall rates (mm/h) for the entire period are largest within 100 km northeast and east of Houston (e.g. the "hypothesized downwind region"). The mean rainfall rate over the Houston urban center is 30.5% larger than the upwind control region. The mean rainfall rate in the downwind region is 34.4% larger than the upwind region. An analysis of a parameter called the urban rainfall ratio (URR) illustrates that 65% (88%) of the satellite-derived rainfall rates in the downwind (upwind control) region are greater (less) than the mean background rainfall rate of the entire study region. When the data is stratified by summer months from 1998 to 2001 (June-August), even greater influence over and downwind of the urban area is observed in the statistics. This result is consistent with published reports of urban-generated rainfall being more prevalent in the warm season. The research demonstrates that the evolving TRMM satellite climatology is a credible way to detect mesoscale precipitation signatures that may be linked to urbanization. Early results also corroborate recent findings on Houston-induced convection/drainfall anomalies. Burian and Shepherd will report on other aspects of the downscaling analysis in future forums, but early rain gauge results are consistent with the satellite-based observations.

Shepherd, J. Marshall↗

The Necessity of Functional Analysis for Space Exploration Programs

As NASA moves toward expanded commercial spaceflight within its human exploration capability, there is increased emphasis on how to allocate responsibilities between government and commercial organizations to achieve coordinated program objectives. The practice of program-level functional analysis offers an opportunity for improved understanding of collaborative functions among heterogeneous partners. Functional analysis is contrasted with the physical analysis more commonly done at the program level, and is shown to provide theoretical performance, risk, and safety advantages beneficial to a government-commercial partnership. Performance advantages include faster convergence to acceptable system solutions; discovery of superior solutions with higher commonality, greater simplicity and greater parallelism by substituting functional for physical redundancy to achieve robustness and safety goals; and greater organizational cohesion around program objectives. Risk advantages include avoidance of rework by revelation of some kinds of architectural and contractual mismatches before systems are specified, designed, constructed, or integrated; avoidance of cost and schedule growth by more complete and precise specifications of cost and schedule estimates; and higher likelihood of successful integration on the first try. Safety advantages include effective delineation of must-work and must-not-work functions for integrated hazard analysis, the ability to formally demonstrate completeness of safety analyses, and provably correct logic for certification of flight readiness. The key mechanism for realizing these benefits is the development of an inter-functional architecture at the program level, which reveals relationships between top-level system requirements that would otherwise be invisible using only a physical architecture. This paper describes the advantages and pitfalls of functional analysis as a means of coordinating the actions of large heterogeneous organizations for space exploration programs.

program management↗

Iterative Discrete Ordinates Solution of the Equation for the Surface-Reflected Radiance

This paper presents a new method of numerical solution of the integral equation for the radiance reflected from an anisotropic surface. The equation relates the radiance at the surface level with BRDF and solutions of the standard radiative transfer problems for a slab with no reflection on its surfaces. It is also shown that the kernel of the equation satisfies the condition of the existence of a unique solution and the convergence of the successive approximations to that solution. The developed method features two basic steps: discretization on a 2D quadrature, and solving the resulting system of algebraic equations with successive over-relaxation method based on the Gauss-Seidel iterative process. Presented numerical examples show good coincidence between the surface-reflected radiance obtained with DISORT and the proposed method. Analysis of contributions of the direct and diffuse (but not yet reflected) parts of the downward radiance to the total solution is performed. Together, they represent a very good initial guess for the iterative process. This fact ensures fast convergence. The numerical evidence is given that the fastest convergence occurs with the relaxation parameter of 1 (no relaxation). An integral equation for BRDF is derived as inversion of the original equation. The potential of this new equation for BRDF retrievals is analyzed. The approach is found not viable as the BRDF equation appears to be an ill-posed problem, and it requires knowledge the surface-reflected radiance on the entire domain of both Sun and viewing zenith angles.

Alexander Radkevich↗

Analysis of a parallel multigrid algorithm

This paper considers the parallel multigrid algorithm of Frederickson and McBryan (1987). This algorithm uses multiple coarse-grid problems (instead of one) in the hope of accelerating convergence, and is found to have a close relationship to traditional multigrid methods. Specifically, the parallel coarse-grid correction operator is identical to a traditional multigrid coarse-grid correction operator, except that the mixing of high and low frequencies caused by aliasing error is removed. Appropriate relaxation operators can be chosen to take advantage of this property. Comparisons between the standard multigrid and the new method are made.

Chan, Tony F.↗

Analysis of a parallel multigrid algorithm

The parallel multigrid algorithm of Frederickson and McBryan (1987) is considered. This algorithm uses multiple coarse-grid problems (instead of one problem) in the hope of accelerating convergence and is found to have a close relationship to traditional multigrid methods. Specifically, the parallel coarse-grid correction operator is identical to a traditional multigrid coarse-grid correction operator, except that the mixing of high and low frequencies caused by aliasing error is removed. Appropriate relaxation operators can be chosen to take advantage of this property. Comparisons between the standard multigrid and the new method are made.

Chan, Tony F.↗

Analysis of a Decision-Directed Receiver with Unknown Priors

A decision-directed receiver (DDR) uses previous outputs (decisions) to estimate unknown parameters and, on the basis of these estimates, modifies the detector structure for subsequent decisions. Although the DDR is less complex to instrument than other adaptive schemes, inherent in the decision-directed approach is the possibility of a runaway. This occurs when the detector commits a sequence of decision errors resulting in a degradation of parameter estimates, which, in turn, results in a further deterioration of detector performance. Because of the dependencies introduced by the learning process, runaway is difficult to analyze. In this paper, a DDR with unknown a priori probabilities is considered. The priors are estimated by the relative frequency of decisions of that event. For binary detection, it is shown that there is a positive probability of a runaway (the estimates converge to 1 or 0), which equals 1 if the signal-to-noise ratio is below a critical value. A tight bound on the probability of a runaway is obtained by approximating the learning process by a random walk with independent increments. The analysis demonstrates that a runaway is quite improbable even for moderate signal-to-noise ratios. The analysis is extended to multiple signals and to the situation where the estimates of the priors are updated continuously through exponential weighting rather than allowed to converge.

Davisson, Lee D.↗

Vibrations of cantilevered shallow cylindrical shells of rectangular planform

A cantilevered, shallow shell of circular cylindrical curvature and rectangular planform exhibits free vibration behavior which differs considerably from that of a cantilevered beam or of a flat plate. Some numerical results can be found for the problem in the previously published literature, mainly obtained by using various finite element methods. The present paper is the first definitive study of the problem, presenting accurate non-dimensional frequency parameters for wide ranges of aspect ratio, shallowness ratio and thickness ratio. The analysis is based upon shallow shell theory. Numerical results are obtained by using the Ritz method, with algebraic polynomial trial functions for the displacements. Convergence is investigated, with attention being given both to the number of terms taken for each co-ordinate direction and for each of the three components of displacement. Accuracy of the results is also established by comparison with finite element results for shallow shells and with other accurate flat plate solutions.

Leissa, A. W.↗

A thermal NO(x) prediction model - Scalar computation module for CFD codes with fluid and kinetic effects

A thermal NO(x) prediction model is developed to interface with a CFD, k-epsilon based code. A converged solution from the CFD code is the input to the postprocessing model for prediction of thermal NO(x). The model uses a decoupled analysis to estimate the equilibrium level of (NO(x))e which is the constant rate limit. This value is used to estimate the flame (NO(x)) and in turn predict the rate of formation at each node using a two-step Zeldovich mechanism. The rate is fixed on the NO(x) production rate plot by estimating the time to reach equilibrium by a differential analysis based on the reaction: O + N2 = NO + N. The rate is integrated in the nonequilibrium time space based on the residence time at each node in the computational domain. The sum of all nodal predictions yields the total NO(x) level.

Mcbeath, Giorgio↗

Steady and unsteady transonic small disturbance analysis of realistic aircraft configurations

A transonic unsteady aerodynamic and aeroelasticity code CAP-TSD (Computational Aeroelasticity Program-Transonic Small Disturbance) has been developed for application to realistic aircraft configurations. The code uses a time-accurate approximate factorization (AF) algorithm for solution of the unsteady transonic small-disturbance potential equation. The paper gives an overview of the CAP-TSD code development effort and reports on recent algorithm modifications. The algorithm modifications include: an Engquist-Osher (E-O) type-dependent switch to treat regions of supersonic flow, extension of the E-O switch for second-order spatial accuracy, nonisentropic effects to treat strong-shock cases, nonreflecting far field boundary conditions for unsteady applications, and several modifications to accelerate convergence to steady state. The modifications have significantly improved the stability of the AF algorithm and hence the reliability of the CAP-TSD code in general. Calculations are also presented from a flutter analysis of a 45-deg sweptback wing which agree well with experimental data. The paper present descriptions of the CAP-TSD code and algorithm details along with results.

Batina, John T.↗

New Flutter Analysis Technique for Time-Domain Computational Aeroelasticity

A new time-domain approach for computing flutter speed is presented. Based on the time-history result of aeroelastic simulation, the unknown unsteady aerodynamics model is estimated using a system identification technique. The full aeroelastic model is generated via coupling the estimated unsteady aerodynamic model with the known linear structure model. The critical dynamic pressure is computed and used in the subsequent simulation until the convergence of the critical dynamic pressure is achieved. The proposed method is applied to a benchmark cantilevered rectangular wing.

HWB aircraft↗

Origin of longitudinal grooving in Tiu Vallis, Mars - Isolation of responsible fluid-types

Fluid instability analysis developed for the growth or decay of large, longitudinal roller vortex perturbations imposed at the base of fluids of exponentially stratified viscosities is presented. The origin of longitudinal grooving and fluting in a region of converging flow in Tiu Vallis, Mars, is analyzed, and the propensity of fluids, such as water, mudflow/debris flow, and lava, to develop longitudinal roll patterns as part of their flow field is tested. It is found that the grooving pattern in Tiu Vallis was probably formed by vortex motion in a mudflow or debris-laden slurry, but exact representation of sediment/debris concentration or sediment transport machanism is not yet possible. It is concluded that longitudinal vortex motion modelled in this analysis can develop in water, lava, or ice, but that the specific flow characteristics of these models which include the vortex motion are unreasonable for flow in Tiu Vallis.

Thompson, D. E.↗