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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1975

Application of a transitional boundary-layer theory in the low hypersonic Mach number regime

An investigation is made to assess the capability of a finite-difference boundary-layer procedure to predict the mean profile development across a transition from laminar to turbulent flow in the low hypersonic Mach-number regime. The boundary-layer procedure uses an integral form of the turbulence kinetic-energy equation to govern the development of the Reynolds apparent shear stress. The present investigation shows the ability of this procedure to predict Stanton number, velocity profiles, and density profiles through the transition region and, in addition, to predict the effect of wall cooling and Mach number on transition Reynolds number. The contribution of the pressure-dilatation term to the energy balance is examined and it is suggested that transition can be initiated by the direct absorption of acoustic energy even if only a small amount (1 per cent) of the incident acoustic energy is absorbed.

Shamroth, S. J.↗

NASA NTRS1975

Calculation of three-dimensional compressible laminar and turbulent boundary layers. Calculation of three-dimensional compressible boundary layers on arbitrary wings

A very general method for calculating compressible three-dimensional laminar and turbulent boundary layers on arbitrary wings is described. The method utilizes a nonorthogonal coordinate system for the boundary-layer calculations and includes a geometry package that represents the wing analytically. In the calculations all the geometric parameters of the coordinate system are accounted for. The Reynolds shear-stress terms are modeled by an eddy-viscosity formulation developed by Cebeci. The governing equations are solved by a very efficient two-point finite-difference method used earlier by Keller and Cebeci for two-dimensional flows and later by Cebeci for three-dimensional flows.

Cebeci, T.↗

NASA NTRS1975

Asymptotic theory of two-dimensional trailing-edge flows

Problems of laminar and turbulent viscous interaction near trailing edges of streamlined bodies are considered. Asymptotic expansions of the Navier-Stokes equations in the limit of large Reynolds numbers are used to describe the local solution near the trailing edge of cusped or nearly cusped airfoils at small angles of attack in compressible flow. A complicated inverse iterative procedure, involving finite-difference solutions of the triple-deck equations coupled with asymptotic solutions of the boundary values, is used to accurately solve the viscous interaction problem. Results are given for the correction to the boundary-layer solution for drag of a finite flat plate at zero angle of attack and for the viscous correction to the lift of an airfoil at incidence. A rational asymptotic theory is developed for treating turbulent interactions near trailing edges and is shown to lead to a multilayer structure of turbulent boundary layers. The flow over most of the boundary layer is described by a Lighthill model of inviscid rotational flow. The main features of the model are discussed and a sample solution for the skin friction is obtained and compared with the data of Schubauer and Klebanoff for a turbulent flow in a moderately large adverse pressure gradient.

Melnik, R. E.↗

NASA NTRS1975

Subgrid or Reynolds stress-modeling for three-dimensional turbulence computations

A review is given of recent advances in two distinct computational methods for evaluating turbulence fields, namely, statistical Reynolds stress modeling and turbulence simulation, where large eddies are followed in time. It is shown that evaluation of the mean Reynolds stresses, rather than use of a scalar eddy viscosity, permits an explanation of streamline curvature effects found in several experiments. Turbulence simulation, with a new volume averaging technique and third-order accurate finite-difference computing is shown to predict the decay of isotropic turbulence in incompressible flow with rather modest computer storage requirements, even at Reynolds numbers of aerodynamic interest.

Rubesin, M. W.↗

NASA NTRS1975

Time-dependent transonic flow solutions for axial turbomachinery

Three-dimensional unsteady transonic flow through an axial turbomachine stage is described in terms of a pair of two-dimensional formulations pertaining to orthogonal surfaces, namely, a blade-to-blade surface and a hub-to-casing surface. The resulting systems of nonlinear, inviscid, compressible equations of motion are solved by an explicit finite-difference technique. The blade-to-blade program includes the periodic interaction between rotor and stator blade rows. Treatment of the boundary conditions and of the blade slipstream motion by a characteristic type procedure is discussed in detail. Harmonic analysis of the acoustic far field produced by the blade row interaction, including an arbitrary initial transient, is outlined. Results from the blade-to-blade program are compared with experimental measurements of the rotating pressure field at the tip of a high-speed fan. The hub-to-casing program determines circumferentially averaged flow properties on a meridional plane. Blade row interactions are neglected in this formulation, but the force distributions over the entire blade surface for both the rotor and stator are obtained. Results from the hub-to-casing program are compared with a relaxation method solution for a subsonic rotor. Results are also presented for a quiet fan stage which includes transonic flow in both the rotor and stator and a normal shock in the stator.

Erdos, J.↗

NASA NTRS1975

A comparison of a shock-capturing technique with experimental data for three-dimensional internal flows

Shock-capturing solutions for an axisymmetric supersonic inlet at small angles of attack are obtained. Good overall agreement between the shock-capturing solutions and experimental data is shown except in regions of strong viscous effects or boundary-layer removal. Although the results indicate a strong potential for the use of shock-capturing or finite-difference solutions for internal flows, improvement in the ability to handle the reflection of strong shockwaves having downstream Mach numbers near 1 is needed.

Presley, L. L.↗

NASA NTRS1975

Internal and external axial corner flows

The inviscid, internal, and external axial corner flows generated by two intersecting wedges traveling supersonically are obtained by use of a second-order shock-capturing, finite-difference approach. The governing equations are solved iteratively in conical coordinates to yield the complicated wave structure of the internal corner and the simple peripheral shock of the external corner. The numerical results for the internal flows compare favorably with existing experimental data.

Kutler, P.↗

NASA NTRS1975

Implicit approximate-factorization schemes for the low-frequency transonic equation

Two- and three-level implicit finite-difference algorithms for the low-frequency transonic small disturbance-equation are constructed using approximate factorization techniques. The schemes are unconditionally stable for the model linear problem. For nonlinear mixed flows, the schemes maintain stability by the use of conservatively switched difference operators for which stability is maintained only if shock propagation is restricted to be less than one spatial grid point per time step. The shock-capturing properties of the schemes were studied for various shock motions that might be encountered in problems of engineering interest. Computed results for a model airfoil problem that produces a flow field similar to that about a helicopter rotor in forward flight show the development of a shock wave and its subsequent propagation upstream off the front of the airfoil.

Ballhaus, W. F.↗

NASA NTRS1976

Historical evolution of vortex-lattice methods

A review of the beginning and some orientation of the vortex-lattice method were given. The historical course of this method was followed in conjunction with its field of computational fluid dynamics, spanning the period from L.F. Richardson's paper in 1910 to 1975. The following landmarks were pointed out: numerical analysis of partial differential equations, lifting-line theory, finite-difference method, 1/4-3/4 rule, block relaxation technique, application of electronic computers, and advanced panel methods.

Deyoung, J.↗

NASA NTRS1975

Internal flow calculations for axisymmetric supersonic inlets at angle of attack

Several solutions for the internal flow in axisymmetric supersonic inlets at angle of attack were obtained. A shock-capturing finite-difference technique was used to obtain the solutions. Good overall agreement between the shock-capturing solutions and experimental data was obtained, except in regions of strong viscous effects or boundary-layer removal. The required centerbody translation with changes in freestream Mach number and angle of attack was obtained and agrees qualitatively with experimental data.

Presley, L. L.↗

NASA NTRS1975

Numerical solutions of the unsteady Navier-Stokes equations for arbitrary bodies using boundary-fitted curvilinear coordinates

A method of automatic body-fitted curvilinear coordinate generation is described and used to construct a finite-difference solution of the full incompressible time-dependent Navier-Stokes equations for the unsteady laminar viscous flow arbitrary two-dimensional airfoils or any other two-dimensional body. A method of controlling the spacing of the coordinate lines encircling the body is developed in order to treat high Reynolds number flows, since the coordinate lines must concentrate near the surface to a greater degree as the Reynolds number increases. Multiple airfoils and submerged hydrofoils are treated as illustrative examples. The solution shows good agreement with the Blasius boundary layer solution for the flow past a semi-infinite flat plate.

Thompson, J. F.↗

NASA NTRS1975

Diagonalization and simultaneous symmetrization of the gas-dynamic matrices

The hyperbolic nature of the unsteady, inviscid, gas-dynamic equations implies the existence of a similarity transformation for diagonalizing an arbitrary linear combination of coefficient matrices. It is shown that the individual matrices are simultaneously symmetrized by the similarity transformation. The transformations and their norms can be applied to the well-posedness of the Cauchy problem, linear stability theory for finite-difference approximations, and simplification of block-tridiagonal systems that arise in implicit time-split algorithms.

Warming, R. F.↗

NASA NTRS1975

A two-level trajectory decomposition algorithm featuring optimal intermediate target selection

A decomposition algorithm is presented that optimizes complex missions by partitioning the trajectory into natural segments such as ascent or entry. Each segment defines a full-rank targeting subproblem. These are solved sequentially using the Newton-Raphson algorithm. The master problem, representing the complete mission, is to determine subproblem targets and master-problem controls that optimize the mission objective subject to intersegment constraints. The gradient projection algorithm solves this problem using derivatives obtained analytically from finite-difference subproblem sensitivities. Thus, the mission is optimized by coordinating the solution of tractible subproblems. Computational results for a synchronous equatorial mission are included.

Petersen, F. M.↗

NASA NTRS1975

Multiple slot skin friction reduction

The paper investigates analytically the effect of multiple slot injection on skin friction for a representative fuselage shape (ogive-cylinder body) and evaluates the potential of slot injection as a drag reduction system in subsonic flow. Typical CTOL cruise flight conditions (Mach number equals 0.82 at altitudes of 11 km) were adopted for a fuselage 67.06 m in length and with maximum diameter of 7.32 m. The numerical method of Price and Harris (1972) was used to calculate the boundary-layer characteristics up to the first slot, while the finite-difference method of Beckwith and Bushnell (1971) was used to calculate the velocity profile downstream of one, three, five, or ten slots. An integral expression is proposed for characterizing skin friction reduction effectiveness, and it is seen that large reductions in viscous drag (50%) are available through slot injection. Skin friction reduction is improved by increasing the number of injection slots but at a diminishing rate.

Howard, F. G.↗

NASA NTRS1975

Computation of three-dimensional, inviscid supersonic flows

The paper sets forth in detail a method for the finite-difference computation of three-dimensional supersonic fields in an Eulerian mesh. First-, second-, and third-order finite difference schemes are examined. Attention is given to proper treatment of the impermeable and permeable boundaries encompassing the computational plane. Numerical results are presented for certain specific configurations: a conical wing-body combination, internal corner flow, a two-dimensional blunt body, an interfering shock problem, and three-dimensional inviscid supersonic flow past a shuttle-orbiter type vehicle.

Kutler, P.↗

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↗

NASA NTRS1950

Direct method of design and stress analysis of rotating disks with temperature gradient

A method is presented for the determination of the contour of disks, typified by those of aircraft gas turbines, to incorporate arbitrary elastic-stress distributions resulting from either centrifugal or combined centrifugal and thermal effects. The specified stress may be radial, tangential, or any combination of the two. Use is made of the finite-difference approach in solving the stress equations, the amount of computation necessary in the evolution of a design being greatly reduced by the judicious selection of point stations by the aid of a design chart. Use of the charts and of a preselected schedule of point stations is also applied to the direct problem of finding the elastic and plastic stress distribution in disks of a given design, thereby effecting a great reduction in the amount of calculation. Illustrative examples are presented to show computational procedures in the determination of a new design and in analyzing an existing design for elastic stress and for stresses resulting from plastic flow.

Manson, S S↗

NASA NTRS1961

Cross-Sectional Deformations of Monocoque Beams and Their Effects on the Natural Vibration Frequencies

The variational principle, differential equations, and boundary conditions governing the cross-sectional distortions due to inertia loading of a two-dimensional model of a thin monocoque wing are shown. A theoretical analysis of this simplified model is made in order to determine the nature of the coupling between the cross-sectional modes and the spanwise deformation modes. General solutions are obtained in finite-difference form for arbitrary cross sections and an exact solution is presented for a parabolic-arc cross section of constant cover thickness. The application of these results in evaluating the coupled frequencies of the actual structure is discussed. Frequencies evaluated for a parabolic-arc monocoque beam show good agreement with experimental values.

Thomson, Robert G.↗