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.
A method using a transonic small disturbance code with successive line over-relaxation is described for treating wing/fuselage configurations with a nacelle/pylon/powered jet. Examples illustrating its use for the NASA transport research model are given. Reasonable test/theory comparisons were obtained.
A power series expansion in the damping parameter, epsilon, of the limit cycle of the free van der Pol equation is constructed and analyzed. Coefficients in the expansion are computed in exact rational arithmetic using the symbolic manipulation system MACSYMA and using a FORTRAN program. The series is analyzed using Pade approximants. The convergence of the series for the maximum amplitude of the limit cycle is limited by two pair of complex conjugate singularities in the complex epsilon-plane. A new expansion parameter is introduced which maps these singularities to infinity and leads to a new expansion for the amplitude which converges for all real values of epsilon. Amplitudes computed from this transformed series agree very well with reported numerical and asymptotic results. For the limit cycle itself, convergence of the series expansion is limited by three pair of complex conjugate branch point singularities. Two pair remain fixed throughout the cycle, and correspond to the singularities found in the maximum amplitude series, while the third pair moves in the epsilon-plane as a function of t from one of the fixed pairs to the other. The limit cycle series is transformed using a new expansion parameter, which leads to a new series that converges for larger values of epsilon.
The application of linear optimal control to the design of systems with integral control action on specified outputs is considered. Using integral terms in a quadratic performance index, an asymptotic analysis is used to determine the effect of variable quadratic weights on the eigenvalues and eigenvectors of the closed loop system. It is shown that for small integral terms the placement of integrator poles and gain calculation can be effectively decoupled from placement of the primary system eigenvalues. This technique is applied to the design of integral controls for a STOL aircraft outer loop guidance system.
A power series expansion in the damping parameter, epsilon, of the limit cycle of the free van der Pol equation is constructed and analyzed. Coefficients in the expansion are computed in exact rational arithmetic using the symbolic manipulation system MACSYMA and using a FORTRAN program. The series is analyzed using Padeapproximants. The convergence of the series for the maximum amplitude of the limit cycle is limited by two pair of complex conjugate singularities in the complex epsilon-plane. A new expansion parameter is introduced which maps these singularities to infinity and leads to a new expansion for the amplitude which converges for all real values of epsilon. Amplitudes computed from this transformed series agree very well with reported numerical and asymptotic results. For the limit cycle itself, convergence of the series expansion is limited by three pair of complex conjugate branch point singularities. Two pair remain fixed throughout the cycle, and correspond to the singularities found in the maximum amplitude series, while the third pair moves in the epsilon-plane as a function of t from one of the fixed pairs to the other. The limit cycle series is transformed using a new expansion parameter, which leads to a new series that converges for larger values of epsilon.
Perturbation analysis for angular orbital elements of Relay 2 satellite
Perturbation solution analysis of vehicle transfer from one planetary circular orbit to circular orbit around another planet using low thrust acceleration
Rounding error bounds by perturbation condition analysis, and conditioning for eigenvalue problem
Satellite orbit perturbation analysis for lunar gravitational field determination using von Zeipel transformation
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.
Analysis of the responses of rising or falling spherical wind sensors to atmospheric wind perturbations on the wind profile in the vertical, using Fourier transform techniques. The linearized equations of motion of a sensor that is subject to drag and gravitational body forces are developed by perturbing a sensor about an equilibrium uniform motion with wind fluctuations which have vertical variations. The wind environment and sensor velocities are decomposed with stochastic Fourier-Stieltjes integrals, and the linearized equations of motion are used to derive the response functions and phase angles of the sensor motions. The results of the analysis are used to analyze the response properties of the Jimsphere balloon wind sensor. It is shown that, in general, the transfer functions associated with the horizontal sensor motions are smaller than the transfer functions associated with the vertical sensor motions in the omega times T product range from zero to infinity, omega being the wind perturbation frequency and T a time constant of the system. Thus, the sensor is more responsive to vertical than to horizontal air motions.
The aerodynamic characteristics of potential flow past an axisymmetric slender body and a thin airfoil are calculated using a uniform perturbation analysis method. The method is based on the superposition of potentials of point singularities distributed inside the body. The strength distribution satisfies a linear integral equation by enforcing the flow tangency condition on the surface of the body. The complete uniform asymptotic expansion of its solution is obtained with respect to the slenderness ratio by modifying and adapting an existing technique. Results calculated by the perturbation analysis method are compared with the existing surface singularity panel method and some available analytical solutions for a number of cases under identical conditions. From these comparisons, it is found that the perturbation analysis method can provide quite accurate results for bodies with small slenderness ratio. The present method is much simpler and requires less memory and computation time than existing surface singularity panel methods of comparable accuracy.
A model of the geopotential field in spherical harmonics to degree and order 30 is obtained from Geos 3 satellite to sea surface altimetry data, terrestrial gravity measurements and satellite perturbation analysis. A general perturbation solution is employed for the calculation of the orbits of 10 satellites based on satellite laser ranging data, and 1 deg x 1 deg surface gravity data are used to compute 550 km x 550 km block anomalies by means of autocovariance analysis. Altimeter-determined sea-surface heights, which are taken as the geoid, are averaged for each 1 deg x 1 deg ocean surface area and treated by autocovariance analysis to obtain 550 x 550 km block undulations. Observation and normal equations are formed from the altimeter and surface gravity data, which together cover 1635 out of 1654 possible surface elements, and are combined with the available satellite-derived normal equations to obtain a solution for the spherical harmonics coefficients. In addition, a value of 6,378,138.23 + or - 1.3 m is obtained for the earth's semimajor axis.
Information on sensitivity analysis in computational aerodynamics is given in outline, graphical, and chart form. The prediction accuracy if the MCAERO program, a perturbation analysis method, is discussed. A procedure for calculating perturbation matrix, baseline wing paneling for perturbation analysis test cases and applications of an inviscid sensitivity matrix are among the topics covered.
The objective of this research effort was to develop a real-time guidance approach for launch vehicles ascent to orbit injection. Various analytical approaches combined with a variety of model order and model complexity reduction have been investigated. Singular perturbation methods were first attempted and found to be unsatisfactory. The second approach based on regular perturbation analysis was subsequently investigated. It also fails because the aerodynamic effects (ignored in the zero order solution) are too large to be treated as perturbations. Therefore, the study demonstrates that perturbation methods alone (both regular and singular perturbations) are inadequate for use in developing a guidance algorithm for the atmospheric flight phase of a launch vehicle. During a second phase of the research effort, a hybrid analytic/numerical approach was developed and evaluated. The approach combines the numerical methods of collocation and the analytical method of regular perturbations. The concept of choosing intelligent interpolating functions is also introduced. Regular perturbation analysis allows the use of a crude representation for the collocation solution, and intelligent interpolating functions further reduce the number of elements without sacrificing the approximation accuracy. As a result, the combined method forms a powerful tool for solving real-time optimal control problems. Details of the approach are illustrated in a fourth order nonlinear example. The hybrid approach is then applied to the launch vehicle problem. The collocation solution is derived from a bilinear tangent steering law, and results in a guidance solution for the entire flight regime that includes both atmospheric and exoatmospheric flight phases.
A method to calculate the dynamic coefficients of seals with arbitrary geometry is presented. The Navier-Stokes equations are used in conjunction with the k-e turbulence model to describe the turbulent flow. These equations are solved by a full 3-dimensional finite-difference procedure instead of the normally used perturbation analysis. The time dependence of the equations is introduced by working with a coordinate system rotating with the precession frequency of the shaft. The results of this theory are compared with coefficients calculated by a perturbation analysis and with experimental results.
This paper proposes a systematic approach for identifying the perturbation parameter in singular perturbation analysis of aircraft optimal guidance, and in particular considers a family of problems related to aircraft energy climbs. The approach, which is based on a nondimensionalization of the equations of motion, is used to evaluatae the appropriateness of forced singular perturbation formulations used in the past for transport and fighter aircraft, and to assess the applicability of energy state approximations and singular perturbation analysis for airbreathing transatmospheric vehicles with hypersonic cruise and orbital capabilities.
The paper considers a thin elastic bidirectional symmetric composite laminate subjected to a uniform axial strain. All interlaminar stresses are determined by a zeroth-order perturbation analysis of the governing equations. The method of solution considers the laminate to be composed of an interior region (removed from the free edge) and a boundary layer region. The interior region is analyzed by dropping terms related to powers of the thickness-to-width ratio in the dimensionless elasticity equations. The boundary layer region is studied by introducing a stretching transformation to the governing equations. A uniform solution is then formed by satisfying Prandtl's matching principle of perturbation theory. Results for a four-layer bidirectional graphite-epoxy laminate are compared with existing finite-difference results. It is shown that the perturbation analysis provides more insight into the stress distribution near the free edge.
The deviation of the solution of the differential equation y' = f(t, y), y(O) = y sub O from the solution of the perturbed system z' = f(t, z) + g(t, z), z(O) = z sub O was investigated for the case where f and g are continuous functions on I x R sup n into R sup n, where I = (o, a) or I = (o, infinity). These functions are assumed to satisfy the Lipschitz condition in the variable z. The space Lip(I) of all such functions with suitable norms forms a Banach space. By introducing a suitable norm in the space of continuous functions C(I), introducing the problem can be reduced to an equivalent problem in terminology of operators in such spaces. A theorem on existence and uniqueness of the solution is presented by means of Banach space technique. Norm estimates on the rate of growth of such solutions are found. As a consequence, estimates of deviation of a solution due to perturbation are obtained. Continuity of the solution on the initial data and on the perturbation is established. A nonlinear perturbation of the harmonic oscillator is considered a perturbation of equations of the restricted three body problem linearized at libration point.