Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear 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 109 records · Page 6

NASA-Ames three-dimensional potential flow analysis system (POTFAN) equation solver code (SOLN) version 1

A computer program known as SOLN was developed as an independent segment of the NASA-Ames three-dimensional potential flow analysis systems of linear algebraic equations. Methods used include: LU decomposition, Householder's method, a partitioning scheme, and a block successive relaxation method. Due to the independent modular nature of the program, it may be used by itself and not necessarily in conjunction with other segments of the POTFAN system.

Davis, J. E.↗

A nonlinear finite-element analysis of wings in steady incompressible flows with wake roll-up

The problem of lifting surfaces and complex aircraft configurations in steady incompressible flow is considered. For lifting surfaces the problem is formulated in terms of an integral equation relating the potential discontinuity on wing and wake to the normal derivative of the potential on the lifting surface. For complex configurations the problem is formulated in terms of an integral equation relating the potential to its normal derivative on the surface of the aircraft. The integral equation is approximated by a system of linear algebraic equations obtained by dividing the surfaces into small quadrilateral elements and by assuming the potential (or the potential discontinuity) and its normal derivative to be constant within each element. The wake geometry is obtained by iteration by satisfying the condition that the velocity be tangent to the surface of the wake and that the potential discontinuity be constant along the streamlines.

Suciu, E. O.↗

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.↗

Nonlinear steady incompressible lifting-surface analysis with wake roll-up

The problem of lifting surfaces in steady incompressible flow is considered in terms of an integral equation relating the potential discontinuity on wing and wake to the normal derivative of the potential (normal wash) on the lifting surface. The integral equation is approximated by a system of linear algebraic equations obtained by dividing the surfaces into small quadrilateral elements and by assuming the potential discontinuity and the normal wash to be constant within each element. The wake geometry is obtained through iteration by satisfying the condition that the velocity be tangent to the surface of the wake and that the potential discontinuity be constant along the streamlines. Numerical results are in good agreement with existing ones

Suciu, E. O.↗

Unsteady flow in a supersonic cascade with strong in-passage shocks

Linearized theory is used to study the unsteady flow in a supersonic cascade with in-passage shock waves. We use the Wiener-Hopf technique to obtain a closed-form analytical solution for the supersonic region. To obtain a solution for the rotational flow in the subsonic region we must solve an infinite set of linear algebraic equations. The analysis shows that it is possible to correlate quantitatively the oscillatory shock motion with the Kutta condition at the trailing edges of the blades. This feature allows us to account for the effect of shock motion on the stability of the cascade. Unlike the theory for a completely supersonic flow, the present study predicts the occurrence of supersonic bending flutter. It therefore provides a possible explanation for the bending flutter that has recently been detected in aircraft-engine compressors at higher blade loadings.

Goldstein, M. E.↗

Calculation of unsteady transonic flows using the integral equation method

The basic integral equations for a harmonically oscillating airfoil in a transonic flow with shock waves are derived; the reduced frequency is assumed to be small. The problems associated with shock wave motion are treated using a strained coordinate system. The integral equation is linear and consists of both line integrals and surface integrals over the flow field which are evaluated by quadrature. This leads to a set of linear algebraic equations that can be solved directly. The shock motion is obtained explicitly by enforcing the condition that the flow is continuous except at a shock wave. Results obtained for both lifting and nonlifting oscillatory flows agree satisfactorily with other accurate results.

Nixon, D.↗

Image-analysis library

MATHPAC image-analysis library is collection of general-purpose mathematical and statistical routines and special-purpose data-analysis and pattern-recognition routines for image analysis. MATHPAC library consists of Linear Algebra, Optimization, Statistical-Summary, Densities and Distribution, Regression, and Statistical-Test packages.

Source record↗

Ball bearing heat analysis program (BABHAP)

The Ball Bearing Heat Analysis Program (BABHAP) is an attempt to assemble a series of equations, some of which are non-linear algebraic systems, in a logical order, which when solved, provide a complex analysis of load distribution among the balls, ball velocities, heat generation resulting from friction, applied load, and ball spinning, minimum lubricant film thickness, and many additional characteristics of ball bearing systems. Although initial design requirements for BABHAP were dictated by the core limitations of the PDP 11/45 computer, (approximately 8K of real words with limited number of instructions) the program dimensions can easily be expanded for large core computers such as the UNIVAC 1108. The PDP version of BABHAP is also operational on the UNIVAC system with the exception that the PDP uses 029 punch and the UNIVAC uses 026. A conversion program was written to allow transfer between machines.

Source record↗

The application of sensitivity analysis to models of large scale physiological systems

A survey of the literature of sensitivity analysis as it applies to biological systems is reported as well as a brief development of sensitivity theory. A simple population model and a more complex thermoregulatory model illustrate the investigatory techniques and interpretation of parameter sensitivity analysis. The role of sensitivity analysis in validating and verifying models, and in identifying relative parameter influence in estimating errors in model behavior due to uncertainty in input data is presented. This analysis is valuable to the simulationist and the experimentalist in allocating resources for data collection. A method for reducing highly complex, nonlinear models to simple linear algebraic models that could be useful for making rapid, first order calculations of system behavior is presented.

Leonard, J. I.↗

An analysis of interface crack in adhesively bonded lap shear joints

This paper presents an analytical study of interface crack problems in adhesively bonded lap shear joints. A convenient and accurate method of analysis is developed on the basis of conservation laws in elasticity and of the fundamental relationships in fracture mechanics of interface crack problems. The problem is eventually reduced to a system of linear algebraic equations of mixed-mode stress intensity factor solutions. One of the salient features of the present approach is that the interface crack-tip stress intensity solutions can be determined accurately and directly by the information extracted from the far field. Solutions of the interface crack problems in lap shear joints have been obtained for the cases of various adhesive/adherend material systems and geometric configurations of the crack and the joint. Of particular interest is the evaluation of the mixed-mode stress intensity factors and the associated energy release rates. They are determined explicitly to illustrate the fundamental nature of the joint fracture problem.

Wang, S. S.↗

A mixed-mode crack analysis of rectilinear anisotropic solids using conservation laws of elasticity

A very simple and convenient method of analysis for studying two-dimensional mixed-mode crack problems in rectilinear anisotropic solids is presented. The analysis is formulated on the basis of conservation laws of anisotropic elasticity and of fundamental relationships in anisotropic fracture mechanics. The problem is reduced to a system of linear algebraic equations in mixed-mode stress intensity factors. One of the salient features of the present approach is that it can determine directly the mixed-mode stress intensity solutions from the conservation integrals evaluated along a path removed from the crack-tip region without the need of solving the corresponding complex near-field boundary value problem. Several examples with solutions available in the literature are solved to ensure the accuracy of the current analysis. This method is further demonstrated to be superior to other approaches in its numerical simplicity and computational efficiency. Solutions of more complicated and practical engineering problems dealing with the crack emanating from a circular hole in composites are presented also to illustrate the capacity of this method.

Wang, S. S.↗

Iterative optical vector-matrix processors (survey of selected achievable operations)

An iterative optical vector-matrix multiplier with a microprocessor-controlled feedback loop capable of performing a wealth of diverse operations was described. A survey and description of many of its operations demonstrates the versatility and flexibility of this class of optical processor and its use in diverse applications. General operations described include: linear difference and differential equations, linear algebraic equations, matrix equations, matrix inversion, nonlinear matrix equations, deconvolution and eigenvalue and eigenvector computations. Engineering applications being addressed for these different operations and for the IOP are: adaptive phased-array radar, time-dependent system modeling, deconvolution and optimal control.

Casasent, D.↗

Analysis of cracks emanating from a circular hole in unidirectional fiber reinforced composites, part 2

An analytical method is developed for cracks emanating from a circular hole in an off-axis unidirectional fiber-reinforced composite. The method which is formulated by using conservation laws of elasticity and fundamental relationships in anisotropic fracture mechanics, provides a convenient and accurate means to examine the complicated crack behavior, when used in conjunction with a suitable numerical scheme such as the finite element method. The formulation is eventually reduced to a system of linear algebraic equations of mixed-mode stress intensity factors. Fracture parameters, describing crack-tip deformation and fracture in the composite, are obtained explicitly. Effects of material anisotropy and crack/hole geometry are examined also. Of particular interest are the energy release rates associated with crack extension; their values are evaluated for various cases. Results show that mixed-mode stress intensity factors and energy release rates associated with the cracks emanating from a hole change very appreciably with fiber orientation in the composite. K sub 1 and G increase monotonically with increasing theta; but K sub 2 reaches its maximum at theta = 45 deg, and then decreases gradually as theta increases further.

Wang, S. S.↗

Analysis of interface cracks in adhesively bonded lap shear joints, part 4

Conservation laws of elasticity for nonhomogeneous materials were developed and were used to study the crack behavior in adhesively bonded lap shear joints. By using these laws and the fundamental relationships in fracture mechanics of interface cracks, the problem is reduced to a pair of linear algebraic equations, and stress intensity solutions can be determined directly by information extracted from the far field. The numerical results obtained show that: (1) in the lap-shear joint with a given adherend, the opening-mode stress intensity factor, (K sub 1) is always larger than that of the shearing-mode (K sub 2); (2) (K sub 1) is not sensitive to adherent thickness abut (K sub 2) increases rapidly with increasing thickness; and (3) (K sub 1) and (K sub 2) increase simultaneously as the interfacial crack length increases.

Wang, S. S.↗

Periodic solutions of the Riccati equation

It is shown that the Riccati equations of control theory are generically Morse-Smale. The existence of periodic solutions is explicitly demonstrated. The principal contribution here is thought to be the presentation of a mathematical structure in which to study the Riccati equation. Mathematical structures already in the literature are systematically exploited. It is also shown that most problems of the Riccati equation can be reduced to problems in linear algebra. A framework is developed for studying the Riccati equation, and the proof that the Morse-Smale conditions are satisfied is outlined.

Hermann, R.↗

Subsonic 3-D surface panel method for rapid analysis of multiple geometry perturbations

For any baseline aircraft configuration, a matrix of partial derivatives of surface velocity potential with respect to geometry coordinates is calculated. Linear extrapolation is then used to analyze the subcritical potential flow corresponding to a series of arbitrary small geometry perturbations. Each perturbation analysis is more than an order of magnitude more efficient than a conventional panel method solution because no influence coefficients are calculated and no large system of linear algebraic equations is solved. Wing and wing-fuselage examples are presented to demonstrate that the predicted pressure distributions are nearly exact for large changes to wing camber, thickness, and leading edge radius.

Bristow, D. R.↗

Interface cracks in adhesively bounded lap-shear joints

A study on the elastic behavior of interface cracks in adhesively bonded lap-shear joints is presented. The problem is investigated by using a recently developed method of analysis based on conservation laws in elasticity for nonhomogeneous solids and fundamental relationships in fracture mechanics of dissimilar materials. The formulation leads to a pair of linear algebraic equations in mixed-mode stress intensity factors. Singular crack-tip stress intensity solutions are determined directly by information extracted from the far field. Stress intensity factors and associated energy release rates are obtained for various cases of interest. Fundamental nature of the interfacial flaw behavior in lap-shear adhesive joints is examined in detail.

Wang, S. S.↗

A tire contact solution technique

An efficient method for calculating the contact boundary and interfacial pressure distribution was developed. This solution technique utilizes the discrete Fourier transform to establish an influence coefficient matrix for the portion of the pressurized tire surface that may be in the contact region. This matrix is used in a linear algebra algorithm to determine the contact boundary and the array of forces within the boundary that are necessary to hold the tire in equilibrium against a specified contact surface. The algorithm also determines the normal and tangential displacements of those points on the tire surface that are included in the influence coefficient matrix. Displacements within and outside the contact region are calculated. The solution technique is implemented with a finite-element tire model that is based on orthotropic, nonlinear shell of revolution elements which can respond to nonaxisymmetric loads. A sample contact solution is presented.

Tielking, J. T.↗