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.
Asymmetrical nonlinear oscillations of control systems analyzed by parameter plane method, considering systems with reference-perturbing and slow-varying input signals
Satellite orbit perturbation analysis for lunar gravitational field determination using von Zeipel transformation
Three body stellar problem libration calculation using nonlinear mechanics methods, and application to lunar satellite perturbation by earth and lunar gravitational effects
The development of improved computer algorithms is considered for calculating earth satellite orbital trajectories by optimum selection of the analytical method that minimizes the number of perturbative acceleration computations for a given accuracy. A variation of parameter algorithm considering the equation of motion for a satellite proved superior for the geosynchronous orbit.
The spherical harmonic coefficients in the gravity field of Mars are determined using a two stage filtering process. In the first stage, a weighted least squares orbit determination processor is applied to Mariner 9 Doppler tracking data to estimate Kepler elements and rates. In the second stage, the Kepler element rates are used as input to another weighted least squares processor that estimates gravity coefficients using the Lagrange perturbation equations. The novelty in this method is the fitting of element rates rather than the elements themselves. This procedure decouples the results from the effects of internal spacecraft dynamics.
Calculations for transonic flow about wing-body combinations at angles of attack of the order of the body thickness-to-length ratio are presented. These calculations show that lift can alter the flow drastically at near-sonic speeds in this angle-of-attack range so that the area rule must be modified. This effect of lift is explained from simple physical considerations. The computational procedure is based on slender-wing theory and a two-variable method of relaxation solution. A small perturbation analysis is used to show that this simple procedure is adequate for these transonic wing-body flow fields.
The techniques considered are related to gravity profiling from Doppler residuals, a solution for a surface mass distribution, and a dynamical approach for the estimation of random forces at each data point. The dynamical approach makes use of a new filtering method which uses sequential estimation theory to determine unknown and somewhat correlated forces which randomly perturb the orbit. It is concluded that lunar data reduction methods would be suitable for applications in connection with the GEOS-C/ATS-F mission.
Nonlinear equations of motion for a cantilever rotor blade are derived for the hovering flight condition. The blade is assumed to have twist, precone, droop, sweep, torque offset and blade root offset, and the elastic axis and the axes of center of mass, tension, and aerodynamic center coincident at the quarter chord. The blade is cantilevered in bending, but has a torsional root spring to simulate pitch link flexibility. Aerodynamic forces acting on the blade are derived from strip theory based on quasi-steady two-dimensional airfoil theory. The equations are hybrid, consisting of one integro-differential equation for root torsion and three integro-partial differential equations for flatwise and chordwise bending and elastic torsion. The equations are specialized for a uniform blade and reduced to nonlinear ordinary differential equations by Galerkin's method. They are linearized for small perturbation motions about the equilibrium operating condition. Modal analysis leads to formulation of a standard eigenvalue problem where the elements of the stability matrix depend on the solution of the equilibrium equations. Two different forms of the root torsion equation are derived that yield virtually identical numerical results. This provides a reasonable check for the accuracy of the equations.
The 'mean frequency' technique, a simple procedure introduced by Bebb and Gold for the approximate evaluation of sums occurring in high-order perturbation theory, represents a useful approximation method. Its predictions compare favorably to exact results obtained by Gontier and Trahin for multiphoton bound-bound transitions in hydrogen. However, the technique can be in error if the 'mean frequency' lies near certain integers.
A two-dimensional nonlinear theory of sound transmission through a nonuniform duct carrying a near-sonic throat flow is described. The mean flow in the duct is treated according to a generalized quasi-one dimensional model, and the analysis of the unsteady perturbations is carried out using the Method of Matched Asymptotic Expansions. The linearized acoustic field in the subsonic regions of the duct is approximated by a Wave Envelope model which is matched asymptotically with a nonlinear inner solution valid in the near-sonic throat region. Numerical results are presented to illustrate the predictions of the theory. It is found that, in general, shock waves develop in the acoustic field during transmission through the near-sonic flow. Unlike previous one-dimensional theories, the current study shows that dispersion can play a major role in the propagation process.
The applicability of current finite-volume CFD algorithms based on the Euler equations to the vortex flow over a helicopter in forward flight is investigated analytically. The general characteristics of the flow are reviewed; existing Euler, Navier-Stokes, perturbation, high-order, and adaptive methods are briefly characterized; and a novel Eulerian/Lagrangian approach with entropy and vorticity corrections is presented in detail. Numerical results for simple convection of a finite-core Lamb vortex moving downstream with its axis perpendicular to the flow are presented in graphs, and the possibility of extending the method to three-dimensional, viscous, and shock flows is discussed.
This paper describes the results of a study to evaluate the orbit determinatioin of Tracking and Data Relay Satellite System (TDRSS) user spacecraft within the dual-Tracking and Data Relay Satellite (TDRS) environment. Dense TDRSS tracking of the Earth Radiation Budget Satellite (ERBS) was acquired for the period August 16 through 22, 1989. This tracking information was processed to evaluate the orbit determination consistency achieved using the Goddard Trajectory Determination System batch least-squares estimator. The effects of the use of the second operational relay spacecraft, of refinements in orbit determination models (geopotentials, polar motion, solid earth tidal gravitational perturbations, ionospheric refraction corrections), and of methods for providing relay spacecraft spacecraft position information were also studied.
The minimization of inviscid fluid drag is studied for aerodynamic shapes satisfying the conditions of linearized theory, and subject to imposed constraints on lift, pitching moment, base area, or volume. The problem is transformed to one of determining two-dimensional potential flows satisfying either Laplace's or Poisson's equations with boundary values fixed by the imposed conditions. A general method for determining integral relations between perturbation velocity components is developed. This analysis is not restricted in application to optimum cases; it may be used for any supersonic wing problem.
A low fidelity simulation approach is used to explore how to place and use noise monitors during X-59 QueSST community tests, where people’s annoyance to the noise produced by the X-59 aircraft will be gathered. Several recommendations are provided including: 1) the desired number of sparsely spaced noise monitor sites within the survey area, 2) whether to group and average measurements across multiple noise monitors located at a site, 3) what spacing should be used if grouped noise monitors are used, 4) an approach to mitigate ambient noise contamination at the measurement sites, 5) a method to combine empirical and predicted dose estimates to provide a single dose estimate for respondents, and 6) assessing how changes in turbulence intensity and array configuration affect dose uncertainty. To make these recommendations, the error that is expected when fitting contrived, smoothly varying sonic boom “reference exposure surfaces” is studied when a spatially sparse and scattered set of samples is used as responses for the fit. The reference exposure surfaces mimic the sonic boom exposure at ground level that might be expected in the X-59 survey area in the absence of atmospheric turbulence, ambient noise, and other localized effects. The spatial extent of these surfaces varies and is representative of the different survey area sizes that might be expected during future X-59 community overflight tests. These contrived reference surfaces are sampled, and those reference samples are then perturbed to mimic atmospheric turbulence, ambient noise and other localized effects that might affect noise monitor measurements within overflown communities. Two different surface fitting methods are investigated when fitting these perturbed samples to approximate the reference surface. The first method uses interpolation between the perturbed data at the scattered sites to compute the fit. The second method fits a polynomial surface model to the perturbed data using ordinary least squares regression analysis. For both fitting methods, the root mean square fit error is computed from the pointwise difference between the fit surface and the reference surface as the count and configuration of the sites is varied while also averaging the error across many different realizations of both the smooth variation of the reference exposure surface and the random, localized perturbations at the sample sites. Different site configurations are compared using this error statistic to make the recommendations noted above. Additionally, the two fitting approaches (interpolation vs linear regression) are compared based on the fit error observed in these simulations. These analyses, comparisons, and recommendations should inform future decisions on the noise monitor placement and the methods used to analyze the noise monitor data that is collected during X-59 community overflights.
We examine the nonlinear stability of the Wisdom-Holman (WH) symplectic mapping applied to the integration of perturbed, highly eccentric (e-0.9) two-body orbits. We find that the method is unstable and introduces artificial chaos into the computed trajectories for this class of problems, unless the step size chosen 1s small enough that PeriaPse is always resolved, in which case the method is generically stable. This 'radial orbit instability' persists even for weakly perturbed systems. Using the Stark problem as a fiducial test case, we investigate the dynamical origin of this instability and argue that the numerical chaos results from the overlap of step-size resonances; interestingly, for the Stark-problem many of these resonances appear to be absolutely stable. We similarly examine the robustness of several alternative integration methods: a time-regularized version of the WH mapping suggested by Mikkola; the potential-splitting (PS) method of Duncan, Levison, Lee; and two original methods incorporating approximations based on Stark motion instead of Keplerian motion. The two fixed point problem and a related, more general problem are used to conduct a comparative test of the various methods for several types of motion. Among the algorithms tested, the time-transformed WH mapping is clearly the most efficient and stable method of integrating eccentric, nearly Keplerian orbits in the absence of close encounters. For test particles subject to both high eccentricities and very close encounters, we find an enhanced version of the PS method-incorporating time regularization, force-center switching, and an improved kernel function-to be both economical and highly versatile. We conclude that Stark-based methods are of marginal utility in N-body type integrations. Additional implications for the symplectic integration of N-body systems are discussed.
In this paper, the development of a digital three-dimensional automatic control law designated to achieve an optimal transition of a large commercial aircraft between the glideslope conditions and the desired final touchdown condition is presented. The linearized equations representing the perturbed motion of the aircraft from the nominal glideslope trajectory are presented. A method of incorporating the spatial, low-level wind shears into these perturbed equations of motion is indicated. It is shown that the system equations then assume the familiar form of the linear regulator problem, acted upon by a constant disturbance. A design procedure is presented to compute a digital, time-invariant, optimal control law for the discrete regulator problem acted upon by a constant disturbance. The response trajectories presented demonstrate the capability of this controller to perform optimal flare maneuvers in the presence of various wind shear conditions.
Brouwers theory of planetary motion for second- order perturbation computation using Fourier series
An analytical solution to the perturbative multiple collision series of a fragmenting HZE ion beam has limited usefulness since the first collision term has several hundred contributions, the second collision term has tens of thousands of contributions, and each successive collision term progresses to unwieldy computational proportions. Our previous work has revealed the multiple collision terms in the straight-ahead approximation to be simple products of a spatially dependent factor times a linear energy-dependent factor of limited domain and unit normalization. The properties of these forms allow the development of the nonperturbative summation of the series to all orders assuming energy-independent nuclear cross sections as matrix products of a scaled Green's function described herein. This nonperturbative Green's function with multiple scattering correction factors compares well with experiments using 670 MeV/u neon-20 ion beams in thick water targets.