Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical schemes”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

An iterative radiative transfer code for ocean-atmosphere systems

The details of an iterative radiative transfer code for computing the intensity and degree of polarization of diffuse radiation in models of the ocean-atmosphere system are described. The present code neglects the upwelling radiation from below the ocean surface and as such can be applied to the part of the spectrum where the absorption by water is strong. To establish the reliability of the numerical scheme and the computer code, the results are compared with those of Fraser and Walker (1968), Dave (1972), and Mullamaa (1964); they are found to be in excellent agreement. The computations also show that both the intensity and the degree of polarization of the upwelling diffuse radiation at the top of the atmosphere vary significantly when the rough ocean at the base of the atmosphere is replaced by a Lambertian surface that reflects the same energy as the rough ocean.

Ahmad, Z.↗
NASA NTRS1981

Numerical analysis of turbulent coaxial flow with internal heat generation

A computational method with which to obtain a physical understanding of the turbulent field of two coaxial jets entering an axisymmetric chamber is developed. Even the laminar field of this flow is quite complicated. This is due to the many different domains which exist in the field especially in the entrance region. Physically, three regions may be identified: the wall region, the initial region near the axis of symmetry and the mixing region. Advancing downstream, these regions change relative size with the ratio of the two jets' mass fluxes as the main parameter. The turbulent field of these flows is much more complicated due to the difference in the effective transport coefficients and turbulence level from region to region. However, being aware beforehand of the complications and the different regions of this field, the appropriate turbulence model and numerical scheme can be adjusted to treat the problem.

Lin, A.↗
NASA NTRS1984

Spectral methods in time for hyperbolic equations

A pseudospectral numerical scheme for solving linear, periodic, hyperbolic problems is described. It has infinite accuracy both in time and in space. The high accuracy in time is achieved without increasing the computational work and memory space which is needed for a regular, one step explicit scheme. The algorithm is shown to be optimal in the sense that among all the explicit algorithms of a certain class it requires the least amount of work to achieve a certain given resolution. The class of algorithms referred to consists of all explicit schemes which may be represented as a polynomial in the spatial operator.

Tal-Ezer, H.↗
NASA NTRS1984

Contact stresses in pin-loaded orthotropic plates

The effects of pin elasticity, friction, and clearance on the stresses near the hole in a pin-loaded orthotropic plate are described. The problem is modeled as a contact elasticity problem using complex variable theory, the pin and the plate being two elastic bodies interacting through contact. This modeling is in contrast to previous works which assumed that the pin is rigid or that it exerts a known cosinusoidal radial traction on the hole boundary. Neither of these approaches explicitly involves a pin. A collocation procedure and iteration were used to obtain numerical results for a variety of plate and pin elastic properties and various levels of friction and clearance. Collocation was used to enforce the boundary and iteration was used to find the contact and no-slip regions on the boundary. Details of the numerical scheme are discussed.

Hyer, M. W.↗
NASA NTRS1984

A numerical study of the controlled flow tunnel for a high lift model

A controlled flow tunnel employs active control of flow through the walls of the wind tunnel so that the model is in approximately free air conditions during the test. This improves the wind tunnel test environment, enhancing the validity of the experimentally obtained test data. This concept is applied to a three dimensional jet flapped wing with full span jet flap. It is shown that a special treatment is required for the high energy wake associated with this and other V/STOL models. An iterative numerical scheme is developed to describe the working of an actual controlled flow tunnel and comparisons are shown with other available results. It is shown that control need be exerted over only part of the tunnel walls to closely approximate free air flow conditions. It is concluded that such a tunnel is able to produce a nearly interference free test environment even with a high lift model in the tunnel.

Parikh, P. C.↗
NASA NTRS1980

Solution of 3-dimensional time-dependent viscous flows. Part 2: Development of the computer code

There is considerable interest in developing a numerical scheme for solving the time dependent viscous compressible three dimensional flow equations to aid in the design of helicopter rotors. The development of a computer code to solve a three dimensional unsteady approximate form of the Navier-Stokes equations employing a linearized block emplicit technique in conjunction with a QR operator scheme is described. Results of calculations of several Cartesian test cases are presented. The computer code can be applied to more complex flow fields such as these encountered on rotating airfoils.

Weinberg, B. C.↗
NASA NTRS1984

Legendre-tau approximation for functional differential equations. Part 2: The linear quadratic optimal control problem

The numerical scheme based on the Legendre-tau approximation is proposed to approximate the feedback solution to the linear quadratic optimal control problem for hereditary differential systems. The convergence property is established using Trotter ideas. The method yields very good approximations at low orders and provides an approximation technique for computing closed-loop eigenvalues of the feedback system. A comparison with existing methods (based on averaging and spline approximations) is made.

Ito, K.↗
NASA NTRS1984

Boundary layer transition effects on flow separation around V/STOL engine inlets at high incidence

Numerical methods for calculating laminar and turbulent boundary layer development around vertical-short take off and landing engine inlets at high incidence angles are investigated. Various transition models were compared and evaluated in calculations off flow separation bound inside the inlet. Results of the transition effects on the boundary layer characteristics at onset of separation for two types of engine inlet geometries are presented. Some of the numerical results are compared with existing wind-tunnel test data for scaled inlet models to demonstrate the effects of transition models in the numerical scheme. The effects of transition modeling on the boundary layer development are illustrated for typical engine operating conditions.

Chou, D. C.↗
NASA NTRS1984

A two-dimensional analysis of laser heat addition in converging nozzles

The two-dimensional equations of motion describing the interaction between a laser beam and a flowing gas are considered. An implicit numerical scheme is used to solve these equations for unchoked flow through a converging-diverging nozzle. Separate grids are used for the fluid dynamics and the radiation equations. The effects of beam focusing and cross-beam intensity profiles are included. The calculations are based upon real gas properties for all quantities except the gas absorptivity, which is taken as a constant. The solutions contain the expected hot central core region with cool gas near the walls. This results in steep temperature gradients in both the streamwise and cross-stream directions. The absorption zone acts as a blockage in the nozzle causing a nonuniform velocity profile at the inlet and an overall decrease in mass flow. The absorption region also forces the streamlines to move away from the axis of symmetry, although this effect is not strong.

Molvik, G. A.↗
NASA NTRS1983

An eigenvalue method for solving transient heat conduction problems

The eigenvalue method, which has been used by researchers in structure mechanics, is applied to problems in heat conduction. Its formulation is decribed in terms of an examination of transient heat conduction in a square slab. Taking advantage of the availability of the exact solution, we compare the accuracy and other numerical properties of the eigenvalue method with those of existing numerical schemes. The comparsion shows that, overall, the eigenvalue method appears to be fairly attractive. Furthermore, only a few dominant eigenvalues and their corresponding eigenvectors need to be computed and retained to yield reasonably high accuracy. Greater savings are attained in the computation time for a transient problem with long time duration and a large computational domain.

Shih, T. M.↗
NASA NTRS1982

Numerical solution of viscous flow in unbounded fluid

Incompressible viscous flow fields induced by initial vorticity distributions with bounded support or exponential decay in the far field are investigated. A numerical scheme for the solution of the vorticity distribution and the velocity field is presented with special emphasis on the treatment of the boundary data. The efficiency of the scheme is demonstrated. The present method has been applied to the study of the merging and collision of vortex rings.

Liu, C. H.↗
NASA NTRS1984

Dynamic response of shock waves in transonic diffuser and supersonic inlet - An analysis with the Navier-Stokes equations and adaptive grid

An existing method which solves the multi-dimensional ensemble-averaged compressible time-dependent Navier-Stokes equations in conjunction with mixing length turbulence model and shock capturing technique has been extended to include the shock-tracking adaptive grid systems. The numerical scheme for solving the governing equations is based on a linearized block implicit approach. The effects of grid-motion and grid-distribution on the calculated flow solutions have been studied in relative detail and this is carried out in the context of physically steady, shocked flows computed with non-stationary grids. Subsequently, the unsteady dynamics of the flows occurring in a supercritically operated transonic diffuser and a mixed compression supersonic inlet have been investigated with the adaptive grid systems by solving the Navier-Stokes equations.

Liu, N.-S.↗
NASA NTRS1984

Thrust chamber performance using Navier-Stokes solution

The viscous, axisymmetric flow in the thrust chamber of the space shuttle main engine (SSME) was computed on the CRAY 205 computer using the general interpolants method (GIM) code. Results show that the Navier-Stokes codes can be used for these flows to study trends and viscous effects as well as determine flow patterns; but further research and development is needed before they can be used as production tools for nozzle performance calculations. The GIM formulation, numerical scheme, and computer code are described. The actual SSME nozzle computation showing grid points, flow contours, and flow parameter plots is discussed. The computer system and run times/costs are detailed.

Chan, J. S.↗
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 NTRS1985

Absorbing boundary conditions for exterior problems

Elliptic and hyperbolic problems in unbounded regions are considered. These problems, when one wants to solve them numerically, have the difficulty of prescribing boundary conditions at infinity. Computationally, one needs a finite region in which to solve these problems. The corresponding conditions at infinity imposed on the finite distance boundaries should dictate the boundary condition at infinity and be accurate with respect to the interior numerical scheme. Such boundary conditions are commonly referred to as absorbing boundary conditions. A treatment is given of these boundary conditions for wave-like equations.

Hariharan, S. I.↗
NASA NTRS1985

A procedure for the calculation of supersonic flows with strong viscous-inviscid interaction

The present paper is concerned with the calculation of strong viscous-inviscid interactions in two-dimensional laminar supersonic flows with and without separation. The equations solved are the so-called parabolized Navier-Stokes equations. The streamwise pressure gradient term is written as a combination of a forward and a backward difference to provide a path for upstream propagation of information. Global iteration is utilized to repeatedly update the pressure field from an initial guess until convergence is achieved. The numerical scheme employed is a new alternating direction explicit (ADE) procedure which is used as an alternative to the more difficult to program multigrid strategy to accelerate convergence. Results are presented for flows past two flat plate related bodies.

Barnett, M.↗
NASA NTRS1985

Contributions to the finite element solution of the fan noise radiation problem

The radiation of fan generated noise to the far field from a nacelle of realistic geometry is investigated using the finite element method. Several innovations have been introduced to minimize the computational requirements and create a highly efficient numerical scheme. The innovations include: (1) formulation of the problem in terms of velocity potential and density in such a way that no inlet mean flow velocity derivatives are required in the field equations, (2) the use of 'wave envelope' elements in an outer region permitting a grid much coarser than would be used for conventional finite elements, (3) the use of a mesh which deforms with an increase of forward flight speed so that mesh lines are always lines of constant phase and rays for a point source, permitting the use of wave envelope elements and simple boundary conditions for any case of forward velocity, (4) an efficient scheme for introducing the noise source via modal amplitude coefficients, and (5) the use of a frontal solution technique which for physically realistic problems drastically reduces the active storage requirements. The finite element scheme is outlined, as are the specific details of the innovations. Results are given for cases where comparable experimental data are available.

Eversman, W.↗
NASA NTRS1985

A three-dimensional boundary-layer analysis including heat-transfer and blade-rotation effects

An analysis has been developed for calculating three-dimensional viscous flows over arbitrary bodies through the solution of complete three-dimensional boundary-layer equations. The governing equations include the effect of blade rotation and are written in nonorthogonal surface-fitted coordinates to allow maximum flexibility. A coupled finite-difference numerical scheme is used for solving the governing equations. The accuracy of this scheme has been assessed through comparisons with known analytical/experimental results for several test cases. The present analysis has then been applied to compute viscous flows over a turbine endwall and airfoil suction surface, and fairly good correlation has been obtained with available experimental data for surface heat transfer and viscous streamline skewing.

Vatsa, V. N.↗