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.

SourceResearch connectionFrom yearThrough year

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

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1974

Lateral conduction effects on heat-transfer data obtained with the phase-change paint technique

A computerized tool, CAPE, (Conduction Analysis Program using Eigenvalues) has been developed to account for lateral heat conduction in wind tunnel models in the data reduction of the phase-change paint technique. The tool also accounts for the effects of finite thickness (thin wings) and surface curvature. A special reduction procedure using just one time of melt is also possible on leading edges. A novel iterative numerical scheme was used, with discretized spatial coordinates but analytic integration in time, to solve the inverse conduction problem involved in the data reduction. A yes-no chart is provided which tells the test engineer when various corrections are large enough so that CAPE should be used. The accuracy of the phase-change paint technique in the presence of finite thickness and lateral conduction is also investigated.

Maise, G.↗

NASA NTRS1974

Effects of mean winds and horizontal temperature gradients on solar and lunar semidiurnal tides in the atmosphere

A description is given of specifications of mean seasonal zonal winds and temperatures. Equations for tidal perturbations of the basic states are discussed along with a numerical scheme which is, essentially, the scheme described by Lindzen and Kuo (1969). Results for the solar semidiurnal tide are considered along with an investigation of the lunar semidiurnal tide.

Lindzen, R. S.↗

NASA NTRS1974

Three dimensional flows around airfoils with shocks

The present work describes a mathematical model and numerical scheme for the computer-aided calculation of two- and three-dimensional transonic flow over an isolated yawed wing with oblique shock waves and a trailing vortex sheet. The flow is modeled by the potential equation for irrotational flow which is hyperbolic at supersonic points and elliptic at subsonic points. A coordinate-invariant difference scheme is used in which retarded difference formulas are constructed to conform with the local flow direction. The resulting 'rotated' difference scheme allows complete flexibility in the choice of a coordinate system. Shock waves are located automatically in the form of compression bands spread over a few mesh widths. The scheme has proven to be stable and convergent throughout the transonic range. Calculations have been performed for Mach number up to 1.2 and yaw angles up to 60 deg, the likely operating range of a yawed-wing transport designed to fly at supersonic speeds. Calculations becomes less accurate towards upper end of range, because the difference scheme is first-order accurate in the supersonic range.

Jameson, A.↗

NASA NTRS1974

Method of Characteristics Calculations and Computer Code for Materials with Arbitrary Equations of State and Using Orthogonal Polynomial Least Square Surface Fits

A numerical scheme using the method of characteristics to calculate the flow properties and pressures behind decaying shock waves for materials under hypervelocity impact is developed. Time-consuming double interpolation subroutines are replaced by a technique based on orthogonal polynomial least square surface fits. Typical calculated results are given and compared with the double interpolation results. The complete computer program is included.

Chang, T. S.↗

NASA NTRS1975

Numerical solutions of supersonic and hypersonic laminar flows over a two-dimensional compression corner

An efficient time-splitting, second-order accurate, numerical scheme is used to solve the complete Navier-Stokes equations for supersonic and hypersonic laminar flow over a two-dimensional compression corner. A fine, exponentially stretched mesh spacing is used in the region near the wall for resolving the viscous layer. Good agreement is obtained between the present computed results and experimental measurement for a Mach number of 14.1, a Reynolds number of 104,000, and wedge angles of 15, 18, and 24 deg. The details of the pressure variation across the boundary layer are given, and a correlation between the leading edge shock and the peaks in surface pressure and heat transfer is observed.

Hung, C. M.↗

NASA NTRS1975

Energy transfer by radiation in non-grey atomic gases in isothermal and non-isothermal slabs

A multiband model for the absorption coefficient of atomic hydrogen-helium plasmas is constructed which includes continuum and line contributions. Emission from 28 stronger lines of 106 that have been screened is considered, of which 21 are from hydrogen and 7 belong to helium, with reabsorption due to line-line, line-continuum overlap accurately accounted for. The model is utilized in the computation of intensities and fluxes from shock-heated slabs of 85% H2-15% He mixtures for slab thicknesses from 1 to 30 cm, temperature from 10,000 to 20,000 K, and for different densities. In conjunction with the multiband model, simple numerical schemes have been devised which provide a quick and comprehensive way of computing radiative energy transfer in nonisothermal and nongrey gases.

Poon, P. T. Y.↗

NASA NTRS1976

Influence of surface roughness and waviness on film thickness and pressure distribution in elastohydrodynamic contacts

The Christensen theory of a stochastic model for hydrodynamic lubrication of rough surfaces was extended to elastohydrodynamic lubrication between two rollers. Solutions for the reduced pressure at the entrance as a function of the ratio of the average nominal film thickness to the rms surface roughness, were obtained numerically. Results were obtained for purely transverse as well as purely longitudinal surface roughness for cases with or without slip. The reduced pressure was shown to decrease slightly by considering longitudinal surface roughness. The same approach was used to study the effect of surface roughness on lubrication between rigid rollers and lubrication of an infinitely wide slider bearing. Using the flow balance concept, the perturbed Reynolds equation, was derived and solved for the perturbed pressure distribution. In addition, Cheng's numerical scheme was modified to incorporate a single two-dimensional elastic asperity on the stationary surface. The perturbed pressures obtained by these three different models were compared.

Chow, L. S. H.↗

NASA NTRS1976

Numerical simulation of supersonic and hypersonic turbulent compression corner flows using relaxation models

Relaxation turbulence eddy viscosity models are incorporated to solve the complete Navier-Stokes equations for supersonic and hypersonic flows over a two-dimensional compression corner. The system of equations is solved by a time-split, second-order accurate numerical scheme. Details of relaxation process are studied, and several relaxation models are tested and compared. Good improvement in the prediction of upstream pressure propagation is obtained for a Mach number of 2.96, and Reynolds number of 10,000,000, with wedge angle of 25 deg. However, the application of the relaxation models to hypersonic flow at Mach number 8.66 and Reynolds number of 22,000,000, with highly cooled wall, shows unfavorable effects on heat transfer and skin friction in the reattachment and recompression regions.

Hung, C. M.↗

NASA NTRS1976

A two-dimensional model for the transport of pollutants in an urban basin

The distribution of an inert chemical species, carbon monoxide (CO), emitted mainly from automobile exhaust is investigated. The transport model, in analogy with the shallow water theory in fluid dynamics, considers variation of all physical quantities in the horizontal direction below the temperature inversion layer. Pollutants are found to be carried primarily by the wind; turbulent diffusion in a normal day plays only a minor role. The concentration of CO predicted by the present model for the entire Los Angeles basin is compared with observed data at nine stations. Accuracy of four different numerical schemes, the effect of turbulent diffusivity, and the source strengths are examined.

Liu, C. Y.↗

NASA NTRS1977

Improved methods for large-eddy simulations of turbulence

By using Fourier transforms for evaluating spatial derivatives, the accuracy of the large-eddy simulation of homogeneous isotropic turbulence is improved. Especially, the treatment of certain terms that arise in filtering the equations is considerably improved in both speed and accuracy. Use of vorticity as the principal variable is shown to be a viable and potentially useful alternative to the primitive variables. A method of deriving conservation properties of numerical schemes is given which is much simpler than previous methods and is widely applicable. The methods are applied to the computation of homogeneous isotropic turbulence, and it is found that the subgrid scale model is improved by using finite differences in place of 'exact' derivatives.

Mansour, N. N.↗

NASA NTRS1977

Numerical solution of supersonic laminar flow over a three-dimensional compression corner

A newly developed, rapid numerical scheme is extended to three dimensions to solve the complete Navier-Stokes equations for a supersonic, laminar flow over a compression corner with sidewall effects. The program is coded so that it can solve for a general curved ramp surface geometry such as found in inlets and fuselage-wing-flap junctions. A test case of Mach 3.0 flow is calculated. In regions where three-dimensional effects are small, good agreement is obtained between the present calculation and previous two-dimensional solutions. In other regions, the results show complex three-dimensional flow-field interactions including shock-shock and shock/boundary-layer interactions

Hung, C. M.↗

NASA NTRS1977

The separated turbulent boundary layer over a wavy wall

A study and application of the fourth order spline collocation procedure, numerical solution of boundary layer like differential equations, is presented. A simple inversion algorithm for the simultaneous solution of the resulting difference equations is given. Particular attention is focused on the boundary condition representation for the spline second derivative approximations. Solutions using the spline procedure, as well as the three point finite difference method, are presented for several model problems in order to assess and improve the spline numerical scheme. Application of the resulting algorithm to the incompressible laminar self similar boundary layer equations is presented.

Polak, A.↗

NASA NTRS1978

Numerical solution of a three-dimensional shock wave and turbulent boundary-layer interaction

A rapid numerical scheme is used to solve the complete mass-averaged Navier-Stokes equations for supersonic turbulent flow over a three-dimensional compression corner. A simple eddy viscosity model is developed, and the interaction of a swept shock wave and a three-dimensional turbulent boundary layer is studied. Good agreement is obtained between the present results and experimental measurements for the case of a wedge with an angle of 6 deg on a flat-plate sidewall. For the case of a 12-deg wedge angle, the computed results do not show the existence of a peak pressure found experimentally. However, the range of interaction, the plateau pressure, and the peak heat transfer are closely predicted for all cases. The high heat transfer near the axial corner is due to the thinning of the boundary layer and inflow of fresh high-momentum fluid. The heat transfer is relieved through pressure reduction and boundary-layer thickening.

Hung, C. M.↗

NASA NTRS1977

A two-dimensional model for the transport of pollutants in an urban basin

The distribution of carbon monoxide emitted mainly from an automobile exhaust is investigated. Carbon monoxide is assumed to be chemically inert. The transport model, in analogy with the shallow-water theory in fluid dynamics, considers variation of all physical quantities in the horizontal direction below the temperature inversion layer. Pollutants are found to be carried primarily by the wind; turbulent diffusion in a normal day plays only a minor role. The concentration of CO predicted by the present model for the entire Los Angeles basin is compared with observed data at nine stations. Accuracy of four different numerical schemes, the effect of turbulent diffusivity, and the source strengths are examined.

Liu, C. Y.↗

NASA NTRS1979

A calculation procedure for viscous flow in turbomachines, volume 1

A method for analyzing the nonadiabatic viscous flow through turbomachine rotors is presented. The field analysis is based upon the numerical integration of the full incompressible stream function vorticity form of the Navier-Stokes equations, together with the energy equation, over the rotor blade-to-blade stream channels. The numerical code used to solve the governing equations employs a nonorthogonal boundary fitted coordinate system that suits the most complicated blade geometries. A numerical scheme is used to carry out the necessary integration of the elliptic governing equations. The flow characteristics within the rotor of a radial inflow turbine are investigated over a wide range of operating conditions. The calculated results are compared to existing experimental data. The flow in a radial compressor is analyzed in order to study the behavior of viscous flow in diffusing cascades. The results are compared qualitatively to known experimental trends. The solution obtained provides insight into the flow phenomena in this type of turbomachine. It is concluded that the method of analysis is quite general and gives a good representation of the actual flow behavior within turbomachine passages.

Khalil, I.↗

NASA NTRS1978

A stochastic analysis of terrain evaluation variables for path selection

A stochastic analysis was performed on the variables associated with the characteristics of the terrain encountered by a roving system with an autonomous navigation system. A laser rangefinder is employed to detect terrain features at ranges up to 75 m. Analytic expressions and a numerical scheme were developed to calculate the variance of data on these four variables: (1) body clearance, (2) in-path slope, (3) tilt slope, and (4) wheel deviation. The variance is due to noise in the range data. It was found that the standard deviation of these terrain variables is large enough to warrant the use of a safety margin to aid the roving vehicle in avoiding high risk areas.

Donohue, J. G.↗

NASA NTRS1979

The application of finite element techniques to acoustic transmission in lined ducts with flow

The finite element method (FEM) is used to analyze the propagation of sound in two-dimensional nonuniform ducts carrying a compressible subsonic mean flow. Galerkin and residual least squares (RLS) methods with natural and forced boundary conditions are considered. The accuracy of FEM results for the eigenvalue and transmission problems is assessed by comparison with alternative numerical schemes for nonuniform ducts. The results presented and those from associated investigations indicate that modal coupling is a significant feature of the acoustic field, especially at high Mach numbers. A multimodal model therefore appears to be essential if any reliable conclusions are to be drawn in the context of turbofan inlet regions. Improvements to the eigenvalue scheme following the implementation of higher-order Hermitian elements indicate a similar modification for the transmission problem.

Astley, R. J.↗

NASA NTRS1979

Numerical computation of aeroelastically corrected transonic loads

A numerical scheme is presented for the computation of transonic aerodynamic loads on flexible wings. The method consists of iteratively applying the loads computed by a 3D transonic aerodynamics code to a structural model to obtain elastic twist, and then recomputing the loads. Because this iteration is performed concurrently with the iterations performed in computing the aerodynamics, flexible loads are obtained in roughly the same amount of computing time as required to obtain rigid loads. Applications of this method to a flexible supercritical transonic transport wing are presented and compared with model test data.

Chipman, R.↗