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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
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.↗
NASA NTRS1979

Numerical study of an explosion in a non-homogeneous medium with and without magnetic fields

A version of the two step Lax-Wendroff difference method with second order accuracy is used to seek solutions of the unsteady magnetohydrodynamic (MHD) equations for the study of an explosion in a non-homogeneous medium with and without magnetic fields. The explosion is initiated by introducing a finite amount of energy in a small volume of gas which leads to an instantaneous increase of gas temperature. Two types of magnetic field configurations (i.e. open and closed) are considered to illustrate the dependence and differences of magnetic and a non-magnetic flow motion. Numerical results show the strong dependence of induced MHD flow field on the magnetic field configuration and strength. The development and decay of the fast and slow MHD shock waves, as well as the ordinary gasdynamic shock waves, are also represented well by the present computer simulations. In conclusion, we have demonstrated that the numerical scheme we have used is a reliable one for studying 2-dimensional MHD problems in non-homogeneous medium.

Han, S. M.↗
NASA NTRS1980

Iterative methods based upon residual averaging

Iterative methods for solving boundary value problems for systems of nonlinear partial differential equations are discussed. The methods involve subtracting an average of residuals from one approximation in order to arrive at a subsequent approximation. Two abstract methods in Hilbert space are given and application of these methods to quasilinear systems to give numerical schemes for such problems is demonstrated. Potential theoretic matters related to the iteration schemes are discussed.

Neuberger, J. W.↗
NASA NTRS1980

Numerical techniques for solving nonlinear instability problems in smokeless tactical solid rocket motors

The selection of a satisfactory numerical method for calculating the propagation of steep fronted shock life waveforms in a solid rocket motor combustion chamber is discussed. A number of different numerical schemes were evaluated by comparing the results obtained for three problems: the shock tube problems; the linear wave equation, and nonlinear wave propagation in a closed tube. The most promising method--a combination of the Lax-Wendroff, Hybrid and Artificial Compression techniques, was incorporated into an existing nonlinear instability program. The capability of the modified program to treat steep fronted wave instabilities in low smoke tactical motors was verified by solving a number of motor test cases with disturbance amplitudes as high as 80% of the mean pressure.

Baum, J. D.↗
NASA NTRS1981

A new payload integration method

This paper introduces a new payload integration method. The system equations of motion are derived in terms of interface-restrained booster and payload modes. These system equations are directly integrated, avoiding the solution of an expensive eigenvalue problem. A modified Newmark-Chan-Beta numerical scheme is used to perform this integration. Due to the unique form of the system equations it is possible to achieve significant savings. A similar savings is feasible for the computation of the internal loads. This new approach does not involve approximations. Although completely general, the presented technique is most effective when the number of interface degrees of freedom is relatively small. Both determinate and indeterminate interfaces are allowed.

Engels, R. C.↗
NASA NTRS1981

Analysis of cracks emanating from a circular hole in unidirectional fiber reinforced composites, part 2

An analytical method is developed for cracks emanating from a circular hole in an off-axis unidirectional fiber-reinforced composite. The method which is formulated by using conservation laws of elasticity and fundamental relationships in anisotropic fracture mechanics, provides a convenient and accurate means to examine the complicated crack behavior, when used in conjunction with a suitable numerical scheme such as the finite element method. The formulation is eventually reduced to a system of linear algebraic equations of mixed-mode stress intensity factors. Fracture parameters, describing crack-tip deformation and fracture in the composite, are obtained explicitly. Effects of material anisotropy and crack/hole geometry are examined also. Of particular interest are the energy release rates associated with crack extension; their values are evaluated for various cases. Results show that mixed-mode stress intensity factors and energy release rates associated with the cracks emanating from a hole change very appreciably with fiber orientation in the composite. K sub 1 and G increase monotonically with increasing theta; but K sub 2 reaches its maximum at theta = 45 deg, and then decreases gradually as theta increases further.

Wang, S. S.↗
NASA NTRS1982

Conjugated heat transfer in a laminar boundary layer with heat source at the wall

This work presents the solution of the temperature field in a two-dimensional laminar incompressible flow over a conducting solid plate with a line heat source located at the fluid-solid interface perpendicular to the flow direction. A numerical scheme was used to obtain the temperature profiles as a function of the source strength, and of the properties of the fluid and the solid. The heat conduction and forced convection in the fluid and the heat conduction in the solid were solved for the case of moderate temperature rise, where the assumption of constant properties applies. The model enables the improvement of an instrument for the detection of boundary layer separation. It was found that for the actual parameters of the separation detector, a distance of 4 to 24 mm between the sensors gives an indication of 70 percent of the maximum temperature difference.

Brosh, A.↗