SEARCH · Engineering Papers
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.
Stochastic differential equations specifying dynamical structure of filters generating posterior probability distribution when inputs are time functions
Method of inverting nonlinear equations of rotational dynamics of rigid body used to design feedback control of orientation of body. Applicable to both direction-cosine and quaternion formulations suitable for large-angle maneuvers. Exploiting some apparently little-known properties of direction cosine and quaternion formulations, method leads to equations for model-follower control system that exhibits exactly linear attitude-error dynamics. Quarternion system more robust in responding to large roll-angle commands.
In recent years, the SUite of Nonlinear and DIfferential/ALgebraic equation Solvers (SUNDIALS) has been redesigned to better enable the use of application-specific and third-party algebraic solvers and data structures. Throughout this work, we have adhered to specific guiding principles that minimized the impact to current users while providing maximum flexibility for later evolution of solvers and data structures. The redesign was done through the addition of new linear and nonlinear solvers classes, enhancements to the vector class, and the creation of modern Fortran interfaces. The vast majority of this work has been performed “behind-the-scenes,” with minimal changes to the user interface and no reduction in solver capabilities or performance. These changes allow SUNDIALS users to more easily utilize external solver libraries and create highly customized solvers, enabling greater flexibility on extreme-scale, heterogeneous computational architectures.
The nonlinear development of the most unstable Gortler vortex mode in boundary layer flows over curved walls is investigated. The most unstable Gortler mode is confined to a viscous wall layer of thickness O(G -1/5) and has spanwise wavelength O(G 11/5); it is, of course, most relevant to flow situations where the Gortler number G is much greater than 1. The nonlinear equations covering the evolution of this mode over an O(G -3/5) streamwise lengthscale are derived and are found to be of a fully nonparallel nature. The solution of these equations is achieved by making use of the numerical scheme used by Hall (1988) for the numerical solution of the nonlinear Gortler equations valid for O(1) Gortler numbers. Thus, the spanwise dependence of the flow is described by a Fourier expansion, whereas the streamwise and normal variations of the flow are dealt with by employing a suitable finite difference discretization of the governing equations. Our calculations demonstrate that, given a suitable initial disturbance, after a brief interval of decay, the energy in all the higher harmonics grows until a singularity is encountered at some downstream position. The structure of the flowfield as this singularity is approached suggests that the singularity is responsible for the vortices, which are initially confined to the thin viscous wall layer, moving away from the wall and into the core of the boundary layer.
The nonlinear development of the most unstable Goertler vortex mode in boundary layer flows over curved walls is investigated. The most unstable Goertler mode is confined to a viscous wall layer of thickness O(G -1/5) and has spanwise wavelength O(G 11/5); it is, of course, most relevant to flow situations where the Goertler number G is much greater than 1. The nonlinear equations covering the evolution of this mode over an O(G -3/5) streamwise lengthscale are derived and are found to be of a fully nonparallel nature. The solution of these equations is achieved by making use of the numerical scheme used by Hall (1988) for the numerical solution of the nonlinear Goertler equations valid for O(1) Goertler numbers. Thus, the spanwise dependence of the flow is described by a Fourier expansion, whereas the streamwise and normal variations of the flow are dealt with by employing a suitable finite difference discretization of the governing equations. Our calculations demonstrate that, given a suitable initial disturbance, after a brief interval of decay, the energy in all the higher harmonics grows until a singularity is encountered at some downstream position. The structure of the flowfield as this singularity is approached suggests that the singularity is responsible for the vortices, which are initially confined to the thin viscous wall layer, moving away from the wall and into the core of the boundary layer.
A general set of finite-element equations is developed for treating the nonlinear dynamic response of membranes consisting of complex nonlinear materials. The approach proposed by Key (1974) for a most complete treatment of the finite-element method to the nonlinear dynamic analysis of solids is used in developing the nonlinear equations of motion for membranes with cubic quadrilateral finite elements. The computations are kept general with respect to material modeling by computing the deformation gradient and requiring a constitutive subroutine to supply the Cauchy stress. An example calculation of the vibration of a square membrane demonstrates the accuracy of the formulation for small deformation theory and shows the nonlinear effect of large deformations.
A variable-skew oblique wing offers a substantial aerodynamic performance advantage for aircraft missions that require both high efficiency in subsonic flight and supersonic dash or cruise. The most obvious characteristic of the oblique-wing concept is the asymmetry associated with wing-skew angle which results in significant aerodynamic and inertial cross-coupling between the aircraft longitudinal and lateral-directional axes. A technique for synthesizing a decoupling controller while providing the desired stability augmentation. The proposed synthesis procedure uses the cncept of explicit model following. Linear quadratic optimization techniques are used to design the linear feedback system. The effectiveness of the control laws developed in achieving the desired decoupling is illustrated for a given flight condition by application to linearized equations of motion, and also to the nonlinear equations of six degrees of freedom of motion with nonlinear aerodynamic data.
A variable-skew oblique wing offers a substantial aerodynamic performance advantage for aircraft missions that require both high efficiency in subsonic flight and supersonic dash or cruise. The most obvious characteristic of the oblique-wing concept is the asymmetry associated with wing-skew angle which results in significant aerodynamic and inertial cross-coupling between the aircraft longitudinal and lateral-directional axes. This paper presents a technique for synthesizing a decoupling controller while providing the desired stability augmentation. The proposed synthesis procedure uses the concept of explicit model following. Linear quadratic optimization techniques are used to design the linear feedback system. The effectiveness of the control laws developed in achieving the desired decoupling is illustrated for a given flight condition by application to linearized equations of motion, and also to the nonlinear equations of six degrees of freedom of motion with nonlinear aerodynamic data.
Parameter-plane method analysis of nonlinear equations of asymmetric oscillations in adaptive feedback control systems
An iterative method is presented to solve the internal and external camera calibration parameters, given model target points and their images from one or more camera locations. The direct linear transform formulation was used to obtain a guess for the iterative method, and herein lies one of the strengths of the present method. In all test cases, the method converged to the correct solution. In general, an overdetermined system of nonlinear equations is solved in the least-squares sense. The iterative method presented is based on Newton-Raphson for solving systems of nonlinear algebraic equations. The Jacobian is analytically derived and the pseudo-inverse of the Jacobian is obtained by singular value decomposition.
A differential correction technique for solving nonlinear minimax problems is presented. The basis of the technique is a linear programing algorithm which solves the linear minimax problem. By linearizing the original nonlinear equations about a nominal solution, both nonlinear approximation and estimation problems using the minimax norm may be solved iteratively. Some consideration is also given to improving convergence and to the treatment of problems with more than one measured quantity. A sample problem is treated with this technique and with the least-squares differential correction method to illustrate the properties of the minimax solution. The results indicate that for the sample approximation problem, the minimax technique provides better estimates than the least-squares method if a sufficient amount of data is used. For the sample estimation problem, the minimax estimates are better if the mathematical model is incomplete.
The paper presents a review of the historical development in unsteady transonic aerodynamics, along with the foundations and accomplishments of several approaches to solve the equations of unsteady transonic flow. The discussion covers the linearized unsteady flow theory, numerical solution of the exact equations for an inviscid compressible gas, nonlinear small disturbance theory of transonic flow and linearization of the unsteady component about the nonlinear solution for the steady state, local linearization solution for unsteady transonic flow, unsteady transonic flow theory for slender wings and bodies, and three-dimensional unsteady transonic flows. The relation between the calculated results and experiment is examined. It is shown that the newly emerging numerical methods are capable of solving the nonlinear equations for two-dimensional flow and can be extended to three-dimensional flows.
The Navier-Stokes equation is transformed into a closed nonlinear equation for the convected curl operator. It is shown how Lie algebraic methods can be employed to solve the obtained double-commutator equation. The proposed operator equation provides a way to the further development of the Navier-Stokes fluid flow theory.
For the determination of aerodynamic characteristics such as lift, drag, and moment coefficients, it is crucial to compute the properties of steady flow past an airfoil. This investigation provides a set of flow dependent grid systems and initial flowfield guesses which substantially improve convergence rates when applied to the Euler equations for flows past an airfoil. The basic equations are examined, taking into account nonlinear equations which are difficult to solve. A good approximation to these equations under certain conditions can be obtained by introducing the so-called 'tangent gas approximation' considered by Woods (1961), in which the isentropic relation between rho and p is replaced by a tangent to a curve. Attention is given to the solution procedure, the analysis (direct) problem, and a comparison of the tangent gas solution with the converged Euler solution.
Two contrasting representations of land surface variability used in general circulation models (GCMs) are compared through an analysis of their corresponding surface energy balance equations. In one representation (the 'mixture' approach), different vegetation types are assumed to be homogeneously mixed over a grid square, so that the GCM atmosphere sees near-surface conditions pertaining to the mixture only. In the second representation (the 'mosaic' approach), different vegetation types are viewed as separate 'tiles' of a grid-square 'mosaic', and each tile interacts with the atmosphere independently. The mosaic approach is computationally simpler and in many ways more flexible than the mixture approach. Analytical solutions to the linearized energy balance equations and numerical solutions to the nonlinear equations both demonstrate that the mixture strategy, when applied to two coexisting vegetation types that differ only in canopy transpiration resistance, promotes both total turbulent flux and latent heat flux relative to the mosaic strategy. The effective differences between the strategies, however, are small over a wide range of conditions. In particular, the strategies are effectively equivalent when the transpiration resistances of the different vegetation types are of the same order of magnitude.
Nonlinear equations of motion of symmetrically laminated anisotropic plates are derived considering von Karman strains and transverse shear effects. Using a single-mode Galerkin procedure, the nonlinear modal equation is obtained. The direct equivalent linearization method is employed for solution of this equation. The effects of transverse shear on large deflection vibration of laminated plates with mixed boundary support conditions under random excitation are studied. Mean-square deflections and mean-square inplane strains are obtained for symmetric graphite-epoxy laminates. Using equilibrium equations and the continuity requirements, the mean-square transverse shear stresses are calculated. The results obtained will be useful in the sonic fatigue design of composite aircraft panels. The analysis is presented in detail for plates with two opposing edges simply supported and the remaining edges clamped.
Some of the basic finite-difference schemes that can be used to solve the nonlinear equations that describe unsteady inviscid and viscous transonic flow are reviewed. Numerical schemes for solving the unsteady Euler and Navier-Stokes, boundary-layer, and nonlinear potential equations are described. Emphasis is given to the elementary ideas used in constructing various numerical procedures, not specific details of any one procedure.
A review is presented of some recently developed numerical methods for the solution of nonlinear equations of mixed type. The methods considered use finite difference approximations to the differential equation. Central difference formulas are employed in the subsonic zone and upwind difference formulas are used in the supersonic zone. The relaxation method for the small disturbance equation is discussed and a description is given of difference schemes for the potential flow equation in quasi-linear form. Attention is also given to difference schemes for the potential flow equation in conservation form, the analysis of relaxation schemes by the time dependent analogy, the accelerated iterative method, and three-dimensional calculations.