Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical solution”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Numerical solutions of the triple-deck equations for laminar trailing-edge stall

The problem of determining the effect of laminar boundary layers on the lift of thin wings in subsonic flow at high Reynolds numbers is considered. The boundary value problem is formulated in the framework of the triple-deck theory of Brown and Stewartson. The resulting fourth-order boundary value was solved by an iterative finite-difference technique. An inverse iteration procedure provides proper treatment of the trailing-edge singularity, and asymptotic far-field expansions and coordinate stretchings are used to deal with the problem of the slow algebraic decay of the solution.

Chow, R.↗
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 NTRS1978

A numerical solution of the Navier-Stokes equations for chemically nonequilibrium, merged stagnation shock layers on spheres and two-dimensional cylinders in air

Results of solving the Navier-Stokes equations for chemically nonequilibrium, merged stagnation shock layers on spheres and two-dimensional cylinders are presented. The effects of wall catalysis and slip are also examined. The thin shock layer assumption is not made, and the thick viscous shock is allowed to develop within the computational domain. The results show good comparison with existing data. Due to the more pronounced merging of shock layer and boundary layer for the sphere, the heating rates for spheres become higher than those for cylinders as the altitude is increased.

Johnston, K. D.↗
NASA NTRS1978

The numerical solution of viscous flows at high Reynolds number

A review is presented of implicit and hybrid methods applicable to solving viscous flows at high Reynolds numbers. Flows within axisymmetric channels containing stationary shock waves, past blunt-nosed lifting airfoils, past sharp-nosed symmetric airfoils with buffet, past three-dimensional compression ramps with side walls, and past ogive- and hemisphere-cylinders at angle of attack have been examined. Reynolds numbers as high as 10 to the 9th power have been used.

Maccormack, R. W.↗
NASA NTRS1979

Numerical solution for perpendicular sonic hydrogen injection into a ducted supersonic airstream

This note discusses a computer program being developed to study the flow field near opposing perpendicular fuel injectors in scramjets. The MacCormack time-split, finite difference relaxation technique was used to solve the full two-dimensional compressible Navier-Stokes equations along with energy and species equations. By using this technique, a program was developed to consider the turbulent nonreacting flow of hydrogen and air in a rectangular duct. A damping term, proportional to the second derivative of pressure and temperature, was used to produce a stable solution behind the hydrogen jet in the neighborhood of the recompression shock. A case using actual conditions encountered in current scramjet design was analyzed, with results agreeing qualitatively with experimental observations.

Drummond, J. P.↗
NASA NTRS1979

Artificial compressibility methods for numerical solutions of transonic full potential equation

New methods for transonic flow computations based on the full potential equation in conservation form are presented. The idea is to modify slightly the density (due to the artificial viscosity in the supersonic region), and solve the resulting elliptic-like problem iteratively. It is shown that standard discretization techniques (central differencing) as well as some standard iterative procedures (SOR, ADI, and explicit methods) are applicable to the modified transonic mixed-type equation. Calculations of transonic flows around cylinders and airfoils are discussed with special emphasis on the explicit methods that are suitable for vector processing on the STAR 100 computer.

Hafez, M.↗
NASA NTRS1979

Numerical solution of supersonic laminar flow over an inclined body of revolution

A mixed explicit-implicit scheme is used to solve the time-dependent thin-layer approximation of the Navier-Stokes equations for a supersonic laminar flow over an inclined body of revolution. Test cases for Mach 2.8 flow over a cylinder with 15 deg flare angle at angles of attack of 0, 1, and 4 deg are calculated. Good agreement is obtained between the present computed results and experimental measurements of surface pressure. A pair of vortices on the leeward and a peak in the normal force distribution near the flared juncture are predicted; the role of circumferential communication is discussed.

Hung, C. M.↗
NASA NTRS1980

Boundary conditions for the numerical solution of elliptic equations in exterior regions

Elliptic equations in exterior regions frequently require a boundary condition at infinity to ensure the well-posedness of the problem. Examples of practical applications include the Helmholtz equation and Laplace's equation. Computational procedures based on a direct discretization of the elliptic problem require the replacement of the condition on a finite artificial surface. Direct imposition of the condition at infinity along the finite boundary results in large errors. A sequence of boundary conditions is developed which provides increasingly accurate approximations to the problem in the infinite domain. Estimates of the error due to the finite boundary are obtained for several cases. Computations are presented which demonstrate the increased accuracy that can be obtained by the use of the higher order boundary conditions. The examples are based on a finite element formulation but finite difference methods can also be used.

Bayliss, A.↗
NASA NTRS1980

Stellar convection. II - A multimode numerical solution for convection in spheres

A stable, equilibrium solution for convection in a self-gravitating sphere of Boussinesq fluid was computed by using a modal analysis in which the theta, phi dependence of the fluid is expanded in a set of 168 spherical harmonics. The temperature, velocity, and convective flux of the fluid as well as the kinetic and thermal energy spectra as functions of wavelength are computed. The spectra are found to be in agreement with both experimental observations and analytic scaling laws. The dynamics of the energy cascade is examined by computing the ratio of the amount of energy dissipation at a particular wavelength to the amount of energy produced at that same wavelength. It is found that there is only a slight cascade of kinetic energy to smaller wavelengths but a large cascade of thermal energy.

Marcus, P. S.↗
NASA NTRS1982

Numerical Solution for Navier-Stokes Equations

Carefully selected blend of computational techniques solves complete set of equations for viscous, unsteady, hypersonic flow in general curvilinear coordinates. New algorithm has tested computation of axially directed flow about blunt body having shape similar to that of such practical bodies as wide-body aircraft or artillery shells. Method offers significant computational advantages because of conservation-law form of equations and because it reduces amount of metric data required.

Warsi, Z. U. A.↗
NASA NTRS1980

Nonlinear grid error effects on numerical solution of partial differential equations

Finite difference solutions of nonlinear partial differential equations require discretizations and consequently grid errors are generated. These errors strongly affect stability and convergence properties of difference models. Previously such errors were analyzed by linearizing the difference equations for solutions. Properties of mappings of decadence were used to analyze nonlinear instabilities. Such an analysis is directly affected by initial/boundary conditions. An algorithm was developed, applied to nonlinear Burgers equations, and verified computationally. A preliminary test shows that Navier-Stokes equations may be treated similarly.

Dey, S. K.↗
NASA NTRS1980

Numerical solutions of transonic flows by parametric differentiation and integral equation techniques

The paper discusses results of an exploratory study of the advantages obtained by combining integral equation methods with the method of parametric differentiation in the treatment of transonic flow problems. In the proposed method, the nonlinear unsteady transonic flow equation for small perturbations is transformed into a linear equation by the use of the method of parametric differentiation. The linear equation is split into a pair of weakly coupled partial differential equations by writing the transformed perturbation potential as the sum of a steady component and an unsteady component. The solution of the steady equation as an integral equation is based on Ogana's treatment (1978). As a test case, the formulation is applied to predict the steady transonic flow over a nonlifting parabolic-arc airfoil.

Sivaneri, N. T.↗
NASA NTRS1981

Numerical solution of supersonic viscous flow over blunt delta wings

A general parabolized Navier-Stokes code has been developed to compute the steady supersonic viscous flow around arbitrary body shapes at high angles of attack. A nonorthogonal three-dimensional coordinate frame permits the code to march with solution surfaces which are the most appropriate. The code has been used to calculate the laminar flow over a slab delta wing with 70 deg sweep at angles of attack up to 41.5 deg and Mach numbers of 6.8 and 9.6. The computed shock shapes, surface pressures and heat transfer coefficients are compared with experiment and show agreement.

Tannehill, J. C.↗
NASA NTRS1981

Numerical solution of transonic stream function equation

The stream function equation, in conservation form, looks similar to the full potential equation and existing methods (e.g. artificial compressibility) can be readily applied. Rotational flows can be calculated once the vorticity (due to shocks or nonuniformity) is evaluated. There are, however, two main difficulties: First, the density is not uniquely determined in terms of the flux (there are two solutions; the subsonic and the supersonic branch with a square root singularity at the sonic point). Methods to overcome this difficulty are studied and results are presented with some remarks on inviscid separation and closed stream lines. Second, the need of two stream functions for three dimensional calculations is briefly discussed.

Hafez, M.↗
NASA NTRS1982

Numerical solution of three-dimensional unsteady transonic flow over wings including inviscid/viscous interactions

A method is presented for computing the unsteady three-dimensional transonic flowfield about twin wings of arbitrary planform and section geometry including the effects of inviscid/viscous interaction. Dominant effects of the shock-boundary-layer interaction are accounted for by a simple empirically defined model, while viscous regions adjacent to the wing surface and along the trailing wake are described by a set of integral equations appropriate for two-dimensional turbulent flows. The validity of the proposed method is verified by comparison with steady three-dimensional experimental data and a solution of more exact equations for a typical transport wing configuration.

Rizzetta, D. P.↗
NASA NTRS1982

Numerical solution of Space Shuttle Orbiter flow field

The supersonic, viscous laminar flow around the Space Shuttle Orbiter forebody has been computed with a parabolized Navier-Stokes code using a generalized coordinate transformation. The initial solution for the nose part of the Orbiter geometry was obtained with a three-dimensional time-dependent Navier-Stokes solver. It was necessary to employ a wind axis oriented coordinate system to obtain the initial solution with the time-dependent code. The generalized PNS technique used in this study allows the solution to be marched from the given initial data surface to any desired surface downstream. A grid point clustering scheme was employed to accurately describe the body shape by clustering points at the wing tip and at the wing body juncture. The computed heat transfer coefficients, pressure coefficients, and shock shapes are compared with the available experimental data for 0 degrees and 30 degrees angle of attack.

Venkatapathy, E.↗