SEARCH · Search NASA
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 study has been made of a dynamical system composed of a pendulum and a harmonic oscillator, in order to show the remarkable resemblance with many classical celestial mechanics problems, in particular, the restricted three-body problem. It is shown that the well-known investigations of periodic orbits can be applied to the present dynamics problem.
This report describes a computer program (program D2180) that has been prepared to implement the analysis described in (N71-10866) for calculating the aerodynamic forces on a class of harmonically oscillating planar lifting surfaces in supersonic potential flow. The planforms treated are the delta and modified-delta (arrowhead) planforms with subsonic leading and supersonic trailing edges, and (essentially) pointed tips. The resulting aerodynamic forces are applied in a Galerkin modal flutter analysis. The required input data are the flow and planform parameters including deflection-mode data, modal frequencies, and generalized masses.
The Newtonian differential equations of motion for the two-body problem can be transformed into four linear harmonic-oscillator equations by simultaneously applying the regularization step dt/ds = r and the Kustaanheimo-Stieffel (KS) transformation. The regularization step changes the independent variable from time to a new variable s, and the KS transformation transforms the position and velocity vectors from Cartesian space into a four-dimensional space. A derivation of a uniform, regular solution for the perturbed two-body problem in the four-dimensional space is presented. The variation-of-parameters technique is used to develop expressions for the derivatives of ten elements (which are constants in the unperturbed motion) for the general case that includes both perturbations which can arise from a potential and perturbations which cannot be derived from a potential. This ten-element solution has mixed secular terms that degrade the long-term accuracy during numerical integration. Therefore, to eliminate these terms, the solution is modified by introducing two additional elements.
The general theory of potential aerodynamic flow around a lifting body having arbitrary shape and motion is presented. By using the Green function method, an integral representation for the potential is obtained for both supersonic and subsonic flow. Under small perturbation assumption, the potential at any point, P, in the field depends only upon the values of the potential and its normal derivative on the surface, sigma, of the body. Hence, if the point P approaches the surface of the body, the representation reduces to an integro-differential equation relating the potential and its normal derivative (which is known from the boundary conditions) on the surface sigma. For the important practical case of small harmonic oscillation around a rest position, the equation reduces to a two-dimensional Fredholm integral equation of second-type. It is shown that this equation reduces properly to the lifting surface theories as well as other classical mathematical formulas. The question of uniqueness is examined and it is shown that, for thin wings, the operator becomes singular as the thickness approaches zero. This fact may yield numerical problems for very thin wings.
The optimal design was investigated of simple structures subjected to dynamic loads, with constraints on the structures' responses. Optimal designs were examined for one dimensional structures excited by harmonically oscillating loads, similar structures excited by white noise, and a wing in the presence of continuous atmospheric turbulence. The first has constraints on the maximum allowable stress while the last two place bounds on the probability of failure of the structure. Approximations were made to replace the time parameter with a frequency parameter. For the first problem, this involved the steady state response, and in the remaining cases, power spectral techniques were employed to find the root mean square values of the responses. Optimal solutions were found by using computer algorithms which combined finite elements methods with optimization techniques based on mathematical programming. It was found that the inertial loads for these dynamic problems result in optimal structures that are radically different from those obtained for structures loaded statically by forces of comparable magnitude.
A semiclassical model proposed by Pearson and Hansen (1974) for computing collision-induced transition probabilities in diatomic molecules is tested by the direct-simulation Monte Carlo method. Specifically, this model is described by point centers of repulsion for collision dynamics, and the resulting classical trajectories are used in conjunction with the Schroedinger equation for a rigid-rotator harmonic oscillator to compute the rotational energy transition probabilities necessary to evaluate the rotation-translation exchange phenomena. It is assumed that a single, average energy spacing exists between the initial state and possible final states for a given collision.
The computational procedure and numerical results are presented for a new method to solve Kuessner's integral equation in the case of subsonic compressible flow about harmonically oscillating planar surfaces with controls. Kuessner's equation is a linear transformation from pressure to normalwash. The unknown pressure is expanded in terms of prescribed basis functions and the unknown basis function coefficients are determined in the usual manner by satisfying the given normalwash distribution either collocationally or in the complex least squares sense. The present method of solution differs from previous ones in that the basis functions are defined in a continuous fashion over a relatively small portion of the aerodynamic surface and are zero elsewhere. This method, termed the local basis function method, combines the smoothness and accuracy of distribution methods with the simplicity and versatility of panel methods. Predictions by the local basis function method for unsteady flow are shown to be in excellent agreement with other methods. Also, potential improvements to the present method and extensions to more general classes of solutions are discussed.
The general behavior of unsteady airloads in the frequency domain is explained. Based on this, a systematic procedure is described whereby the airloads, produced by completely arbitrary, small, time-dependent motions of a thin lifting surface in an airstream, can be predicted. This scheme employs as raw materials any of the unsteady linearized theories that have been mechanized for simple harmonic oscillations. Each desired aerodynamic transfer function is approximated by means of an appropriate Pade approximant, that is, a rational function of finite degree polynomials in the Laplace transform variable. Although these approximations have many uses, they are proving especially valuable in the design of automatic control systems intended to modify aeroelastic behavior.
The paper examines some questions relating to the optimal design of continuous one-dimensional structures driven by harmonically oscillating loads. Optimal-control methods are applied to a cantilever bar driven sinusoidally by an axial force at its tip to illustrate the minimum-weight design of one-dimensional structures under dynamic excitation. Realistic constraints are imposed during the optimizations, including a maximum allowable stress amplitude at any point along the bar and a minimum cross-sectional area. It is shown that in the absence of damping, the design space may contain many disjoint feasible regions, and multiple optima can exist. Detailed solutions are obtained for continuous bars with an excitation frequency less than, and then greater than, the fundamental free-vibration frequency. It is found that above a certain excitation frequency, two or more arcs with different constraints characterize the optimal designs. It is concluded that when more than two different constrained arcs characterize the optimal solution, the continuum approach may be impractical, and finite-element approximations may offer the only alternative.
A general theory of finite state modeling of aerodynamic loads on thin airfoils and lifting surfaces performing completely arbitrary, small, time-dependent motions in an airstream is developed and presented. The nature of the behavior of the unsteady airloads in the frequency domain is explained, using as raw materials any of the unsteady linearized theories that have been mechanized for simple harmonic oscillations. Each desired aerodynamic transfer function is approximated by means of an appropriate Pade approximant, that is, a rational function of finite degree polynomials in the Laplace transform variable. The modeling technique is applied to several two dimensional and three dimensional airfoils. Circular, elliptic, rectangular and tapered planforms are considered as examples. Identical functions are also obtained for control surfaces for two and three dimensional airfoils.
Results indicating that unsteady aerodynamic loads derived under the assumption of simple harmonic motions executed by airfoil or wing can be extended to arbitrary motions are summarized. The generalized Theodorsen (1953) function referable to loads due to simple harmonic oscillations of a wing section in incompressible flow, the Laplace inversion integral for unsteady aerodynamic loads, calculations of root loci of aeroelastic loads, and analysis of generalized compressible transient airloads are discussed.
The time-optimal control of a spin-stabilized spacecraft with a movable telescoping appendage (boom) is considered analytically and numerically. The motion of a control mass at the end of the boom is determined such that the terminal time will be minimized for two-axis control of a symmetric spacecraft. The equations of rotational motion are linearized about the desired state of spin about the symmetry axis. The equations for the transverse angular velocity components have the form of a coupled two dimensional harmonic oscillator with boom motion as a control force. The control function which brings the system to the desired state is known to be a series of positive and negative pulses. If the initial state is such that the system can be driven to rest in a single switch, the responses, switching and final times, and required boom motion may be determined analytically. Some typical numerical results based on these solutions are discussed.
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.
The limitations of the detectability of extremely weak signals (gravitational radiation for instance) imposed by Heisenberg's uncertainty principle on the sequential determination of those signals have been explored recently. A variety of schemes have been proposed to circumvent these limitations. Although all of the earlier attempts have been proven fruitless a recent proposal seems to be quite promising. The scheme, consisting of two harmonic oscillators interacting with each other in a peculiar way, allows for an exact analytical solution which is derived here. If it can be assumed that the expectation value of one of the canonical variables of the total system suffices to monitor the weak signal it can be shown that, in the absence of thermal noise, arbitrarily weak signals can in principle be measured without interference from the uncertainty principle.
A two dimensional cascade of harmonically oscillating airfoils was designed to model a near tip section from a rotor which was known to have experienced supersonic torsional flutter. This five bladed cascade had a solidity of 1.17 and a setting angle of 1.07 rad. Graphite epoxy airfoils were fabricated to achieve the realistically high reduced frequency level of 0.44. The cascade was tested over a range of static pressure ratios approximating the blade element operating conditions of the rotor along a constant speed line which penetrated the flutter boundary. The time-steady and time-unsteady flow field surrounding the center cascade airfoil were investigated. The effects of reduced solidity and decreased setting angle on the flow field were also evaluated.
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.
An advanced method is presented for solving the linear integral equations for subsonic unsteady flow in three dimensions. The method is applicable to flows about arbitrary, nonplanar boundary surfaces undergoing small amplitude harmonic oscillations about their steady mean locations. The problem is formulated with a wake model wherein unsteady vorticity can be convected by the steady mean component of flow. The geometric location of the unsteady source and doublet distributions can be located on the actual surfaces of thick bodies in their steady mean locations. The method is an outgrowth of a recently developed steady flow panel method and employs the linear source and quadratic doublet splines of that method.
Use has been made of flow graphs to deduce Feenberg perturbation expansions for radiative interactions. It is demonstrated that these expansions can in certain cases be summed to provide closed form expressions for the molecular response. In particular, it is shown that the coherent state response can be obtained by the summation of a continued fraction perturbation expansion for the harmonic oscillator. Anharmonicity in the lower levels is treated and its shown to introduce Rabi flopping identifiable with multiphoton transitions among isolated tightly coupled subsystems of levels. Relevance to laser induced multiphoton excitation and energy level shift calculations in the presence of a strong field are also discussed.