Engineering Papers⌕ Search

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.

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

At least 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

Relaxation solutions for inviscid axisymmetric transonic flow over blunt or pointed bodies.

A finite-difference relaxation method is presented for numerical solution of the full potential equation and exact boundary conditions for general axisymmetric bodies is inviscid, steady transonic flow. Body-normal coordinates are used in the nose region and sheared cylindrical coordinates are used on the afterbody to accommodate corners such as boattails and flares. An improved difference scheme is used which does not require that the flow be nearly alined with a coordinate direction in supersonic regions, and which treats either subsonic or supersonic free streams. Numerical results are illustrated for some simple classical shapes such as spheres and ellipsoids, and for more practical shapes like tangent-ogives with boattails. Special attention is given to bodies which have been studied for area-rule applications. Agreement with available experimental results is good in cases where viscous effects and wind-tunnel wall interference are not important.

South, J. C., Jr.↗
NASA NTRS1973

Viscous flow over spinning cones at angle of attack.

A numerical finite-difference method is developed for evaluating the Magnus coefficients on spinning cones in laminar flow. The merged layer, the strong interaction region, and the downstream boundary layer are all considered. The numerical method is a predictor-corrector scheme developed for three-dimensional flows with or without crossflow diffusion. This method is particularly useful in problems in which a symmetry plane does not exist. Several contributions to the Magnus force and moments are considered. These include asymmetries in displacement thickness, centrifugal force and crossflow shear, and the effects of crossflow separation and vortex formation. Comparisons are made with experimental data and other analyses.

Lin, T. C.↗
NASA NTRS1973

Numerical solution for the inviscid supersonic flow in the corner formed by two intersecting wedges.

The inviscid, interference corner flow generated by two intersecting wedges immersed in a supersonic stream is obtained by use of a second-order, shock-capturing, finite-difference approach. The governing equations are solved iteratively in conical coordinates to yield the flow structure consisting of multiple shock and slip surfaces. The numerical results for shock wave and slip surface position and structure, pitot pressure traverses, and surface pressure distributions are compared with experimental data obtained over a wide range of Reynolds numbers. The comparisons show the best agreement with the high Reynolds number (greater than 3,000,000) results for which the boundary layer is turbulent.

Kutler, P.↗
NASA NTRS1973

Apollo 15 measurement of lunar surface brightness temperatures - Thermal conductivity of the upper 1.5 meters of regolith.

In situ measurements of lunar surface brightness temperatures made as a part of the Apollo Lunar Surface Experiments Package at the Apollo 15 Hadley Rille landing site are reported. Data derived from five thermocouples of the Heat Flow Experiment, which are lying on or just above the surface, are used to examine the thermal properties of the upper 15 cm of the lunar regolith using eclipse and nighttime cool-down temperatures. Application of finite-difference techniques in modeling the lunar soil shows that the thermocouple data are best fit by a model consisting of a low-density and low-thermal conductivity surface layer approximately 2 cm thick overlying a region increasing in conductivity and density with depth. Conductivities on the order of 0.00001 W per cm per deg K are postulated for the upper layer, with conductivity increasing to the order of 0.0001 W per cm per deg K at depths exceeding 20 cm. An increase in mean temperature with depth indicates that the ratio of radiative to conductive transfer at 350 K is 2.7 for at least the upper few centimeters of lunar soil; this value is nearly twice that measured for returned lunar fines.

Keihm, S. J.↗
NASA NTRS1974

Computer program for calculating laminar, transitional, and turbulent boundary layers for a compressible axisymmetric flow

Finite-difference computer program calculates viscous compressible boundary layer flow over either planar or axisymmetric surfaces. Flow may be initially laminar and progress through transitional zone to fully turbulent flow, or it may remain laminar, depending on imposed boundary conditions, laws of viscosity, and numerical solution of momentum and energy equations.

Albers, J. A.↗
NASA NTRS1972

Thermochemical nonequilibrium in atomic hydrogen at elevated temperatures

A numerical study of the nonequilibrium flow of atomic hydrogen in a cascade arc was performed to obtain insight into the physics of the hydrogen cascade arc. A rigorous mathematical model of the flow problem was formulated, incorporating the important nonequilibrium transport phenomena and atomic processes which occur in atomic hydrogen. Realistic boundary conditions, including consideration of the wall electrostatic sheath phenomenon, were included in the model. The governing equations of the asymptotic region of the cascade arc were obtained by writing conservation of mass and energy equations for the electron subgas, an energy conservation equation for heavy particles and an equation of state. Finite-difference operators for variable grid spacing were applied to the governing equations and the resulting system of strongly coupled, stiff equations were solved numerically by the Newton-Raphson method.

Scott, R. K.↗
NASA NTRS1973

View factor computer program (program view) user's manual

The purpose of program VIEW is to compute view factors between specified surfaces and to be compatible with level 15.5 of the NASTRAN structural analysis program. Program VIEW is a modification of a (finite element) view factor computation program called RAVFAC. VIEW is designed to run on an IBM System/360 operating under OS (operating system), with a minimum region size of 110 K bytes. The actual computation of view factors is still performed exactly as it was in the original version of RAVFAC. In developing VIEW, RAVFAC was modified to satisfy the following compatibility requirements: (1) accept finite element input which can also be used as input to NASTRAN, (2) produce output (view factors) in a format which can be used as input to NASTRAN, and (3) follow NASTRAN program design so that in the future VIEW can be incorporated into NASTRAN as a subroutine. The VIEW program permits computation of the view factors between surfaces, taking into account the presence of any intermediate surfaces. VIEW also computes these view factors either by contour integration or by finite difference (double summation) methods.

Puccinelli, E. F.↗
NASA NTRS1974

Decay of far flow field in trailing vortices

A finite difference machine code is used in the wake vortex problems in the quasi-cylindrical boundary layer approximation. A turbulent energy model containing new features is developed that accounts for the major effects disclosed by more advanced models in which the parameters are not yet established. Several puzzles that arose in previous theoretical investigations of wake vortices are resolved.

Baldwin, B. S.↗
NASA NTRS1974

A rapid numerical solution to subsonic flow over planar and axisymmetric profiles at an angle of attack of 0 deg

A numerical solution is presented for the incompressible flow over thin planar and axisymmetric profiles at an angle of attack of 0 degrees. The method uses a finite-difference field solution to the governing equation with a Gauss-Seidel successive overrelaxation scheme. However, the use of a simple Cartesian grid system restricts this method to slender profiles. Results are presented for a cambered airfoil, airfoil in wall effect (two-dimensional flowthrough inlet), body of revolution, and flowthrough nacelle. A computer program is presented which can be used for any of the previously mentioned cases with simple input changes. Results for compressible flow are available with the use of the appropriate two-dimensional or axisymmetric compressibility corrections. Computational time for a typical field calculation of 3000 grid points and 200 cycles through the field is less than 1 minute with less than 50,000 octal storage on the Control Data Corporation 6600 computing system.

Mascitti, V. R.↗
NASA NTRS1974

The computation of steady nozzle flow by a time-dependent method

The steady flow in two-dimensional and axisymmetric nozzles was computed using a time-dependent method. In this method the interior mesh points were computed using the MacCormack finite-difference scheme, while a characteristic scheme was used to calculate the boundary mesh points. No explicit artificial viscosity term was included. The fluid was assumed to be a perfect gas. This method was used to compute the flow in a 45 deg - 15 deg conical, converging-diverging nozzle, a 15 deg conical, converging nozzle, and a 10 deg conical, plug nozzle. Good agreement between the numerical solution and experimental data was found. In contrast to previous time-dependent methods, the computational times were less than one minute on a CDC 6600 computer.

Cline, M. C.↗
NASA NTRS1974

Numerical calculations of two dimensional, unsteady transonic flows with circulation

The feasibility of obtaining two-dimensional, unsteady transonic aerodynamic data by numerically integrating the Euler equations is investigated. An explicit, third-order-accurate, noncentered, finite-difference scheme is used to compute unsteady flows about airfoils. Solutions for lifting and nonlifting airfoils are presented and compared with subsonic linear theory. The applicability and efficiency of the numerical indicial function method are outlined. Numerically computed subsonic and transonic oscillatory aerodynamic coefficients are presented and compared with those obtained from subsonic linear theory and transonic wind-tunnel data.

Beam, R. M.↗
NASA NTRS1974

Estimation of deviation angle for axial-flow compressor blade sections using inviscid-flow solutions

Development of a method of estimating deviation angles by analytical procedures was begun. Solutions for inviscid, irrotational flow in the blade-to-blade plane were obtained with a finite-difference calculation method. Deviation angles for a plane cascade with a rounded trailing edge were estimated by using the inviscid-flow solutions and three trailing-edge hypotheses. The estimated deviation angles were compared with existing experimental data over a range of incidence angles at inlet flow angles of 30 deg and 60 deg. The results indicate that deviation angles can be estimated accurately (within 1 deg) by using one of the three trailing-edge hypotheses, but only when pressure losses are low. A new trailing-edge hypotheses is presented which is suitable (for the cascade considered) for both low- and high-loss operating points.

Miller, M. J.↗
NASA NTRS1974

Two-demensional analysis of heat and mass transfer in porous media using the strongly implicit procedure

Numerical results of the heat and mass transfer in a porous matrix are presented. The coupled, nonlinear partial differential equations describing this physical phenomenon are solved in finite difference form for two dimensions, using a new iterative technique (the strongly implicit procedure). The influence of the external environment conditions (heating and pressure) is shown to produce two-dimensional flow in the porous matrix. Typical fluid and solid temperature distributions in the porous matrix and internal pressure distributions are presented.

Curry, D. M.↗
NASA NTRS1974

Input description for Jameson's three-dimensional transonic airfoil analysis program

The input parameters are presented for a computer program which performs calculations for inviscid isentropic transonic flow over three dimensional airfoils with straight leading edges. The free stream Mach number is restricted only by the isentropic assumption. Weak shock waves are automatically located where they occur in the flow. The finite difference form of the full equation for the velocity potential is solved by the method of relaxation, after the flow exterior to the airfoil is mapped to the upper half plane.

Newman, P. A.↗
NASA NTRS1972

Hydrodynamic water impact

The hydrodynamic impact of a falling body upon a viscous incompressible fluid was investigated by numerically solving the equations of motion. Initially the mathematical model simulated the axisymmetric impact of a rigid right circular cylinder upon the initially quiescent free surface of a fluid. A compressible air layer exists between the falling cylinder and the liquid free surface. The mathematical model was developed by applying the Navier-Stokes equations to the incompressible air layer and the incompressible fluid. Assuming the flow to be one dimensional within the air layer, the average velocity, pressure and density distributions were calculated. The liquid free surface was allowed to deform as the air pressure acting on it increases. For the liquid the normalized equations were expressed in two-dimensional cylindrical coordinates. The governing equations for the air layer and the liquid were expressed in finite difference form and solved numerically. For the liquid a modified version of the Marker-and-Cell method was used. The mathematical model has been reexamined and a new approach has recently been initiated. Essentially this consists of examining the impact of an inclined plate onto a quiesent water surface with the equations now formulated in cartesian coordinates.

Kettleborough, C. F.↗
NASA NTRS1973

Derivation of a three dimensional numerical water quality model for estuary and continental shelf application

A derivation is given for a three dimensional mass transport equation which is appropriate for numerical modeling of estuary and continental shelf water quality variations for both the time dependent and steady state cases. A finite difference approximation to the derived equation is presented and a solution scheme for the resulting equations outlined. Preliminary results are obtained using the model for the extremely simple problems which have analytical solutions. The numerical model, as presented, will provide a scheme to study water quality problems in coastal waters for both steady state and time dependent cases.

Spaulding, M.↗
NASA NTRS1974

Numerical solution of the time-dependent compressible Navier-Stokes equations in inlet regions

The results of a study to determine the effects of compressibility on the viscous flow through channels that have straight, parallel walls are presented. Two channel configurations are considered, the flow between two semi-infinite flat plates with uniform flow prescribed at the inlet plane and a cascade of semi-infinite flat plates with uniform flow introduced upstream. The flow field is modeled by using the time dependent, compressible Navier-Stokes equations. Time dependent solutions are obtained by using an explicit finite difference technique which advances the pressure on near field subsonic boundaries such that accurate steady state solutions are obtained. Steady state results at Reynolds number 20 and 150 are presented for Mach numbers between 0.09 and 0.36 and compared with the incompressible solutions of previous studies.

Olson, L. E.↗