Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solvers”

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 1,531 records · Page 85

Design and experimental evaluation of compact radial-inflow turbines

The application of a multistage 3D Euler solver to the aerodynamic design of two compact radial-inflow turbines is presented, along with experimental results evaluating and validating the designs. The objectives of the program were to design, fabricate, and rig test compact radial-inflow turbines with equal or better efficiency relative to conventional designs, while having 40 percent less rotor length than current traditionally-sized radial turbines. The approach to achieving these objectives was to apply a calibrated 3D multistage Euler code to accurately predict and control the high rotor flow passage velocities and high aerodynamic loadings resulting from the reduction in rotor length. A comparison of the advanced compact designs to current state-of-the-art configurations is presented.

Fredmonski, A. J.↗

Solution-adaptive program SADAP3D

A generic solution-adaptive grid program based on an error equi-distribution concept has been developed for use in complex multidimensional flows. The capability of a generic, user-friendly, multidimensional algorithm is demonstrated using several complex 3D flows in high-speed regimes. A simple scheme enforces the continuity of grid spacing along each currently adapted grid line as the computation proceeds. Calculation of normal vectors at each point of a grid plane is enhanced. The grid movement is automatically controlled by reinstating grid-spacing continuity locally in the mainstream of the unidirectional computation. It is concluded that the scheme can be readily coupled with many flow solvers to enhance accuracy and efficiency of solutions.

Djomehri, M. J.↗

Space Station resource node flow field analysis

An analysis of the flow field within the Space Station Freedom resource node with operational intermodule ventilation and temperature/humidity control ventilation systems has been conducted. The INS3D code, an incompressible, steady-state Navier-Stokes solver has been used to assess the design of the ventilation system via quantification of the level of fluid mixing and identification of 'dead air' regions and short-circuit ventilation. Numerical results indicate significant short-circuit ventilation in the forward and midsections of the node and insufficient fluid mixing is found to exist in the aft node section. These results as well as results from a solution grid dependence study are presented.

Kania, Lee↗

Aerodynamic analysis of complex configurations using unstructured grids

The purpose of this paper is to assess the accuracy and utility of a new unstructured, inviscid, upwind flow solver for the aerodynamic analysis of two aircraft configurations. The two configurations consist of a low-wing transport with nacelle/pylon on and off, and a generic high-speed civil transport. Computations are made at subsonic and transonic Mach numbers for the low-wing transport and at transonic and low-supersonic speeds for the high-speed civil transport. The results include an assessment of grid sensitivity and provide comparisons with experimental data.

Frink, Neal T.↗

Wing design code using three-dimensional Euler equations and optimization

This paper describes a new wing design code which is based on the Euler equations and a constrained numerical optimization technique. The geometry modification is based on a set of fundamental modes define on the unit interval. A design example involving a high speed civil transport wing is presented to demonstrate the usefulness of the design code. It is shown that the use of an Euler solver in the direct numerical optimization procedures is affordable on the current generation of supercomputers.

Chang, I-Chung↗

Recent literature on structural modeling, identification, and analysis

The literature on the mathematical modeling of large space structures is first reviewed, with attention given to continuum models, model order reduction, substructuring, and computational techniques. System identification and mode verification are then discussed with reference to the verification of mathematical models of large space structures. In connection with analysis, the paper surveys recent research on eigensolvers and dynamic response solvers for large-order finite-element-based models.

Craig, Roy R., Jr.↗

An assessment of computational fluid dynamic techniques in the analysis and design of turbomachinery - The 1990 Freeman Scholar Lecture

Various computational fluid dynamic techniques are reviewed focusing on the Euler and Navier-Stokes solvers with a brief assessment of boundary layer solutions, and quasi-3D and quasi-viscous techniques. Particular attention is given to a pressure-based method, explicit and implicit time marching techniques, a pseudocompressibility technique for incompressible flow, and zonal techniques. Recommendations are presented with regard to the most appropriate technique for various flow regimes and types of turbomachinery, incompressible and compressible flows, cascades, rotors, stators, liquid-handling, and gas-handling turbomachinery.

Lakshminarayana, B.↗

An implementation of the look-ahead Lanczos algorithm for non-Hermitian matrices, part 2

It is shown how the look-ahead Lanczos process (combined with a quasi-minimal residual QMR) approach) can be used to develop a robust black box solver for large sparse non-Hermitian linear systems. Details of an implementation of the resulting QMR algorithm are presented. It is demonstrated that the QMR method is closely related to the biconjugate gradient (BCG) algorithm; however, unlike BCG, the QMR algorithm has smooth convergence curves and good numerical properties. We report numerical experiments with our implementation of the look-ahead Lanczos algorithm, both for eigenvalue problem and linear systems. Also, program listings of FORTRAN implementations of the look-ahead algorithm and the QMR method are included.

Freund, Roland W.↗

An efficient and robust algorithm for two dimensional time dependent incompressible Navier-Stokes equations: High Reynolds number flows

An algorithm is presented for unsteady two-dimensional incompressible Navier-Stokes calculations. This algorithm is based on the fourth order partial differential equation for incompressible fluid flow which uses the streamfunction as the only dependent variable. The algorithm is second order accurate in both time and space. It uses a multigrid solver at each time step. It is extremely efficient with respect to the use of both CPU time and physical memory. It is extremely robust with respect to Reynolds number.

Goodrich, John W.↗

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↗

Optimal parallel solution of sparse triangular systems

A method for the parallel solution of triangular sets of equations is described that is appropriate when there are many right-handed sides. By preprocessing, the method can reduce the number of parallel steps required to solve Lx = b compared to parallel forward or backsolve. Applications are to iterative solvers with triangular preconditioners, to structural analysis, or to power systems applications, where there may be many right-handed sides (not all available a priori). The inverse of L is represented as a product of sparse triangular factors. The problem is to find a factored representation of this inverse of L with the smallest number of factors (or partitions), subject to the requirement that no new nonzero elements be created in the formation of these inverse factors. A method from an earlier reference is shown to solve this problem. This method is improved upon by constructing a permutation of the rows and columns of L that preserves triangularity and allow for the best possible such partition. A number of practical examples and algorithmic details are presented. The parallelism attainable is illustrated by means of elimination trees and clique trees.

Alvarado, Fernando L.↗

Distributed memory compiler methods for irregular problems: Data copy reuse and runtime partitioning

Outlined here are two methods which we believe will play an important role in any distributed memory compiler able to handle sparse and unstructured problems. We describe how to link runtime partitioners to distributed memory compilers. In our scheme, programmers can implicitly specify how data and loop iterations are to be distributed between processors. This insulates users from having to deal explicitly with potentially complex algorithms that carry out work and data partitioning. We also describe a viable mechanism for tracking and reusing copies of off-processor data. In many programs, several loops access the same off-processor memory locations. As long as it can be verified that the values assigned to off-processor memory locations remain unmodified, we show that we can effectively reuse stored off-processor data. We present experimental data from a 3-D unstructured Euler solver run on iPSC/860 to demonstrate the usefulness of our methods.

Das, Raja↗

Nonequilibrium turbulence modeling effects on transonic vortical flows about delta wings

The Johnson-King turbulence model that is a viable method for calculating two dimensional transonic separated flows was extended into three dimensions. The implementation was done for Navier Stokes flow solvers written in general curvilinear coordinates. The present approach used in turbulence modeling is based on streamwise integration of an ordinary differential equation (o.d.e.) that governs the maximum Reynolds shear stress behavior. Streamwise integration of the o.d.e. approach was found to offer great mathematical simplicity and economy for three dimensional Navier Stokes methods. Thus, the new method is quick, simple, and very cheap. The new method was first checked against the data of a well known transonic axisymmetric bump experiment, and a good agreement was obtained. Later, the new method was used to compute the flow around a low aspect ratio wing in a transonic wind tunnel. Finally it was employed to study the nonequilibrium turbulence effects on the transonic vortical flows about a 65 deg sweep round leading edge delta wing.

Unver Kaynak↗

Navier-Stokes calculation of transonic flow past the NTF 65 deg delta wing

Models of four delta wings were built and tested in the 8 x 8 ft Transonic Wind Tunnel at LaRC. The wings are identical in planform shape with a swept-back angle of 65 degrees, but bear different leading edge profiles. The models were tested under pressurized and cryogenic conditions to simulate true flight Reynolds numbers. Data on the aerodynamic forces and pressure distributions at various locations on the wings were taken at various flight Mach numbers and angles of attack. Effects of high Reynolds number and leading edge radius on the aerodynamic characteristics of the wings are being accessed. To thoroughly understand the turbulent, vortical flow around the wings, an effort to perform computational aerodynamic analysis is being made. The objective of the analysis is to supplement and validate the experimental data and explain the high Reynolds number and leading edge effects. GRIDGEN, a software developed by General Dynamics, is being used to generate the grid topology for the flow field around the wings. The flow solver to be used is CFL3D, a computation fluid dynamics code developed at LaRC. Based on the geometric description of the wings, a FORTRAN program called WINGSURF was written to generate the databases defining the surface geometry of the wings. To match the true geometry of the models in the wind tunnel for realistic comparison with experimental results, databases defining the sting support for the wing models were also created by two other FORTRAN programs, STINJOIN and STINREAR. Listing of the programs are attached and the geometry of a typical wing model with the sting support is shown. Other aspects of the investigation are discussed.

Wu, Chivey C.↗

Application of direct inverse analogy method (DIVA) and viscous design optimization techniques

A direct-inverse approach to the transonic design problem was presented in its initial state at the First International Conference on Inverse Design Concepts and Optimization in Engineering Sciences (ICIDES-1). Further applications of the direct inverse analogy (DIVA) method to the design of airfoils and incremental wing improvements and experimental verification are reported. First results of a new viscous design code also from the residual correction type with semi-inverse boundary layer coupling are compared with DIVA which may enhance the accuracy of trailing edge design for highly loaded airfoils. Finally, the capabilities of an optimization routine coupled with the two viscous full potential solvers are investigated in comparison to the inverse method.

Greff, E.↗

An inverse method for the aerodynamic design of three-dimensional aircraft engine nacelles

A fast, efficient and user friendly inverse design system for 3-D nacelles was developed. The system is a product of a 2-D inverse design method originally developed at NASA-Langley and the CFL3D analysis code which was also developed at NASA-Langley and modified for nacelle analysis. The design system uses a predictor/corrector design approach in which an analysis code is used to calculate the flow field for an initial geometry, the geometry is then modified based on the difference between the calculated and target pressures. A detailed discussion of the design method, the process of linking it to the modified CFL3D solver and its extension to 3-D is presented. This is followed by a number of examples of the use of the design system for the design of both axisymmetric and 3-D nacelles.

Bell, R. A.↗

Airfoil optimization with efficient gradient calculations

The viscous airfoil design analysis code XFOIL was extended to allow optimization using conformal mapping coefficients as design variables. The optimization technique used was the Steepest Descent method applied to a Penalty Function. The gradients of the aerodynamic variables with respect to the design variables were cheaply calculated as by-products of XFOIL's integral boundary layer Newton solver. The speed of the optimization process further increased by updating the Newton system boundary layer variables after each optimization step using the available gradient information. Two examples are presented.

Sorensen, Thomas↗

Optimization of transistor design including large signal device/circuit interactions at extremely high frequencies (20-100+GHz)

Transistor design for extremely high frequency applications requires consideration of the interaction between the device and the circuit to which it is connected. Traditional analytical transistor models are to approximate at some of these frequencies and may not account for variations of dopants and semiconductor materials (especially some of the newer materials) within the device. Physically based models of device performance are required. These are based on coupled systems of partial differential equations and typically require 20 minutes of Cray computer time for a single AC operating point. A technique is presented to extract parameters from a few partial differential equation solutions for the device to create a nonlinear equivalent circuit model which runs in approximately 1 second of personal computer time. This nonlinear equivalent circuit model accurately replicates the contact current properties of the device as computed by the partial differential solver on which it is based. Using the nonlinear equivalent circuit model of the device, optimization of systems design can be performed based on device/circuit interactions.

Levy, Ralph↗