Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iteration 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 973 records · Page 54

Infrared band extinctions and complex refractive indices of crystalline C2H2 and C4H2

Thermal IR absorption intensities are obtained for thin films of crystalline C2H2 and C4H2 at 70 K, and their n and k complex refractive indices are ascertained by separating true film absorption from interface reflection on the basis of an analysis of the transmission spectrum ratio for two sample thicknesses. This method significantly simplifies the n and k iteration process. The n and k values determined in the laboratory will in most cases reproduce a given sample thickness' observed transmission to within + or - 5 percent.

Khanna, R. K.↗

Chemical nonequilibrium and viscous flow computation for conic aerobrake bodies

Three-dimensional analyses are presented for the viscous, reactive flow over a complete entry-body configuration with a wide-angle conic surface. The predictive method uses a split approach that solves iteratively the Navier-Stokes and the continuity equations of chemical species. The finite-difference formulation and the computational grid are adapted to the bow shock and the conformally mapped body such that the velocity components are in the computational spherical-polar space. Combinations of several conic forebody and afterbody configurations have been studied using wind-tunnel, Space Shuttle, and aerobraking orbital transfer vehicle (AOTV) entry conditions. The effects of the borebody bluntness and of finite-rate chemical reactions on the shock layer, the wall catalycity on the boundary layer, the shear-layer impingement on the afterbody, and the base-flow environment are discussed.

Li, C. P.↗

Efficient real gas Navier-Stokes computations of high speed flows using an LU scheme

An efficient method to account for the chemically frozen thermodynamic and transport properties of air in three dimensional Navier-Stokes calculations was demonstrated. This approach uses an explicitly specified equation of state (EOS) so that the fluid pressure, temperature and transport properties are directly related to the flow variables. Since the pressure is explicitly known as a general function of the flow variables no assumptions are made regarding the pressure derivatives in the construction of the flux Jacobians. The method is efficient since no sub-iterations are required to deduce the pressure and temperature from the flux variables and allows different equations of state to be easily supplied to the code. The flexibility of the EOS approach is demonstrated by implementing a high order TVD upwinding scheme based upon flux differencing and Van Leer's flux vector splitting. The EOS approach is demonstrated by computing the hypersonic flow through the corner region of two mutually perpendicular flat plates and through a simplified model of a scramjet module gap-seal configuration.

Coirier, William J.↗

An application generator for rapid prototyping of Ada real-time control software

The need to increase engineering productivity and decrease software life cycle costs in real-time system development establishes a motivation for a method of rapid prototyping. The design by iterative rapid prototyping technique is described. A tool which facilitates such a design methodology for the generation of embedded control software is described.

Johnson, Jim↗

Efficient real gas Navier-Stokes computations of high speed flows using an LU scheme

An efficient method to account for the chemically frozen thermodynamic and transport properties of air in three dimensional Navier-Stokes calculations was demonstrated. This approach uses an explicitly specified equation of state (EOS) so that the fluid pressure, temperature and transport properties are directly related to the flow variables. Since the pressure is explicitly known as a general function of the flow variables no assumptions are made regarding the pressure derivatives in the construction of the flux Jacobians. The method is efficient since no sub-iterations are required to deduce the pressure and temperature from the flux variables and allows different equations of state to be easily supplied to the code. The flexibility of the EOS approach is demonstrated by implementing a high order TVD upwinding scheme based upon flux differencing and Van Leer's flux vector splitting. The EOS approach is demonstrated by computing the hypersonic flow through the corner region of two mutually perpendicular flat plates and through a simplified model of a scramjet module gap-seal configuration.

Coirier, William J.↗

A tool for interactive exploration of a hierarchical segmentation

Hierarchical segmentation is discussed as a form of region growing in which the sequence of merges is controlled by a 'best merge first' principle, and a record of the region merging sequence is often retained for later analysis. In the author's massively parallel implementation of hierarchical segmentation, which he calls iterative parallel region growing (IPRG), a set of directional edge maps are used to store the region merging sequence information. An iteractive tool is described that allows an analyst to fully explore a hierarchical segmentation, potentially producing an image segmentation that is a combination of segmentations produced at several different iterations of the IPRG algorithm. A potential method for automatically producing such segmentations is also discussed.

Tilton, James C.↗

Symmetric Galerkin boundary formulations employing curved elements

Accounts of the symmetric Galerkin approach to boundary element analysis (BEA) have recently been published. This paper attempts to add to the understanding of this method by addressing a series of fundamental issues associated with its potential computational efficiency. A new symmetric Galerkin theoretical formulation for both the (harmonic) heat conduction and the (biharmonic) elasticity problem that employs regularized singular and hypersingular boundary integral equations (BIEs) is presented. The novel use of regularized BIEs in the Galerkin context is shown to allow straightforward incorporation of curved, isoparametric elements. A symmetric reusable intrinsic sample point (RISP) numerical integration algorithm is shown to produce a Galerkin (i.e., double) integration strategy that is competitive with its counterpart (i.e., singular) integration procedure in the collocation BEA approach when the time saved in the symmetric equation solution phase is also taken into account. This new formulation is shown to be capable of employing hypersingular BIEs while obviating the requirement of C 1 continuity, a fact that allows the employment of the popular continuous element technology. The behavior of the symmetric Galerkin BEA method with regard to both direct and iterative equation solution operations is also addressed. A series of example problems are presented to quantify the performance of this symmetric approach, relative to the more conventional unsymmetric BEA, in terms of both accuracy and efficiency. It is concluded that appropriate implementations of the symmetric Galerkin approach to BEA indeed have the potential to be competitive with, if not superior to, collocation-based BEA, for large-scale problems.

Kane, J. H.↗

Transonic Flow Computations Using Nonlinear Potential Methods

This presentation describes the state of transonic flow simulation using nonlinear potential methods for external aerodynamic applications. The presentation begins with a review of the various potential equation forms (with emphasis on the full potential equation) and includes a discussion of pertinent mathematical characteristics and all derivation assumptions. Impact of the derivation assumptions on simulation accuracy, especially with respect to shock wave capture, is discussed. Key characteristics of all numerical algorithm types used for solving nonlinear potential equations, including steady, unsteady, space marching, and design methods, are described. Both spatial discretization and iteration scheme characteristics are examined. Numerical results for various aerodynamic applications are included throughout the presentation to highlight key discussion points. The presentation ends with concluding remarks and recommendations for future work. Overall. nonlinear potential solvers are efficient, highly developed and routinely used in the aerodynamic design environment for cruise conditions. Published by Elsevier Science Ltd. All rights reserved.

Holst, Terry L.↗

Flutter Prediction for Aircraft Conceptual Design

Flutter prediction is usually a knowledge-based analysis process that aims to reduce the cost of aeroelastic stability margin certification. However, early detection of flutter problems is beneficial in the development of unconventional aircraft. The recently developed automation tool ConceptFEA for structural sizing of aircraft concepts paves the way for rapid physics-based flutter prediction of aircraft concepts. A match-point iteration procedure using the p-k method is implemented for ConceptFEA with minimum user input requirements to generate flutter boundary points. A subsonic business jet concept and its high aspect-ratio wing variant are used to demonstrate how the newly developed flutter prediction capability can be used during aircraft conceptual design. Sized structures, flutter boundary curves, and flutter sensitivity analysis results are generated for these two concepts using ConceptFEA. The relevant equivalent plate theory is provided to show the quantitative relationships between a stiffened panel and its equivalent NASTRAN PSHELL panel. The rapid flutter prediction capability of ConceptFEA makes multidisciplinary collaborations between systems analysts and aeroelasticity experts feasible in practice.

Li, Wu↗

NASA Kennedy Space Center Swamp Works: Innovation Methods for Advanced Research & Technology Development

Kennedy Space Center’s (KSC) Swamp Works, provides government and commercial space ventures with the technologies required for working and living on the surfaces of the Moon or other planets and bodies in our solar system. The Swamp Works team establishes agile, innovative and cost effective exploration mission solutions through leveraging of partnerships across NASA, industry and academia. Concepts start small and build up efficiently, with lean development processes and a hands on approach. Testing is performed in early stages to drive design improvements and progressively increase Technology Readiness Levels (TRL). Swamp Works provides concepts, architecture studies and trades, designs, data, technology development, technology demonstration hardware, flight hardware, testing, flight support and knowledge in support of the development of surface systems. This paper will summarize the state of the art in innovation methods, and it will compare the various other methods to how innovation is accomplished at the NASA KSC Swamp Works. Since 2013, Swamp Works has developed a wide variety of technologies, and has demonstrated viable methods for innovation while also collaborating with other government agencies, industry and academia to support space exploration. These methods focus on agile design, test, build iterations where each cycle of prototyping results in improved knowledge and designs. Failures are tole rated in the early stages of these iterations, since the overall development cycle is shorter more efficient and cost effective by retiring risk early in the cycle. Momentum is maintained and the development teams remain highly motivated by the high level of team work, technical excellence, engagement and creative freedom.

Robotic Mining↗

Nested Krylov methods and preserving the orthogonality

Recently the GMRESR inner-outer iteraction scheme for the solution of linear systems of equations was proposed by Van der Vorst and Vuik. Similar methods have been proposed by Axelsson and Vassilevski and Saad (FGMRES). The outer iteration is GCR, which minimizes the residual over a given set of direction vectors. The inner iteration is GMRES, which at each step computes a new direction vector by approximately solving the residual equation. However, the optimality of the approximation over the space of outer search directions is ignored in the inner GMRES iteration. This leads to suboptimal corrections to the solution in the outer iteration, as components of the outer iteration directions may reenter in the inner iteration process. Therefore we propose to preserve the orthogonality relations of GCR in the inner GMRES iteration. This gives optimal corrections; however, it involves working with a singular, non-symmetric operator. We will discuss some important properties, and we will show by experiments that, in terms of matrix vector products, this modification (almost) always leads to better convergence. However, because we do more orthogonalizations, it does not always give an improved performance in CPU-time. Furthermore, we will discuss efficient implementations as well as the truncation possibilities of the outer GCR process. The experimental results indicate that for such methods it is advantageous to preserve the orthogonality in the inner iteration. Of course we can also use iteration schemes other than GMRES as the inner method; methods with short recurrences like GICGSTAB are of interest.

Desturler, Eric↗

Comparative study of the convergence rates of two numerical techniques

The paper examines the applicability of the three-step Stetter (1968) method to the problem of hypersonic viscous flow over a blunt axisymmetric body used for planetary entry probes at zero angle of attack. The flow-field results using the two-step finite-difference MacCormack (1969) method are reported by Kumar and Graves (1977). Only the computational efficiency of Stetter's method is compared with that of MacCormack's in terms of the iterative time steps and computing time required for the steady-state solution. Advantages of Stetter's method over MacCormack's are established.

Kumar, A.↗

Analysis of energy states in modulation doped multiquantum well heterostructures

A precise and effective numerical procedure to model the band diagram of modulation doped multiquantum well heterostructures is presented. This method is based on a self-consistent iterative solution of the Schroedinger equation and the Poisson equation. It can be used rather easily in any arbitrary modulation-doped structure. In addition to confined energy subbands, the unconfined states can be calculated as well. Examples on realistic device structures are given to demonstrate capabilities of this procedure. The numerical results are in good agreement with experiments. With the aid of this method the transitions involving both the confined and unconfined conduction subbands in a modulation doped AlGaAs/GaAs superlattice, and in a strained layer InGaAs/GaAs superlattice are identified. These results represent the first observation of unconfined transitions in modulation doped multiquantum well structures.

Ji, G.↗

Temperature retrievals from TIROS-N

A nonstatistical, nonlinear iterative scheme, has been developed, based on the methods of Chahine (1970), to retrieve surface temperature, vertical temperature profile and cloud distribution from TIROS-N HIRS 2/MSU observations. Fundamental to the technique is the ability to accurately calculate clear column radiances given atmospheric and surface characteristics. The cloud filtering algorithm, making simultaneous use of IR 15 micron channels and microwave channels, allows accurate infrared retrievals in the presence of up to 80 percent cloud cover. Details for the iterative scheme are presented together with preliminary results from analysis of TIROS-N data.

Susskind, J.↗

Ames Optimized TCA Configuration

Configuration design at Ames was carried out with the SYN87-SB (single block) Euler code using a 193 x 49 x 65 C-H grid. The Euler solver is coupled to the constrained (NPSOL) and the unconstrained (QNMDIF) optimization packages. Since the single block grid is able to model only wing-body configurations, the nacelle/diverter effects were included in the optimization process by SYN87's option to superimpose the nacelle/diverter interference pressures on the wing. These interference pressures were calculated using the AIRPLANE code. AIRPLANE is an Euler solver that uses a unstructured tetrahedral mesh and is capable of computations about arbitrary complete configurations. In addition, the buoyancy effects of the nacelle/diverters were also included in the design process by imposing the pressure field obtained during the design process onto the triangulated surfaces of the nacelle/diverter mesh generated by AIRPLANE. The interference pressures and nacelle buoyancy effects are added to the final forces after each flow field calculation. Full details of the (recently enhanced) ghost nacelle capability are given in a related talk. The pseudo nacelle corrections were greatly improved during this design cycle. During the Ref H and Cycle 1 design activities, the nacelles were only translated and pitched. In the cycle 2 design effort the nacelles can translate vertically, and pitch to accommodate the changes in the lower surface geometry. The diverter heights (between their leading and trailing edges) were modified during design as the shape of the lower wing changed, with the drag of the diverter changing accordingly. Both adjoint and finite difference gradients were used during optimization. The adjoint-based gradients were found to give good direction in the design space for configurations near the starting point, but as the design approached a minimum, the finite difference gradients were found to be more accurate. Use of finite difference gradients was limited by the CPU time limit available on the Cray machines. A typical optimization run using finite difference gradients can use only 30 to 40 design variables and one optimization iteration within the 8 hour queue limit for the chosen grid size and convergence level. The efficiency afforded by the adjoint method allowed for 50-120 design variables and 5-10 optimization iterations in the 8 hour queue. Geometric perturbations to the wing and fuselage were made using the Hicks/Henne (HH) shape functions. The HH functions were distributed uniformly along the chords of the wing defining sections and lofted linearly. During single-surface design, constraints on thickness and volume at selected wing stations were imposed. Both fuselage camber and cross-sectional area distributions were permitted to change during design. The major disadvantage to the use of these functions is the inherent surface waviness produced by repeated use of such functions. Many smoothing operations were required following optimization runs to produce a configuration with reasonable smoothness. Wagner functions were also used on the wing sections but were never used on the fuselage. The Wagner functions are a family of increasingly oscillatory functions that have also been used extensively in airfoil design. The leading and trailing edge regions of the wing were designed by use of polynomial and monomial functions respectively. Twist was attempted but was abandoned because of little performance improvement available from changing the baseline twist.

Cliff, Susan E.↗

Streamlined Convergence Acceleration for CFD Codes

Enigma, a simplified interface to the PETSc library, is shown to enable the rapid solution of discrete partial differential equations. Two CFD codes, LAURA and HyperSolve, use Enigma to compute steady solutions of the Navier-Stokes equations. Using PETSc, Enigma is shown to provide a Jacobian-Free Newton-Krylov method (JFNK), globalized with pseudotransient continuation, that improves efficiency over the point-implicit relaxation method traditionally used by LAURA. It is shown that iterative error has a large impact on surface heat transfer predicted by LAURA on an axisymmetric sphere-cone geometry. Also, the convergence rate of HyperSolve simulating subsonic flow over a delta wing geometry with the JFNK method is shown to be more efficient than employing a defect correction method as the nonlinear solver.

Thompson, Kyle B.↗

Effects of nonlinear aerodynamics and static aeroelasticity on mission performance calculations for a fighter aircraft

During conceptual design studies of advanced aircraft, the usual practice is to use linear theory to calculate the aerodynamic characteristics of candidate rigid (nonflexible) geometric external shapes. Recent developments and improvements in computational methods, especially computational fluid dynamics (CFD), provide significantly improved capability to generate detailed analysis data for the use of all disciplines involved in the evaluation of a proposed aircraft design. A multidisciplinary application of such analysis methods to calculate the effects of nonlinear aerodynamics and static aeroelasticity on the mission performance of a fighter aircraft concept is described. The aircraft configuration selected for study was defined in a previous study using linear aerodynamics and rigid geometry. The results from the previous study are used as a basis of comparison for the data generated herein. Aerodynamic characteristics are calculated using two different nonlinear theories, potential flow and rotational (Euler) flow. The aerodynamic calculations are performed in an iterative procedure with an equivalent plate structural analysis method to obtain lift and drag data for a flexible (nonrigid) aircraft. These static aeroelastic data are then used in calculating the combat and mission performance characteristics of the aircraft.

Giles, Gary L.↗

VIPER Lunar Rover Agile Mission Systems

Agile development methods, which have gone from outlier to mainstream in software development, are poised to expand into all aspects of space mission development. Modern software development operates on a principle of continuous deployment, where progress is verified not with conventional metrics, but with a continuous build, available to key stakeholders, enabling direct examination of the state of the code base, and assessment of progress through demonstration of capability. Delivery times are measured in weeks, not months. Stakeholders are part of the process on an ongoing basis. The cost of change is comparatively low and requirements, which often are not precisely defined at the start of a project, may be iteratively refined in a series of agile development cycles. Agile methods are compatible with traditional system engineering methods and may be tailored to the space operations environment. The low cost of change and iterative development cycles of agile enable requirements to be defined as outcomes and constraints, with design details to be refined during the development cycle. We are now at a point where agile methods may be extended beyond software, to Mission Systems, including the Mission Operations System and the Ground Data System. For NASA’s VIPER Lunar Rover Mission, scheduled to land at a lunar pole in late 2023, we are developing the Mission System using agile methods. As in agile software, where the measure of progress is working code, in agile mission system development, the measure of capability is what we can demonstrate. Demonstrations over presentations. We demonstrate mission system capability using simulations. The concept of operations, from commanding, to driving the rover, to how we downlink images for evaluation for a near-real time command cycle, will be tested and proven in simulation, years before we begin the traditional simulation cycle for training. “Say it then simulate it.” We develop and refine our designs using simulations, with an emphasis on new components of the system that are not well known early. For example, the required duration of a mission planning cycle for a lunar surface asset such as VIPER, that operates twenty-four hours a day, seven days a week, with continuous communications and a unique set of constraints based on the physics of the lunar poles and the line of site to Earth, is a unique problem in mission planning that is unlikely to be solved in a series of meetings. A small number of requirements specifying the outcomes may serve as the jumping off point to an agile development cycle, with demonstration in simulations. We have already demonstrated this process with simulations of rover driver decision time. VIPER is driven using near-real time command and control to waypoints. The driver decision time between waypoints is a fundamental enabling unit of productivity to accomplish the mission timeline. We have validated driver decision time in simulations of rover driving at the lunar South Pole, using the prototype mission tools for driving, command and control. The capability to develop and refine designs using simulations as part of agile Mission System development cycle changes the nature of team interactions, creating a focus on doing, rather than analyzing and documenting. Waterfall development cycles were, in part, a product of the significant cost of change in the early days of spaceflight. When the cost of change is high, it is vital to get your requirements right at the outset, because the system will be built to those specifications, and, when change is expensive, you better get it right early. However, modern technology has greatly lowered the cost of change, enabling iterative, rapid development cycles, in which key operations concepts may be tested and refined during development. Extending agile development to the Mission System for VIPER is a significant step in moving agile development methods for space operations beyond software, to the Mission System.

Agile↗