Assessment of a sponge layer as a non-reflective boundary treatment with highly accurate gust–airfoil interaction results

This paper deals with the numerical performance of a sponge layer as a non-reflective boundary condition. This technique is well known and widely adopted, but only recently have the reasons for a sponge failure been recognised, in analysis by Mani. For multidimensional problems, the ineffectiveness of the method is due to the self-reflections of the sponge occurring when it interacts with an oblique acoustic wave. Based on his theoretical investigations, Mani gives some useful guidelines for implementing effective sponge layers. However, in our opinion, some practical indications are still missing from the current literature. Here, an extensive numerical study of the performance of this technique is presented. Moreover, we analyse a reduced sponge implementation characterised by undamped partial differential equations for the velocity components. The main aim of this paper relies on the determination of the minimal width of the layer, as well as of the corresponding strength, required to obtain a reflection error of no more than a few per cent of that observed when solving the same problem on the same grid, but without employing the sponge layer term. For this purpose, a test case of computational aeroacoustics, the single airfoil gust response problem, has been addressed in several configurations. As a direct consequence of our investigation, we present a well documented and highly validated reference solution for the far-field acoustic intensity, a result that is not well established in the literature. Lastly, the proof of the accuracy of an algorithm for coupling sub-domains solved by the linear and non-liner Euler governing equations is given. This result is here exploited to adopt a linear-based sponge layer even in a non-linear computation.     

Correlation study in shock wave�Cturbulent boundary layer interaction

Shock wave–turbulent boundary layer interaction is a critical problem in aircraft design. Therefore, a thorough understanding of the processes occurring in such flows is necessary. The most important task is to study the unsteady phenomena, in particular, the low-frequency ones, for this interaction. An experimental study of separated flow has been performed in the zone of interaction of the incident oblique shock wave with a turbulent boundary layer at Mach 2. Two-point correlation data in the separation zone and the upstream flow were obtained and showed that low-frequency oscillations of the reflected shock waves are related to pulsations in the inflow turbulent boundary layer.     

Bubble–wall interaction and two-phase flow parameters on a full-scale boat boundary layer

Full scale bubbly flow experiments were performed on a 6 m flat bottom survey boat, measuring the void fraction, bubble velocity and size distributions as the bubbles naturally entrained at the bow of the boat interact with the boat’s boundary layer. Double-tip sapphire optical probes capable of measuring bubbles down to 50 μm in diameter were specifically designed and built for this experiment. The probes were positioned under the hull at the bow near the bubble entrainment region and at the stern at the exit of the bottom flat plate. Motorized positioners were used to vary the probe distance to the wall from 0 to 50 mm. The experiments were performed in fresh water (Coralville Lake, IA) and salt water (Panama City Beach, FL), at varying velocities with most data analysis performed at 10, 14 and 18 knots. The results indicate that the bubbles interact significantly with the boundary layer. At low velocity in fresh water, bubble accumulation under the hull and coalescence are evident by the presence of large bubbles at the stern. At high speeds bubble breakup dominates and very small bubbles are produced near the wall. It is also observed that salt water inhibits coalescence, even at low boat speeds. The void fraction increases with speed beyond 10 knots and peaks near the wall. Bubble velocities show slip with the wall at all speeds and exhibit large RMS fluctuations, increasing near the wall.     

A study of boundary layer–wake interaction in shallow open channel flows

This paper describes an experimental investigation of the interaction between the boundary layer on a horizontal floor of a shallow open channel flow and the wake of a thin flat plate mounted vertically on the floor of the channel. The nominal thickness of the flat plate was limited to 2 mm in order to minimize the effect of the flume side walls on the generated wake. Two flat plate chord to thickness ratios (10 and 25) and two depths of flow (50 and 80 mm) were considered. The boundary layer thickness of the approaching flow was comparable with the depth of flow. The recovery of the boundary layer is then studied by observing the characteristics of the velocity profile downstream of the flat plate along the wake axis. The results indicate that the recovery process is slow, and that it is clearly non-monotonic. When compared with the approaching flow, the streamwise turbulence intensity values increase in the near-wake region followed by a gradual but slow decrease with increasing axial distance. Neither mean nor higher-order moments indicate a complete recovery even at large distances from the wake generator. The present results also indicate that the inner region appears to develop more quickly than the outer flow. Since the development of the quasi-two-dimensional wake is also of interest, velocity measurements are also presented along the wake cross-section. These velocity profiles indicate that the wake effects are still prevalent at 200 plate widths downstream of the wake generator. Through a proper choice of scaling variables, the mean velocity profiles across the wake can be collapsed onto a single curve, indicating a sense of similarity. Received: 23 September 1999/Accepted: 30 August 2000     

“Non-free” interaction of plane shock waves and the boundary layer in the neighborhood of the leading edge of a plate with slip

The interaction between a normally impinging shock wave and the boundary layer on a plate with slip is studied in the neighborhood of the leading edge using various experimental methods, including special laser technology, to visualize the supersonic conical gas flows. It is found that in the “non-free” interaction, when the leading edge impedes the propagation of the boundary layer separation line upstream, the structure of the disturbed flow is largely identical to that in the developed “free” interaction, but with higher parameter values and gradients in the leading part of the separation zone. The fundamental property of developed separation flows, namely, coincidence of the values of the pressure “plateau” in the separation zone and the pressure behind the oblique shock above the separation zone of the turbulent boundary layer, is conserved. Moscow. e-mail: ostap@inmech.msu.su. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 57–69, May–June, 2000. The work was carried out with financial support from the Russian Foundation for Basic Research (project No. 97-01-00099).     

Experimental study of vortex–structure interaction noise radiated from rod–airfoil configurations

Vortex–structure interaction noise radiated from an airfoil embedded in the wake of a rod is investigated experimentally in an anechoic wind tunnel by means of a phased microphone array for acoustic tests and particle image velocimetry (PIV) for the flow field measurements. The rod–airfoil configuration is varied by changing the rod diameter (D), adjusting the cross-stream position (Y) of the rod and the streamwise gap (L) between the rod and the airfoil leading edge. Two noise control concepts, including “air blowing” on the upstream rod and a soft-vane leading edge on the airfoil, are applied to control the vortex–structure interaction noise. The motivation behind this study is to investigate the effects of the three parameters on the characteristics of the radiated noise and then explore the influences of the noise control concepts. Both the vortex–structure interaction noise and the rod vortex shedding tonal noise are analysed. The acoustic test results show that both the position and magnitude of the dominant noise source of the rod–airfoil model are highly dependent on the parameters considered. In the case where the vortex–structure interaction noise is dominant, the application of the air blowing and the soft vane can effectively attenuate the interaction noise. Flow field measurements suggest that the intensity of the vortex–structure interaction and the flow impingement on the airfoil leading edge are suppressed by the control methods, giving a reduction in noise.     

Thermodynamic fluctuations in canonical shock–turbulence interaction: effect of shock strength

In this work, we use numerical simulation and linear inviscid theory to study the thermodynamic field generated by the interaction of a shock wave with homogeneous isotropic turbulence. Fluctuations in density, pressure, temperature and entropy can play an important role in shock-induced mixing, combustion and energy transfer processes. Data from shock-captured direct numerical simulations (scDNS) are used to investigate the variation of thermodynamic fluctuations for varying shock strengths, and the results are compared with linear interaction analysis (LIA). The density, pressure and temperature variances attain large values at the shock, followed by, in general, a rapid decay in the downstream flow. The rapid variation behind the shock makes it difficult to compare numerical results with theoretical predictions. A threshold method based on instantaneous shock dilatation is used to overcome this problem, and it gives excellent match between scDNS and LIA. We find cases with non-monotonic variation with Mach number as well as local peaks in density fluctuations behind the shock. These are explained in terms of the contribution of the post-shock acoustic and entropy modes in the LIA solution and their cross-correlation. Budget of the transport equations reveals interesting insight into the physics governing the thermodynamic field behind the shock wave. It is found that the variances are primarily determined by the competing effects of dilatational and dissipation mechanisms. The dominant mechanisms are identified for a range of conditions, and their implication for developing predictive models is highlighted.     

“Thermal layer” effect in MHD: Interaction of a parallel shock with a layer of decreased density

Interaction of a parallel fast MHD shock with a layer of decreased density is discussed using ideal MHD approach. This is an extrapolation of gas dynamic thermal layer effect on ideal MHD. Computer simulations show that a magnetic field of a moderate intensity ( 1) may change the character of the flow for intermediate Mach numbers (M 5) and a new raking regime may occur which is not observed in the absence of a magnetic field. Self similar precursor analogous to that in gas dynamics may develop in the case of highM and low density in the layer but magnetic forces essentially decrease its growth rate. This problem appears in connection with cosmical shock propagation where planetary magnetic tails play the role of the thermal layer, and it may also be observed in the laboratory when the shock is strong enough to heat the walls ahead of it.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Experimental study of the shock–bubble interaction with reshock

The interaction of a planar shock wave with a spherical density inhomogeneity is studied experimentally under reshock conditions. Reshock occurs when the incident shock wave, which has already accelerated the spherical bubble, reflects off the tube end wall and reaccelerates the inhomogeneity for a second time. These experiments are performed at the Wisconsin Shock Tube Laboratory, in a 9m-long vertical shock tube with a large square cross section (25.4×25.4 cm²). The bubble is prepared on a pneumatically retracted injector and released into a state of free fall. Planar diagnostic methods are used to study the bubble morphology after reshock. Data are presented for experiments involving two Atwood numbers (A = 0.17 and 0.68) and three Mach numbers (1.35 < M < 2.33). For the low Atwood number case, a secondary vortex ring appears immediately after reshock which is not observed for the larger Atwood number. The post-reshock vortex velocity is shown to be proportional to the incident Mach number, M, the initial Atwood number, A, and the incident shock wave speed, W _i.     

Modeling of the unsteady force for shock–particle interaction

The interaction between a particle and a shock wave leads to unsteady forces that can be an order of magnitude larger than the quasi-steady force in the flow field behind the shock wave. Simple models for the unsteady force have so far not been proposed because of the complicated flow field during the interaction. Here, a simple model is presented based on the work of Parmar et al. (Phil Trans R Soc A 366:2161–2175, 2008). Comparisons with experimental and computational data for both stationary spheres and spheres set in motion by shock waves show good agreement in terms of the magnitude of the peak and the duration of the unsteady force.      

Relaminarization of a turbulent boundary layer with a mach number M∞=4

Results of experimental studies are presented on relaminarization of a supersonic turbulent boundary layer behind an expansion fan for a freestream Mach number M=4 within a range of Reynolds numbers Re₁=8·10⁶ – 26·10⁶ m^–1. Experimental data on distributions of the mean velocity and massflow fluctuations and the skin friction force are obtained. Partial relaminarization of the boundary layer is reached in the experiments. The calculations of relaminarization criteria show that they can be used to predict the onset of the relaminarization process at high supersonic flow velocities.     

Time–space characteristics of stratified shear layer from UVP measurements

The stratified shear layer flow pattern involves a fresh water layer flowing over a salted water one. The instabilities arising due to velocity gradients are mainly convective and thus evolve in time and space during their downstream propagation. This work was carried out within the framework of a study on the interaction between fresh water and seawater in the estuaries of rivers, where the main part of the physicochemical and biological phenomena occurs under the control of hydrodynamic conditions. The spatio-temporal velocity profile measurement by the ultrasonic Doppler method is well adapted to the stability analysis of such flows. It allows a comparison between the experimental wave properties and the theoretical results given by a linear temporal approach extended to the spatial point of view.     

The effect of G?rtler instability on hypersonic boundary layer transition

The evolution of G?rtler vortices and its interaction with other instabilities are investigated in this paper.Both the Mack mode and the G?rtler mode exist in hypersonic boundary-layer flows over concave surfaces, and their interactions are crucially important in boundary layer transition. We carry out a direct numerical simulation to explore the interaction between the G?rtler and the oblique Mack mode.The results indicate that the interaction between the forced G?rtler mode and the oblique Mack mode promotes the onset of the transition. The forced oblique Mack mode is susceptible to nonlinear interaction.Because of the development of the G?rtler mode, the forced Mack mode and other harmonic modes are excited.     

Application of the PIV technique to measurements around and inside a forming drop in a liquid–liquid system

A particle image velocimetry (PIV) method has been developed to measure the velocity field inside and around a forming drop with a final diameter of 1 mm. The system, including a microscope, was used to image silicon oil drops forming in a continuous phase of water and glycerol. Fluorescent particles with a diameter of 1 μm were used as seeding particles. The oil was forced through a 200 μm diameter glass capillary into a laminar cross-flow in a rectangular channel. The velocity field was computed with a double-frame cross-correlation function down to a spatial resolution of 21 × 21 μm. The method can be used to calculate the shear stress induced at the interface by the cross-flow of the continuous phase and the main forces involved in the drop formation process.     

Generation of the tollmien–schlichting wave in a supersonic boundary layer by two sinusoidal acoustic waves

The problem is solved using parabolized equations of stability for threedimensional perturbations of a compressible boundary layer on a flat plate. Nonlinearity is taken into account by quadratic terms that are most significant in estimates of the viscous critical layer of the stability theory. These terms are determined by the total field of two acoustic perturbations, and the equations become linear and inhomogeneous. The calculations are performed for one acoustic wave being stationary in the reference system fitted to the plate for Mach numbers M=2 and 5. Solutions are presented, which are identified very accurately with Tollmien–Schlichting waves at a rather large distance from the plate edge.     

Effect of shock–wave loading on the properties of cryogenic TiNi — a shape–memory alloy

Results of an experimental study of the shock–wave deformation of TiNi and its effect on the crystallographic structure and temperature of austenite–martensite transformations are given. It is found that, for pressures of up to 2 GPa, shock–wave loading changes the defect structure and parameters of the lattice; however, this does not lead to a noticeable change in the temperature of the austenite–martensite transformation and the manifestation of the shapeNdash;memory effect.     

Numerical simulation of shock–vortex interaction in Schardin’s problem

Shock diffraction over a two-dimensional wedge and subsequent shock–vortex interaction have been numerically simulated using the AUSM $+$ + scheme. After the passage of the incident shock over the wedge, the generated tip vortex interacts with a reflected shock. The resulting shock pattern has been captured well. It matches the existing experimental and numerical results reported in the literature. We solve the Navier–Stokes equations using high accuracy schemes and extend the existing results by focussing on the Kelvin–Helmholtz instability generated vortices which follow a spiral path to the vortex core and on their way interact with shock waves embedded within the vortex. Vortex detection algorithms have been used to visualize the spiral structure of the initial vortex and its final breakdown into a turbulent state. Plotting the dilatation field we notice a new source of diverging acoustic waves and a lambda shock at the wedge tip.