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 271 records · Page 15

Calculation of unsteady two-dimensional laminar and turbulent boundary layers with fluctuations in external velocity

A numerical method is proposed for calculating unsteady two-dimensional laminar and turbulent boundary layers with allowance for fluctuations in external flow velocity. An eddy-viscosity formulation is used to model the Reynolds shear stress term appropriate to turbulent flow, and the governing boundary layer equations are solved by a two-point finite-difference scheme. For illustration, the method is applied to the calculation of the phase angles between the shear stress at the wall and an oscillating external laminar boundary layer on a flat plate. The effectiveness of the method (for engineering purposes) is demonstrated by comparing calculations for turbulent flows with experimental data.

Cebeci, T.↗

Prospects for computational aerodynamics

The integral representations approach, for the solution of the Navier-Stokes equations is discussed as well as experience in its development and in applying available finite-difference and finite-element techniques to the treatment of three-dimensional problems, and the computation of turbulent flow. The magnitude of efforts required to develop turbulence models and three-dimensional algorithms indicates that the computational fluid dynamics research must have a broad base. Broader access to modern computing facilities that are in existence within NASA should be promoted for active researchers not directly affiliated with that agency.

Wu, J. C.↗

Comparison of two computer programs by predicting turbulent mixing of helium in a ducted supersonic airstream

Two parabolic flow computer programs, SHIP (a finite-difference program) and COMOC (a finite-element program), are used for predicting three-dimensional turbulent reacting flow fields in supersonic combustors. The theoretical foundation of the two computer programs are described, and then the programs are applied to a three-dimensional turbulent mixing experiment. The cold (nonreacting) flow experiment was performed to study the mixing of helium jets with a supersonic airstream in a rectangular duct. Surveys of the flow field at an upstream were used as the initial data by programs; surveys at a downstream station provided comparison to assess program accuracy. Both computer programs predicted the experimental results and data trends reasonably well. However, the comparison between the computations from the two programs indicated that SHIP was more accurate in computation and more efficient in both computer storage and computing time than COMOC.

Pan, Y. S.↗

Computational design of the basic dynamical processes of the UCLA general circulation model

The 12-layer UCLA general circulation model encompassing troposphere and stratosphere (and superjacent 'sponge layer') is described. Prognostic variables are: surface pressure, horizontal velocity, temperature, water vapor and ozone in each layer, planetary boundary layer (PBL) depth, temperature, moisture and momentum discontinuities at PBL top, ground temperature and water storage, and mass of snow on ground. Selection of space finite-difference schemes for homogeneous incompressible flow, with/without a free surface, nonlinear two-dimensional nondivergent flow, enstrophy conserving schemes, momentum advection schemes, vertical and horizontal difference schemes, and time differencing schemes are discussed.

Arakawa, A.↗

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

A decomposition algorithm is presented which 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.↗

A perturbation solution for interlaminar stresses in bidirectional laminates

The paper considers a thin elastic bidirectional symmetric composite laminate subjected to a uniform axial strain. All interlaminar stresses are determined by a zeroth-order perturbation analysis of the governing equations. The method of solution considers the laminate to be composed of an interior region (removed from the free edge) and a boundary layer region. The interior region is analyzed by dropping terms related to powers of the thickness-to-width ratio in the dimensionless elasticity equations. The boundary layer region is studied by introducing a stretching transformation to the governing equations. A uniform solution is then formed by satisfying Prandtl's matching principle of perturbation theory. Results for a four-layer bidirectional graphite-epoxy laminate are compared with existing finite-difference results. It is shown that the perturbation analysis provides more insight into the stress distribution near the free edge.

Hsu, P. W.↗

Calculation of transonic potential flowfields about complex, three-dimensional configurations

Methods for extending iterative, finite-difference calculations of transonic potential flowfields to complex three-dimensional configurations are discussed. One particularly attractive approach is to use relatively simple conformal mappings in combination with shearing transformations to generate computational domains that are nearly-conformally mapped from the physical space in one family of coordinate surfaces, and which map the complex boundaries to grid surfaces. The application of such a method to a general wing-body combination or to a multi-bladed fan is discussed. A transformation to map the wing-fuselage or fan-hub combination to a convenient computational domain is proposed. The transformation is useful in its own right for treating the two-dimensional problems of flow past a profile in a wind tunnel or through a cascade. Some results of preliminary calculations are presented.

Caughey, D. A.↗

Displacement thickness distributions in transonic flows about 3-D wings

The transonic 3-D inviscid small-perturbation solution of Bailey and Ballhaus is combined with a finite-difference solution for Prandtl's boundary-layer equations in order to include viscous effects. The inviscid-viscous interaction is modeled by means of the displacement surface, which can be thought of as the effective body surface seen by the inviscid flow. Displacement thickness, lift, and pressure distributions resulting from the combined solution are presented for transonic flows about the RAE 101 A wing and a Lockheed transport wing, both at small angles of attack. The influence of changing arbitrarily the start of transition on the displacement surface and lift is discussed for the RAE wing flow.

Kordulla, W.↗

Meteorite impact ejecta - Dependence of mass and energy lost on planetary escape velocity

The amounts of material and energy which escape a planet in a meteorite impact event is calculated as functions of impact and escape velocities. Results are obtained from the computed flow induced by the impact of iron and gabbroic anorthosite spheres onto a half-space of anorthosite at impact velocities of 5 to 45 km/sec. The impact-induced flows were determined by a numerical method using the mass, momentum, and energy conservation relations in finite-difference approximation, within an Eulerian computational grid. The impact velocities at which ejecta losses equal meteorite mass gains are found to be approximately 20, 35, and 45 km/sec for anorthosite objects and approximately 25, 35, and 40 km/sec for iron objects striking anorthosite surfaces for the gravity fields of the moon, Mercury and Mars.

Okeefe, J. D.↗

Comparison of numerical and experimental 'conical' flow fields in supersonic corners with compression and/or expansion

The flow field produced by the intersection of two plane solid surfaces in a supersonic stream is a complex interference flow. These flows can be fully compressive, fully expansive, or of mixed compression-expansion nature. This paper presents a comparison of the experimentally obtained flow-field structure in an axial corner with that predicted numerically by using a shock-capturing finite-difference method. The effect of sweep and surface deflection are evaluated, and the general influence of each is presented for the three classes of corner flow. The results show that the numerical method is a valuable aid in understanding the flow structure for simple configurations. In addition, confidence in the numerical method is gained for use in solving more general three-dimensional configurations where the flow is nonconical and several wave interaction may be presented.

Anderson, D. A.↗

Navier-Stokes solution for the laminar and turbulent supersonic flow near the trailing edge of a flat plate

The present paper deals with problems of two-dimensional laminar and turbulent flow of a viscous incompressible fluid near the trailing edge of a thin flat plate. The complete system of Navier-Stokes equations is solved by means of MacCormack's finite-difference scheme. In the laminar case, the results are compared with the triple deck solution of Daniels (1973). The comparisons indicate that the asymptotic triple deck theory for supersonic trailing-edge flow is accurate within five percent for Reynolds numbers greater than 1000. In the turbulent case, the Prandtl-Van Driest-Clausner algebraic eddy viscosity model is used, showing that the region of upstream influence is roughly of the order of the boundary layer thickness.

Hussaini, M. Y.↗

Impact-induced energy partitioning, melting, and vaporization on terrestrial planets

The reported investigation was conducted in connection with the desirability to obtain information regarding the response of a solidified planetary surface to the primary impact of a meteoroid whose dimensions are small compared to that of the planet. The investigation is an extension of studies conducted by Ahrens and O'Keefe (1977) concerned with a calculation via finite-difference techniques of the entire flow field, the energy partitioning, and the quantity of melt and vapor produced upon impact of an iron object on a gabbroic anorthosite half-space at 15 km/sec. The assumptions used in these calculations are summarized. Calculations based on an equation of state for the low- and high-pressure mineral assemblage of gabbroic anorthosite, initial density 2.94 g/cu cm are considered. An algorithm is formulated for the description of the phase transformation upon shock loading. It is shown that the relative volume of melt and the equivalent volume of shock-induced vapor is proportional to the meteoroid kinetic energy (at high speeds) using the similarity scaling of projectile and target properties as suggested by Dienes and Walsh (1970).

Okeefe, J. D.↗

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↗

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↗

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

Advances in Application of Fast Semidirect Computational Methods in Transonic Flow

This paper is intended as a review and summary of the advances made in a recently developed approach for rapid numerical solution of the equations of inviscid transonic aerodynamics. The investigation has been limited to two-dimensional, steady, inviscid flow over airfoils in a subsonic free stream, with emphasis on development of a rapid computational technique, rather than on generality of application. The approach uses finite-difference algorithms called "fast direct elliptic solvers" within an iteration scheme. "Direct" means that the entire computation field is solved at once, rather than in successive traverses over the field as in a point- or line-relaxation method. Such an iterative method is referred to as "semidirect." The iterative convergence can be faster than in other relaxation methods because changes are felt simultaneously at all points in each succeeding iteration. Direct elliptic solvers and semidirect methods have restrictions, but these are gradually being removed. Direct solvers were first developed for solving Poisson's equation on a rectangle without interior boundaries. A method to treat first-order systems, a direct Cauchy-Riemann solver has also been developed. Numerical treatment of part of a system of nonlinear equations by a Poisson solver has been reported. Also Poisson solvers in semidirect methods were used for nonseparable elliptic equations. The semidirect method was extended to the solution of a problem of mixed type, where the improved Murman-Cole transonic small-disturbance difference equations were solved. A slightly supercritical flow over a biconvex airfoil was treated successfully, but the iterations did not converge for more strongly supercritical conditions In another work the addition of terms ot both sides of the difference equations stabilized the iteration for supercritical conditions with large supersonic zones. For this, the Cauchy-Riemann solver was revised to incl,ude the needed terms. Most recently, the evaluation of parameters for rapid convergence and comparisons, with Murman's line-relaxation method was described. The method was extended to full second order accuracy in a fully conservative formulation in another work.

Martin, E. Dale↗