Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite difference”

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.

At least 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

The application of Riegels' rule and time-like damping to transonic flow calculations

Finite difference relaxation solutions of the nonlinear small perturbation equations have proven reliable and successful in determining the transonic flowfields about thin airfoils. However, application of the small perturbation approach to thick airfoils usually results in an accuracy less than desirable. The incorporation of Riegels' Rule and time-like damping into the small perturbation approach and their application to thick and thin airfoils in transonic flow are discussed. Studies for thick and thin airfoils are presented. It is concluded that Riegels' Rule and damping should both be included in small perturbation transonic flow calculations.

Herron, D. L.↗
NASA NTRS1973

Survey of computational methods for three-dimensional supersonic inviscid flows with shocks

Hyperbolic systems of partial differential equations governing supersonic inviscid flows are discussed and analyzed. Finite-difference analogues for integrating these systems in the interior of fluid domains are described from two points of view: a differential form approach and an integral form approach. The algorithms presented are analyzed for stability and accuracy. The concept of time splitting is discussed and supplied to these methods to achieve increased numerical efficiency. Techniques for treating conditions at the boundaries of the fluid domain and shock-wave discontinuities at surfaces within the domain are described.

Maccormack, R. W.↗
NASA NTRS1974

Computer program for calculating laminar, transitional, and turbulent boundary layers for a compressible axisymmetric flow

A finite-difference program is described for calculating the viscous compressible boundary layer flow over either planar or axisymmetric surfaces. The flow may be initially laminar and progress through a transitional zone to fully turbulent flow, or it may remain laminar, depending on the imposed boundary conditions, laws of viscosity, and numerical solution of the momentum and energy equations. The flow may also be forced into a turbulent flow at a chosen spot by the data input. The input may contain the factors of arbitrary Reynolds number, free-stream Mach number, free-stream turbulence, wall heating or cooling, longitudinal wall curvature, wall suction or blowing, and wall roughness. The solution may start from an initial Falkner-Skan similarity profile, an approximate equilibrium turbulent profile, or an initial arbitrary input profile.

Albers, J. A.↗
NASA NTRS1974

Numerical computation of two dimensional viscous blunt body flows with an impinging shock, part 2

Two-dimensional viscous blunt body flows with an impinging shock have been computed using a time-dependent finite-difference method which solves the complete set of Navier-Stokes equations for a compressible flow. For low Reynolds number flows, the entire flow field, including the bow shock and impinging shock, has been captured in the computation. For higher Reynolds number flows, the bow shock is treated as a discontinuity across which the Rankine-Hugoniot equations are applied, while the boundary layer and interaction regions are captured as before. Using this latter shock-fitting approach, a Type III shock interaction flow field has been computed with flow conditions corresponding to the space shuttle orbiter freestream conditions at 61 km (200,000 ft).

Holst, T. L.↗
NASA NTRS1974

Numerical solution of the steady-state Navier-Stokes equations for hypersonic flow about blunt axisymmetric bodies

The steady-state Navier-Stokes equations are solved for hypersonic flow about blunt axisymmetric bodies. The equations of motion are solved by successive approximations using an implicit finite-difference scheme. The results are compared with viscous shock-layer theory, experimental data, and time-dependent solutions of the Navier-Stokes equations. It is demonstrated that viscous shock-layer theory is sufficiently accurate for the range of flight conditions normally encountered by entry vehicles.

Anderson, E. C.↗
NASA NTRS1974

Separated flow over bodies of revolution using an unsteady discrete-vorticity cross wake. Part 2: Computer program description

A method is developed to determine the flow field of a body of revolution in separated flow. The computer was used to integrate various solutions and solution properties of the sub-flow fields which made up the entire flow field without resorting to a finite difference solution to the complete Navier-Stokes equations. The technique entails the use of the unsteady cross flow analogy and a new solution to the two-dimensional unsteady separated flow problem based upon an unsteady, discrete-vorticity wake. Data for the forces and moments on aerodynamic bodies at low speeds and high angle of attack (outside the range of linear inviscid theories) such that the flow is substantially separated are produced which compare well with experimental data. In addition, three dimensional steady separated regions and wake vortex patterns are determined. The computer program developed to perform the numerical calculations is described.

Marshall, F. J.↗
NASA NTRS1974

The numerical calculation of laminar boundary-layer separation

Iterative finite-difference techniques are developed for integrating the boundary-layer equations, without approximation, through a region of reversed flow. The numerical procedures are used to calculate incompressible laminar separated flows and to investigate the conditions for regular behavior at the point of separation. Regular flows are shown to be characterized by an integrable saddle-type singularity that makes it difficult to obtain numerical solutions which pass continuously into the separated region. The singularity is removed and continuous solutions ensured by specifying the wall shear distribution and computing the pressure gradient as part of the solution. Calculated results are presented for several separated flows and the accuracy of the method is verified. A computer program listing and complete solution case are included.

Klineberg, J. M.↗
NASA NTRS1974

Separated flow over bodies of revolution using an unsteady discrete-vorticity cross wake. Part 1: Theory and application

A method is developed to determine the flow field of a body of revolution in separated flow. The technique employed is the use of the computer to integrate various solutions and solution properties of the sub-flow fields which made up the entire flow field without resorting to a finite difference solution to the complete Navier-Stokes equations. The technique entails the use of the unsteady cross flow analogy and a new solution to the required two-dimensional unsteady separated flow problem based upon an unsteady, discrete-vorticity wake. Data for the forces and moments on aerodynamic bodies at low speeds and high angle of attack (outside the range of linear inviscid theories) such that the flow is substantially separated are produced which compare well with experimental data. In addition, three dimensional steady separation regions and wake vortex patterns are determined.

Marshall, F. J.↗
NASA NTRS1974

Analysis of sound propagation in ducts using the wave envelope concept

A finite difference formulation is presented for sound propagation in a rectangular two-dimensional duct without steady flow for plane wave input. Before the difference equations are formulated, the governing Helmholtz equation is first transformed to a form whose solution does not oscillate along the length of the duct. This transformation reduces the required number of grid points by an order of magnitude, and the number of grid points becomes independent of the sound frequency. Physically, the transformed pressure represents the amplitude of the conventional sound wave. Example solutions are presented for sound propagation in a one-dimensional straight hard-wall duct and in a two-dimensional straight soft-wall duct without steady flow. The numerical solutions show evidence of the existence along the duct wall of a developing acoustic pressure diffusion boundary layer which is similar in nature to the conventional viscous flow boundary layer. In order to better illustrate this concept, the wave equation and boundary conditions are written such that the frequency no longer appears explicitly in them. The frequency effects in duct propagation can be visualized solely as an expansion and stretching of the suppressor duct.

Baumeister, K. J.↗
NASA NTRS1974

Transonic conical flow

The problem of inviscid, steady transonic conical flow, formulated in terms of the small disturbance theory, is studied. The small disturbance equation and similarity rules are presented, and a boundary value problem is formulated for the case of a supersonic freestream Mach number. The equation for the perturbation potential is solved numerically using an elliptic finite difference system. The difference equations are solved with a point relaxation algorithm that is also capable of capturing the shock wave during the iteration procedure by using the boundary conditions at the shock. Numerical calculations, for shock location, pressure distribution and drag coefficient, are presented for a family of nonlifting conical wings. The theory of slender wings is also presented and analytical results for pressure and drag coefficients are obtained.

Agopian, K. G.↗
NASA NTRS1974

Comparison of several methods for predicting separation in a compressible turbulent boundary layer

Several methods for predicting the separation point for a compressible turbulent boundary layer were applied to the flow over a bump on a wind-tunnel wall. Measured pressure distributions were used as input. Two integral boundary-layer methods, three finite-difference boundary-layer methods, and three simple methods were applied at five free-stream Mach numbers ranging from 0.354 to 0.7325. Each of the boundary-layer methods failed to explicitly predict separation. However, by relaxing the theoretical separation criteria, several boundary-layer methods were made to yield reasonable separation predictions, but none of the methods accurately predicted the important boundary-layer parameters at separation. Only one of the simple methods consistently predicted separation with reasonable accuracy in a manner consistent with the theory. The other methods either indicated several possible separation locations or only sometimes predicted separation.

Gerhart, P. M.↗
NASA NTRS1974

Flutter analysis of swept-wing subsonic aircraft with parameter studies of composite wings

A computer program is presented for the flutter analysis, including the effects of rigid-body roll, pitch, and plunge of swept-wing subsonic aircraft with a flexible fuselage and engines mounted on flexible pylons. The program utilizes a direct flutter solution in which the flutter determinant is derived by using finite differences, and the root locus branches of the determinant are searched for the lowest flutter speed. In addition, a preprocessing subroutine is included which evaluates the variable bending and twisting stiffness properties of the wing by using a laminated, balanced ply, filamentary composite plate theory. The program has been substantiated by comparisons with existing flutter solutions. The program has been applied to parameter studies which examine the effect of filament orientation upon the flutter behavior of wings belonging to the following three classes: wings having different angles of sweep, wings having different mass ratios, and wings having variable skin thicknesses. These studies demonstrated that the program can perform a complete parameter study in one computer run. The program is designed to detect abrupt changes in the lowest flutter speed and mode shape as the parameters are varied.

Housner, J. M.↗
NASA NTRS1974

Boundary-layer transition on a plate subjected to simultaneous spanwise and chordwise pressure gradients

The boundary-layer transition on a short plate was studied by means of the china-clay visual technique. The plate model was mounted in a wind tunnel so that it was subjected to small simultaneous spanwise and chordwise pressure gradients. Results of the experimental study, which was performed at three subsonic velocities, indicated that the transition pattern was appreciably curved in the spanwise direction but quite smooth and well behaved. Reasonable comparisons between predictions of transition and experiment were obtained from two finite-difference two-dimensional boundary-layer calculation methods which incorporated transition models based on the concept of a transition intermittency factor.

Boldman, D. R.↗
NASA NTRS1974

Analysis of dynamic stall using unsteady boundary-layer theory

The unsteady turbulent boundary layer and potential flow about a pitching airfoil are analyzed using numerical methods to determine the effect of pitch rate on the delay in forward movement of the rear flow reversal point. An explicit finite difference scheme is used to integrate the unsteady boundary layer equations, which are coupled at each instant of time to a fully unsteady and nonlinear potential flow analysis. A substantial delay in forward movement of the reversal point is demonstrated with increasing pitch rate, and it is shown that the delay results partly from the alleviation of the gradients in the potential flow, and partly from the effects of unsteadiness in the boundary layer itself. The predicted delay in flow-reversal onset, and its variation with pitch rate, are shown to be in reasonable agreement with experimental data relating to the delay in dynamic stall. From the comparisons it can be concluded (a) that the effects of time-dependence are sufficient to explain the failure of the boundary layer to separate during the dynamic overshoot, and (b) that there may be some link between forward movement of the reversal point and dynamic stall.

Scruggs, R. M.↗
NASA NTRS1973

Stability analysis of cylinders with circular cutouts

The stability of axially compressed cylinders with circular cutouts is analyzed numerically. An extension of the finite-difference method is used which removes the requirement that displacement components be defined in the directions of the grid lines. The results of this nonlinear analysis are found to be in good agreement with earlier experimental results.

Almroth, B. O.↗
NASA NTRS1973

Initiation and structure of axisymmetric eddies in a rotating stream

Axisymmetric flow of a rotating stream is examined numerically to determine conditions under which an isolated eddy will form on the axis of rotation. An explicit finite difference procedure is used to integrate the time dependent transport equations. Solutions provide details of the flow structure and are presented for a range of Reynolds numbers and swirl ratios. Calculated results are interpreted in terms of recent physical experiments insofar as is possible.

Kopecky, R. M.↗
NASA NTRS1973

Heating of foods in space-vehicle environments

In extended space missions, foods will be heated to enhance the psychological as well as the physiological well-being of the crew. In the low-gravity space environment natural convection is essentially absent so that the heat transfer within the food is by conduction alone. To prevent boiling in reduced pressure environments the maximum temperature of the heating system is severely limited. The Skylab food-heating system utilizes a tray with receptables for the food containers. The walls of the receptacles are lined with thermally controlled, electrical-resistance, blanket-type heating elements. A finite difference model is employed to perform parametric studies on the food-heating system. The effects on heating time of the (1) thermophysical properties of the food, (2) heater power level, (3) initial food temperatures, (4) container geometry, and (5) heater control temperature are presented graphically. The optimal heater power level and container geometry are determined.

Bannerot, R. B.↗