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.
The status of fourth-order conformal Weyl gravity theory as an alternative theory of gravity to standard second-order Einstein theory is briefly reviewed. The questions remaining open are addressed.
Solutions have been computed and results are presented for Problem 1, the case of Mach 9 transitional flow past a 7 deg half-angle cone at zero incidence. The solutions were computed using a code developed for the integration of the parabolized Navier-Stokes equations. The algorithm employed in the code is based on a Roe-type flux-difference-splitting scheme applied following a finite-volume approach. The basic algorithm has been modified to make it implicit and second-order accurate in the crossflow directions. Results are presented in terms of surface pressure and heat transfer as well as boundary layer profiles of pitot pressure, Mach number, and tangential velocity. The case was recalculated several times in an effort to determine sensitivities to such parameters as grid density, wall temperature, turbulence model parameters, as well as freestream expansion. Comparisons with the experimental data are presented and discussed.
Program LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) is an upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for 3D viscous hypersonic flows in chemical and thermal nonequilibrium. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with a modified Roe's averaging and with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. The code has been applied to Problem 8.2 of this workshop for the case of thermochemical nonequilibrium flow through a nozzle. Chemical reaction rates are defined with the model of Park (1987). Thermal nonequilibrium is modeled using a two-temperature approximation in which the vibrational energies of all molecules are assumed to be in equilibrium at a single temperature which is generally different from the translational-rotational temperature. Two grids were used to define the flow for the original problem, with a stagnation temperature of 6500 K. A third case with a stagnation temperature of 10,000 K is also presented. The solution domain includes the converging nozzle, subsonic flow domain in which the gas is substantially in thermochemical equilibrium and the diverging nozzle, hypersonic flow domain in which the gas is substantially in thermochemical nonequilibrium.
A dissipative compact two-four scheme (second-order time, fourth-order space) similar to the original MacCormack scheme has been developed, that exhibits greater accuracy than conventional fourth-order schemes. The dissipative nature of the scheme allows it to resolve weak discontinuities without artificial damping. A derivation of the scheme is presented, as well as the theoretical stability characteristics. The temporal scheme is then generalized into a steady-state formulation which achieves fourth-order spatial accuracy at steady-state. Several test problems are used to show that the scheme is more accurate than the traditional MacCormack scheme, and is nearly as efficient.
An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for three-dimensional, viscous, hypersonic flows in chemical and thermal nonequilibrium is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. The relaxation strategy is well suited for computers employing either vector or parallel architectures, and the relation between computer architecture and algorithm is emphasized. It is also well suited to the numerical solution of the governing equations on unstructured grids. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large. block tri-diagonal systems. A single relaxation step depends only on information from nearest neighbors. Predictions for pressure distributions, surface heating, and aerodynamic coefficients compare well with experimental data for Mach 10 flow over a blunt body. Predictions for the hypersonic flow of air in chemical and thermal nonequilibrium (velocity = 8917 m/s, altitude = 78 km.) over the Aeroassist Flight Experiment (AFE) configuration obtained on a multi-domain grid are discussed.
A laminate strength and life prediction methodology has been postulated for laminated composites which accounts for the progressive development of microstructural damage to structural failure. A damage dependent constitutive model predicts the stress redistribution in an average sense that accompanies damage development in laminates. Each mode of microstructural damage is represented by a second-order tensor valued internal state variable which is a strain like quantity. The mechanics framework together with the global-local strategy for predicting laminate strength and life is presented in the paper. The kinematic effects of damage are represented by effective engineering moduli in the global analysis and the results of the global analysis provide the boundary conditions for the local ply level stress analysis. Damage evolution laws are based on experimental results.
The objective of this work is to develop a space-time accurate numerical method for the solution of incompressible Navier-Stokes equations in generalized coordinates. The resulting code is to be used for direct and large-eddy simulation of turbulence in complex geometries. In a previous paper, the system of Navier-Stokes equations in general curvilinear coordinates was solved by a second-order accurate finite-difference scheme. Satisfactory results were obtained for several flows in two and three dimensions. The system of Navier-Stokes for the fluxes are given in Orlandi (1989). The main deficiency of the numerical scheme was the large CPU time required for the solution of the Poisson equation for the 'pressure' field. The point SOR relaxation, in conjunction with a multigrid scheme, was used for the Poisson equation. In some cases, particularly with very fine grids, it was impossible to obtain a divergent-free flow. A preliminary attempt is made to compute the spatially evolving flow of Swearingen & Blackwelder. To reduce the streamwise distance, the inflow was at a distance x = 60 cm from the leading edge.
Afterbody pressure distributions and wave drag were calculated using a second-order theory for a variety of conical boattails at zero angle of attack. Results are presented for Mach numbers from 1.5 to 4.5, area ratios from 0.200 to 0.800, and boattail angle from 3 degrees to 11 degrees. The results indicate that for a given boattail angle, the wave drag decreases with increasing Mach number and area ratio. The wave drag, for a constant area ratio, increases with increasing boattail angle. For a specific Mach number, area ratio, and fineness ratio, a comparison of the wave-drag coefficients for conical, tangent-parabolic, and secant-parabolic boattails showed the conical boattail to have the smallest wave drag.
Pressure distributions and forces for a series of four bodies of revolution having nose-fineness ratios varying from 4 to 10 have been obtained and compared with theory for a Mach number of 3.12, a Reynolds number range of 2x10(sup)6 to 14x10(sup)6, and angles of attack from zero to 9 degrees. In general, a comparison of the experimental data with a second-order theory showed good agreement for the range of variables investigated.
Pressure distributions for a series of four boattailed bodies of revolution were obtained and compared with theory for a Mach number of 3.12, a Reynolds number range of 2 x 10 to 6th power to 14 x 10 to the 6th power, and angles of attack from zero to 9 degrees. Second-order theory adequately predicted the pressure distribution for regions free of the effects of cross-flow separation.
Shock position was assumed to be controlled by air bleed through a diffuser-exit bypass door operated by a second-order servomotor. Pressure disturbance resulting from changes in bypass-door position had whereas a disturbance in engine speed had a large effect on diffuser pressure and could impose severe requirements on diffuser control response. The controlled diffuser response improved with faster servomotor response and smaller diffuser dead time and lag. The effect of diffuser dead time on diffuser response was greater than the effect of diffuser lag.
A discontinuous variation of coefficients of the differential equation describing the linear control system before nonlinear elements are added is studied in detail. The nonlinear feedback is applied to a second-order system. Simulation techniques are used to study performance of the nonlinear control system and to compare it with the linear system for a wide variety of inputs. A detailed quantitative study of the influence of relay delays and of a transport delay is presented.
A new mathematical technique, due to Milne-Thomson, is used to obtain an improved form of the method of Poggi for calculating the effect of compressibility on the subsonic flow past an obstacle. By means of this new method, the difficult surface integrals of the original Poggi method can be replaced by line integrals. These line integrals are then solved by the use of residue theory. In this way an equation is obtained giving the second-order effect of compressibility on the velocity of the fluid. The method is practicable for obtaining the higher-order effects of compressibility on the velocity field. As an illustration of the general result, the flow past an elliptic cylinder is discussed.
Methods are presented that use general correlative time-response input and output data for a linear system to determine the frequency-response function of that system. These methods give an exact description of any linear system for which such transient data are available. Examples are shown of application of a method to both an underdamped and a critically damped exact second-order system, and to an exact first-order system with and without dead time. Experimental data for a turbine-propeller engine showing the response of engine speed to change in propeller-blade angle are presented and analyzed.
The particular integrals of the second-order and third-order Prandtl-Busemann iteration equations for the flow of a compressible fluid are obtained by means of the method in which the complex conjugate variables are utilized as the independent variables of the analysis. The assumption is made that the Prandtl-Glauert solution of the linearized or first-order iteration equation for the two-dimensional flow of a compressible fluid is known. The forms of the particular integrals, derived for subsonic flow, are readily adapted to supersonic flows with only a change in sign of one of the parameters of the problem.
Several methods of obtaining the time response of Linear systems to either a unit impulse or an arbitrary input from frequency-response data are described and compared. Comparisons indicate that all the methods give good accuracy when applied to a second-order system; the main difference is the required computing time. The methods generally classified as inverse Laplace transform methods were found to be most effective in determining the response to a unit impulse from frequency-response data of higher order systems. Some discussion and examples are given of the use of such methods as flight-data-analysis techniques in predicting loads and motions of a flexible aircraft on the basis of simple calculations when the aircraft frequency response is known.
In controller design for flexible structures, certain system modes are extremely important for the overall performance of the structure. A reduced-order model (ROM) based control focuses on these modes, providing a viable, active control algorithm for large systems. Unfortunately, unmodeled structure dynamics can interact with the ROM controller (CSI) and cause crippling deterioration of system performance, possibly to the point that system stability is lost. A residual model filter (RMF) eliminates one channel of control structure interaction (CSI), while adding only a simple, second-order filter to the control loop. Thus, the ROM controller can be designed independently, based strictly on performance criteria, and residual mode filters can then be selected to compensate for CSI. A flexible robot manipulator is used for preliminary experimentation with the ROM/RMF design methodology. Since the controller was to be implemented both with, and without compensation for CSI, the ROM control gains are carefully chosen such that closed loop stability is never compromised. In this way, RMF effectiveness is easily evaluated in terms of the improvement in system performance resulting from CSI compensation.
A solution method, fourth-order accurate in the body-normal direction and second-order accurate in the stream surface directions, to solve the compressible 3-D boundary layer equations is presented. The transformation used, the discretization details, and the solution procedure are described. Ten validation cases of varying complexity are presented and results of calculation given. The results range from subsonic flow to supersonic flow and involve 2-D or 3-D geometries. Applications to laminar flow past wing and fuselage-type bodies are discussed. An interface procedure is used to solve the surface Euler equations with the inviscid flow pressure field as the input to assure accurate boundary conditions at the boundary layer edge. Complete details of the computer program used and information necessary to run each of the test cases are given in the Appendix.