Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence analysis”

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 145 records · Page 8

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

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

Rasmussen, John↗

Convergence acceleration of an aeroelastic Navier-Stokes solver

New capabilities have been added to a Navier-Stokes solver to perform steady-state simulations more efficiently. The flow solver for solving the Navier-Stokes equations is completely rewritten with a combination of the LU-SGS (Lower-Upper factored Symmetric Gauss-Seidel) implicit method and the modified HLLE (Harten-Lax-van Leer-Einfeldt) upwind scheme. A pseudo-time marching method is used for the directly coupled structural equations to improve overall convergence rates for static aeroelastic analysis. Results are demonstrated for transonic flows over rigid and flexible wings.

Obayashi, S.↗

An Optimal Order Nonnested Mixed Multigrid Method for Generalized Stokes Problems

A multigrid algorithm is developed and analyzed for generalized Stokes problems discretized by various nonnested mixed finite elements within a unified framework. It is abstractly proved by an element-independent analysis that the multigrid algorithm converges with an optimal order if there exists a 'good' prolongation operator. A technique to construct a 'good' prolongation operator for nonnested multilevel finite element spaces is proposed. Its basic idea is to introduce a sequence of auxiliary nested multilevel finite element spaces and define a prolongation operator as a composite operator of two single grid level operators. This makes not only the construction of a prolongation operator much easier (the final explicit forms of such prolongation operators are fairly simple), but the verification of the approximate properties for prolongation operators is also simplified. Finally, as an application, the framework and technique is applied to seven typical nonnested mixed finite elements.

Deng, Qingping↗

A Non-Linear Approach to Spacecraft Trajectory Control in the Vicinity of a Libration Point

An expanding interest in mission design strategies that exploit libration point regions demands the continued development of enhanced, efficient, control algorithms for station-keeping and formation maintenance. This paper discusses the development of a non-linear, station-keeping, control algorithm for trajectories in the vicinity of a libration point. The control law guarantees exponential convergence, based on a Lyaponov analysis. Controller performance is evaluated using FreeFlyer(R) and MATLAB(R) for a spacecraft stationed near the L2 libration point in the Earth-Moon system, tracking a pre-defined reference trajectory. Evaluation metrics are fuel usage and tracking accuracy. Simulation results are compared with a linear-based controller for a spacecraft tracking the same reference trajectory. Although the analysis is framed in the context of station-keeping, the control algorithm is equally applicable to a formation flying problem with an appropriate definition of the reference trajectory.

Luquette, Richard J.↗

A Non-Linear Approach to Spacecraft Formation Control in the Vicinity of a Collinear Libration Point

An expanding interest in mission design strategies that exploit libration point regions demands the continued development of enhanced, efficient, control algorithms for station-keeping and formation maintenance. This paper discusses the development of a non-linear, formation maintenance, control algorithm for trajectories in the vicinity of a libration point. However, the formulation holds for any trajectory governed by the equations of motion for the restricted three body problem. The control law guarantees exponential convergence, based on a Lyaponov analysis. FreeFlyer and MATLAB provide the simulation environment for controller performance evaluation. The simulation, modeled after the MAXIM Pathfinder mission, maintains the relative position of a "follower" spacecraft with respect to a "leader" spacecraft, stationed near the L2 libration point in the Sun-Earth system. Evaluation metrics are fuel usage and tracking accuracy.

Luquette, Richard J.↗

Numerical Solution of Multi-Dimensional Hyperbolic Conservation Laws on Unstructured Meshes

The lecture material will discuss the application of one-dimensional approximate Riemann solutions and high order accurate data reconstruction as building blocks for solving multi-dimensional hyperbolic equations. This building block procedure is well-documented in the nationally available literature. The relevant stability and convergence theory using positive operator analysis will also be presented. All participants in the minisymposium will be asked to solve one or more generic test problems so that a critical comparison of accuracy can be made among differing approaches.

Barth, Timothy J.↗

Guidelines for Computing Longitudinal Dynamic Stability Characteristics of a Subsonic Transport

A systematic study is presented to guide the selection of a numerical solution strategy for URANS computation of a subsonic transport configuration undergoing simulated forced oscillation about its pitch axis. Forced oscillation is central to the prevalent wind tunnel methodology for quantifying aircraft dynamic stability derivatives from force and moment coefficients, which is the ultimate goal for the computational simulations. Extensive computations are performed that lead in key insights of the critical numerical parameters affecting solution convergence. A preliminary linear harmonic analysis is included to demonstrate the potential of extracting dynamic stability derivatives from computational solutions.

Thompson, Joseph R.↗

A Coupled Nonlinear Spacecraft Attitude Controller and Observer with an Unknown Gyro Misalignment and Gyro Bias

A nonlinear control scheme for attitude control of a spacecraft is combined with a nonlinear gyro misalignment and bias observer for the case of constant gyro misalignment and bias. A persistency of excitation analysis shows the observer gyro bias estimates converge to the true bias values exponentially fast. The convergence of the misalignment estimates is also presented. Then; the resulting coupled, closed loop dynamics are proven by a Lyapunov analysis to be globally stable, with asymptotically perfect tracking. The analysis is extended to consider the effects of noise in addition to the gyro misalignment and bias. A simulation of the proposed observer-controller design is given for a rigid spacecraft tracking a specified, time-varying attitude sequence to illustrate the theoretical claims.

Thienel, Julie↗

Some flow characteristics of conventional and tapered high pressure drop simulated seals

The leak rates through shaft seals with large pressure drops were simulated using gaseous hydrogen, or nitrogen flowing through an annulus with a nonrotating centerbody. The flows were choked. For concentric or eccentric position of the rotor and parallel or convergent tapered flow passages, data and analysis revealed that mass flux or leak rate can be determined from a relation whose normalizing parameters depend on the thermodynamic critical constants of the working fluid and an average flow area expressed in terms of the inlet and exit cross-sectional areas. Using these normalized relations, the flow data for parallel and three convergent tapered shaft seal configurations are in good agreement. Generalization to any simple gas or gas mixture is implied and demonstrated.

Hendricks, R. C.↗

Some flow characteristics of conventional and tapered high-pressure-drop simulated seals

The leak rates through shaft seals with large pressure drops were simulated using gaseous hydrogen, or nitrogen flowing through an annulus with a nonrotating centerbody. The flows were choked. For concentric or eccentric position of the rotor and parallel or convergent tapered flow passages, data and analysis revealed that mass flux or leak rate can be determined from a relation whose normalizing parameters depend on the thermodynamic critical constants of the working fluid and an average flow area expressed in terms of the inlet and exit cross-sectional areas. Using these normalized relations, the flow data for parallel and three convergent, tapered, shaft-seal configurations are in good agreement. Generalization to any simple gas or gas mixtures is implied and demonstrated in part.

Hendricks, R. C.↗

Finite element analysis of transonic flows in cascades: Importance of computational grids in improving accuracy and convergence

The finite element method is applied for the solution of transonic potential flows through a cascade of airfoils. Convergence characteristics of the solution scheme are discussed. Accuracy of the numerical solutions is investigated for various flow regions in the transonic flow configuration. The design of an efficient finite element computational grid is discussed for improving accuracy and convergence.

Ecer, A.↗

The multigrid preconditioned conjugate gradient method

A multigrid preconditioned conjugate gradient method (MGCG method), which uses the multigrid method as a preconditioner of the PCG method, is proposed. The multigrid method has inherent high parallelism and improves convergence of long wavelength components, which is important in iterative methods. By using this method as a preconditioner of the PCG method, an efficient method with high parallelism and fast convergence is obtained. First, it is considered a necessary condition of the multigrid preconditioner in order to satisfy requirements of a preconditioner of the PCG method. Next numerical experiments show a behavior of the MGCG method and that the MGCG method is superior to both the ICCG method and the multigrid method in point of fast convergence and high parallelism. This fast convergence is understood in terms of the eigenvalue analysis of the preconditioned matrix. From this observation of the multigrid preconditioner, it is realized that the MGCG method converges in very few iterations and the multigrid preconditioner is a desirable preconditioner of the conjugate gradient method.

Tatebe, Osamu↗

Material nonlinear analysis via mixed-iterative finite element method

The performance of elastic-plastic mixed-iterative analysis is examined through a set of convergence studies. Membrane and bending behaviors are tested using 4-node quadrilateral finite elements. The membrane result is excellent, which indicates the implementation of elastic-plastic mixed-iterative analysis is appropriate. On the other hand, further research to improve bending performance of the method seems to be warranted.

Sutjahjo, Edhi↗

Dual Extended Kalman Filter for the Identification of Time-Varying Human Manual Control Behavior

A Dual Extended Kalman Filter was implemented for the identification of time-varying human manual control behavior. Two filters that run concurrently were used, a state filter that estimates the equalization dynamics, and a parameter filter that estimates the neuromuscular parameters and time delay. Time-varying parameters were modeled as a random walk. The filter successfully estimated time-varying human control behavior in both simulated and experimental data. Simple guidelines are proposed for the tuning of the process and measurement covariance matrices and the initial parameter estimates. The tuning was performed on simulation data, and when applied on experimental data, only an increase in measurement process noise power was required in order for the filter to converge and estimate all parameters. A sensitivity analysis to initial parameter estimates showed that the filter is more sensitive to poor initial choices of neuromuscular parameters than equalization parameters, and bad choices for initial parameters can result in divergence, slow convergence, or parameter estimates that do not have a real physical interpretation. The promising results when applied to experimental data, together with its simple tuning and low dimension of the state-space, make the use of the Dual Extended Kalman Filter a viable option for identifying time-varying human control parameters in manual tracking tasks, which could be used in real-time human state monitoring and adaptive human-vehicle haptic interfaces.

manual control↗

On Bi-Grid Local Mode Analysis of Solution Techniques for 3-D Euler and Navier-Stokes Equations

A procedure is presented for utilizing a bi-grid stability analysis as a practical tool for predicting multigrid performance in a range of numerical methods for solving Euler and Navier-Stokes equations. Model problems based on the convection, diffusion and Burger's equation are used to illustrate the superiority of the bi-grid analysis as a predictive tool for multigrid performance in comparison to the smoothing factor derived from conventional von Neumann analysis. For the Euler equations, bi-grid analysis is presented for three upwind difference based factorizations, namely Spatial, Eigenvalue and Combination splits, and two central difference based factorizations, namely LU and ADI methods. In the former, both the Steger-Warming and van Leer flux-vector splitting methods are considered. For the Navier-Stokes equations, only the Beam-Warming (ADI) central difference scheme is considered. In each case, estimates of multigrid convergence rates from the bi-grid analysis are compared to smoothing factors obtained from single-grid stability analysis. Effects of grid aspect ratio and flow skewness are examined. Both predictions are compared with practical multigrid convergence rates for 2-D Euler and Navier-Stokes solutions based on the Beam-Warming central scheme.

Ibraheem, S. O.↗

Elastic-plastic mixed-iterative finite element analysis: Implementation and performance assessment

An elastic-plastic algorithm based on Von Mises and associative flow criteria is implemented in MHOST-a mixed iterative finite element analysis computer program developed by NASA Lewis Research Center. The performance of the resulting elastic-plastic mixed-iterative analysis is examined through a set of convergence studies. Membrane and bending behaviors of 4-node quadrilateral shell finite elements are tested for elastic-plastic performance. Generally, the membrane results are excellent, indicating the implementation of elastic-plastic mixed-iterative analysis is appropriate.

Sutjahjo, Edhi↗

Consistent linearization of the element-independent corotational formulation for the structural analysis of general shells

A consistent linearization is provided for the element-dependent corotational formulation, providing the proper first and second variation of the strain energy. As a result, the warping problem that has plagued flat elements has been overcome, with beneficial effects carried over to linear solutions. True Newton quadratic convergence has been restored to the Structural Analysis of General Shells (STAGS) code for conservative loading using the full corotational implementation. Some implications for general finite element analysis are discussed, including what effect the automatic frame invariance provided by this work might have on the development of new, improved elements.

Rankin, C. C.↗

Fourier analysis of the SOR iteration

The SOR iteration for solving linear systems of equations depends upon an overrelaxation factor omega. It is shown that for the standard model problem of Poisson's equation on a rectangle, the optimal omega and corresponding convergence rate can be rigorously obtained by Fourier analysis. The trick is to tilt the space-time grid so that the SOR stencil becomes symmetrical. The tilted grid also gives insight into the relation between convergence rates of several variants.

Leveque, R. J.↗