Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differential algebra”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

A brief survey of constrained mechanics and variational problems in terms of differential forms

There has been considerable interest recently in constrained mechanics and variational problems. This is in part due to applied interests (such as 'non-holonomic mechanics in robotics') and in other part due to the fact that several schools of 'pure' mathematics have found that this classical subject is of importance for what they are trying to do. I have made various attempts at developing these subjects since my Lincoln lab days of the late 1950's. In this Chapter, I will sketch a Unified point of view, using Cartan's approach with differential forms. This has the advantage from the C-O-R viewpoint being developed in this Volume that the extension from 'smooth' to 'generalized' data is very systematic and algebraic. (I will only deal with the 'smooth' point of view in this Chapter; I will develop the 'generalized function' material at a later point.) The material presented briefly here about Variational Calculus and Constrained Mechanics can be found in more detail in my books, 'Differential Geometry and the Calculus of Variations', 'Lie Algebras and Quantum Mechanics', and 'Geometry, Physics and Systems'.

Hermann, Robert↗
NASA NTRS1993

Large-scale computation of incompressible viscous flow by least-squares finite element method

The least-squares finite element method (LSFEM) based on the velocity-pressure-vorticity formulation is applied to large-scale/three-dimensional steady incompressible Navier-Stokes problems. This method can accommodate equal-order interpolations and results in symmetric, positive definite algebraic system which can be solved effectively by simple iterative methods. The first-order velocity-Bernoulli function-vorticity formulation for incompressible viscous flows is also tested. For three-dimensional cases, an additional compatibility equation, i.e., the divergence of the vorticity vector should be zero, is included to make the first-order system elliptic. The simple substitution of the Newton's method is employed to linearize the partial differential equations, the LSFEM is used to obtain discretized equations, and the system of algebraic equations is solved using the Jacobi preconditioned conjugate gradient method which avoids formation of either element or global matrices (matrix-free) to achieve high efficiency. To show the validity of this scheme for large-scale computation, we give numerical results for 2D driven cavity problem at Re = 10000 with 408 x 400 bilinear elements. The flow in a 3D cavity is calculated at Re = 100, 400, and 1,000 with 50 x 50 x 50 trilinear elements. The Taylor-Goertler-like vortices are observed for Re = 1,000.

Jiang, Bo-Nan↗
NASA NTRS1977

A finite element algorithm for sound propagation in axisymmetric ducts containing compressible mean flow

The described investigation is concerned with he development of a finite element scheme which can be used in a study of the acoustics of aircraft-engine ducts. In the absence of suitable variational principles for acoustic fluctuations within an aircraft fan engine, an acoustic analysis must proceed directly from the differential equations which describe compressible flow. The derived equations cannot be solved algebraically. The numerical technique used for solving them makes use of a linear rectangular element of a type considered by Zienkiewicz (1971). Attention is given to aspects of element derivation, the global matrix assembly, the solution of the matrix equation, questions of acoustic attenuation, and illustrations of the potential of the current model in duct optimization.

Abrahamson, A. L.↗
NASA NTRS1981

A Galerkin type finite element method for rotary-wing aeroelasticity in hover and forward flight

A Galerkin finite element method for the spatial discretization of the nonlinear, nonselfadjoint, partial differential equations governing rotary-wing aeroelasticity is presented. This method reduces algebraic manipulative labor significantly when compared to the global Galerkin method based on assumed modes. Furthermore, the Galerkin finite element method is ideally suited to treat rotor blades with discontinuous mass and stiffness distribution and structurally redundant configurations as they appear in bearingless rotors. Implementation of the method is illustrated for the coupled flap-lag aeroelastic problem of hingeless rotor blades in hover and forward flight. Numerical results for stability and response illustrate the numerical properties and convergence behavior of the method. It is concluded that the Galerkin finite element method is a practical tool for solving rotary-wing aeroelastic stability and response problems.

Straub, F. K.↗
NASA NTRS1988

Investigating a hybrid perturbation-Galerkin technique using computer algebra

A two-step hybrid perturbation-Galerkin method is presented for the solution of a variety of differential equations type problems which involve a scalar parameter. The resulting (approximate) solution has the form of a sum where each term consists of the product of two functions. The first function is a function of the independent field variable(s) x, and the second is a function of the parameter lambda. In step one the functions of x are determined by forming a perturbation expansion in lambda. In step two the functions of lambda are determined through the use of the classical Bubnov-Gelerkin method. The resulting hybrid method has the potential of overcoming some of the drawbacks of the perturbation and Bubnov-Galerkin methods applied separately, while combining some of the good features of each. In particular, the results can be useful well beyond the radius of convergence associated with the perturbation expansion. The hybrid method is applied with the aid of computer algebra to a simple two-point boundary value problem where the radius of convergence is finite and to a quantum eigenvalue problem where the radius of convergence is zero. For both problems the hybrid method apparently converges for an infinite range of the parameter lambda. The results obtained from the hybrid method are compared with approximate solutions obtained by other methods, and the applicability of the hybrid method to broader problem areas is discussed.

Andersen, Carl M.↗
NASA NTRS1986

On substructuring algorithms and solution techniques for the numerical approximation of partial differential equations

Substructuring methods are in common use in mechanics problems where typically the associated linear systems of algebraic equations are positive definite. Here these methods are extended to problems which lead to nonpositive definite, nonsymmetric matrices. The extension is based on an algorithm which carries out the block Gauss elimination procedure without the need for interchanges even when a pivot matrix is singular. Examples are provided wherein the method is used in connection with finite element solutions of the stationary Stokes equations and the Helmholtz equation, and dual methods for second-order elliptic equations.

Gunzburger, M. D.↗
NASA NTRS1980

Research on the application of a decoupling algorithm for structure analysis

The mathematical theory for decoupling mth-order matrix differential equations is presented. It is shown that the decoupling precedure can be developed from the algebraic theory of matrix polynomials. The role of eigenprojectors and latent projectors in the decoupling process is discussed and the mathematical relationships between eigenvalues, eigenvectors, latent roots, and latent vectors are developed. It is shown that the eigenvectors of the companion form of a matrix contains the latent vectors as a subset. The spectral decomposition of a matrix and the application to differential equations is given.

Denman, E. D.↗
NASA NTRS1989

Prediction of separated transonic wing flows with a non-equilibrium algebraic model

A nonequilibrium algebraic turbulence model, which is based on the turbulence closure scheme of Johnson and King (1985), is proposed to predict separated transonic wing flows. The influence of history effects are modeled by solving a partial differential equation for the maximum total Reynolds shear stress, which is then used to scale the eddy viscosity of an algebraic model. The turbulence model is implemented in a three-dimensional, Reynolds-averaged Navier-Stokes code. Comparisons with experimental data are presented which show clearly that the nonequilibrium type of turbulence model is essential for accurate prediction of transonic separated flows.

Abid, Ridha↗
NASA NTRS1983

Symmetric linear systems - An application of algebraic systems theory

Dynamical systems which contain several identical subsystems occur in a variety of applications ranging from command and control systems and discretization of partial differential equations, to the stability augmentation of pairs of helicopters lifting a large mass. Linear models for such systems display certain obvious symmetries. In this paper, we discuss how these symmetries can be incorporated into a mathematical model that utilizes the modern theory of algebraic systems. Such systems are inherently related to the representation theory of algebras over fields. We will show that any control scheme which respects the dynamical structure either implicitly or explicitly uses the underlying algebra.

Hazewinkel, M.↗
NASA NTRS1987

Quadrupole source in prediction of the noise of rotating blades - A new source description

The aim of this paper is to perform a theoretical study of the quadrupole term of the Ffowcs Williams-Hawkings (FW-H) equation to obtain practical results for applications to rotating blades. The quadrupole term of the FW-H equation is algebraically manipulated into volume, surface and line sources using generalized function theory and differential geometry. The volume source is of the type in Lighthill's jet noise theory. The surface sources are on the blade and shock surfaces and the line source is at the trailing edge. It is shown that contribution of volume sources in the boundary layer and wakes can be written in the form of surface integrals. It is argued that the surface and line sources and the part of the volume sources in the boundary layer, wakes and vortices near the blades should be sufficient in calculation of the noise of high speed rotating blades. The integrals correspoding to the various sources appearing in the formula for calculation of the acoustic pressure are briefly derived.

Farassat, F.↗
NASA NTRS1999

A Textbook for a First Course in Computational Fluid Dynamics

This paper describes and discusses the textbook, Fundamentals of Computational Fluid Dynamics by Lomax, Pulliam, and Zingg, which is intended for a graduate level first course in computational fluid dynamics. This textbook emphasizes fundamental concepts in developing, analyzing, and understanding numerical methods for the partial differential equations governing the physics of fluid flow. Its underlying philosophy is that the theory of linear algebra and the attendant eigenanalysis of linear systems provides a mathematical framework to describe and unify most numerical methods in common use in the field of fluid dynamics. Two linear model equations, the linear convection and diffusion equations, are used to illustrate concepts throughout. Emphasis is on the semi-discrete approach, in which the governing partial differential equations (PDE's) are reduced to systems of ordinary differential equations (ODE's) through a discretization of the spatial derivatives. The ordinary differential equations are then reduced to ordinary difference equations (O(Delta)E's) using a time-marching method. This methodology, using the progression from PDE through ODE's to O(Delta)E's, together with the use of the eigensystems of tridiagonal matrices and the theory of O(Delta)E's, gives the book its distinctiveness and provides a sound basis for a deep understanding of fundamental concepts in computational fluid dynamics.

Zingg, D. W.↗
NASA NTRS1976

Finite difference procedure for boundary layers including effects of longitudinal and transverse curvatures

A second order viscous layer solution procedure has been developed that does in a consistent way include curvature effects and the corresponding normal pressure gradients. In the present system, the normal momentum equation is retained. The parabolic system of nonlinear partial differential equations is converted by linear finite differencing procedures to a system of linear algebraic equations and solved in primitive coordinates. The solutions have been shown to give smooth stable distributions for all the variables, most particularly the normal velocity which plays an important role in the interaction procedure. An algorithm for matching the viscous layer solution with a rotational characteristics outer solution has been developed.

Tassa, Y.↗
NASA NTRS1986

The algebraic criteria for the stability of control systems

This paper critically examines the standard algebraic criteria for the stability of linear control systems and their proofs, reveals important previously unnoticed connections, and presents new representations. Algebraic stability criteria have also acquired significance for stability studies of non-linear differential equation systems by the Krylov-Bogoljubov-Magnus Method, and allow realization conditions to be determined for classes of broken rational functions as frequency characteristics of electrical network.

Cremer, H.↗
NASA NTRS1971

Similarity analysis of differential equations by Lie group.

Methods for transforming partial differential equations into forms more suitable for analysis and solution are investigated. The idea of Lie's infinitesimal contact transformation group is introduced to develop a systematic method which involves mostly algebraic manipulations. A thorough presentation of the application of this general method to the problem of similarity analysis in a broader sense - namely, the similarity between partial and ordinary differential equations, boundary value and initial value problems, and nonlinear and linear equations - is given with new and very general methods evolved for deriving the possible groups of transformations.

Na, T. Y.↗
NASA NTRS1972

Application of a one-strip integral method to the unsteady supersonic aerodynamics of an inclined flat surface

The method of integral relations is applied in a one-strip approximation to the perturbation equations governing small motions of an inclined, sharp-edged, flat surface about the mean supersonic steady flow. Algebraic expressions for low reduced-frequency aerodynamics are obtained and a set of ordinary differential equations are obtained for general oscillatory motion. Results are presented for low reduced-frequency aerodynamics and for the variation of the unsteady forces with frequency. The method gives accurate results for the aerodynamic forces at low reduced frequency which are in good agreement with available experimental data. However, for cases in which the aerodynamic forces vary rapidly with frequency, the results are qualitatively correct, but of limited accuracy. Calculations indicate that for a range of inclination angles near shock detachment such that the flow in the shock layer is low supersonic, the aerodynamic forces vary rapidly both with inclination angle and with reduced frequency.

Bennett, R. M.↗