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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
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

Internal hypersonic flow

An approach for studying hypersonic internal flow with the aid of a thin-shock-layer approximation is discussed, giving attention to a comparison of thin-shock-layer results with the data obtained on the basis of the imposition theory or a finite-difference integration of the Euler equations. Relations in the case of strong interaction are considered together with questions of pressure distribution and aspects of the boundary-layer solution.

Lin, T. C.↗
NASA NTRS1974

Stress, stability and vibration of complex, branched shells of revolution

A comprehensive computer program, designated BOSOR4, for analysis of the stress, stability, and vibration of segmented, ring-stiffened, branched shells of revolution and prismatic shells and panels is described. The program performs large-deflection axisymmetric stress analysis, small-deflection nonsymmetric stress analysis, modal vibration analysis with axisymmetric nonlinear prestress included, and buckling analysis with axisymmetric or nonsymmetric prestress. One of the main advantages of the code is the provision for realistic engineering details such as eccentric load paths, internal supports, arbitrary branching conditions, and a 'library' of wall constructions. The program is based on the finite-difference energy method, which is very rapidly convergent with increasing numbers of mesh points. The organization of the program is briefly described with the flow of calculations charted for each of the types of analysis. Overlay charts and core storage requirements are given for the CDC 6600, IBM 370/165, and UNIVAC 1108 versions of BOSOR4.

Bushnell, D.↗
NASA NTRS1974

Deviation angle estimation for axial-flow compressors using inviscid flow solutions

Development of a more satisfactory method of estimating deviation angles was attempted using analytical procedures. Solutions for inviscid, irrotational flow in the blade-to-blade plane were obtained with a finite-difference calculation method. Deviation angles for three cascades with rounded trailing edges were estimated using the inviscid flow solutions and four existing trailing edge hypotheses. The estimated deviation angles were compared with experimental values for a range of incidence angle and inlet flow angle. None of the hypotheses were satisfactory at all incidence angles for all the cascades considered. An alternative hypothesis is presented which is suitable for one of the cascade profiles over a range of incidence and blade setting angle.

Miller, M. J.↗
NASA NTRS1974

Liquefaction dynamics of n-octadecane in cylindrical coordinates

The study presents one method of mathematically modeling the effects of gravity-induced free convection upon the melting phenomena of a finite paraffin slab. The theoretical model is based on the numerical computer solution of the solid phase and liquid phase energy equations, coupled with an imposed ideal-viscous flow approximation. An implicit alternating-direction technique is used as the finite-difference approximation for the basic energy equations in the numerical solution. The theoretical predictions are compared to experimental data from an annular test cell, in which n-octadecane, C18H38, was used as the test material. The results show that the model is a good first approximation to the solution of solid-liquid phase change involving gravity-induced free convection.

Bain, R. L.↗
NASA NTRS1974

Influence of ablation impurities on blunt body re-entry ionization

Electron density profiles which include the effect of an ablated sodium impurity were computed for the boundary layer on a blunt-nosed body re-entering the atmosphere at 7.62 km/sec. Profiles are computed from the nose to a distance of four diameters along the RAM C-payload. A finite-difference, laminar, nonequilibrium chemistry boundary-layer program was used. Comparison of theory with S-band diagnostic antenna results, electron concentration deduced from X- and C-band attenuation data, and Langmuir probe data at several different aft body locations show that agreement is good at high altitude. At the lower altitudes there is disagreement between theory and S-band antenna data where the apparent discrepancy is attributed to the three-body recombination rate constant used for deionization of sodium coupled with the effect of angle of attack.

Schexnayder, C. J., Jr.↗
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

On the shock-on-shock interaction problem

The unsteady, inviscid flowfield that results when a supersonic vehicle strikes a planar oblique shock wave, though difficult to simulate experimentally, is quite easy to model and compute numerically. The complicated flowfield, which contains multiple shock wave interactions, is determined using a second-order, shock-capturing, finite-difference approach which solves the time-dependent Euler equations under a self-similar transformation. A series of numerical results for a simple two-dimensional wedge is presented which describes the entire disturbed region, including the wave structure, and shows good agreement with the available two- and three-dimensional experimental data.-

Kutler, P.↗
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

Numerical simulation of high Reynolds number transonic flows

A code has been developed for simulating high Reynolds number transonic flow fields of arbitrary configuration. An explicit finite-difference method with time splitting is used to solve the time-dependent equations for compressible turbulent flow. A nonorthogonal computational mesh of arbitrary configuration facilitates the description of the flow field. The code is applied to simulate the flow over an 18 percent thick circular-arc biconvex airfoil at zero angle of attack and free-stream Mach number of 0.775. A simple mixing-length model is used to describe the turbulence and chord Reynolds numbers of 1, 2, 4, and 10 million are considered. The solution describes in sufficient detail both the shock-induced and trailing-edge separation regions, and provides the profile and friction drag.

Deiwert, G. S.↗
NASA NTRS1974

Investigation of the three-dimensional turbulent flow downstream of swept slot injection in hypersonic flow

Results of an experimental and numerical investigation of tangential swept slot injection (sweep angles of 22.5 and 45 deg) into a thick turbulent boundary layer at Mach 6 are presented. Film cooling effectiveness, skin friction, and flow structure downstream of the swept slot injection are investigated. The data are compared to that for unswept slots, and it is found that cooling effectiveness and skin-friction reductions are not significantly affected by sweeping the slot. Predictions of cooling effectiveness and skin friction obtained by a numerical finite-difference technique agree reasonably well with experimental surface variables. As in previous supersonic two-dimensional slot research, reduced mixing was found downstream of the slot lip in the present three-dimensional case.

Hefner, J. N.↗
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

Diffusive ignition and combustion in a wall jet

Hydrogen is injected from a downstream facing step in a wall into a high-temperature stream. Temperature and hydroxyl radical concentration are measured downstream of the injection plane by means of ultraviolet absorption spectroscopy. The experimental results are compared with theory which is based on a finite-difference solution of boundary-layer equations. Finite-rate kinetics equations are included in the analysis. The analytic predictions are also compared with previously obtained experimental results which are based on probe measurements. Comparison is made between calculated and observed ignition distances.

Kurkov, A. P.↗
NASA NTRS1974

A comparison of numerical solutions of the advective equation

Second- and third-order finite-difference methods recently applied to problems in high-speed fluid flow are applied to the model advection equation cast in conservative form. The differencing methods considered use forward time differencing with both centered and preferential space differences employed in the predictor-corrector sequences. The free parameter required for stability in the third-order method is adjusted to cause the solution to be either minimum dispersive or minimum dissipative in nature. Results indicate that the third-order method using minimum dispersion is the most accurate method tested. Computer time requirements are approximately twice those needed for second-order techniques.

Anderson, D.↗
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.↗