Experimental study of vortex structures in a hypersonic shock layer on a flat plate

The characteristics of travelling perturbations of density in a hypersonic shock layer on a flat plate for the Mach number M_∞=21 and unit Reynolds numberRe ₁=6·10⁵ m⁻¹ were experimentally studied by the method of electron-beam fluorescence. The perturbations were generated by interaction of the shock layer behind an oblique gas-dynamic whistle and the leading edge of the plate. The cases of unsteady and quasi-steady interaction were considered. In both cases, vortex disturbances of finite amplitude were generated. The measurements were performed at the fundamental frequency F=0.6·10⁻⁴ and at the harmonic; the streamwise phase velocities, the growth rates of the disturbances, and the angles of wave propagation were obtained. The measurement results are compared with some experimental data for subsonic flows, some particular results of the linear stability theory for compressible flows, and the results obtained on the basis of a simple model of the nonlinear stage of disturbance evolution in a hypersonic boundary layer. Institute of Theoretical and Applied Mechanics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 6, pp. 41–47, November–December, 1999.     

Wave processes in the shock layer on a flat plate at an angle of attack

A numerical and experimental study of receptivity of the viscous shock layer on a flat plate aligned at an angle of attack to external acoustic perturbations is performed. Density and pressure fluctuations are measured in experiments at the free-stream Mach number M _∞ = 21 and Reynolds number Re ₁ = 6·10 ⁵ m ⁻¹ . Direct numerical simulations of receptivity of the viscous shock layer to external acoustic perturbations in wide ranges of the governing parameters are performed by solving the Navier-Stokes equations with the use of high-order shock-capturing schemes. The calculated intensities of density and pressure fluctuations are found to be in good agreement with experimental data. Results of the study show that entropy-vortex disturbances dominate in the shock layer at small angles of attack, whereas acoustic perturbations prevail at angles of attack above 20°.     

Experimental study of evolution of disturbances in a supersonic boundary layer on a swept-wing model under controlled conditions

Experimental data on stability of a three-dimensional supersonic boundary layer on a swept wing are presented. Evolution of artificial wave trains was studied. The experiments were conducted for Mach numberM=2.0 and unit Reynolds numberRe ₁=6.6·10⁶m⁻¹ on a swept-wing model with a lenticular profile and a40° sweep angle of the leading edge at zero incidence. Excitation of high-frequency disturbances caused by secondary-flow instability at a high initial amplitude was observed. It is shown that the evolution of disturbances at frequencies of10, 20, and30 kHz is similar to the development of travelling waves for the case of subsonic velocities. Institute of Theoretical and Applied Mechanics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 1, pp. 50–56, January–February, 2000.     

Experimental study and direct numerical simulation of the evolution of disturbances in a viscous shock layer on a flat plate

The evolution of disturbances in a hypersonic viscous shock layer on a flat plate excited by slow-mode acoustic waves is considered numerically and experimentally. The parameters measured in the experiments performed with a free-stream Mach number M _∞ = 21 and Reynolds number Re _L = 1.44 · 10⁵ are the transverse profiles of the mean density and Mach number, the spectra of density fluctuations, and growth rates of natural disturbances. Direct numerical simulation of propagation of disturbances is performed by solving the Navier-Stokes equations with a high-order shock-capturing scheme. The numerical and experimental data characterizing the mean flow field, intensity of density fluctuations, and their growth rates are found to be in good agreement. Possible mechanisms of disturbance generation and evolution in the shock layer at hypersonic velocities are discussed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 5, pp. 3–15, September–October, 2006.     

The ultra-low Reynolds number airfoil wake

Lift force and the near wake of an NACA 0012 airfoil were measured over the angle (α) of attack of 0°–90° and the chord Reynolds number (Re _c ), 5.3 × 10³–5.1 × 10⁴, with a view to understand thoroughly the near wake of the airfoil at low- to ultra-low Re _c . While the lift force is measured using a load cell, the detailed flow structure is captured using laser-Doppler anemometry, particle image velocimetry, and laser-induced fluorescence flow visualization. It has been found that the stall of an airfoil, characterized by a drop in the lift force, occurs at Re _c  ≥ 1.05 × 10⁴ but is absent at Re _c  = 5.3 × 10³. The observation is connected to the presence of the separation bubble at high Re _c but absence of the bubble at ultra-low Re _c , as evidenced in our wake measurements. The near-wake characteristics are examined and discussed in detail, including the vortex formation length, wake width, spanwise vorticity, wake bubble size, wavelength of K–H vortices, Strouhal numbers, and their dependence on α and Re _c .     

Influence of unsteady wake on a turbulent separation bubble

 An experimental study was made of turbulent separated and reattaching flow over a blunt body, where unsteady wake was generated by a spoked-wheel type of wake generator with cylindrical rods. The influence of unsteady wake was scrutinized by altering the rotation direction (clockwise and counter-clockwise) and the normalized passing frequency (0 ≤ St _H  ≤ 0.20). The Reynolds number based on the cylindrical rod was Re _d =375. A phase-averaging technique was employed to characterize the unsteady wake. The effect of different rotation directions, which gave a significant reduction of x _R , was examined in detail. The wall pressure fluctuations on the blunt body were analyzed in terms of the spectrum and the coherence. Received: 15 January 2001 / Accepted: 17 July 2001     

Study of Density Effects in Turbulent Buoyant Jets Using Large-Eddy Simulation

Large-Eddy simulations (LES) of spatially evolving turbulent buoyant round jets have been carried out with two different density ratios. The numerical method used is based on a low-Mach-number version of the Navier–Stokes equations for weakly compressible flow using a second-order centre-difference scheme for spatial discretization in Cartesian coordinates and an Adams–Bashforth scheme for temporal discretization. The simulations reproduce the typical temporal and spatial development of turbulent buoyant jets. The near-field dynamic phenomenon of puffing associated with the formation of large vortex structures near the plume base with a varicose mode of instability and the far-field random motions of small-scale eddies are well captured. The pulsation frequencies of the buoyant plumes compare reasonably well with the experimental results of Cetegen (1997) under different density ratios, and the underlying mechanism of the pulsation instability is analysed by examining the vorticity transport equation where it is found that the baroclinic torque, buoyancy force and volumetric expansion are the dominant terms. The roll-up of the vortices is broken down by a secondary instability mechanism which leads to strong turbulent mixing and a subsequent jet spreading. The transition from laminar to turbulence occurs at around four diameters when random disturbances with a 5% level of forcing are imposed to a top-hat velocity profile at the inflow plane and the transition from jet-like to plume-like behaviour occurs further downstream. The energy-spectrum for the temperature fluctuations show both −5/3 and −3 power laws, characteristic of buoyancy-dominated flows. Comparisons are conducted between LES results and experimental measurements, and good agreement has been achieved for the mean and turbulence quantities. The decay of the centreline mean velocity is proportional to x ^−1/3 in the plume-like region consistent with the experimental observation, but is different from the x ⁻¹ law for a non-buoyant jet, where x is the streamwise location. The distributions of the mean velocity, temperature and their fluctuations in the near-field strongly depend upon the ratio of the ambient density to plume density ρ_a/ρ₀. The increase of ρ_a/ρ₀ under buoyancy forcing causes an increase in the self-similar turbulent intensities and turbulent fluxes and an increase in the spatial growth rate. Budgets of the mean momentum, energy, temperature variance and turbulent kinetic energy are analysed and it is found that the production of turbulence kinetic energy by buoyancy relative to the production by shear is increased with the increase of ρ_a/ρ₀. Received 16 June 2000 and accepted 26 June 2001     

Laminar-turbulent transition in the boundary layer on cones in a hypersonic flow at high Reynolds numbers per meter

The laminar-turbulent transition is experimentally studied in boundary-layer flows on cones with a rectangular axisymmetric step in the base part of the cone and without the step. The experiments are performed in an A-1 two-step piston-driven gas-dynamic facility with adiabatic compression of the working gas with Mach numbers at the nozzle exit M _∞ = 12–14 and pressures in the settling chamber P₀ = 60–600 MPa. These values of parameters allow obtaining Reynolds numbers per meter near the cone surface equal to Re _1e = (53–200) · 10⁶ m ⁻¹. The transition occurs at Reynolds numbers Re _tr = (2.3–5.7) · 10⁶. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 3, pp. 76–83, May–June, 2007.     

Time resolved PIV analysis of flow over a NACA 0015 airfoil with Gurney flap

A NACA 0015 airfoil with and without a Gurney flap was studied in a wind tunnel with Re _c = 2.0 × 10⁵ in order to examine the evolving flow structure of the wake through time-resolved PIV and to correlate this structure with time-averaged measurements of the lift coefficient. The Gurney flap, a tab of small length (1–4% of the airfoil chord) that protrudes perpendicular to the chord at the trailing edge, yields a significant and relatively constant lift increment through the linear range of the C _L versus α curve. Two distinct vortex shedding modes were found to exist and interact in the wake downstream of flapped airfoils. The dominant mode resembles a Kàrmàn vortex street shedding behind an asymmetric bluff body. The second mode, which was caused by the intermittent shedding of fluid recirculating in the cavity upstream of the flap, becomes more coherent with increasing angle of attack. For a 4% Gurney flap at α = 8°, the first and second modes corresponded with Strouhal numbers based on flap height of 0.18 and 0.13. Comparison of flow around ‘filled’ and ‘open’ flap configurations suggested that the second shedding mode was responsible for a significant portion of the overall lift increment.     

Aerodynamic modification of flow over bluff objects by plasma actuation

Particle image velocimetry and smoke visualization are used to study the alteration of the flow field in the wake of a bluff body by use of an alternating current (AC) surface dielectric barrier discharge. Staggered, surface, and buried electrodes were positioned on the downstream side of circular cylinders at conditions of Re _D = 1 × 10⁴−4 × 10⁴ configured to impose a force due to the ion drift that is either along or counter to the free-stream flow direction. Smoke visualization and Particle Image Velocimetry (PIV) in the wake of the flow confirms that the configuration of the surface electrodes and operation of the discharge significantly alters the location of the flow separation point and the time-averaged velocity profiles in the near and distant wake. Measurements of the vibrational and the rotational temperature using optical emission spectroscopy on the N₂ second positive system (C³Π_u–B³Π_g) indicates that the resulting plasma is highly non-equilibrium and discounts the possibility of a thermal effect on the flow separation process. The mechanism responsible for reduction or enhancement of flow separation is attributed to the streamwise force generated by the asymmetric ion wind—the direction of which is established by the electrode geometry and the local surface charge accumulated on AC cycles.     

Thermal characteristics of the wake shear layers from a slightly heated circular cylinder

In this paper, we investigate the thermal characteristics of wake shear layers generated by a slightly heated circular cylinder. Measurements of the fluctuating temperature were made in the region x/d = 0.6 to x/d = 3 (where x is the downstream distance from the cylinder axis and d is the cylinder diameter) using a single cold-wire probe. The Reynolds number Re was varied in the range 2,600–8,600. For Re = 5,500, simultaneous measurements were made with a rake of 16 cold wires, aligned in the direction of the mean shear, at x/d = 2 and 3. The results indicate that the passive temperature can be an effective marker of various instabilities of the wake shear layers, including the Kelvin–Helmholtz (KH) instability. The temperature data have confirmed the approximate Re ^m dependence of the KH instability frequency (f _KH) with different values of m over different ranges of Re, as reported previously in the literature. However, it is found that this power-law dependence is not exact, and a third-order polynomial dependence appears to fit the data well over the full range of Re. Importantly, it is found that the wake shear-layer instabilities can be grouped into three categories: (1) one with frequencies much smaller than the Bénard–Kármán-vortex shedding frequency, (2) one associated with the vortex shedding and (3) one related to the KH instability. The low-frequency shear-layer instabilities from both sides of the cylinder are in-phase, in contrast to the anti-phase high-frequency KH instabilities. Finally, the observed streamwise decrease in the mean KH frequency provides strong support for the occurrence of vortex pairing in wake shear layers from a circular cylinder, thus implying that both the wake shear layer and a mixing layer develop in similar fashion.     

Stochastic estimation of flow near the trailing edge of a NACA0012 airfoil

A stochastic estimation technique has been applied to simultaneously acquired data of velocity and surface pressure as a tool to identify the sources of wall-pressure fluctuations. The measurements have been done on a NACA0012 airfoil at a Reynolds number of Re _c  = 2 × 10⁵, based on the chord of the airfoil, where a separated laminar boundary layer was present. By performing simultaneous measurements of the surface pressure fluctuations and of the velocity field in the boundary layer and wake of the airfoil, the wall-pressure sources near the trailing edge (TE) have been studied. The mechanisms and flow structures associated with the generation of the surface pressure have been investigated. The “quasi-instantaneous” velocity field resulting from the application of the technique has led to a picture of the evolution in time of the convecting surface pressure generating flow structures and revealed information about the sources of the wall-pressure fluctuations, their nature and variability. These sources are closely related to those of the radiated noise from the TE of an airfoil and to the vibration issues encountered in ship hulls for example. The NACA0012 airfoil had a 30 cm chord and aspect ratio of 1.     

Dependence of the wake on inclination of a stationary cylinder

Three-dimensional vorticity in the wake of an inclined stationary circular cylinder was measured simultaneously using a multi-hot wire vorticity probe over a streamwise range of x/d = 10–40. The study aimed to examine the dependence of the wake characteristics on cylinder inclination angle α (=0°–45°). The validity of the independence principle (IP) for vortex shedding was also examined. It was found that the spanwise mean velocity which represents the three-dimensionality of the wake flow, increases monotonically with α. The root-mean-square (rms) values of the streamwise (u) and spanwise (w) velocities and the three vorticity components decrease significantly with the increase of α, whereas the transverse velocity (v) does not follow the same trend. The vortex shedding frequency decreases with the increase of α. The Strouhal number (St _N), obtained by using the velocity component normal to the cylinder axis, remains approximately a constant within the experimental uncertainty (±8%) when α is smaller than about 40°. The autocorrelation coefficients ρ _u and ρ _v of the u and v velocity signals show apparent periodicity for all inclination angles. With increasing α, ρ _u and ρ _v decrease and approach zero quickly. In contrast, the autocorrelation coefficient ρ _w of w increases with α in the near wake, implying an enhanced three-dimensionality of the wake.     

Scale-adaptive simulation of flow past wavy cylinders at a subcritical Reynolds number

The Xu & Yan scale-adaptive simulation (XYSAS) model is employed to simulate the flows past wavy cylinders at Reynolds number 8 × 10 3.This approach yields results in good agreement with experimental measurements.The mean flow field and near wake vortex structure are replicated and compared with that of a corresponding circular cylinder.The effects of wavelength ratios λ/D m from 3 to 7,together with the amplitude ratios a /D m of 0.091 and 0.25,are fully investigated.Owing to the wavy configuration,a maximum reduction of Strouhal number and root-meansquare (r.m.s) fluctuating lift coefficients are up to 50% and 92%,respectively,which means the vortex induced vibration (VIV) could be effectively alleviated at certain larger values of λ/D m and a /D m.Also,the drag coefficients can be reduced by 30%.It is found that the flow field presents contrary patterns with the increase of λ/D m.The free shear layer becomes much more stable and rolls up into mature vortex only further downstream when λ/D m falls in the range of 5-7.The amplitude ratio a /D m greatly changes the separation line,and subsequently influences the wake structures.     

Self-consistent modeling of entangled network strands and linear dangling structures in a single-strand mean-field slip-link model

Linear viscoelastic (LVE) measurements as well as non-linear elongation measurements have been performed on stoichiometrically imbalanced polymeric networks to gain insight into the structural influence on the rheological response (Jensen et al., Rheol Acta 49(1):1–13, 2010). In particular, we seek knowledge about the effect of dangling ends and soluble structures. To interpret our recent experimental results, we exploit a molecular model that can predict LVE data and non-linear stress–strain data. The slip-link model has proven to be a robust tool for both LVE and non-linear stress–strain predictions for linear chains (Khaliullin and Schieber, Phys Rev Lett 100(18):188302–188304, 2008, Macromolecules 42(19):7504–7517, 2009; Schieber, J Chem Phys 118(11):5162–5166, 2003), and it is thus used to analyze the experimental results. Initially, we consider a stoichiometrically balanced network, i.e., all strands in the ensemble are attached to the network in both ends. Next we add dangling strands to the network representing the stoichiometric imbalance, or imperfections during curing. By considering monodisperse network strands without dangling ends, we find that the relative low-frequency plateau, G₀/G_N⁰G_0/G_N^0, decreases linearly with the average number of entanglements. The decrease from G_N⁰G_N^0 to G ₀ is a result of monomer fluctuations between entanglements, which is similar to “longitudinal modes” in tube theory. It is found that the slope of G′ is dependent on the fraction of network strands and the structural distribution of the network. The power-law behavior of G ^″ is not yet captured quantitatively by the model, but our results suggest that it is a result of polydisperse dangling and soluble structures.     

Numerical study of transient laminar natural convection heat transfer over a sphere subjected to a constant heat flux

Transient laminar natural convection over a sphere which is subjected to a constant heat flux has been studied numerically for high Grashof numbers (10⁵ ≤ Gr ≤ 10⁹) and a wide range of Prandtl numbers (Pr = 0.02, 0.7, 7, and 100). A plume with a mushroom-shaped cap forms above the sphere and drifts upward continuously with time. The size and the level of temperature of the transient cap and plume stem decrease with increasing Gr and Pr. Flow separation and an associated vortex may appear in the wake of the sphere depending on the magnitude of Gr and Pr. A recirculation vortex which appears and grows until “steady state” is attained was found only for the very high Grashof numbers (10⁵ ≤ Gr ≤ 10⁹) and the lowest Prandtl number considered (Pr = 0.02). The appearance and subsequent disappearance of a vortex was observed for Gr = 10⁹ and Pr = 0.7. Over the lower hemisphere, the thickness of both the hydrodynamic (δ_H) and the thermal (δ_T) boundary layers remain nearly constant and the sphere surface is nearly isothermal. The surface temperature presents a local maximum in the wake of the sphere whenever a vortex is established in the wake of the sphere. The surface pressure recovery in the wake of the sphere increases with decreasing Pr and with increasing Gr. For very small Pr, unlike forced convection, the ratio δ_T/δ_H remains close to unity. The results are in good agreement with experimental data and in excellent agreement with numerical results available in the literature. A correlation has also been presented for the overall Nusselt number as a function of Gr and Pr.     

Effect of sound-absorbing materials on intensity of disturbances in the shock layer on a flat plate aligned at an angle of attack

Results of a numerical and experimental study of characteristics of disturbances in a hypersonic shock layer on a flat plate covered by a sound-absorbing coating and aligned at an angle of attack are presented. Experiments and computations are performed for the free-stream Mach number M _∞ = 21 and Reynolds number Re _L = 6 · 10⁴. A possibility of suppressing pressure fluctuations in the shock layer at frequencies of 20–40 kHz with the use of tubular and porous materials incorporated into the plate surface is demonstrated. Results of numerical simulations are found to be in good agreement with experimental data.     

Shear layers of a circular cylinder with rotary oscillation

The behavior of the separated shear layers and the near wake of a circular cylinder with small-amplitude rotary oscillations (Ω₁ = 0.05−0.15 for f _f/f _o ≤ 1.25) were investigated experimentally at Re = 3,700. Measurements of an unforced cylinder were also made for 2,000 ≤ Re ≤ 10,000 to better understand the effects of rotary oscillations. The results show that the shear-layer vortices formed closer to the cylinder and the distance separating them was found to decrease with cylinder oscillations. The shear-layer frequency, however, increased with increasing forcing frequency f _f. The formation-region length l _f decreased significantly with increasing f _f while decreased to a lesser extent with increasing normalized oscillation amplitude Ω₁. The shear layer also diffused to a length L _d larger than that of an unforced cylinder, while the l _f-L _d-Strouhal frequency offsetting mechanism was generally maintained. The near wake was of lower momentum compared to an unforced cylinder, and the transverse velocity fluctuations associated with the unforced vortex-shedding frequency f _o always presented a local peak at f _f/f _o = 0.5, regardless of Ω₁ tested.     

Mach reflection of a plane shock wave passing a mountain of 45° inclination

The transition from regular reflection (RR) to Mach reflection (MR) as a plane shock wave diffracts around a triangular mountain of 45° inclination is analysed in this paper, both by optical measurement in a shock tube and by numerical simulation the numerical method developed by Li Yingfan^[1] is of the FLIC type with triangular mesh. The dependence of the critical transition point Lk ofRR→MR on shock Mach numberM _i is analyzed and the variations of the incidence angle ω _i of the impinging shock and the reflection angle ω _r with the distanceL ^* are investigated. Our experimental and numerical results agree well with the theoretical results of Iton and Italya.     

Production of nanomaterials by vaporizing ceramic targets irradiated by a moderate-power continuous-wave CO<Subscript>2</Subscript> laser

The efficiency of utilization of CO ₂ laser energy for vaporization of Al ₂ O ₃ ceramics is evaluated using a mathematical model for the interaction of laser radiation with materials. It is shown that the calculated efficiency of radiation-energy utilization is not higher than 15% at a radiation power density of 10⁵ W/cm ² on the target. On the experimental facility designed for the synthesis of nanopowders, a vaporization rate of 1 g/h was achieved for Al ₂ O ₃, which corresponds to a 3% efficiency of radiation-energy utilization. The dependence of the characteristic particle size of a zirconium oxide nanopowder on helium pressure in the range of 0.01–1.00 atm was studied. Results of experiments on vaporization of multicomponent materials (LaNiO ₃ and the Tsarev meteorite) are given. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 2, pp. 172–184, March–April, 2007.