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 343 records · Page 19

Equatorial wave-mean flow interaction - A numerical study of the role of latitudinal shear

A time-dependent primitive equation model for an equatorial channel is used to assess the interaction of equatorial Kelvin and mixed Rossby-gravity waves with the mean flow. The proposed model involves a semiimplicit time-differencing scheme and a finite-difference grid in the meridional plane. It is shown that forced equatorial waves interact with mean flow to produce equatorial jets characterized by downward-moving westerly (Kelvin wave forcing) and easterly (mixed Rossby-gravity wave forcing) shear zones, respectively. For parameters characteristic of the observed waves in the equatorial stratosphere, the wave-mean flow interaction process always reduces the amplitude of any initial cross-equatorial mean wind shear. The mean flow profile tends to become symmetric about the equator as the interaction process continues.

Holton, J. R.↗

A new boundary-layer inviscid iteration technique for separated flow

A new viscous-inviscid iteration procedure is presented which is applicable to separated flows. The new procedure is simple, converges rapidly, and does not require numerical smoothing, at least in the cases computed thus far. Calculations are presented for both incompressible and compressible separated flow in the juncture region between an axisymmetric body and sting. The viscous computation is obtained with a recently developed inverse boundary-layer procedure which is no more complicated than a standard direct finite-difference boundary-layer technique. The inviscid computation is made with an axisymmetric transonic code called RAXBOD, in which surface injection is used to represent the viscous effects.

Carter, J. E.↗

A fast, conservative algorithm for solving the transonic full-potential equation

A fast, fully implicit approximate factorization (AF) algorithm designed to solve the conservative transonic full-potential equation in either two or three dimensions is described. The algorithm uses an upwind bias of the density coefficient for stability in supersonic regions. This provides an effective upwind difference of the streamwise terms for any orientation of the velocity vector (i.e., 'rotated differencing'), and thereby greatly enhances the reliability of the present algorithm. A numerical transformation is used to establish an arbitrary body-fitted finite-difference mesh. Computed results for both airfoils and simplified wings demonstrate substantial improvement in convergence speed for the new algorithm relative to standard successive-line overrelaxation algorithms.

Holst, T. L.↗

Evaluation of finite element formulations for transient conduction forced-convection analysis

Numerical studies clarifying the advantages and disavantages of conventional versus upwind convective finite elements are presented along with lumped versus consistent formulations for practical conduction forced-convection analysis. A finite-element procedure for treatment of negligible capacitance fluid nodes is presented. The procedure is based on procedures used in finite-element structural dynamics to treat nodes with negligible structural mass. Two finite-element programs and a finite-difference lumped-parameter program used in the studies are discussed. Evaluation studies utilizing three convection and two combined conduction-convection problems are then presented and discussed. Additionally, the computational time saving offered by the finite element procedure is considered for a practical combined conduction-convection problem.

Thornton, E. A.↗

Application of numerical methods to heat transfer and thermal stress analysis of aerospace vehicles

The paper describes a thermal-structural design analysis study of a fuel-injection strut for a hydrogen-cooled scramjet engine for a supersonic transport, utilizing finite-element methodology. Applications of finite-element and finite-difference codes to the thermal-structural design-analysis of space transports and structures are discussed. The interaction between the thermal and structural analyses has led to development of finite-element thermal methodology to improve the integration between these two disciplines. The integrated thermal-structural analysis capability developed within the framework of a computer code is outlined.

Wieting, A. R.↗

Analysis of Combustion Instability in Liquid Fuel Rocket Motors

The development of a technique to be used in the solution of nonlinear velocity-sensitive combustion instability problems is described. The orthogonal collocation method was investigated. It found that the results are heavily dependent on the location of the collocation points and characteristics of the equations, so the method was rejected as unreliabile. The Galerkin method, which has proved to be very successful in analysis of the pressure sensitive combustion instability was found to work very well. It was found that the pressure wave forms exhibit a strong second harmonic distortion and a variety of behaviors are possible depending on the nature of the combustion process and the parametric values involved. A one-dimensional model provides further insight into the problem by allowing a comparison of Galerkin solutions with more exact finite-difference computations.

Wong, K. W.↗

Calculation of two-dimensional inlet flow fields in a supersonic free stream by an implicit marching code with nonorthogonal mesh generation: User's manual

An implicit, shock-capturing finite-difference code which is used to calculate two-dimensional inlet flow fields in a supersonic free stream is explained. The Euler equations are subjected to general nonorthogonal transformation and a body-fitted coordinate system is employed. The mathematical formulation of the problem is given along with the numerical algorithm. Initial and boundary conditions, numerical stability, program limitations, and accuracy is discussed. An overall program logic as well as instructions for program use and operation are also furnished.

Biringen, S. H.↗

Amplitude effects on the dynamic performance of hydrostatic gas thrust bearings

A strip gas film bearing with inherently compensated inlets is analyzed to determine the effect of disturbance amplitude on its dynamic performance. The governing Reynolds' equation is solved using finite-difference techniques. The time dependent load capacity is represented by a Fourier series up to and including the third harmonics. For the range of amplitudes investigated the linear stiffness was independent of the amplitude, and the linear damping was inversely proportional to (1 - epsilon-squared) to the 1.5 power where epsilon is the amplitude relative to the film thickness.

Stiffler, A. K.↗

Finite-element numerical modeling of atmospheric turbulent boundary layer

A dynamic turbulent boundary-layer model in the neutral atmosphere is constructed, using a dynamic turbulent equation of the eddy viscosity coefficient for momentum derived from the relationship among the turbulent dissipation rate, the turbulent kinetic energy and the eddy viscosity coefficient, with aid of the turbulent second-order closure scheme. A finite-element technique was used for the numerical integration. In preliminary results, the behavior of the neutral planetary boundary layer agrees well with the available data and with the existing elaborate turbulent models, using a finite-difference scheme. The proposed dynamic formulation of the eddy viscosity coefficient for momentum is particularly attractive and can provide a viable alternative approach to study atmospheric turbulence, diffusion and air pollution.

Lee, H. N.↗

Laminar and turbulent flows over a spherically blunted cone with massive surface blowing

Numerical solutions are presented for the flow over a spherically blunted cone with massive surface blowing. Time-dependent viscous shock layer equations are used to describe the flowfield. The boundary conditions on the body surface include a prescribed blowing rate distribution. The governing equations are solved by a time-asymptotic finite-difference method. Results presented here are only for a perfect gas-type flow at zero angle of attack. Both laminar and turbulent flow solutions are obtained. It is found that the surface blowing smooths out the effect of the curvature discontinuity at the sphere-cone juncture point on the laminar flowfield and results in a negative pressure gradient over the body. The shock slope increases on the downstream portion of the body as the surface blowing rate is increased. The turbulent flow with surface blowing is found to redevelop a boundary-layer-like region near the surface. The effects of this boundary-layer-like region on the flowfield and heating rates are discussed.

Kumar, A.↗

An efficient explicit-implicit-characteristic method for solving the compressible Navier-Stokes equations

Explicit, implicit, and characteristic finite-difference methods are applied to solve model equations representative of the compressible Navier-Stokes equations. An approach is then formulated for solving the Navier-Stokes equation at high Reynolds numbers. The approach has drastically reduced the computation time required to obtain viscous flow solutions. Computational results for shock wave separated flows are presented.

Maccormack, R. W.↗

Computational simulation of corotational vortex merger using 0, 1, and 2 equation turbulence models

A computational study has been conducted to gain a better understanding of the role of turbulence in the merger of corotating vortices and in the decay of a single vortex. The flow field is simulated by using an inviscid flow model and zero-, one-, and two-equation turbulent-flow models. An implicit finite-difference procedure is employed to integrate the unsteady, two-dimensional, incompressible flow equations in a cross-plane. In the zero- and one-equation models, a spatially varying mixing-length model that simulates the streamline curvature effect is used to define an eddy viscosity. The computational merging distance criterion is in good agreement with the experimental one. In the case of the two-equation model, the turbulent kinetic energy and a modified rate of dissipation equation incorporating a curvature correction term are solved. A comparative study of the computational results corresponding to the various flow models and the available experimental data is presented.

Raj, P.↗

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↗