Search NASA⌕ Search

SEARCH · Search NASA

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

Solutions of the Navier-Stokes equations for vortex breakdown

Steady solutions of the Navier-Stokes equations, in terms of velocity and pressure, for breakdown in an unconfined viscous vortex are obtained numerically using the artificial compressibility technique of Chorin combined with an ADI finite-difference scheme. Axisymmetry is assumed and boundary conditions are carefully applied at the boundaries of a large finite region in an axial plane while resolution near the axis is maintained by a coordinate transformation. The solutions, which are obtained for Reynolds numbers up to 200 based on the free-stream axial velocity and a characteristic core radius, show that breakdown results from the diffusion and convection of vorticity away from the vortex core which, because of the strong coupling between the circumferential and axial velocity fields in strongly swirling flows, can lead to stagnation and reversal of the axial flow near the axis.

Grabowski, W. J.↗
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

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 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

An improved numerical process for solution of solid mechanics problems

This paper gives an overview of the development and status of an improved numerical process for the solution of solid mechanics problems. The proposed process uses a mixed formulation with the fundamental unknowns consisting of both stress and displacement parameters. The problem is formulated either by means of first-order partial differential equations or in a variational form by using a Hellinger-Reissner-type mixed variational principle. For presentation purposes, the components of a numerical process are characterized and the criteria for an ideal process are outlined. Commonly used finite-difference and finite-element procedures are examined in the light of these criteria and it is shown that they fall short in a number of ways. The proposed numerical process, on the other hand, satisfies most of the optimality criteria and appears to be particularly suited for use with the forthcoming generation computers.

Noor, A. K.↗
NASA NTRS1973

Calculation of separated flows at subsonic and transonic speeds

A boundary-layer integral approach is combined with a finite-difference relaxation method to calculate viscous interactions between separated flows at subsonic and transonic velocities. Results are obtained for separated laminar flows on circular-arc airfoils at zero angle of attack and are compared with data of Collins (1972). Inviscid and viscous flows are covered.

Klineberg, J. M.↗
NASA NTRS1974

Thermoacoustic convection of fluids in low gravity

The heat flow in a confined perfect gas in low gravity is investigated, including the effects of conduction and thermal convection. Buoyancy-driven flow is neglected, due to the low-gravity environment, but the effect of thermoacoustic motion due to fluid compressibility is included. One-dimensional mathematical models are constructed from the conservation equations for a compressible, viscous, heat-conducting fluid. A conservative, time-dependent finite-difference method is used to generate numerical solutions on a digital computer. Problems for flat plates and cylindrical segments are solved for specified thermal boundary conditions. Numerical results are given which indicate that thermoacoustic convection can significantly increase the transient heat flow over conduction model predictions for cases where a confined gas is rapidly heated.

Spradley, L. W.↗
NASA NTRS1974

Hypersonic chemically reacting viscous shock layers over sphere-cones and cylinder-wedges

Hypersonic, nonequilibrium viscous flow over nonanalytic blunt bodies is considered. The equations which govern the viscous shock-layer flow are presented and the method by which the equations are solved is discussed. The predictions of the present finite-difference method are compared with other numerical predictions as well as with experimental data. Three flow conditions are considered; the experimental, wind tunnel conditions of Pappas and Lee for a 7.5 deg sphere-cone at Mach 13 and two cases considered by Kang and Dunn, a 9 deg sphere-cone at 233,000 ft and a 20 deg sphere-cone at 280 and 310 Kft. The predictions of the present method agreed well with the experimental heat-transfer data, but substantial differences were found between the present predictions and the more approximate predictions of Kang and Dunn for heat-transfer distributions and temperature profiles.

Miner, E. W.↗
NASA NTRS1973

Computer simulation of shock wave boundary layer interaction

Development of a boundary layer correction for an analysis program previously used by Bauer et al. (1972) to study two-dimensional flow past transonic airfoils at off-design conditions where shocks which interact with the boundary layer appear. Using a method which combines conformal mapping with a finite-difference scheme due to Murman and Cole (1971), an analysis is made of the interaction between shock waves on a supercritical wing section and a turbulent boundary layer, assuming that the shocks are weak and that separation is insignificant. It is shown that relatively simple shock wave/boundary layer interactions of the kind envisioned in this study can be treated by determining the displacement thickness iteratively together with the flow and altering the profile accordingly. It is concluded that the proposed boundary layer correction appears to furnish an essential improvement in cases where both pronounced aft loading and a large enclosed supersonic region of flow are present.

Bauer, F.↗
NASA NTRS1974

Fast direct numerical solution of the nonhomogeneous Cauchy-Riemann equations

A fast direct (noniterative) 'Cauchy-Riemann Solver' is developed for solving the finite-difference equations representing systems of first-order elliptic partial differential equations in the form of the nonhomogeneous Cauchy-Riemann equations. The method is second-order accurate and requires approximately the same computer time as a fast cyclic-reduction Poisson solver. The accuracy and efficiency of the direct solver are demonstrated in an application to solving an example problem in aerodynamics: subsonic inviscid flow over a biconvex airfoil. The analytical small-perturbation solution contains singularities, which are captured well by the computational technique. The algorithm is expected to be useful in nonlinear subsonic and transonic aerodynamics.

Lomax, H.↗
NASA NTRS1974

Inviscid drag at transonic speeds

A systematic asymptotic expansion procedure is introduced into the Euler equations to obtain the transonic, small-disturbance potential equation and corresponding relation for the drag. The relation consists of an integral around any contour enclosing the body and an integral along all shocks contained within the contour. Interpretations are given for the origin of the pressure drag due to shock waves, lift, and wind-tunnel boundaries. Numerical procedures are discussed for evaluating these quantities from finite-difference relaxation calculations. Computed results are included for various contours, and results are compared with some experimental data.

Murman, E. M.↗
NASA NTRS1974

Viscous shock-layer flows for the space shuttle windward plane of symmetry

Predictions are presented for hypersonic, nonequilibrium viscous flows over nonanalytic blunt bodies at reentry conditions. Emphasis is placed on predictions of the viscous flowfield for the windward plane of symmetry of the space shuttle and other axisymmetric bodies which approximate the space shuttle geometry. The equations which govern the viscous shock-layer flow are presented, and the method by which the equations are solved is discussed. The predictions of the present finite-difference method are compared with predictions of other numerical methods as well as with experimental data.

Miner, E. W.↗
NASA NTRS1974

Using the NASTRAN Thermal Analyzer to simulate a flight scientific instrument package

The NASTRAN Thermal Analyzer has proven to be a unique and useful tool for thermal analyses involving large and complex structures where small, thermally induced deformations are critical. Among its major advantages are direct grid point-to-grid point compatibility with large structural models; plots of the model that may be generated for both conduction and boundary elements; versatility of applying transient thermal loads especially to repeat orbital cycles; on-line printer plotting of temperatures and rate of temperature changes as a function of time; and direct matrix input to solve linear differential equations on-line. These features provide a flexibility far beyond that available in most finite-difference thermal analysis computer programs.

Lee, H.-P.↗
NASA NTRS1974

Development of solution techniques for nonlinear structural analysis

Nonlinear structural solution methods in the current research literature are classified according to order of the solution scheme, and it is shown that the analytical tools for these methods are uniformly derivable by perturbation techniques. A new perturbation formulation is developed for treating an arbitrary nonlinear material, in terms of a finite-difference generated stress-strain expansion. Nonlinear geometric effects are included in an explicit manner by appropriate definition of an applicable strain tensor. A new finite-element pilot computer program PANES (Program for Analysis of Nonlinear Equilibrium and Stability) is presented for treatment of problems involving material and geometric nonlinearities, as well as certain forms on nonconservative loading.

Vos, R. G.↗
NASA NTRS1972

Steady-state and dynamic performance of a gas-lubricated seal

Steady-state and dynamic performance of a gas-lubricated, self-acting face seal was determined using numerical methods based on a variable grid, finite-difference, time-transient procedure. Results were obtained for a gas turbine main shaft seal operating at 206.9 newton per square centimeter (300 psi) sealed air pressure and 152.4 meters per second (500 ft/sec) sliding velocity. Analysis of the seal dynamics revealed that the response of the seal nosepiece to runout of the seat face is markedly affected by secondary seal friction and by nosepiece inertia. The nosepiece response was determined for various levels of secondary seal friction and seat face runout magnitudes.

Colsher, R.↗
NASA NTRS1974

The influence of the computational mesh on accuracy for initial value problems with discontinuous or nonunique solutions

Discontinuous, or weak, solutions of the wave equation, the inviscid form of Burgers equation, and the time-dependent, two-dimensional Euler equations are studied. A numerical method of second-order accuracy in two forms, differential and integral, is used to calculate the weak solutions of these equations for several initial value problems, including supersonic flow past a wedge, a double symmetric wedge, and a sphere. The effect of the computational mesh on the accuracy of computed weak solutions including shock waves and expansion phenomena is studied. Modifications to the finite-difference method are presented which aid in obtaining desired solutions for initial value problems in which the solutions are nonunique.

Maccormack, R. W.↗
NASA NTRS1947

Determination of elastic stresses in gas-turbine disks

A method is presented for the calculation of elastic stresses in symmetrical disks typical of those of a high-temperature gas turbine. The method is essentially a finite-difference solution of the equilibrium and compatibility equations for elastic stresses in a symmetrical disk. Account can be taken of point-to-point variations in disk thickness, in temperature, in elastic modulus, in coefficient of thermal expansion, in material density, and in Poisson's ratio. No numerical integration or trial-and-error procedures are involved and the computations can be performed in rapid and routine fashion by nontechnical computers with little engineering supervision. Checks on problems for which exact mathematical solutions are known indicate that the method yields results of high accuracy. Illustrative examples are presented to show the manner of treating solid disks, disks with central holes, and disks constructed either of a single material or two or more welded materials. The effect of shrink fitting is taken into account by a very simple device.

Manson, S S↗