Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shear Layers”

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 397 records · Page 22

Analysis of Nozzle Jet Plume Effects on Sonic Boom Signature

An axisymmetric full Navier-Stokes computational fluid dynamics (CFD) study was conducted to examine nozzle exhaust jet plume effects on the sonic boom signature of a supersonic aircraft. A simplified axisymmetric nozzle geometry, representative of the nozzle on the NASA Dryden NF-15B Lift and Nozzle Change Effects on Tail Shock (LaNCETS) research airplane, was considered. The highly underexpanded nozzle flow is found to provide significantly more reduction in the tail shock strength in the sonic boom N-wave pressure signature than perfectly expanded and overexpanded nozzle flows. A tail shock train in the sonic boom signature, similar to what was observed in the LaNCETS flight data, is observed for the highly underexpanded nozzle flow. The CFD results provide a detailed description of the nozzle flow physics involved in the LaNCETS nozzle at different nozzle expansion conditions and help in interpreting LaNCETS flight data as well as in the eventual CFD analysis of a full LaNCETS aircraft. The current study also provided important information on proper modeling of the LaNCETS aircraft nozzle. The primary objective of the current CFD research effort was to support the LaNCETS flight research data analysis effort by studying the detailed nozzle exhaust jet plume s imperfect expansion effects on the sonic boom signature of a supersonic aircraft. Figure 1 illustrates the primary flow physics present in the interaction between the exhaust jet plume shock and the sonic boom coming off of an axisymmetric body in supersonic flight. The steeper tail shock from highly expanded jet plume reduces the dip of the sonic boom N-wave signature. A structured finite-volume compressible full Navier-Stokes CFD code was used in the current study. This approach is not limited by the simplifying assumptions inherent in previous sonic boom analysis efforts. Also, this study was the first known jet plume sonic boom CFD study in which the full viscous nozzle flow field was modeled, without coupling to a sonic boom propagation analysis code, from the stagnation chamber of the nozzle to the far field external flow, taking into account all nonisentropic effects in the shocks, boundary layers, and free shear layers, and their interactions at distances up to 30 times the nozzle exit diameter from the jet centerline. A CFD solution is shown in Figure 2. The flow field is very complicated and multi-dimensional, with shock-shock and shockplume interactions. At the time of this reporting, a full three-dimensional CFD study was being conducted to evaluate the effects of nozzle vectoring on the aircraft tail shock strength.

Bui, Trong↗

The instability of sheared liquid layers

The critical values of the Reynolds number (Rc) above which unstable travelling waves exist in a steady shear flow are calculated with the inclusion of a term in the normal stress boundary condition omitted by Miles (1960). Results show a value of Rc = 34.2, compared to Miles' value of 203, for the case of a zero-surface tension interface. The critical wavenumber, critical phase speed, and Rc are calculated as functions of a nondimensional surface tension, and the mechanism of the instability is investigated by examining the disturbance energy equation. A long wave instability is found at small values of the Reynolds number when the shear flow has a parabolic velocity profile. In addition, these results are extended to larger values of the wavenumber by numerical methods, and the relation between this long wave instability and the profile curvature is examined.

Smith, M. K.↗

A numerical study of wave propagation in a confined mixing layer by eigenfunction expansions

It is well known that the growth rate of instability waves of a two-dimensional free shear layer is reduced greatly at supersonic convective Mach numbers. In previous works, it has been shown that new wave modes exist when the shear layers are bounded by a channel due to the coupling effect between the acoustic wave modes and the motion of the mixing layer. The present work studies the simultaneous propagation of multiple stability waves using numerical simulation. It is shown here that the coexistence of two wave modes in the flow field can lead to an oscillatory growth of disturbance energy with each individual wave mode propagating linearly. This is particularly important when the growth rates of the unstable waves are small. It is also shown here that the propagation of two neutrally stable wave modes can lead to a stationary periodic structure of rms fluctuations. In the numerical simulations presented here the forced wave modes are propagating at same frequency, but with different phase velocities. In order to track the growth of each wave mode as it propagates downstream, a numerical method that can effectively detect and separate the contribution of the individual wave is given. It is demonstrated that by a least square fitting of the disturbance field with eigenfunctions the amplitude of each wave mode can be found. Satisfactory results as compared to linear theory are obtained.

Hu, Fang Q.↗

LES/RANS Simulation of a Supersonic Reacting Wall Jet

This work presents results from large-eddy / Reynolds-averaged Navier-Stokes (LES/RANS) simulations of the well-known Burrows-Kurkov supersonic reacting wall-jet experiment. Generally good agreement with experimental mole fraction, stagnation temperature, and Pitot pressure profiles is obtained for non-reactive mixing of the hydrogen jet with a non-vitiated air stream. A lifted flame, stabilized between 10 and 22 cm downstream of the hydrogen jet, is formed for hydrogen injected into a vitiated air stream. Flame stabilization occurs closer to the hydrogen injection location when a three-dimensional combustor geometry (with boundary layer development resolved on all walls) is considered. Volumetric expansion of the reactive shear layer is accompanied by the formation of large eddies which interact strongly with the reaction zone. Time averaged predictions of the reaction zone structure show an under-prediction of the peak water concentration and stagnation temperature, relative to experimental data and to results from a Reynolds-averaged Navier-Stokes calculation. If the experimental data can be considered as being accurate, this result indicates that the present LES/RANS method does not correctly capture the cascade of turbulence scales that should be resolvable on the present mesh. Instead, energy is concentrated in the very largest scales, which provide an over-mixing effect that excessively cools and strains the flame. Predictions improve with the use of a low-dissipation version of the baseline piecewise parabolic advection scheme, which captures the formation of smaller-scale structures superimposed on larger structures of the order of the shear-layer width.

Edwards, Jack R.↗

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↗

Propulsive re-entry aerodynamics Interim report, 1 Jun. 1967 - 1 Jun. 1968

This report discusses the progress accomplished during the first year of a two-year contract to study Propulsive Re-Entry Aerodynamics. Analytical methods were developed during the first year to predict the flow field in the immediate vicinity of a planetary entry system composed of an aeroshell with a single retrorocket which exhausts into a subsonic or supersonic counterflowing planetary atmosphere. In the second year a digital computer program will be developed for the analytical solutions, and the analytical model will be tested by performing an experimental test program and comparing the experimental results with the analytical predictions. The analytical investigation considers the flow field due to a single supersonic retrorocket exhausting from an aeroshell into oncoming subsonic, transonic, or supersonic streams in the large thrusting coefficient regime, 0.5 < C(sub T) = T/q(sub ∞) A(sub m) < 15, which is the range of thrusting coefficients encountered during the terminal landing phase for planetary landing on Mars or Venus. Analyses are provided for all the geometrical features of the flow field such as the location of the terminal shock, jet boundary, and the profile for the interface between the jet and atmospheric gases. In addition, analyses are developed which describe the growth of the mixing layers along the interface and at the edge of the jet boundary and dead-air region and the reattachment of the shear layer to the aeroshell or the recompression of the shear layer in wake type flow. A comprehensive discussion of the experimental wind tunnel test program to be accomplished in Phase II is presented.

Philip O Jarvinen↗

An integral turbulent kinetic energy analysis of free shear flows

Mixing of coaxial streams is analyzed by application of integral techniques. An integrated turbulent kinetic energy (TKE) equation is solved simultaneously with the integral equations for the mean flow. Normalized TKE profile shapes are obtained from incompressible jet and shear layer experiments and are assumed to be applicable to all free turbulent flows. The shear stress at the midpoint of the mixing zone is assumed to be directly proportional to the local TKE, and dissipation is treated with a generalization of the model developed for isotropic turbulence. Although the analysis was developed for ducted flows, constant-pressure flows were approximated with the duct much larger than the jet. The axisymmetric flows under consideration were predicted with reasonable accuracy. Fairly good results were also obtained for the fully developed two-dimensional shear layers, which were computed as thin layers at the boundary of a large circular jet.

Peters, C. E.↗

Plume effects on the flow around a blunted cone at hypersonic speeds

Tests at M = 8.2 show that a simulated rocket plume at the base of a blunted cone can cause large areas of separated flow, with dramatic effects on the heat transfer rate distribution. The plume was simulated by solid discs of varying sizes or by an annular jet of gas. Flow over the cone without a plume is fully laminar and attached. Using a large disc, the boundary layer is laminar at separation at the test Reynolds number. Transition occurs along the separated shear layer and the boundary layer quickly becomes turbulent. The reduction in heat transfer associated with a laminar separated region is followed by rising values as transition occurs and the heat transfer rates towards the rear of the cone substantially exceed the values obtained without a plume. With the annular jet or a small disc, separation occurs much further aft, so that heat transfer rates at the front of the cone are comparable with those found without a plume. Downstream of separation the shear layer now remains laminar and the heat transfer rates to the surface are significantly lower than the attached flow values.

Atcliffe, P.↗

Far-field noise of a subsonic jet under controlled excitation

The role of large-scale coherent structures in broadband jet noise suppression and amplification under controlled excitation are studied experimentally. It is found that suppression occurs only at low Reynolds and Mach numbers, when the boundary layer at the jet exit is laminar; the optimum St(Theta) (St, based on the initial shear-layer momentum thickness) value is 0.017, at which the excitation results in a quick roll-up and transition of the laminar shear-layer vortices, yielding coherent structures similar to those existing at high speeds. At the asymptotic level the broadband jet noise can only be amplified by the excitation. The amplification is inferred to be maximum for excitation in the St(D) (St based on the jet diameter) range of 0.65-0.85. Finally, it is noted that the pairing process induced by the excitation is at the origin of the broadband noise amplification.

Zaman, K. B. M. Q.↗

Filtered Rayleigh scattering based measurements in compressible mixing layers

Results are presented of experiments in which a Rayleigh-scattering-based technique was used to make quantitative planar measurements in the compressible free shear layers. To this end, the absorption characteristics of the iodine molecular filter were investigated, and initial planar density measurements in M(c) = 0.51 shear layer were performed. The preliminary results presented here demonstrate the validity of the measurement technique.

Elliott, G. S.↗

Cold-Flow Study of Low Frequency Pressure Instability in Hybrid Rocket Motors

Past experience with hybrid rockets has shown that certain motor operating conditions are conducive to the formation of low frequency pressure oscillations, or flow instabilities, within the motor. Both past and present work in the hybrid propulsion community acknowledges deficiencies in the understanding of such behavior, though it seems probable that the answer lies in an interaction between the flow dynamics and the combustion heat release. Knowledge of the fundamental flow dynamics is essential to the basic understanding of the overall stability problem. A first step in this direction was a study conducted at NASA Marshall Space Flight Center (MSFC), centered around a laboratory-scale two dimensional water flow model of a hybrid rocket motor. Principal objectives included: (1) visualization of flow and measurement of flow velocity distributions: (2) assessment of the importance of shear layer instabilities in driving motor pressure oscillations; (3) determination of the interactions between flow induced shear layers with the mainstream flow, the secondary (wall) throughflow, and solid boundaries; (4) investigation of the interactions between wall flow oscillations and the mainstream flow pressure distribution.

Jenkins, Rhonald M.↗

Performance of Low Dissipative High Order Shock-Capturing Schemes for Shock-Turbulence Interactions

Accurate and efficient direct numerical simulation of turbulence in the presence of shock waves represents a significant challenge for numerical methods. The objective of this paper is to evaluate the performance of high order compact and non-compact central spatial differencing employing total variation diminishing (TVD) shock-capturing dissipations as characteristic based filters for two model problems combining shock wave and shear layer phenomena. A vortex pairing model evaluates the ability of the schemes to cope with shear layer instability and eddy shock waves, while a shock wave impingement on a spatially-evolving mixing layer model studies the accuracy of computation of vortices passing through a sequence of shock and expansion waves. A drastic increase in accuracy is observed if a suitable artificial compression formulation is applied to the TVD dissipations. With this modification to the filter step the fourth-order non-compact scheme shows improved results in comparison to second-order methods, while retaining the good shock resolution of the basic TVD scheme. For this characteristic based filter approach, however, the benefits of compact schemes or schemes with higher than fourth order are not sufficient to justify the higher complexity near the boundary and/or the additional computational cost.

Sandham, N. D.↗

Performance of Low Dissipative High Order Shock-Capturing Schemes for Shock-Turbulence Interactions

Abstract Accurate and efficient direct numerical simulation of turbulence in the presence of shock waves represents a significant challenge for numerical methods. The objective of this paper is to evaluate the performance of high order compact and non-compact central spatial differencing employing total variation diminishing (TVD) shock-capturing dissipations as characteristic based filters for two model problems combining shock wave and shear layer phenomena. A vortex pairing model evaluates the ability of the schemes to cope with shear layer instability and eddy shock waves, while a shock wave impingement on a spatially-evolving mixing layer model studies the accuracy of computation of vortices passing through a sequence of shock and expansion waves. A drastic increase in accuracy is observed if a suitable artificial compression formulation is applied to the TVD dissipations. With this modification to the filter step the fourth-order non-compact scheme shows improved results in comparison to second-order methods, while retaining the good shock resolution of the basic TVD scheme. For this characteristic based filter approach, however, the benefits of compact schemes or schemes with higher than fourth order are not sufficient to justify the higher complexity near the boundary and/or the additional computational cost.

Sandham, N. D.↗

Transitional Bubble in Periodic Flow Phase Shift

One particular characteristic observed in unsteady shear layers is the phase shift relative to the main flow. In attached boundary layers this will have an effect both on the instantaneous skin friction and heat transfer. In separation bubbles the contribution to the drag is dominated by the pressure distribution. However, the most significant effect appears to be the phase shift on the transition process. Unsteady transition behaviour may determine the bursting of the bubble resulting in an un-recoverable full separation. An early analysis of the phase shift was performed by Stokes for the incompressible boundary layer of an oscillating wall and an oscillating main flow. An amplitude overshoot within the shear layer as well as a phase shift were observed that can be attributed to the relatively slow diffusion of viscous stresses compared to the fast change of pressure. Experiments in a low speed facility with the boundary layer of a flat plate were evaluated in respect to phase shift. A pressure distribution similar to that on the suction surface of a turbomachinery aerofoil was superimposed generating a typical transitional separation bubble. A periodically unsteady main flow in the suction type wind tunnel was introduced via a rotating flap downstream of the test section. The experiments covered a range of the three similarity parameters of momentum-loss-thickness Reynolds-number of 92 to 226 and Strouhal-number (reduced frequency) of 0.0001 to 0.0004 at the separation point, and an amplitude range up to 19 %. The free stream turbulence level was less than 1% .Upstream of the separation point the phase shift in the laminar boundary layer does not appear to be affected significantly bay either of the three parameters. The trend perpendicular to the wall is similar to the Stokes analysis. The problem scales well with the wave velocity introduced by Stokes, however, the lag of the main flow near the wall is less than indicated analytically. The separation point comes closest to the Stokes analysis but the phase is still 20 degrees lower at the wall.

Talan, M.↗

Zonal Two Equation Kappa-Omega Turbulence Models for Aerodynamic Flows

Two new versions of the kappa-omega two-equation turbulence model will be presented. The new Baseline (BSL) model is designed to give results similar to those of the original kappa-omega model of Wilcox, but without its strong dependency on arbitrary freestream values. The BSL model is identical to the Wilcox model in the inner 50% of the boundary-layer but changes gradually to the standard kappa-epsilon model (in a kappa- omega formulation) towards the boundary-layer edge. The free shear layers. The second version of the model is called Shear-Stress Transport (SST) model. It is a variation of the BSL model with the additional ability to account for the transport of the principal turbulent shear stress in adverse pressure gradient boundary-layers. The model is based on Bradshaw's assumption that the principal shear-stress is proportional to the turbulent kinetic energy, which is introduced into the definition of the eddy-viscosity. Both models are tested for a large number of different flowfields. The results of the BSL model are similar to those of the original kappa-omega model, but without the undesirable freestream dependency. The predictions of the SST model are also independent of the freestream values but show better agreement with experimental data for adverse pressure gradient boundary-layer flows.

Menter, Florian R.↗

Investigation of supersonic modes and three-dimensionality in bounded, free shear flows

It is found through two-dimensional temporal simulations of high-speed free shear layers that mean flow distortion is significantly increased when supersonic disturbances are introduced as initial conditions. The shear layer exhibits no subharmonic growth or roll-up, but rather a spectral broadening as energy is distributed into higher harmonics. Increasing the velocity of one side of the mixing layer (u2) to roughly 1/5 the speed of the high speed side (u1), allows a slight subharmonic growth at a very slow rate for two-dimensional modes. A first look at three-dimensional free shear flows is also presented for M = 2. No effect is seen for incompressible flow; however, stabilization is seen with respect to maximum temporal growth rates as the transverse velocity is increased. A much stronger, but similar effect is seen if u2 is increased. The wave direction of maximum growth for u2 is greater than 0.4 is found to be the direction of the faster stream (u1) over a broad range of transverse velocities.

Macaraeg, Michele G.↗

Physical description of boundary-layer transition: Experimental evidence

The problems of understanding the origins of turbulent flow and transition to turbulent flow are the most important unsolved problems of fluid mechanics and aerodynamics. It is well known that the stability, transition, and turbulent characteristics of bounded shear layers are fundamentally different from those of free shear layers. Likewise, the stability, transition, and turbulent characteristics of open systems are fundamentally different from those of closed systems. Because of the influence of indigenous disturbances, surface geometry and roughness, sound, heat transfer, and ablation, it is not possible to develop general prediction schemes for transition location and the nature of turbulent structures in boundary-layer flows. At the present time no mathematical model exists that can predict the transition Reynolds number on a flat plate. The recent progress in this area is encouraging, in that a number of distinct transition mechanisms have been found experimentally. The theoretical work finds them to be amplitude and Reynolds-number dependent. The theory remains rather incomplete with regard to predicting transition. Amplitude and spectral characteristics of the disturbances inside the laminar viscous layer strongly influence which type of transition occurs. The major need in this area is to understand how freestream disturbances are entrained into the boundary layer, i.e., to answer the question of receptivity. We refer receptivity to the mechanism(s) that cause freestream disturbances to enter the boundary layer and create the initial amplitudes for unstable waves.

Saric, William S.↗

Experimental Investigation of Separated and Transitional Boundary Layers Under Low-Pressure Turbine Airfoil Conditions

Modern low-pressure turbine airfoils are subject to increasingly stronger pressure gradients as designers impose higher loading in an effort to improve efficiency and to reduce part count. The adverse pressure gradients on the suction side of these airfoils can lead to boundary-layer separation, particularly under cruise conditions. Separation bubbles, notably those which fail to reattach, can result in a significant degradation of engine efficiency. Accurate prediction of separation and reattachment is hence crucial to improved turbine design. This requires an improved understanding of the transition flow physics. Transition may begin before or after separation, depending on the Reynolds number and other flow conditions, has a strong influence on subsequent reattachment, and may even eliminate separation. Further complicating the problem are the high free-stream turbulence levels in a real engine environment, the strong pressure gradients along the airfoils, the curvature of the airfoils, and the unsteadiness associated with wake passing from upstream stages. Because of the complicated flow situation, transition in these devices can take many paths that can coexist, vary in importance, and possibly also interact, at different locations and instances in time. The present work was carried out in an attempt to systematically sort out some of these issues. Detailed velocity measurements were made along a flat plate subject to the same nominal dimensionless pressure gradient as the suction side of a modern low-pressure turbine airfoil ('Pak-B'). The Reynolds number based on wetted plate length and nominal exit velocity, Re, was varied from 50;000 to 300; 000, covering cruise to takeoff conditions. Low, 0.2%, and high, 7%, inlet free-stream turbulence intensities were set using passive grids. These turbulence levels correspond to about 0.2% and 2.5% turbulence intensity in the test section when normalized with the exit velocity. The Reynolds number and free-stream turbulence level do not have a significant effect on the location of boundary-layer separation unless they are high enough to induce transition upstream of separation. The location and extent of the transition zone, in contrast, depend strongly on Re and TI. The beginning of reattachment closely follows the onset of transition. Under low free-stream turbulence conditions the boundary layer is laminar at separation and then begins to exhibit fluctuations in a finite frequency band in the shear layer over the separation bubble. These fluctuations are due to instability waves. The fluctuations grow in magnitude, higher harmonics are generated, and finally lead to a breakdown to turbulence. Transition begins in the shear layer, but quickly spreads to the near wall region and causes the boundary layer to reattach. The transition is rapid and the resulting turbulence contains a full range of high and low frequencies. Under high free-stream turbulence conditions, slowly growing low-frequency fluctuations are induced in the pretransitional boundary layer by the free-stream. The separation bubbles are considerably thinner than in the low TI cases, resulting in thinner boundary layers at the end of the test wall. At Re=50,000 and 100,000, the pre-transitional boundary layer separates at about the same location as in the low TI cases. Transition occurs through a bypass mode, begins upstream of the corresponding low-TI location, and proceeds in a manner similar to that of an attached boundary layer. Under high TI at Re=200,000 and 300,000, transition begins before separation. The boundary layer may separate, but if it does the separation bubble is very short and does not significantly affect the downstream development of the boundary layer. A comparison is made to previous work in a simulated cascade.

Hultgren, Lennart S.↗