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 1,765 records · Page 98

A method of limit point calculation in finite element structural analysis

An approach is presented for the calculation of limit points for structures described by discrete coordinates, and whose governing equations derive from finite element concepts. The nonlinear load-displacement path of the imperfect structure is first traced by use of a direct iteration scheme and the determinant of the governing algebraic equations is calculated at each solution point. The limit point is then established by extrapolation and imposition of the condition of zero slope of the plot of load vs. determinant. Three problems are solved in illustration of the approach and in comparison with alternative procedures and test data.

Gallagher, R. H.↗

Application of Contraction Mappings to the Control of Nonlinear Systems

The theoretical and applied aspects of successive approximation techniques are considered for the determination of controls for nonlinear dynamical systems. Particular emphasis is placed upon the methods of contraction mappings and modified contraction mappings. It is shown that application of the Pontryagin principle to the optimal nonlinear regulator problem results in necessary conditions for optimality in the form of a two point boundary value problem (TPBVP). The TPBVP is represented by an operator equation and functional analytic results on the iterative solution of operator equations are applied. The general convergence theorems are translated and applied to those operators arising from the optimal regulation of nonlinear systems. It is shown that simply structured matrices and similarity transformations may be used to facilitate the calculation of the matrix Green functions and the evaluation of the convergence criteria. A controllability theory based on the integral representation of TPBVP's, the implicit function theorem, and contraction mappings is developed for nonlinear dynamical systems. Contraction mappings are theoretically and practically applied to a nonlinear control problem with bounded input control and the Lipschitz norm is used to prove convergence for the nondifferentiable operator. A dynamic model representing community drug usage is developed and the contraction mappings method is used to study the optimal regulation of the nonlinear system.

Killingsworth, W. R., Jr.↗

Methods for design and evaluation of integrated hardware-software systems for concurrent computation

Research activities and publications are briefly summarized. The major tasks reviewed are: (1) VAX implementation of the PISCES parallel programming environment; (2) Apollo workstation network implementation of the PISCES environment; (3) FLEX implementation of the PISCES environment; (4) sparse matrix iterative solver in PSICES Fortran; (5) image processing application of PISCES; and (6) a formal model of concurrent computation being developed.

Pratt, T. W.↗

Application of an advanced computerized structural design system to an arrow-wing supersonic cruise aircraft

A structural design study of an arrow-wing supersonic cruise aircraft has been made using the integrated design system, ATLAS, and a relatively large analytical finite-element model containing 8500 degrees of freedom. This paper focuses on structural design methods developed and used in support of the study with emphasis on aeroelasticity. The use of ATLAS permitted (1) automatic resizing of the wing structure for multiple load conditions, (2) rapid evaluation of aeroelastic effects, and (3) an iterative approach to the correction of flutter deficiencies. The significant results of the study are discussed along with the advantages derived from the use of an advanced structural design system in preliminary design studies.

Robinson, J. C.↗

Combustion instability with partial length acoustic liners

The development is discussed of analytical techniques and numerical methods for the prediction of the stability behavior of liquid propellant rocket combustors with partial length acoustic liners, injector face baffles, and nonlinear combustion responses. Three primary objectives included: (1) extension of the Green's Function integral-iteration technique previously developed so that nonlinear three dimensional wave propagation effects could be included in stability predictions; (2) development of an analytical technique and a computer program for the prediction of the damping effects caused by injector face baffles; and (3) inclusion of a nonlinear, wave distortation dependent, combustion response model in the global stability analysis and determination of its effect on combustor stability.

Mitchell, C. E.↗

Simplified contact analysis

A simple and straightforward method for calculating the elastic deformation at the center of contact between two solids with different radii of curvature is proposed instead of the conventional contact deformation analysis for ball bearings, gears, and cams, which usually involves tedious iterative procedures or the use of design charts. Ellipticity of the contact is approximated from a least-squares power fit. Simplified expressions are derived which allow rapid calculation of deformation from the material properties and geometry of the contacting elements.

Brewe, D. E.↗

Full potential integral solution for transonic flows with and without embedded Euler domains

Two methods are presented to solve for the transonic airfoil flow problems. The first method is based on the integral equation solution of the full-potential equation in terms of the velocity field, and a Shock Capturing-Shock Fitting (SCSF) scheme has been developed. The SCSF-scheme consists of a shock-capturing part and a shock-fitting part in which shock panels are introduced at the shock location. The sock panels are fitted and crossed by using the Rankine-Hugoniot relations. The second method is based on coupling the integral equation of the full-potential equation with the pseudo time integration of Euler equations in a small embedded region around the shock. The integral solution provides the initial and boundary conditions for the Euler domain. This scheme is named as the Integral Equation-Embedded Euler (IEEE) scheme. The two methods are applied to NACA 0012 and NACA 64A010A over a wide range of Mach numbers, and the reults are in good agreement with the experimental data and other computational results. The schemes converge within a number of iterations which is one-order of magnitude less than the finite-difference schemes.

Kandil, Osama A.↗

Thermostructural tailoring of fiber composite structures

A significant area of interest in design of complex structures involves the study of multidisciplined problems. The coordination of several different intricate areas of study to obtain a particular design of a structure is a new and pressing area of research. In the past, each discipline would perform its task consecutively using the appropriate inputs from the other disciplines. This process usually required several time-consuming iterations to obtain a satisfactory design. The alternative pursued here is combining various participating disciplines and specified design requirements into a formal structural computer code. The main focus of this research is to develop a multidiscipline structural tailoring method for select composite structures and to demonstrate its application to specific areas. The development of an integrated computer program involves the coupling of three independent computer programs using an excutive module. This module will be the foundation for integrating a structural optimizer, a composites analyzer and a thermal analyzer. With the completion of the executive module, the first step was taken toward the evolution of multidiscipline software in the field of composite mechanics. Through the use of an array of cases involving a variety of objective functions/constraints and thermal-mechanical load conditions, it became evident that simple composite structures can be designed to a combined loads environment.

Acquaviva, Thomas H.↗

A Parallel Genetic Algorithm for Automated Electronic Circuit Design

Parallelized versions of genetic algorithms (GAs) are popular primarily for three reasons: the GA is an inherently parallel algorithm, typical GA applications are very compute intensive, and powerful computing platforms, especially Beowulf-style computing clusters, are becoming more affordable and easier to implement. In addition, the low communication bandwidth required allows the use of inexpensive networking hardware such as standard office ethernet. In this paper we describe a parallel GA and its use in automated high-level circuit design. Genetic algorithms are a type of trial-and-error search technique that are guided by principles of Darwinian evolution. Just as the genetic material of two living organisms can intermix to produce offspring that are better adapted to their environment, GAs expose genetic material, frequently strings of 1s and Os, to the forces of artificial evolution: selection, mutation, recombination, etc. GAs start with a pool of randomly-generated candidate solutions which are then tested and scored with respect to their utility. Solutions are then bred by probabilistically selecting high quality parents and recombining their genetic representations to produce offspring solutions. Offspring are typically subjected to a small amount of random mutation. After a pool of offspring is produced, this process iterates until a satisfactory solution is found or an iteration limit is reached. Genetic algorithms have been applied to a wide variety of problems in many fields, including chemistry, biology, and many engineering disciplines. There are many styles of parallelism used in implementing parallel GAs. One such method is called the master-slave or processor farm approach. In this technique, slave nodes are used solely to compute fitness evaluations (the most time consuming part). The master processor collects fitness scores from the nodes and performs the genetic operators (selection, reproduction, variation, etc.). Because of dependency issues in the GA, it is possible to have idle processors. However, as long as the load at each processing node is similar, the processors are kept busy nearly all of the time. In applying GAs to circuit design, a suitable genetic representation 'is that of a circuit-construction program. We discuss one such circuit-construction programming language and show how evolution can generate useful analog circuit designs. This language has the desirable property that virtually all sets of combinations of primitives result in valid circuit graphs. Our system allows circuit size (number of devices), circuit topology, and device values to be evolved. Using a parallel genetic algorithm and circuit simulation software, we present experimental results as applied to three analog filter and two amplifier design tasks. For example, a figure shows an 85 dB amplifier design evolved by our system, and another figure shows the performance of that circuit (gain and frequency response). In all tasks, our system is able to generate circuits that achieve the target specifications.

Long, Jason D.↗

Automating CPM-GOMS

CPM-GOMS is a modeling method that combines the task decomposition of a GOMS analysis with a model of human resource usage at the level of cognitive, perceptual, and motor operations. CPM-GOMS models have made accurate predictions about skilled user behavior in routine tasks, but developing such models is tedious and error-prone. We describe a process for automatically generating CPM-GOMS models from a hierarchical task decomposition expressed in a cognitive modeling tool called Apex. Resource scheduling in Apex automates the difficult task of interleaving the cognitive, perceptual, and motor resources underlying common task operators (e.g. mouse move-and-click). Apex's UI automatically generates PERT charts, which allow modelers to visualize a model's complex parallel behavior. Because interleaving and visualization is now automated, it is feasible to construct arbitrarily long sequences of behavior. To demonstrate the process, we present a model of automated teller interactions in Apex and discuss implications for user modeling. available to model human users, the Goals, Operators, Methods, and Selection (GOMS) method [6, 21] has been the most widely used, providing accurate, often zero-parameter, predictions of the routine performance of skilled users in a wide range of procedural tasks [6, 13, 15, 27, 28]. GOMS is meant to model routine behavior. The user is assumed to have methods that apply sequences of operators and to achieve a goal. Selection rules are applied when there is more than one method to achieve a goal. Many routine tasks lend themselves well to such decomposition. Decomposition produces a representation of the task as a set of nested goal states that include an initial state and a final state. The iterative decomposition into goals and nested subgoals can terminate in primitives of any desired granularity, the choice of level of detail dependent on the predictions required. Although GOMS has proven useful in HCI, tools to support the construction of GOMS models have not yet come into general use.

GOMS↗

A Method for Combining Experimentation and Molecular Dynamics Simulation to Improve Cohesive Zone Models for Metallic Microstructures

Fracture processes within a material begin at the nanometer length scale at which the formation, propagation, and interaction of fundamental damage mechanisms occur. Physics-based modeling of these atomic processes quickly becomes computationally intractable as the system size increases. Thus, a multiscale modeling method, based on the aggregation of fundamental damage processes occurring at the nanoscale within a cohesive zone model, is under development and will enable computationally feasible and physically meaningful microscale fracture simulation in polycrystalline metals. This method employs atomistic simulation to provide an optimization loop with an initial prediction of a cohesive zone model (CZM). This initial CZM is then applied at the crack front region within a finite element model. The optimization procedure iterates upon the CZM until the finite element model acceptably reproduces the near-crack-front displacement fields obtained from experimental observation. With this approach, a comparison can be made between the original CZM predicted by atomistic simulation and the converged CZM that is based on experimental observation. Comparison of the two CZMs gives insight into how atomistic simulation scales.

Hochhalter, J. D.↗

Squash-Box Feasibility Driven Differential Dynamic Programming

Recently, Differential Dynamic Programming (DDP) and other similar algorithms have become the solvers of choice when performing non-linear Model Predictive Control (nMPC) with modern robotic devices. The reason is that they have a lower computational cost per iteration when compared with off-the-shelf Non-Linear Programming (NLP) solvers, which enables its online operation. However, they cannot handle constraints, and are known to have poor convergence capabilities. In this paper, we propose a method to solve the optimal control problem with control bounds through a squashing function (i.e., a sigmoid, which is bounded by construction). It has been shown that a naive use of squashing functions damage the convergence rate. To tackle this, we first propose to add a quadratic barrier that avoids the difficulty of the plateau produced by the sigmoid. Second, we add an outer loop that adapts both the sigmoid and the barrier; it makes the optimal control problem with the squashing function converge to the original control-bounded problem. To validate our method, we present simulation results for different types of platforms including a multi-rotor, a biped, a quadruped and a humanoid robot.

Navarro, Angel Santamaria↗

Methods for design and evaluation of integrated hardware/software systems for concurrent computation

Two testbed programming environments to support the evaluation of a large range of parallel architectures have been implemented under the program Parallel Implementation of Scientific Computing Environments (PISCES). The PISCES 1 environment was applied to two areas of aerospace interest: a sparse matrix iterative equation solver and a dynamic scene analysis system. Currently, the NICE/SPAR testbed system for structural analysis is being modified for parallel operation under PISCES 2; the PISCES 1 applications are also being adapted for PISCES 2. A new formal model of concurrent computation has been developed, based on the mathematical system known as H graph semantics together with a timed Petri net model of the parallel aspects of a system.

Pratt, Terrence W.↗

Canonical-variables multigrid method for steady-state Euler equations

A novel approach for the solution of inviscid flow problems for subsonic compressible flows is described. The approach is based on canonical forms of the equations, in which subsystems governed by hyperbolic operators are separated from those governed by elliptic ones. The discretizations used as well as the iterative techniques for the different subsystems are inherently different. Hyperbolic parts, which describe, in general, propagation phenomena, are discretized using upwind schemes and are solved by marching techniques. Elliptic parts, which are directionally unbiased, are discretized using h-elliptic central discretizations, and are solved by pointwise relaxations together with coarse grid acceleration. The resulting discretization schemes introduce artificial viscosity only for the hyperbolic parts of the system; thus a smaller total artificial viscosity is used, while the multigrid solvers used are much more efficient. Solutions of the subsonic compressible Euler equations are achieved at the same efficiency as the full potential equation.

Taasan, Shlomo↗

A two dimensional model of viscous slot flow

The effects of viscosity on flows through slots are investigated. Using Briley's formulation of the Alternating Direction Implicit method together with the Thompson transformation, numerical solution of viscous incompressible flow through a slot is simulated, for Reynolds number 1000 and openness ratio one tenth. The problem is formulated in terms of stream function and vorticity. A wall iteration of Israeli type is employed to iteratively induce the no slip condition on velocity.

Cooke, C. H.↗

A theoretical study of a two-wavelength lidar technique for the measurement of atmospheric temperature profiles

The theory of differential absorption lidar measurements for lines with a Voigt profile is given and applied to a two-wavelength technique for measuring the atmospheric temperature profile using a high J line in the oxygen A band. Explicit expressions for the temperature and pressure dependence of the absorption coefficient are developed for lines with a Voigt profile. An iteration procedure for calculating the temperature for narrow laser bandwidths is described which has an accuracy better than 0.2 K for bandwidths less than 0.01/cm. To reduce the errors in lidar measurements due to uncertainties in pressure, a method for estimating the pressure from the temperature profile is described. A procedure for extending the differential absorption technique to the case of finite laser bandwidth with good accuracy is also described. Simulation results show that a knowledge of the laser frequency is needed to 0.005/cm for accurate temperature measurements. Evaluation of the sensitivity for both ground- and Shuttle-based measurements shows accuracies generally better than 1 K. This technique allows up to an order of magnitude improvement in sensitivity compared to other differential absorption lidar techniques.

Korb, C. L.↗

Primitive numerical simulation of circular Couette flow - Carrousel wind tunnel nonturbulent solutions

The azimuthal-invariant, three-dimensional cylindrical, incompressible Navier-Stokes equations are solved to steady state for a finite-length, physically realistic model. The numerical method relies on an alternating-direction implicit scheme that is formally second-order accurate in space and first-order accurate in time. The equations are linearized and uncoupled by evaluating variable coefficients at the previous time iteration. Wall grid clustering is provided by a Roberts transformation in radial and axial directions. A vorticity-velocity formulation is found to be preferable to a vorticity-streamfunction approach. Subject to no-slip, Dirichlet boundary conditions, except for the inner cylinder rotation velocity (impulsive start-up) and zero-flow initial conditions, nonturbulent solutions are obtained for sub- and supercritical Reynolds numbers of 100 to 400 for a finite geometry where R(outer)/R(inner) = 1.5, H/R(inner) = 0.73, and H/Delta-R = 1.5. An axially-stretched model solution is shown to asymptotically approach the one-dimensional analytic Couette solution at the cylinder midheight. Flowfield change from laminar to Taylor-vortex flow is discussed as a function of Reynolds number. Three-dimensional velocities, vorticity, and streamfunction are presented via two-dimensional graphs and three-dimensional surface and contour plots.

Hasiuk, Jan↗

On 3D inelastic analysis methods for hot section components

Accomplishments are described for the 2-year program, to develop advanced 3-D inelastic structural stress analysis methods and solution strategies for more accurate and cost effective analysis of combustors, turbine blades and vanes. The approach was to develop a matrix of formulation elements and constitutive models. Three constitutive models were developed in conjunction with optimized iterating techniques, accelerators, and convergence criteria within a framework of dynamic time incrementing. Three formulations models were developed; an eight-noded mid-surface shell element, a nine-noded mid-surface shell element and a twenty-noded isoparametric solid element. A separate computer program was developed for each combination of constitutive model-formulation model. Each program provides a functional stand alone capability for performing cyclic nonlinear structural analysis. In addition, the analysis capabilities incorporated into each program can be abstracted in subroutine form for incorporation into other codes or to form new combinations.

Mcknight, R. L.↗