Investigation of the Flow Through a Single-stage Two-dimensional Nozzle in the Langley 11-inch Hypersonic Tunnel
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
The present paper presents a numerical procedure to calculate supersonic and hypersonic flows using the compressible law of the wall. The turbulence models under consideration include the Launder-Reece-Rodi-Gibson Reynolds-stress model and the k-epsilon model. The models coupled with the proposed wall function technique have been tested in both separated and unseparated flows. The flows include (1) an insulated flat plate flow over a range of Mach numbers, (2) a Mach 5 flat plate flow with cold wall conditions, (3) a two dimensional supersonic compression corner flow, (4) a hypersonic flow over an axisymmetric flare, and (5) a hypersonic flow over a 2-D compression corner. Results indicate that the wall function technique gives improved predictions of skin friction and heat transfer in separated flows compared with models using wall dampers. Predictions of the extent of separation are not improved over the wall damper models except with the Reynolds-stress model for the supersonic compression corner flow case.
A computational study of the separated flow through a 2-D asymmetric subsonic diffuser has been performed. The Wind Computational Fluid Dynamics code is used to predict the separation and reattachment behavior for an incompressible diffuser flow. The diffuser inlet flow is a two-dimensional, turbulent, and fully-developed channel flow with a Reynolds number of 20,000 based on the centerline velocity and the channel height. Wind solutions computed with the Menter SST, Chien k-epsilon, Spalart-Allmaras and Explicit Algebraic Reynolds Stress turbulence models are compared with experimentally measured velocity profiles and skin friction along the upper and lower walls. In addition to the turbulence model study, the effects of grid resolution and use of wall functions were investigated. The grid studies varied the number of grid points across the diffuser and varied the initial wall spacing from y(sup +) = 0.2 to 60. The wall function study assessed the applicability of wall functions for analysis of separated flow. The SST and Explicit Algebraic Stress models provide the best agreement with experimental data, and it is recommended wall functions should only be used with a high level of caution.
The objective was to study diamond film growth and its properties in order to enhance the laser damage threshold of substrate materials. Calculations were performed to evaluate laser induced thermal stress parameter, R(sub T) of diamond. It is found that diamond has several orders of magnitude higher in value for R(sub T) compared to other materials. Thus, the laser induced damage threshold (LIDT) of diamond is much higher. Diamond films were grown using a microwave plasma enhanced chemical vapor deposition (MPECVD) system at various conditions of gas composition, pressure, temperature, and substrate materials. A 0.5 percent CH4 in H2 at 20 torr were ideal conditions for growing of high quality diamond films on substrates maintained at 900 C. The diamond films were polycrystalline which were characterized by scanning electron microscopy (SEM) and Raman scattering spectroscopy. The top surface of the growing film is always rough due to the facets of polycrystalline film while the back surface of the film replicates the substrate surface. An analytical model based on two dimensional periodic heat flow was developed to calculate the effective in-plane (face parallel) diffusivity of a two layer system. The effective diffusivity of diamond/silicon samples was measured using a laser pulse technique. The thermal conductivity of the films was measured to be 13.5 W/cm K, which is better than that of a type Ia natural diamond. Laser induced damage experiments were performed on bare Si substrates, diamond film coated Si, and diamond film windows. Significant improvements in the LIDT were obtained for diamond film coated Si compared to the bare Si.
The laminar base flow field of a two-dimensional reentry body has been studied by Telenin's method. The flow domain was divided into strips along the x-axis, and the flow variations were represented by Lagrange interpolation polynomials in the transformed vertical coordinate. The complete Navier-Stokes equations were used in the near wake region, and the boundary layer equations were applied elsewhere. The boundary conditions consisted of the flat plate thermal boundary layer in the forebody region and the near wake profile in the downstream region. The resulting two-point boundary value problem of 33 ordinary differential equations was then solved by the multiple shooting method. The detailed flow field and thermal environment in the base region are presented in the form of temperature contours, Mach number contours, velocity vectors, pressure distributions, and heat transfer coefficients on the base surface. The maximum heating rate was found on the centerline, and the two-dimensional stagnation point flow solution was adquate to estimate the maximum heating rate so long as the local Reynolds number could be obtained.
A new method has been developed for two- and three-dimensional computations of viscous supersonic flows with embedded subsonic regions adjacent to solid boundaries. The approach employs a reduced form of the Navier-Stokes equations which allows solution as an initial-boundary value problem in space, using an efficient noniterative forward marching algorithm. Numerical instability associated with forward marching algorithms for flows with embedded subsonic regions is avoided by approximation of the reduced form of the Navier-Stokes equations in the subsonic regions of the boundary layers. Supersonic and subsonic portions of the flow field are simultaneously calculated by a consistently split linearized block implicit computational algorithm. The results of computations for a series of test cases relevant to internal supersonic flow is presented and compared with data. Comparison between data and computation are in general excellent thus indicating that the computational technique has great promise as a tool for calculating supersonic flow with embedded subsonic regions. Finally, a User's Manual is presented for the computer code used to perform the calculations.
The complexity of codes for solutions of real aerodynamic problems has progressed from simple two-dimensional models to three-dimensional inviscid and viscous models. As the algorithms used in the codes increased in accuracy, speed and robustness, the codes were steadily incorporated into standard design processes. The highly sophisticated codes, which provide solutions to the truly complex flows, require computers with large memory and high computational speed. The advent of high-speed supercomputers, such that the solutions of these complex flows become more practical, permits the introduction of the codes into the design system at an earlier stage. The results of several codes which either were already introduced into the design process or are rapidly in the process of becoming so, are presented. The codes fall into the area of turbomachinery aerodynamics and hypersonic propulsion. In the former category, results are presented for three-dimensional inviscid and viscous flows through nozzle and unducted fan bladerows. In the latter category, results are presented for two-dimensional inviscid and viscous flows for hypersonic vehicle forebodies and engine inlets.
Nonradial hydrodynamic flow can be generated or amplified during plasma compression by various mechanisms, including the compression itself. In certain circumstances, the plasma may reach a viscous state; for example, in compression experiments seeking fusion, the fuel plasma may reach a viscous state late in the compression due in part to the rising fuel temperature. Here, we consider viscous dissipation of nonradial flow in the case of initially isotropic, three-dimensional (3D), turbulent flow fields compressed at constant velocity in two dimensions. Prior work in the case of 3D compressions has shown the possibility of effective viscous dissipation of nonradial flow under compression. We show that, theoretically, complete viscous dissipation of the nonradial flow should still occur in the 2D case when the plasma heating is adiabatic and the viscosity has the (strong) Braginskii temperature dependence (μ~ T 5/2 ). Futhermore, in the general case, the amount of compression required is very large even for modest initial Reynolds numbers, with the compression reaching an intermediate state dominated by variations only in the noncompressed direction. We show that both the nonlinearity and boundary conditions can play important roles in setting the characteristics and ease of the viscous dissipation.
The prediction of boundary layer deceleration and separation in regions of adverse pressure gradients are discussed as an unsolved problem of fundamental practical importance in aerodynamics. A three layered theoretical flow model of boundary layer separation including viscous-inviscid interaction was conceived for the case of subsonic two dimensional steady laminar flow.
A flow model has been devised to deal with the viscid-inviscid interaction of a class of two-dimensional incompressible separated flow problems. It is suggested that the corresponding inviscid flow of these problems is described by the free streamline theory with few unspecified parameters and their values are, in turn, determined by the viscous flow considerations. The problem of a flow past a backward facing step is selected for study in detail. The viscous flow components of turbulent jet mixing, recompression, and reattachment are delineated and studied individually. When they are later combined, it is found that the point of reattachment behaves as a saddle-point-type singularity in the system of differential equations describing the viscous flow process. This feature is employed to the determination of the aforementioned free parameters and thus the establishment of the overall corresponding inviscid flow field. The resulting base pressure coefficient for the specific case agrees reasonably well with the available experimental data.
A variety of computational methods were developed to obtain shockless or near shockless flow past two-dimensional airfoils. The approach used was the method of complex characteristics, which determines smooth solutions to the transonic flow equations based on an input speed distribution. General results from fluid mechanics are presented. An account of the method of complex characteristics is given including a description of the particular spaces and coordinates, conformal transformations, and numerical procedures that are used. The operation of the computer program COMPRES is presented along with examples of blade sections designed with the code. A user manual is included with a glossary to provide additional information which may be helpful. The computer program in Fortran, including numerous comment cards is listed.
Unsteady two- and three-dimensional flow structure at leading and trailing edges of bodies can be characterized effectively using recently developed techniques for acquisition and interpretation of flow visualization. The techniques addressed here include: flow image/surface pressure correlations; 3-D reconstruction of flow structure from flow images; and interactive interpretation of flow images with theoretical simulations. These techniques can be employed in conjunction with: visual correlation and ensemble-averaging, both within a given image and between images; recognition of patterns from images; and estimates of velocity eigenfunctions from images.
New error estimation techniques are described for two-dimensional convection-dominated flows. Examples of applications to linear convection problems are given. The methods are based on discontinuous finite-element techniques. Results show that conventional interpolation methods may lead to gross underestimates of actual error and be inadequate for use in adaptive methods.
The dynamic behavior of the vapor flow in heat pipes is investigated at startup and during operational transients. The vapor is modeled as two-dimensional, compressible viscous flow in an enclosure with inflow and outflow boundary conditions. For steady-state and operating transients, the SIMPLER method is used. In this method a control volume approach is employed on a staggered grid which makes the scheme very stable. It is shown that for relatively low input heat fluxes the compressibility of the vapor flow is low and the SIMPLER scheme is suitable for the study of transient vapor dynamics. When the input heat flux is high or the process under a startup operation starts at very low pressures and temperatures, the vapor is highly compressible and a shock wave is created in the evaporator. It is shown that for a wide range of input heat fluxes, the standard methods, including the SIMPLER scheme, are not suitable. A nonlinear filtering technique, along with the centered difference scheme, are then used for shock capturing as well as for the solution of the cell Reynolds-number problem. For high heat flux, the startup transient phase involves multiple shock reflections in the evaporator region. Each shock reflection causes a significant increase in the local pressure and a large pressure drop along the heat pipe. Furthermore, shock reflections cause flow reversal in the evaporation region and flow circulations in the adiabatic region. The maximum and maximum-averaged pressure drops in different sections of the heat pipe oscillate periodically with time because of multiple shock reflections. The pressure drop converges to a constant value at steady state. However, it is significantly higher than its steady-state value at the initiation of the startup transient. The time for the vapor core to reach steady-state condition depends on the input heat flux, the heat pipe geometry, the working fluid, and the condenser conditions. However, the vapor transient time, for an Na-filled heat pipe is on the order of seconds. Depending on the time constant for the overall system, the vapor transient time may be very short. Therefore, the vapor core may be assumed to be quasi-steady in the transient analysis of a heat pipe operation.
Two-step flow control adjusts to dimensional changes in purged equipment. Two pressure actuated valves provide two-step control. When turbopump starts, helium flows through both valves. As pump operates its parts expand, clearance decrease and back pressure builds up on valves. At preset back pressure value, one valve closes, reducing helium flow.
Forced and naturally occurring, self-sustaining oscillations of transonic flows in two-dimensional diffusers were computed using MacCormack's hybrid method. Depending upon the shock strengths and the area ratios, the flow was fully attached or separated by either the shock or the adverse pressure gradient associated with the enlarging diffuser area. In the case of forced oscillations, a sinusoidal plane pressure wave at frequency 300 Hz was prescribed at the exit. A sufficiently large amount of data were acquired and Fourier analyzed. The distrbutions of time-mean pressures, the power spectral density, and the amplitude with phase angle along the top wall and in the core region were determined. Comparison with experimental results for the forced oscillation generally gave very good agreement; some success was achieved for the case of self-sustaining oscillation despite substantial three-dimensionality in the test. An observation of the sequence of self-sustaining oscillations was given.
Methods are considered for reducing a non-random fatigue loading history to a concise description and then for reconstructing a time history similar to the original. In particular, three methods of reconstruction based on a rain-flow cycle counting matrix are presented. A rain-flow matrix consists of the numbers of cycles at various peak and valley combinations. Two methods are based on a two dimensional rain-flow matrix, and the third on a three dimensional rain-flow matrix. Histories reconstructed by any of these methods produce a rain-flow matrix identical to that of the original history, and as a result the resulting time history is expected to produce a fatigue life similar to that for the original. The procedures described allow lengthy loading histories to be stored in compact form.
Report discusses application of lower-upper (LU) implicit numerical-integration scheme, together with interactive grid-generation method, to simulation of inviscid, compressible flow in two-dimensional model of turbine cascade. Study part of continuing effort to develop techniques for computation of complicated three-dimensional flows in turbomachinery.