Search NASA⌕ Search

SEARCH · Search NASA

Results for “iterative method”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

Propulsion-Airframe Integration for Conceptual Redesign of a Low-Boom Supersonic Transport

A low-boom supersonic transport was generated in a previous low-boom multidisciplinary optimization (MDO) study that used computational fluid dynamics (CFD) off-body pressure to compute the undertrack sonic boom ground signature and calibrated low-fidelity aerodynamics analyses to compute the mission performance metrics. This low-boom aircraft, referred to as the Mach 1.7 40-PAX concept, can carry 40 passengers for a low-boom overland mission with cruise Mach 1.7 for the airport pairs over the continental US and has the potential to achieve an undertrack sonic boom ground noise level below 70 perceived level of decibels at the start of overland cruise (SOC). The engine for the Mach 1.7 40-PAX concept was designed using the Numerical Propulsion System Simulation (NPSS) for the mission analysis and modeled as a flow-through nacelle for CFD-based sonic boom analysis. To understand how the engine plume affects the undertrack ground signature, the Mach 1.7 40-PAX concept is redesigned in this paper after replacing the flow-through nacelles by CFD engines for sonic boom analysis using aeropropulsive CFD simulation. A process is developed for approximation of the NPSS engine at SOC by a CFD engine for aeropropulsive CFD simulation. The generated CFD engine has the identical nozzle boundary conditions and approximately the same mass flow and thrust as those of the NPSS engine at SOC. Then, the outer mold line of the configuration with the CFD engines is redesigned to approximately restore the low-boom characteristics of the Mach 1.7 40-PAX concept at SOC. Finally, the CFD simulation data for the redesigned concept with the CFD engines is used to calibrate the low-fidelity aerodynamic analyses for mission analysis of this concept. The cyclic dependency of the involved disciplinary analyses is resolved using an iteration method for a consistent coupling of the mission analysis, NPSS engine analysis, and low-boom redesign using the aeropropulsive CFD simulation. This low-boom redesign study is used as an example to demonstrate how the propulsion-airframe integration could be implemented for conceptual design of supersonic transports that satisfy both the low-boom and mission performance requirements.

Low-boom supersonic transport↗
NASA NTRSDate not supplied

Standalone Hazard Evaluation and Refinement From Instrument Findings (S.H.E.R.I.F)

SHERIF is a novel algorithm which defines a unique, non-iterative method for the combination of LiDAR scans into a Digital Elevation Map (DEM) and the simultaneous evaluation of that DEM for Safe Site Selection (SSL) and Hazard Detection (HD). This creates a sensor independent method for a spacecraft conducting Entry Descent and Landing (EDL) operations to synthesize terrain scans and generate a more complete and evolving understanding of the landing site topography as well as the associated hazards. SHERIF aims to increase safety and mission success probability by more effectively using the data from a single sensor to increase situational awareness, knowledge of the landing site and its hazards, as well as informing a more streamlined transition from Terrain Relative Navigation (TRN) to Hazard Relative Navigation (HRN) during descent.

Hazard Detection↗
NASA NTRS1976

The Davidon-Fletcher-Powell penalty function method: A generalized iterative technique for solving parameter optimization problems

The Fletcher-Powell version of the Davidon variable metric unconstrained minimization technique is described. Equations that have been used successfully with the Davidon-Fletcher-Powell penalty function technique for solving constrained minimization problems and the advantages and disadvantages of using them are discussed. The experience gained in the behavior of the method while iterating is also related.

Johnson, I. L., Jr.↗
NASA NTRS1990

Fast secant methods for the iterative solution of large nonsymmetric linear systems

A family of secant methods based on general rank-1 updates was revisited in view of the construction of iterative solvers for large non-Hermitian linear systems. As it turns out, both Broyden's good and bad update techniques play a special role, but should be associated with two different line search principles. For Broyden's bad update technique, a minimum residual principle is natural, thus making it theoretically comparable with a series of well known algorithms like GMRES. Broyden's good update technique, however, is shown to be naturally linked with a minimum next correction principle, which asymptotically mimics a minimum error principle. The two minimization principles differ significantly for sufficiently large system dimension. Numerical experiments on discretized partial differential equations of convection diffusion type in 2-D with integral layers give a first impression of the possible power of the derived good Broyden variant.

Deuflhard, Peter↗
NASA NTRS1990

Linear iterative solvers for implicit ODE methods

The numerical solution of stiff initial value problems, which lead to the problem of solving large systems of mildly nonlinear equations are considered. For many problems derived from engineering and science, a solution is possible only with methods derived from iterative linear equation solvers. A common approach to solving the nonlinear equations is to employ an approximate solution obtained from an explicit method. The error is examined to determine how it is distributed among the stiff and non-stiff components, which bears on the choice of an iterative method. The conclusion is that error is (roughly) uniformly distributed, a fact that suggests the Chebyshev method (and the accompanying Manteuffel adaptive parameter algorithm). This method is described, also commenting on Richardson's method and its advantages for large problems. Richardson's method and the Chebyshev method with the Mantueffel algorithm are applied to the solution of the nonlinear equations by Newton's method.

Saylor, Paul E.↗
DOE OSTI2025

An iterative dynamic chemical stiffness removal method for reacting flow simulations

Abstract An iterative dynamic chemical stiffness removal method (IDCSR) based on quasi-steady-state approximation (QSSA) is proposed. The IDCSR method is built on a previously developed non-iterative method which has proved to work well for small timestep sizes. A novel iterative procedure is designed in IDCSR to enable explicit time integration of stiff chemistry at relatively large timestep sizes relevant to practical reacting flow simulations. The effectiveness of the iterative procedure is first demonstrated with a toy problem and homogeneous auto-ignition with fixed integration step sizes, showing that larger timestep sizes can be allowed for explicit time integration using IDCSR compared with the previous non-iterative method. IDCSR is then compared with existing explicit chemistry solvers for simulations of homogeneous auto-ignition and shows similar or lower computational cost but significantly higher accuracy across a wide range of timestep sizes. IDCSR is further combined with an automatic adaptive time-stepping scheme for simulations of 0-D homogeneous auto-ignition and a 2-D laminar lifted n -dodecane jet flame. For the 0-D auto-ignition simulations, IDCSR is shown to reduce both the error (by 43%–90%) and computational cost (by 6–15 times) compared with existing explicit solvers, while achieving speed-up factors of up to 400 compared with VODE for a wide range of timestep sizes and reaction mechanisms. For the 2-D jet flame simulations, speed-up factors of 15 and 31 for chemistry integration, and 5 and 9 for overall simulation, are achieved by IDCSR compared with CVODE with and without analytic Jacobian, respectively.

Xu, Chao (ORCID:0000000153074159)↗
NASA NTRS1972

A review of nonstatistical techniques for the estimation of vertical atmospheric structure from remote infrared measurements

Mathematical techniques used in atmospheric profile inversion are reviewed, with the discussion limited to methods that do not require extensive information on the statistical properties of the profiles. The basic problem of the calculation of vertical profiles from noise contaminated infrared spectral measurements is examined. Three specific examples of methods of solution are given: the iterative minimum estimation method, a second linear iterative method, and the algorithm of Chahine. Satellite data are used in an empirical comparison of the methods. Several factors are considered, including the stability of solutions, convergence behavior of the iterative methods, and dependence of the solutions on the first guess, a factor that is particularly important in the treatment of data from other planets.

Conrath, B. J.↗
NASA NTRS1990

Accumulated approximation: A new method for structural optimization by iterative improvement

A new method for the solution of non-linear mathematical programming problems in the field of structural optimization is presented. It is an iterative scheme which for each iteration refines the approximation of objective and constraint functions by accumulating the function values of previously visited design points. The method has proven to be competitive for a number of well-known examples of which one is presented here. Furthermore because of the accumulation strategy, the method produces convergence even when the sensitivity analysis is inaccurate.

Rasmussen, John↗
NASA NTRS1996

Multigrid Methods for EHL Problems

In many bearings and contacts, forces are transmitted through thin continuous fluid films which separate two contacting elements. Objects in contact are normally subjected to friction and wear which can be reduced effectively by using lubricants. If the lubricant film is sufficiently thin to prevent the opposing solids from coming into contact and carries the entire load, then we have hydrodynamic lubrication, where the lubricant film is determined by the motion and geometry of the solids. However, for loaded contacts of low geometrical conformity, such as gears, rolling contact bearings and cams, this is not the case due to high pressures and this is referred to as Elasto-Hydrodynamic Lubrication (EHL) In EHL, elastic deformation of the contacting elements and the increase in fluid viscosity with pressure are very significant and cannot be ignored. Since the deformation results in changing the geometry of the lubricating film, which in turn determines the pressure distribution, an EHL mathematical model must simultaneously satisfy the complex elasticity (integral) and the Reynolds lubrication (differential) equations. The nonlinear and coupled nature of the two equations makes numerical calculations computationally intensive. This is especially true for highly loaded problems found in practice. One novel feature of these problems is that the solution may exhibit sharp pressure spikes in the outlet region. To this date both finite element and finite difference methods have been used to solve EHL problems with perhaps greater emphasis on the use of the finite difference approach. In both cases, a major computational difficulty is ensuring convergence of the nonlinear equations solver to a steady state solution. Two successful methods for achieving this are direct iteration and multigrid methods. Direct iteration methods (e.g Gauss Seidel) have long been used in conjunction with finite difference discretizations on regular meshes. Perhaps one of the best examples of the application of such methods is the recent Effective Influence Method of Dowson and Wang. Multigrid methods have also been used with great success by Venner and Venner and Lubrecht with a good summary being given by Venner. As both these finite difference discretization based approaches appear to provide an efficient way of solving EHL problems, it is important to understand their relative merits. This paper is a first attempt at providing such an understanding in the context of EHL point contact problem, (contact of two spheres), in which the contact zone is a point and an ellipse or circle for unloaded and loaded dry contacts respectively. Since the film thickness and the contact width are generally small compared to the local radius of curvature of the two surfaces, the reduced geometry of the surfaces in the contact area can be accurately approximated to the contact between a paraboloid and a flat surface. The layout of the remainder of this paper is as follows. In section 2 we introduce the form of the equations to be solved. The Effective Influence Newton Method is described in Section 3 while Section 4 describes the Multigrid method to be used. Sections 5 and 6 describe the test problems to be used in the comparison between the two methods and compare the performance of the two methods. Section 7 concludes the paper with an argument of the two methods and suggests some future research directions.

Nurgat, Elyas↗
NASA NTRS1982

Computer tomography of flows external to test models

Computer tomographic techniques for reconstruction of three-dimensional aerodynamic density fields, from interferograms recorded from several different viewing directions were studied. Emphasis is on the case in which an opaque object such as a test model in a wind tunnel obscures significant regions of the interferograms (projection data). A method called the Iterative Convolution Method (ICM), existing methods in which the field is represented by a series expansions, and analysis of real experimental data in the form of aerodynamic interferograms are discussed.

Prikryl, I.↗
NASA NTRS1988

Domain decomposition preconditioners for the spectral collocation method

Several block iteration preconditioners are proposed and analyzed for the solution of elliptic problems by spectral collocation methods in a region partitioned into several rectangles. It is shown that convergence is achieved with a rate which does not depend on the polynomial degree of the spectral solution. The iterative methods here presented can be effectively implemented on multiprocessor systems due to their high degree of parallelism.

Quarteroni, Alfio↗
NASA NTRS1984

A statistical correlation method for the retrieval of atmospheric moisture profiles by microwave radiometry

A statistical correlation technique is applied to the retrieval of vertical moisture profiles under clear-sky conditions from down-looking radiometric measurements of atmospheric radiation at microwave wavelengths. For a given set of channels, the method selects the optimum radiometric channels for estimating water vapor at specific pressure levels between the surface and 300 mb. The water vapor mixing ratio at these pressure levels is then calculated from a linear combination of the selected channel brightness temperatures. To test its validity the algorithm was applied, in a numerical experiment, to fifty independent tropical radiosondes. The rms absolute deviation of the estimated moisture profiles from the actual profiles was comparable to that obtained using an iterative retrieval method reported earlier. The statistical method, however, requires several orders of magnitude less computer time than the iterative method; it is suitable for high speed processing of large amounts of data.

Kakar, R. K.↗
NASA NTRS1987

A new art code for tomographic interferometry

A new algebraic reconstruction technique (ART) code based on the iterative refinement method of least squares solution for tomographic reconstruction is presented. Accuracy and the convergence of the technique is evaluated through the application of numerically generated interferometric data. It was found that, in general, the accuracy of the results was superior to other reported techniques. The iterative method unconditionally converged to a solution for which the residual was minimum. The effects of increased data were studied. The inversion error was found to be a function of the input data error only. The convergence rate, on the other hand, was affected by all three parameters. Finally, the technique was applied to experimental data, and the results are reported.

Tan, H.↗
NASA NTRS1996

Singularity Preserving Numerical Methods for Boundary Integral Equations

In the past twelve months (May 8, 1995 - May 8, 1996), under the cooperative agreement with Division of Multidisciplinary Optimization at NASA Langley, we have accomplished the following five projects: a note on the finite element method with singular basis functions; numerical quadrature for weakly singular integrals; superconvergence of degenerate kernel method; superconvergence of the iterated collocation method for Hammersteion equations; and singularity preserving Galerkin method for Hammerstein equations with logarithmic kernel. This final report consists of five papers describing these projects. Each project is preceeded by a brief abstract.

Kaneko, Hideaki↗
NASA NTRS1973

Rapid iterative reanalysis for automated design

A method for iterative reanalysis in automated structural design is presented for a finite-element analysis using the direct stiffness approach. A basic feature of the method is that the generalized stiffness and inertia matrices are expressed as functions of structural design parameters, and these generalized matrices are expanded in Taylor series about the initial design. Only the linear terms are retained in the expansions. The method is approximate because it uses static condensation, modal reduction, and the linear Taylor series expansions. The exact linear representation of the expansions of the generalized matrices is also described and a basis for the present method is established. Results of applications of the present method to the recalculation of the natural frequencies of two simple platelike structural models are presented and compared with results obtained by using a commonly applied analysis procedure used as a reference. In general, the results are in good agreement. A comparison of the computer times required for the use of the present method and the reference method indicated that the present method required substantially less time for reanalysis. Although the results presented are for relatively small-order problems, the present method will become more efficient relative to the reference method as the problem size increases. An extension of the present method to static reanalysis is described, ana a basis for unifying the static and dynamic reanalysis procedures is presented.

Bhatia, K. G.↗
NASA NTRS1992

Transonic airfoil and wing design using Navier-Stokes codes

An iterative design method has been implemented into 2D and 3D Navier-Stokes codes for the design of airfoils or wings with given target pressure distributions. The method begins with the analysis of an initial geometry, and obtains the analysis pressure distributions of that geometry. The differences between analysis pressures and target pressures are used to drive geometry changes through the use of a streamline curvature method. This paper describes the procedure that makes the iterative design method work for Navier-Stokes codes. Examples of 2D airfoil design, and 3D wing design are included. It is demonstrated that the method is highly effective for airfoil or wing design at flow conditions where no substantial separation occurs. Problems encountered in the airfoil design with shock induced flow separations are discussed.

Yu, N. J.↗
NASA NTRS1990

A transient response analysis of the space shuttle vehicle during liftoff

A proposed transient response method is formulated for the liftoff analysis of the space shuttle vehicles. It uses a power series approximation with unknown coefficients for the interface forces between the space shuttle and mobile launch platform. This allows the equation of motion of the two structures to be solved separately with the unknown coefficients at the end of each step. These coefficients are obtained by enforcing the interface compatibility conditions between the two structures. Once the unknown coefficients are determined, the total response is computed for that time step. The method is validated by a numerical example of a cantilevered beam and by the liftoff analysis of the space shuttle vehicles. The proposed method is compared to an iterative transient response analysis method used by Martin Marietta for their space shuttle liftoff analysis. It is shown that the proposed method uses less computer time than the iterative method and does not require as small a time step for integration. The space shuttle vehicle model is reduced using two different types of component mode synthesis (CMS) methods, the Lanczos method and the Craig and Bampton CMS method. By varying the cutoff frequency in the Craig and Bampton method it was shown that the space shuttle interface loads can be computed with reasonable accuracy. Both the Lanczos CMS method and Craig and Bampton CMS method give similar results. A substantial amount of computer time is saved using the Lanczos CMS method over that of the Craig and Bampton method. However, when trying to compute a large number of Lanczos vectors, input/output computer time increased and increased the overall computer time. The application of several liftoff release mechanisms that can be adapted to the proposed method are discussed.

Brunty, J. A.↗