An improved method for calculating supersonic pressure fields about bodies of revolution
Development of near-field method for determining supersonic flow properties about bodies of revolution
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.
Development of near-field method for determining supersonic flow properties about bodies of revolution
The present study evaluates a higher-order modal method proposed by Leung for transient structural analysis entitled the force-derivative method. This method repeatedly integrates by parts with respect to time the convolution-integral form of the structural response to produce successively better approximations to the contribution of the higher modes which are neglected in the modal summation. Comparisons are made of the force-derivative, the mode-displacement, and the mode-acceleration methods for several numerical example problems for various times, levels of damping, and forcing functions. The example problems include a tip-loaded cantilevered beam and a simply-supported multispan beam. The force-derivative method is shown to converge to an accurate solution in fewer modes than either the mode-displacement or the mode-acceleration methods. In addition, for problems in which there are a large number of closely-spaced frequencies whose mode shapes have a negligible contribution to the response, the force derivative method is very effective in representing the effect of the important, but otherwise neglected, higher modes.
It is pointed out that edge delamination is a serious failure mechanism for laminated composite materials. Various numerical methods have been utilized in attempts to calculate the interlaminar stress components which precede delamination in a laminate. There are, however, discrepancies regarding the results provided by different methods, taking into account a finite-difference procedure, a perturbation procedure, and finite element approaches. The present investigation has the objective to assess the capacity of a finite difference method to predict the character and magnitude of the interlaminar stress distributions near an interface corner. A second purpose of the investigation is to determine if predictions by finite element method in-plane, interlaminar stress components near an interface corner represent actual laminate behavior.
A method is presented for modifying finite difference solutions of the potential equation to include the calculation of non-planar vortex wake features. The approach is an adaptation of Baker's 'cloud in cell' algorithm developed for the stream function-vorticity equations. The vortex wake is tracked in a Lagrangian frame of reference as a group of discrete vortex filaments. These are distributed to the Eulerian mesh system on which the velocity is calculated by a finite difference solution of the potential equation. An artificial viscosity introduced by the finite difference equations removes the singular nature of the vortex filaments. Computed examples are given for the two-dimensional time dependent roll-up of vortex wakes generated by wings with different spanwise loading distributions.
We report on the first part of a study of electron-hydrogen scattering, using a method which allows for the ab initio calculation of total and elastic cross sections at higher energies. In its general form the method uses complex 'radial' correlation functions, in a (Kohn) T-matrix formalism. The titled method, abbreviated Complex Correlation Kohn T (CCKT) method, is reviewed, in the context of electron-hydrogen scattering, including the derivation of the equation for the (complex) scattering function, and the extraction of the scattering information from the latter. The calculation reported here is restricted to S-waves in the elastic region, where the correlation functions can be taken, without loss of generality, to be real. Phase shifts are calculated using Hylleraas-type correlation functions with up to 95 terms. Results are rigorous lower bounds; they are in general agreement with those of Schwartz, but they are more accurate and outside his error bounds at a couple of energies,
The methods of NACA reports 865 and 1090 have been applied to the calculation of the rolling-moment and yawing-moment coefficients due to rolling for unswept wings with or without flaps or ailerons. The methods are based on lifting-line theory and allow the use of nonlinear section lift data. The method presented in this report permits calculations to be made somewhat beyond maximum lift for wings having no twist or continuous twist and employing airfoil sections which do not display large discontinuities in the lift curves. Calculations can be made up to maximum lift for wings with discontinuous twist such as that produced by partial-span flaps or ailerons, or both. Two calculated examples are presented in simplified computing forms in order to illustrate the procedures involve.
A nonsteady vortex-lattice method is introduced for predicting the dynamic stability derivatives of a delta wing undergoing an oscillatory motion. The analysis is applied to several types of small oscillations in pitch. The angle of attack varied between + or - 1 deg, with the mean held at 0 deg when the flow was assumed to be attached and between + or - 1 deg and the mean held at 15 deg when both leading-edge separation and wake roll-up were included. The computed results for damping in pitch are compared with several other methods and with experiments, and are found to be consistent and in good agreement.
A method is described for the analysis of the unsteady, incompressible potential flow associated with a helicopter rotor and it's wake in forward flight. This method is particularly useful in low advance ratio flight due to the major contribution, in the near field, of the deformed wake. The rotor geometry is prescribed and the unsteady wake geometry is computed from the local flow perturbation velocities. The wake is modeled as a full vortex lattice. The rotor geometry is arbitrary and several rotor blades can be represented. The unsteady airloads on the rotor blades are computed in the presence of the deformed rotor wake by a time-stepping technique. Solution for the load distribution on the blade surfaces is found by prescribing boundary conditions in a reference system which rotates with the blade tips. Transformation tensors are used to describe the contribution of the wake in the inertial system to the rotor in the rotating reference system. The effects of blade cyclic pitch variation are computed using a rotation tensor. The deformation of the wake is computed in the inertial frame. The wake is started impulsively from rest, allowing a natural convection of the wake with time.
The use of the so-called 'nose method' for computing critical cooling rates for glass formation is examined and compared with other methods, presenting data for the glass-forming systems SiO2, GeO2, and P2O5. It is shown that, for homogeneous crystallization, the nose-method will give an overestimate of Rc, a conclusion which was drawn after assessing the enfluence of a range of values for the parameters which control crystal growth and nucleation. The paper also proposes an alternative simple procedure (termed the 'cutoff method') for computing critical cooling rates from T-T-T diagrams, which was shown in the SiO2 and GeO2 systems to be superior to the nose method.
This paper describes a method which has been developed for estimating the safe fatigue life of compact, highly-stressed and inaccessible components for aeroplanes and helicopters of the Royal Air Force. It is explained why the Design Requirements for British Military Aircraft do not favor the use of a damage-tolerance approach in these circumstances.
A new theory of microwave absorption by a lossy dielectric sphere in a single mode rectangular cavity has been recently developed. The absorption was treated in the framework of an electromagnetic scattering problem. That theory is summarized here and calculated results that bear on optimizing the processing of materials are illustrated. Methods for observing power absorption and other results predicted by the scattering model are discussed. Cavity perturbation theory provides a bridge between theoretical calculations and experimental observations, and a special problem that arises when an established version of cavity perturbation theory is applied to spheres is identified, analyzed, and resolved. The direct problem of predicting shifts in frequency and Q from model calculations is discussed for a sphere in a cavity when the sphere's complex dielectric constant is known. Also, the inverse problem of determining the complex dielectric constants from measured values of those shifts is considered. The small sphere limit, where an electrostatic or quasistatic model is valid, is treated in detail, and planned work on parallel problems for larger spheres is described.
The methods of NACA Reports 865 and 1090 have been applied to the calculation of the rolling- and yawing-moment coefficients due to rolling for unswept wings with or without flaps or ailerons. The methods allow the use of nonlinear section lift data together with lifting-line theory. Two calculated examples are presented in simplified computing forms in order to illustrate the procedures involved.
Method which can be applied to the problem of an incompressible flow past an arbitrarily shaped cylindrical body in the presence of a magnetic field normal to the axis of the cylinder
An approximate solution technique was developed for 3-D inviscid, hypersonic flows. The method employs Maslen's explicit pressure equation in addition to the assumption of approximate stream surfaces in the shock layer. This approximation represents a simplification to Maslen's asymmetric method. The present method presents a tractable procedure for computing the inviscid flow over 3-D surfaces at angle of attack. The solution procedure involves iteratively changing the shock shape in the subsonic-transonic region until the correct body shape is obtained. Beyond this region, the shock surface is determined using a marching procedure. Results are presented for a spherically blunted cone, paraboloid, and elliptic cone at angle of attack. The calculated surface pressures are compared with experimental data and finite difference solutions of the Euler equations. Shock shapes and profiles of pressure are also examined. Comparisons indicate the method adequately predicts shock layer properties on blunt bodies in hypersonic flow. The speed of the calculations makes the procedure attractive for engineering design applications.
This report has been prepared to provide a practical method for determining the chordwise distribution of the rate of heat transfer from the surface of a wing or body of revolution to air. The method is limited in use to the determination of heat transfer from the forward section of such bodies when the flow is laminar. A comparison of the calculated average heat-transfer coefficient for the nose section of the wing of a Lockheed 12-A airplane with that experimentally determined shows a satisfactory agreement. A sample calculation is appended.
Many overarching standards, regulations, or other requirements necessitate that the critical temperature (T crit ) of an energetic material be known or estimate prior to operations such as heating. However, they rarely, if ever, provide a specific method for doing so. While other methods exist to calculate T crit , such as the Frank-Kamenetskii equation, computer simulation, etc., these cannot be readily utilized for energetics that are physical mixtures (i.e. most pyrotechnics) or materials where detailed material properties required for such calculations are lacking. Here, in this study, a COTS (commercial off-the-shelf) SBAT (Simulated Bulk Auto-ignition Test) apparatus is modified to perform Henkin cookoff tests. This creates a simple, efficient, and cost-effective solution to estimate required T crit values for energetic materials. After a historical overview of the evolution of the Henkin test, several common energetic materials (PETN, RDX, etc.) for which T crit has been readily calculated, and historical data is available, were analyzed for comparison and verification purposes. This was followed by a variety of pyrotechnic mixtures and other materials to where calculation methods cannot be readily used. The modified apparatus, as well as an updated sample shell sealing method, produced results that aligned well with historical data. It also produced reasonable T crit estimates for those materials where the Frank-Kamenetskii equation cannot be applied.
An approximate method for the calculation of the compressible laminar boundary layer with heat transfer and arbitrary pressure gradient, based on Thwaites' correlation concept, is presented. With the definition of dimensionless shear and heat-transfer parameters and an assumed correlation of these parameters in terms of a momentum parameter, a complete system of relations for calculating skin friction and heat transfer results. Knowledge of velocity or temperature profiles is not necessary in using this calculation method. When the method is applied to a convergent-divergent, axially symmetric rocket nozzle, it shows that high rates of heat transfer are obtained at the initial stagnation point and at the throat of the nozzle. Also indicated are negative displacement thicknesses in the convergent portion of the nozzle; these occur because of the high density within the lower portions of the cooled boundary layer. (author)
Woodward's panel method for subsonic and supersonic flow was improved by employing control points determined by exactly matching two-dimensional pressure at a finite number of points. The results show great improvement in the predicted pressure distribution of a flapped airfoil. With the paneling scheme of cosine law in both chordwise and spanwise directions, the method is shown to accurately predict leading edge and side edge suction forces of various configurations in subsonic and supersonic flow.