Supersonic flow separation in planar nozzles

We present experimental results on separation of supersonic flow inside a convergent–divergent (CD) nozzle. The study is motivated by the occurrence of mixing enhancement outside CD nozzles operated at low pressure ratio. A novel apparatus allows investigation of many nozzle geometries with large optical access and measurement of wall and centerline pressures. The nozzle area ratio ranged from 1.0 to 1.6 and the pressure ratio ranged from 1.2 to 1.8. At the low end of these ranges, the shock is nearly straight. As the area ratio and pressure ratio increase, the shock acquires two lambda feet. Towards the high end of the ranges, one lambda foot is consistently larger than the other and flow separation occurs asymmetrically. Downstream of the shock, flow accelerates to supersonic speed and then recompresses. The shock is unsteady, however, there is no evidence of resonant tones. The separation shear layer on the side of the large lambda foot exhibits intense instability that grows into large eddies near the nozzle exit. Time-resolved wall pressure measurements indicate that the shock oscillates in a piston-like manner and most of the energy of the oscillations is at low frequency.      

Experimental flow study over a blunt-nosed axisymmetric body at incidence

The flowfield over a blunt-nosed cylinder was examined experimentally at a low subsonic speed for Re=1.88×10⁵ and angles of attack up to 40°. Velocity measurements were carried out (employing a seven-hole Pitot tube) as well as wall static pressure and wall shear-stress measurements. Surface flow visualization was applied using liquid crystals and a mixture of oil–TiO₂. For all the examined cases no flow asymmetries were found. For high angles of attack (20° and above) a separation “bubble” appears at the leeside of the nose area (streamwise flow separation). The basic feature of the circumferential pressure distribution at the after body area for these angles of attack is a plateau close to the suction peak and a fast recovery next to it. One streamwise vortex on each side of the symmetry plane is formed as well as a separation bubble about 90° far from this plane, where the cross-flow primary separation line is located. Each cross-flow primary separation line starts at the leeside nose area and moves towards the windward side along the cylindrical after body. The space between the two primary separation lines close to the wall is characterized by high flow fluctuations on the leeside, compared to the low fluctuations of the windward side.     

Experimental investigation of the wall dynamics of the A-pillar vortex flow

The influence of the A-pillar vortex on the wall flow of the side window of a car is investigated experimentally using a 30° dihedron model. The measurement of the unsteady pressure at the wall provides a map of the pressure fluctuation intensity, and a spectral analysis is performed to track the dominant frequencies of the wall pressure fluctuations. The wall flow generated by the vortex structure that develops parallel to the side-wall is characterised by means of particle image velocimetry (PIV). Its structure is analysed and compared to cross-sections of the A-pillar vortex in order to identify the different separation and reattachment lines. A comparison of the field of turbulent kinetic energy obtained by PIV with the map of the pressure fluctuations shows a correlation between the structure of the A-pillar vortex and the pressure fluctuations. It is found that the dominant wall pressure fluctuations are located at the secondary separation line, whereas the primary reattachment line does not show any significant pressure variations, that the A-pillar vortex will not naturally break down and that discrete vortices may be associated with the pressure fluctuations.     

Transient three-dimensional startup side load analysis of a regeneratively cooled nozzle

The objective of this effort is to develop a computational methodology to capture the side load physics and to anchor the computed aerodynamic side loads with the available data by simulating the startup transient of a regeneratively cooled, high-aspect-ratio nozzle, hot-fired at sea level. The computational methodology is based on an unstructured-grid, pressure-based, reacting flow computational fluid dynamics and heat transfer formulation, and a transient inlet history based on an engine system simulation. Emphases were put on the effects of regenerative cooling on shock formation inside the nozzle, and ramp rate on side load reduction. The results show that three types of asymmetric shock physics incur strong side loads: the generation of combustion wave, shock transitions, and shock pulsations across the nozzle lip, albeit the combustion wave can be avoided with sparklers during hot-firing. Results from both regenerative cooled and adiabatic wall boundary conditions capture the early shock transitions with corresponding side loads matching the measured secondary side load. It is theorized that the first transition from free-shock separation to restricted-shock separation is caused by the Coanda effect. After which the regeneratively cooled wall enhances the Coanda effect such that the supersonic jet stays attached, while the hot adiabatic wall fights off the Coanda effect, and the supersonic jet becomes detached most of the time. As a result, the computed peak side load and dominant frequency due to shock pulsation across the nozzle lip associated with the regeneratively cooled wall boundary condition match those of the test, while those associated with the adiabatic wall boundary condition are much too low. Moreover, shorter ramp time results show that higher ramp rate has the potential in reducing the nozzle side loads.

---

     

A method for separating shock wave motion and turbulence in LDV measurements

 Two-component laser Doppler velocimetry (LDV) measurements were made in a planar, two-dimensional flow containing an unsteady oblique shock wave formed by the convergence of two supersonic streams past a thick plate. High-speed wall pressure measurements locate the shock wave and, consequently, allow separation of the effects of shock wave motion from the turbulence fluctuations in the LDV measurements of the shock-separated free shear layer. In the current flow isolating the large-scale changes in the position of the shock from the turbulence primarily reduces the experimental scatter rather than significantly changing the shapes or magnitudes of the turbulent stress profiles. Changes in the direction of shock motion do not significantly affect the mean velocity, but do affect the turbulent stresses. Received: 11 August 1997/Accepted: 30 September 1998     

Experimental study of boundary layer separation in truncated ideal contour nozzles

DLR Lampoldshausen carried out a cold flow test series to study the boundary layer separation and the related flow field in a truncated ideal contour nozzle. A special focus was set on low nozzle pressure ratios to identify the origin of a locally re-attached flow condition that was detected in previous test campaigns. A convex shaped Mach disc was found for nozzle pressure ratios less than 10 and a slight concave one for nozzle pressure ratios more than 20. Due to boundary layer transition at low nozzle pressure ratios the convex Mach disc is temporary tilted and redirects the flow towards the nozzle wall. A simple separation criterion for turbulent nozzle flows is presented that fits well for both cold and hot flows. It is shown that the oblique separation shock recompresses the flow to 90% of the ambience. The separation zone of the presented film cooled nozzle is compared with a conventional one around 40% longer. Furthermore a relation between shear layer shape and forced side loads is described.      

Origin of flow asymmetry in planar nozzles with separation

An experimental investigation was conducted to study the mechanisms that lead to the origin of flow asymmetry in overexpanded planar nozzles, especially at low nozzle pressure ratios. Three Mach 2 planar nozzles with different divergent wall angles but same area-ratio were tested. For all three nozzles, a large portion of the dimensional pressure rise data across the separation shock shows the nature of boundary layer to be in the laminar/transitional state. Depending upon the local flow conditions, the flow can, therefore, experience either an early or a delayed separation on either wall. This can result in a free or a restricted shock separation condition on either wall which can initiate the beginning of flow asymmetry in nozzles at low nozzle pressure ratio. However, a higher nozzle wall angle was observed to prevent initiation of such a flow asymmetry. The present tests, therefore, indicate that in addition to the state of the boundary layer along the nozzle wall, the proximity of the separated shear layer to the nozzle walls also seems to play a dominant role in initiating conditions that favor the origin of flow asymmetry in nozzles. A significant drop in the shock unsteadiness levels was also indicated by increasing the wall angle.     

Shock unsteadiness in a thrust optimized parabolic nozzle

This paper discusses the nature of shock unsteadiness, in an overexpanded thrust optimized parabolic nozzle, prevalent in various flow separation modes experienced during start up and shut down sequences. The results are based on simultaneously acquired data from real-time wall pressure measurements using Kulite pressure transducers, high-speed schlieren (2 kHz) of the exhaust flow-field and from strain-gauges installed on the nozzle bending tube. Shock unsteadiness in the separation region is seen to increase significantly just before the onset of each flow transition, even during steady nozzle operation. The intensity of this measure (rms level) is seen to be strongly influenced by relative locations of normal and overexpansion shock, the decrease in radial size of re-circulation zone in the back-flow region, and finally, the local nozzle wall contour. During restricted shock separation, the pressure fluctuations in separation region exhibit periodic characteristics rather than the usually observed characteristics of intermittent separation. The possible physical mechanisms responsible for the generation of flow unsteadiness in various separation modes are discussed. The results are from an experimental study conducted in P6.2 cold-gas subscale test facility using a thrust optimized parabolic nozzle of area-ratio 30.

---

     

The Effect of the Mach Number on a Turbulent Backward-Facing Step Flow

The flow around a backward-facing step in the sub-, trans- and supersonic regimes was investigated at the Trisonic Wind Tunnel Munich with particle image velocimetry and dynamic pressure measurements. These two techniques were combined to simultaneously measure and correlate the velocity fluctuations in a streamwise vertical plane with the pressure fluctuations on the reattachment surface. The results show that the dynamic loads on the reattachment surface increase from subsonic up to the transonic regime while the mean reattachment location moves downstream. As soon as the flow becomes locally supersonic aft of the backward-facing step, the mean reattachment location suddenly moves upstream while the normalized dynamic loads drastically decrease. By correlating the velocity and the dynamic pressure data, it was shown that a clear separation between outer flow and the flow close to the surface aft of the step is responsible for the drastic load reduction. Due to the large difference in pressure/density, the disturbances from the locally supersonic flow do not have an effect on the flow close to the surface. This is also reflected in the power spectral densities of the pressure fluctuations on the surface, showing that at supersonic free-stream Mach numbers a low-frequency pumping motion of the locally subsonic flow is the dominant mode, while in sub-/transonic flow Kelvin-Helmholtz instabilities and a cross-pumping motion of the shear layer dominate the dynamic loads.     

Wall pressure fluctuations and topology in separated flows over a forward-facing step

Laboratory measurements of wall pressure fluctuations and aerodynamic fields were made in separated flows over a forward facing step (h = 30, 40 and 50 mm with U _e = 15–40 m/s). An array of 16 off-set pressure probes extending in the streamwise and the spanwise directions was especially developed for sensing the wall pressure fluctuations. The flow field was also investigated by wall flow visualizations and PIV to analyze the flow topology in an open section wind tunnel. The results show a different behavior of the flow depending on the aspect ratio l/h and δ/h for high Reynolds numbers. The space time correlations between the wall pressure and the velocity fields were highlighted. The results show that high levels of these correlations are located at the top of the recirculation bubble, mainly in the shear layer and are extended downstream of the re-attachment point. Indeed, the results indicate that the flapping motion at the separation is important in the flow organization at the re-attachment point.     

Wall temperature effects on shock unsteadiness in a reattaching compressible turbulent shear layer

In this study, the effect of heat transfer on the compressible turbulent shear layer and shockwave interaction in a scramjet has been investigated. To this end, highly resolved Large Eddy Simulations (LES) are performed to explore the effect of wall thermal conditions on the behavior of a reattaching free shear layer interacting with an oblique shock in compressible turbulent flows. Various wall-to-recovery temperature ratios are considered, and results are compared to the adiabatic wall. It is found that the wall temperature affects the reattachment location and the shock behavior in the interaction region. Furthermore, fluctuating heat flux exhibits a strong intermittent behavior with severe heat transfer compared to the mean, characterized by scattered spots. The distribution of the Stanton number shows a strong heat transfer and complex pattern within the interaction, with the maximum thermal (heat transfer rates) and dynamic loads (root-mean-square wall pressure) found for the case of the cold wall. The analysis of LES data reveals that the thermal boundary condition can significantly impact the wall pressure fluctuations level. The primary mechanism for changes in the flow unsteadiness due to the wall thermal condition is linked to the reattaching shear layer, which agrees with the compressible turbulent boundary layer theory.     

Wall pressure fluctuations in the reattachment region of a supersonic free shear layer

A study was made of the wall pressure fluctuations in the reattachment region of a supersonic free shear layer. The free shear layer was formed by the separation of a Mach 2.9 turbulent boundary layer from a backward facing step. Reattachment occurred on a 20° ramp. By adjusting the position of the ramp, the base pressure at the step was set equal to the freestream pressure, and the free shear layer formed in the absence of any turning. An array of flush-mounted, miniature, high-frequency pressure transducers was used in the vicinity of the reattachment region to make multichannel measurements of the fluctuating wall pressure. Contrary to previous observations of this flow, the reattachment region was found to be highly unsteady, and the pressure fluctuations were found to be large. The overall behavior of the wall pressure loading is similar in scale and magnitude to the unsteadiness of the wall pressure field in compression ramp flows at the same Mach number. Rayleigh scattering was used to visualize the instantaneous shock structure in the streamwise and spanwise direction. Spanwise wrinkles on the order of half the boundary layer thickness were observed on the shock sheet.     

Dual-element directional shock wave attenuators

This study explores the interaction between shock waves and dual-element porous plates used to ameliorate the hazardous effects of these waves. Tests were performed in a shock tube to determine the effects that a pair of porous plates with directional resistance properties had on the initial peak pressure and impulse amelioration experienced by an end wall. Mild steel test specimens, ranging in porosity values from 6.6% to 41.1%, were mounted two at a time at different spacings in the test section. Each plate had directional properties, i.e. resistance to flow was different for flow coming from either side. Four plates were used, and 48 plate configurations were tested. Side wall and end wall pressure measurements and schlieren photographs were taken of the interactions. Tests were run at Mach numbers of 1.23, 1.35 and 1.42. The separation distances between the plate specimens were varied between 30 mm and 60 mm; however, the distance between the downstream plate and the end wall was kept constant at 140 mm for all tests.Both the initial peak pressure and impulse amelioration values were found to be dependant on the plate combination porosity. As the porosity of the combination increased, the amelioration values decreased. Complementary plate combinations produced differing results as different wave interactions occur when plate positions were interchanged. The porosity of the combined plates was found to have an overriding influence on the end wall peak pressure and impulse amelioration values when compared to the effect that plate arrangement (i.e. geometrical influences) had. Impulse amelioration values were found to increase as the separation distance between plates were increased. The amplitude of the end wall pressure trace was found to increase as the incident Mach number was increased. The significant attenuation of the incident shock wave obtained during this study is attributed to the system of multiple reflected and transmitted waves that are produced by the presence of the plate specimens in series. This increases the frequency of shock wave and barrier interactions, when compared to just using a single barrier.     

翼型跨声速气动特性的不确定性及全局灵敏度分析

针对马赫数和仰角的随机不确定性会导致气动性能波动的现象, 采用非嵌入式的混沌多项式方法对绕NACA0012 翼型跨声速随机气动特性进行不确定性及全局灵敏度分析. 具体分析了飞行状态的不确定性对气动载荷分布、流场及气动力系数的影响并通过全局灵敏度分析找出重要因素. 不确定性分析结果表明翼型上表面的激波以及激波后分离泡是造成气动性能剧烈波动的主要原因. 灵敏度分析结果表明在跨声速区域马赫数对激波处气动性能影响最大, 此外, 虽然马赫数和仰角相互耦合作用对气动力系数贡献比较小, 但对于激波位置处的流场, 这种互耦合作用不可忽略.      

Cold-gas experiments to study the flow separation characteristics of a dual-bell nozzle during its transition modes

An experimental investigation has been carried out to study the effect of test environment on transition characteristics and the flow unsteadiness associated with the transition modes of a dual-bell nozzle. Cold-gas tests using gaseous nitrogen were carried out in (i) a horizontal test-rig with nozzle exhausting into atmospheric conditions and, (ii) a high altitude simulation chamber with nozzle operation under self-evacuation mode. Transient tests indicate that increasing δP ₀/δt (the rate of stagnation chamber pressure change) reduces the amplitude of pressure fluctuations of the separation shock at the wall inflection point. This is preferable from the viewpoint of lowering the possible risk of any structural failure during the transition mode. Sea-level tests show 15–17% decrease in the transition nozzle pressure ratio (NPR) during subsequent tests in a single run primarily due to frost formation in the nozzle extension up to the wall inflection location. Frost reduces the wall inflection angle and hence, the transition NPR. However, tests inside the altitude chamber show nearly constant NPR value during subsequent runs primarily due to decrease in back temperature with decrease in back pressure that prevents any frost formation.     

Computational analysis of impinging shock-wave boundary layer interaction under conditions of incipient separation

The interaction of an oblique shock wave with a turbulent boundary layer under conditions of incipient separation is analyzed by means of large-eddy simulation (LES) and Reynolds-averaged Navier–Stokes (RANS) turbulence models, with the objective to explore their predictive capabilities, in particular with respect to the unsteady features of the interaction. Consistent with earlier direct numerical simulations, we have found that the flow dynamics in the interaction zone is characterized by strong intermittency associated with the formation of scattered spots of flow reversal near the nominal position of the reflected shock. Comparison with experimental results (at much larger Reynolds number) show that the qualitative features of the interaction are predicted reasonably well by both LES and RANS models. RANS models supplemented with a semi-empirical closure are also found to provide reasonable estimate of the fluctuating pressure loads at the wall.     

Numerical calculation of FSS/RSS transition in highly overexpanded rocket nozzle flows

Turbulent flow separation in over-expanded rocket nozzles is investigated numerically in a sub-scale parabolic nozzle fed with cold nitrogen. Depending upon the feeding to ambient pressure ratio either a free shock separation or a restricted shock separation is computed, with a significant hysteresis between these two flow regimes. This hysteresis was also found in experimental tests with the same nozzle geometry. The present study is mainly focused on the transition between the two shock separation patterns. The analysis of the numerical solutions aims to provide clues for the explanation of the hysteresis cycle.     

Experimental evaluation of side-loads in LE-7A prototype engine nozzle

During development tests of the LE-7A prototype engine, severe side-loads were observed. The side-load peaks appeared only in certain limited conditions during start-up and shut-down transients. To investigate phenomena causing those severe side-loads observed in the LE-7A prototype engine nozzle, series of cold-flow tests and hot-firing tests as well as CFD analyses were conducted. As a result of the hot-firing tests, two different phenomena were found to cause severe side-loads in the LE-7A prototype engine nozzle. One was a restricted shock separation (RSS) flow structure and the other was a phenomenon termed “separation jump,” the rapid movement of the separation location in the vicinity of the step. A step was installed in the LE-7A prototype to supply film-cooling gas. Hot-firing test results showed that RSS can occur for a limited mixture ratio. Detailed flow structure of RSS on the nozzle surface was revealed by the cold-flow tests. Measured pressures and visualized images of cold-flow tests clarified the mechanism causing the separation jump. The key phenomenon ruling the separation jump was found to be the base flow behind the step. Based on the results of the present study, the latest LE-7A engine nozzle design has been changed to eliminate the severe side-load.

---

     

Delayed detached eddy simulation of the end-effect regime and side-loads in an overexpanded nozzle flow

The separated flow in an overexpanded nozzle featuring a restricted shock separation is investigated numerically using delayed detached eddy simulation and compared with the experimental data of Nguyen et al. (Int J Flow Turbul Combust 71(1):161–181, 2003). First, the enormous cost of a Large Eddy Simulation for such a nozzle flow is assessed before being performed to motivate the practical need for using an hybrid RANS/LES method. The calculation is then used to investigate the “end-effect” regime which involves a strong global unsteadiness with very large amplitude fluctuations of about 15–20% of nozzle divergent length. The flow regime is characterized by high wall pressure fluctuations which are hopefully nearly axisymmetric. The main properties (rms levels, amplitude of displacement of the separation) of the motion are rather well reproduced by DDES compared to the experiment. However, a major difference lies in the frequency of the computed motion which is higher than in the experiment. This major discrepancy is currently not explained by the author. The properties of the side-loads are also briefly discussed.      

Assessment of the organization of a turbulent separated and reattaching flow by measuring wall pressure fluctuations

The effect of local forcing on the organization of a turbulent separated and reattaching flow was assessed by measuring wall pressure fluctuations. Multi-arrayed microphones were installed on the surface to measure the simultaneous spatial and temporal wall pressure fluctuations. Local forcing at the separation edge was applied to the separated flow over a backward-facing step through a thin slit. The organization of the separated and reattaching flow was found to be greatest at the effective forcing frequency. The flow structure was diagnosed by analyzing several characteristics of the wall pressure fluctuations: the wall pressure fluctuation coefficients, wall pressure spectrum, wavenumber-frequency spectrum, coherence, cross-correlation, and multi-resolution autocorrelations of pressure fluctuations using the maximum overlap discrete wavelet transform and continuous wavelet transform. Features indicative of the amalgamation of vortices under the local forcing were observed; this amalgamation process accounted for the observed reduction of the reattachment length. Examination of the wall pressure fluctuations revealed that introduction of local forcing enhanced flapping motion as well as the streamwise and spanwise dispersions of vortical structures.