Analysis of the transient flow of hydrogen in heated passages in the laminar-turbulent transition region
Hydrogen transient flow analysis in heated passages in laminar-turbulent transition region
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.
Hydrogen transient flow analysis in heated passages in laminar-turbulent transition region
A prediction of whether a point on a computer-generated surface is adjacent to laminar or turbulent flow is made using a transition prediction technique. A plurality of boundary-layer properties at the point are obtained from a steady-state solution of a fluid flow in a region adjacent to the point. A plurality of instability modes are obtained, each defined by one or more mode parameters. A vector of regressor weights is obtained for the known instability growth rates in a training dataset. For each instability mode in the plurality of instability modes, a covariance vector is determined, which is the covariance of a predicted local growth rate with the known instability growth rates. Each covariance vector is used with the vector of regressor weights to determine a predicted local growth rate at the point. Based on the predicted local growth rates, an n-factor envelope at the point is determined.
Computations were completed of transition curves in the conventional annulus, including hysteresis effect. The model GEOSIM was used to compute the transition between axisymmetric flow and baroclinic wave flow in the conventional annulus experiments. Thorough testing and documentation of the GEOSIM code were also completed. The Spacelab 3 results from the Geophysical Fluid Flow Cell (GFFC) were reviewed and numerical modeling was performed of many of the cases with horizontal temperature gradients as well as heating from below, with different rates of rotation. A numerical study of the lower transition to axisymmetric flow in the baroclinic annulus was performed using GEOSIM.
The problem of one-dimensional gas flow through a sinusoidal gravitational potential with a series of equally spaced shock fronts is considered. For time-independent and spatially periodic flows, an integral equation relates the flow velocity to the gravitational potential and source functions of energy and momentum densities. It is suggested that this problem simulates some of the dynamical effects of the azimuthal flow on a parcel of gas at a fixed radius in a galactic disk that supports a spiral density wave structure. It omits the radial motions. In this context, a typical azimuthal velocity across a spiral arm is specified and the resultant source functions are obtained. The implications of this calculation suggest heating of the gaseous disk by O and B stars behind the shock front, and subsequent interarm cooling is a plausible mechanism for maintaining the periodic cycling of the gaseous component of the galactic disk.
Two-phase gas-liquid flows occur in a wide variety of situations. In addition to normal gravity applications, such flows may occur in space operations such as active thermal control systems, power cycles, and storage and transfer of cryogenic fluids. Various flow patterns exhibiting characteristic spatial and temporal distribution of the two phases are observed in two-phase flows. The magnitude and orientation of gravity with respect to the flow has a strong impact on the flow patterns observed and on their boundaries. The identification of the flow pattern of a flow is somewhat subjective. The same two-phase flow (especially near a flow pattern transition boundary) may be categorized differently by different researchers. Two-phase flow patterns are somewhat simplified in microgravity, where only three flow patterns (bubble, slug and annular) have been observed. Annular flow is obtained for a wide range of gas and liquid flow rates, and it is expected to occur in many situations under microgravity conditions. Slug flow needs to be avoided, because vibrations caused by slugs result in unwanted accelerations. Therefore, it is important to be able to accurately predict the flow pattern which exists under given operating conditions. It is known that the wavy liquid film in annular flow has a profound influence on the transfer of momentum and heat between the phases. Thus, an understanding of the characteristics of the wavy film is essential for developing accurate correlations. In this work, we review our recent results on flow pattern transitions and wavy films in microgravity.
Some of the recent work at NASA Langley on transition, turbulence, and reacting flows is summarized. Much of this effort has been motivated by outstanding technological problems in high-speed flow. A class of numerical algorithms for these problems has been developed and a variety of physical problems have been simulated. Descriptions are provided of the basic mathematical models, the nature of the numerical methods, and some of the recent simulations.
The ongoing effort at NASA Langley Research Center aimed at developing transitional zone models for Reynolds-averaged Navier-Stokes calculations is described. Results from three models that we have examined - two zero-equation models and a two-equation model of the k-omega type - are presented for high-speed flows. The use of analytical and computational tools to provide guidance and to provide the transitional flow-field data bases necessary to develop such models is discussed. Some results from validation of these tools are given; qualitative comparisons between modeled and theoretical/computational representations of the transitional zone are provided.
A structural performance and resizing finite element thermal analysis computer program was used in the reentry heat transfer analysis of the space shuttle. Two typical wing cross sections and a midfuselage cross section were selected for the analysis. The surface heat inputs to the thermal models were obtained from aerodynamic heating analyses, which assumed a purely turbulent boundary layer, a purely laminar boundary layer, separated flow, and transition from laminar to turbulent flow. The effect of internal radiation was found to be quite significant. With the effect of the internal radiation considered, the wing lower skin temperature became about 39 C (70 F) lower. The results were compared with fight data for space transportation system, trajectory 1. The calculated and measured temperatures compared well for the wing if laminar flow was assumed for the lower surface and bay one upper surface and if separated flow was assumed for the upper surfaces of bays other than bay one. For the fuselage, good agreement between the calculated and measured data was obtained if laminar flow was assumed for the bottom surface. The structural temperatures were found to reach their peak values shortly before touchdown. In addition, the finite element solutions were compared with those obtained from the conventional finite difference solutions.
The 'classical' polar wind is an ambipolar outflow of thermal plasma from the terrestrial ionosphere at high latitudes. As the plasma escapes along diverging geomagnetic flux tubes, it undergoes four major transitions, including a transition from chemical to diffusion dominance, a transition from subsonic to supersonic flow, a transition from collision-dominated to collisionless regimes, and a transition from a heavy to a light ion. A further complication arises because of horizontal convection of the flux tubes owing to magnetospheric electric fields. Recent modelling predictions indicate that the polar wind has the following characteristics: (1) the ion and electron distributions are anisotropic and asymmetric in the collisionless regime; (2) elevated electron temperatures act to produce significant escape fluxes of suprathermal O(+) ions; (3) the interaction of the hot magnetospheric and cold ionospheric electron population leads to a localized (double layer) electric field which accelerates the polar wind ions; (4) a time-dependent expansion produces suprathermal ions; and (5) large perturbations lead to the formation of forward and reverse shocks. These and other results are reviewed.
Space Shuttle preflight hypergol oxidizer and fuel loading were monitored using a nonintrusive flow measurement system (NFMS) during preflight operation. A pair of 4-megahertz ultrasonic transducers measured the flow rate using a transit-time flow measurement technique. Using wireless cellular phone technology, flow-rate data was monitored and recorded. Excellent correlation was observed between the flow profiles measured using the proposed nonintrusive ultrasonic flowmeters (UFM's) and those made from conventional intrusive turbine meters (TFM). Based on the preliminary tests, it is concluded that the nonintrusive method of flow measurement has the same or higher accuracy, is simpler, and costs less than the existing TFM. Benefits of UFM's include a highly flexible, cost-effective, reliable, hazard free, and streamlined hypergol operation. Redundant installation of ultrasonic flowmeters was recommended for a few more launches prior to the replacement of the existing TFM.
A prediction of whether a point on a computer-generated surface is adjacent to laminar or turbulent flow is made using a transition prediction technique. A plurality of instability modes are obtained, each defined by one or more mode parameters. A vector of regressor weights is obtained for the known instability growth rates in a training dataset. For an instability mode in the plurality of instability modes, a covariance vector is determined. A predicted local instability growth rate at the point is determined using the covariance vector and the vector of regressor weights. Based on the predicted local instability growth rate, an n-factor envelope at the point is determined.
An experiment on a 45-degree swept wing was conducted to study three-dimensional boundary-layer characteristics using surface-mounted, micro-thin, multi-element hot-film sensors. Cross-flow vortex structure and boundary-layer transition were measured from the simultaneously acquired signals of the hot films. Spanwise variation of the root-mean-square (RMS) hot-film signal show a local minima and maxima. The distance between two minima corresponds to the stationary cross-flow vortex wavelength and agrees with naphthalene flow-visualization results. The chordwise and spanwise variation of amplified traveling (nonstationary) cross-flow disturbance characteristics were measured as Reynolds number was varied. The frequency of the most amplified cross-flow disturbances agrees with linear stability theory.
This paper describes the creation, by temporal direct numerical simulation and the analysis based on the Reynolds stress transport equations, of a high-quality data set that represents the laminar-turbulent transition of a high-speed boundary-layer flow. Following Pruett and Zang (1992), and with the help of algorithmic refinements, the evolution of an axial, Mach 4.5 boundary-layer flow along a hollow cylinder is simulated numerically. Favre-averaged Reynolds stress transport equations are derived in generalized curvilinear coordinates and are then specialized to the cylindrical geometry at hand. Reynolds stresses and various turbulence quantities, such as turbulent kinetic energy and turbulent Mach number, are calculated from the numerical data at various stages of the transition process. The kinetic energy 'budgets' are constructed from the transport equations. Various contributing terms for the evolution of kinetic energy, like the rates of production and dissipation, transport, and diffusion, are presented. The compressible dissipation rate is small in comparison with the solenoidal dissipation rate for all times. The pressure-dilatation term is of the same order of magnitude as the compressible dissipation rate.
The possibility of global-scale transitions between the Hadley and Rossby atmospheric regimes is investiated using a simple three-dimensional rotating spherical model without a boundary layer structure, bottom topography, or cumulus friction, and the expected occurrence of the Hadley to Rossby transition is demonstrated. It is shown that a transition from Hadley flow to wavenumber-5 Rossby flow is preferred, in agreement with standard baroclinic instability results. This result gives a reasonable Rossby wave bifurcation from the Hadley solution. For the cases examined, it was found that the upper symmetric Hadley regime does not exist and that the Hadley to Rossby transition depends on the values of the eddy viscosities.
The role of secondary instabilities in the transition to turbulence in three-dimensional boundary layers (especially in swept-wing flows) is characterized, reviewing the results of recent theoretical and experimental investigations. Consideration is given to cross-flow vortices, boundary-layer profile measurements, flow visualization, spanwise wavelength determinations, interactions, problems involving rotating disks and cones, and leading-edge contamination. All of the flows investigated are shown to exhibit streamwise vorticity and to depend strongly on initial conditions.
The stability of three-dimensional rotating disk flow is investigated, including the effects of Coriolis forces and streamline curvature. The numerical results show that the critical Reynolds number for establishment of stationary vortex flow is 287. These vortices spiral outward at an angle of about 11.2 deg, and transition to turbulence occurs when their total amplification is about e to the 11th. New experimental results are also reported on the spatial growth rates of the stationary vortices. It is shown that the analysis gives growth rates that compare much better with the experimental results than do results obtained using the Orr-Sommerfeld equation. The experimental results tend to support the numerical prediction that the number of stationary vortices varies directly with the Reynolds number. The calculations also indicate the existence of weakly unstable propagating (type II) modes at low Reynolds numbers (Critical Reynolds Number being approximately equal to 49).
Circular-to-rectangular transition ducts are used as exhaust system components of aircraft with rectangular exhaust nozzles. Often, the incoming flow of these transition ducts includes a swirling velocity component remaining from the gas turbine engine. Previous transition duct studies have either not included inlet swirl or when inlet swirl was considered, only overall performance parameters were evaluated. Circular-to-rectangular transition duct flows with and without inlet swirl were explored in order to understand the effect of inlet swirl on the transition duct flow field and to provide detailed duct flow data for comparison with numerical code predictions. A method was devised to create a swirling, solid body rotational flow with minimal associated disturbances. Coefficients based on velocities and total and static pressures measured incross stream planes at four axial locations within the transition duct, along with surface static pressure measurements and surface oil film visualization, are presented for both nonswirling and swirling incoming flow. In both cases the inlet centerline Mach number was 0.35. The Reynolds number based on the inlet centerline velocity and duct inlet diameter was 1,547,000 for nonswirling and 1,366,000 for swirling flow. The maximum swirl angle was 15.6 deg. Two pair of counter-rotating side wall vortices appeared in the duct flow without inlet swirl. These vortices were absent in the swirling incoming flow cases.
A low Reynolds number k-epsilon turbulence model and conditioned momentum, energy and turbulence equations were used to predict bypass transition heat transfer on a flat plate in a high-disturbance environment with zero pressure gradient. The use of conditioned equations was demonstrated to be an improvement over the use of the global-time-averaged equations for the calculation of velocity profiles and turbulence intensity profiles in the transition region of a boundary layer. The approach of conditioned equations is extended to include heat transfer and a modeling of transition events is used to predict transition onset and the extent of transition on a flat plate. The events, which describe the boundary layer at the leading edge, result in boundary-layer regions consisting of: (1) the laminar, (2) pseudolaminar, (3) transitional, and (4) turbulent boundary layers. The modeled transition events were incorporated into the TEXSTAN 2-D boundary-layer code which is used to numerically predict the heat transfer. The numerical predictions in general compared well with the experimental data and revealed areas where additional experimental information is needed.