Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear equations”

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 199 records · Page 11

Leapfrog variants of iterative methods for linear algebra equations

Two iterative methods are considered, Richardson's method and a general second order method. For both methods, a variant of the method is derived for which only even numbered iterates are computed. The variant is called a leapfrog method. Comparisons between the conventional form of the methods and the leapfrog form are made under the assumption that the number of unknowns is large. In the case of Richardson's method, it is possible to express the final iterate in terms of only the initial approximation, a variant of the iteration called the grand-leap method. In the case of the grand-leap variant, a set of parameters is required. An algorithm is presented to compute these parameters that is related to algorithms to compute the weights and abscissas for Gaussian quadrature. General algorithms to implement the leapfrog and grand-leap methods are presented. Algorithms for the important special case of the Chebyshev method are also given.

Saylor, Paul E.↗

Bounds on the Fourier coefficients for the periodic solutions of non-linear oscillator equations

The differential equations describing nonlinear oscillations (as seen in mechanical vibrations, electronic oscillators, chemical and biochemical reactions, acoustic systems, stellar pulsations, etc.) are investigated analytically. The boundedness of the Fourier coefficients for periodic solutions is demonstrated for two special cases, and the extrapolation of the results to higher-dimensionsal systems is briefly considered.

Mickens, R. E.↗

On the Coupling Between a Supersonic Turbulent Boundary Layer and a Flexible Structure

A mathematical model and a computer code have been developed to fully couple the vibration of an aircraft fuselage panel to the surrounding flow field, turbulent boundary layer and acoustic fluid. The turbulent boundary layer model is derived using a triple decomposition of the flow variables and applying a conditional averaging to the resulting equations. Linearized panel and acoustic equations are used. Results from this model are in good agreement with existing experimental and numerical data. It is shown that in the supersonic regime, full coupling of the flexible panel leads to lower response and radiation from the panel. This is believed to be due to an increase in acoustic damping on the panel in this regime. Increasing the Mach number increases the acoustic damping, which is in agreement with earlier work.

Frendi, Abdelkader↗

An approximately factored incremental strategy for calculating consistent discrete aerodynamic sensitivity derivatives

An incremental strategy is presented for iteratively solving very large systems of linear equations, which are associated with aerodynamic sensitivity derivatives for advanced CFD codes. It is shown that the left-hand side matrix operator and the well-known factorization algorithm used to solve the nonlinear flow equations can also be used to efficiently solve the linear sensitivity equations. Two airfoil problems are considered as an example: subsonic low Reynolds number laminar flow and transonic high Reynolds number turbulent flow.

Korivi, V. M.↗

Preprocessing Inconsistent Linear System for a Meaningful Least Squares Solution

Mathematical models of many physical/statistical problems are systems of linear equations. Due to measurement and possible human errors/mistakes in modeling/data, as well as due to certain assumptions to reduce complexity, inconsistency (contradiction) is injected into the model, viz. the linear system. While any inconsistent system irrespective of the degree of inconsistency has always a least-squares solution, one needs to check whether an equation is too much inconsistent or, equivalently too much contradictory. Such an equation will affect/distort the least-squares solution to such an extent that renders it unacceptable/unfit to be used in a real-world application. We propose an algorithm which (i) prunes numerically redundant linear equations from the system as these do not add any new information to the model, (ii) detects contradictory linear equations along with their degree of contradiction (inconsistency index), (iii) removes those equations presumed to be too contradictory, and then (iv) obtain the minimum norm least-squares solution of the acceptably inconsistent reduced linear system. The algorithm presented in Matlab reduces the computational and storage complexities and also improves the accuracy of the solution. It also provides the necessary warning about the existence of too much contradiction in the model. In addition, we suggest a thorough relook into the mathematical modeling to determine the reason why unacceptable contradiction has occurred thus prompting us to make necessary corrections/modifications to the models - both mathematical and, if necessary, physical.

Sen, Syamal K.↗

Preprocessing in Matlab Inconsistent Linear System for a Meaningful Least Squares Solution

Mathematical models of many physical/statistical problems are systems of linear equations~ Due to measurement and possible human errors/mistakes in modeling/data, as well as due to certain assumptions to reduce complexity, inconsistency (contradiction) is injected into the model, viz. the linear system. While any inconsistent system irrespective of the degree of inconsistency has always a least-squares solution, one needs to check whether an equation is too much inconsistent or, equivalently too much contradictory. Such an equation will affect/distort the least-squares solution to such an extent that renders it unacceptable/unfit to be used in a real-world application. We propose an algorithm which (i) prunes numerically redundant linear equations from the system as these do not add any new information to the model, (ii) detects contradictory linear equations along with their degree of contradiction (inconsistency index), (iii) removes those equations presumed to be too contradictory, and then (iv) obtain the . minimum norm least-squares solution of the acceptably inconsistent reduced linear system. The algorithm presented in Matlab reduces the computational and storage complexities and also improves the accuracy of the solution. It also provides the necessary warning about the existence of too much contradiction in the model. In addition, we suggest a thorough relook into the mathematical modeling to determine the reason why unacceptable contradiction has occurred thus prompting us to make necessary corrections/modifications to the models - both mathematical and, if necessary, physical.

Sen, Symal K.↗

Coupled rotor-body equations of motion hover flight

A set of linearized equations of motion to predict the linearized dynamic response of a single rotor helicopter in a hover trim condition to cyclic pitch control inputs is described. The equations of motion assume four fuselage degrees of freedom: lateral and longitudinal translation, roll angle, pitch angle: four rotor degrees of freedom: flapping (lateral and longitudinal tilt of the tip path plane), lagging (lateral and longitudinal displacement of the rotor plane center of mass); and dynamic inflow (harmonic components). These ten degrees of freedom correspond to a system with eighteen dynamic states. In addition to examination of the full system dynamics, the computer code supplied with this report permits the examination of various reduced order models. The code is presented in a specific form such that the dynamic response of a helicopter in flight can be investigated. With minor modifications to the code the dynamics of a rotor mounted on a flexible support can also be studied.

Curtiss, H. C., Jr.↗

A linear analysis of the transition curve for the baroclinic annulus

A primitive equation linear model is shown to be capable of quantitatively predicting the neutral curve for a baroclinic annulus. Flows unstable to small nonaxisymmetric perturbations are assumed to lie on the nonaxisymmetric side of a transition curve, while the stable flows are on the axisymmetric side. The stability in the presence of perturbations is tested with linearized hydrostatic primitive equations. Attention is given to the fastest growing eigenmodes, demonstrating that the linear and nonlinear wave have similar structures at a particular point in the nonaxisymmetric regime. All the waves are similar to the Eddy wave except at the extreme to the transition curve, where the waves maintain little structure with height. The eddy kinetic energy is baroclinic with a potential energy suppresed by the large static stability of the basic state. Eddy dissipation is noted to be significant near all boundaries.

Miller, T. L.↗

Analysis of discretization errors in LES

All numerical simulations of turbulence (DNS or LES) involve some discretization errors. The integrity of such simulations therefore depend on our ability to quantify and control such errors. In the classical literature on analysis of errors in partial differential equations, one typically studies simple linear equations (such as the wave equation or Laplace's equation). The qualitative insight gained from studying such simple situations is then used to design numerical methods for more complex problems such as the Navier-Stokes equations. Though such an approach may seem reasonable as a first approximation, it should be recognized that strongly nonlinear problems, such as turbulence, have a feature that is absent in linear problems. This feature is the simultaneous presence of a continuum of space and time scales. Thus, in an analysis of errors in the one dimensional wave equation, one may, without loss of generality, rescale the equations so that the dependent variable is always of order unity. This is not possible in the turbulence problem since the amplitudes of the Fourier modes of the velocity field have a continuous distribution. The objective of the present research is to provide some quantitative measures of numerical errors in such situations. Though the focus of this work is LES, the methods introduced here can be just as easily applied to DNS. Errors due to discretization of the time-variable are neglected for the purpose of this analysis.

Ghosal, Sandip↗

Transient Growth Analysis of Compressible Boundary Layers with Parabolized Stability Equations

The linear form of parabolized linear stability equations (PSE) is used in a variational approach to extend the previous body of results for the optimal, non-modal disturbance growth in boundary layer flows. This methodology includes the non-parallel effects associated with the spatial development of boundary layer flows. As noted in literature, the optimal initial disturbances correspond to steady counter-rotating stream-wise vortices, which subsequently lead to the formation of stream-wise-elongated structures, i.e., streaks, via a lift-up effect. The parameter space for optimal growth is extended to the hypersonic Mach number regime without any high enthalpy effects, and the effect of wall cooling is studied with particular emphasis on the role of the initial disturbance location and the value of the span-wise wavenumber that leads to the maximum energy growth up to a specified location. Unlike previous predictions that used a basic state obtained from a self-similar solution to the boundary layer equations, mean flow solutions based on the full Navier-Stokes (NS) equations are used in select cases to help account for the viscous-inviscid interaction near the leading edge of the plate and also for the weak shock wave emanating from that region. These differences in the base flow lead to an increasing reduction with Mach number in the magnitude of optimal growth relative to the predictions based on self-similar mean-flow approximation. Finally, the maximum optimal energy gain for the favorable pressure gradient boundary layer near a planar stagnation point is found to be substantially weaker than that in a zero pressure gradient Blasius boundary layer.

Compressible boundary layer↗

Linearizing the BPS equations with vector and tensor multiplets

We analyse the BPS equations of N = (1, 0) supergravity theory in six dimensions coupled to a vector and tensor multiplet. We show how these BPS equations can be reduced to a set of linear differential equations. This system is triangular in that each layer of equations, while linear, is quadratically sourced by the solutions of the previous layers. We examine several explicit examples and discuss the construction of new families of microstate geometries. We expect that the result presented here will open up new branches of superstrata in which the momentum is encoded in a new class of charge carriers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Longitudinal stability and control in wind shear with energy height rate feedback

The longitudinal linearized equations of motion in wind shear were derived for the NASA Terminal Configured Vehicle, a modified Boeing 737 airplane. In addition to the apparent acceleration terms resulting from wind shear, the equations included altitude dependent stability derivatives. A linear analysis of these equations indicates a first order divergence type of instability due to wind shear in which head wind decreased with altitude. Furthermore, this instability cannot be stabilized by attitude control alone. However, attitude control used in combination with an addition feedback loop which consisted of the energy height rate feedback to the throttle proved to be effective in suppressing instability due to wind shear. A brief piloted, real time, nonlinear simulation indicated the desirability of using a display based on the rate of change of energy height rate and of commanded thrust.

Gera, J.↗

Identification and control of structures in space

The derivation of the equations of motion for the Spacecraft Control Laboratory Experiment (SCOLE) is reported and the equations of motion of a similar structure orbiting the earth are also derived. The structure is assumed to undergo large rigid-body maneuvers and small elastic deformations. A perturbation approach is proposed whereby the quantities defining the rigid-body maneuver are assumed to be relatively large, with the elastic deformations and deviations from the rigid-body maneuver being relatively small. The perturbation equations have the form of linear equations with time-dependent coefficients. An active control technique can then be formulated to permit maneuvering of the spacecraft and simultaneously suppressing the elastic vibration.

Meirovitch, L.↗

Second-order non-iterative ADI solution of non-linear partial differential equations

A new method for the solution of non-linear partial differential equations by an ADI procedure is described. Although the method is second order accurate in time, it does not require either iterations or predictor corrector methods to overcome the nonlinearity of the equations. Thus the computational effort required for the solution of the non-linear problem becomes similar to that required for the linear case. The method is applied to a two-dimensional 'extended Burgers equation'. Linear stability is studied, and some numerical solutions obtained. The improved accuracy obtained by the 2nd order truncation error is clearly manifested.

Wolfshtein, M.↗