SEARCH · Engineering Papers

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 1,729 records · Page 96

All sources NASA NTRS DOE OSTI

NASA NTRS1976

Impact ejecta on the moon

The response of a lunar-sized object to the impact of meteoroids no more than about 100 km in radius is studied by means of a numerical model. The partitioning of impact energy into the kinetic and internal energy of the ejecta is obtained by using the conservation of mass, momentum, and energy conservation equations in finite-difference form within an Eulerian framework with approximate equations of state. The calculations are performed for a 15 km/sec impact of an iron object 5 cm in radius on a gabbroic anorthosite surface. Ejecta ballistic analysis is then performed. Most of the material lost escaping the moon is lunar crust material. Only 0.2% of the meteoroid escapes, all in the vapor phase.

Okeefe, J. D.↗

NASA NTRS1977

An implicit factored scheme for the compressible Navier-Stokes equations

An implicit finite difference scheme is developed for the numerical solution of the compressible Navier-Stokes equations in conservation-law form. The algorithm is second-order-time accurate, noniterative, and spatially factored. In order to obtain an efficient factored algorithm, the spatial cross-derivatives are evaluated explicitly. However, the algorithm is unconditionally stable and, although a three-time-level scheme, requires only two-time-levels of data storage. The algorithm is constructed in a 'delta' form (i.e., increments of the conserved variables and fluxes) that provides a direct derivation of the scheme and leads to an efficient computational algorithm. In addition, the delta form has the advantageous property of a steady-state (if one exists) independent of the size of the time step. Numerical results are presented for a two-dimensional shock boundary-layer interaction problem.

Beam, R. M.↗

NASA NTRS1977

A numerical study of the unsteady leading edge separation bubble on an oscillating airfoil

A preliminary study of the unsteady viscous flow in the region of an airfoil leading edge was performed, in which the interaction between the viscous and inviscid flow fields is neglected. The solution method uses the finite difference form of the governing equations throughout the separated flow field and incorporates a transition model based on the integral turbulence kinetic energy equation. The validity of the numerical procedure is verified by making comparisons with analytical solutions to several test problems, including unsteady flow over a plate oscillating in its own plane. The method was then applied to the problem of unsteady viscous flow over a NACA 0012 airfoil oscillating sinusoidally in pitch. The flow field characteristics were in qualitative agreement with experimental results. The bubble moved forward on the airfoil and decreased in size as incidence was increased. Viscous flow in the leading edge region was found to be quasi-steady, while bubble height varied inversely with Reynolds number.

Kreskovsky, J. P.↗

NASA NTRS1977

Turbulent plane Couette flow using probability distribution functions

A numerical scheme employing a combination of the discrete ordinate method and finite differences is developed for solving the one-dimensional form of Lundgren's (1967) model equation for turbulent plane Couette flow. The approach used requires no a priori assumption about the form of the turbulent distribution function, and the numerical solution is obtained directly from the governing differential equations. Two different types of boundary conditions (zero-gradient and Chapman-Enskog) for the distribution function are evaluated by comparing the numerical results with experimental data. It is found that: (1) the present approach gives convergent and stable results over a wide range of Reynolds numbers; (2) Lundgren's equation yields results that compare well with experimental data for mean velocity and skin friction in the case of simple Couette flow; (3) the zero-gradient boundary condition leads to a logarithmic flow profile; and (4) the Chapman-Enskog boundary condition provides very good agreement with experimental data when applied within the near-wall region.

Srinivasan, R.↗

NASA NTRS1977

Compressible laminar boundary layers with suction on swept and tapered wings

In this paper we present a numerical method for solving the compressible laminar boundary-layer equations with suction on swept and tapered wings. The method employs an efficient two-point finite-difference method to solve the governing equations, and a very convenient similarity transformation which removes the wall normal velocity as a boundary condition and places it into the governing equations as a parameter. In this way the awkward nonlinear boundary condition which couples all the variables is avoided. To test and demonstrate the method, we present a sample calculation for a typical laminar-flow-control (LFC) wing.

Kaups, K.↗

NASA NTRS1976

Numerical solution of the steady-state Navier-Stokes equations for hypersonic flow about blunt axisymmetric bodies

The steady-state Navier-Stokes equations are solved for hypersonic flow about blunt axisymmetric bodies. The equations of motion are solved by successive approximations using an implicit finite-difference scheme. The results are compared with viscous shock-layer theory, experimental data, and time-dependent solutions of the Navier-Stokes equations. It is demonstrated that viscous shock-layer theory is sufficiently accurate for the range of flight conditions normally encountered by entry vehicles.

Anderson, E. C.↗

NASA NTRS1977

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

NASA NTRS1977

Polynomial interpolation methods for viscous flow calculations

Higher-order collocation procedures which result in block-tridiagonal matrix systems are derived from (1) Taylor series expansions and from (2) polynomial interpolation, and the relationships between the two formulations, called respectively Hermite and spline collocation, are investigated. A Hermite block-tridiagonal system for a nonuniform mesh is derived, and the Hermite approach is extended in order to develop a variable-mesh sixth-order block-tridiagonal procedure. It is shown that all results obtained by Hermite development can be recovered by appropriate spline polynomial interpolation. The additional boundary conditions required for these higher-order procedures are also given. Comparative solutions using second-order accurate finite difference and spline and Hermite formulations are presented for the boundary layer on a flat plate, boundary layers with uniform and variable mass transfer, and the viscous incompressible Navier-Stokes equations describing flow in a driven cavity.

Rubin, S. G.↗

NASA NTRS1977

TOMCAT - A code for numerical generation of boundary-fitted curvilinear coordinate systems on fields containing any number of arbitrary two-dimensional bodies

A method for automatic generation of boundary-fitted curvilinear coordinate systems, where the transformed coordinates are solutions of an elliptic differential system in the physical plane, and where the coordinate lines are coincident with all boundaries of a general multiply-connected, two-dimensional region containing any number of arbitrarily shaped bodies, and is described along with a suitable computer code for implementing the method. Any partial differential system can be solved on the boundary-fitted coordinate system by appropriate transformations. The transformed equations are approximated by finite differences and solved numerically in the transformed plane. All computations, whether for generating coordinate system or then solving the transformed equations, can be done on a rectangular field with square mesh with no interpolation required on the boundaries. The physical boundaries may even be time-dependent.

Thompson, J. F.↗

NASA NTRS1976

Analysis of transonic flow over lifting and oscillating airfoils

Finite element procedures are presented for the numerical solution of transonic flow over lifting and oscillating airfoils based on small disturbance theory. The algorithm is simple and effective in handling the mixed flow problem and, for airfoil in free air, treats the entire flow domain by making use of the far field asymptotic solution. Numerical results are presented and compared with experimental data and those obtained by finite difference relaxation techniques.

Chan, S. T. K.↗

NASA NTRS1978

Flow velocities and streamlines

Computer program calculates subsonic or transonic flow on hubshroud, midchannel, stream surface of single-blade row of turbomachine. Program uses finite-different and quasi-orthogonal (velocity-gradient) methods. Program is reported in two volumes: Part I is User's Manual, Part II is Programmer's Manual.

Katsanis, T.↗

NASA NTRS1978

Internal and external 2-d boundary layer flows

Computer program computes general two dimensional turbulent boundary-layer flow using finite-difference techniques. Structure allows for user modification to accommodate unique problems. Program should prove useful in many applications where accurate boundary-layer flow calculations are required.

Crawford, M. E.↗

NASA NTRS1977

Analysis of simple 2-D and 3-D metal structures subjected to fragment impact

Theoretical methods were developed for predicting the large-deflection elastic-plastic transient structural responses of metal containment or deflector (C/D) structures to cope with rotor burst fragment impact attack. For two-dimensional C/D structures both, finite element and finite difference analysis methods were employed to analyze structural response produced by either prescribed transient loads or fragment impact. For the latter category, two time-wise step-by-step analysis procedures were devised to predict the structural responses resulting from a succession of fragment impacts: the collision force method (CFM) which utilizes an approximate prediction of the force applied to the attacked structure during fragment impact, and the collision imparted velocity method (CIVM) in which the impact-induced velocity increment acquired by a region of the impacted structure near the impact point is computed. The merits and limitations of these approaches are discussed. For the analysis of 3-d responses of C/D structures, only the CIVM approach was investigated.

Witmer, E. A.↗

NASA NTRS1977

An investigation of temperature distribution in cooled guide vanes

A numerical study to determine the temperature distribution in the guide vane blades of a radial inflow turbine is presented. A computer program was developed which permits the temperature distribution to be calculated when the blade is cooled internally using a combination of impingement and film cooling techniques. The study is based on the use of the finite difference method in a two dimensional heat conduction problem. The results are then compared to determine the best cooling configuration for a certain coolant to primary mass flow ratio.

Kotwal, R.↗

NASA NTRS1977

HCMM energy budget data as a model input for assessing regions of high potential groundwater pollution

The author has identified the following significant results. The finite difference model was used to calculate the differences in surface temperature between two hypothetical sites which result from a temperature difference at 50 cm due to the presence of shallow ground water at one of the sites. Although qualitative results of the model seemed consistant with experimental results, further evaluation showed a need for taking account of differences in thermal conductivity due to different moisture profiles at the two sites considered.

Moore, D. G.↗

NASA NTRS1977

Neutrally stable atmospheric flow over a two-dimensional rectangular block

The phenomena of atmospheric flow over a two dimensional surface obstruction such as a building modeled as a rectangular block are analyzed by an approach using the Navier-Stokes equations with a two equation model of turbulence. The partial differential equations for the vorticity, stream function, turbulence kinetic energy, and turbulence length scale are solved by a finite difference technique. The predicted results are in agreement with the limited experimental data available. Current computed results show that the separation bubble originates from the upper front corner of the block and extends approximately 11.5 block heights behind the block. The decay of the mean velocity along the wake center line coincides almost perfectly with the experimental data. The vertical profiles of the mean velocity defect are also in reasonable agreement with wind tunnel results. Velocity profiles in the mixing region are shown to agree with the error function profile typically found in the shear layer. Details of the behavior of the turbulence kinetic energy and the turbulence length scale are also discussed.

Shieh, C. F.↗

NASA NTRS1977

A solution to the Navier-Stokes equations based upon the Newton Kantorovich method

An implicit finite difference scheme based on the Newton-Kantorovich technique was developed for the numerical solution of the nonsteady, incompressible, two-dimensional Navier-Stokes equations in conservation-law form. The algorithm was second-order-time accurate, noniterative with regard to the nonlinear terms in the vorticity transport equation except at the earliest few time steps, and spatially factored. Numerical results were obtained with the technique for a circular cylinder at Reynolds number 15. Results indicate that the technique is in excellent agreement with other numerical techniques for all geometries and Reynolds numbers investigated, and indicates a potential for significant reduction in computation time over current iterative techniques.

Davis, J. E.↗

NASA NTRS1977

An automated procedure for calculating system matrices from perturbation data generated by an EAI Pacer and 100 hybrid computer system

Techniques are presented for determining the elements of the A, B, C, and D state variable matrices for systems simulated on an EAI Pacer 100 hybrid computer. An automated procedure systematically generates disturbance data necessary to linearize the simulation model and stores these data on a floppy disk. A separate digital program verifies this data, calculates the elements of the system matrices, and prints these matrices appropriately labeled. The partial derivatives forming the elements of the state variable matrices are approximated by finite difference calculations.

Milner, E. J.↗