SEARCH · Engineering Papers

Results for “Finite-differences”

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 361 records · Page 20

All sources NASA NTRS DOE OSTI

NASA NTRS1981

Computation of simple three-dimensional wing-vortex interaction in transonic flow

An implicit delta form finite-difference algorithm for Euler equations in conservation law form has been used in preliminary calculations of three-dimensional wing-vortex interactions. Both steady and unsteady transonic flow wing-vortex interactions are computed. The computations themselves are meant to guide upcoming wind tunnel experiments of the same flow field. Various modifications to the numerical method that are intended to improve computational efficiency are also described and tested in both two- and three-dimensions.

Srinivasan, G. R.↗

NASA NTRS1981

A numerical simulation of two-dimensional inlet flow fields

Inlet flow fields for airbreathing missiles are calculated by the adaptation of a two-dimensional computational method developed for the flow around airfoils. A supersonic free stream is assumed to allow the forebody calculation to be uncoupled from the inlet calculation. The inlet calculation employs an implicit, time-marching finite-difference procedure to solve the thin-layer Navier-Stokes equations formulated in body-fitted coordinates. Because the method can be used for a flow field with both subsonic and supersonic regions, it is applicable to subcritical as well as supercritical inlet operation.

Biringen, S.↗

NASA NTRS1981

A split explicit reformulation of the regional numerical weather prediction model of the Japan Meteorological Agency

The split explicit integration scheme for numerical weather prediction models is employed in a version of the regional numerical weather prediction model of the Japan Meteorological Agency. The finite-difference scheme of the model is designed in the manner proposed by Okamura (1975). The horizontal advection terms in the governing equations are integrated with a time step limited by the wind speed while the terms which describe inertial-gravity oscillations are integrated in a succession of shorter time steps. The physical processes included within the model are precipitation, small-scale convection, surface exchanges of sensible and latent heat, and radiative heating and cooling. An example of a surface pressure forecast over Europe is shown for initial data observed at 0000 GMT 29 December 1979. Quantitative precipitation forecasts over Europe and North America for the 24 h period beginning at 0000 GMT 30 December 1979 are also shown. It is concluded that the model is capable of realistically depicting the evolution of synoptic-scale systems.

Duffy, D. G.↗

NASA NTRS1981

A modified Dodge algorithm for the parabolized Navier-Stokes equation and compressible duct flows

A revised version of Dodge's split-velocity method for numerical calculation of compressible duct flow was developed. The revision incorporates balancing of mass flow rates on each marching step in order to maintain front-to-back continuity during the calculation. The (checkerboard) zebra algorithm is applied to solution of the three dimensional continuity equation in conservative form. A second-order A-stable linear multistep method is employed in effecting a marching solution of the parabolized momentum equations. A checkerboard iteration is used to solve the resulting implicit nonlinear systems of finite-difference equations which govern stepwise transition. Qualitive agreement with analytical predictions and experimental results was obtained for some flows with well-known solutions.

Cooke, C. H.↗

NASA NTRS1982

Guidelines for developing vectorizable computer programs

Some fundamental principles for developing computer programs which are compatible with array-oriented computers are presented. The emphasis is on basic techniques for structuring computer codes which are applicable in FORTRAN and do not require a special programming language or exact a significant penalty on a scalar computer. Researchers who are using numerical techniques to solve problems in engineering can apply these basic principles and thus develop transportable computer programs (in FORTRAN) which contain much vectorizable code. The vector architecture of the ASC is discussed so that the requirements of array processing can be better appreciated. The "vectorization" of a finite-difference viscous shock-layer code is used as an example to illustrate the benefits and some of the difficulties involved. Increases in computing speed with vectorization are illustrated with results from the viscous shock-layer code and from a finite-element shock tube code. The applicability of these principles was substantiated through running programs on other computers with array-associated computing characteristics, such as the Hewlett-Packard (H-P) 1000-F.

Miner, E. W.↗

NASA NTRS1980

A vectorized code for calculating laminar and turbulent hypersonic flows about blunt axisymmetric bodies at zero and small angles of attack

A user's guide is provided for a computer code which calculates the laminar and turbulent hypersonic flows about blunt axisymmetric bodies, such as spherically blunted cones, hyperboloids, etc., at zero and small angles of attack. The code is written in STAR FORTRAN language for the CDC-STAR-100 computer. Time-dependent, viscous-shock-layer-type equations are used to describe the flow field. These equations are solved by an explicit, two-step, time asymptotic, finite-difference method. For the turbulent flow, a two-layer, eddy-viscosity model is used. The code provides complete flow-field properties including shock location, surface pressure distribution, surface heating rates, and skin-friction coefficients. This report contains descriptions of the input and output, the listing of the program, and a sample flow-field solution.

Kumar, A.↗

NASA NTRS1980

User's guide for the computer code COLTS for calculating the coupled laminar and turbulent flow over a Jovian entry probe

A user's guide for a computer code 'COLTS' (Coupled Laminar and Turbulent Solutions) is provided which calculates the laminar and turbulent hypersonic flows with radiation and coupled ablation injection past a Jovian entry probe. Time-dependent viscous-shock-layer equations are used to describe the flow field. These equations are solved by an explicit, two-step, time-asymptotic finite-difference method. Eddy viscosity in the turbulent flow is approximated by a two-layer model. In all, 19 chemical species are used to describe the injection of carbon-phenolic ablator in the hydrogen-helium gas mixture. The equilibrium composition of the mixture is determined by a free-energy minimization technique. A detailed frequency dependence of the absorption coefficient for various species is considered to obtain the radiative flux. The code is written for a CDC-CYBER-203 computer and is capable of providing solutions for ablated probe shapes also.

Kumar, A.↗

NASA NTRS1981

Computation of airfoil buffet boundaries

The ILLIAC IV computer has been programmed with an implicit, finite-difference code for solving the thin layer compressible Navier-Stokes equation. Results presented for the case of the buffet boundaries of a conventional and a supercritical airfoil section at high Reynolds numbers are found to be in agreement with experimentally determined buffet boundaries, especially at the higher freestream Mach numbers and lower lift coefficients where the onset of unsteady flows is associated with shock wave-induced boundary layer separation.

Levy, L. L., Jr.↗

NASA NTRS1981

Flow aerodynamics modeling of an MHD swirl combustor - Calculations and experimental verification

The paper describes a computer code for calculating the flow dynamics of a constant-density flow in the second-stage trumpet shaped nozzle section of a two-stage MHD swirl combustor for application to a disk generator. The primitive pressure-velocity variable, finite-difference computer code has been developed for the computation of inert nonreacting turbulent swirling flows in an axisymmetric MHD model swirl combustor. The method and program involve a staggered grid system for axial and radial velocities, and a line relaxation technique for the efficient solution of the equations. The code produces as output the flow field map of the nondimensional stream function, axial and swirl velocity. It was found that the best location for seed injection to obtain a uniform distribution at the combustor exit is in the central location for seed injected at the entrance to the second stage combustor.

Gupta, A. K.↗

NASA NTRS1981

Test problems for inviscid transonic flow

The paper discusses some results obtained in the solving of test problems for inviscid transonic flow with shock waves using the TRANDES program. The method used employs the full inviscid perturbation-potential flow equation in a Cartesian grid system that is stretched to infinity. The equation is represented by a nonconservative system of finite-difference equations that includes at supersonic points a rotated difference scheme and is solved by column relaxation. Except for the NACA 0012 case with the freestream Mach number equal to 0.95, all test problems were solved straightforwardly and appeared to be converged or close to convergence. The only difficulty was some sensitivity to grid placement, which is typical of the Cartesian formulation

Carlson, L. A.↗

NASA NTRS1982

A numerical solution of the supersonic flow over a rearward facing step with transverse non-reacting hydrogen injection

The present investigation represents an application of computational fluid dynamics to a problem associated with the flow in the combustor region of a supersonic combustion ramjet engine (scramjet). The governing equations are considered, taking into account the Navier-Stokes equations, a molecular viscosity calculation, the molecular thermal conductivity, molecular diffusion, and a turbulence model. The employed numerical solution is patterned after the explicit, time-dependent, unsplit, predictor-corrector, finite-difference method given by MacCormack (1969). The calculation is concerned with the supersonic flow over a rearward-facing step with transverse H2 injection at conditions germane to the combustor region of a scramjet engine. The H2 jet acts as an effective body which essentially shields the primary flow from the rearward-facing step, thus substantially changing the wave pattern in the primary flow.

Berman, H. A.↗

NASA NTRS1982

Application of two-point difference schemes to the conservative Euler equations for one-dimensional flows

An implicit finite-difference method is presented for obtaining steady-state solutions to the time-dependent, conservative Euler equations for flows containing shocks. The method uses a two-point central difference scheme with dissipation added at supersonic points via the retarded density concept. Application of the method to the one-dimensional nozzle flow equations for various combinations of subsonic and supersonic boundary conditions show the method to be very efficient. Residuals are typically reduced to machine zero in approximately 35 time steps for 50 mesh points. It is shown that the scheme offers certain advantages over the more widely-used three-point schemes, especially in regard to application of boundary conditions.

Wornom, S. F.↗

NASA NTRS1981

Use of interactive graphics to analyze QUICK-geometry

The QUICK InterActive Graphics Analysis (QUIAGA) program and its advantages for displaying aircraft QUICK-geometry to aid in detection and analysis of errors are described. The QUICK-geometry system generates a completely analytical aircraft geometry description for use by finite-difference flow codes. The QUIAGA program was developed to exercise the QUICK-geometry subroutines to examine the analytic definition of a configuration by plotting cross sections and body lines on a graphics terminal. A number of options are available, including multiple cross-section views, hidden-line removal, and display of control point locations. Use of these options for the detection and analysis of errors in the QUICK-geometry definition can be of great assistance in speedily arriving at a correct analytical geometry description for flow-field computation. The QUIAGA program has been used in developing a QUICK-geometry model of the NASA Space Shuttle Orbiter, and examples from this experience are given to show some of the program's features. Details of program usage and an example session are given in the appendixes.

Townsend, J. C.↗

NASA NTRS1981

Mesh generation by conformal and quasiconformal mappings

It is pointed out that many recent advances in the finite-difference solution of elliptic equations have been limited to regions whose boundary contours coincide with coordinate lines of the Cartesian coordinate system. The reason for this is related to the fact that in the case of an arbitrary curvilinear coordinate system the original equation becomes much more complex. However, there is no added complexity if an orthogonal coordinate system is generated from a conformal mapping. In the present investigation, a finite difference method developed for the construction of conformal mappings has been generalized to construct quasi-conformal mappings. It is expected that the use of more sophisticated numerical algorithms could lead to improvements in both speed and accuracy. Quasi-conformal mappings have applications not only in the solution of elliptic equations but also in other areas such as orthogonal mesh generation on surfaces and the solution of certain fluid flow problems.

Mastin, C. W.↗

NASA NTRS1981

Viscous-inviscid interaction method including wake effects for three-dimensional wing-body configurations

A viscous-inviscid interaction method has been developed by using a three-dimensional integral boundary-layer method which produces results in good agreement with a finite-difference method in a fraction of the computer time. The integral method is stable and robust and incorporates a model for computation in a small region of streamwise separation. A locally two-dimensional wake model, accounting for thickness and curvature effects, is also included in the interaction procedure. Computation time spent in converging an interacted result is, many times, only slightly greater than that required to converge an inviscid calculation. Results are shown from the interaction method, run at experimental angle of attack, Reynolds number, and Mach number, on a wing-body test case for which viscous effects are large. Agreement with experiment is good; in particular, the present wake model improves prediction of the spanwise lift distribution and lower surface cove pressure.

Streett, C. L.↗

NASA NTRS1982

Design principles for high efficiency small-grain polysilicon solar cells, with supporting experimental studies

Design principles suggested here aim toward high conversion efficiency (greater than 15 percent) in polysilicon cells. The principles seek to decrease the liabilities of both intragranular and grain-boundary-surface defects. The advantages of a phosphorus atom concentration gradient in a thin (less than 50 microns) base of a p(+)/n(x)/n(+) drift-field solar cell, which produces favorable gradients in chemical potential, minority-carrier mobility and diffusivity, and recombination lifetime (via phosphorus gettering) are suggested. The degrading effects of grain boundaries are reduced by these three gradients and by substituting atoms (P, H, F or Li) for vacancies on the grain-boundary surface. From recent experiments comes support for the benefits of P diffusion down grain boundaries and, for quasi-grain-boundary-free and related structures. New analytic solutions for the n(x)-base include the effect of a power-law dependence between P concentration and lifetime. These provide an upper-bound estimate on the open circuit voltage. Finite-difference numerical solutions of the six Shockley equations furnish complete information about all solar-cell parameters and add insight concerning design.

Lindholm, F. A.↗

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↗