SEARCH · Engineering Papers
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.
The use of the finite element and finite difference methods often leads to the problem of solving large, sparse, positive definite systems of linear equations. MACSYMA plays a major role in the generation of formulas representing the time required for execution of the dissection algorithms. The use of MACSYMA in the generation of those formulas is described.
The free surface shapes of circular liquid jets impinging normal to sharp-edged disks in zero gravity are determined. Zero gravity drop tower experiments yielded three distinct flow patterns that were classified in terms of the relative effects of surface tension and inertial forces. An order of magnitude analysis was conducted that indicated regions where viscous forces were not significant in the computation of free surface shapes. The free surface analysis was simplified by transforming the governing potential flow equations and boundary conditions into the inverse plane, where the stream function and velocity potential became the coordinates. The resulting nonlinear equations were solved by standard finite difference methods, and comparisons were made with the experimental data for the inertia dominated regime.
An implicit finite-difference method has been developed for computing the flow in the near field of a fuel injector as part of a broader study of the effects of fuel injector geometry on fuel-air mixing and combustion. Detailed numerical results have been obtained for cases of laminar and turbulent flow without base injection, corresponding to the supersonic base flow problem. These numerical results indicated that the method is stable and convergent, and that significant savings in computer time can be achieved, compared with explicit methods.
Hypersonic viscous flow over spherically blunted cones of large half angle is computed at small angles of attack in the plane of symmetry of the flow field. Time-dependent viscous shock-layer equations in body-oriented coordinates are used to describe the flow field. The shock wave is treated as a discontinuity, across which the Rankine-Hugoniot relations are used to compute the flow conditions behind the shock. A time-marching second-order finite-difference method is used to solve the equations for a perfect gas. The local CFL (Courant-Friedrich-Lewy) time increment is used to advance the solution in time at each grid point. A fourth-order damping is used to damp the oscillations in the flow quantities. The numerical results of the present analysis for quantities such as shock standoff distance, surface-pressure distribution, and heating rates compare well with the existing theoretical and experimental results.
A numerical study to determine the temperature distribution in the guide vanes of a radial inflow turbine is presented. A computer program has been developed to calculate the temperature distribution when the vanes are cooled internally using a combination of impingement and film cooling techniques. The study is based on the use of the finite difference method in a two dimensional heat conduction problem. The results are then compared to determine the best cooling configuration for a certain coolant to primary mass flow ratio.
Turbulent separating flows over axisymmetric afterbody-boattail configurations with solid sting plume simulators are computed with a time-dependent finite-difference method to solve the compressible Navier-Stokes equations. The Reynolds stress terms are replaced with a two-layer eddy viscosity model including a relaxation formula to model the nonequilibrium effects of the separated flow. The mesh is alined with the boattail body through an analytic transformation which accommodates a wide variety of boattail geometries. Numerical results for a series of boattail geometries over a wide range of Reynolds number (140,000 to 140 million) are presented and, when possible, compared with experimental data or independent numerical results.
A numerical study to determine the temperature distribution in the guide vane blades of a radial inflow turbine is presented. A computer program was developed which permits the temperature distribution to be calculated when the blade is cooled internally using a combination of impingement and film cooling techniques. The study is based on the use of the finite difference method in a two dimensional heat conduction problem. The results are then compared to determine the best cooling configuration for a certain coolant to primary mass flow ratio.
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.
A computer program is described which provides solutions of two dimensional equations appropriate to laminar and turbulent boundary layers for boundary conditions with an external flow which fluctuates in magnitude. The program is based on the numerical solution of the governing boundary layer equations by an efficient two point finite difference method. An eddy viscosity formulation was used to model the Reynolds shear stress term. The main features of the method are briefly described and instructions for the computer program with a listing are provided. Sample calculations to demonstrate its usage and capabilities for laminar and turbulent unsteady boundary layers with an external flow which fluctuated in magnitude are presented.
The behavior of unsteady aerodynamic loadings on airfoils oscillating in transonic flow has been investigated numerically with particular attention given to supercritical airfoil sections. A previously developed finite difference method, which is based on the full potential equation and which uses a quasi-conservative scheme for proper capture of a shock wave motion, was employed for the present study. The unsteady aerodynamic pressure and load distributions on several different airfoil sections are presented with particular emphasis on the effects of free-stream Mach number, reduced frequency, and mean angle of attack. These parameters are demonstrated to have a significant effect on the behavior of the unsteady aerodynamic loadings. Comparisons of the present calculations with the exact inviscid solution and with the experimental results are also presented.
For the problem of predicting one-dimensional heat transfer between conducting and radiating mediums by an implicit finite difference method, four different formulations were used to approximate the surface radiation boundary condition while retaining an implicit formulation for the interior temperature nodes. These formulations are an explicit boundary condition, a linearized boundary condition, an iterative boundary condition, and a semi-iterative boundary method. The results of these methods in predicting surface temperature on the space shuttle orbiter thermal protection system model under a variety of heating rates were compared. The iterative technique caused the surface temperature to be bounded at each step. While the linearized and explicit methods were generally more efficient, the iterative and semi-iterative techniques provided a realistic surface temperature response without requiring step size control techniques.
The flow field produced by the intersection of two plane solid surfaces in a supersonic stream is a complex interference flow. These flows can be fully compressive, fully expansive, or of mixed compression-expansion nature. This paper presents a comparison of the experimentally obtained flow-field structure in an axial corner with that predicted numerically by using a shock-capturing finite-difference method. The effect of sweep and surface deflection are evaluated, and the general influence of each is presented for the three classes of corner flow. The results show that the numerical method is a valuable aid in understanding the flow structure for simple configurations. In addition, confidence in the numerical method is gained for use in solving more general three-dimensional configurations where the flow is nonconical and several wave interaction may be presented.
The equations governing the laminar and turbulent flows of reacting gas mixtures in chemical equilibrium over axially symmetric blunt bodies at small angles of attack are developed and presented in the unsteady conservative form. Solutions are obtained in the planes of symmetry of the flow field for the typical conditions encountered by a probe entering the Jovian atmosphere at small angles of attack. The eddy-viscosity is approximated by a two-layer model. The shock wave is treated as a discontinuity across which the shock relations are used to compute the flow conditions behind the shock. A time-asymptotic finite-difference method is used to solve the equations. The zero angle of attack results for a hyperboloid are compared with the existing results but neither experimental nor computational results are available for comparison with the present results at angle of attack.
The use of grooved evaporator surfaces in heat pipes has increased in popularity in the past few years primarily due to the reproducibility achievable with grooved walls and the relatively low costs of the threading or extrusion processes involved in their production. The present study combines both analyses and experiments on square groove geometries, with special emphasis on overcoming the limitations of earlier analyses with finite-difference methods and groove-fillet hydrodynamic simplifications. The groove fillet, which has in previous analyses been assumed constant in radius of curvature, is permitted to change in thickness and curvature consistent with hydrodynamics and heat loss from the groove. A model is developed for accurate determination of the effect of constriction resistance on groove performance. The grooved-surface tests to be conducted are briefly described which will provide data under closely controlled operation to allow comparison and verification of the analyses.
An implicit finite-difference method has been developed to compute two-dimensional, turbulent, blunt body flows with an impinging shock wave. The full time-averaged Navier-Stokes equations are solved with algebraic eddy viscosity and turbulent Prandtl number models employed for shear stress and heat flux. The irregular-shaped bow shock is treated as a discontinuity across which the Rankine-Hugoniot equations are applied. A Type III turbulent shock interference flow field has been computed and the numerical results compare favorably with existing experimental data. In addition, comparisons are made between the present implicit code and a previous explicit code.
A new trigonometric approach to the finite difference calculus was applied to the problem of beam buckling as represented by virtual work and equilibrium equations. The trigonometric functions were varied by adjusting a wavelength parameter in the approximating Fourier series. Values of the critical force obtained from the modified approach for beams with a variety of boundary conditions were compared to results using the conventional finite difference method. The trigonometric approach produced significantly more accurate approximations for the critical force than the conventional approach for a relatively wide range in values of the wavelength parameter; and the optimizing value of the wavelength parameter corresponded to the half-wavelength of the buckled mode shape. It was found from a modal analysis that the most accurate solutions are obtained when the approximating function closely represents the actual displacement function and matches the actual boundary conditions.
This code solves the two-dimensional, transonic, small-disturbance equations for flow past lifting airfoils in both free air and various wind-tunnel environments by using a variant of the finite-difference method. A description of the theoretical and numerical basis of the code is provided, together with complete operating instructions and sample cases for the general user. In addition, a programmer's manual is also presented to assist the user interested in modifying the code. Included in the programmer's manual are a dictionary of subroutine variables in common and a detailed description of each subroutine.
Numerical solutions are presented for the flow over a spherically blunted cone and hyperboloid with massive surface blowing. Time-dependent viscous shock-layer equations are used to describe the flow field. The boundary conditions on the body surface include a prescribed blowing-rate distribution. The governing equations are solved by a time-asymptotic finite-difference method. Results presented here are only for a perfect gas-type flow at zero angle of attack. Both laminar and turbulent flow solutions are obtained. It is found that the effect of the surface blowing on the laminar flow field is to smooth out the curvature discontinuity at the sphere-cone juncture point, which results in a positive pressure gradient over the body. The shock slope increases on the downstream portion of the body as the surface blowing rate is increased. The turbulent flow with surface blowing is found to redevelop a boundary-layer-like region near the surface. The effects of this boundary-layer region on the flow field and heating rates are discussed.