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 modified stepwise regression is applied to flight data from a light research air-plane operating at high angles at attack. The well-known phenomenon referred to as buckling or porpoising is analyzed and modeled using both power series and spline expansions of the aerodynamic force and moment coefficients associated with the longitudinal equations of motion.
Analysis of global magnetic anomaly maps derived from satellite data is facilitated by inversion to the equivalent magnetization in a constant thickness magnetic crust or, equivalently, by reduction to the pole. Previous inversions have proven unstable near the geomagnetic equator. The instability results from magnetic moment distributions which are admissible in the inversion solution but which make only small contribution to the computed values of anomaly field. Their admissibility in the solution could result from noisy or incomplete data or from small poorly resolved anomalies. The resulting magnetic moments are unrealistically large and oscillatory. Application of the method of principal components (e.g. eigenvalue decomposition and selective elimination of less significant eigenvectors) is proposed as a way of overcoming the instability and the method is demonstrated by applying it to the region around the Bangui anomaly in Central Africa.
Analytical and numerical methods of calculation of the energy and three dimensional EPS characteristics are reported. The angular and lateral functions of electrons in EPS have been obtained by the Landau and small angle approximations A and B and compared with earlier data. A numerical method of solution of cascade equations for the EPS distribution function moments has been constructed. Considering the equilibrium rms angle as an example, errors appearing when approximating the elementary process cross sections by their asymptotic expressions are analyzed.
A FORTRAN coded computer program and method for evaluation of the rigid body disturbing moments for a launch vehicle first stage based on post-flight measurements is described. The technique is a straightforward deterministic approach. Residual moments are computed to satisfy the equations of motion. Residuals are expressed in terms of altered vehicle characteristics; the aerodynamic coefficients, thrust misalignment, and control effectiveness. This method was used on the Scout launch vehicle and uncovered several significant differences between flight data and wind tunnel data. The computer program is written in FORTRAN IV for a CDC CYBER 173 computer system.
Mutual coupling between rectangular apertures in a finite antenna array, in an infinite ground plane, is analyzed using the vector potential approach. The method of moments is used to solve the equations that result from setting the tangential magnetic fields across each aperture equal. The approximation uses a set of vector potential model functions to solve for equivalent magnetic currents. A computer program was written to carry out this analysis and the resulting currents were used to determine the co- and cross-polarized far zone radiation patterns. Numerical results for various arrays using several modes in the approximation are presented. Results for one and two aperture arrays are compared against published data to check on the agreement of this model with previous work. Computer derived results are also compared against experimental results to test the accuracy of the model. These tests of the accuracy of the program showed that it yields valid data.
The paper discusses a design for an attitude control system for the Space Station which produces fast response, with minimal overshoot and cross-coupling with the use of Control Moment Gyros (CMG). The rigid body equations of motion are linearized and discretized and a Linear Quadratic Regulator (LQR) design and analysis study is performed. The resulting design is then modified such that integral and differential terms are added to the state equations to enhance response characteristics. Methods for reduction of computation time through channelization are discussed as well as the reduction of initial torque requirements.
Theoretical models are developed and numerical studies conducted on various types of flows including both elliptic and parabolic. The purpose of this study is to find better higher order closure models for the computations of complex flows. This report summarizes three new achievements: (1) completion of the Reynolds-stress closure by developing a new pressure-strain correlation; (2) development of a parabolic code to compute jets and wakes; and, (3) application to a flow through a 180 deg turnaround duct by adopting a boundary fitted coordinate system. In the above mentioned models near-wall models are developed for pressure-strain correlation and third-moment, and incorporated into the transport equations. This addition improved the results considerably and is recommended for future computations. A new parabolic code to solve shear flows without coordinate tranformations is developed and incorporated in this study. This code uses the structure of the finite volume method to solve the governing equations implicitly. The code was validated with the experimental results available in the literature.
Theoretical research in two areas was performed. The aerothermodynamics research focused on the hard-body and rocket plume flows. Analytical real gas models to describe finite rate chemistry were developed and incorporated into the three-dimensional flow codes. New numerical algorithms capable of treating multi-species reacting gas equations and treating flows with large gradients were also developed. The computational chemistry research focused on the determination of spectral radiative intensity factors, transport properties and reaction rates. Ab initio solutions to the Schrodinger equation provided potential energy curves transition moments (radiative probabilities and strengths) and potential energy surfaces. These surfaces were then coupled with classical particle reactive trajectories to compute reaction cross-sections and rates.
Due to the requirement of increased performance and maneuverability, the flight envelope of a modern fighter is frequently extended to the high angle-of-attack regime. Vehicles maneuvering in this regime are subjected to nonlinear aerodynamic loads. The nonlinearities are due mainly to three-dimensional separated flow and concentrated vortex flow that occur at large angles of attack. Accurate prediction of these nonlinear airloads is of great importance in the analysis of a vehicle's flight motion and in the design of its flight control system. A satisfactory evaluation of the performance envelope of the aircraft may require a large number of coupled computations, one for each change in initial conditions. To avoid the disadvantage of solving the coupled flow-field equations and aircraft's motion equations, an alternate approach is to use a mathematical modeling to describe the steady and unsteady aerodynamics for the aircraft equations of motion. Aerodynamic forces and moments acting on a rapidly maneuvering aircraft are, in general, nonlinear functions of motion variables, their time rate of change, and the history of maneuvering. A numerical method was developed to analyze the nonlinear and time-dependent aerodynamic response to establish the generalized indicial function in terms of motion variables and their time rates of change.
A selection of successes and failures of Computational Fluid Dynamics (CFD) is discussed. Experiment/CFD correlations involving full potential and Euler computations of the aerodynamic characteristics of four commercial transport wings and two low aspect ratio, delta wing configurations are shown. The examples consist of experiment/CFD comparisons for aerodynamic forces, moments, and pressures. Navier-Stokes equations are not considered.
Discussed here is the future direction of various levels of engineering turbulence modeling related to computational fluid dynamics (CFD) computations for propulsion. For each level of computation, there are a few turbulence models which represent the state-of-the-art for that level. However, it is important to know their capabilities as well as their deficiencies in order to help engineers select and implement the appropriate models in their real world engineering calculations. This will also help turbulence modelers perceive the future directions for improving turbulence models. The focus is on one-point closure models (i.e., from algebraic models to higher order moment closure schemes and partial differential equation methods) which can be applied to CFD computations. However, other schemes helpful in developing one-point closure models, are also discussed.
Adaptive linear predictors have been used extensively in practice in a wide variety of forms. In the main, their theoretical development is based upon the assumption of stationarity of the signals involved, particularly with respect to the second order statistics. On this basis, the well-known normal equations can be formulated. If high- order statistical stationarity is assumed, then the equivalent normal equations involve high-order signal moments. In either case, the cross moments (second or higher) are needed. This renders the adaptive prediction procedure non-blind. A novel procedure for blind adaptive prediction has been proposed and considerable implementation has been made in our contributions in the past year. The approach is based upon a suitable interpretation of blind equalization methods that satisfy the constant modulus property and offers significant deviations from the standard prediction methods. These blind adaptive algorithms are derived by formulating Lagrange equivalents from mechanisms of constrained optimization. In this report, other new update algorithms are derived from the fundamental concepts of advanced system identification to carry out the proposed blind adaptive prediction. The results of the work can be extended to a number of control-related problems, such as disturbance identification. The basic principles are outlined in this report and differences from other existing methods are discussed. The applications implemented are speech processing, such as coding and synthesis. Simulations are included to verify the novel modelling method.
Aromatic polyimides are used as matrix resins in advanced composites, as high strength films, and as high-temperature adhesives, owing in part to their unusual thermal and chemical stability. The polyimides' desirable qualities of very high softening temperatures and negligibly small solubilities in and low permeabilities by most solvents have limited the kinds of fundamental studies that can be performed on these systems. Consequently, relationships between the molecular structure of polyimides and their bulk properties are not as well understood as might be expected given their widespread applications. In particular, the intermolecular forces in polyimides that play a critical role determining their densities, solubilities, viscosities, moduli, glass transitions, etc. are less well characterized for polyimides than for other widely used polymeric materials. The purpose of the present study is to obtain experimental data for establishing parameters that characterize the intermolecular forces in polyimides. We report here our studies on tractable low molecular-weight imides that contain the same structural features that are present in polyimide materials. We have measured equation-of-state properties and dipole moments for a variety of such systems in the liquid state. Both pure compounds and binary mixtures have been studied.
In radiometer applications, it is required to design antennas that meet low average sidelobe levels and low average return loss over a specified frequency bandwidth. It is a challenge to meet such specifications over a frequency range when one uses resonant elements such as waveguide feed slots. In addition to their inherent narrow frequency band performance, the problem is exacerbated due to modeling errors and manufacturing tolerances. There was a need to develop a design methodology to solve the problem. An iterative design procedure was developed by starting with an array architecture, lattice spacing, aperture distribution, waveguide dimensions, etc. The array was designed using Elliott s technique with appropriate values of the total slot conductance in each radiating waveguide, and the total resistance in each feed waveguide. Subsequently, the array performance was analyzed by the full wave method of moments solution to the pertinent integral equations. Monte Carlo simulations were also carried out to account for amplitude and phase errors introduced for the aperture distribution due to modeling errors as well as manufacturing tolerances. If the design margins for the average sidelobe level and the average return loss were not adequate, array architecture, lattice spacing, aperture distribution, and waveguide dimensions were varied in subsequent iterations. Once the design margins were found to be adequate, the iteration was stopped and a good design was achieved. A symmetric array architecture was found to meet the design specification with adequate margin. The specifications were near 40 dB for angular regions beyond 30 degrees from broadside. Separable Taylor distribution with nbar=4 and 35 dB sidelobe specification was chosen for each principal plane. A non-separable distribution obtained by the genetic algorithm was found to have similar characteristics. The element spacing was obtained to provide the required beamwidth and close to a null in the E-plane end-fire direction. Because of the alternating slot offsets, grating lobes called butterfly lobes are produced in non-principal planes close to the H-plane. An attempt to reduce the influence of such grating lobes resulted in a symmetric design.
Governing equations of motion are derived for a bulk-flow model of the leakage path between an impeller shroud and a pump housing. The governing equations consist of a path-momentum, a circumferential - momentum, and a continuity equation. The fluid annulus between the impeller shroud and pump housing is assumed to be circumferentially symmetric when the impeller is centered; i.e., the clearance can vary along the pump axis but does not vary in the circumferential direction. A perturbation expansion of the governing equations in the eccentricity ratio yields a set of zeroth and first-order governing equations. The zeroth-order equations define the leaking rate and the circumferential and path velocity distributions and pressure distributions for a centered impeller position. The first-order equations define the perturbations in the velocity and pressure distributions due to either a radial-displacement perturbation or a tilt perturbation of the impeller. Integration of the perturbed pressure and shear-stress distribution acting on the rotor yields the reaction forces and moments acting on the impeller face.
The accurate modeling of electromagnetic penetration is an important topic in computational electromagnetics. Electromagnetic penetration occurs through intentional or inadvertent openings in an otherwise closed electromagnetic scatterer, which prevent the contents from being fully shielded from external fields. To efficiently model electromagnetic penetration, aperture or slot models can be used with surface integral equations to solve Maxwell's equations. A necessary step towards establishing the credibility of these models is to assess the correctness of the implementation of the underlying numerical methods through code verification. Surface integral equations and slot models yield multiple interacting sources of numerical error and other challenges, which render traditional code-verification approaches ineffective. In this paper, we provide approaches to separately measure the numerical errors arising from these different error sources for the method-of-moments implementation of the electric-field integral equation with a slot model. Finally, we demonstrate the effectiveness of these approaches for a variety of cases.
A general method for finding the steady flow velocity relative to a body in plane curvilinear motion, whence the pressure is found by Bernoulli's energy principle is described. Integration of the pressure supplies basic formulas for the zonal forces and moments on the revolving body. The application of the steady flow method for calculating the velocity and pressure at all points of the flow inside and outside an ellipsoid and some of its limiting forms is presented and graphs those quantities for the latter forms. In some useful cases experimental pressures are plotted for comparison with theoretical. The pressure, and thence the zonal force and moment, on hulls in plane curvilinear flight are calculated. General equations for the resultant fluid forces and moments on trisymmetrical bodies moving through a perfect fluid are derived. Formulas for potential coefficients and inertia coefficients for an ellipsoid and its limiting forms are presented.
An implementation of the Asymptotic Waveform Evaluation (AWE) technique is presented for obtaining the frequency response of the Radar Cross Section (RCS) of arbitrarily shaped, three-dimensional perfect electric conductor (PEC) bodies. An Electric Field Integral Equation (EFIE) is solved using the Method of Moments (MoM) to compute the RCS. The electric current, thus obtained, is expanded in a Taylor series around the frequency of interest. The coefficients of the Taylor series (called 'moments') are obtained using the frequency derivatives of the EFIE. Using the moments, the electric current on the PEC body is obtained over a frequency band. Using the electric current at different frequencies, RCS of the PEC body is obtained over a wide frequency band. Numerical results for a square plate, a cube, and a sphere are presented over a bandwidth. A good agreement between AWE and the exact solution over the bandwidth is observed.