Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical computation”

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

Numerical computation of transient coaxial entry tube flows

A numerical program was developed to compute transient laminar flows in two dimensions including multicomponent mixing and chemical reaction. The program can compute both incompressible flows and compressible flows at all speeds, and it is applied to describe transient and steady state solutions for low subsonic, coaxial entry, tue flows. Single component, nonreacting flows comprise most of the solutions, but one steady state solution is presented for trace concentration constituents engaging in a second order reaction. Numerical stability was obtained by adding at each calculation point a correction for numerical diffusion errors caused by truncation of the Taylor series used to finite difference the conservation equations. Transient computations were made for fluids initially at rest, then subjected to step velocity inputs that were uniform across each region of the entry plane and were held constant throughout the computation period. For center tube to annulus velocity ratios of 0.5 and 2.0, the bulk fluid in the tube initially moved in plug flow, but strong radial flows developed near the injection plane which moved the fluid into the high shear region between the jets and away from the tube wall.

Wieber, P. R.↗
NASA NTRS1990

Numerical computation of transition to turbulence

Current applications of numerical simulations, including studies of mildly nonlinear and strongly nonlinear effects, and incipient turbulence are reviewed, and it is observed that one mildly nonlinear application for which direct simulations need to take a leading part is in establishing confidence in the dominant instability of a flow, while direct simulations are the primary tool when nonlinear effects are stronger than can be handled by existing asymptotic methods. The characteristics of existing algorithms for transition simulations are discussed, along with algorithmic improvements which would extend the range of problems amenable to numerical simulations. Current and future issues and problems are outlined.

Zang, Thomas A.↗
NASA NTRS1976

Numerical computation of transonic flows with shock waves

Some recent developments in numerical methods for calculating solutions to the transonic potential flow equation are reviewed, including (1) the construction of the stable coordinate independent difference schemes; (2) the use of conservation form to insure proper shock jump conditions; (3) analysis of the relaxation method by the time dependent analogy; (4) accelerated iterative schemes.

Jameson, A.↗
NASA NTRS1983

Numerical computation of transonic flow governed by the full-potential equation

Numerical solution techniques for solving transonic flow fields governed by the full potential equation are discussed. In a general sense relaxation schemes suitable for the numerical solution of elliptic partial differential equations are presented and discussed with emphasis on transonic flow applications. The presentation can be divided into two general categories: An introductory treatment of the basic concepts associated with the numerical solution of elliptic partial differential equations and a more advanced treatment of current procedures used to solve the full potential equation for transonic flow fields. The introductory material is presented for completeness and includes a brief introduction (Chapter 1), governing equations (Chapter 2), classical relaxation schemes (Chapter 3), and early concepts regarding transonic full potential equation algorithms (Chapter 4).

Holst, T. L.↗
NASA NTRS1977

Ray trajectories in a torus: An application of MACSYMA to a complex numerical computation

The study of ray trajectories of plasma waves in a torodial geometry using MACSYMA is an example of how symbolic, numerical, and graphical facilities can be used in concert to accomplish a complex computational goal. Computational features of this study which are of particular significance include: the derivation of code (i.e. writing functions to generate program fragments), the use of array functions to simplify the specification of a numerical iteration scheme, and the graphical presentation of the results. Mathematically, this study originates in the solution of a linear inhomogeneous partial differential equation in 3 dimensions by the method of characteristics. While it is possible to describe this equation compactly by using vector notation, and by specifying the spatial variation of the coefficients in terms of intermediate parameters, the transformation of the equation into a form amenable to solution is very tedious. A MACSYMA program is presented for obtaining description of the rf field structure excited by a waveguide located at the edge of a toroidal plasma confinement device.

Kulp, J. L.↗
NASA NTRS1975

Numerical computation of viscous flows on the lee side of blunt shapes flying at supersonic speeds

A numerical method for solving the parabolic approximation to the steady-state compressible Navier-Stokes equations is examined. The approximation neglects only the streamwise gradients of shear stress. An implicit finite difference method is used which advances the solution downstream from an initial data surface and determines the complete viscous-inviscid flow between the body and bow shock wave. It is necessary that the inviscid portion of the flow field be supersonic. Crossflow separation is determined as part of the solution. The method is applied to a 15 deg sphere-cone at 15 deg angle of attack, and the results are compared with an inviscid method-of-characteristics calculation.

Rakich, J. V.↗
NASA NTRS1983

Numerical computation of transonic flow governed by the full-potential equation

Basic concepts associated with the numerical solution of elliptic partial differential equations are introduced, and procedures used to solve the full potential equation for transonic flow fields are discussed. Governing equations, classical relaxation schemes and concepts regarding transonic, full potential equation algorithms are covered. The equation transformation and grid generation procedures; full potential spatial differencing schemes; full potential iteration schemes, emphasizing convergence acceleration; and three dimensional applications are presented.

Holst, T. L.↗
NASA NTRS1983

Numerical computation of transonic flow about wing-fuselage configurations on a vector computer

The transonic wing analysis code TWING, which uses the AF2 relaxation algorithm, has been vectorized to run on the Cray-1S computer. Vectorization of this code improved computational efficiency over that of the CDC 7600 computer by factors of 11 to 13. The improvement compares favorably with the prediction of a theoretical performance model. A convenient generalization now permits the treatment of rudimentary wing-fuselage combinations. Flow predictions for a transport configuration in both isolated-wing and wing-fuselage modes show the expected trends in shock strength and position when compared with wind-tunnel results. An isolated fighter wing is examined in terms of execution time on three different computers and in comparison with experimental data. The computational fluid dynamics code produced during this study is a careful union of an efficient three-dimensional, transonic, numerical algorithm and the vector features presently available on modern computers.

Thomas, S. D.↗
NASA NTRS1989

Numerical computation of shock wave-turbulent boundary layer interaction in transonic flow over an axisymmetric curved hill

A control-volume based finite difference computation of a turbulent transonic flow over an axisymmetric curved hill is presented. The numerical method is based on the SIMPLE algorithm, and hence the conservation of mass equation is replaced by a pressure correction equation for compressible flows. The turbulence is described by a k-epsilon turbulence model supplemented by a near-wall turbulence model. In the method, the dissipation rate in the region very close to the wall is obtained from an algebraic equation and that for the rest of the flow domain is obtained by solving a partial differential equation for the dissipation rate. The other flow equations are integrated up to the wall. It is shown that the present turbulence model yields the correct location of the compression shock. The other computational results are also in good agreement with experimental data.

Kim, S.-W.↗
NASA NTRS1994

Numerical Computation of Flame Spread over a Thin Solid in Forced Concurrent Flow with Gas-phase Radiation

Excerpts from a paper describing the numerical examination of concurrent-flow flame spread over a thin solid in purely forced flow with gas-phase radiation are presented. The computational model solves the two-dimensional, elliptic, steady, and laminar conservation equations for mass, momentum, energy, and chemical species. Gas-phase combustion is modeled via a one-step, second order finite rate Arrhenius reaction. Gas-phase radiation considering gray non-scattering medium is solved by a S-N discrete ordinates method. A simplified solid phase treatment assumes a zeroth order pyrolysis relation and includes radiative interaction between the surface and the gas phase.

Jiang, Ching-Biau↗
NASA NTRS1975

Numerical computation of three-dimensional blunt body flow fields with an impinging shock

A time-marching finite-difference method was used to solve the compressible Navier-Stokes equations for the three-dimensional wing-leading-edge shock impingement problem. The bow shock was treated as a discontinuity across which the exact shock jump conditions were applied. All interior shock layer detail such as shear layers, shock waves, jets, and the wall boundary layer were automatically captured in the solution. The impinging shock was introduced by discontinuously changing the freestream conditions across the intersection line at the bow shock. A special storage-saving procedure for sweeping through the finite-difference mesh was developed which reduces the required amount of computer storage by at least a factor of two without sacrificing the execution time. Numerical results are presented for infinite cylinder blunt body cases as well as the three-dimensional shock impingement case. The numerical results are compared with existing experimental and theoretical results.

Holst, T. L.↗
NASA NTRS1992

Numerical computation of steady-state acoustic disturbances in flow

Two time domain methods for computing two dimensional steady-state acoustic disturbances propagating through internal subsonic viscous flow fields in the presence of variable area are investigated. The first method solves the Navier-Stokes equations for the combined steady and acoustic field together and subtracts the steady flow to obtain the acoustic field. The second method solves a system of perturbation equations to obtain the acoustic disturbances, making use of a separate steady flow computation as input to the system. In each case the periodic steady-state acoustic fluctuations are obtained numerically on a supercomputer using a second order unsplit explicit MacCormack predictor-corrector method. Results show that the first method is not very effective for computing acoustic disturbances of even moderate amplitude. It appears that more accurate steady flow algorithms are required for this method to succeed. On the other hand, linear and nonlinear acoustic disturbances extracted from the perturbation approach are shown to exhibit expected behavior for the problems considered. It is also found that inflow boundary conditions for an equivalent uniform duct can be successfully applied to a nonuniform duct to obtain steady-state acoustic disturbances.

Watson, W. R.↗
NASA NTRS1990

On the numerical computation of nonlinear force-free magnetic fields

An algorithm has been developed to extrapolate nonlinear force-free magnetic fields from the photosphere, given the proper boundary conditions. This paper presents the results of this work, describing the mathematical formalism that was developed, the numerical techniques employed, and comments on the stability criteria and accuracy developed for these numerical schemes. An analytical solution is used for a benchmark test; the results show that the computational accuracy for the case of a nonlinear force-free magnetic field was on the order of a few percent (less than 5 percent). This newly developed scheme was applied to analyze a solar vector magnetogram, and the results were compared with the results deduced from the classical potential field method. The comparison shows that additional physical features of the vector magnetogram were revealed in the nonlinear force-free case.

Wu, S. T.↗
NASA NTRS1985

On the numerical computation of nonlinear force-free magnetic fields

An algorithm has been developed to extrapolate nonlinear force-free magnetic fields from a source surface, given the proper boundary conditions. The results of this work; describing the mathematical formalism that was developed, the numerical techniques employed, and the stability criteria developed for these numerical schemes are presented. An analytical solution is used for a test case; the results show that the computational accuracy for the case of a nonlinear force-free magnetic field was on the order of a few percent ( 5%).

Wu, S. T.↗
NASA NTRS1978

On the use of co-ordinate stretching in the numeral computation of high frequency scattering

The scattering of the sound of a jet engine by an airplane fuselage is modeled by solving the axially symmetric Helmholtz equation exterior to a long thin ellipsoid. The integral equation method based on the single layer potential formulation is used. A family of coordinate systems on the body is introduced and an algorithm is presented to determine the optimal coordinate system. Numerical results verify that the optimal choice enables the solution to be computed with a grid that is coarse relative to the wavelength.

Bayliss, A.↗
NASA NTRS1976

Laser beam scintillation beyond the turbulent atmosphere A numerical computation

The extended Huygens-Fresnel formulation for propagation through turbulence is used to examine scintillation of a finite laser beam. The method is demonstrated analytically for propagation beyond a weak Gaussian phase screen. A numerical integration technique is used to extend the results to a more realistic turbulence model. Results are compared with existing Gaussian beam propagation theory.

Bufton, J. L.↗
NASA NTRS1996

Active Control of Fan Noise: Feasibility Study: Numerical Computation of Acoustic Mode Reflection Coefficients for an Unflanged Cylindrical Duct - Volume 5

A computational method to predict modal reflection coefficients in cylindrical ducts has been developed based on the work of Homicz, Lordi, and Rehm, which uses the Wiener-Hopf method to account for the boundary conditions at the termination of a thin cylindrical pipe. The purpose of this study is to develop a computational routine to predict the reflection coefficients of higher order acoustic modes impinging on the unflanged termination of a cylindrical duct. This effort was conducted wider Task Order 5 of the NASA Lewis LET Program, Active Noise Control of aircraft Engines: Feasibility Study, and will be used as part of the development of an integrated source noise, acoustic propagation, ANC actuator coupling, and control system algorithm simulation. The reflection coefficient prediction will be incorporated into an existing cylindrical duct modal analysis to account for the reflection of modes from the duct termination. This will provide a more accurate, rapid computation design tool for evaluating the effect of reflected waves on active noise control systems mounted in the duct, as well as providing a tool for the design of acoustic treatment in inlet ducts. As an active noise control system design tool, the method can be used preliminary to more accurate but more numerically intensive acoustic propagation models such as finite element methods. The resulting computer program has been shown to give reasonable results, some examples of which are presented. Reliable data to use for comparison is scarce, so complete checkout is difficult, and further checkout is needed over a wider range of system parameters. In future efforts the method will be adapted as a subroutine to the GEAE segmented cylindrical duct modal analysis program.

Kraft, R. E.↗
NASA NTRS1984

On accuracy conditions for the numerical computation of waves

The Helmholtz equation (Delta + K(2)n(2))u = f with a variable index of refraction n, and a suitable radiation condition at infinity serves as a model for a wide variety of wave propagation problems. Such problems can be solved numerically by first truncating the given unbounded domain and imposing a suitable outgoing radiation condition on an artificial boundary and then solving the resulting problem on the bounded domain by direct discretization (for example, using a finite element method). In practical applications, the mesh size h and the wave number K, are not independent but are constrained by the accuracy of the desired computation. It will be shown that the number of points per wavelength, measured by (Kh)(-1), is not sufficient to determine the accuracy of a given discretization. For example, the quantity K(3)h(2) is shown to determine the accuracy in the L(2) norm for a second-order discretization method applied to several propagation models.

Bayliss, A.↗