Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “precondition”

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 289 records · Page 16

Discrete sensitivity derivatives of the Navier-Stokes equations with a parallel Krylov solver

This paper solves an 'incremental' form of the sensitivity equations derived by differentiating the discretized thin-layer Navier Stokes equations with respect to certain design variables of interest. The equations are solved with a parallel, preconditioned Generalized Minimal RESidual (GMRES) solver on a distributed-memory architecture. The 'serial' sensitivity analysis code is parallelized by using the Single Program Multiple Data (SPMD) programming model, domain decomposition techniques, and message-passing tools. Sensitivity derivatives are computed for low and high Reynolds number flows over a NACA 1406 airfoil on a 32-processor Intel Hypercube, and found to be identical to those computed on a single-processor Cray Y-MP. It is estimated that the parallel sensitivity analysis code has to be run on 40-50 processors of the Intel Hypercube in order to match the single-processor processing time of a Cray Y-MP.

Ajmani, Kumud↗

Parallel implicit unstructured grid Euler solvers

A mesh-vertex finite volume scheme for solving the Euler equations on triangular unstructured meshes is implemented on an MIMD (multiple instruction/multiple data stream) parallel computer. An explicit four-stage Runge-Kutta scheme is used to solve two-dimensional flow problems. A family of implicit schemes is also developed to solve these problems, where the linear system that arises at each time step is solved by a preconditioned GMRES algorithm. Two partitioning strategies are employed, one that partitions triangles and the other that partitions vertices. The choice of the preconditioner in a distributed memory setting is discussed. All the methods are compared both in terms of elapsed times and convergence rates. It is shown that the implicit schemes offer adequate parallelism at the expense of minimal sequential overhead. The use of a global coarse grid to further minimize this overhead is also investigated. The schemes are implemented on a distributed memory parallel computer, the iPSC/860.

Venkatakrishnan, V.↗

CFD analyses for advanced pump design

As one of the activities of the NASA/MSFC Pump Stage Technology Team, the present effort was focused on using CFD in the design and analysis of high performance rocket engine pumps. Under this effort, a three-dimensional Navier-Stokes code was used for various inducer and impeller flow field calculations. An existing algebraic grid generation procedure was-extended to allow for nonzero blade thickness, splitter blades, and hub/shroud cavities upstream or downstream of the (main) blades. This resulted in a fast, robust inducer/impeller geometry/grid generation package. Problems associated with running a compressible flow code to simulate an incompressible flow were resolved; related aspects of the numerical algorithm (viz., the matrix preconditioning, the artificial dissipation, and the treatment of low Mach number flows) were addressed. As shown by the calculations performed under the present effort, the resulting code, in conjunction with the grid generation package, is an effective tool for the rapid solution of three-dimensional viscous inducer and impeller flows.

Dejong, F. J.↗

Process modelling for materials preparation experiments

The main goals of the research under this grant consist of the development of mathematical tools and measurement techniques for transport properties necessary for high fidelity modelling of crystal growth from the melt and solution. Of the tasks described in detail in the original proposal, two remain to be worked on: development of a spectral code for moving boundary problems, and development of an expedient diffusivity measurement technique for concentrated and supersaturated solutions. We have focused on developing a code to solve for interface shape, heat and species transport during directional solidification. The work involved the computation of heat, mass and momentum transfer during Bridgman-Stockbarger solidification of compound semiconductors. Domain decomposition techniques and preconditioning methods were used in conjunction with Chebyshev spectral methods to accelerate convergence while retaining the high-order spectral accuracy. During the report period we have further improved our experimental setup. These improvements include: temperature control of the measurement cell to 0.1 C between 10 and 60 C; enclosure of the optical measurement path outside the ZYGO interferometer in a metal housing that is temperature controlled to the same temperature setting as the measurement cell; simultaneous dispensing and partial removal of the lower concentration (lighter) solution above the higher concentration (heavier) solution through independently motor-driven syringes; three-fold increase in data resolution by orientation of the interferometer with respect to diffusion direction; and increase of the optical path length in the solution cell to 12 mm.

Rosenberger, Franz↗

Thermal and hydraulic considerations regarding the fate of water discharged by the outflow channels to the Martian northern plains

The identification of possible shorelines in the Martian northern plains suggests that the water discharged by the circum-Chryse outflow channels may have led to the formation of transient seas, or possibly even an ocean, covering as much as one-third of the planet. Speculations regarding the possible fate of this water have included local ponding and reinfiltration into the crust; freezing, sublimation, and eventual cold-trapping at higher latitudes; or the in situ survival of this now frozen water to the present day -- perhaps aided by burial beneath a protective cover of eolian sediment or lavas. Although neither cold-trapping at higher latitudes nor the subsequent freezing and burial of flood waters can be ruled out, thermal and hydraulic considerations effectively eliminate the possibility that any significant reassimilation of this water by local infiltration has occurred given climatic conditions resembling those of today. The arguments against the local infiltration of flood water into the northern plains are two-fold. First, given the climatic and geothermal conditions that are thought to have prevailed on Mars during the Late Hesperian (the period of peak outflow channel activity in the northern plains), the thickness of the cryosphere in Chryse Planitia is likely to have exceeded 1 km. A necessary precondition for the widespread occurrence of groundwater is that the thermodynamic sink represented by the cryosphere must already be saturated with ice. For this reason, the ice-saturated cryosphere acts as an impermeable barrier that effectively precludes the local resupply of subpermafrost groundwater by the infiltration of water discharged to the surface by catastraphic floods. Note that the problem of local infiltration is not significantly improved even if the cryosphere were initially dry, for as water attempts to infiltrate the cold, dry crust, it will quickly freeze, creating a seal that prevents any further infiltration from the ponded water above. The second argument against the local infiltration of flood water in the northern plains is based on hydraulic considerations. Repeated impacts have likely brecciated the Martian crust down to a depth of roughly 10 km. Given a value of permeability no greater than that inferred for the top 10 km of the Earth's crust (approximately 10(exp -2) darcies), a timescale as much as a billion years or more for the Martian groundwater system to achieve hydrostatic equilibrium, and the approximately 2-4 km elevation difference between the outflow channel source regions and the northern plains, the water confined beneath the frozen crust of the northern plains should have been under a significant hydraulic head. Thus, the existence of a hydraulic pathway between the ponded flood waters above the northern plains and the confined aquifer lying beneath it would not have led to the infiltration of flood water back into the crust, but rather the additional expulsion of groundwater onto the surface.

Clifford, S. M.↗

The least-squares finite element method for low-mach-number compressible viscous flows

The present paper reports the development of the Least-Squares Finite Element Method (LSFEM) for simulating compressible viscous flows at low Mach numbers in which the incompressible flows pose as an extreme. Conventional approach requires special treatments for low-speed flows calculations: finite difference and finite volume methods are based on the use of the staggered grid or the preconditioning technique; and, finite element methods rely on the mixed method and the operator-splitting method. In this paper, however, we show that such difficulty does not exist for the LSFEM and no special treatment is needed. The LSFEM always leads to a symmetric, positive-definite matrix through which the compressible flow equations can be effectively solved. Two numerical examples are included to demonstrate the method: first, driven cavity flows at various Reynolds numbers; and, buoyancy-driven flows with significant density variation. Both examples are calculated by using full compressible flow equations.

Yu, Sheng-Tao↗

Multigrid one shot methods for optimal control problems: Infinite dimensional control

The multigrid one shot method for optimal control problems, governed by elliptic systems, is introduced for the infinite dimensional control space. ln this case, the control variable is a function whose discrete representation involves_an increasing number of variables with grid refinement. The minimization algorithm uses Lagrange multipliers to calculate sensitivity gradients. A preconditioned gradient descent algorithm is accelerated by a set of coarse grids. It optimizes for different scales in the representation of the control variable on different discretization levels. An analysis which reduces the problem to the boundary is introduced. It is used to approximate the two level asymptotic convergence rate, to determine the amplitude of the minimization steps, and the choice of a high pass filter to be used when necessary. The effectiveness of the method is demonstrated on a series of test problems. The new method enables the solutions of optimal control problems at the same cost of solving the corresponding analysis problems just a few times.

Arian, Eyal↗

Integrating deliberative planning in a robot architecture

The role of planning and reactive control in an architecture for autonomous agents is discussed. The postulated architecture seperates the general robot intelligence problem into three interacting pieces: (1) robot reactive skills, i.e., grasping, object tracking, etc.; (2) a sequencing capability to differentially ativate the reactive skills; and (3) a delibrative planning capability to reason in depth about goals, preconditions, resources, and timing constraints. Within the sequencing module, caching techniques are used for handling routine activities. The planning system then builds on these cached solutions to routine tasks to build larger grain sized primitives. This eliminates large numbers of essentially linear planning problems. The architecture will be used in the future to incorporate in robots cognitive capabilites normally associated with intelligent behavior.

Elsaesser, Chris↗

Investigation of ramp injectors for supersonic mixing enhancement

A comparative study of wall mounted swept ramp injectors fitted with injector nozzles of different shape has been conducted in a constant area duct to explore mixing enhancement techniques for scramjet combustors. Six different injector nozzle inserts, all having equal exit and throat areas, were tested to explore the interaction between the preconditioned fuel jet and the vortical flowfield produced by the ramp: circular nozzle (baseline), nozzle with three downstream facing steps, nozzle with four vortex generators, elliptical nozzle, tapered-slot nozzle, and trapezoidal nozzle. The main flow was air at Mach 2, and the fuel was simulated by air injected at Mach 1.63 or by helium injected at Mach 1.7. Pressure and temperature surveys, combined with Mie and Rayleigh scattering visualization, were used to investigate the flow field. The experiments were compared with three dimensional Navier-Stokes computations. The results indicate that the mixing process is dominated by the streamwise vorticity generated by the ramp, the injectors' inner geometry having a minor effect. It was also found that the injectant/air mixing in the far-field is nearly independent of the injector geometry, molecular weight of the injectant, and the initial convective Mach number.

Haimovitch, Y.↗

Characterization of path effects on the properties of balloon films

The effect of stress and temperature path on the mechanical properties of a balloon film material (Astrofilm E2) have been studied using two techniques. The first technique is based on uniaxial tensile measurements of prestrained strips of the film as a function of temperature. Data are analyzed in terms of two temperature coefficients, defined here as a stiffness coefficient and a strength coefficient; and by comparing the areas under the stress-strain curves. The second technique investigated the effect of stress and temperature path on inflated cylinders of the same material. The changes in the mechanical properties due to preconditioning in the biaxial stress state were evaluated using the ball burst test. Preliminary findings indicate that the material is stress and temperature path dependent. Changes in stiffness coefficient, draw strength and toughness were measured and discussed. The above techniques may be suitable to document path dependent changes, and to discriminate among films of different manufacture.

Said, Magdi A.↗

Parallel implicit unstructured grid Euler solvers

A mesh-vertex finite volume scheme for solving the Euler equations on triangular unstructured meshes is implemented on a multiple-instruction/multiple-data stream parallel computer. An explicit four-stage Runge-Kutta scheme is used to solve two-dimensional flow problems. A family of implicit schemes is also developed to solve these problems, where the linear system that arises at each time step is solved by a preconditioned GMRES algorithm. Two partitioning strategies are employed: one that partitions triangles and the other that partitions vertices. The choice of the preconditioner in a distributed memory setting is discussed. All of the methods are compared both in terms of elapsed times and convergence rates. It is shown that the implicit schemes offer adequate parallelism at the expense of minimal sequential overhead. The use of a global coarse grid to further minimize this overhead is also investigated. The schemes are implemented on a distributed memory parallel computer, the Intel iPSC/860.

TRT-THEORETICAL↗

Fast structural design and analysis via hybrid domain decomposition on massively parallel processors

A hybrid domain decomposition framework for static, transient and eigen finite element analyses of structural mechanics problems is presented. Its basic ingredients include physical substructuring and /or automatic mesh partitioning, mapping algorithms, 'gluing' approximations for fast design modifications and evaluations, and fast direct and preconditioned iterative solvers for local and interface subproblems. The overall methodology is illustrated with the structural design of a solar viewing payload that is scheduled to fly in March 1993. This payload has been entirely designed and validated by a group of undergraduate students at the University of Colorado using the proposed hybrid domain decomposition approach on a massively parallel processor. Performance results are reported on the CRAY Y-MP/8 and the iPSC-860/64 Touchstone systems, which represent both extreme parallel architectures. The hybrid domain decomposition methodology is shown to outperform leading solution algorithms and to exhibit an excellent parallel scalability.

Farhat, Charbel↗

Cognitive factors associated with immersion in virtual environments

Immersion into the dataspace provided by a computer, and the feeling of really being there or 'presence', are commonly acknowledged as the uniquely important features of virtual reality environments. How immersed one feels appears to be determined by a complex set of physical components and affordances of the environment, and as yet poorly understood psychological processes. Pimentel and Teixeira say that the experience of being immersed in a computer-generated world involves the same mental shift of 'suspending your disbelief for a period of time' as 'when you get wrapped up in a good novel or become absorbed in playing a computer game'. That sounds as if it could be right, but it would be good to get some evidence for these important conclusions. It might be even better to try to connect these statements with theoretical positions that try to do justice to complex cognitive processes. The basic precondition for understanding Virtual Reality (VR) is understanding the spatial representation systems that localize our bodies or egocenters in space. The effort to understand these cognitive processes is being driven with new energy by the pragmatic demands of successful virtual reality environments, but the literature is largely sparse and anecdotal.

Psotka, Joseph↗

Three-dimensional aerodynamic shape optimization using discrete sensitivity analysis

An aerodynamic shape optimization procedure based on discrete sensitivity analysis is extended to treat three-dimensional geometries. The function of sensitivity analysis is to directly couple computational fluid dynamics (CFD) with numerical optimization techniques, which facilitates the construction of efficient direct-design methods. The development of a practical three-dimensional design procedures entails many challenges, such as: (1) the demand for significant efficiency improvements over current design methods; (2) a general and flexible three-dimensional surface representation; and (3) the efficient solution of very large systems of linear algebraic equations. It is demonstrated that each of these challenges is overcome by: (1) employing fully implicit (Newton) methods for the CFD analyses; (2) adopting a Bezier-Bernstein polynomial parameterization of two- and three-dimensional surfaces; and (3) using preconditioned conjugate gradient-like linear system solvers. Whereas each of these extensions independently yields an improvement in computational efficiency, the combined effect of implementing all the extensions simultaneously results in a significant factor of 50 decrease in computational time and a factor of eight reduction in memory over the most efficient design strategies in current use. The new aerodynamic shape optimization procedure is demonstrated in the design of both two- and three-dimensional inviscid aerodynamic problems including a two-dimensional supersonic internal/external nozzle, two-dimensional transonic airfoils (resulting in supercritical shapes), three-dimensional transport wings, and three-dimensional supersonic delta wings. Each design application results in realistic and useful optimized shapes.

Burgreen, Gregory W.↗

Design of 3-D Nacelle near Flat-Plate Wing Using Multiblock Sensitivity Analysis (ADOS)

One of the major design tasks involved in reducing aircraft drag is the integration of the engine nacelles and airframe. With this impetus, nacelle shapes with and without the presence of a flat-plate wing nearby were optimized. This also served as a demonstration of the 3-D version of the recently developed aerodynamic design optimization methodology using sensitivity analysis, ADOS. The required flow analyses were obtained by solving the three-dimensional, compressible, thin-layer Navier-Stokes equations using an implicit, upwind-biased, finite volume scheme. The sensitivity analyses were performed using the preconditioned version of the SADD scheme (sensitivity analysis on domain decomposition). In addition to demonstrating the present method's capability for automatic optimization, the results offered some insight into two important issues related to optimizing the shapes of multicomponent configurations in close proximity. First, inclusion of the mutual interference between the components resulted in a different shape as opposed to shaping an isolated component. Secondly, exclusion of the viscous effects compromised not only the flow physics but also the optimized shapes even for isolated components.

Eleshaky, Mohamed E.↗

Coupled ice-ocean variability in the Greenland Sea

A major surface feature of the Greenland Sea during winter is the frequent eastward extension of sea ice south of 75 deg N and an associated embayment to the north. These features are nominally connected with the East Greenland Current, and both the promontory) and the embayment are readily apparent on climatic ice charts. However, there are significant changes in these features on time-scales as short as a few days. Using a combination of satellite microwave images (SSM/I) of ice cover, meteorological data and in situ velocity, temperature and salinity records, we relate the ice distribution and its changes to the developing structure and circulation of the upper ocean during winter 1988-1989. Our measurements illustrate the preconditioning that leads to convective overturn, which in turn brings warmer water to the surface and results in the rapid disappearance of ice. In particular, the surface was cooled to the freezing point by early December and the salinity then increased through ice formation (about 0.016 m/d) and brine rejection. Once the vertical density gradient was sufficiently eroded, a period of high heat flux (greater than 300 W/sq m) in late January provided enough buoyancy loss to convectively mix the upper water column to at least 200 m. We estimate vertical velocities at about 3 cm/s downward during the initial sinking. The deepening of the thermocline raised surface temperatures by over 1 C resulting in nearly 1.5 x 105 km of ice-melt within two days. Average rates of ice retreat are about 2 km/d southwestward, generally consistent with a wind-driven flow. Comparison of hydrographic surveys from before and after the overturning indicate the fresh water was advected out of the area, possibly to the south and east of our moorings.

Roach, A. T.↗

Numerical Modeling of Spray Combustion with an Unstructured-Grid Method

The present unstructured-grid method follows strictly the basic finite volume forms of the conservation laws of the governing equations for the entire flow domain. High-order spatially accurate formulation has been employed for the numerical solutions of the Navier-Stokes equations. A two-equation k-epsilon turbulence model is also incorporated in the unstructured-grid solver. The convergence of the resulted linear algebraic equation is accelerated with preconditioned Conjugate Gradient method. A statistical spray combustion model has been incorporated into the present unstructured-grid solver. In this model, spray is represented by discrete particles, rather than by continuous distributions. A finite number of computational particles are used to predict a sample of total population of particles. Particle trajectories are integrated using their momentum and motion equations and particles exchange mass, momentum and energy with the gas within the computational cell in which they are located. The interaction calculations are performed simultaneously and eliminate global iteration for the two-phase momentum exchange. A transient spray flame in a high pressure combustion chamber is predicted and then the solution of liquid-fuel combusting flow with a rotating cup atomizer is presented and compared with the experimental data. The major conclusion of this investigation is that the unstructured-grid method can be employed to study very complicated flow fields of turbulent spray combustion. Grid adaptation can be easily achieved in any flow domain such as droplet evaporation and combustion zone. Future applications of the present model can be found in the full three-dimensional study of flow fields of gas turbine and liquid propulsion engine combustion chambers with multi-injectors.

Shang, H. M.↗

An Integrated Planning Representation Using Macros, Abstractions, and Cases

Planning will be an essential part of future autonomous robots and integrated intelligent systems. This paper focuses on learning problem solving knowledge in planning systems. The system is based on a common representation for macros, abstractions, and cases. Therefore, it is able to exploit both classical and case based techniques. The general operators in a successful plan derivation would be assessed for their potential usefulness, and some stored. The feasibility of this approach was studied through the implementation of a learning system for abstraction. New macros are motivated by trying to improve the operatorset. One heuristic used to improve the operator set is generating operators with more general preconditions than existing ones. This heuristic leads naturally to abstraction hierarchies. This investigation showed promising results on the towers of Hanoi problem. The paper concludes by describing methods for learning other problem solving knowledge. This knowledge can be represented by allowing operators at different levels of abstraction in a refinement.

Baltes, Jacky↗