Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transitional flow”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

A spectral multi-domain technique applied to high-speed chemically reacting flows

The first applications of a spectral multidomain method for viscous compressible flow is presented. The method imposes a global flux balance condition at the interface so that high-order continuity of the solution is preserved. The global flux balance is imposed in terms of a spectral integral of the discrete equations across adjoining domains. Since the discretized equations interior to each domain solved are uncoupled from each other, and since the interface relation has a block structure, the solution scheme can be adapted to the particular requirements of each subdomain. The spectral multidomain technique presented is well-suited for the multiple scales associated with the chemically reacting and transition flows in hypersonic research. A nonstaggered multidomain discretization is used for the chemically reacting flow calculation, and the first implementation of a staggered multidomain mesh is presented for accurately solving the stability equation for a viscous compressible fluid.

Macaraeg, Michele G.↗
NASA NTRS1980

Calculation of skin-friction coefficients for low Reynolds number turbulent boundary layer flows

An analysis is presented of the reliability of various generally accepted empirical expressions for the prediction of the skin-friction coefficient C/sub f/ of turbulent boundary layers at low Reynolds numbers in zero-pressure-gradient flows on a smooth flat plate. The skin-friction coefficients predicted from these expressions were compared to the skin-friction coefficients of experimental profiles that were determined from a graphical method formulated from the law of the wall. These expressions are found to predict values that are consistently different than those obtained from the graphical method over the range 600 Re/sub theta 2000. A curve-fitted empirical relationship was developed from the present data and yields a better estimated value of C/sub f/ in this range. The data, covering the range 200 Re/sub theta 7000, provide insight into the nature of transitional flows. They show that fully developed turbulent boundary layers occur at Reynolds numbers Re/sub theta/ down to 425. Below this level there appears to be a well-ordered evolutionary process from the laminar to the turbulent profiles. These profiles clearly display the development of the turbulent core region and the shrinking of the laminar sublayer with increasing values of Re/sub theta/.

Barr, P. K.↗
NASA NTRS2002

Effect of Favorable Pressure Gradients on Turbine Blade Pressure Surface Heat Transfer

Recent measurements on a turbine rotor showed significant relaminarization effects. These effects were evident on the pressure surface heat transfer measurements. The character of the heat transfer varied with Reynolds number. Data were obtained for exit Reynolds numbers between 500,000 and 880,000. Tests were done with a high level of inlet turbulence, 7.5%. At lower Reynolds numbers the heat transfer was similar to that for laminar flow, but at a level higher than for laminar flow. At higher Reynolds numbers the heat transfer was similar to turbulent flow, when the acceleration parameter, K, was sufficiently small. The proposed paper discusses the experimental results, and also discusses approaches to calculating the surface heat transfer for the blade surface. Calculations were done using a three-dimensional Navier-Stokes CFD analysis. The results of these tests, when compared with previous blade tests in the same facility, illustrate modeling difficulties that were encountered in CFD predictions. The two blades were in many ways similar. However, the degree of agreement between the same analysis and the experimental data was significantly different. These differences are highlighted to illustrate where improvements in modeling approaches are needed for transitional flows.

Boyle, Robert J.↗
NASA NTRS1992

Crystallization, sublimation, and gas release in the interior of a porous comet nucleus

A numerical code is developed for evolutionary calculations of the thermal structure and composition of a porous comet nucleus made of water ice, in amorphous or crystalline form, other volatiles, dust, and gases trapped in amorphous ice. Bulk evaporation, crystallization, gas release, and free (Knudsen) flow of gases through the pores are taken into account. The numerical scheme yields exact conservation laws for mass and energy. The code is used to study the effect of bulk evaporation of ice in the interior of a comet nucleus during crystallization. It is found that evaporation controls the temperature distribution; the vapor prevents cooling of the crystallized layer of ice, by recondensation and release of latent heat. Thus high temperatures are maintained below the surface of the nucleus and down to depths of tens or hundreds of meters, even at large heliocentric distances, as long as crystallization goes on. Gas trapped in the ice and released during the phase transition flows both toward the interior and toward the surface and out of the nucleus. The progress of crystallization is largely determined by the contribution of gas fluxes to heat transfer.

Prialnik, Dina↗
NASA NTRSDate not supplied

Identification and Study of Validation Level Test Cases for Computational Modeling of Non-Charring Ablators

Computational modeling of Thermal Protection System (TPS) materials, used for aerospace applications, provides numerous advantages in preliminary selection and design of a heatshield material and shape for atmospheric entry vehicles. However, to serve as a reliable tool for prediction of material thermal and ablative behavior, the modeling approach needs to be validated against real experimental and flight data, preferably at a range of applied conditions. The validation study is typically very complex as it requires reliable measured data not only for the material thermal response and surface recession, but also well characterized environmental conditions. The validation problem becomes even more complex when the material thermal response is dictated by multi-physics effects such as solid conduction, in-depth thermal decomposition, pyrolysis gas flow and chemical reactions. The multi-physics effects complicate not only the modeling effort, but also the experimental measurement for validation of various aspects of the highly coupled problem. In this study, an attempt is made to identify suitable experimental data that could serve as a source for validation of material thermal response modeling tools. To reduce the computational complexity, this study focuses only on non-charring ablators, where the material thermal response could be modeled with a single governing equation for solid conduction and the ablation is limited only to the surface of the material. With a well characterized and publicly available experimental data being sparse, the study is limited in presenting test cases for only three materials: camphor, graphite and FiberForm® in the sequence of increased modeling complexity. Graphite is a commonly used TPS material for aerospace applications, both for leading edges of high-speed vehicles and internal insulation of solid rocket motors. FiberForm® is a porous carbon pre-form used in preparation of the well known PICA material Tran et al. [1996]. Inclusion of camphor into the list is conditioned with the relative simplicity in modeling the material thermal and chemical response and the low-enthalpy flow environment. In addition, camphor has been used as a simple test material for study of flow transition behavior by Stock and Ginoux [1973] and assessment of a heatshield shape change at flight relevant conditions by Rotondi et al. [2022]. In this work, the identified experimental data was extracted from the public literature and test cases that yet have been published. As it was found from the review, not a single test case contains an exhaustive set of data that would validate every aspect of the material physics. However, in the data collected, various aspects of the material behavior can be still validated, such as surface and in-depth temperature, amount of recession and a shape change. The identified experimental data for each case is accompanied with a characterized flow environment and simulated boundary conditions predicted by a Data-Parallel Line Relaxation (DPLR) code Wright et al. [1998]. In addition, material thermal response numerical simulations in each test case are performed with Kentucky Aerothermodynamics and Thermal Response System (KATS-MR) Zibitsker et al. [2022] providing a comparative study and a sanity check for the proposed validation data. Sample results from the performed numerical study are shown below. Figure 1 shows distribution of surface heat flux and pressure values on a hemi-cylinder model made of FiberForm® and tested in HyMETS arc-jet facility. The results are shown for the high pressure condition among the two tests. Flow simulation was performed with DPLR code on a quarter of original geometry. In the figure, the quarter shape was mirrored across zx and xy planes to show the complete distribution. Figure 2 shows the material response results for the high pressure case (7500 Pa), simulated with KATS-MR and a comparison to the experimental data for the surface temperature and shape shape. The simulation was performed on a 2-D slice, extracted in the xy plane at the middle of the sample. Figure 3 shows the material response simulation for the low pressure case (3500 Pa) and a comparison to the experimental data for surface temperature and shape change.

ablation↗
NASA NTRS1995

Direct numerical simulation of transitional and turbulent flow over a heated flat plate using finite-difference schemes

This report deals with the direct numerical simulation of transitional and turbulent flow at low Mach numbers using high-order-accurate finite-difference techniques. A computation of transition to turbulence of the spatially-evolving boundary layer on a heated flat plate in the presence of relatively high freestream turbulence was performed. The geometry and flow conditions were chosen to match earlier experiments. The development of the momentum and thermal boundary layers was documented. Velocity and temperature profiles, as well as distributions of skin friction, surface heat transfer rate, Reynolds shear stress, and turbulent heat flux, were shown to compare well with experiment. The results indicate that the essential features of the transition process have been captured. The numerical method used here can be applied to complex geometries in a straightforward manner.

Madavan, Nateri K.↗
NASA NTRS1994

Asymptotic behavior of solutions of the renormalization group K-epsilon turbulence model

Presently, the only efficient way to calculate turbulent flows in complex geometries of engineering interest is to use Reynolds-average Navier-Stokes (RANS) equations. As compared to the original Navier-Stokes problem, these RANS equations posses much more complicated nonlinear structure and may exhibit far more complex nonlinear behavior. In certain cases, the asymptotic behavior of such models can be studied analytically which, aside from being an interesting fundamental problem, is important for better understanding of the internal structure of the models as well as to improve their performances. The renormalization group (RNG) K-epsilon turbulence model, derived directly from the incompresible Navier-Stokes equations, is analyzed. It has already been used to calculate a variety of turbulent and transitional flows in complex geometries. For large values of the RNG viscosity parameter, the model may exhibit singular behavior. In the form of the RNG K-epsilon model that avoids the use of explicit wall functions, a = 1, so the RNG viscosity parameter must be smaller than 23.62 to avoid singularities.

Yakhot, A.↗
NASA NTRS1993

Burnett solutions along the stagnation line of a cooled cylinder in low-density hypersonic flows

The Burnett equations in the cylindrical coordinate system have been applied to the hypersonic flow over a cylinder in low density environments. Through a local similarity transformation, the Burnett equations are degenerated into nonlinear third-order ordinary differential equations which are only valid along the stagnation line. These equations are solved by the multi-staged Runge-Kutta integration technique. Steady state solutions are asymptotically reached by the time-marching procedure. A newly-developed computer code is used to solve both the Navier-Stokes and the Burnett equations within the shock layer. The computed heat transfer coefficients are consistent with wind-tunnel and STS flight data. Comparisons with the DSMC solutions and the Boltzmann solutions show that the Burnett equations are more appropriate as the governing equations than the Navier-Stokes equations in the transitional flow regime.

Liaw, G. S.↗
NASA NTRS1983

Stability and Interaction of Coherent Structure in Supersonic Reactive Wakes

A theoretical formulation and analysis is presented for a study of the stability and interaction of coherent structure in reacting free shear layers. The physical problem under investigation is a premixed hydrogen-oxygen reacting shear layer in the wake of a thin flat plate. The coherent structure is modeled as a periodic disturbance and its stability is determined by the application of linearized hydrodynamic stability theory which results in a generalized eigenvalue problem for reactive flows. Detailed stability analysis of the reactive wake for neutral, symmetrical and antisymmetrical disturbance is presented. Reactive stability criteria is shown to be quite different from classical non-reactive stability. The interaction between the mean flow, coherent structure and fine-scale turbulence is theoretically formulated using the von-Kaman integral technique. Both time-averaging and conditional phase averaging are necessary to separate the three types of motion. The resulting integro-differential equations can then be solved subject to initial conditions with appropriate shape functions. In the laminar flow transition region of interest, the spatial interaction between the mean motion and coherent structure is calculated for both non-reactive and reactive conditions and compared with experimental data wherever available. The fine-scale turbulent motion determined by the application of integral analysis to the fluctuation equations. Since at present this turbulence model is still untested, turbulence is modeled in the interaction problem by a simple algebraic eddy viscosity model. The applicability of the integral turbulence model formulated here is studied parametrically by integrating these equations for the simple case of self-similar mean motion with assumed shape functions. The effect of the motion of the coherent structure is studied and very good agreement is obtained with previous experimental and theoretical works for non-reactive flow. For the reactive case, lack of experimental data made direct comparison difficult. It was determined that the growth rate of the disturbance amplitude is lower for reactive case. The results indicate that the reactive flow stability is in qualitative agreement with experimental observation.

Menon, Suresh↗
NASA NTRS1994

Aerothermodynamics of Pyrolizing Surfaces in Hypersonic Rarefied Flows

Direct simulation Monte Carlo (DSMC) calculations of rarefied flows about entry bodies typically employ a fixed surface temperature or a radiative-equilibrium energy balance to compute that temperature. Such boundary conditions neglect any effects of heat capacitance and heat conduction in the spacecraft heat shield and, therefore, provide an upper bound for the surface temperature. Such calculations also neglect pyrolysis from the heat shield which can be significant for a high-energy incident flow at very low densities. Accurate prediction of both heating and aerodynamic forces requires including pyrolysis and surface heat transfer in the models for gas-surface interaction employed in DSMC methods. Although these physical models have long appeared in various continuum flow calculation codes, they have only recently appeared in DSMC codes which are required to simulate rarefied flows during entry at high altitudes. In the current implementation, routines from the widely distributed Charring Material Thermal Response and Ablation (CMA) program are coupled into a DSMC code to calculate the one-dimensional heat transfer into the carbon phenolic heat shield at each point on a vehicle surface. Temperature-dependent material properties, surface re-radiation, and in-depth pyrolysis were included in the calculation, but surface ablation was neglected. Sample calculations for entry of the Galileo probe into the atmosphere of Jupiter demonstrate that including pyrolysis in the model leads to significant differences in predicted aerodynamics. Granted, the drag coefficient does not depend strongly on the surface temperature which can itself be significantly below the radiative equilibrium value during entry. However, the surface mass flux due to pyrolysis of the material is significant once the probe drops to altitudes characterized by transition flow. This leads to a noticeable increase in drag and a decrease in heating compared to a body without pyrolysis.

Haas, Brian L.↗
NASA NTRS1992

Transition of oscillatory flow in tubes - An empirical model for application to Stirling engines

Attention is given to an empirical model for transition to turbulence in oscillatory flows in straight tubes. Designed after a correlation for transition of a boundary layer on a flat plate, the model yields the laminar flow momentum thickness Reynolds number that must be met before transition to turbulence will occur. The transition point is located by comparing this to the actual momentum thickness Reynolds number. A scheme is proposed for estimating the momentum thickness Reynolds number in terms of the position within the cycle, the maximum value of the diameter Reynolds within the cycle, Re(max), and the dimensionless frequency, Valensi number. Results from an experimental study of oscillatory flow in a tube are employed to develop the model. When the flow is determined to be turbulent, it is proposed that a fully-developed, steady flow friction coefficient be applied. When the flow is laminar, the assumption of fully developed flow cannot be made; thus, a method is suggested for estimating the friction factor.

Simon, T. W.↗
NASA NTRS1972

Spreading of a turbulent disturbance.

The presented paper shows the effect of local Mach number on the turbulent disturbance spreading angle relative to the wall and on lateral disturbance spreading angles. Almost all the disturbances angles relative to the wall were determined from investigations where hot-wire contours or hot-film surveys of a 'laminar' boundary layer were obtained. Lateral disturbance spreading angles were obtained from investigations of various conditions including turbulent bursts, reported observations of transverse contamination, and observed transitional flow. It is noted that the disturbance spreading angle relative to the wall seems to remain essentially invariant with Mach number, while the lateral spreading angle decreases sharply with increasing Mach number up to about 6. The good agreement between lateral disturbance spreading angle data and data for the variation of turbulent jet spreading angle with Mach number implies that in the lateral dimension, turbulence in a boundary layer may develop essentially free of wall constraints (similar to a free shear layer).

Fischer, M. C.↗
NASA NTRS1976

Molecular beam mass spectrometer development

An analytical model, based on the kinetics theory of a drifting Maxwellian gas is used to determine the nonequilibrium molecular density distribution within a hemispherical shell open aft with its axis parallel to its velocity. The concept of a molecular shield in terrestrial orbit above 200 km is also analyzed using the kinetic theory of a drifting Maxwellian gas. Data are presented for the components of the gas density within the shield due to the free stream atmosphere, outgassing from the shield and enclosed experiments, and atmospheric gas scattered off a shield orbiter system. A description is given of a FORTRAN program for computating the three dimensional transition flow regime past the space shuttle orbiter that employs the Monte Carlo simulation method to model real flow by some thousands of simulated molecules.

Brock, F. J.↗
NASA NTRS1976

Discharge stabilization studies of CO laser gas mixtures in quasi-steady supersonic flow

Experiments were conducted to study the applicability of a double discharge stabilization scheme in conditions appropriate for high energy CO lasers in supersonic flows. A Ludwieg tube impulse flow facility and a ballasted capacitor bank provided essentially steady flow and discharge conditions (d.c.) for times longer than ten electrode length-flow transit times. Steady, arc-free, volume discharges were produced in a Mach 3 test cavity using an auxiliary discharge to stabilize the main discharge in N2 and He/CO mixture. A signigicant result is the lack of observed plasma E/N changes in response to auxiliary discharge current changes. Also, where glow discharges were obtained, the energy loading achieved was very much less than the threshold level required for laser operation.

Srinivasan, G.↗