Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear systems”

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 793 records · Page 44

Discretized Newton-relaxation solution of high resolution flux-difference split schemes

A method is presented for the implicit solution of high resolution flux-difference split schemes for which it is impractical to obtain the Jacobian matrix analytically. The method is used frequently in numerical schemes for solving nonlinear systems of equations, but apparently has not been used to date for approximating the Jacobian of the Roe numerical flux vectors. The method consists of obtaining a discretized Jacobian by using simple finite-differences of the flux vectors, and then using this discretized Jacobian in a Newton-relaxation solution method. The method is applied to both compressible and incompressible flows. Numerical results are presented that demonstrate the method works for both inviscid and viscous flows.

Whitfield, David L.↗

Experiments with conjugate gradient algorithms for homotopy curve tracking

There are algorithms for finding zeros or fixed points of nonlinear systems of equations that are globally convergent for almost all starting points, i.e., with probability one. The essence of all such algorithms is the construction of an appropriate homotopy map and then tracking some smooth curve in the zero set of this homotopy map. HOMPACK is a mathematical software package implementing globally convergent homotopy algorithms with three different techniques for tracking a homotopy zero curve, and has separate routines for dense and sparse Jacobian matrices. The HOMPACK algorithms for sparse Jacobian matrices use a preconditioned conjugate gradient algorithm for the computation of the kernel of the homotopy Jacobian matrix, a required linear algebra step for homotopy curve tracking. Here, variants of the conjugate gradient algorithm are implemented in the context of homotopy curve tracking and compared with Craig's preconditioned conjugate gradient method used in HOMPACK. The test problems used include actual large scale, sparse structural mechanics problems.

Irani, Kashmira M.↗

Resonant frequency calculations using a hybrid perturbation-Galerkin technique

A two-step hybrid perturbation Galerkin technique is applied to the problem of determining the resonant frequencies of one or several degree of freedom nonlinear systems involving a parameter. In one step, the Lindstedt-Poincare method is used to determine perturbation solutions which are formally valid about one or more special values of the parameter (e.g., for large or small values of the parameter). In step two, a subset of the perturbation coordinate functions determined in step one is used in Galerkin type approximation. The technique is illustrated for several one degree of freedom systems, including the Duffing and van der Pol oscillators, as well as for the compound pendulum. For all of the examples considered, it is shown that the frequencies obtained by the hybrid technique using only a few terms from the perturbation solutions are significantly more accurate than the perturbation results on which they are based, and they compare very well with frequencies obtained by purely numerical methods.

Geer, James F.↗

A real-time recursive filter for the attitude determination of the Spacelab instrument pointing subsystem

A real-time estimation filter which reduces sensitivity to system variations and reduces the amount of preflight computation is developed for the instrument pointing subsystem (IPS). The IPS is a three-axis stabilized platform developed to point various astronomical observation instruments aboard the shuttle. Currently, the IPS utilizes a linearized Kalman filter (LKF), with premission defined gains, to compensate for system drifts and accumulated attitude errors. Since the a priori gains are generated for an expected system, variations result in a suboptimal estimation process. This report compares the performance of three real-time estimation filters with the current LKF implementation. An extended Kalman filter and a second-order Kalman filter are developed to account for the system nonlinearities, while a linear Kalman filter implementation assumes that the nonlinearities are negligible. The performance of each of the four estimation filters are compared with respect to accuracy, stability, settling time, robustness, and computational requirements. It is shown, that for the current IPS pointing requirements, the linear Kalman filter provides improved robustness over the LKF with less computational requirements than the two real-time nonlinear estimation filters.

West, M. E.↗

The symbolic computation of series solutions to ordinary differential equations using trees (extended abstract)

Algorithms previously developed by the author give formulas which can be used for the efficient symbolic computation of series expansions to solutions of nonlinear systems of ordinary differential equations. As a by product of this analysis, formulas are derived which relate to trees to the coefficients of the series expansions, similar to the work of Leroux and Viennot, and Lamnabhi, Leroux and Viennot.

Grossman, Robert↗

A global nonlinear tracking problem

The tracking problem for a globally minimum phase nonlinear system is considered. It is assumed that the signal to be tracked is slowly varying and a priori bounds on its magnitude are known. It is shown that if the system has bounded derivatives and exponentially stable zero dynamics then the system admits an output feedback controller which solves the tracking problem.

Dayawansa, W. P.↗

Note on unit tangent vector computation for homotopy curve tracking on a hypercube

Probability-one homotopy methods are a class of methods for solving nonlinear systems of equations that are globally convergent from an arbitrary starting point. The essence of all such algorithms is the construction of an appropriate homotopy map and subsequent tracking of some smooth curve in the zero set of the homotopy map. Tracking a homotopy curve involves finding the unit tangent vector at different points along the zero curve, which amounts to calculating the kernel of the n x (n + 1) Jacobian matrix. While computing the tangent vector is just one part of the curve tracking algorithm, it can require a significant percentage of the total tracking time. This note presents computational results showing the performance of several different parallel orthogonal factorization/triangular system solving algorithms for the tangent vector computation on a hypercube.

Chakraborty, A.↗

Lava flows are fractals

Results are presented of a preliminary investigation of the fractal nature of the plan-view shapes of lava flows in Hawaii (based on field measurements and aerial photographs), as well as in Idaho and the Galapagos Islands (using aerial photographs only). The shapes of the lava flow margins are found to be fractals: lava flow shape is scale-invariant. This observation suggests that nonlinear forces are operating in them because nonlinear systems frequently produce fractals. A'a and pahoehoe flows can be distinguished by their fractal dimensions (D). The majority of the a'a flows measured have D between 1.05 and 1.09, whereas the pahoehoe flows generally have higher D (1.14-1.23). The analysis is extended to other planetary bodies by measuring flows from orbital images of Venus, Mars, and the moon. All are fractal and have D consistent with the range of terrestrial a'a and have D consistent with the range of terrestrial a'a and pahoehoe values.

Bruno, B. C.↗

Nonlinear control via approximate input-output linearization - The ball and beam example

A study is made of approximate input-output linearization of nonlinear systems which fail to have a well defined relative degree. For such systems, a method is provided for constructing approximate systems that are input-output linearizable. The analysis presented in this note is motivated through its application to a common undergraduate control laboratory experiment, the ball and beam system, where it is shown to be more effective for trajectory tracking than the standard Jacobian linearization.

Hauser, John↗

Resonant frequency calculations using a hybrid perturbation-Galerkin technique

A two-step hybrid perturbation Galerkin technique is applied to the problem of determining the resonant frequencies of one or several degrees of freedom nonlinear systems involving a parameter. In one step, the Lindstedt-Poincare method is used to determine perturbation solutions which are formally valid about one or more special values of the parameter (e.g., for large or small values of the parameter). In step two, a subset of the perturbation coordinate functions determined in step one is used in Galerkin type approximation. The technique is illustrated for several one degree of freedom systems, including the Duffing and van der Pol oscillators, as well as for the compound pendulum. For all of the examples considered, it is shown that the frequencies obtained by the hybrid technique using only a few terms from the perturbation solutions are significantly more accurate than the perturbation results on which they are based, and they compare very well with frequencies obtained by purely numerical methods.

Geer, James F.↗

Effects of structural offset, axial shortening, and gravitational torque on the slewing of a flexible spacecraft

The large-angle maneuvering of a Shuttle-beam-reflector spacecraft in the plane of a circular earth orbit is examined by considering the effects of the structural offset connection, the axial shortening, and the gravitational torque on the slewing motion. The offset effect is analyzed by changing the attachment point of the reflector to the beam. As the attachment point is moved away from the mass center of the reflector, the responses of the nonlinear system deviate from those of the linearized system. The axial geometric shortening effect induced by the deformation of the beam contributes to the system equations through second order terms in the modal amplitudes and rates. The gravitational torque effect is relatively small.

Li, Feiyue↗

Recognition of partially occluded threat objects using the annealed Hopefield network

Recognition of partially occluded objects has been an important issue to airport security because occlusion causes significant problems in identifying and locating objects during baggage inspection. The neural network approach is suitable for the problems in the sense that the inherent parallelism of neural networks pursues many hypotheses in parallel resulting in high computation rates. Moreover, they provide a greater degree of robustness or fault tolerance than conventional computers. The annealed Hopfield network which is derived from the mean field annealing (MFA) has been developed to find global solutions of a nonlinear system. In the study, it has been proven that the system temperature of MFA is equivalent to the gain of the sigmoid function of a Hopfield network. In our early work, we developed the hybrid Hopfield network (HHN) for fast and reliable matching. However, HHN doesn't guarantee global solutions and yields false matching under heavily occluded conditions because HHN is dependent on initial states by its nature. In this paper, we present the annealed Hopfield network (AHN) for occluded object matching problems. In AHN, the mean field theory is applied to the hybird Hopfield network in order to improve computational complexity of the annealed Hopfield network and provide reliable matching under heavily occluded conditions. AHN is slower than HHN. However, AHN provides near global solutions without initial restrictions and provides less false matching than HHN. In conclusion, a new algorithm based upon a neural network approach was developed to demonstrate the feasibility of the automated inspection of threat objects from x-ray images. The robustness of the algorithm is proved by identifying occluded target objects with large tolerance of their features.

Kim, Jung H.↗

A new flux conserving Newton's method scheme for the two-dimensional, steady Navier-Stokes equations

A new numerical method is developed for the solution of the two-dimensional, steady Navier-Stokes equations. The method that is presented differs in significant ways from the established numerical methods for solving the Navier-Stokes equations. The major differences are described. First, the focus of the present method is on satisfying flux conservation in an integral formulation, rather than on simulating conservation laws in their differential form. Second, the present approach provides a unified treatment of the dependent variables and their unknown derivatives. All are treated as unknowns together to be solved for through simulating local and global flux conservation. Third, fluxes are balanced at cell interfaces without the use of interpolation or flux limiters. Fourth, flux conservation is achieved through the use of discrete regions known as conservation elements and solution elements. These elements are not the same as the standard control volumes used in the finite volume method. Fifth, the discrete approximation obtained on each solution element is a functional solution of both the integral and differential form of the Navier-Stokes equations. Finally, the method that is presented is a highly localized approach in which the coupling to nearby cells is only in one direction for each spatial coordinate, and involves only the immediately adjacent cells. A general third-order formulation for the steady, compressible Navier-Stokes equations is presented, and then a Newton's method scheme is developed for the solution of incompressible, low Reynolds number channel flow. It is shown that the Jacobian matrix is nearly block diagonal if the nonlinear system of discrete equations is arranged approximately and a proper pivoting strategy is used. Numerical results are presented for Reynolds numbers of 100, 1000, and 2000. Finally, it is shown that the present scheme can resolve the developing channel flow boundary layer using as few as six to ten cells per channel width, depending on the Reynolds number.

Scott, James R.↗

Comparison between the performance of two classes of fuzzy controllers

This paper presents an application comparison between two classes of fuzzy controllers: the Clearness Transformation Fuzzy Controller (CTFC) and the CRI-based Fuzzy Controller. The comparison is performed by studying the application of the controllers to simulation examples of nonlinear systems. The CTFC is a new approach for the organization of fuzzy controllers based on a cognitive model of parameter driven control, the notion of fuzzy patterns to represent fuzzy knowledge and the Clearness Transformation Rule of Inference (CTRI) for approximate reasoning. The approach facilitates the implementation of the basic modules of the controller: the fuzzifier, defuzzifier, and the control protocol in a rule-based architecture. The CTRI scheme for approximate reasoning does not require the formation of fuzzy relation matrices yielding improved performance in comparison with the traditional organization of fuzzy controllers.

Janabi, T. H.↗

Error-based adaptive non-linear control and regions of feasibility

The nonlinear adaptive algorithm of Kanellakopoulos et al. (1991) was modified to produce an error-based algorithm. This permits global stabilizability for a large subset of pure-feedback nonlinear systems. The algorithm was demonstrated on the single-input stabilization problem, but extends easily to the multiple input tracking problems.

Teel, Andrew R.↗

Pinching solutions of slender cylindrical jets

Simplified equations for slender jets are derived for a circular jet of one fluid flowing into an ambient second fluid, the flow being confined in a circular tank. Inviscid flows are studied which include both surface tension effects and Kelvin-Helmholtz instability. For slender jets a coupled nonlinear system of equations is found for the jet shape and the axial velocity jump across it. The equations can break down after a finite time and similarity solutions are constructed, and studied analytically and numerically. The break-ups found pertain to the jet pinching after a finite time, without violation of the slender jet ansatz. The system is conservative and admissible singular solutions are those which conserve the total energy, mass, and momentum. Such solutions are constructed analytically and numerically, and in the case of vortex sheets with no surface tension certain solutions are given in closed form.

Papageorgiou, Demetrios T.↗

Delay analysis of combinations of pass transistors and classical logic gates

Combinations of pass transistors and logic gates driving nonlinear capacitive loads are analyzed for the presence of characteristics that will permit easier and more accurate digital logic simulation. Accurate time delay models are developed by studying the nature of the response of simplified circuit models to variations of input waveform rise and fall times and output loading. The nonlinear effects of the CMOS logic devices are minimized to permit easier interpretation of the influence of nonlinear capacitive loads. The performance of a CMOS inverter with a complex nonlinear load consisting of a pass transistor that separates a range of capacitances is compared to the same inverter circuit with a linear capacitive load to develop an understanding of the unique requirements of modeling a nonlinear system. Several methods of modeling the delay of CMOS circuits are reviewed, and a multi-parameter linear model is described. General guidelines for designing CMOS circuits with complex load circuits are developed, emphasizing that the circuit output rise delays and fall delays must be separately analyzed.

Gibson, Jonathan S.↗

The block adaptive multigrid method applied to the solution of the Euler equations

In the present study, a scheme capable of solving very fast and robust complex nonlinear systems of equations is presented. The Block Adaptive Multigrid (BAM) solution method offers multigrid acceleration and adaptive grid refinement based on the prediction of the solution error. The proposed solution method was used with an implicit upwind Euler solver for the solution of complex transonic flows around airfoils. Very fast results were obtained (18-fold acceleration of the solution) using one fourth of the volumes of a global grid with the same solution accuracy for two test cases.

Pantelelis, Nikos↗