Keith, Theo G., Jr. — Page 5 | Engineering Papers

NASA NTRS1990

Analysis and design of optimized truncated scarfed nozzles subject to external flow effects

Rao's method for computing optimum thrust nozzles is modified to study the effects of external flow on the performance of a class of exhaust nozzles. Members of this class are termed scarfed nozzles. These are two-dimensional, nonsymmetric nozzles with a flat lower wall. The lower wall (the cowl) is truncated in order to save weight. Results from a parametric investigation are presented to show the effects of the external flowfield on performance.

Shyne, Rickey J.↗

NASA NTRS1990

Numerical simulation of rarefied gas flow through a slit

Two different approaches, the finite-difference method coupled with the discrete-ordinate method (FDDO), and the direct-simulation Monte Carlo (DSMC) method, are used in the analysis of the flow of a rarefied gas from one reservoir to another through a two-dimensional slit. The cases considered are for hard vacuum downstream pressure, finite pressure ratios, and isobaric pressure with thermal diffusion, which are not well established in spite of the simplicity of the flow field. In the FDDO analysis, by employing the discrete-ordinate method, the Boltzmann equation simplified by a model collision integral is transformed to a set of partial differential equations which are continuous in physical space but are point functions in molecular velocity space. The set of partial differential equations are solved by means of a finite-difference approximation. In the DSMC analysis, three kinds of collision sampling techniques, the time counter (TC) method, the null collision (NC) method, and the no time counter (NTC) method, are used.

Keith, Theo G., Jr.↗

NASA NTRS1990

A numerical simulation of the flow in the diffuser of the NASA Lewis icing research tunnel

The flow in the diffuser section of the Icing Research Tunnel at the NASA Lewis Research Center is numerically investigated. To accomplish this, an existing computer code is utilized. The code, known as PARC3D, is based on the Beam-Warming algorithm applied to the strong conservation law form of the complete Navier-Stokes equations. The first portion of the paper consists of a brief description of the diffuser and its current flow characteristics. A brief discussion of the code work follows. Predicted velocity patterns are then compared with the measured values.

Addy, Harold E., Jr.↗

NASA NTRS1990

Analysis and design of optimized truncated scarfed nozzles subject to external flow effects

Rao's method for computing optimum thrust nozzles is modified to study the effects of external flow on the performance of a class of exhaust nozzles. Members of this class are termed scarfed nozzles. These are two-dimensional, nonsymmetric nozzles with a flat lower wall. The lower wall (the cowl) is truncated in order to save weight. Results from a parametric investigation are presented to show the effects of the external flowfield on performance.

Shyne, Rickey J.↗

NASA NTRS1990

A numerical simulation of the flow in the diffuser of the NASA Lewis Icing Research Tunnel

The flow in the diffuser section of the Icing Research Tunnel at the NASA Lewis Research Center is numerically investigated. To accomplish this, an existing computer code is utilized. The code, known as PARC3D, is based on the Beam-Warming algorithm applied to the strong conservation law form of the complete Navier-Stokes equations. The first portion of the paper consists of a brief description of the diffuser and its current flow characteristics. A brief discussion of the code work follows. Predicted velocity patterns are then compared with the measured values.

Addy, Harold E., Jr.↗

NASA NTRS1990

Time domain flutter analysis of cascades using a full-potential solver

A time domain approach is used to determine the dynamic aeroelastic stability of a cascade of blades. The structural model for each blade is a typical section with two degrees of freedom. The aerodynamic model is the unsteady, two-dimensional, full-potential flow through the cascade of airfoils. The unsteady equations of motion for the structure and the fluid are integrated simultaneously in time starting with the steady flowfield and a small initial disturbance applied to the airfoils. The motion of each blade is analyzed to determine the aeroelastic stability of the cascade. The effect of interblade phase angle is included in the analysis by allowing each blade to have an independent motion and considering a number of blade passages. Calculations are made using an airfoil section and structural parameters that are representative of a propfan. The results are compared with those from a separate frequency domain analysis. Good agreement between the results is observed. With the time domain approach, it is possible to consider nonlinear structural models and nonlinear force-displacement relations. The method allows a realistic simulation of the motion of the fluid and the cascade blades for a better physical understanding and it also has the potential for saving computational time when compared to the frequency domain approach for the flutter analysis of cascades.

Bakhle, Milind A.↗

NASA NTRS1989

Numerical modeling of a vortex stabilized arcjet thruster

A numerical method to solve the equations governing a vortex stabilized arcjet thruster is being developed. The model will allow the effects of swirling flow and geometry on arcjet thruster performance to be determined. The propellant flow equations, i.e., the axisymmetric, thin layer, Navier-Stokes equations, are solved using a Gauss-Siedel line-relaxation procedure. An implicit FTCS method is used to solve the electromagnetic field equations. A grid generation scheme was developed for an arbitrary arcjet geometry. The model will allow the radial and axial components of velocity and current distributions to be determined from a region upstream of the cathode, through the constrictor, to the exit plane of the nozzle.

Keith, Theo G., Jr.↗

NASA NTRS1989

The effect of test-cell pressure on resistojet nozzle flow

Previous experimental work has shown that measured resistojet thrust decreases from that obtained at hard vacuum conditions as the test-cell pressure rises above 0.001 torr. Thrust losses have been observed for both cold and heated flow conditions, and the most significant losses have been experienced using thrusters with low Reynolds number flow and high area ratio nozzles. In order to further investigate nozzle flow characteristics, a pressure probe having four degrees of freedom has been used to obtain stagnation pressure surveys across the nozzle exit planes of four resistojets. The surveys show a change in the ratio of the supersonic core to the viscous boundary layer flow areas as the test-cell pressure increases. The surveys are also used for detecting whether an oblique shock is present in an overexpanded nozzle flow. Thruster temperature measurements and nozzle exit plane pressure surveys indicate that thrust losses are the combined result of convection heat losses and nozzle flow momentum effects.

Penko, Paul F.↗

NASA NTRS1989

Rarefied gas flow through two-dimensional nozzles

A kinetic theory analysis is made of the flow of a rarefied gas from one reservoir to another through two-dimensional nozzles with arbitrary curvature. The Boltzmann equation simplified by a model collision integral is solved by means of finite-difference approximations with the discrete ordinate method. The physical space is transformed by a general grid generation technique and the velocity space is transformed to a polar coordinate system. A numerical code is developed which can be applied to any two-dimensional passage of complicated geometry for the flow regimes from free-molecular to slip. Numerical values of flow quantities can be calculated for the entire physical space including both inside the nozzle and in the outside plume. Predictions are made for the case of parallel slots and compared with existing literature data. Also, results for the cases of convergent or divergent slots and two-dimensional nozzles with arbitrary curvature at arbitrary knudsen number are presented.

De Witt, Kenneth J.↗

NASA NTRS1989

Numerical simulation of axisymmetric turbulent flow in combustors and diffusers

Numerical studies of turbulent flow in an axisymmetric 45-deg-expansion combustor and bifurcated diffuser are presented. The Navier-Stokes equations incorporating a k-epsilon model were solved in a nonorthogonal curvilinear coordinate system. A zonal-grid method, where the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. Boundary overlap and interpolating techniques were used, and an adjustment of the flow variables was required to assure conservation of mass flux. Three finite-differencing methods (hybrid, quadratic upwind, and skew upwind) were used to represent the convection terms. Results were compared with existing experimental data. In general, good agreement between predicted and measured values was obtained.

Yung, Chain-Nan↗

NASA NTRS1989

Investigation of the flow in the diffuser section of the NASA Lewis icing research tunnel

The flow in the diffuser section of the Icing Research Wind Tunnel at NASA Lewis Research Center is investigated using both tunnel calibration measurements and numerical simulation techniques. Local pressure and temperature measurements are made to establish velocity and temperature profiles in the diffuser of the tunnel. These profiles are compared with similar measurements made prior to renovating the equipment which generates the tunnel's icing cloud. This comparison indicates the manner in which this change affected the flow. The measured data were also compared with a numerical simulation of the flow to help understand how such changes may favorably alter the tunnel flow.

Addy, Harold E., Jr.↗

NASA NTRS1989

Thermal analysis of engine inlet anti-icing systems

A hot air anti-icing system of a gas turbine engine inlet is analyzed numerically. A three-dimensional potential flow code, which accounts for compressibility effects, is used to determine the flowfield in and around the inlet. A particle trajectory code is developed using a local linearization technique. The trajectory code is used to calculate local water impingement rates. Energy balances are performed on both the surface runback water and the metallic skin to determine their temperature distributions. A variety of test cases are considered in order to validate the various numerical components of the process as well as to demonstrate the procedure.

Keith, Theo G., Jr.↗

NASA NTRS1989

Investigation of the flow in the diffuser section of the NASA Lewis Icing Research Tunnel

The flow in the diffuser section of the Icing Research Wind Tunnel at NASA Lewis Research Center is investigated using both tunnel calibration measurements and numerical simulation techniques. Local pressure and temperature measurements are made to establish velocity and temperature profiles in the diffuser of the tunnel. These profiles are compared with similar measurements made prior to renovating the equipment which generates the tunnel's icing cloud. This comparison indicates the manner in which this change affected the flow. The measured data were also compared with a numerical simulation of the flow to help understand how such changes may favorably alter the tunnel flow.

Addy, Harold E., Jr.↗

NASA NTRS1989

Application of a full-potential solver to bending-torsion flutter in cascades

A two-dimensional, unsteady, full-potential cascade aerodynamic code has been used for flutter calculations. The code, which uses a time-marching implicit scheme with internal Newton-iterations at each time step, allows oscillatory blade motions with phase lag between adjacent blades to be modeled using a single inter-blade passage. This aerodynamic code has been used with a two-degrees-of-freedom typical section structural model to determine the aeroelastic stability of a bladed disk. The unsteady pressure results calculated from the aerodynamic code have been verified by comparison with results from the linear and linearized potential theories. Flutter results have been presented for two cases. A comparison with linear theory shows identical trends; however, the present calculations show a lower flutter frequency.

Bakhle, Milind A.↗

NASA NTRS1989

Flow of rarefied gases over two-dimensional bodies

A kinetic-theory analysis is made of the flow of rarefied gases over two-dimensional bodies of arbitrary curvature. The Boltzmann equation simplified by a model collision integral is written in an arbitrary orthogonal curvilinear coordinate system, and solved by means of finite-difference approximation with the discrete ordinate method. A numerical code is developed which can be applied to any two-dimensional submerged body of arbitrary curvature for the flow regimes from free-molecular to slip at transonic Mach numbers. Predictions are made for the case of a right circular cylinder.

Jeng, Duen-Ren↗

NASA NTRS1988

Applications of an exponential finite difference technique

An exponential finite difference scheme first presented by Bhattacharya for one dimensional unsteady heat conduction problems in Cartesian coordinates was extended. The finite difference algorithm developed was used to solve the unsteady diffusion equation in one dimensional cylindrical coordinates and was applied to two and three dimensional conduction problems in Cartesian coordinates. Heat conduction involving variable thermal conductivity was also investigated. The method was used to solve nonlinear partial differential equations in one and two dimensional Cartesian coordinates. Predicted results are compared to exact solutions where available or to results obtained by other numerical methods.

Handschuh, Robert F.↗