Search NASA⌕ Search

SEARCH · Search NASA

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 523 records · Page 29

Frequency domain compensation of a DYNGEN turbofan engine model

Following Rosenbrock's ideas regarding the advantages of dominance in linear multivariable control systems, a new graphical technique is used for the design of compensators that achieve dominance. The technique is illustrated with an application to the problem of designing compensators for a linear turbofan-engine model. The resulting design is put into perspective by examining it in the light of two other multivariable frequency-domain methods. One, MacFarlane's method of characteristic loci, is used to realize a final design for stability and low interaction. The other is a direct technique based upon the algebraic expansion of the determinant of the return difference in terms of it's elements. Results from simulations carried out on the NASA DYNGEN software are included.

Schafer, R. M.↗

An algebraic structure of discrete-time biaffine systems

New results on the realization of finite-dimensional, discrete-time, internally biaffine systems are presented in this paper. The external behavior of such systems is described by multiaffine functions and the state space is constructed via Nerode equivalence relations. We prove that the state space is an affine space. An algorithm which amounts to choosing a frame for the affine space is presented. Our algorithm reduces in the linear and bilinear case to a generalization of algorithms existing in the literature. Explicit existence criteria for span-canonical realizations as well as an affine isomorphism theorem are given.

Tarn, T.-J.↗

Semilinear (topological) spaces and applications

Semivector spaces are defined and some of their algebraic aspects are developed including some structure theory. These spaces are then topologized to obtain semilinear topological spaces for which a hierarchy of local convexity axioms is identified. A number of fixed point and minmax theorems for spaces with various local convexity properties are established. The spaces of concern arise naturally as various hyperspaces of linear and semilinear (topological) spaces. It is indicated briefly how all this can be applied in socio-economic analysis and optimization.

Prakash, P.↗

Nonlinear adhesive behavior effects in a cracked orthotropic sheet stiffened by a semi-infinite orthotropic sheet

The stress-intensity factors are determined for a cracked orthotropic sheet adhesively bonded to an orthotropic stringer where the adhesive layer is modeled with a nonlinear stress-strain curve. By the use of Green's functions and the complex variable theory of orthotropic elasticity, a set of integral equations is obtained. The integral equations are replaced by an equivalent set of algebraic equations, which are solved to obtain the shear stress distribution in the adhesive layer, with which the crack-tip stress-intensity factors are found. When the adhesive was modeled with a nonlinear stress-strain curve, the peak shear stresses in the adhesive were considerably reduced in comparison to the solution for the linear elastic adhesive. This resulted in increases in the stress-intensity factors for the nonlinear adhesive solution compared to the linear adhesive solution. The nonlinear adhesive has no significant effect on the stress-intensity factor unless the near crack tip is beneath the stringer. It is assumed that the adhesive bond remains intact and it is predicted that onset of adhesive failure occurs at decreasing levels of applied stress as the crack propagates beneath the stringer.

Bigelow, C. A.↗

Nonlinear adhesive behavior effects in a cracked orthotropic sheet stiffened by a semi-infinite orthotropic sheet

The stress-intensity factors are determined for a cracked orthotropic sheet adhesively bonded to an orthotropic stringer where the adhesive layer is modeled with a nonlinear stress-strain curve. By the use of Green's functions and the complex variable theory of orthotropic elasticity, a set of integral equations is obtained. The integral equations are replaced by an equivalent set of algebraic equations, which are solved to obtain the shear stress distribution in the adhesive layer, with which the crack-tip stress-intensity factors are found. When the adhesive was modeled with a nonlinear stress-strain curve, the peak shear stresses in the adhesive were considerably reduced in comparison to the solution for the linear elastic adhesive. This resulted in increases in the stress-intensity factors for the nonlinear adhesive solution compared to the linear adhesive solution. The nonlinear adhesive has no significant effect on the stress-intensity factor unless the near crack tip is beneath the stringer. It is assumed that the adhesive bond remains intact and it is predicted that onset of adhesive failure occurs at decreasing levels of applied stress as the crack propagates beneath the stringer.

Bigelow, C. A.↗

Output feedback and generic stabilizability

Attention is given to questions of pole placement and stabilization for generic linear systems with prescribed state, input and output dimension for the case where the controller must be implemented by linear memoryless output feedback. A criterion is presented in terms of a special pole placement property for generic stabilizability, and this criterion is used to describe constraints on the dimension which are consistent with generic stabilizability. Also treated are the rationality and solvability by radicals of stabilizing or pole positioning gains, and the way in which decision algebra can theoretically handle existence questions for generic systems is described.

Byrnes, C. I.↗

Instability and sound emission from a flow over a curved surface

The growth and decay of a wavepacket convecting in a boundary layer over a concave-convex surface is studied numerically using direct computations of the Navier-Stokes equations. The resulting sound radiation is computed using the linearized Euler equations with the pressure from the Navier-Stokes solution as a time-dependent boundary condition. It is shown that on the concave portion the amplitude of the wavepacket increases and its bandwidth broadens while on the convex portion some of the components in the packet are stabilized. The pressure field decays exponentially away from the surface and then algebraically exhibits a decay characteristic of acoustic waves in two dimensions. The far-field acoustic pressure exhibits a peak at a frequency corresponding to the inflow instability frequency.

Maestrello, L.↗

Application of sound and temperature to control boundary-layer transition

The growth and decay of a wave packet convecting in a boundary layer over a concave-convex surface and its active control by localized surface heating are studied numerically using direct computations of the Navier-Stokes equations. The resulting sound radiations are computed using linearized Euler equations with the pressure from the Navier-Stokes solution as a time-dependent boundary condition. It is shown that on the concave portion the amplitude of the wave packet increases and its bandwidth broadens while on the convex portion some of the components in the packet are stabilized. The pressure field decays exponentially away from the surface and then algebraically, exhibiting a decay characteristic of acoustic waves in two dimensions. The far-field acoustic behavior exhibits a super-directivity type of behavior with a beaming downstream. Active control by surface heating is shown to reduce the growth of the wave packet but have little effect on acoustic far field behavior for the cases considered. Active control by sound emanating from the surface of an airfoil in the vicinity of the leading edge is experimentally investigated. The purpose is to control the separated region at high angles of attack. The results show that injection of sound at shedding frequency of the flow is effective in an increase of lift and reduction of drag.

Maestrello, Lucio↗

Inverse methods for assessing ship-of-opportunity networks and estimating circulation and winds from tropical expendable bathythermograph data

Inverse methods for estimating the surface ciculation of the equatorial Pacific by combining a linear reduced-gravity shallow-water model with the Tropical Ocean-Global Atmosphere ship-of-opportunity expendable bathythermograph (TOGA SOP XBT) observing program are examined. It is demonstrated that a simple linear model of the upper circulation of the equatorial Pacific can be successfully used as a weak constraint when smoothing the TOGA SOP XBT data. A circulation is sought as the weighted least squares fit to the dynamics and the data. The solution method is an expansion in representer functions, and the generalized inverse problem is thereby reduced from a functional problem to an algebraic problem for the coefficients of the representer. A specific inverse calculation using synthetic forcing and data is presented.

Bennett, Andrew F.↗

An extended structure-based model based on a stochastic eddy-axis evolution equation

We have proposed and implemented an extension of the structure-based model for weak deformations. It was shown that the extended model will correctly reduce to the form of standard k-e models for the case of equilibrium under weak mean strain. The realizability of the extended model is guaranteed by the method of its construction. The predictions of the proposed model were very good for rotating homogeneous shear flows and for irrotational axisymmetric contraction, but were seriously deficient in the case of plane strain and axisymmetric expansion. We have concluded that the problem behind these difficulties lies in the algebraic constitutive equation relating the Reynolds stresses to the structure parameters rather than in the slow model developed here. In its present form, this equation assumes that under irrotational strain the principal axes of the Reynolds stresses remain locked onto those of the eddy-axis tensor. This is correct in the RDT limit, but inappropriate under weaker mean strains, when the non-linear eddy-eddy interactions tend to misalign the two sets of principal axes and create some non-zero theta and gamma.

Kassinos, S. C.↗

Nonlinear, nonbinary cyclic group codes

New cyclic group codes of length 2(exp m) - 1 over (m - j)-bit symbols are introduced. These codes can be systematically encoded and decoded algebraically. The code rates are very close to Reed-Solomon (RS) codes and are much better than Bose-Chaudhuri-Hocquenghem (BCH) codes (a former alternative). The binary (m - j)-tuples are identified with a subgroup of the binary m-tuples which represents the field GF(2 exp m). Encoding is systematic and involves a two-stage procedure consisting of the usual linear feedback register (using the division or check polynomial) and a small table lookup. For low rates, a second shift-register encoding operation may be invoked. Decoding uses the RS error-correcting procedures for the m-tuple codes for m = 4, 5, and 6.

Solomon, G.↗

Numerical Solution of Incompressible Navier-Stokes Equations Using a Fractional-Step Approach

A fractional step method for the solution of steady and unsteady incompressible Navier-Stokes equations is outlined. The method is based on a finite volume formulation and uses the pressure in the cell center and the mass fluxes across the faces of each cell as dependent variables. Implicit treatment of convective and viscous terms in the momentum equations enables the numerical stability restrictions to be relaxed. The linearization error in the implicit solution of momentum equations is reduced by using three subiterations in order to achieve second order temporal accuracy for time-accurate calculations. In spatial discretizations of the momentum equations, a high-order (3rd and 5th) flux-difference splitting for the convective terms and a second-order central difference for the viscous terms are used. The resulting algebraic equations are solved with a line-relaxation scheme which allows the use of large time step. A four color ZEBRA scheme is employed after the line-relaxation procedure in the solution of the Poisson equation for pressure. This procedure is applied to a Couette flow problem using a distorted computational grid to show that the method minimizes grid effects. Additional benchmark cases include the unsteady laminar flow over a circular cylinder for Reynolds Numbers of 200, and a 3-D, steady, turbulent wingtip vortex wake propagation study. The solution algorithm does a very good job in resolving the vortex core when 5th-order upwind differencing and a modified production term in the Baldwin-Barth one-equation turbulence model are used with adequate grid resolution.

Kiris, Cetin↗

The effect on engine performance of change in jacket-water outlet temperature

Tests made on a Curtiss D-12 engine in the Altitude Laboratory at the Bureau of Standards show the following effects on engine performance of change in jacket-water outlet temperature: 1) Friction at all altitudes is a linear function of the jacket-water temperature, decreasing with increasing temperature. 2) The brake horsepower below an altitude of about 9,000 feet decreases, and at higher altitudes increases, with jacket-water temperature. 3) The brake specific fuel consumption tends to decrease, at all altitudes, with increasing jacket-water temperature. 4) The percentage change in brake power output is roughly equal to the algebraic sum of the percentage change in volumetric efficiency and mechanical efficiency.

Garlock, E A↗

A Method to Solve Interior and Exterior Camera Calibration Parameters for Image Resection

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.

Samtaney, Ravi↗

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 1 is the Analysis Description, and describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Towne, Charles E.↗

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 2: User's guide

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 2 is the User's Guide, and describes the program's general features, the input and output, the procedure for setting up initial conditions, the computer resource requirements, the diagnostic messages that may be generated, the job control language used to run the program, and several test cases.

Towne, Charles E.↗

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 3: Programmer's reference

A new computer code was developed to solve the 2-D or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating-direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 3 is the Programmer's Reference, and describes the program structure, the FORTRAN variables stored in common blocks, and the details of each subprogram.

Towne, Charles E.↗

A Nonlinear Modal Aeroelastic Solver for FUN3D

A nonlinear structural solver has been implemented internally within the NASA FUN3D computational fluid dynamics code, allowing for some new aeroelastic capabilities. Using a modal representation of the structure, a set of differential or differential-algebraic equations are derived for general thin structures with geometric nonlinearities. ODEPACK and LAPACK routines are linked with FUN3D, and the nonlinear equations are solved at each CFD time step. The existing predictor-corrector method is retained, whereby the structural solution is updated after mesh deformation. The nonlinear solver is validated using a test case for a flexible aeroshell at transonic, supersonic, and hypersonic flow conditions. Agreement with linear theory is seen for the static aeroelastic solutions at relatively low dynamic pressures, but structural nonlinearities limit deformation amplitudes at high dynamic pressures. No flutter was found at any of the tested trajectory points, though LCO may be possible in the transonic regime.

Goldman, Benjamin D.↗